The Gravitational wave (GW) has opened a new window to the universe beyond the electromagnetic spectrum. Since 2015, dozens of GW events have been caught by ground based GW detectors through laser interferometry. However, both the ground based and future space-based detectors are two-dimensional, which have very limited directional response to the GW. Here we propose a simple three-dimensional (3-D) laser interferometer in the shape of a regular triangular pyramid. Such a 3-D detector can provide much stronger direction constraints for GW sources and wider directional response than the two-dimensional ones, which is more efficient in the research jointing GW and electromagnetic emission. The most sensitive band of such a detector is around kilo-Hz, which is suitable for measuring the post-merger signal from the binary neutron stars. Such a detection will be of importance in the restriction of the equation of state of the neutron stars.

Since the discovery of solar cycle related with magnetic field in 1908, deep seated oscillatory dynamo has been studied extensively. However, there are still open questions on the solar dynamo, e.g., asymmetric conversion between large-scale poloidal and toroidal field as well as physics underlying the butterfly pattern of sunspots. Here we report a new generation of large-scale magnetic field and process of energy release. The inductive action of fluid motions pervading the solar interior is represented by a RLC circuit in which the toroidal field built up through twisting of poloidal field, so called $\omega$-effect, plays the role of a capacitor. Such a RLC circuit not only provides a self-sustained oscillatory system avoiding Cowling's antidynamo theorem, but also site of rapid magnetic reconnection which reproduces quadrupole magnetic field interpreting the behavior of sunspots and moving of foot-point in solar activities. Moreover, parameters of the circuit and the Sun are well consistent with the 22-year solar cycle.

Weak Equivalent Principle (WEP) can be tested through the parameterized post-Newtonian parameter $\gamma$, representing the space curvature produced by unit rest mass. The parameter $\gamma$ in turn has been constrained by comparing the arrival times of photons originating in distant transient events, such as gamma-ray bursts, fast radio bursts as well as giant pulses of pulsars. Those measurements normally correspond to an individual burst event with very limited energy bands and signal-to-noise ratio (S/N). In this letter, the discrepancy in the pulse arrival times of the Crab Pulsar between different energy bands is obtained by the phase difference between corresponding pulse profiles. This allows us to compare the pulse arrival times at the largest energy band differences, between radio and optical, radio and X-ray, radio and gamma-ray respectively. As the pulse profiles are generated by phase-folding thousands of individual pulses, the time discrepancies between two energy bands are actually measured from thousands of events at each energy band, which corresponds to much higher S/N. The upper limit of the $\gamma$ discrepancy set by such an extensively-observed and well-modeled source is as follows: $\gamma_{radio}-\gamma_{\gamma-ray}<3.28\times10^{-9}$ at the energy difference of $E_{\gamma-ray}/E_{radio}\sim10^{13}$, $\gamma_{radio}-\gamma_{X-ray}<4.01\times10^{-9}$ at the energy difference of $E_{X-ray}/E_{radio}\sim10^{9}$, $\gamma_{radio}-\gamma_{optical}<2.63\times10^{-9}$ at $E_{optical}/E_{radio}\sim10^{5}$, and $\gamma_{optical}-\gamma_{\gamma-ray}<3.03\times10^{-10}$ at $E_{\gamma-ray}/E_{optical}\sim10^8$. This actually measures the arrival times of freely-falling photons in the gravitational field of the Milky Way with the largest amount of events and with data of the highest S/N, which tests WEP at the energy band differences that has never been reached before.

The huge derivative of orbital period observed in binary pulsar PSR B1259-63, the torque reversal displaying on low mass X-ray binary, 4U1627-67 and the long term change of orbital period of PSR J2051-0827, seem totally unrelated phenomena occurring at totally different pulsar systems. In this paper, they are simply interpreted by the same mechanism, residual Doppler shift. In a binary system with periodic signals sending to an observer, the drift of the signal frequency actually changes with the varying orbital velocity, projected to line of sight at different phases of orbit. And it has been taken for granted that the net red-shift and blue-shift of an full orbit circle be cancelled out, so that the effect of Doppler shift to the signal in binary motion cannot be accumulated over the orbital period. However, taking the propagation time at each velocity state into account, the symmetry of the velocity distribution over the orbital phase is broken. Consequently, the net Doppler shift left in an orbit is non-zero. Understanding this Newtonian second Doppler effect not only makes pulsars better laboratory in the test of gravitational effects, but also allows us to extract the angular momentum of the pulsar of PSR J2051-0827, $\leq 2\times 10^{43}gcm^2$; and the accretion disc of 4U 1627-67, $7\times 10^{50} gcm^2/s$, respectively, which are of importance in the study of structure of neutron stars and the physics of accretion disc of X-ray binaries.

Twin kilo-Hertz Quasi-periodic oscillation with ratio 3/2 has been found in some compact sources, which is believed to be related with their innermost regions of accretion disks, and hence carrying information of gravity in strong regime. However, more complicated phenomena have been revealed, e.g., the twin kilo-Hertz Quasi-periodic oscillation of 4U 1820-20 start increasing at certain spectral state and then saturate upon reaching certain level. Moreover, such quasi-periodic oscillation is not uniform and has multiple peaks, which displays random feature. This letter suggests that these challenges could be easy to understand if the quasi-periodic oscillation originates in jets of compact objects. With a seed periodicity originating in either neutron star spin or accretion, shell collision develops in jet. The corresponding twin kilo-Hertz quasi-periodic oscillations automatically carry random feature, vary in frequency; and couple with luminosity, spectral properties, which well account for the observation of 4U 1820-30. New quasi-periodic oscillation of 4U 1820-30 is predicted, which can test the validity of this model. And the scenario is applicable to other compact object like AGN.

The quasi-periodic feature of 1-10 years exhibited in pulsar timing noise has not been well understood since 1980. The recently demonstrated correlation between timing noise and variation of pulse profile motivates us to further investigate its origins. We suggest that the quasi-periodicity feature of timing noise, with rapid oscillations lying on lower frequency structure, comes from the geodetic precession of an unseen binary system, which induces additional motion of the pulsar spin axis. The resultant change of azimuth and latitude at which the observer's line of sight crosses the emission beam is responsible for the variation of timing noise and pulse profile respectively. The first numerical simulation to both timing noise and pulse profile variation are thus performed, from which the orbital periods of these pulsars are of 1-35 minutes. Considering the existence of the ultracompact binary white dwarf of orbital period of 5.4 minute, HM Cancri, such orbital periods to pulsar binaries are not strange. The change of latitude of the magnetic moment exceeding the range of emission beam of a pulsar results in the intermittency, which explains the behavior of PSR B1931+24. Therefore, it provides not only a mechanism of the quasi-periodic feature displayed on some ``singular" pulsars, like timing noise, variation of pulse profile and intermittency; but also a new approach of searching ultra-compact binaries, possibly pulsar-black hole binary systems.

The sprinkler pivots on a bearing on top of its threaded attachment nut. It is driven in a circular motion by a spring-loaded arm pushed back by the water stream which returns to "impact" the stream. The water stream can thus rotate around a fix axis. Analogously in our universe, the outflow or jet formed by the relativistic plasma corresponds to the water stream of a sprinkler; and the baryons in the tilted outer accretion disc or torus play the role of "impact arm". Then the jet aligning with the inner parts of a warped disc can directly "touch" the outer region of the disc. The resultant collision between such rapid leptons and slow baryons automatically accounts for the main features of broad-line region of active galactic nuclei. Moreover, it naturally provides a channel of dissipating the angular momentum of an accretion disc, which has long been a problem in theory of accretion disc.

The model of non-ballistic jet motion proposed in 2008 provides a simple explanation to the inward jet motion and bent jet. Recently, evidences of such a non-radial motion increase rapidly, and more complicated morphologies appear. On the other hand, the ballistic plus precession model likely holds in majority samples of jet motion. This paper discusses the relationship between the ballistic and non-ballistic model of jet motion, which suggests that the interaction of ejectors with ambient matter can produce knots at different stages of evolution and hence different separations to the core. And as a jet precesses, knots produced between the core and the deceleration radius result in spiral pattern expected by the model of ballistic plus precession; and knots generated at the deceleration radius display non-radial motion such as bent jet or oscillation of ridge-line. This paper develops the first non-ballistic model in four aspects. Firstly, it provides a numerical simulation to the production of multi-knot for a precessing jet. Secondly, it fits the precession behavior of multi-knot and interprets the oscillation of ridge lines like S5 1803+784. Thirdly, it gives an unified interpretation to the bent jet applicable to both multi-knot and single knot. And fourthly, the problem of very large numbers of observed outward motions as opposed to the inward ones is addressed in a new scope.

Short precession periods like 164d of SS433 can be well determined by observations of time scales longer or much longer than the precession period. However, it doesn't work for sources with precession periods of millions of years. This paper utilizes the particular morphologies of X-shaped sources, so that the 3 dimension kinematics of lobes can be obtained. Thus, for the first time, the million-year precession period of X-shaped sources by observer on the Earth can be determined elegantly: $(6.1\pm 1.5)$Myr, $(1.8\pm 0.5)$Myr, and $(3.2\pm 1.2)$Myr for 3C52, 3C223.1 and 4C12.03 respectively. The result naturally explains the asymmetry displayed in the morphology of these sources, and the effect of propagation time on the diversity of morphologies is well demonstrated. The precession period may originate from long-term effects of a binary super-massive black hole system, which is a potential source of gravitational wave radiation.

Context. We present evidence that parsec-scale jets in BL Lac objects may be significantly distinct in kinematics from their counterparts in quasars. We argued this previously for the BL lac sources 1803+784 and 0716+714, report here a similar pattern for another well-known BL Lac object, PKS 0735+178, whose nuclear jet is found to exhibit kinematics atypical of quasars. Aims. A detailed study of the jet components' motion reveals that the standard AGN paradigm of apparent superluminal motion does not always describe the kinematics in BL Lac objects. We study 0735+178 here to augment and improve the understanding of the peculiar motions in the jets of BL Lac objects as a class. Methods. We analyzed 15 GHz VLBA (Very Long Baseline Array) observations (2cm/MOJAVE survey) performed at 23 epochs between 1995.27 and 2008.91. Results. We found a drastic structural mode change in the VLBI jet of 0735+178, between 2000.4 and 2001.8 when its twice sharply bent trajectory turned into a linear shape.We further found that this jet had undergone a similar transition sometime between December 1981 and June 1983. A mode change, occurring in the reverse direction (between mid-1992 and mid-1995) has already been reported in the literature. These structural mode changes are found to be reflected in changed kinematical behavior of the nuclear jet, manifested as an apparent superluminal motion and stationarity of the radio knots. In addition, we found the individual mode changes to correlate in time with the maxima in the optical light curve. The last two transitions occurred before a (modest) radio flare. The behavior of this pc-scale jet appears to favor a scenario involving non-ballistic motions of the radio knots, produced by the precession of a continuous jet within the ambient medium.

Faster-than-light or superluminal motion was originally predicted as a relativistic illusion of ballistic moving ejecta, and confirmed in a few tens of sources observationally. However, the recent results of the long-term multi-epoch observations of quasars, active galaxies, tracing the structure further along the jets and following the motion of individual features for longer time, rise questions that are difficult to understand by the standard ballistic model. I.e., the ejecta are aligned with the local jet direction, instead of the core; and within individual jets apparently inward-moving features are observed. Here we show that these unexpected phenomena, although only a small fraction among large samples, indicate the existence of non-ballistic jet motion, in which a continuous jet produces a discrete hot spot. And the precession of such a hot spot in the plane of the sky appears superluminal. Therefore, an unified and simple interpretation to the new results is obtained, which can be further tested through its predictions on the evolution of ejecta. The study is of importance in the understanding of the nature of superluminal motion, the interaction of jets and surrounding materials, as well as the common physics underlying quasars and microquasars.

An extreme Kerr black hole (BH) surrounded by a precessing disk is invoked to explain the light curves of gamma-ray bursts (GRBs) based on the coexistence of the Blandford-Znajek (BZ) and the magnetic coupling (MC) processes. The overall shape of the light curves and the duration of GRBs are interpreted by the evolution of the half-opening angle of the magnetic flux on the BH horizon, and the complex temporal structures are modulated by the precession and nutation of the jet powered by the BZ process. The time profile of the emission exhibits a fast rise and a slow decay due to the effect of the evolution of the half-opening angle. The light curves of several GRBs are well fitted by this model with only six free parameters.

Superluminal motion has been interpreted as relativistically moving out flow. The increasing observations of such sources provide chance to investigate their mechanism in detail. This paper proposes that superluminal motion may be jet precession induced blob motion. This model can interpret a number of observational phenomena which are not well understood under the scenario of bulk motion, such as the absence of the receding blob after the outburst and the receding blob becomes brighter than the approaching one in XTE J1550-564 . Moreover the compton scattering of photos from a warped accretion disk may form a structured jet. And the precession of the structured jet may lead to the complicated state transitions observed in microquasars like GRS 1915+105. Applying to AGNs, the model may explain the scales of length and time of the phenomena proportioning to the mass of the black hole. The new model can be tested easily on microquasars, i.e., XTE J1550-564, on which the previous receding blob should move towards the core instead of away from it after 2002, and the previous approaching blob shall disappear. %For GRS 1915+105 the vibration of the two blobs should be observed in both %the time scales of year and days.

The strange timing property of X-ray pulsar 1E 1207.4-5209 can be explained by the hypothesis that it is a member of an ultra-compact binary system. This paper confronts the ultra-compact assumption with the observed properties of this pulsar. The gravitational potential well of an ultra-compact binary can enlarge the corotation radius and thus make it possible for accreting material to reach the surface of the NS in the low accretion rate case. Thus the generation of the absorption features should be similar to the case of accreting pulsars. The close equality of the energy loss by fast cooling of the postsupernova neutron star and the energy dissipation needed for a wide binary evolving to an ultra-compact binary demonstrates that the ultra-compact binary may be formed in 10-100yr after the second supernova explosion. Moreover, the ultra-compact binary hypothesis can well explain the the absence of optical counterpart and the observed two black body emissions. We suggest a simple method which can test the binary nature directly with XMM-Newton and Chandra observations. We further predict that the temperature of the two black bodies should vary at different pulse periods.

There are four puzzles on 1E 1207.4-5209: (1) the characteristic age of the pulsar is much higher than the estimated age of the supernova remnant; (2) the magnetic field inferred from spin-down is significantly different from the value obtained from the cyclotron absorption lines; (3) the spinning down of the pulsar is non-monotonic; (4) the magnitude of the frequency's first derivative varies significantly and its sign is also variable. The third puzzle can be explained by a wide binary system, with orbital period from 0.2 to 6 yr. This letter proposes that all four puzzles can be explained naturally by an ultra-compact binary with orbital period of between 0.5 and 3.3 min. With the shortest orbital period and a close distance of 2 kpc, the characteristic amplitude of gravitational waves is $h\sim3\times10^{-21}$. It would be an excellent source of gravitational-wave detectors such as the Laser Interferometer Space Antenna.

Unexpected sign, significant magnitude and time-variation of frequency second derivative exist not only in singular radio pulsars but also in Soft Gamma repeaters (SGRs) and Anomalous X-ray pulsars (AXPs). This paper shows that the these phenomena are related, and can be interpreted by a simple unified model, long-term orbital effect. Thus many of previous ``isolated'' pulsars may be binary pulsars, i.e., orbital period $P_b\approx(47, 72)$min for AXP 1E 2259+586, and $P_b\approx (20, 34)$min for PSR J1614-5047, $P_b\approx (3.6, 6.4)$min for SGR 1900+14, and $P_b\approx (1.5, 5.8)$min for SGR 1806-20. In addition to X-ray pulsar, 1E 1207.4-5209, these two SGRs may be the new ultra-compact binary pulsars with orbital period of a few minutes. There might be more binary pulsars that suitable for gravitational wave detection than we had thought.

There are two kinds of expressions on the precession of orbital plane of a binary pulsar system, which are given by Barker $&$ O'Connell (1975) and Apostolatos et al. (1994), Kidder (1995) respectively. This paper points out that these two kinds of orbital precession velocities are actually obtained by the same Lagrangian under different degrees of freedom. Correspondingly the former expression is not consistent with the conservation of the total angular momentum vector; whereas the latter one is. Damour $&$ Sch\"afer (1988) and Wex $&$ Kopeikin (1999) have applied Barker $&$ O'Connell's orbital precession velocity in pulsar timing measurement. This paper applies Apostolatos et al. $&$ Kidder's orbital precession velocity in pulsar timing measurement. We analyze that Damour $&$ Sch\"afer's treatment corresponds to negligible Spin-Orbit induced precession of periastron. Whereas the effects corresponding to Wex $&$ Kopeikin and this paper are both significant (however they are not equivalent). The observational data of two typical binary pulsars, PSR J2051-0827 and PSR J1713+0747 apparently support significant Spin-Orbit coupling effect. Further more, the discrepancies between Wex $&$ Kopeikin and this paper can be tested on specific binary pulsars with orbital plane nearly edge on. If the orbital period derivative of double-pulsar system PSRs J0737-3039 A and B, with orbital inclination angle $i=87.7_{-29}^{+17}$deg, is much larger than that of the gravitational radiation induced one, then the expression of this paper is supported, otherwise Wex $&$ Kopeikin's expression is supported.

PSR B1937+21 was the first millisecond pulsar ever measured. Which has been appeared as a singular pulsar. The high-precision observation of this pulsar shows systemic long-term variation in the residual of time of arrivals. This paper modelled the secular variation by the orbital precession induced time delay of a binary pulsar system. The fitting requires that as a binary pulsar, PSR B1937+21 should have small companion star, $m_2\sim 10^{-2}M_{\odot}$, and projected semi-major axis, $x\sim 10^{-4}$s to $10^{-3}$s. Which corresponds to ignorable radial velocity of the pulsar to the line of sight. This might explain why it has been measured as a singular pulsar instead of a binary pulsar.

The measurability of the spin--orbit (S--L) coupling induced orbital effect is dependent on the orientation and magnitude of the orbital precession velocity, ${\bf \Omega}_0$. This paper derives ${\bf \Omega}_0$ in the case that both spins in the binary system contribute to the spin--orbit (S--L) coupling, which is suitable for the most popular binary pulsars, Neutron star--White Dwarf star (NS--WD) binaries (as well as for NS--NS binaries). This paper shows that from two constraints, the conservation of the total angular momentum and the triangle formed by the orbital angular momentum, ${\bf L}$, the sum the spin angular momenta of the two stars, ${\bf S}$, and the total angular momentum, ${\bf J}$, the orbital precession velocity, ${\bf \Omega}_0$, along ${\bf J}$ is inevitable. Moreover, by the relation, $S/L\ll 1$, which is satisfied for a general binary pulsar, a significant ${\bf \Omega}_0$ (in magnitude) is inevitable, 1.5 Post Newtonian order (PN). Which are similar to the case of one spin as discussed by many authors. However unlike the one spin case, the magnitude of the precession velocity of ${\bf \Omega}_0$ varies significantly due to the variation of the sum the spin angular momenta of the two stars, ${\bf S}$, which can lead to significant secular variabilities in binary pulsars.

Variations of dispersion measure (DM) have been observed in some binary pulsars, which can not be well explained by the propagation effects, such as turbulence of the interstellar media (ISM) between the Earth and the pulsar. This paper provides an alternative interpretation of the phenomena, the geodetic precession of the orbit plane of a binary pulsar system. The dynamic model can naturally avoid the difficulties of propagation explanations. Moreover the additional time delay represented by the DM variation of two binary pulsars can be fitted numerically, through which some interesting parameters of the binary pulsar system, i.e., the moment of inertia of pulsars can be obtained, $I_1=(2.0\pm0.6)\times 10^{45}$g cm$^{2}$. The elimination of the additional time delay by the dynamic effect means that ISM between the these pulsars and the Earth might also be stable, like some other binary pulsars.

More and more binary pulsars show significant secular variations, in which the measured projected semi-major axis, $\dot{x}^{obs}$, and the first derivative of orbital period, $\dot{P}_{b}^{obs}$, are several order of magnitude larger than the prediction of general relativity (GR). This paper shows that the geodetic precession induced orbital effects can explain both $\dot{x}$ and $\dot{P}_{b}$ measured in binary pulsars. Moreover, by this model we can automatically estimate the magnitude of the spin angular momenta of the pulsar and its companion star, and therefore the moment of inertia ($10^{44}$g cm$^2$ to $10^{45}$g cm$^2$) of pulsar of binary pulsar systems, which agrees well with theoretical predictions. In other words, the contamination (residual represented by $\dot{x}$ and $\dot{P}_{b}$) in pulsar timing measurements might be caused by geodetic precession, an interesting gravitational effect we have been seeking for.