We investigate a simplified freeze-in dark-matter model in which the dark matter only interacts with the standard-model neutrinos via a light scalar. The extremely small coupling for the freeze-in mechanism is naturally realized in several neutrino-portal scenarios with the secret neutrino interactions. We study possible evolution history of the dark sector: the dark sector would undergo pure freeze-in production if the interactions between the dark-sector particles are negligible, while thermal equilibrium within the dark sector could occur if the reannihilation of the dark matter and the scalar mediator is rapid enough. We investigate the relic abundance in the freeze-in and dark freeze-out regimes, calculate evolution of the dark temperature, and study its phenomenological aspects on BBN and CMB constraints, the indirect-detection signature, as well as the potential to solve the small scale structure problem.

Fast radio bursts (FRBs) are short pulses observed in radio band from cosmological distances, some of which emit repeating bursts. The physical origins of these mysterious events have been subject to wide speculations and heated debates. One class of models invoke soft gamma-ray repeaters (SGRs), or magnetars, as the sources of FRBs. Magnetars are rotating neutron stars with extremely strong magnetic field and can sporadically emit bursts from X-ray ($\sim$keV) to soft gamma-ray ($\sim$sub-MeV) with duration from $10^{-2}$~s to $10^2$~s. However, even though some bright radio bursts have been observed from some magnetars, no FRB-like events had been detected to be associated with any magnetar burst, including one giant flare, and no radio burst has been associated with any X-ray event from any magnetar. Therefore, there is still no observational evidence for magnetar-FRB association up to today. Recently, a pair of FRB-like bursts (FRB~200428 hereafter) separated by 30 milliseconds (ms) were detected from the general direction of the Galactic magnetar SGR~J1935+2154. Here we report the detection of a non-thermal X-ray burst in the 1-250 keV energy band with the Insight-HXMT satellite, which we identify as emitted from SGR~J1935+2154. The burst showed two hard peaks with a separation of $\sim30$ ms, consistent with the separation between the two bursts in FRB~200428. The delay time between the double radio and X-ray peaks is $\sim8.57$ s, fully consistent with the dispersion delay of FRB~200428. We thus identify the non-thermal X-ray burst is associated with FRB~200428 whose high energy counterpart is the two hard peaks in X-ray. Our results suggest that the non-thermal X-ray burst and FRB~200428 share the same physical origin in an explosive event from SGR~J1935+2154.

We use the In data collected during the 2019 outburst from X-ray pulsar 4U 1901+03 to complement the orbital parameters reported by Fermi/GBM. Using the Insight-HXMT, we examine the correlation between the derivative of the intrinsic spin frequency and bolometric flux based on accretion torque models. It was found that the pulse profiles significantly evolve during the outburst. The existence of two types of the profile's pattern discovered in the Insight-HXMT data indicates that this source experienced transition between a super-critical and a sub-critical accretion regime during its 2019 outburst. Based on the evolution of the pulse profiles and the torque model, we derive the distance to 4U 1901+03 as 12.4+-0.2 kpc.

We study the evolution of the temporal properties of MAXI 1820+070 during the 2018 outburst in its hard state from MJD 58190 to 58289 with Insight-HXMT in a broad energy band 1-150 keV. We find different behaviors of the hardness ratio, the fractional rms and time lag before and after MJD 58257, suggesting a transition occurred at around this point. The observed time lags between the soft photons in the 1-5 keV band and the hard photons in higher energy bands, up to 150 keV, are frequency-dependent: the time lags in the low-frequency range, 2-10 mHz, are both soft and hard lags with a timescale of dozens of seconds but without a clear trend along the outburst; the time lags in the high-frequency range, 1-10 Hz, are only hard lags with a timescale of tens of milliseconds; first increase until around MJD 58257 and decrease after this date. The high-frequency time lags are significantly correlated to the photon index derived from the fit to the quasi-simultaneous NICER spectrum in the 1-10 keV band. This result is qualitatively consistent with a model in which the high-frequency time lags are produced by Comptonization in a jet.

We study the effect of compact extra dimensions on the gravitational wave luminosity and waveform. We consider a toy model, with a compactified fifth dimension, and matter confined on a brane. We work in the context of five dimensional ($5d$) general relativity, though we do make connections with the corresponding Kaluza-Klein effective $4d$ theory. We show that the luminosity of gravitational waves emitted in $5d$ gravity by a binary with the same characteristics (same masses and separation distance) as a $4d$ binary is 20.8\% less relative to the $4d$ case, to leading post-Newtonian order. The phase of the gravitational waveform differs by 26\% relative to the $4d$ case, to leading post-Newtonian order. Such a correction arises mainly due to the coupling between matter and dilaton field in the effective $4d$ picture and agrees with previous calculations when we set black holes' scalar charges to be those computed from the Kaluza-Klein reduction. The above correction is inconsistent with the recent gravitational-wave observations and it thus effectively rules out the possibility of such a simple compactified higher dimensions scenario. We also comment on how our results change if there are several compactified extra dimensions, and show that the discrepancy with $4d$ general relativity only increases.

Glitches correspond to sudden jumps of rotation frequency ($\nu$) and its derivative ($\dot{\nu}$) of pulsars, the origin of which remains not well understood yet, partly because the jump processes of most glitches are not well time-resolved. There are three large glitches of the Crab pulsar, detected in 1989, 1996 and 2017, which were found to have delayed spin-up processes before the normal recovery processes. Here we report two additional glitches of the Crab pulsar occurred in 2004 and 2011 for which we discovered delayed spin up processes, and present refined parameters of the largest glitch occurred in 2017. The initial rising time of the glitch is determined as $<0.48$ hour. We also carried out a statistical study of these five glitches with observed spin-up processes. The two glitches occurred in 2004 and 2011 have delayed spin-up time scales ($\tau_{1}$) of $1.7\pm0.8$\,days and $1.6\pm0.4$\,days, respectively. We find that the $\Delta{\nu}$ vs. $|\Delta{\dot\nu}|$ relation of these five glitches is similar to those with no detected delayed spin-up process, indicating that they are similar to the others in nature except that they have larger amplitudes. For these five glitches, the amplitudes of the delayed spin-up process ($|\Delta{\nu}_{\rm d1}|$) and recovery process ($\Delta{\nu}_{\rm d2}$), their time scales ($\tau_{1}$, $\tau_{2}$), and permanent changes in spin frequency ($\Delta{\nu}_{\rm p}$) and total frequency step ($\Delta{\nu}_{\rm g}$) have positive correlations. From these correlations, we suggest that the delayed spin-up processes are common for all glitches, but are too short and thus difficult to be detected for most glitches.

No robust detection of prompt electromagnetic counterparts to fast radio bursts (FRBs) has yet been obtained, in spite of several multi-wavelength searches carried out so far. Specifically, X/gamma-ray counterparts are predicted by some models. We planned on searching for prompt gamma-ray counterparts in the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) data, taking advantage of the unique combination of large effective area in the keV-MeV energy range and of sub-ms time resolution. We selected 39 FRBs that were promptly visible from the High-Energy (HE) instrument aboard Insight-HXMT. After calculating the expected arrival times at the location of the spacecraft, we searched for a significant excess in both individual and cumulative time profiles over a wide range of time resolutions, from several seconds down to sub-ms scales. Using the dispersion measures in excess of the Galactic terms, we estimated the upper limits on the redshifts. No convincing signal was found and for each FRB we constrained the gamma-ray isotropic-equivalent luminosity and the released energy as a function of emission timescale. For the nearest FRB source, the periodic repeater FRB180916.J0158+65, we find $L_{\gamma,iso}<5.5\times 10^{47}$ erg/s over 1 s, whereas $L_{\gamma,iso}<10^{49}-10^{51}$ erg/s for the bulk of FRBs. The same values scale up by a factor of ~100 for a ms-long emission. Even on a timescale comparable with that of the radio pulse itself no keV-MeV emission is observed. A systematic association with either long or short GRBs is ruled out with high confidence, except for subluminous events, as is the case for core-collapse of massive stars (long) or binary neutron star mergers (short) viewed off axis. Only giant flares from extra-galactic magnetars at least ten times more energetic than Galactic siblings are ruled out for the nearest FRB.

The Medium Energy X-ray Telescope (ME) is one of the main payloads of the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT). The background of Insight-HXMT/ME is mainly caused by the environmental charged particles and the background intensity is modulated remarkably by the geomagnetic field, as well as the geographical location. At the same geographical location, the background spectral shape is stable but the intensity varies with the level of the environmental charged particles. In this paper, we develop a model to estimate the ME background based on the ME database that is established with the two-year blank sky observations of the high Galactic latitude. In this model, the entire geographical area covered by Insight-HXMT is divided into grids of $5^{\circ}\times5^{\circ}$ in geographical coordinate system. For each grid, the background spectral shape can be obtained from the background database and the intensity can be corrected by the contemporary count rate of the blind FOV detectors. Thus the background spectrum can be obtained by accumulating the background of all the grids passed by Insight-HXMT during the effective observational time. The model test with the blank sky observations shows that the systematic error of the background estimation in $8.9-44.0$ keV is $\sim1.3\%$ for a pointing observation with an average exposure $\sim5.5$ ks. We also find that the systematic error is anti-correlated with the exposure, which indicates the systematic error is partly contributed by the statistical error of count rate measured by the blind FOV detectors.

We report on our analysis of the 2019 outburst of the X-ray accreting pulsar 4U 1901+03 observed with Insight-HXMT and NICER. Both spectra and pulse profiles evolve significantly in the decaying phase of the outburst. Dozens of flares are observed throughout the outburst. They are more frequent and brighter at the outburst peak. We find that the flares, which have a duration from tens to hundreds of seconds, are generally brighter than the persistent emission by a factor of $\sim$ 1.5. The pulse profile shape during the flares can be significantly different than that of the persistent emission. In particular, a phase shift is clearly observed in many cases. We interpret these findings as direct evidence of changes of the pulsed beam pattern, due to transitions between the sub- and super-critical accretion regimes on a short time scale. We also observe that at comparable luminosities the flares' pulse profiles are rather similar to those of the persistent emission. This indicates that the accretion on the polar cap of the neutron star is mainly determined by the luminosity, i.e., the mass accretion rate.

A new black hole X-ray binary (BHXRB) MAXI J1535-571 was discovered by MAXI during its outburst in 2017. Using observations taken by the first Chinese X-ray satellite, the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT), we perform a joint spectral analysis (2-150 keV) in both energy and time domains. The energy spectra provide the essential input for probing the intrinsic Quasi-Periodic Oscillation (QPO) fractional rms spectra (FRS). Our results show that during the intermediate state, the energy spectra are in general consistent with those reported by Swift/XRT and NuSTAR. However, the QPO FRS become harder and the FRS residuals may suggest the presence of either an additional power-law component in the energy spectrum or a turn-over in the intrinsic QPO FRS at high energies.

We study the 2018 outburst of Aql X-1 via the monitor of all sky X-ray image (MAXI) data. We show that the outburst starting in February 2018 is a member of short-low class in the frame of outburst duration and the peak count rate although the outburst morphology is slightly different from the other fast-rise-exponential-decay (FRED) type outbursts with a milder rising stage. We study the partial accretion in the weak propeller stage of Aql X-1 via the MAXI data of the 2018 outburst. We report on the spectral analysis of 3 observations of Aquila X-1 obtained by Insight - hard X-ray modulation telescope (Insight-HXMT) during the late decay stage of the 2018 outburst. We discuss that the data taken by Insight-HXMT is just after the transition to the weak propeller stage. Our analysis shows the necessity of a comptonization component to take into account the existence of an electron cloud resulting photons partly up-scattered.

As China's first X-ray astronomical satellite, the Hard X-ray Modulation Telescope (HXMT), which was dubbed as Insight-HXMT after the launch on June 15, 2017, is a wide-band (1-250 keV) slat-collimator-based X-ray astronomy satellite with the capability of all-sky monitoring in 0.2-3 MeV. It was designed to perform pointing, scanning and gamma-ray burst (GRB) observations and, based on the Direct Demodulation Method (DDM), the image of the scanned sky region can be reconstructed. Here we give an overview of the mission and its progresses, including payload, core sciences, ground calibration/facility, ground segment, data archive, software, in-orbit performance, calibration, background model, observations and some preliminary results.

We present a detailed timing study of the brightest persistent X-ray source Sco X-1 using the data collected by the Hard X-ray Modulation Telescope ($Insight$-HXMT) from July 2017 to August 2018. A complete $Z$-track hardness-intensity diagram (HID) is obtained. The normal branch oscillations (NBOs) at $\sim$ 6 Hz in the lower part of the normal branch (NB) and the flare branch oscillations (FBOs) at $\sim$ 16 Hz in the beginning part of the flaring branch (FB) are found in observations with the Low Energy X-ray Telescope (LE) and the Medium Energy X-ray Telescope (ME) of $Insight$-HXMT, while the horizontal branch oscillations (HBOs) at $\sim$ 40 Hz and the kilohertz quasi-periodic oscillations (kHz QPOs) at $\sim$ 800 Hz are found simultaneously up to 60 keV for the first time on the horizontal branch (HB) by the High Energy X-ray Telescope (HE) and ME. We find that for all types of the observed QPOs, the centroid frequencies are independent of energy, while the root mean square (rms) increases with energy; the centroid frequencies of both the HBOs and kHz QPOs increase along the $Z$-track from the top to the bottom of the HB; and the NBOs show soft phase lags increasing with energy. A continuous QPO transition from the FB to NB in $\sim$ 200 s are also detected. Our results indicate that the non-thermal emission is the origin of all types of QPOs, the innermost region of the accretion disk is non-thermal in nature, and the corona is nonhomogeneous geometrically.

The evidences for the influence of thermonuclear (type-I) X-ray bursts upon the surrounding environments in neutron star low-mass X-ray binaries (LMXB) were detected previously via spectral and timing analyses. Benefitting from a broad energy coverage of Insight-HXMT, we analyze one photospheric radius expansion (PRE) burst, and find an emission excess at soft X-rays. Our spectral analysis shows that, such an excess is not likely relevant to the disk reflection induced by the burst emission and can be attributed to an enhanced pre-burst/persistent emission. We find that the burst and enhanced persistent emissions sum up to exceed Eddington luminosity by $\sim$ 40 percentages. We speculate that the enhanced emission is from a region beyond the PRE radius, or through the Comptonization of the corona.

Corona cooling was detected previously from stacking a series of short type-I bursts occurred during the low/had state of atoll outburst. Type-I bursts are hence regarded as sharp probe to our better understanding on the basic property of the corona. The launch of the first Chinese X-ray satellite Insight-HXMT has large detection area at hard X-rays which provide almost unique chance to move further in this research field. We report the first detection of the corona cooling by Insight-HXMT from single short type-I burst showing up during {\bf flare} of 4U 1636-536. This type-I X-ray burst has a duration of $\sim$13 seconds and hard X-ray shortage is detected with significance 6.2~$\sigma$ in 40-70 keV. A cross-correlation analysis between the lightcurves of soft and hard X-ray band, shows that the corona shortage lag the burst emission by 1.6 $\pm$1.2~s. These results are consistent with those derived previously from stacking a large amount of bursts detected by RXTE/PCA within a series of {\bf flares} of 4U 1636-536. Moreover, the broad bandwidth of Insight-HXMT allows as well for the first time to infer the burst influence upon the continuum spectrum via performing the spectral fitting of the burst, which ends up with the finding that hard X-ray shortage appears at around 40 keV in the continuum spectrum. These results suggest that the evolution of the corona along with the outburst{\bf /flare} of NS XRB may be traced via looking into a series of embedded type-I bursts by using Insight-HXMT.

The Hard X-ray Modulation Telescope (HXMT) named Insight is China's first X-ray astronomical satellite. The Low Energy X-ray Telescope (LE) is one of its main payloads onboard. The detectors of LE adopt swept charge device CCD236 with L-shaped transfer electrodes. Charges in detection area are read out continuously along specific paths, which leads to a time response distribution of photons readout time. We designed a long exposure readout mode to measure the time response distribution. In this mode, CCD236 firstly performs exposure without readout, then all charges generated in preceding exposure phase are read out completely. Through analysis of the photons readout time in this mode, we obtained the probability distribution of photons readout time.

The Medium Energy X-ray telescope (ME) is one of the three main telescopes on board the Insight Hard X-ray Modulation Telescope (Insight-HXMT) astronomy satellite. ME contains 1728 pixels of Si-PIN detectors sensitive in 5-30 keV with a total geometrical area of 952 cm2. Application Specific Integrated Circuit (ASIC) chips, VA32TA6, is used to achieve low power consumption and low readout noise. The collimators define three kinds of field of views (FOVs) for the telescope, 1{\deg}{\times}4{\deg}, 4{\deg}{\times}4{\deg}, and blocked ones. Combination of such FOVs can be used to estimate the in-orbit X-ray and particle background components. The energy resolution of ME is ~3 keV at 17.8 keV (FWHM) and the time resolution is 255 {\mu}s. In this paper, we introduce the design and performance of ME.

We present a study of timing properties of the accreting pulsar 2S 1417-624 observed during its 2018 outburst, based on Swift/BAT, Fermi/GBM, Insight-HXMT and NICER observations. We report a dramatic change of the pulse profiles with luminosity. The morphology of the profile in the range 0.2-10.0keV switches from double to triple peaks at $\sim2.5$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$ and from triple to quadruple peaks at $\sim7$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$. The profile at high energies (25-100keV) shows significant evolutions as well. We explain this phenomenon according to existing theoretical models. We argue that the first change is related to the transition from the sub to the super-critical accretion regime, while the second to the transition of the accretion disc from the gas-dominated to the radiation pressure-dominated state. Considering the spin-up as well due to the accretion torque, this interpretation allows to estimate the magnetic field self-consistently at $\sim7\times 10^{12}$G.

The long-term evolution of the centroid energy of the CRSF in Her X-1 is still a mystery. We report a new measurement from a campaign between {\sl Insight}-HXMT and {\sl NuSTAR} performed in February 2018. Generally, the two satellites show well consistent results of timing and spectral properties. The joint spectral analysis confirms that the previously observed long decay phase has ended, and that the line energy instead keeps constant around 37.5 keV after flux correction.

We report on analysis of observations of the bright transient X-ray pulsar \src obtained during its 2017-2018 giant outburst with Insight-HXMT, \emph{NuSTAR}, and \textit{Swift} observatories. We focus on the discovery of a sharp state transition of the timing and spectral properties of the source at super-Eddington accretion rates, which we associate with the transition of the accretion disk to a radiation pressure dominated (RPD) state, the first ever directly observed for magnetized neutron star. This transition occurs at slightly higher luminosity compared to already reported transition of the source from sub- to super-critical accretion regime associate with onset of an accretion column. We argue that this scenario can only be realized for comparatively weakly magnetized neutron star, not dissimilar to other ultra-luminous X-ray pulsars (ULPs), which accrete at similar rates. Further evidence for this conclusion is provided by the non-detection of the transition to the propeller state in quiescence which strongly implies compact magnetosphere and thus rules out magnetar-like fields.

In this work, we report the in-orbit demonstration of X-ray pulsar navigation with Insight-Hard X-ray Modulation Telescope (Insight-HXMT), which was launched on Jun. 15th, 2017. The new pulsar navigation method 'Significance Enhancement of Pulse-profile with Orbit-dynamics' (SEPO) is adopted to determine the orbit with observations of only one pulsar. In this test, the Crab pulsar is chosen and observed by Insight-HXMT from Aug. 31th to Sept. 5th in 2017. Using the 5-day-long observation data, the orbit of Insight-HXMT is determined successfully with the three telescopes onboard - High Energy X-ray Telescope (HE), Medium Energy X-ray Telescope (ME) and Low Energy X-ray Telescope (LE) - respectively. Combining all the data, the position and velocity of the Insight-HXMT are pinpointed to within 10 km (3 sigma) and 10 m/s (3 sigma), respectively.

The recently discovered neutron star transient Swift J0243.6+6124 has been monitored by {\it the Hard X-ray Modulation Telescope} ({\it Insight-\rm HXMT). Based on the obtained data, we investigate the broadband spectrum of the source throughout the outburst. We estimate the broadband flux of the source and search for possible cyclotron line in the broadband spectrum. No evidence of line-like features is, however, found up to $\rm 150~keV$. In the absence of any cyclotron line in its energy spectrum, we estimate the magnetic field of the source based on the observed spin evolution of the neutron star by applying two accretion torque models. In both cases, we get consistent results with $B\rm \sim 10^{13}~G$, $D\rm \sim 6~kpc$ and peak luminosity of $\rm >10^{39}~erg~s^{-1}$ which makes the source the first Galactic ultraluminous X-ray source hosting a neutron star.

In this paper we present the enhanced X-ray Timing and Polarimetry mission - eXTP. eXTP is a space science mission designed to study fundamental physics under extreme conditions of density, gravity and magnetism. The mission aims at determining the equation of state of matter at supra-nuclear density, measuring effects of QED, and understanding the dynamics of matter in strong-field gravity. In addition to investigating fundamental physics, eXTP will be a very powerful observatory for astrophysics that will provide observations of unprecedented quality on a variety of galactic and extragalactic objects. In particular, its wide field monitoring capabilities will be highly instrumental to detect the electro-magnetic counterparts of gravitational wave sources. The paper provides a detailed description of: (1) the technological and technical aspects, and the expected performance of the instruments of the scientific payload; (2) the elements and functions of the mission, from the spacecraft to the ground segment.

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.

Finding the electromagnetic (EM) counterpart of binary compact star merger, especially the binary neutron star (BNS) merger, is critically important for gravitational wave (GW) astronomy, cosmology and fundamental physics. On Aug. 17, 2017, Advanced LIGO and \textit{Fermi}/GBM independently triggered the first BNS merger, GW170817, and its high energy EM counterpart, GRB 170817A, respectively, resulting in a global observation campaign covering gamma-ray, X-ray, UV, optical, IR, radio as well as neutrinos. The High Energy X-ray telescope (HE) onboard \textit{Insight}-HXMT (Hard X-ray Modulation Telescope) is the unique high-energy gamma-ray telescope that monitored the entire GW localization area and especially the optical counterpart (SSS17a/AT2017gfo) with very large collection area ($\sim$1000 cm$^2$) and microsecond time resolution in 0.2-5 MeV. In addition, \textit{Insight}-HXMT quickly implemented a Target of Opportunity (ToO) observation to scan the GW localization area for potential X-ray emission from the GW source. Although it did not detect any significant high energy (0.2-5 MeV) radiation from GW170817, its observation helped to confirm the unexpected weak and soft nature of GRB 170817A. Meanwhile, \textit{Insight}-HXMT/HE provides one of the most stringent constraints (~10$^{-7}$ to 10$^{-6}$ erg/cm$^2$/s) for both GRB170817A and any other possible precursor or extended emissions in 0.2-5 MeV, which help us to better understand the properties of EM radiation from this BNS merger. Therefore the observation of \textit{Insight}-HXMT constitutes an important chapter in the full context of multi-wavelength and multi-messenger observation of this historical GW event.

We study the X-ray phases of the Crab pulsar utilizing the 11-year observations from the Rossi X-ray Timing Explorer, 6-year radio observations from the Nanshan Telescope, and the ephemeris from Jodrell Bank Observatory. It is found that the X-ray phases in different energy bands and the radio phases from Nanshan Telescope show similar behaviors, including long-time evolution and short-time variations. Such strong correlations between the X-ray and radio phases imply that the radio and X-ray timing noises are both generated from the pulsar spin that cannot be well described by the the monthly ephemeris from the Jodrell Bank observatory. When using the Nanshan phases as references to study the X-ray timing noise, it has a significantly smaller variation amplitude and shows no long-time evolution, with a change rate of $(-1.1\pm1.1)\times10^{-7}$ periods per day. These results show that the distance of the X-ray and radio emission regions on the Crab pulsar has no detectable secular change, and it is unlikely that timing-noises resulted from any unique physical processes in the radio or X-ray emitting regions. The similar behaviors of the X-ray and radio timing noises also imply that the variation of the interstellar medium is not the origin of the Crab pulsar's timing noises, which is consistent with the results obtained from the multi-frequency radio observations of PSR B1540-06.

Integrated pulse profiles at 8.6~GHz obtained with the Shanghai Tian Ma Radio Telescope (TMRT) are presented for a sample of 26 pulsars. Mean flux densities and pulse width parameters of these pulsars are estimated. For eleven pulsars these are the first high-frequency observations and for a further four, our observations have a better signal-to-noise ratio than previous observations. For one (PSR J0742-2822) the 8.6~GHz profiles differs from previously observed profiles. A comparison of 19 profiles with those at other frequencies shows that in nine cases the separation between the outmost leading and trailing components decreases with frequency, roughly in agreement with radius-to-frequency mapping, whereas in the other ten the separation is nearly constant. Different spectral indices of profile components lead to the variation of integrated pulse profile shapes with frequency. In seven pulsars with multi-component profiles, the spectral indices of the central components are steeper than those of the outer components. For the 12 pulsars with multi-component profiles in the high-frequency sample, we estimate the core width using gaussian fitting and discuss the width-period relationship.

The emission geometries, e.g. the emission region height, the beam shape, and radius-to-frequency mapping, are important predictions of pulsar radiation model. The multi-band radio observations carry such valuable information. In this paper, we study two bright pulsars, (PSRs B0329+54 and B1642-03) and observe them in high frequency (2.5 GHz, 5 GHz, and 8 GHz). The newly acquired data together with historical archive provide an atlas of multi-frequency profiles spanning from 100 MHz to 10 GHz. We study the frequency evolution of pulse profiles and the radiation regions with the these data. We firstly fit the pulse profiles with Gaussian functions to determine the phase of each component, and then calculate the radiation altitudes of different emission components and the radiation regions. We find that the inverse Compton scattering (ICS) model can reproduce the radiation geometry of these two pulsars. But for PSR B0329+54 the radiation can be generated in either annular gap (AG) or core gap (CG), while the radiation of PSR B1642-03 can only be generated in the CG. This difference is caused by the inclination angle and the impact angle of these two pulsars. The relation of beaming angle (the angle between the radiation direction and the magnetic axis) and the radiation altitudes versus frequency is also presented by modelling the beam-frequency evolution in the ICS model. The multi-band pulse profiles of these two pulsars can be described well by the ICS model combined with the CG and AG.

Astrophysical techniques have pioneered the discovery of neutrino mass properties. Current cosmological observations give an upper bound on neutrino masses by attempting to disentangle the small neutrino contribution from the sum of all matter using precise theoretical models. We discover the differential neutrino condensation effect in our TianNu N-body simulation. Neutrino masses can be inferred using this effect by comparing galaxy properties in regions of the universe with different neutrino relative abundance (i.e. the local neutrino to cold dark matter density ratio). In "neutrino-rich"' regions, more neutrinos can be captured by massive halos compared to "neutrino-poor" regions. This effect differentially skews the halo mass function and opens up the path to independent neutrino mass measurements in current or future galaxy surveys.

The pulse profile of PSR B1133+16 is usually regarded as a conal-double structure. However, its multifrequency profiles cannot simply be fitted with two Gaussian functions, and a third component is always needed to fit the bridge region (between two peaks). This would introduce additional, redundant parameters. In this paper, through a comparison of five fitting functions (Gaussian, von Mises, hyperbolic secant, square hyperbolic secant, and Lorentz), it is found that the square hyperbolic secant function can best reproduce the profile, yielding an improved fit. Moreover, a symmetric 2D radiation beam function, instead of a simple 1D Gaussian function, is used to fit the profile. Each profile with either well-resolved or not-so-well-resolved peaks could be fitted adequately using this beam function, and the bridge emission between the two peaks does not need to be a new component. Adopting inclination and impact angles based on polarization measurements, the opening angle ({\theta}_{\mu}0) of the radiation beam in a certain frequency band is derived from beam-function fitting. The corresponding radiation altitudes are then calculated. Based on multi-frequency profiles, we also computed the Lorentz factors of the particles and their dispersion at those locations in both the curvature-radiation (CR) and inverse-Compton-scattering (ICS) models. We found that the Lorentz factors of the particles decrease rapidly as the radiation altitude increases. Besides, the radiation prefers to be generated in annular region rather than core region, and this needs further validation.

Due to the lack of long term pulsed emission in quiescence and the strong timing noise, it is impossible to directly measure the braking index $n$ of a magnetar. Based on the estimated ages of their potentially associated supernova remnants (SNRs), we estimate the values of the mean braking indices of eight magnetars with SNRs, and find that they cluster in a range of $1\sim$42. Five magnetars have smaller mean braking indices of $1<n<3$, and we interpret them within a combination of magneto-dipole radiation and wind aided braking, while the larger mean braking indices of $n>3$ for other three magnetars are attributed to the decay of external braking torque, which might be caused by magnetic field decay. We estimate the possible wind luminosities for the magnetars with $1<n<3$, and the dipolar magnetic field decay rates for the magnetars with $n>3$ within the updated magneto-thermal evolution models. Although the constrained range of the magnetars' braking indices is tentative, due to the uncertainties in the SNR ages, which come from distance uncertainties and the unknown conditions of the expanding shells, our method provides an effective way to constrain the magnetars' braking indices if the measurements of the SNRs' ages are reliable, which can be improved by future observations.

PSR B1821$-$24 is a solitary millisecond pulsar (MSP) which radiates multi-wavelength pulsed photons. It has complex radio, X-ray and $\gamma$-ray pulse profiles with distinct peak phase-separations that challenge the traditional caustic emission models. Using the single-pole annular gap model with suitable magnetic inclination angle ($\alpha=40^\circ$) and viewing angle ($\zeta=75^\circ$), we managed to reproduce its pulse profiles of three wavebands. It is found that the middle radio peak is originated from the core gap region at high altitudes, and the other two radio peaks are originated from the annular gap region at relatively low altitudes. Two peaks of both X-ray and $\gamma$-ray wavebands are fundamentally originated from annular gap region, while the $\gamma$-ray emission generated from the core gap region contributes somewhat to the first $\gamma$-ray peak. Precisely reproducing the multi-wavelength pulse profiles of PSR B1821$-$24 enables us to understand emission regions of distinct wavebands and justify pulsar emission models.

PSR J0218$+$4232 is a millisecond pulsar (MSP) with a flux density $\sim$ 0.9 mJy at 1.4 GHz. It is very bright in the high-energy X-ray and $\gamma$-ray domains. We conducted an astrometric program using the European VLBI Network (EVN) at 1.6 GHz to measure its proper motion and parallax. A model-independent distance would also help constrain its $\gamma$-ray luminosity. We achieved a detection of signal-to-noise ratio S/N > 37 for the weak pulsar in all five epochs. Using an extragalactic radio source lying 20 arcmin away from the pulsar, we estimate the pulsar's proper motion to be $\mu_{\alpha}\cos\delta=5.35\pm0.05$ mas yr$^{-1}$ and $\mu_{\delta}=-3.74\pm 0.12$ mas yr$^{-1}$, and a parallax of $\pi=0.16\pm0.09$ mas. The very long baseline interferometry (VLBI) proper motion has significantly improved upon the estimates from long-term pulsar timing observations. The VLBI parallax provides the first model-independent distance constraints: $d=6.3^{+8.0}_{-2.3}$ kpc, with a corresponding $3\sigma$ lower-limit of $d=2.3$ kpc. This is the first pulsar trigonometric parallax measurement based solely on EVN observations. Using the derived distance, we believe that PSR J0218$+$4232 is the most energetic $\gamma$-ray MSP known to date. The luminosity based on even our 3$\sigma$ lower-limit distance is high enough to pose challenges to the conventional outer gap and slot gap models.

Pulsed gamma-ray emission from millisecond pulsars (MSPs) has been detected by the sensitive Fermi, which sheds light on studies of the emission region and mechanism. In particular, the specific patterns of radio and gamma-ray emission from PSR J0101-6422 challenge the popular pulsar models, e.g. outer gap and two-pole caustic (TPC) models. Using the three dimension (3D) annular gap model, we have jointly simulated radio and gamma-ray light curves for three representative MSPs (PSR J0034-0534, PSR J0101-6422 and PSR J0437-4715) with distinct radio phase lags and present the best simulated results for these MSPs, particularly for PSR J0101-6422 with complex radio and gamma-ray pulse profiles and for PSR J0437-4715 with a radio interpulse. It is found that both the gamma-ray and radio emission originate from the annular gap region located in only one magnetic pole, and the radio emission region is not primarily lower than the gamma-ray region in most cases. In addition, the annular gap model with a small magnetic inclination angle instead of "orthogonal rotator" can account for MSPs' radio interpulse with a large phase separation from the main pulse. The annular gap model is a self-consistent model not only for young pulsars but also MSPs, and multi-wavelength light curves can be fundamentally explained by this model.

The concept of a "magnetar" was proposed mainly because of two factors. First, the X-ray luminosity of Anomalous X-ray Pulsars (AXPs) and Soft Gamma-Ray Repeaters (SGRs) is larger than the rotational energy loss rate, and second, the magnetic field strength calculated from "normal method" is super strong. It is proposed that the radiation energy of magnetar comes from its magnetic fields. Here it is argued that the magnetic field strength calculated through the normal method is incorrect when X-ray luminosity is larger than rotational energy loss rate, because the wind braking is not taken into account. Besides, the "anti-magnetar" and some other X-ray and radio observations are difficult to understand with a magnetar model. Instead of the magnetar, we propose a "quarctar", which is a crusted quark star in an accretion disk, to explain the observations. In this model, the persistent X-ray emission, burst luminosity, spectrum of AXPs and SGRs can be understood naturally. The radio-emitting AXPs, which are challenging the magnetar, can also be explained by the quarctar model.

The multi-wavelength pulsed emission from young pulsars and millisecond pulsars can be well modeled with the single-pole 3-dimension annular gap and core gap model. To distinguish our single magnetic pole model from two-pole models (e.g. outer gap model and two-pole caustic model), the convincing values of the magnetic inclination angle and the viewing angle will play a key role.

In the framework of the three-dimensional (3D) annular gap model with reasonable parameters (the magnetic inclination angle \alpha = 45 deg and the view angle \zeta = 63 deg), we first use the latest hight energy data to self-consistently calculate radio, X-ray, gamma-ray and TeV (MAGIC and VERITAS) light curves, phase-averaged spectrum and phase-resolved spectra for the Crab pulsar. It is found that the acceleration electric field and potential in the annular gap and core gap are huge enough in the several tens of neutron star radii. The pulsed emission of radio, X-ray, gamma-ray and TeV are mainly generated from the emission of primary particles or secondary particles with different emission mechanisms in the nearly similar region of the annular gap located in the only one magnetic pole, which leads to the nearly "phase-aligned" multi-wavelength light curves. The emission of peak 1 (P1) and peak 2 (P2) is originated from the annular gap region near the null charge surface, while the emission of bridge is mainly originated from the core gap region. The phase-averaged spectrum and phase-resolved spectra of the Crab pulsar from soft X-ray to TeV band are produced by four components: synchrotron radiation from CR-induced and ICS-induced pairs dominates the X-ray band to soft gamma-ray band (100 eV to 10 MeV); curvature radiation and synchrotron radiation from the primary particles mainly contribute to gamma-ray band (10 MeV to \sim 20 GeV); ICS from the pairs significantly contributes to the TeV gamma-ray band (\sim 20 GeV to 400 GeV). The multi-wavelength pulsed emission from the Crab pulsar has been well modeled with the annular gap and core gap model. To distinguish our single magnetic pole model from two-pole models, the convincing values of the magnetic inclination angle and the viewing angle will play a key role.

Recently, type Ia supernovae data appear to support a dark energy whose equation of state $w$ crosses -1, which is a much more amazing problem than the acceleration of the universe. We show that it is possible for the equation of state to cross the phantom divide by a scalar field in the gravity with an additional inverse power-law term of Ricci scalar in the Lagrangian. The necessary and sufficient condition for a universe in which the dark energy can cross the phantom divide is obtained. Some analytical solutions with $w<-1$ or $w>-1$ are obtained. A minimal coupled scalar with different potentials, including quadratic, cubic, quantic, exponential and logarithmic potentials are investigated via numerical methods, respectively. All these potentials lead to the crossing behavior. We show that it is a robust result which is hardly dependent on the concrete form of the potential of the scalar.

The Vela pulsar represents a distinct group of {\gamma}-ray pulsars. Fermi {\gamma}-ray observations reveal that it has two sharp peaks (P1 and P2) in the light curve with a phase separation of 0.42 and a third peak (P3) in the bridge. The location and intensity of P3 are energy-dependent. We use the 3D magnetospheric model for the annular gap and core gap to simulate the {\gamma}-ray light curves, phase-averaged and phase-resolved spectra. We found that the acceleration electric field along a field line in the annular gap region decreases with heights. The emission at high energy GeV band is originated from the curvature radiation of accelerated primary particles, while the synchrotron radiation from secondary particles have some contributions to low energy {\gamma}-ray band (0.1 - 0.3 GeV). The {\gamma}-ray light curve peaks P1 and P2 are generated in the annular gap region near the altitude of null charge surface, whereas P3 and the bridge emission is generated in the core gap region. The intensity and location of P3 at different energy bands depend on the emission altitudes. The radio emission from the Vela pulsar should be generated in a high-altitude narrow regions of the annular gap, which leads to a radio phase lag of ~ 0.13 prior to the first {\gamma}-ray peak.

Determining radiation location observationally plays a very important role in testing the pulsar radiation models. One-photon pair production in the strong magnetic field, $\gamma-e^{+}e^{1}$, is one of the important physical processes in pulsar radiation mechanisms. Photons near pulsar surface with sufficient energy will be absorbed in the magnetosphere and the absorption optical depth for these GeV $\gamma$-ray photons is usually large. In this paper, we include the aberrational, rotational and general relativistic effects and calculate the $\gamma$-B optical depth for $\gamma$-ray photons. Then we use the derived optical depth to determine the radiation altitude lower bounds for photons with given energies. As a case study, we calculate the lower bounds of radiation altitudes of Crab pulsar for photons with energy from 5 GeV to 1 TeV.

Anomalous X-ray Pulsars and Soft Gamma-Ray Repeaters have been generally recognized as neutron stars with super strong magnetic fields, namely "magnetars". The "magnetars" manifest that the luminosity in X-ray band are larger than the rotational energy loss rate, i.e. $L_{X}>\dot {E}_{\rm rot}$, and then the radiation energy is coming from the energy of magnetic field. Here it is argued that magnetars may not really exist. Some X-ray and radio observational results are contradicted with the magnetar model. (1) The X-ray luminosity of PSR J1852+0040 is much larger than the rotational energy loss rate ($L_{X}/\dot {E}_{\rm rot}\simeq 18)$, but the magnetic field is just $3.1\times 10^{11}$ G. Does this X-ray radiation energy come from the magnetic field? (2) In contrast to the above, the magnetic fields of radio pulsars J1847-0130 and PSR J1718-3718 are higher than that of AXP 1E 2259+586, why is the radiation energy of those two radio pulsars still coming from rotational energy? Furthermore, the magnetic field of the newly discovered SGR 0418+5729 with the lowest magnetic field is 3.0 e13 G, lower than the critical magnetic field $B_{\rm C}=4.414$ e13 G) (Esposito et al. 2010). (3) Some "magnetars" also emit normal transient radio pulses, what is the essential difference between radio pulsars and the "magnetars"? The observational fact arguments will be presented at first, then we discuss in what situation the conventional method to obtain magnetic field could not be correct.

Pulsed high energy radiation from pulsars is not yet completely understood. In this paper, we use the 3D self-consistent annular gap model to study light curves for both young and millisecond pulsars observed by the Fermi Gamma-ray Space Telescope. The annular gap can generate high energy emission for short-period pulsars. The annular gap regions are so large that they have enough electric potential drop to accelerate charged particles to produce gamma-ray photons. For young pulsars, the emission region is from the neutron star surface to about half of the light cylinder radius, and the peak emissivity is in the vicinity of the null charge surface. The emission region for the millisecond pulsars is located much lower than that of the young pulsars. The higher energy gamma-ray emission comes from higher altitudes in the magnetosphere. We present the simulated light curves for three young pulsars (the Crab, the Vela, the Geminga) and three millisecond pulsars (PSR J0030+0451, PSR J0218+4232, PSR J0437-3715) using the annular gap model. Our simulations can reproduce the main properties of observed light curves.

Pulsars have been recognized as normal neutron stars, but sometimes argued as quark stars. {\it Sub-millisecond pulsars, if detected, would play an essential and important role in distinguishing quark stars from neutron stars.} We focus on the formation of such sub-millisecond pulsars in this paper. A new approach to form a sub-millisecond pulsar (quark star) via accretion induced collapse (AIC) of a white dwarf is investigated here. Under this AIC process, we found that: (1) almost all the newborn quark stars could have an initial spin period of $\sim 0.1$ ms; (2) the nascent quark stars (even with a low mass) have sufficiently high spin-down luminosity and satisfy the conditions for pair production and sparking process to be as sub-millisecond radio pulsars; (3) in most cases, the timescales of newborn quark stars in the phase of spin period $< 1$ (or $<0.5$) ms can be long enough to be detected. As a comparison, an accretion spin-up process (for both neutron and quark stars) is also investigated. It is found that, quark stars formed through AIC process can have shorter periods ($\leq$ 0.5 ms); while the periods of neutron stars formed in accretion spin-up process must be longer than 0.5ms. Thus if a pulsar with a period less than 0.5 ms can be identified in the future, it should be a quark star.

Pulsars have been recognized as normal neutron stars or quark stars. Sub-millisecond pulsars, if detected, would play an essential and important role in distinguishing quark stars from neutron stars. A key question is how sub-millisecond pulsars could form. Both sub-Keplerian (for neutron and quark stars) and super-Keplerian cases (only for quark stars, which are bound additionally by strong interaction) have been discussed in this paper in order to investigate possible ways of forming sub-millisecond pulsars. In the sub-Keplerian case, the equilibrium periods of both neutron and quark stars could be as low as ~0.4 ms when they are spun up through accretion in binary systems. In the super-Keplerian case, pulsars could very likely have an initial period of ~0.1 ms if quark stars with different mass could be formed from accretion-induced collapse (AIC) of white dwarfs. The timescale for a sub-millisecond pulsar to keep its period <1 ms is restricted due to gravitational wave radiation. We found that the timescales of neutron stars within sub-millisecond periods are approximately ~10 yr, but the timescales estimated (>10 Myr) for low-mass quark stars could be long enough for us to detect.

The star HE 1305-0007 is a metal-poor double-enhanced star with metallicity [Fe/H] $=-2.0$, which is just at the upper limit of the metallicity for the observed double-enhanced stars. Using a parametric model, we find that almost all s-elements were made in a single neutron exposure. This star should be a member of a post-common-envelope binary. After the s-process material has experienced only one neutron exposure in the nucleosynthesis region and is dredged-up to its envelope, the AGB evolution is terminated by the onset of common-envelope evolution. Based on the high radial-velocity of HE 1305-0007, we speculate that the star could be a runaway star from a binary system, in which the AIC event has occurred and produced the r-process elements.

We present the AGAPE astrometric and photometric catalogue of 1579 variable stars in a 14'x10' field centred on M31. This work is the first survey devoted to variable stars in the bulge of M31. The R magnitudes of the objects and the B-R colours suggest that our sample is dominated by red long-period variable stars (LPV), with a possible overlap with Cepheid-like type II stars. Twelve nova candidates are identified. Correlations with other catalogues suggest that 2 novae could be recurrent novae and provide possible optical counterparts for 2 supersoft X-ray sources candidates observed with Chandra.

Automatic classification of variability is now possible with tools like neural networks. Here, we present two neural networks for the identification of microlensing events -- the first discriminates against variable stars and the second against supernovae. The inputs to the networks include parameters describing the shape and the size of the lightcurve, together with colour of the event. The network computes the posterior probability of microlensing, together with an estimate of the likely error. An algorithm is devised for direct calculation of the microlensing rate from the output of the neural networks. We present a new analysis of the microlensing candidates towards the Large Magellanic Cloud (LMC). The neural networks confirm the microlensing nature of only 7 of the possible 17 events identified by the MACHO experiment. This suggests that earlier estimates of the microlensing optical depth towards the LMC may have been overestimated. A smaller number of events is consistent with the assumption that all the microlensing events are caused by the known stellar populations in the outer Galaxy/LMC.

We have carried out a survey of the Andromeda galaxy for unresolved microlensing (pixel lensing). We present a subset of four short timescale, high signal-to-noise microlensing candidates found by imposing severe selection criteria: the source flux variation exceeds the flux of an R=21 magnitude star and the full width at half maximum timescale is less than 25 days. Remarkably, in three out of four cases, we have been able to measure or strongly constrain the Einstein crossing time of the event. One event, which lies projected on the M31 bulge, is almost certainly due to a stellar lens in the bulge of M31. The other three candidates can be explained either by stars in M31 and M32 or by MACHOs.

This paper pioneers the use of neural networks to provide a fast and automatic way to classify lightcurves in massive photometric datasets. As an example, we provide a working neural network that can distinguish microlensing lightcurves from other forms of variability, such as eruptive, pulsating, cataclysmic and eclipsing variable stars. The network has five input neurons, a hidden layer of five neurons and one output neuron. The five input variables for the network are extracted by spectral analysis from the lightcurve datapoints and are optimised for the identification of a single, symmetric, microlensing bump. The output of the network is the posterior probability of microlensing. The committee of neural networks successfully passes tests on noisy data taken by the MACHO collaboration. When used to process 5000 lightcurves on a typical tile towards the bulge, the network cleanly identifies the single microlensing event. When fed with a sub-sample of 36 lightcurves identified by the MACHO collaboration as microlensing, the network corroborates this verdict in the case of 27 events, but classifies the remaining 9 events as other forms of variability. For some of these discrepant events, it looks as though there are secondary bumps or the bump is noisy or not properly contained. Neural networks naturally allow for the possibility of novelty detection -- that is, new or unexpected phenomena which we may want to follow up. The advantages of neural networks for microlensing rate calculations, as well as the future developments of massive variability surveys, are both briefly discussed.

We report the discovery of a microlensing candidate projected 2'54" from the center of M32, on the side closest to M31. The blue color (R-I= 0.00 +/- 0.14) of the source argues strongly that it lies in the disk of M31, while the proximity of the line of sight to M32 implies that this galaxy is the most likely host of the lens. If this interpretation is correct, it would confirm previous arguments that M32 lies in front of M31. We estimate that of order one such event or less should be present in the POINT-AGAPE data base. If more events are discovered in this direction in a dedicated experiment, they could be used to measure the mass function of M32 up to an unknown scale factor. By combining microlensing observations of a binary-lens event with a measurement of the M31-M32 relative proper motion using the astrometric satellites SIM or GAIA, it will be possible to measure the physical separation of M31 and M32, the last of the six phase-space coordinates needed to assign M32 an orbit.

We present the first results of the analysis of data collected during the 1998-99 observational campaign at the 1.3 meter McGraw-Hill Telescope, towards the Andromeda galaxy (M31), aimed to the detection of gravitational microlensing effects as a probe of the presence of dark matter in our and in M31 halo. The analysis is performed using the pixel lensing technique, which consists in the study of flux variations of unresolved sources and has been proposed and implemented by the AGAPE collaboration. We carry out a shape analysis by demanding that the detected flux variations be achromatic and compatible with a Paczynski light curve. We apply the Durbin-Watson hypothesis test to the residuals. Furthermore, we consider the background of variables sources. Finally five candidate microlensing events emerge from our selection. Comparing with the predictions of a Monte Carlo simulation, assuming a standard spherical model for the M31 and Galactic haloes, and typical values for the MACHO mass, we find that our events are only marginally consistent with the distribution of observable parameters predicted by the simulation.

We report the discovery of a short-duration microlensing candidate in the northern field of the POINT-AGAPE pixel lensing survey towards M31. The full-width half-maximum timescale is very short, just 1.8 days. Almost certainly, the source star has been identified on Hubble Space Telescope archival images, allowing us to infer an Einstein crossing time of 10.4 days, a maximum magnification of about 18, and a lens-source proper motion greater than 0.3 microarcsec/day. The event lies projected at 8' from the center of M31, which is beyond the bulk of the stellar lens population. The lens is likely to reside in one of three locations. It may be a star in the M31 disk, or a massive compact halo object (Macho) in either M31 or the Milky Way. The most probable mass is 0.06 solar masses for an M31 Macho, 0.02 solar masses for a Milky Way Macho and 0.2 solar masses for an M31 stellar lens. Whilst the stellar interpretation is plausible, the Macho interpretation is the most probable for halo fractions above 20%.

POINT-AGAPE is an Anglo-French collaboration which is employing the Isaac Newton Telescope (INT) to conduct a pixel-lensing survey towards M31. In this paper we investigate what we can learn from pixel-lensing observables about the MACHO mass and fractional contribution in M31 and the Galaxy for the case of spherically-symmetric near-isothermal haloes. We employ detailed pixel-lensing simulations which include many of the factors which affect the observables. For a maximum MACHO halo we predict an event rate in V of up to 100 per season for M31 and 40 per season for the Galaxy. However, the Einstein radius crossing time is generally not measurable and the observed full-width half-maximum duration provides only a weak tracer of lens mass. Nonetheless, we find that the near-far asymmetry in the spatial distribution of M31 MACHOs provides significant information on their mass and density contribution. We present a likelihood estimator for measuring the fractional contribution and mass of both M31 and Galaxy MACHOs which permits an unbiased determination to be made of MACHO parameters, even from data-sets strongly contaminated by variable stars. If M31 does not have a significant population of MACHOs in the mass range 0.001-1 Solar masses strong limits will result from the first season of INT observations. Simulations based on currently favoured density and mass values indicate that, after three seasons, the M31 MACHO parameters should be constrained to within a factor four uncertainty in halo fraction and an order of magnitude uncertainty in mass (90% confidence). Interesting constraints on Galaxy MACHOs may also be possible. For a campaign lasting ten years, comparable to the lifetime of current LMC surveys, reliable estimates of MACHO parameters in both galaxies should be possible. (Abridged)

AgapeZ1 is the brightest and the shortest duration microlensing candidate event found in the Agape data. It occured only 42" from the center of M31. Our photometry shows that the half intensity duration of the event6 is 4.8 days and at maximum brightness we measure a stellar magnitude of R=18.0 with B-R=0.80 mag color. A search on HST archives produced a single resolved star within the projected event position error box. Its magnitude is R=22.

M31 is a very tempting target for a microlensing search of compact objects in galactic haloes. It is the nearest large galaxy, it probably has its own dark halo, and its tilted position with respect to the line of sight provides an unmistakable signature of microlensing. However most stars of M31 are not resolved and one has to use the ``pixel method'': monitor the pixels of the image rather than the stars. AGAPE is the implementation of this idea. Data have been collected and treated during two autumns of observation at the 2 metre telescope of Pic du Midi. The process of geometric and photometric alignment, which must be performed before constructing pixel light curves, is described. Seeing variations are minimised by working with large super-pixels (2.1) compared with the average seeing. A high level of stability of pixel fluxes, crucial to the approach, is reached. Fluctuations of super-pixels do not exceed 1.7 times the photon noise which is 0.1% of the intensity for the brightest ones. With such stable data, 10 microlensing events are expected for a full ``standard halo''. With a larger field, a regular and short time sampling and a long lever arm in time, the pixel method will be a very efficient tool to explore the halo of M31.

AGAPE is an observational search of massive compact halo objects (MACHOs) in the direction of M31 by means of a novel method: the gravitational microlensing of unresolved stars. The search consists in examining CCD pixel light curves for microlensing features. The high level of temporal stability necessary to detect microlensing events has been achieved, with quiet pixels stable within a factor of two of the photon noise (the brightest ones down to a level of 0.001 mag). The data analysis is still in progress, but hundreds of variable objects (cepheids, novae, ...) have already been found. Among them there are several lightcurves that resemble microlensing events.

A status report of the microlensing search by the pixel method in the direction of M31, on the 2 meter telescope at Pic du Midi is given. Pixels are stable to a level better than 0.5%. Pixel variations as small as 0.02 magnitude can clearly be detected.