Galactic transient black hole candidate (BHC) MAXI~J1535-571 was discovered on 2017 September 02 simultaneously by {\it MAXI}/GSC and {\it Swift}/BAT instruments. It has also been observed by India's first multi-wavelength astronomy-mission satellite {\it AstroSat}, during the rising phase of its 2017-18 outburst. We make both the spectral and the temporal analysis of the source during 2017 September 12-17 using data of {\it AstroSat}'s Large Area X-ray Proportional Counter (LAXPC) in the energy range of $3-40$~keV to infer the accretion flow properties of the source. Spectral analysis is done with the physical two-component advective flow (TCAF) solution-based {\it fits} file. From the nature of the variation of the TCAF model fitted physical flow parameters, we conclude and confirm that during our analysis period source was in the soft-intermediate spectral state. We observe sharp type-C quasi-periodic oscillations (QPOs) in the frequency range of $\sim 1.74-2.81$~Hz, which are not common in the soft-intermediate spectral state. For a better understanding of the nature and evolution of these type-C QPOs, a dynamic study of the power density spectra is done. We also investigate the origin of these QPOs from the shock oscillation model. We find that as is expected in the soft-Intermediate spectral states, non-satisfaction of Rankine-Hugoniot conditions for non-dissipative shocks and not their resonance oscillations may be the cause of the origin of type-C QPOs.

Instruments flown on CubeSats are small. Meaningful applications of CubeSats in astronomical observations rely on the choice of a particular subject that is feasible for CubeSats. Here we report the result of a feasibility study for observing gamma-ray polarization from Cygnus X-1 using a small Compton polarimeter on board a 3U CubeSat. Silicon detectors and cerium bromide scintillators were employed in the instrument models that we discussed in this study. Through Monte Carlo simulations with Geant4-based MEGAlib package, we found that, with a 10-Ms on-axis, zenith-direction observation in a low inclination, low altitude earth orbit radiation background environment, the minimum detectable polarization degree can be down to about 10\% in 160 - 250 keV, 20\% in 250 - 400 keV, and 65\% in 400 - 2000 keV. A 3U CubeSat dedicated to observing Cygnus X-1 can therefore yield useful information on the polarization state of gamma-ray emissions from the brightest persistent X-ray black-hole binary in the sky.

The 2015 Outburst of V404 Cygni is an unusual one with several X-ray and radio flares and rapid variation in the spectral and timing properties. The outburst occurred after $26$ years of inactivity of the black hole. We study the accretion flow properties of the source during its initial phase of the outburst using {\it Swift}/XRT and {\it Swift}/BAT data in the energy range of $0.5-150$ keV. We have done spectral analysis with the two-component advective flow (TCAF) model fits file. Several flow parameters such as two types of accretion rates (Keplerian disk and sub-Keplerian halo), shock parameters (location and compression ratio) are extracted to understand the accretion flow dynamics. We calculated equipartition magnetic field $B$ for the outburst and found that the highest $B \sim 900$~Gauss. Power density spectra (PDS) showed no break, which indicates no or very less contribution of the Keplerian disk component, which is also seen from the result of the spectral analysis. No signature of prominent quasi-periodic oscillations (QPOs) is observed in the PDS. This is due to the mismatch of the cooling timescale and infall timescale of the post-shock matter.

Application of cubesats in astronomical observations has been getting more and more mature in recent years. Here we report a concept study of a small Compton polarimeter to fly on a cubesat for observing polarization of soft gamma-rays from a black-hole X-ray binary, Cygnus X-1. Polarization states provide very useful diagnostics on the emission mechanism and the origin of those gamma rays. In our study, we conducted Monte Carlo simulations to decide the basic design of this small polarimeter. Silicon detectors and cerium bromide scintillators were employed in this study. We estimated its on-axis Compton efficiency at different energies and its data telemetry requirement when flying in a low earth orbit. Our results indicate that it is feasible to achieve high signal-to-noise ratio for observing Cyg X-1 with such a small instrument. Based on this study, we will proceed to have a more realistic design and look for opportunities of a cubesat space mission.

In this Astro2020 APC White Paper, we describe a Small Explorer (SMEX) mission concept called the Compton Spectrometer and Imager. COSI is a Compton telescope that covers the bandpass often referred to as the "MeV Gap" because it is the least explored region of the whole electromagnetic spectrum. COSI provides a significant improvement in sensitivity along with high-resolution spectroscopy, enabling studies of 511 keV electron-positron annihilation emission and measurements of several radioactive elements that trace the Galactic history of supernovae. COSI also measures polarization of gamma-ray bursts (GRBs), accreting black holes, and pulsars as well as detecting and localizing multimessenger sources. In the following, we describe the COSI science, the instrument, and its capabilities. We highlight many Astro2020 science WPs that describe the COSI science in depth.

We build HI absorption spectra towards Supernova Remnant (SNRs) G16.7+0.1 and G15.9+0.2 using the THOR survey data. With the absorption spectra, we give a new distance range of 7 to 16 kpc for G15.9+0.2. We also resolve the near/far-side distance ambiguity of G16,7+0.1 and confirm its kinematic distance of about 14 kpc. In addition, we analyze the CO (J=3-2) spectra towards G16.7+0.1 and find obvious CO emission at the 20 kms$^{-1}$ OH 1720 MHz maser site. This supports Reynoso and Mangum (2000)'s suggestions that the velocity difference between the maser and southern molecular cloud is caused by the shock acceleration. We discuss the impact of the distances on other physical parameters of the two SNRs.

We present spectral and timing properties of the newly discovered X-ray transient source, MAXI J1535-571, which is believed to be a Galactic X-ray binary containing a black hole candidate (BHC) as the primary object. After its discovery on 2017 Sep. 2, it has been monitored regularly in multi-wavelength bands by several satellites. We use archival data of Swift (XRT and BAT) and MAXI (GSC) satellite instruments to study accretion flow dynamics of the source during the outburst. During its outburst, the source became very bright in the sky with a maximum observed flux of $5$~Crab in the $2-10$~keV GSC band. Similar to other transient BHCs, it also shows signatures of low frequency quasi-periodic oscillations (QPOs) during the outburst. Spectral data of different instruments are fitted with the transonic flow solution based two-component advective flow (TCAF) model fits file to find the direct accretion flow parameters. Evolution of spectral states and their transitions are understood from the model fitted physical flow parameters and nature of QPOs. We also estimate probable mass of the black hole from our spectral analysis as $7.9-9.9~M_\odot$ or $8.9\pm1.0~M_\odot$.

The composition of the astrophysical relativistic jets remains uncertain. By kinetic particle-in-cell simulations, we show that the baryon component in the jet, or the so-called baryon loading effect (BLE), heavily affects relativistic jets transport dynamics in the interstellar medium. On the one hand, with the BLE, relativistic jets can transport in a much longer distance, because jet electrons draw a significant amount of energy from jet baryons via the Buneman-induced electrostatic waves and the Weibel-mediated collisionless shock; on the other hand, the jet electron phase space distribution may transform from a bottom-wide-single-peak structure to a center-wide-multiple-peak one by increasing the BLE, which largely influences the observed jet morphology. Implications for related astrophysical studies are also discussed.

We present the X-ray timing results of the new black hole candidate (BHC) MAXI J1535-571 during its 2017 outburst from Hard X-ray Modulation Telescope (\emph{Insight}-HXMT) observations taken from 2017 September 6 to 23. Following the definitions given by \citet{Belloni2010}, we find that the source exhibits state transitions from Low/Hard state (LHS) to Hard Intermediate state (HIMS) and eventually to Soft Intermediate state (SIMS). Quasi-periodic oscillations (QPOs) are found in the intermediate states, which suggest different types of QPOs. With the large effective area of \emph{Insight}-HXMT at high energies, we are able to present the energy dependence of the QPO amplitude and centroid frequency up to 100 keV which is rarely explored by previous satellites. We also find that the phase lag at the type-C QPOs centroid frequency is negative (soft lags) and strongly correlated with the centroid frequency. By assuming a geometrical origin of type-C QPOs, the source is consistent with being a high inclination system.

We present X-ray observations of PWN G16.73+0.08/SNR G16.7+0.1 using archival data of {\it Chandra} ACIS. The X-ray emission peak location of this pulsar wind nebula is found to be offset by 24 arcsec from the centre of the 1.4-GHz emission of this nebula. The X-ray nebula is elongated in the direction from the X-ray peak to the 1.4-GHz emission centre. This offset suggests that G16.73+0.08 is an evolved pulsar wind nebula interacting with the supernova remnant reverse shock. We identify a point source, CXO J182058.16-142001.5, near the location of the X-ray peak. The spectrum of the X-ray nebula can be described by an absorbed power law of photon index $0.98^{+0.79}_{-0.71}$ and hydrogen column density $N_{\rm H}=4.99^{+2.75}_{-2.28}\times 10^{22}$ cm$^{-2}$. CXO J182058.16-142001.5 is likely a pulsar. We estimate its spin-down power to be about $2.6\times 10 ^{36}$ erg s$^{-1}$. Assuming its age at 3000 and 10,000 years, its dipole magnetic field strength at the polar surface is estimated to be about $4.2 \times 10^{13}$ G and $1.3 \times 10^{13}$ G, respectively.

Astrophysical polarization measurements in the soft gamma-ray band are becoming more feasible as detectors with high position and energy resolution are deployed. Previous work has shown that the minimum detectable polarization (MDP) of an ideal Compton polarimeter can be improved by $\sim 21\%$ when an unbinned, maximum likelihood method is used instead of the standard approach of fitting a sinusoid to a histogram of azimuthal scattering angles. Here we outline a procedure for implementing this maximum likelihood approach for real, non-ideal polarimeters. As an example, we use the recent observation of GRB 160530A with the Compton Spectrometer and Imager. We find that the MDP for this observation is reduced by $20\%$ when the maximum likelihood method is used instead of the standard method.

A long duration gamma-ray burst, GRB 160530A, was detected by the Compton Spectrometer and Imager (COSI) during the 2016 COSI Super Pressure Balloon campaign. As a Compton telescope, COSI is inherently sensitive to the polarization of gamma-ray sources in the energy range 0.2-5.0 MeV. We measured the polarization of GRB 160530A using 1) a standard method (SM) based on fitting the distribution of azimuthal scattering angles with a modulation curve, and 2) an unbinned, maximum likelihood method (MLM). In both cases, the measured polarization level was below the $99\%$ confidence minimum detectable polarization levels of $72.3 \pm 0.8\%$ (SM) and $57.5 \pm 0.8\%$ (MLM). Therefore, COSI did not detect polarized gamma-ray emission from this burst. Our most constraining $90\%$ confidence upper limit on the polarization level was $46\%$ (MLM).

We search for the gamma-ray counterparts of stellar-mass black holes using long-term Fermi archive to investigate the electrostatic acceleration of electrons and positrons in the vicinity of the event horizon, by applying the pulsar outer-gap model to their magnetosphere. When a black hole transient (BHT) is in a low-hard or quiescent state, the radiatively inefficient accretion flow cannot emit enough MeV photons that are required to sustain the force-free magnetosphere in the polar funnel via two-photon collisions. In this charge-starved gap region, an electric field arises along the magnetic field lines to accelerate electrons and positrons into ultra-relativistic energies. These relativistic leptons emit copious gamma-rays via the curvature and inverse-Compton (IC) processes. It is found that these gamma-ray emissions exhibit a flaring activity when the plasma accretion rate stays typically between 0.01 and 0.005 percent of the Eddington value for rapidly rotating, stellar-mass black holes. By analyzing the detection limit determined from archival Fermi/LAT data, we find that the 7-year averaged duty cycle of such flaring activities should be less than 5% and 10% for XTE J1118+480 and 1A 0620-00, respectively, and that the detection limit is comparable to the theoretical prediction for V404 Cyg. It is predicted that the gap emission can be discriminated from the jet emission, if we investigate the high-energy spectral behaviour or observe nearby BHTs during deep quiescence simultaneously in infrared wavelength and very-high energies.

We investigate the acceleration of electrons and positrons by magnetic-field-aligned electric fields in the polar funnel of an accreting black hole (BH). Applying the pulsar outer-gap theory to BH magnetospheres, we find that such a lepton accelerator arises in the immediate vicinity of the event horizon due to frame-dragging, and that their gamma-ray luminosity increases with decreasing accretion rate. Furthermore, we demonstrate that the gamma-ray flux is enhanced along the rotation axis by more than an order of magnitude if the BH spin increases from $a=0.90M$ to $a=0.9999M$. As a result, if a ten-solar-mass, almost-maximally rotating BH is located within 3 kpc, when its accretion rate is between 0.005% and 0.01% of the Eddington rate, its high-energy flare becomes detectable with the Fermi/Large Area Telescope, provided that the flare lasts longer than 1.2 months and that we view the source nearly along the rotation axis. In addition, its very-high-energy flux is marginally detectable with the Cherenkov Telescope Array, provided that the flare lasts longer than a night and that our viewing angle is about 45 degrees with respect to the rotation axis.

Around a rapidly rotating black hole (BH), when the plasma accretion rate is much less than the Eddington rate, the radiatively inefficient accretion flow (RIAF) cannot supply enough MeV photons that are capable of materializing as pairs. In such a charge-starved BH magnetosphere, the force-free condition breaks down in the polar funnels. Applying the pulsar outer-magnetospheric lepton accelerator theory to super-massive BHs, we demonstrate that a strong electric field arises along the magnetic field lines in the direct vicinity of the event horizon in the funnels, that the electrons and positrons are accelerated up to 100~TeV in this vacuum gap, and that these leptons emit copious photons via inverse-Compton (IC) process between 0.1~TeV and 30~TeV for a distant observer. It is found that these IC fluxes will be detectable with Imaging Atmospheric Cherenkov Telescopes, provided that a low-luminosity active galactic nucleus is located within 1~Mpc for a million-solar-mass central BH or within 30~Mpc for a billion-solar-mass central BH. These very-high-energy fluxes are beamed in a relatively small solid angle around the rotation axis because of the inhomogeneous and anisotropic distribution of the RIAF photon field, and show an anti-correlation with the RIAF submillimeter fluxes. The gap luminosity little depends on the three-dimensional magnetic-field configuration, because the Goldreich-Julian charge density, and hence the exerted electric field is essentially governed by the frame-dragging effect, not by the magnetic field configuration.

The Compton Spectrometer and Imager (COSI) is a balloon-borne gamma-ray (0.2-5 MeV) telescope with inherent sensitivity to polarization. COSI's main goal is to study astrophysical sources such as $\gamma$-ray bursts, positron annihilation, Galactic nucleosynthesis, and compact objects. COSI employs a compact Compton telescope design utilizing 12 high-purity cross strip germanium detectors (size: $8\times8\times1.5$ cm$^3$, 2 mm strip pitch). We require well-benchmarked simulations to simulate the full instrument response used for data analysis, to optimize our analysis algorithms, and to better understand our instrument and the in-flight performance. In order to achieve a reasonable agreement, we have built a comprehensive mass model of the instrument and developed a detailed detector effects engine, which takes into account the individual performance of each strip as well as the characteristics of the overall detector system. We performed detailed Monte-Carlo simulations with Cosima/Geant4, applied the detector effects engine, and benchmarked the results with pre-flight calibrations using radioactive sources. After applying the detector effects engine, the simulations closely resemble the measurements, and the standard calibration, event reconstruction, and imaging pipeline used for measurements can also be applied to the simulations. In this manuscript, we will describe the detector effects engine, the benchmarking tests with calibrations, and the application to preliminary results from COSI's 46-day balloon flight in 2016.

The Compton Spectrometer and Imager (COSI) is a balloon-borne, soft-gamma ray imager, spectrometer, and polarimeter with sensitivity from 0.2 to 5 MeV. Utilizing a compact Compton telescope design with twelve cross-strip, high-purity germanium detectors, COSI has three main science goals: study the 511 keV positron annihilation line from the Galactic plane, image diffuse emission from stellar nuclear lines, and perform polarization studies of gamma-ray bursts and other extreme astrophysical environments. COSI has just completed a successful 46-day flight on NASA's new Super Pressure Balloon, launched from Wanaka, New Zealand, in May 2016. We present an overview of the instrument and the 2016 flight, and discuss COSI's main science goals, predicted performance, and preliminary results.

We examine a quintessence model with a modified exponential potential given by $V(\phi) = V_0(1+e^{-\lambda \phi})$. Unlike quintessence with a standard exponential potential, our model can yield an acceptable accelerated expansion at late times, while producing a distinct "early dark energy" signature at high redshift. We determine the evolution of the equation of state parameter, $w_\phi$, and the density parameter, $\Omega_\phi$, as a function of the scale factor. The strongest constraints on the model come from cosmic microwave background observations rather than supernova data. The former give the limit $\lambda > 13$. This model predicts a value of the effective neutrino number during Big Bang nucleosynthesis larger than the standard model value. It also provides a partial solution to the coincidence problem, in that the ratio of the quintessence energy density is always within a few orders of magnitude of the background radiation or matter density from the early universe up to the present, but it does not explain why the accelerated expansion is beginning near the present day, suggesting that these two different ways of characterizing the coincidence problem are not entirely equivalent.

We investigate the electrostatic acceleration of electrons and positrons in the vicinity of the event horizon, applying the pulsar outer-gap model to black hole magnetospheres. During a low accretion phase, the radiatively inefficient accretion flow (RIAF) cannot emit enough MeV photons that are needed to sustain the force-free magnetosphere via two-photon collisions. In such a charge-starved region (or a gap), an electric field arises along the magnetic field lines to accelerate electrons and positrons into ultra-relativistic energies. These relativistic leptons emit copious gamma-rays via curvature and inverse-Compton (IC) processes. Some of such gamma-rays collide with the submillimeter-IR photons emitted from the RIAF to materialize as pairs, which polarize to partially screen the original acceleration electric field. It is found that the gap gamma-ray luminosity increases with decreasing accretion rate. However, if the accretion rate decreases too much, the diminished RIAF soft photon field can no longer sustain a stationary pair production within the gap. As long as a stationary gap is formed, the magnetosphere becomes force-free outside the gap by the cascaded pairs, irrespective of the BH mass. If a nearby stellar-mass black hole (BH) is in quiescence, or if a galactic intermediate-mass BH is in a very low accretion state, its curvature and IC emissions are found to be detectable with Fermi/LAT and imaging atmospheric Cherenkov telescopes (IACT). If a low-luminosity active galactic nucleus is located within a few tens of Mpc, the IC emission from its super-massive BH is marginally detectable with IACT.

A sample of the LAMOST spectra of the early type M0-M3 dwarfs is compared with the Kepler observations. It is found that M dwarfs with strong chromospheric emission in $H_{\alpha}$ have large flare activity in general. The rotational periods derived from the Kepler measurements have close correlations with the sizes of the flares, the power-law distribution index and the equivalent widths of the $H_{\alpha}$ emission. A clear trend exists for higher magnetic activities being detected in faster rotating M dwarfs (rotation periods $<$ 20 day).

Using Gaussian Mixture Model (GMM) and Expectation Maximization Algorithm, we perform an analysis of time duration ($T_{90}$) for \textit{CGRO}/BATSE, \textit{Swift}/BAT and \textit{Fermi}/GBM Gamma-Ray Bursts. The $T_{90}$ distributions of 298 redshift-known \textit{Swift}/BAT GRBs have also been studied in both observer and rest frames. Bayesian Information Criterion has been used to compare between different GMM models. We find that two Gaussian components are better to describe the \textit{CGRO}/BATSE and \textit{Fermi}/GBM GRBs in the observer frame. Also, we caution that two groups are expected for the \textit{Swift}/BAT bursts in the rest frame, which is consistent with some previous results. However, \textit{Swift} GRBs in the observer frame seem to show a trimodal distribution, of which the superficial intermediate class may result from the selection effect of \textit{Swift}/BAT.

Using all the RXTE archival data of Sco X-1 and GX 5-1, which amount to about 1.6 mega seconds in total, we searched for possible occultation events caused by Oort Cloud Objects. The detection efficiency of our searching approach was studied with simulation. Our search is sensitive to object size of about 300 m in the inner Oort Cloud, taking 4000 AU as a representative distance, and of 900 m in the outer Oort Cloud, taking 36000 AU as the representative distance. No occultation events were found in the 1.6 Ms data. We derived upper limits to the number of Oort Cloud Objects, which are about three orders of magnitude higher than the highest theoretical estimates in the literature for the inner Oort Cloud, and about six orders higher for the outer Oort Cloud. Although these upper limits are not constraining enough, they are the first obtained observationally, without making any model assumptions about comet injection. They also provide guidance to such serendipitous occultation event search in the future.

Ground-based observations at terahertz (THz) frequencies are a newly explorable area of astronomy for the next ten years. We discuss science cases for a first-generation 10-m class THz telescope, focusing on the Greenland Telescope as an example of such a facility. We propose science cases and provide quantitative estimates for each case. The largest advantage of ground-based THz telescopes is their higher angular resolution (~ 4 arcsec for a 10-m dish), as compared to space or airborne THz telescopes. Thus, high-resolution mapping is an important scientific argument. In particular, we can isolate zones of interest for Galactic and extragalactic star-forming regions. The THz windows are suitable for observations of high-excitation CO lines and [N II] 205 um lines, which are scientifically relevant tracers of star formation and stellar feedback. Those lines are the brightest lines in the THz windows, so that they are suitable for the initiation of ground-based THz observations. THz polarization of star-forming regions can also be explored since it traces the dust population contributing to the THz spectral peak. For survey-type observations, we focus on ``sub-THz'' extragalactic surveys, whose uniqueness is to detect galaxies at redshifts z ~ 1--2, where the dust emission per comoving volume is the largest in the history of the Universe. Finally we explore possibilities of flexible time scheduling, which enables us to monitor active galactic nuclei, and to target gamma-ray burst afterglows. For these objects, THz and submillimeter wavelength ranges have not yet been explored.

In this paper, we present the time evolution of a rotationally axisymmetric gas ring around a non rotating black hole using two dimensional grid-based hydrodynamic simulation. We show the way in which angular momentum transport is included in simulations of non-self-gravitating accretion of matter towards a black hole. We use the Shakura-Sunyaev {\alpha} viscosity prescription to estimate the turbulent viscosity for all major viscous stress tensors. We investigate how a gas ring which is initially assumed to rotate with Keplerian angular velocity is accreted on to a black hole and hence forms accretion disc in the presence of turbulent viscosity. We show that the centrifugal pressure supported sub-Keplerian flow with shocks forms when the ring starts to disperse with inclusion of relatively small amount of viscosity. But, if the viscosity is above the critical value, the shock disappears altogether and the whole disc becomes Kepleiran which is subsonic everywhere except in a region close to the horizon, where it supersonically enters to the black hole. We discovered a multiple valued Mach number solution and the corresponding density distributions that connects matter (a) from the initial Keplerian gas ring to a sub-Keplerian disc with shocks in presence of small amount of viscosity and (b) from the sub-Keplerian flow to a Keplerian disc in presence of huge amount of viscosity. We calculate the temporal variations of the magnitude of various time scales which ensure us about the stability of the flow.

Considering the contribution of the emission from the host galaxies of gamma-ray bursts (GRBs) to the radio afterglows, we investigate the effect of host galaxies on observations statistically. For the three types of events, e.g. low-luminosity, standard and high-luminosity GRBs, it is found that a tight correlation exists between the ratio of the radio flux (RRF) of host galaxy to the total radio peak emission and the observational frequency. Especially, toward lower frequencies, the contribution from the host increases significantly. The correlation can be used to get a useful estimate for the radio brightness of those host galaxies which only have very limited radio afterglow data. Using this prediction, we re-considered the theoretical radio afterglow light curves for four kinds of events, i.e. high-luminosity, low-luminosity, standard and failed GRBs, taking into account the contribution from the host galaxies and aiming at exploring the detectability of these events by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Lying at a typical redshift of $z=1$, most of the events can be detected by FAST easily. For the less fierce low-luminosity GRBs, their radio afterglows are not strong enough to exceed the sensitivity limit of FAST at such distances. However, since a large number of low luminosity bursts actually happen very near to us, it is expected that FAST will still be able to detect many of them.

Pulsars, especially millisecond pulsars, are intrinsically very stable celestial clocks, and their great pulse period stability open up a wide range of potential applications to astronomical phenomena, such as a natural detector for very low frequency ($10^{-7}-10^{-9}$ Hz) gravitational waves (GWs) background from supermassive black hole binaries. Double neutron star (DNS) binary systems, containing one or two radio pulsars, lose orbital energy by gravitational radiation, which leads to the orbital shrink. As a result, two neutron stars get closer and closer, during which it contributes to the emission of high frequency GWs of $1-10^4$ Hz. In this paper, we investigate the frequency shift of pulse signal for radio pulsars in DNS system that is induced by the emission of GWs from the system. We point out that the pulse frequency shift of radio signal in these systems can be a potential tool to hunt for the high-frequency GWs, with massive gravitons, from DNS systems, which resorts to a temporal shift of gravitational constant $\delta G(t)/G$. The sensitivity to high-frequency GWs from DNS by radio pulse frequency shift is discussed. The correlation between timing residuals of pulsar pair in double radio pulsars, such as the system PSR 0737-3039 A(B), is also considered.

Soft gamma repeaters (SGRs) and anomalous X-ray pulsars (AXPs) are two groups of enigmatic objects, which have been extensively investigated in past few decades. Based on the ample information about their timing behaviors, spectra, and variability properties, it was proposed that SGRs/AXPs are isolated neutron stars (NSs) with extremely strong magnetic fields, the so-called magnetars. Nonetheless, some alternative models are probably equally convincing such as those proposing that they are accreting NSs with a fall-back disk or rotation-powered magnetized and massive white dwarfs. The nature and nurture of SGRs/AXPs remain controversial. In this paper, we propose that SGRs/AXPs can originate from normal NSs in binary systems. SGRs are a class of objects containing a neutron-drip core and denser crust with a stiffer equation of state, which is formed from re-explosion of normal NSs in binary systems. It is the accretion onto the denser crust that contributes to the observed hard emissions. AXPs are a group of NSs which reserve huge thermal energy in an insulating layer above the surface. This structure accretes material from a massive disk, manifesting themselves via release of deposited energy. The spin-period clustering is due to either the brake of a slowly rotating envelope or the frictional drag during the common-envelope phase.

A close high-mass binary system consisting of a neutron star (NS) and a massive OB supergiant companion is expected to lead to a TZO structure, which consists of a NS core and a stellar envelope. We use the scenario machine program to calculate the formation tracks of TZOs in close high mass NS binaries and their subsequent evolution. We propose and demonstrate that the explosion and instant contraction of a TZO structure leave its stellar remnant as a soft gamma ray repeater and an anomalous X-ray pulsar respectively.

We report results obtained during the characterization of a commercial front-illuminated progressive scan interline transfer CCD camera. We demonstrate that the unmodified camera operates successfully in temperature and pressure conditions (-40C, 4mBar) representative of a high altitude balloon mission. We further demonstrate that the centroid of a well-sampled star can be determined to better than 2% of a pixel, even though the CCD is equipped with a microlens array. This device has been selected for use in a closed-loop star-guiding and tip-tilt correction system in the BIT-STABLE balloon mission.

In accreting neutron star (NS) low-mass X-ray binary (LMXB) systems, NS accretes material from its low-mass companion via a Keplerian disk. In a viscous accretion disk, inflows orbit the NS and spiral in due to dissipative processes, such as the viscous process and collisions of elements. The dynamics of accretion flows in the inner region of an accretion disk is significantly affected by the rotation of NS. The rotation makes NS, thus the space-time metric, deviate from the originally spherical symmetry, and leads to gravitational quadrupole, on one hand. On the other hand, a rotating NS drags the local inertial frame in its vicinity, which is known as the rotational frame-dragging effect. In this paper, {\bf we investigate the orbital motion of accretion flows of accreting NS/LMXBs and demonstrate that the rotational effects of NS result in a band of quasi-quantized structure in the inner region of the accretion disk, which is different, in nature, from the scenario in the strong gravity of black hole arising from the resonance for frequencies related to epicyclic and orbital motions. We also demonstrate that such a disk structure may account for frequencies seen in X-ray variability, such as quasi-periodic oscillations, and can be a potential promising tool of investigation for photon polarization.

We examine models in which the accelerated expansion of the universe is driven by a scalar field rolling near an inflection point in the potential. For the simplest such models, in which the potential is of the form V(\phi) = V_0 + V_3 (\phi-\phi_0)^3, the scalar field can either evolve toward \phi = \phi_0 at late times, yielding an asymptotic de Sitter expansion, or it can transition through the inflection point, producing a transient period of acceleration. We determine the parameter ranges which produce each of these two possibilities and also map out the region in parameter space for which the equation of state of the scalar field is close to -1 at all times up to the present, mimicking \LambdaCDM. We show that the latter can be consistent with either eternal or transient acceleration. More complicated inflection point models are also investigated.

We explore the detection condition of a wide-separation planet through the perturbation induced by the planetary caustic for various microlensing parameters, especially for the size of the source stars. By constructing the fractional deviation maps at various positions in the space of microlensing parameters, we find that the pattern of the fractional deviation depends on the ratio of the source radius to the caustic size, and the ratio satisfying the observational threshold varies with the star-planet separation. We have also obtained the upper limits of the source size that allows the detection of the signature of the host star as a function of the separation for given observational threshold. It is shown that this relation further leads one to a simple analytic condition for the star-planet separation to detect the boundness of wide-separation planets as a function of the mass ratio and the source radius. For example, when 5% of the detection threshold is assumed, for a source star with the radius of ~1 R_sun, an Earth-mass planet and a Jupiter-mass planet can be recognized of its boundness when it is within the separation range of ~10 AU and ~30 AU, respectively. We also compare the separation ranges of detection by the planetary caustic with those by the central caustic. It is found that when the microlensing light curve caused by the planetary caustic happens to be analyzed, one may afford to support the boundness of the wide-separation planet farther than when that caused by the central caustic is analyzed. Finally, we conclude by briefly discussing the implication of our findings on the next-generation microlensing experiments.

To study the population properties of small, remote objects beyond Neptune's orbit in the outer solar system, of kilometer size or smaller, serendipitous occultation search is so far the only way. For hectometer-sized Trans-Neptunian Objects (TNOs), optical shadows actually disappear because of diffraction. Observations at shorter wave lengths are needed. Here we report the effort of TNO occultation search in X-rays using RXTE/PCA data of Sco X-1 taken from June 2007 to October 2011. No definite TNO occultation events were found in the 334 ks data. We investigate the detection efficiency dependence on the TNO size to better define the sensible size range of our approach and suggest upper limits to the TNO size distribution in the size range from 30 m to 300 m. A list of X-ray sources suitable for future larger facilities to observe is proposed.

We report the energy dependence of normal branch oscillations (NBOs) in Scorpius X-1, a low-mass X-ray binary Z-source. Three characteristic quantities (centroid frequency, quality factor, and fractional root-mean-squared (rms) amplitude) of a quasi-periodic oscillation signal as functions of photon energy are investigated. We found that, although it is not yet statistically well established, there is a signature indicating that the NBO centroid frequency decreases with increasing photon energy when it is below 6-8 keV, which turns out to be positively correlated with the photon energy at the higher energy side. In addition, the rms amplitude increases significantly with the photon energy below 13 keV and then decreases in the energy band of 13-20 keV. There is no clear dependence on photon energy for the quality factor. Based on these results, we suggest that the NBO originates mainly in the transition layer.

We investigate the extraction of the rotational energy of a black hole under different accreting environment. When the accretion rate is moderate, the accretion disk consists of an outer thin disk and an inner advection-dominated accretion flow. In such a combined disk, the outer thin disk can sustain a magnetic field with moderate strength at the event horizon, leading to the formation of relativistic jets with moderate luminosity and speed via the magnetohrodynamic Penrose process. When the accretion rate increases enough, on the other hand, the disk becomes geometrically thin near the horizon. In this slim disk, the denser plasmas can sustain a stronger magnetic field than that in a combined disk, leading to the formation of jets with greater luminosity and speed via the Blandford-Znajek processs. It is discussed that the former jets are associated with the Fanaroff-Riley (FR) I galaxies and the latter with FR II galaxies.

It is known that the Z and atoll sources are two typical types of neutron-star sources in low mass X-ray binaries (LMXBs), which present very different Q-$\nu$ relations of lower kHz QPOs. We propose that the Z and atoll sources are two different phases in the evolutionary track of neutron star in LMXBs, instead of two types of distinct sources.

We take the recently published data of twin kHz quasi-period oscillations (QPOs) in neutron star (NS) lowmass X-ray binaries (LMXBs) as the samples, and investigate the morphology of the samples, which focuses on the quality factor, peak frequency of kHz QPOs, and try to infer their physical mechanism. We notice that: (1) The quality factors of upper kHz QPOs are low (2 ~ 20 in general) and increase with the kHz QPO peak frequencies for both Z and Atoll sources. (2) The distribution of quality factor versus frequency for the lower kHz QPOs are quite different between Z and Atoll sources. For most Z source samples, the quality factors of lower kHz QPOs are low (usually lower than 15) and rise steadily with the peak frequencies except for Sco X-1, which drop abruptly at the frequency of about 750 Hz. While for most Atoll sources, the quality factors of lower kHz QPOs are very high (from 2 to 200) and usually have a rising part, a maximum and an abrupt drop. (3) There are three Atoll sources (4U 1728-34, 4U 1636-53 and 4U 1608-52) of displaying very high quality factors for lower kHz QPOs. These three sources have been detected with the spin frequencies and sidebands, in which the source with higher spin frequency presents higher quality factor of lower kHz QPOs and lower difference between sideband frequency and lower kHz QPO frequency.

A new class of high-mass X-ray binaries (HMXBs) --- supergiant fast X-ray transients (SFXTs) --- are discovered by INTEGRAL, which are associated with OB supergiants and present long spin periods. Observational evidence indicates that some accreting neutron stars in HMXBs display accretion reversals. It has been suggested that the inverted torque can lead to a very slow rotator. According to three characteristic radii in wind-fed accretion, we developed a retrograde accretion scenario and divided the accretion phase into three regimes, to interpret the formation of the long spin period of SFXTs. The accretion regime in some SFXT systems can be determined by their spin and orbital periods.

Relativistic jets are launched from black hole (BH) X-ray binaries and active galactic nuclei when the disk accretion rate is below a certain limit (i.e., when the ratio of the accretion rate to the Eddingtion accretion rate, $\dot{m}$, is below about 0.01) but quenched when above. We propose a new paradigm to explain this observed coupling between the jet and the accretion disk by investigating the extraction of the rotational energy of a BH when it is surrounded by different types of accretion disk. At low accretion rates (e.g., when $\dot{m}\lesssim0.1$), the accretion near the event horizon is quasi-spherical. The accreting plasmas fall onto the event horizon in a wide range of latitudes, breaking down the force-free approximation near the horizon. To incorporate the plasma inertia effect, we consider the magnetohydrodynamical (MHD) extraction of the rotational energy from BHs by the accreting MHD fluid, as described by the MHD Penrose process. It is found that the energy extraction operates, and hence a relativistic jet is launched, preferentially when the accretion disk consists of an outer Shakura-Sunyaev disk (SSD) and an inner advection-dominated accretion flow. When the entire accretion disk type changes into an SSD, the jet is quenched because the plasmas brings more rest-mass energy than what is extracted from the hole electromagnetically to stop the extraction. Several other observed BH disk-jet couplings, such as why the radio luminosity increases with increasing X-ray luminosity until the radio emission drops, are also explained.

We report on an 105 ks Suzaku observation of the supernova remnant CTA 1 (G119.5+10.2). The Suzaku soft X-ray observation was carried out with both timing mode and imaging mode. A ~ 10' extended feature, which is interpreted as a bow-shock component of the pulsar wind nebula (PWN), is revealed in this deep observation for the first time. The nebular spectrum can be modelled by a power-law with a photon index of ~ 1.8 which suggests a slow synchrotron cooling scenario. The photon index is approximately constant across this extended feature. We compare and discuss our observations of this complex nebula with previous X-ray investigations. We do not obtain any significant pulsation from the central pulsar in the soft (0.2-12 keV) and hard (10-60 keV) X-ray data. The non-detection is mainly due to the loss of the precise imaging ability to accurately determine the source contribution. The spectra of XIS and HXD can be directly connected without a significant spectral break according to our analysis. Future observations of NuSTAR and Astro-H would be able to resolve the contamination and provide an accurate hard X-ray measurement of CTA 1.

We examine cyclic phantom models for the universe, in which the universe is dominated sequentially by radiation, matter, and a phantom dark energy field, followed by a standard inflationary phase. Since this cycle repeats endlessly, the Universe spends a substantial portion of its lifetime in a state for which the matter and dark energy densities have comparable magnitudes, thus ameliorating the coincidence problem. We calculate the fraction of time that the universe spends in such a coincidental state and find that it is nearly the same as in the case of a phantom model with a future big rip. In the limit where the dark energy equation of state parameter, w, is close to -1, we show that the fraction of time, f, for which the ratio of the dark energy density to the matter density lies between r_1 and r_2, is f = -(1+w) ln [(\sqrt{r_2} + \sqrt{1+r_2})/(\sqrt{r_1} + \sqrt{1+r_1})].

Chapter from the book "Astrophysics", p.287-308. Edited by Ibrahim Kucuk, ISBN 978-953-51-0473-5, InTech, March 3, 2012.

This work makes the first ever attempt to understand the influence of the black hole background space-time in determining the fundamental properties of the embedded relativistic acoustic geometry. To accomplish such task, the role of the spin angular momentum of the astrophysical black hole (the Kerr parameter $a$ -- a representative feature of the background black hole metric) in estimating the value of the acoustic surface gravity (the representative feature of the corresponding analogue space time) has been investigated for axially symmetric inflow of hydrodynamic fluid onto a rotating black hole. Since almost all astrophysical black holes are supposed to posses some degree of intrinsic rotation, the influence of the Kerr parameter on classical analogue models is very important to understand. For certain values of the initial boundary conditions describing the aforementioned flow, more than one acoustic horizons, namely two black hole type and one white hole type, may form, where the surface gravity may become formally infinite at the acoustic white hole. The connection between the corresponding analogue Hawking temperature with astrophysically relevant observables associated with the spectral signature has been discussed.

Using the recently published data of twin kHz quasi-period oscillations (QPOs) in neutron star low-mass X-ray binaries (LMXBs), we study the different profiles between bright Z sources and less luminous Atoll sources in a statistical way. We find the quality factors of upper kHz QPOs show a narrow distribution both for Z sources and Atoll sources, which concentrate at 7.98 and 9.75 respectively, the quality factors of lower kHz QPOs show a narrow distribution for Z sources and a broader distribution for Atoll sources, which concentrate at 5.25 and 86.22 respectively. In order to investigate the relation between the quality factor and the peakfrequency of kHz QPOs, we fit the data with power-law, linear and exponential functions, respectively. There is an obvious trend that the quality factors increase with the peak frequencies both for upper and lower QPOs. The implications of our results are discussed.

Millisecond radio pulsars have long been proposed to form from a spin-up recycling process in a binary system. In this paper we demonstrate that the accretion-induced field-decay and spin-up model for recycled pulsars can indeed produce those millisecond pulsars with relatively weak magnetic fields of $10^8-10^9$ G and short spin periods of a few milliseconds. Our results also suggest that the value of the currently observed highest spin frequency of millisecond pulsars may simply be constrained by the amount of mass available for accretion.

We study the latitudinal distribution of sunspots observed from 1874 to 2009 using the center-of-latitude (COL). We calculate COL by taking the area-weighted mean latitude of sunspots for each calendar month. We then form the latitudinal distribution of COL for the sunspots appearing in the northern and southern hemispheres separately, and in both hemispheres with unsigned and signed latitudes, respectively. We repeat the analysis with subsets which are divided based on the criterion of which hemisphere is dominant for a given solar cycle. Our primary findings are as follows: (1) COL is not monotonically decreasing with time in each cycle. Small humps can be seen (or short plateaus) around every solar maxima. (2) The distribution of COL resulting from each hemisphere is bimodal, which can well be represented by the double Gaussian function. (3) As far as the primary component of the double Gaussian function is concerned, for a given data subset, the distributions due to the sunspots appearing in two different hemispheres are alike. Regardless of which hemisphere is magnetically dominant, the primary component of the double Gaussian function seems relatively unchanged. (4) When the northern (southern) hemisphere is dominant the width of the secondary component of the double Gaussian function in the northern (southern) hemisphere case is about twice as wide as that in the southern (northern) hemisphere. (5) For the distribution of the COL averaged with signed latitude, whose distribution is basically described by a single Gaussian function, it is shifted to the positive (negative) side when the northern (southern) hemisphere is dominant. Finally, we conclude by briefly discussing the implications of these findings on the variations in the solar activity.

We explore a possibility that the daily sea-level pressure (SLP) over South Korea responds to the high-speed solar wind event. This is of interest in two aspects: First, if there is a statistical association this can be another piece of evidence showing that various meteorological observables indeed respond to variations in the interplanetary environment. Second, this can be a very crucial observational constraint since most models proposed so far are expected to preferentially work in higher latitude regions than the low latitude region studied here. We have examined daily solar wind speed ${\rm V}$, daily SLP difference ${\rm \Delta SLP}$, and daily ${\rm \log(BV^{2})}$ using the superposed epoch analysis in which the key date is set such that the daily solar wind speed exceeds 800 ${\rm kms^{-1}}$. We find that the daily ${\rm \Delta SLP}$ averaged out of 12 events reaches its peak at day +1 and gradually decreases back to its normal level. The amount of positive deviation of ${\rm \Delta SLP}$ is +2.5 hPa. The duration of deviation is a few days. We also find that ${\rm \Delta SLP}$ is well correlated with both the speed of solar wind and ${\rm \log(BV^{2})}$. The obtained linear correlation coefficients and chance probabilities with one-day lag for two cases are $r \simeq 0.81$ with $P> 99.9%$, and $r \simeq 0.84$ with $P> 99.9%$, respectively. We conclude by briefly discussing future direction to pursue.

The Nuclear Compton Telescope (NCT) is a balloon-borne Compton telescope designed for the study of astrophysical sources in the soft gamma-ray regime (200 keV--20 MeV). NCT's ten high-purity germanium crossed-strip detectors measure the deposited energies and three-dimensional positions of gamma-ray interactions in the sensitive volume, and this information is used to restrict the initial photon to a circle on the sky using the Compton scatter technique. Thus NCT is able to perform spectroscopy, imaging, and polarization analysis on soft gamma-ray sources. NCT is one of the next generation of Compton telescopes --- so-called compact Compton telescopes (CCTs) --- which can achieve effective areas comparable to COMPTEL's with an instrument that is a fraction of the size. The Crab Nebula was the primary target for the second flight of the NCT instrument, which occurred on 17--18 May 2009 in Fort Sumner, New Mexico. Analysis of 29.3 ks of data from the flight reveals an image of the Crab at a significance of 4-sigma. This is the first reported detection of an astrophysical source by a CCT.

We investigate the dependence of the amount of the observed galactic cosmic ray (GCR) influx on the solar North-South asymmetry using the neutron count rates obtained from four stations and sunspot data in archives spanning six solar cycles from 1953 to 2008. We find that the observed GCR influxes at Moscow, Kiel, Climax and Huancayo stations are more suppressed when the solar activity in the southern hemisphere is dominant compared with when the solar activity in the northern hemisphere is dominant. Its reduction rates at four stations are all larger than those of the suppression due to other factors including the solar polarity effect on the GCR influx. We perform the student's t-test to see how significant these suppressions are. It is found that suppressions due to the solar North-South asymmetry as well as the solar polarity are significant and yet the suppressions associated with the former are larger and more significant.

Magneto-solid-mechanical model of two-component, core-crust, paramagnetic neutron star responding to quake-induced perturbation by differentially rotational, torsional, oscillations of crustal electron-nuclear solid-state plasma about axis of magnetic field frozen in the immobile paramagnetic core is developed. Particular attention is given to the node-free torsional crust-against-core vibrations under combined action of Lorentz magnetic and Hooke's elastic forces; the damping is attributed to Newtonian force of shear viscose stresses in crustal solid-state plasma. The spectral formulae for the frequency and lifetime of this toroidal mode are derived in analytic form and discussed in the context of quasi-periodic oscillations of the X-ray outburst flux from quaking magnetars. The application of obtained theoretical spectra to modal analysis of available data on frequencies of oscillating outburst emission suggests that detected variability is the manifestation of crustal Alfven's seismic vibrations restored by Lorentz force of magnetic field stresses.

In view of the assumption that any planetary system is likely to be composed of more than one planet, and the multiple planet system with a large mass planet has more chance of detailed follow-up observations, the multiple planet system may be an efficient way to search for sub-Jovian planets. We compare the magnification pattern of the triple lens system with that of a best-fitted binary system composed of a star and a Jovian mass planet, and check the probability in detecting the low-mass secondary planet whose signature will be superposed on that of the primary Jovian mass planet. Detection probabilities of the low-mass planet in the triple lens system are quite similar to the probability in detecting such a low-mass planet in a binary system with a star and only a low-mass planet, which shows that the signature of a low-mass planet can be effectively detected even when it is concurrent with the signature of the more massive planet, implying that the binary superposition approximation works over a relatively broad range of planet mass ratio and separations, and the inaccuracies thereof do not significantly affect the detection probability of the lower mass secondary planet. Since the signature of the Jovian mass planet will be larger and lasts longer, thereby warranting more intensive follow-up observations, the actual detection rate of the low-mass planet in a triple system with a Jovian mass can be significantly higher than that in a binary system with a low-mass planet only. We conclude that it may be worthwhile to develop an efficient algorithm to search for `super-Earth' planets in the paradigm of the triple lens model for high-magnification microlensing events.

Serendipitous stellar occultation search is so far the only way to detect the existence of very small, very dim, remote objects in the solar system. To date, however, there are only very few reported detections for trans-Neptunian objects (TNOs) in optical bands. In the X-ray band, with the RXTE/PCA data of Sco X-1 taken from June 2007 to October 2009, we found one possible X-ray occultation event. We discuss the veracity and properties of this event, and suggest upper limits to the size distribution of TNOs at hectometer size and of Main-Belt Asteroids (MBAs) at decameter size.

We investigate an asteroseismic model of non-rotating paramagnetic neutron star with core-crust stratification of interior pervaded by homogeneous internal and dipolar external magnetic field, presuming that neutron degenerate Fermi-matter of the star core is in the state of Pauli's paramagnetic permanent magnetization caused by polarization of spin magnetic moments of neutrons along the axis of magnetic field of collapsed massive progenitor. The magnetic cohesion between metal-like crust and permanent-magnet-like core is considered as playing a main part in the dynamics of starquake. Focus is laid on the post-quake relaxation of the star by node-free torsional vibrations of highly conducting crustal solid-state plasma, composed of nuclei embedded in the degenerate Fermi-gas of relativistic electrons, about axis of magnetic field frozen in the immobile paramagnetic core. Two scenarios of these axisymmetric seismic vibrations are examined, in first of which these are considered as maintained by combined action of Lorentz magnetic and Hooke's elastic forces and in second one by solely Lorentz force. Based on the energy variational method of magneto-solid-mechanical theory of elastic continuous medium, the spectral formulae for the frequency and lifetime of this toroidal mode are obtained and discussed in the context of theoretical treatment of recently discovered quasi-periodic oscillations of the X-ray outburst flux from SGR 1806-20 and SGR 1900+14 as being produced by above seismic vibrations.

In the context of the white dwarf asteroseismology, we investigate vibrational properties of a non-convective solid star with an axisymmetric purely toroidal intrinsic magnetic field of two different shapes. Focus is laid on regime of node-free global Lorentz-force-driven vibrations about symmetry axis at which material displacements have one and the same form as those for nodeless spheroidal and torsional vibrations restored by Hooke's force of elastic shear stresses. Particular attention is given to the even-parity poloidal Alfven modes whose frequency spectra are computed in analytic form showing how the purely toroidal magnetic fields completely buried beneath the star surface can manifest itself in seismic vibrations of non-magnetic white dwarfs. The obtained spectral formulae are discussed in juxtaposition with those for Alfven modes in the solid star model with the poloidal, homogeneous internal and dipolar external, magnetic field whose inferences are relevant to Alfven vibrations in magnetic white dwarfs.

We explore conventional Neptune migration model with one additional planet of mass at 0.1-2.0 Me. This planet inhabited in the 3:2 mean motion resonance with Neptune during planet migration epoch, and then escaped from the Kuiper belt when Jovian planets parked near the present orbits. Adding this extra planet and assuming the primordial disk truncated at about 45 AU in the conventional Neptune migration model, it is able to explain the complex structure of the observed Kuiper belt better than the usual Neptune migration model did in several respects. However, numerical experiments imply that this model is a low-probability event. In addition to the low probability, two features produced by this model may be inconsistent with the observations. They are small number of low-inclination particles in the classical belt, and the production of a remnant population with near-circular and low-inclination orbit within a = 50-52 AU. According to our present study, including one extra planet in the conventional Neptune migration model as the scenario we explored here may be unsuitable because of the low probability, and the two drawbacks mentioned above, although this model can explain better several features which is hard to produce by the conventional Neptune migration model. The issues of low-probability event and the lack of low-inclination KBOs in the classical belt are interesting and may be studied further under a more realistic consideration.

Using the energy variational method of magneto-solid-mechanical theory of a perfectly conducting elastic medium threaded by magnetic field, the frequency spectrum of Lorentz-force-driven global torsional nodeless vibrations of a neutron star with Ferraro's form of axisymmetric poloidal nonhomogeneous internal and dipole-like external magnetic field is obtained and compared with that for this toroidal Alfv\'en mode in a neutron star with homogeneous internal and dipolar external magnetic field. The relevance of considered asteroseismic models to quasi-periodic oscillations of the X-ray flux during the ultra powerful outbursts of SGR 1806-20 and SGR 1900+14 is discussed.

The eigenfrequency problem of fundamental vibrational mode in a highly inhomogeneous star, modeled by self-gravitating mass of viscous liquid with singular density at the center, is considered in juxtaposition with that for Kelvin fundamental mode in the liquid star model with uniform density. Particular attention is given to the difference between spectral equations for the frequency and lifetime of f-mode in the singular and homogeneous star models. The newly obtained results are discussed in the context of theoretical asteroseismology of pre-white dwarf stars and stellar cocoons -- spherical gas-dust clouds with dense star-forming core at the center.

The two-component, core-crust, model of a neutron star with homogenous internal and dipolar external magnetic field is studied responding to quake-induced perturbation by substantially nodeless differentially rotational Alfv\'en oscillations of the perfectly conducting crustal matter about axis of fossil magnetic field frozen in the immobile core. The energy variational method of the magneto-solid-mechanical theory of a viscoelastic perfectly conducting medium pervaded by magnetic field is utilized to compute the frequency and lifetime of nodeless torsional vibrations of crustal solid-state plasma about the dipole magnetic-moment axis of the star. It is found that obtained two-parametric spectral formula for the frequency of this toroidal Alfven mode provides fairly accurate account of rapid oscillations of the X-ray flux during the flare of SGR 1806-20 and SGR 1900+14, supporting the investigated conjecture that these quasi-periodic oscillations owe its origin to axisymmetric torsional oscillations predominately driven by Lorentz force of magnetic field stresses in the finite-depth crustal region of the above magnetars.

We discuss $\gamma$-ray emissions from the outer gap accelerators of middle-aged pulsars for part of the series of our studies. A two-dimensional electrodynamic model is used to solve the distribution of accelerating electric fields with electron and positron pair creation and radiation processes in the magnetic meridional plane. We compute the curvature radiation and the synchrotron radiation by solving the evolution of the Lorentz factor and the pitch angle. The calculated spectra are compared with observed phase-averaged spectra. We also use a three-dimensional geometrical model to discuss the pulse profiles. We argue that the outer gap of middle-aged pulsars occupies the whole region between the last-open field lines and the critical magnetic field lines, which are perpendicular to the rotational axis at the light cylinder. We assume that there is no outer gap accelerator inside the light cylinder between the rotational axis and the critical magnetic field lines. For the Geminga pulsar, we demonstrate that the outward curvature radiation dominates in the spectrum above 10 MeV, while the inward synchrotron radiation dominates below 10 MeV. We find that the computed spectrum is consistent with the observations in X-ray through $\gamma$-ray bands. With the pulse morphology of the $\gamma$-ray emissions, we argue that the inclination angle and the viewing angle for the Geminga pulsar are $\alpha\sim 50^{\circ}$ and $\xi\sim 90^{\circ}$, respectively.

Within the framework of Newtonian magneto-solid-mechanics, relying on equations appropriate for a perfectly conducting elastic continuous medium threaded by a uniform magnetic field, the asteroseismic model of a neutron star undergoing axisymmetric global torsional nodeless vibrations under the combined action of Hooke's elastic and Lorentz magnetic forces is considered with emphasis on a toroidal Alfv\'en mode of differentially rotational vibrations about the dipole magnetic moment axis of the star. The obtained spectral equation for frequency is applied to $\ell$-pole identification of quasi-periodic oscillations (QPOs) of X-ray flux during the giant flares of SGR 1806-20 and SGR 1900+14. Our calculations suggest that detected QPOs can be consistently interpreted, within the framework of this model, as produced by global torsional nodeless vibrations of quaking magnetar if they are considered to be restored by the joint action of bulk forces of shear elastic and magnetic field stresses.

More reliable constraints on the microlensing optical depth comes from a better understanding of the Galactic model. Based on well-constrained Galactic bulge and disk models constructed from survey observations, such as, HST, 2MASS, and SDSS, we calculate the microlensing optical depths toward the Galactic bulge fields, and compare them with recent results of microlensing surveys. We test chi^2 statistics of microlensing optical depths expected from those models, as well as previously proposed models, using two types of data: optical depth map in (l, b) and averaged optical depth over the Galactic longitude l as a function of the latitude b. From this analysis, we find that the Galactic bulge models of 2MASS, Han & Gould (2003), and G2 of Stanek et al. (1997) show a good agreement with the microlensing optical depth profiles for all the microlensing observations, compared with E2 of Stanek et al. (1997). We find, on the other hand, that models involving an SDSS disk model produce relatively higher chi^2 values. It should be noted that modeled microlensing optical depths diverge in the low Galactic latitude, |b| < 2 arcdeg. Therefore, we suggest the microlensing observation toward much closer to central regions of the Galaxy to further test the proposed Galactic models, if it is more technically feasible than waiting for large data set of microlensing events.

Due to the high efficiency of planet detections, current microlensing planet searches focus on high-magnification events. High-magnification events are sensitive to remote binary companions as well and thus a sample of wide-separation binaries are expected to be collected as a byproduct. In this paper, we show that characterizing binaries for a portion of this sample will be difficult due to the degeneracy of the binary-lensing parameters. This degeneracy arises because the perturbation induced by the binary companion is well approximated by the Chang-Refsdal lensing for binaries with separations greater than a certain limit. For binaries composed of equal mass lenses, we find that the lens binarity can be noticed up to the separations of $\sim 60$ times of the Einstein radius corresponding to the mass of each lens. Among these binaries, however, we find that the lensing parameters can be determined only for a portion of binaries with separations less than $\sim 20$ times of the Einstein radius.

Millisecond dips in the RXTE/PCA archival data of Sco X-1 taken from 1996 to 2002 were reported recently. Those dips were found to be most likely caused by instrumental dead time but may also contain some true astronomical events, which were interpreted as the occultation of X-rays from Sco X-1 by Trans-Neptunian Objects (TNO) of 100-m size. Here we report the results of search for millisecond dip events with the new RXTE/PCA data of Sco X-1 taken in year 2007. Adopting the same selection criteria as that in the previous study, we found only 3 dip events in 72-ks data, much fewer than the 107 events found in the 560-ks data taken from 1996 to 2002 reported earlier. The new data provides more detailed information of individual `very large events' (VLEs), which is not available in the old archival data. Although the number of VLEs does not obviously increase during the occurrence of dip events, all the 3 dip events are coincident in time with VLEs that have no flags set for any of the propane or the 6 main xenon anodes. It is a strong indication of instrumental effects. No significant dips which might be real occultation by 60 -- 100 m TNOs were observed. With only 72-ks data, however, the previously proposed possibility that about 10 percent of the dip events might not be instrumental still cannot be strictly excluded. Using the absence of those anomalous VLEs as the criterion for identifying non-instrumental dip events, we found, at a lower confidence level, 4 dip events of duration 8 - 10 ms in the 72-ks data. Upper limits to the size distribution of TNOs at the small size end are suggested.

Periodicity search in gamma-ray data is usually difficult because of the small number of detected photons. A periodicity in the timing signal at other energy bands from the counterpart to the gamma-ray source may help to establish the periodicity in the gamma-ray emission and strengthen the identification of the source in different energy bands. It may, however, still be difficult to find the period directly from X-ray data because of limited exposure. We developed a procedure, by cross-checking two X-ray data sets, to find candidate periods for X-ray sources which are possible counterparts to gamma-ray pulsar candidates. Here we report the results of this method obtained with all the currently available X-ray data of 8 X-ray sources. Some tempting periodicity features were found. Those candidate periods can serve as the target periods for future search when new data become available so that a blind search with a huge number of trials can be avoided.

The Newtonian solid-mechanical theory of non-compressional spheroidal and torsional nodeless elastic vibrations in the homogenous crust model of a quaking neutron star is developed and applied to the modal classification of the quasi-periodic oscillations (QPOs) of X-ray luminosity in the aftermath of giant flares in SGR 1806-20 and SGR 1900+14. Particular attention is given to the low-frequency QPOs in the data for SGR 1806-20 whose physical origin has been called into question. Our calculations suggest that unspecified QPOs are due to nodeless dipole torsional and dipole spheroidal elastic shear vibrations.

A two-dimensional electrodynamic model is used to study particle acceleration and non-thermal emission mechanisms in the pulsar magnetospheres. We solve distribution of the accelerating electric field with the emission process and the pair-creation process in meridional plane, which includes the rotational axis and the magnetic axis. By solving the evolutions of the Lorentz factor, and of the pitch angle, we calculate spectrum in optical through $\gamma$-ray bands with the curvature radiation, synchrotron radiation, and inverse-Compton process not only for outgoing particles, but also for ingoing particles, which were ignored in previous studies. We apply the theory to the Vela pulsar. We find that the curvature radiation from the outgoing particles is the major emission process above 10 MeV bands. In soft $\gamma$-ray to hard X-ray bands, the synchrotron radiation from the ingoing primary particles in the gap dominates in the spectrum. Below hard X-ray bands, the synchrotron emissions from both outgoing and ingoing particles contribute to the calculated spectrum. The calculated spectrum is consistent with the observed phase-averaged spectrum of the Vela pulsar. We show that the observed five-peak pulse profile in the X-ray bands of the Vela pulsar is reproduced by the inward and outward emissions, and the observed double-peak pulse profile in $\gamma$-ray bands is explained by the outward emissions.

The Newtonian solid-mechanical theory of nodeless spheroidal and torsional of elastic seismic vibrations trapped in the crust of a quaking neutron star is outlined and applied to the modal classification of the quasi-periodic oscillations (QPOs) of X-ray luminosity in the aftermath of giant flares in SGR 1806-20 and SGR 1900+14. The presented analysis relies heavily on the Samuelsson-Andersson identification of the QPOs frequency from the range 30-200 Hz with those for torsional nodeless vibrations of multipole degree in the interval 2-12. Based on this identification, which is used to fix the input parameters entering the obtained spectral formulae, we compute frequency spectrum of nodeless spheroidal elastic vibrations. Particular attention is given to the low-frequency QPOs in the data for SGR 1806-20 whose physical origin has been called into question. Our calculations suggest that these unspecified QPOs are due to nodeless dipole torsional (18 Hz) and dipole spheroidal (26 Hz) elastic shear vibrations.

We present several exact solutions of the eigenfrequency problem for torsional shear vibrations in homogeneous and non-homogeneous models of the neutron star crust governed by canonical equation of solid mechanics with a restoring force of Hookean elasticity. Particular attention is given to regime of large lengthscale nodeless axisymmetric differentially rotational oscillations which are treated in spherical polar coordinates reflecting real geometry of the neutron star crust. Highlighted is the distinction between analytic forms and numerical estimates of the frequency, computed as a function of multipole degree of nodeless torsional oscillations and fractional depth of the crust, caused by different boundary conditions imposed on the toroidal field of material displacements. The relevance of considered models to quasiperiodic oscillations, recently detected during the flare of SGR 1806-20 and SGR 1900+14, is discussed.

The past three decades of investigation on nuclear physics and pulsar astrophysics have seen gradual recognition that elastodynamic approach to the continuum mechanics of nuclear matter provides proper account of macroscopic motions of degenerate Fermi-matter constituting interior of the nuclear material objects, the densest of all known today. This paper focuses on one theoretical issue of this development which is concerned with oscillatory behavior of a viscoelastic solid globe in the regime of quasistatic, force-free, non-compressional oscillations less investigated in the literature compared to oscillations in the regime of standing shear waves. We show that in this case the problem of computing frequency and lifetime of spheroidal and torsional modes of non-radial shear vibrations damped by viscosity can be unambiguously resolved by working from the energy balance equation and taking advantage of the Rayleigh's variational method. The efficiency of this method is demonstrated by solid globe models of nuclear physics and pulsar astrophysics dealing with oscillations of a spherical mass of a viscoelastic Fermi-solid with homogeneous and non-homogeneous profiles of the bulk density, the shear modulus, and the shear viscosity.

We discuss non-thermal emission mechanism of the Crab-like pulsars with both a two-dimensional electrodynamical study and a three-dimensional model. We investigate the emission process in the outer gap accelerator. In the two-dimensional electrodynamical study, we solve the Poisson equation of the accelerating electric field in the outer gap and the equation of motion of the primary particles with the synchrotron and the curvature radiation process and the pair-creation process. We show a solved gap structure which produces a consistent gamma-ray spectrum with EGRET observation. Based on the two-dimensional model, we conduct a three-dimensional emission model to calculate the synchrotron and the inverse-Compton processes of the secondary pairs produced outside the outer gap. We calculate the pulse profiles, the phase-resolved spectra and the polarization characteristics in optical to $\gamma$-ray bands to compare the observation of the Crab pulsar and PSR B0540-69. For the Crab pulsar, we find that the outer gap geometry extending from near the stellar surface to near the light cylinder produces a complex morphology change of the pulse profiles as a function of the photon energy. This predicted morphology change is quite similar with that of the observations. The calculated phase-resolved spectra are consistent with the data through optical to the $\gamma$-ray bands. We demonstrate that the 10$\sim$20 % of the polarization degree in the optical emissions from the Crab pulsar and the Vela pulsar are explained by the synchrotron emissions with the particle gyration motion.

We present a modified outer gap model to study the phase-resolved spectra of the Crab pulsar. A theoretical double peak profile of the light curve containing the whole phase is shown to be consistent with the observed light curve of the Crab pulsar by shifting the inner boundary of the outer gap inwardly to $\sim 10$ stellar radii above the neutron star surface. In this model, the radial distances of the photons corresponding to different phases can be determined in the numerical calculation. Also the local electrodynamics, such as the accelerating electric field, the curvature radius of the magnetic field line and the soft photon energy, are sensitive to the radial distances to the neutron star. Using a synchrotron self-Compton mechanism, the phase-resolved spectra with the energy range from 100 eV to 3 GeV of the Crab pulsar can also be explained.

In this paper, we improve the previous work on the MHD Alfven wave oscillation model for the neutron star (NS) kHz quasi-periodic oscillations (QPOs), and compare the model with the updated twin kHz QPO data. For the 17 NS X-ray sources with the simultaneously detected twin kHz QPO frequencies, the stellar mass M and radius R constraints are given by means of the derived parameter A in the model, which is associated with the averaged mass density of star as ${\rm \lan \rho \ran = 3M/(4\pi R^{3}) \simeq 2.4\times 10^{14} (g/cm^{3}) (A/0.7)^{2}}$, and we also compare the M-R constraints with the star equations of states. Moreover, we also discuss the theoretical maximum kHz QPO frequency and maximum twin peak separation, and some expectations on SAX J1808.4-3658 are mentioned, such as its highest kHz QPO frequency $\sim$~870 (Hz), which is about 1.4-1.5 times less than those of the other known kHz QPO sources. The estimated magnetic fields for both Z sources (about Eddington accretion rate $\med$) and Atoll sources ($ 1% \med$) are approximately $10^9$ G and $10^8$ G respectively.

Millisecond dips in the RXTE/PCA light curve of Sco X-1 were reported recently (Chang et al. 2006), which were interpreted as the occultation of X-rays from Sco X-1 caused by Trans-Neptunian Objects (TNO) of hundred-meter size. Inconclusive signatures of possible instrumental effects in many of these dip events related to high-energy cosmic rays were later found (Jones et al. 2006) and the TNO interpretation became shaky. Here we report more detailed analysis aiming at distinguishing true occultation events from those related to cosmic rays. Based on some indicative criteria derived from housekeeping data and two-channel spectral information, we suggest that about 10% of the dips are probable events of occultation. The total number of TNOs of size from 60 m to 100 m is estiamted to be about 10^{15} accordingly. Limited by the coarser time resolution of standard data modes of RXTE/PCA, however, definite results cannot be obtained. Adequately configured observations with RXTE or other new instruments in the future are very much desired.

We investigate polarization of high-energy emissions from the Crab pulsar in the frame work of the outer gap accelerator. The recent version of the outer gap, which extends from inside the null charge surface to the light cylinder, is used for examining the light curve, the spectrum and the polarization characteristics, simultaneously. The polarization position angle curve and the polarization degree are calculated to compare with the Crab optical data. We show that the outer gap model explains the general features of the observed light curve, the spectrum and the polarization by taking into account the emissions from inside of the null charge surface and from tertiary pairs, which were produced by the high-energy photons from the secondary pairs. For the Crab pulsar, the polarization position angle curve indicates that the viewing angle of the observer measured from the rotational axis is greater than 90 degrees.

To derive physical properties of the neutron star surface with observed spectra, a realistic model spectrum of neutron star surface emission is essential. Limited by computing resources, a full computation of the radiative transfer equations without the diffusion approximation has been conducted up to date only for the case of local magnetic fields being perpendicular to the stellar surface. In this paper we report the full-computation result for an arbitrary field direction. For comparison we also compute the radiative transfer equation using the diffusion approximation. For a given effective temperature, the computed spectrum with the diffusion approximation is always softer than that of a full computation at a non-negligible level. It leads to an over-estimate of the effective temperature if the diffusion approximation spectrum is employed in the spectral fitting. Other characteristics for different magnetic field orientations, such as the beaming pattern of the two polarization modes and the structure of the atmosphere, are also discussed.

We investigate a pair creation cascade in the magnetosphere of a rapidly rotating neutron star. We solve the set of the Poisson equation for the electro-static potential and the Boltzmann equations for electrons, positrons, and gamma-ray photons simultaneously. In this paper, we first examine the time-dependent nature of particle accelerators by solving the non-stationary Boltzmann equations on the two-dimensional poloidal plane in which both the rotational and magnetic axes reside. Evaluating the temperature of the heated polar cap surface, which is located near the magnetic pole, by the bombardment of gap-accelerated particles, and applying the scheme to millisecond pulsar parameters, we demonstrate that the solution converges to a stationary solution of which pair-creation cascade is maintained by the heated polar-cap emission, in a wide range of three-dimensional parameter space (period, period derivative, magnetic inclination angle). We also present the deathlines of millisecond pulsars.

We propose to determine the mass and the radius of a neutron star (NS) using three measurable mass-radius relationships, namely the ``apparent'' radius inferred from neutron star thermal emission, the gravitational redshift inferred from the absorption lines, as well as the averaged stellar mass density inferred from the orbital Keplerian frequency derived from the kilohertz quasi periodic oscillation (kHz QPO) data. We apply the method to constrain the NS mass and the radius of the X-ray sources, 1E 1207.4-5209, Aql X-1 and EXO 0748-676.

We investigate polarization of high-energy emissions from the Crab pulsar in the frame work of the outer gap accelerator, following previous works of Cheng and coworkers. The recent version of the outer gap, which extends from inside the null charge surface to the light cylinder, is used for examining the synchrotron radiations from the secondary and the tertiary pairs, which are produced outside the gap. We calculate the light curve, the spectrum and the polarization characteristics, simultaneously, by taking into account gyration motion of the particles. The polarization position angle curve and the polarization degree are calculated to compare with the Crab optical data. We demonstrate that the radiations from inside the null charge surface make outer-wing and off-pulse emissions in the light curve, and the tertiary pairs contribute to bridge emissions. The emissions from the secondary pairs explain the main features of the observed light curve and spectrum. On the other hand, both emissions from inside the null charge surface and from the tertiary pairs are required to explain the optical polarization behavior of the Crab pulsar. The energy dependence of the polarization features is expected by the present model. For the Crab pulsar, the polarization position angle curve indicates that the viewing angle of the observer measured from the rotational axis is greater than $90^{\circ}$.

We propose a method that can determine the ellipticities of source stars of microlensing events produced by binary lenses. The method is based on the fact that the products of the caustic-crossing timescale, $\Delta t$, and the cosine of the caustic incidence angle of the source trajectory, $\kappa$, of the individual caustic crossings are different for events involving an elliptical source, while the products are the same for events associated with a circular source. The product $\Delta t_\perp =\Delta t \cos\kappa$ corresponds to the caustic-crossing timescale when the incidence angle of the source trajectory is $\kappa=0$. For the unique determination of the source ellipticity, resolutions of at least three caustic crossings are required. Although this requirement is difficult to achieve under the current observational setup based on alert/follow-up mode, it will be possible with the advent of future lensing experiments that will survey wide fields continuously at high cadence. For typical Galactic bulge events, the difference in $\Delta t_\perp$ between caustic crossings is of the order of minutes depending on the source orientations and ellipticities. Considering the monitoring frequency of the future lensing surveys of $\sim 6$ times/hr and the improved photometry especially of the proposed space-based survey, we predict that ellipticity determinations by the proposed method will be possible for a significant fraction of multiple caustic-crossing binary lens events involving source stars having non-negligible ellipticities.

A two-dimensional electrodynamical model is used to study particle acceleration in the outer magnetosphere of a pulsar. The charge depletion from the Goldreich-Julian charge density causes a large electric field along the magnetic field lines. The charge particles are accelerated by the electric field and emit $\gamma$-rays via the curvature process. Some of the emitted $\gamma$-rays may collide with $X$-ray photons to make new pairs, which are accelerated again on the different field lines in the gap and proceed similar processes. We simulate the pair creation cascade in the meridional plane using the pair creation mean-free path, in which the $X$-ray photon number density is proportional to inverse square of radial distance. With the space charge density determined by the pair creation simulation, we solve the electric structure of the outer gap in the meridional plane and calculate the curvature spectrum. Because the two-dimensional model can link both gap width along the magnetic field line and trans-field thickness with the spectral cut-off energy and flux, we can diagnose both the current through the gap and inclination angle between the rotational and magnetic axes. We apply the theory to the Vela pulsar. By comparing the results with the $EGRET$ data, we rule out any cases that have a large particle injection at the outer boundary. We also suggest the inclination angle of $\alpha_{inc}\geq65^{\circ}$. The present model predicts the outer gap starting from near the conventional null charge surface for the Vela pulsar.

Recently, unprecedentedly accurate data on the orbital motion of stars in the vicinity of the Sagittarius $ {\rm A^*}$ have been available. Such information can be used not only to constrain the mass of the supermassive black hole (SMBH) in the Galactic center but also to study the source of the radio emission. Two major competing explanations of the radio spectrum of the Sagittarius $ {\rm A^*}$ are based on two different models, that is, hot accretion disk and jet. Hence, independent observational constraints are required to resolve related issues. We explore the possibility of using the observational data of the star S2, currently closest to the Galactic center, to distinguish physical models for the radio emission of the Sagittarius $ {\rm A^*}$, by applying the stellar cooling model to the Sagittarius $ {\rm A^*}$ with the orbital parameters derived from the observation. The relative difference in the electron temperature due to the stellar cooling from the star S2 at the pericenter is a few parts of a thousand and the consequent relative radio luminosity difference is order of $10^{-4}$. Therefore, one could possibly expect to observe the radio flux variation with a periodic or quasi-periodic modulation in the frequency range at $\nu \la {\rm 100 MHz}$ if the radiatively inefficient hot accretion flows are indeed responsible for the radio emission, contrary to the case of a jet. According to our findings, even though no periodic radio flux variations have been reported up to date a radiatively inefficient hot accretion disk model cannot be conclusively ruled out.

Planetary companions to the source stars of a caustic-crossing binary microlensing events can be detected via the deviation from the parent light curves created when the caustic magnifies the star light reflecting off the atmosphere or surface of the planets. The magnitude of the deviation is delta_p e_p rho_p^{-1/2}, where e_p is the fraction of starlight reflected by the planet and rho_p is the angular radius of the planet in units of angular Einstein ring radius. Due to the extraordinarily high resolution achieved during the caustic crossing, the detailed shapes of these perturbations are sensitive to fine structures on and around the planets. We consider the signatures of rings, satellites, and atmospheric features on caustic-crossing microlensing light curves. We find that, for reasonable assumptions, rings produce deviations of order 10% delta_p, whereas satellites, spots, and zonal bands produce deviations of order 1% delta_p. We consider the detectability of these features using current and future telescopes, and find that, with very large apertures (>30m), ring systems may be detectable, whereas spots, satellites, and zonal bands will generally be difficult to detect. We also present a short discussion of the stability of rings around close-in planets, noting that rings are likely to be lost to Poynting-Robertson drag on a timescale of order 10^5 years, unless they are composed of large (>>1 cm) particles, or are stabilized by satellites.

Recently, from the Hubble Space Telescope (HST) images of one of the Large Magellanic Cloud (LMC) events taken 6.3 years after the original lensing measurement, Alcock et al. were able to directly image the lens. Although the first resolved lens was identified for an LMC event, much more numerous lenses are expected to be resolved for Galactic bulge events. In this paper, we estimate the fraction of Galactic bulge events whose lenses can be directly imaged under the assumption that all bulge events are caused by normal stars. For this determination, we compute the distribution of lens proper motions of the currently detected Galactic bulge events based on standard models of the geometrical and kinematical distributions of lenses and their mass function. We then apply realistic criteria for lens resolution, and the result is presented as a function of the time elapsed after an original lensing measurement, $\Delta t$. If followup observations are performed by using an instrument with a resolving power of \theta_{PSF}=0.1'', which corresponds to that of HST equipped with the new Advanced Camera for Surveys, we estimate that lenses can be resolved for ~ 3% and 22% of disk-bulge events and for ~0.3% and 6% of bulge self-lensing events after $\Delta t=10$ and 20 years, respectively. The fraction increases substantially with the increase of the resolving power. If the instrument has a resolution of \theta_{PSF}=0.05'', which can be achieved by the {\it Next Generation Space Telescope}, we estimate that lenses can be resolved for $\sim 22%$ and 45% of disk-bulge events and for ~6% and 23% of bulge self-lensing events after \Delta t=10 and 20 years, respectively.

Optically thin and geometrically thick accretion flows are known to be responsible for the observed radio/X-ray luminosity relation of the X-ray bright galactic nuclei. It has also been suggested that supermassive black hole masses can be estimated from measurements of the core radio luminosity and the X-ray luminosity by using the advection-dominated accretion flow (ADAF) model. In this study we increase the number of data available by compiling the radio/X-ray fluxes and the mass in published literatures, and compare the observed ratio of the luminosities with predictions from various models of optically thin accretion flows. Semi-analytically derived relations of the luminosities are presented in cases of the standard ADAF model and modified ADAF models, in which a truncation of inner parts of the flows and winds causing a reduction of the infalling matter are included. We show that the observed relation can be used indeed to estimate the supermassive black hole mass, provided that properties of such accretion flows are known. Having investigated sensitivities of the method on modifications of the 'standard' ADAF model, we find that a general trend of model predictions from the 'standard' ADAF, the truncated ADAF and the 'windy' ADAF are somewhat indistinguishable. We also find, however, that the extreme case of the windy model is inconsistent with currently available observational data, unless microphysics parameters are to be substantially changed. High resolution radio observations, however, are required to avoid the contamination of non-disk components, such as, a jet component, which, otherwise, results in the over-estimated SMBH mass.

The principal aim of observational helioseismology is to determine mode parameters of the solar oscillations as accurately as possible. Yet estimates of the mode parameters are subject to many sources of noise, including inevitable gaps in data, stochastic nature of excitation processes, interferences among modes. Another uncertainty in frequency determination results from a modulation of the oscillation frequency in itself. It is well known that variation in the mean strength of the solar magnetic field modulates the frequency. We consider effects of the solar magnetic field variation on the spectral behavior of the power spectrum using a simple model. We show that the solar magnetic field variation can cause unwanted sidelobes beside the main peak in the power spectrum, and consequently show that the effect of the frequency modulation due to the solar magnetic field variation can lead a bias in frequency determinations. This effect should be considered seriously particularly when the $ l \ne 0$ $p$-mode multiplets are analyzed to measure the rotational splitting since the separations could be comparable. In addition to a bias in frequency estimates, the line width is also likely to be overestimated due to a line broadening. We therefore suggest that determination of the mode parameters should be derived with due care considering the magnetic field variation when a long data set is analyzed. The frequency modulation should be taken into account in the analysis to make a solid conclusion on any subtle dependence of the mode parameters, such as, that on the solar cycle. We conclude by pointing out that a new method is required to accommodate the stochastic force and the phase variation.

We have studied the long-term X-ray light curve (2$-$10 keV) of the luminous Seyfert 1 galaxy MCG-2-58-22 by compiling data, from various X-ray satellites, which together cover more than 20 years. We have found two distinct types of time variations in the light curve. One is a gradual and secular decrease of the X-ray flux, and the other is the episodic increase of X-ray flux (or flare) by a factor of 2$-$4 compared with the level expected from the secular variation. We detected 3 such flares in total; a representative duration for the flares is $\sim $2 years, with intervening quiescent intervals lasting $\sim 6-8$ years. We discuss a few possible origins for these variabilities. Though a standard disk instability theory may explain the displayed time variability in the X-ray light curve, the subsequent accretions of stellar debris, from a tidal disruption event caused by a supermassive black hole in MCG-2-58-22, cannot be ruled out as an alternative explanation.

We investigate the pattern of anomalies in the light curves of caustic-crossing binary microlensing events induced by spot(s) on the lensed source star. For this purpose, we perform simulations of events with various models of spots. From these simulations, we find that the spot-induced anomalies take various forms depending on the physical state of spots, which is characterized by the surface brightness contrast, the size, the number, the umbra/penumbra structure, the shape, and the orientation with respect to the sweeping caustic. We also examine the feasibility of distinguishing the two possibly degenerate types of anomalies caused by a spot and a transiting planet and find that the degeneracy in many cases can be broken from the characteristic multiple deviation feature in the spot-induced anomaly pattern caused by the multiplicity of spots.

We investigate consequences of a continuously energy-injecting central engine of gamma-ray burst (GRB) afterglow emission, assuming that a highly magnetized pulsar is left beaming in the core of a GRB progenitor. Beaming and continuous energy-injection are natural consequences of the pulsar origin of GRB afterglows. Whereas previous studies have considered continuous energy-injection from a new-born pulsar to interpret the deviation of afterglow light curves of GRBs from those with the simple power law behavior, a beaming effect, which is one of the most important aspects of pulsar emissions, is ignored in earlier investigations. We explicitly include the beaming effect and consider a change of the beaming with time due to a dynamical evolution of a new-born pulsar. We show that the magnitude of the afterglow from this fireball indeed first decreases with time, subsequently rises, and declines again. One of the most peculiar optical afterglows light curve of GRB 970508 can be accounted for by continuous energy injection with beaming due to a highly magnetized new-born pulsar. We discuss implications on such observational evidence for a pulsar.

We propose a redshift estimator for the long ($T_{90} > 20$) gamma-ray bursts (GRBs) observed by the BATSE. It is based on an empirical relation between the redshift and the power-law index of power density spectra (PDSs) of the observed GRBs. This relation is constructed by using the fact that the power-law index is dependent upon a characteristic timescale of GRB light curves which are inevitably stretched by cosmological time dilation. We construct the empirical relation using both individual PDSs and averaged PDSs. An error estimates of $z$ are 1.09 and 1.11 for the empirical relation by individual PDS fits, 1.72 and 1.56 by averaged PDS fits, for the least squares fit and the maximum likelihood fit, respectively. We attempt to determine the spatial distribution of the GRBs observed by the BATSE as a function of redshifts on the basis of the resulting redshift estimator. We find that the obtained spatial distribution of the observed GRBs seems consistent with that of the GRBs whose redshifts are reported, even though the estimated errors are not very accurate. The GRBs observed by the BATSE seem distributed within $z \sim 5-6$. This result has implications on theoretical calculations of stellar formations at high redshifts and beaming geometry via a statistical study of the observed GRBs involving beaming-induced luminosity functions. We discuss such implications, and possible uncertainties of the suggested method.

In this paper, we show that the pattern of microlensing light curve anomalies induced by multiple planets are well described by the superposition of those of the single-planet systems where the individual planet-primary binary pairs act as independent lens systems. Since the outer deviation regions around the planetary caustics of the individual planets occur in general at different locations, we find that the pattern of anomalies in these regions are hardly affected by the existence of other planet(s). This implies that even if an event is caused by a multiple planetary system, a simple single-planet lensing model is good enough for the description of most anomalies caused by the source passage of the outer deviation regions. Detection of the anomalies resulting from the source trajectory passing both the outer deviation regions caused by more than two planets will provide a new channel of detecting multiple planets.

We study the power density spectrum (PDS) of light curves of the observed gamma-ray bursts (GRBs) to search for a direct signature for cosmological time dilation in the PDS statistics with the GRBs whose redshifts $z$'s are known. The anticorrelation of a timescale measure and a brightness measure is indirect evidence of its effect. On the other hand, we directly demonstrate that a time dilation effect can be seen in GRB light curves. We find that timescales tend to be shorter in bursts with small redshift, as expected from cosmological time-dilation effects, and we also find that there may be non-cosmological effects constituting to this correlation. We discuss its implication on interpretations of the PDS analysis results. We put forward another caution to this kind of analysis when we statistically exercise with GRBs whose $z$ is unknown.

Gaudi & Gould (1997) showed that close companions of remote binary systems can be efficiently detected by using gravitational microlensing via the deviations in the lensing light curves induced by the existence of the lens companions. In this paper, we introduce another channel to detect faint close-in binary companions by using microlensing. This method utilizes a caustic-crossing binary lens event with a source also composed of binary stars, where the companion is a faint star. Detection of the companion is possible because the flux of the companion can be highly amplified when it crosses the lens caustic. The detection is facilitated since the companion is more amplified than the primary because it, in general, has a smaller size than the primary, and thus experiences less finite source effect. The method is extension of the previous one suggested to detect close-in giant planets by Graff & Gaudi (2000) and Lewis & Ibata (2000) and further developed by Ashton & Lewis (2001). From the simulations of realistic Galactic bulge events, we find that companions of K-type main sequence or brighter can be efficiently detected from the current type microlensing followup observations by using the proposed method. We also find that compared to the method of detecting lens companions for which the efficiency drops significantly for binaries with separations $\lesssim 0.2$ of the angular Einstein ring radius, the proposed method has an important advantage of being able to detect companions with substantially smaller separations down to order of 0.01 \theta_E.

We study the power density spectrum (PDS) of artificial light curves of observed gamma-ray bursts (GRBs). We investigate statistical properties of GRB light curves by comparing the reported characteristics in the PDSs of the observed GRBs with those that we model, and discuss implications on interpretations of the PDS analysis results. Results of PDS analysis of observed GRBs suggest that the averaged PDS of GRBs follows a power law over about two decades of frequency with the power law index, -5/3, and the distribution of individual power follows an exponential distribution. Though an attempt to identify the most sensitive physical parameter has been made on the basis of the internal shock model, we demonstrate that conclusions of this kind of approach should be derived with due care. It is indicative that the physical information extracted from the slope can be misleading. We show that the reported slope and the distribution can be reproduced by adjusting the sampling interval in the time domain for a given decaying timescale of individual pulse in a specific form of GRB light curves. In particular, given that the temporal feature is modeled by a two-sided exponential function, the power law behavior with the index of -5/3 and the exponential distribution of the observed PDS is recovered at the 64 ms trigger time scale when the decaying timescale of individual pulse is $\sim 1$ second, provided that the pulse sharply rises. Another way of using the PDS analysis is an application of the same method to individual long bursts in order to examine a possible evolution of the decaying timescale in a single burst.

We study the effects of the beaming-induced luminosity function on statistics of observed GRBs, assuming the cosmological scenario. We select and divide the BATSE 4B data into 588 long bursts (T$_{90}>2.5$ sec) and 149 short bursts (T$_{90}<2.5$ sec), and compare the statistics calculated in each subgroup. The $<V/V_{\rm max}>$ of the long bursts is $ 0.2901\pm 0.0113$, and that of the short bursts is $0.4178\pm 0.0239$, which is a Euclidean value. For luminosity function models, we consider a cylindrical-beam and a conic-beam. We take into account the spatial distribution of GRB sources as well. A broad luminosity function is naturally produced when one introduces beaming of GRBs. We calculate the maximum detectable redshift of GRBs, $z_{\rm max}$. The estimated $z_{\rm max}$ for the cylindrical-beam case is as high as $\sim 14$ for the long bursts and $\sim 3$ for the short bursts. The large $z_{\rm max}$ value for the short bursts is rather surprising in that the $<V/V_{\rm max}>$ for this subgroup is close to the so-called Euclidean value, 0.5. We calculate the fraction of bursts whose redshifts are larger than a certain redshift $z'$, i.e. $f_{\rm > z'}$. When we take $z'=3.42$ and apply the luminosity function derived for the cylindrical-beam, the expected $f_{\rm > z'}$ is $\sim 75 %$ for long bursts. When we increase the opening angle of the conic beam to $\Delta \theta =3^\circ.0$, $f_{\rm > z'}$ decreases to $\sim 20 %$ at $ {\rm z'=3.42}$. We conclude that the beaming-induced luminosity functions are compatible with the redshift distribution of observed GRBs and that the apparent Euclidean value of $<V/V_{\rm max}>$ may not be due to the Euclidean space distribution but to the luminosity distribution.

If a supermassive black hole resides in the centers of galaxies, there are several effects expected to be observed. It is likely that accretion disks are around the supermassive black hole. Stellar interactions with the accretion flows around the supermassive black hole play a role in that an flying-by star may cool a hot accretion disk as a result of Comptonization. It is shown that the Comptonization of the stellar emission will take place in a hot accretion disk such as the ADAFs around the supermassive black hole and become a relatively important source of the accretion disk cooling when the mass accretion rate is low, and stars are passing outer parts of the hot accretion disk. We suggest that such a stellar cooling can be observed in the radio frequency regime since synchrotron luminosity of the ADAFs depends strongly on the electron temperature and occurs much more frequently than a tidal disruption event.

It has been suggested that optically thin and geometrically thick accretion flows are responsible for the observed radio/X-ray luminosity relation of the X-ray bright galactic nuclei. If this is the case then central supermassive black hole masses can be estimated directly from measurements of the core radio luminosity and the X-ray luminosity, provided that properties of such accretion flows are known. Calculated ratios of the luminosities are presented in cases of the standard ADAF model and modified ADAF models, in which a truncation of inner parts of the flows and winds causing a reduction of the infalling matter are included. We compare the observed ratio of the luminosities with predictions from models of optically thin accretion flows. We also discuss the possible effects of the convection in ADAFs. We confirm that the supermassive black hole (SMBH) mass estimate is possible with the radio/X-ray luminosity relation due to ADAF models in the absence of a radio jet. We find that observational data are insufficient to distinguish the standard ADAF model from its modified models. However, the ADAF model with convection is inconsistent with observations, unless microphysics parameters are to be substantially changed. High resolution radio observations are required to avoid the contamination of other components, such as, a jet component. Otherwise, the SMBH mass is inclined to be over-estimated.

We study the statistics of gamma-ray bursts, assuming that gamma-ray bursts are cosmological and they are beamed in the form of a conical jet with a large bulk Lorentz factor $\sim 100$. In such a conic beam, the relativistic ejecta may have a spatial variation in the bulk Lorentz factor and the density distribution of gamma-ray emitting jet material. An apparent luminosity function arises because the axis of the cone is randomly oriented with respect to the observer's line of sight. The width and the shape of the luminosity function are determined by the ratio of the beam opening angle of the conical jet to the inverse of the bulk Lorentz factor, when the bulk Lorentz factor and the jet material density is uniform on the photon emitting jet surface. We calculate effects of spatial variation of the Lorentz factor and the spatial density fluctuations within the cone on the luminosity function and the statistics of gamma-ray bursts. In particular, we focus on the redshift distribution of the observed gamma-ray bursts. The maximum distance to and the average redshift of the gamma-ray bursts are strongly affected by the beaming-induced luminosity function. The bursts with the angle-dependent Lorentz factor which peaks at the center of the cone have substantially higher average gamma-ray burst redshifts. When both the jet material density and the Lorentz factor are inhomogeneous in the conical beam, the average redshift of the bursts could be 5 times higher than that of the case in which relativistic jet is completely homogeneous and structureless. Even the simplest models for the gamma-ray burst jets and their apparent luminosity distributions have a significant effect on the redshift distribution of the gamma-ray bursts.

PSR B0656+14 was observed by the {\it Rossi X-ray Timing Explorer (RXTE)} with the proportional counter array (PCA) and the high-energy x-ray timing experiment (HEXTE) for 160 ksec during August 22 -- September 3, 1997. No pulsation was firmly found in the timing analysis, in which the contemporaneous radio ephemeris and various statistical tests were applied for searching evidence of pulsation. A marginal detection of pulsation at a confidence level of 95.5% based on the $H$-test was found with data in the whole HEXTE energy band. In the energy band of 2-10 keV the {\it RXTE} PCA upper limits are about one order of magnitude lower than that from {\it ASCA} GIS data. If the {\it CGRO} EGRET detection of this pulsar is real, considering the common trait that most EGRET-detected pulsars have a cooling spectrum in hard x-ray and gamma ray energy bands, the estimated {\it RXTE} upper limits indicate a deviation (low-energy turn-over) from a cooling spectrum starting from 20 keV or higher. It in turn suggests an outer-magnetospheric synchrotron-radiation origin for high-energy emissions from PSR B0656+14. The {\it RXTE} PCA upper limits also suggest that a reported power-law component based on {\it ASCA} SIS data in 1-10 keV fitted jointly with {\it ROSAT} data, if real, should be mainly unpulsed.

We report results of RXTE observations of PSR B1951+32 using the PCA instrument for 19k seconds during 1996 March 24th. We applied the contemporaneous radio ephemeris and various statistical tests to search for evidence of pulsation. These analyses yield intriguing yet inconclusive evidence for the presence of the pulsation in the time series: confidence level for the presence of pulsation is 94.5% in the 2.0-4.8 keV band and 97.6% in the 4.8-6.3 keV band based on the H-test. Under the premise of non-detection of pulsation, we derive estimated 2-sigma upper limits for the pulsed flux to be $3.3\times 10^{-6} cm^{-2} s^{-1} keV^{-1}$ in the 2.0-4.8 keV band, $2.7\times 10^{-6} cm^{-2} s^{-1} keV^{-1}$ in the 4.8-8.5 keV band, and $2.0\times 10^{-6} cm^{-2} s^{-1} keV^{-1}$ in the 8.5-13.0 keV band. These upper limits are consistent with the trend of spectral turn-over from high-energy gamma-ray emission as suggested by the OSSE upper limit. Such turn-over strongly suggests the outer magnetosphere as the emission site for pulsed gamma-rays. These RXTE upper limits for X-ray pulsation are, on the other hand, not consistent with the extrapolation of reported power-law spectra from the point source observed by ROSAT in the 0.1-2.4 keV band, assuming a constant pulse fraction: The pulsed soft X-ray emission detected by ROSAT must follow a much softer spectrum than that of the overall point source.