Measuring the dark matter (DM) signals via electron recoil provides an important means for direct detection of light DM particles. The recent XENON1T anomaly with electron recoil energy around $\,E_R^{}\!=\!(2-3)$\,keV can be naturally explained by inelastic scattering of DM which injects energy to the recoiled electrons and gives a narrow peak structure in the recoil spectrum. We present an effective field theory (EFT) analysis of exothermic inelastic DM signals for the Xenon electron recoil detection. For relatively heavy mediator, we fairly formulate the DM-lepton interactions by effective contact operators with two DM fields $(X,\,X')$ and two leptons. Using the XENON1T data, we fit the electron recoil spectrum and constrain the allowed scalar DM mass-splitting as $\,2.1\,\text{keV}\!<\!\Delta m\!<\!3.3\,\text{keV}$ (95\%\,C.L.), with the best fit $\,\Delta{m} \!=\! 2.8\,$keV. We further analyze the relic density produced by such effective DM-electron contact interaction. To provide both the DM relic density and the XENON1T excess, we further derive constraints on the DM mass and the UV cutoff scale of the DM effective interaction. Finally, we study models of UV completion for the effective DM-lepton contact interactions.

Context: Instrumental profile (IP) is the basic property of a spectrograph. Accurate IP characterisation is the prerequisite of accurate wavelength solution. It also facilitates new spectral acquisition methods such as the forward modeling and deconvolution. Aims: We investigate an IP modeling method for the fibre-fed echelle spectrograph with the emission lines of the ThAr lamp, and explore the method to evaluate the accuracy of IP characterisation. Methods: The backbone-residual (BR) model is put forward and tested on the fibre-fed High Resolution Spectrograph (HRS) at the Chinese Xinglong 2.16-m Telescope, which is the sum of the backbone function and the residual function. The backbone function is a bell-shaped function to describe the main component and the spatial variation of IP. The residual function, which is expressed as the cubic spline function, accounts for the difference between the bell-shaped function and the actual IP. The method of evaluating the accuracy of IP characterisation is based on the spectral reconstruction and Monte Carlo simulation. Results: The IP of HRS is characterised with the BR model, and the accuracy of the characterised IP reaches 0.006 of the peak value of the backbone function. This result demonstrates that the accurate IP characterisation has been achieved on HRS with the BR model, and the BR model is an excellent choice for accurate IP characterisation of fibre-fed echelle spectrographs.

We investigate the relationship between the cosmic curvature and the model of dark energy (hereafter DE) with the recent Type Ia supernovae (hereafter SNe Ia) data, i.e., the Pantheon sample including 1048 SNe Ia with $0.01 < z < 2.3$. We obtain the measurements of the dimensionless spatial curvature density today, i.e., $\Omega_{k0} = -0.062^{+0.189}_{-0.169}, -0.004^{+0.228}_{-0.134}, 0.127^{+0.280}_{-0.276}$ and $0.422^{+0.213}_{-0.338}$ at 68\% confidence level (CL), respectively, in the scenarios of $\Lambda$CDM, $\phi$CDM (i.e., scalar field dark energy), $\omega$CDM and $\omega_0\omega_a$CDM models. In the scenario of $\Lambda$CDM model, a closed universe is preferred by the Pantheon sample, which is consistent with that from the Planck CMB spectra. However, the uncertainty of $\Omega_{k0}$ from the Pantheon SNe sample is about 8 times larger than that from the Planck data, so the former one supports a closed universe at a much lower CL than that from the latter one. An open unverse is supported by the Pantheon sample at $\sim$32\% and $\sim$78\% CLs, respectively, in the $\omega$CDM and $\omega_0\omega_a$CDM models. Among these models, the $\phi$CDM model is the one which supports the flat universe most strongly. It shows that $\Omega_{k0}$ is significantly dependent on the adopted model of dark energy, and there is a negative correlation between $\Omega_{k0}$ and the equation of state of DE.

Self-interacting neutrinos that begin to free-stream at close to matter-radiation equality can reduce the physical size of photon sound horizon at last scattering surface. This can be the reason why standard $\Lambda$CDM cosmology sees a lower value of Hubble constant than local measurements from distance ladder. We propose a new realization of self-interacting Dirac neutrinos (SID$\nu$) with light dark photon mediator for a viable interaction mechanism. Our model is UV completed by a Dirac seesaw with anomaly-free dark $U(1)_X^{}$ gauge group which charges the right-handed neutrinos. This naturally generates small masses for Dirac neutrinos and induces self-scattering of right-handed neutrinos. The scattering with left-handed neutrinos is suppressed by a chirality-flip mass insertion when the neutrino energy is much larger than its mass. The resultant neutrino self-scattering is not operative for $E_{\nu}\gtrsim O(\text{keV})$, which avoids the cosmological and laboratory constraints. By evolving Boltzmann equations for left- and right-handed neutrino number densities, we show that about $2/3$ of the left-handed neutrinos are converted into right-handed neutrinos in a short epoch between the Big-Bang Nucleosynthesis and the recombination, and interact with each other efficiently afterwards. The resultant neutrino non-free-streaming is the key ingredient to shrink down the comoving sound horizon at drag epoch, which can reconcile the Hubble tension between early and late time measurements.

We here report a spectroscopic identification of two new changing-look AGNs (CL-AGNs): SDSS\,J104705.16+544405.8 and SDSS\,J120447.91+170256.8 both with a "turn-off" type transition from type 1 to type 1.8/1.9. The identification is arrived by a follow-up spectroscopic observation of the five changing-look AGN (CL-AGN) candidates that are extracted from the sample recently released in Macleod et al. The candidates are extract by the authors from the Sloan Digit Sky Survey Data Release 7 spectroscopically confirmed quasars with large amplitude variability. By compiling a sample of 26 previously identified CL-AGNs, we confirm the claim in Macleod et al. that CL-AGNs tend to be biased against low Eddington ratio, and identify an overlap between the CL-AGNs at their dim state and the so-called intermediate-type AGNs. The overlap implies that there two populations of the intermediate-type AGNs with different origins. One is due to the torus orientation effect, and the another the intrinsic change of the accretion rate of the central supermassive blackholes.

In this paper, three ensemble methods: Random Forest, XGBoost, and a Hybrid Ensemble method were implemented to classify imbalanced pulsar candidates. To assist these methods, tree models were used to select features among 30 features of pulsar candidates from references. The skewness of the integrated pulse profile, chi-squared value for sine-squared fit to amended profile and best S/N value play important roles in Random Forest, while the skewness of the integrated pulse profile is one of the most significant features in XGBoost. More than 20 features were selected by their relative scores and then applied in three ensemble methods. In the Hybrid Ensemble method, we combined Random Forest and XGBoost with EasyEnsemble. By changing thresholds, we tried to make a trade-off between Recall and Precision to make them approximately equal and as high as possible. Experiments on HTRU 1 and HTRU 2 datasets show that the Hybrid Ensemble method achieves higher Recall than the other two algorithms. In HTRU 1 dataset, Recall, Precision, and F-Score of the Hybrid Ensemble method are $0.967$, $0.971$, and $0.969$, respectively. In HTRU 2 dataset, the three values of that are $0.920$, $0.917$, and $0.918$, respectively.

As performance of dedicated facilities continually improved, massive pulsar candidates are being received, which makes selecting valuable pulsar signals from candidates challenging. In this paper, we designed a deep convolutional neural network (CNN) with 11 layers for classifying pulsar candidates. Compared to artificial designed features, CNN chose sub-integrations plot and sub-bands plot in each candidate as inputs without carrying biases. To address the imbalanced problem, data augmentation method based on synthetic minority samples is proposed according to characteristics of pulsars. The maximum pulses of pulsar candidates were first translated to the same position, then new samples were generated by adding up multiple subplots of pulsars. The data augmentation method is simple and effective for obtaining varied and representative samples which keep pulsar characteristics. In the experiments on HTRU 1 dataset, it shows that this model can achieve recall as 0.962 while precision as 0.963.

We present the revised and high quality CCD time series observations in 2017 and 2018 of the W UMa-type contact binary BH~Cas. The combination of our new determinations of minimum times and those collected from the literatures reveal a long-term period increase with a rate of $\textit{dP/dt}\,=\,+3.27 \times 10^{-7} \textrm{d}\,\textrm{yr}^{-1}$, which can be accounted for by angular momentum transfer as the less-massive component loses its mass to the more massive one. One cyclical oscillation ($A'=0.00300$~days and $P'=20.09$~yrs) superimposed on the long-term period increase tendency is attributed to one additional object that orbits BH Cas. The low-resolution spectra of this source are observed at phase $\sim$0.0 and $\sim$0.5 in 2018 and 2019, and the spectral types are identified as K3V$\,\pm\,$1 and G8V$\,\pm\,$2. According to the luminosity contributions and spectra at these phases, the spectral type of the primary star is close to K3V$\,\pm\,$1. We infer a cool spot on the hotter and less-massive component to fit the asymmetry of light curves. In addition, changes in the location, temperature, and area of spots on the secondary star in different observations may indicate that the magnetic field activity of the secondary is more active than that of the primary. Such strong spot activities and optical flare are both detected at the same time from a W UMa-type star for the first time. With the \emph{XMM-Newton} observations taken in 2003 and 2012, the X-ray light curve is obtained and the X-ray spectra can be well described with a black-body or a bremsstrahlung model. The soft X-ray luminosity of BH Cas is estimated as 8.8\,--\,9.7$\,\times\,10^{30}$\,erg\,s$^{-1}$.

We present a new, model-independent method to reconstruct the temporal evolution of the speed of light $c(z)$ using astronomical observations. After validating our pipeline using mock datasets, we apply our method to the latest BAO and supernovae observations, and reconstruct $c(z)$ in the redshift range of $z\in[0,1.5]$. We find no evidence of a varying speed of light, although we see some interesting features of $\Delta c(z)$, the fractional difference between $c(z)$ and $c_0$ (the speed of light in SI), \eg, $\Delta c(z)<0$ and $\Delta c(z)>0$ at $0.2\lesssim z\lesssim0.5$ and $0.8\lesssim z\lesssim1.3$, respectively, although the significance of these features is currently far below statistical importance.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

We propose a model-independent way to determine the cosmic curvature using the Hubble parameter $H(z)$ measured with cosmic chronometers and the comoving distance $D(z)$ inferred from HII galaxies. We employ Gaussian processes to smooth the measure of distance and match it to $30$ values of $H(z)$. The curvature parameter $\Omega_k$ may be obtained individually for each such pair. The weighted average for the complete sample is $\Omega_k=-0.0013\pm0.0004$, suggesting a bias towards negative values. The accuracy of the curvature measurement improves with increased redshift, however, given possible systematic effects associated with local inhomogeneities. We therefore also analyze a high-redshift ($z>1.5$) sub-sample on its own, which is more likely to reflect the geometry of the Universe on large, smooth scales. We find for this set of data that $\Omega_k=-0.0111\pm0.0416$, consistent with zero to better than $1\sigma$. This result is in agreement with the spatially flat universe inferred from the cosmic microwave background observations. We expect this method to yield even tighter constraints on the curvature parameter with future, more accurate observations of HII galaxies at high $z$.

We report first results from the CHinese Exoplanet Searching Program from Antarctica (CHESPA)---a wide-field high-resolution photometric survey for transiting exoplanets carried out using telescopes of the AST3 (Antarctic Survey Telescopes times 3) project. There are now three telescopes (AST3-I, AST3-II, and CSTAR-II) operating at Dome A---the highest point on the Antarctic Plateau---in a fully automatic and remote mode to exploit the superb observing conditions of the site, and its long and uninterrupted polar nights. The search for transiting exoplanets is one of the key projects for AST3. During the Austral winters of 2016 and 2017 we used the AST3-II telescope to survey a set of target fields near the southern ecliptic pole, falling within the continuous viewing zone of the TESS mission \citep{Ricker10}. The first data release of the 2016 data, including images, catalogs and lightcurves of 26578 bright stars ($7.5\le i \le15$) was presented in \citet{Zhang18}. The best precision, as measured by the RMS of the lightcurves at the optimum magnitude of the survey ($i=10$), is around 2\,mmag. We detect 222 objects with plausible transit signals from these data, 116 of which are plausible transiting exoplanet candidates according to their stellar properties as given by the TESS Input Catalog \citep{Stassun17}, Gaia DR2 \citep{Gaia18} and TESS-HERMES spectroscopy \citep{Sharma18}. With the first data release from TESS expected in late 2018, this candidate list will be a timely for improving the rejection of potential false positives.

Located at Dome A, the highest point of the Antarctic plateau, the Chinese Kunlun station is considered to be one of the best ground-based photometric sites because of its extremely cold, dry, and stable atmosphere(Saunders et al. 2009). A target can be monitored from there for over 40 days without diurnal interruption during a polar winter. This makes Kunlun station a perfect site to search for short-period transiting exoplanets. Since 2008, an observatory has been built at Kunlun station and three telescopes are working there. Using these telescopes, the AST3 project has been carried out over the last six years with a search for transiting exoplanets as one of its key programs (CHESPA). In the austral winters of 2016 and 2017, a set of target fields in the Southern CVZ of TESS (Ricker et al. 2009) were monitored by the AST3-II telescope. In this paper, we introduce the CHESPA and present the first data release containing photometry of 26,578 bright stars (m_i < 15). The best photometric precision at the optimum magnitude for the survey is around 2 mmag. To demonstrate the data quality, we also present a catalog of 221 variables with a brightness variation greater than 5 mmag from the 2016 data. Among these variables, 179 are newly identified periodic variables not listed in the AAVSO databasea), and 67 are listed in the Candidate Target List(Stassun et al. 2017). These variables will require careful attention to avoid false-positive signals when searching for transiting exoplanets. Dozens of new transiting exoplanet candidates will be also released in a subsequent paper(Zhang et al. 2018b).

We perform a tomographic structure growth and expansion rate analysis using the monopole, quadrupole and hexadecapole of the redshift-space galaxy power spectrum derived from the Sloan Digital Sky Survey (SDSS-III) Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12 combined sample, which covers the redshift range of $0.20<z<0.75$. By allowing for overlap between neighbouring redshift slices in order to extract information on the light-cone, we successfully obtain joint BAO and RSD constraints with a precision of $2-3\%$ for $D_A$, $3-10\%$ for $H$ and $9-12\%$ for $f\sigma_8$ with a redshift resolution of $\Delta z\sim0.04$. Our measurement is consistent with that presented in arXiv:1709.05173, where the analysis is performed in configuration space. We apply our measurement to constrain the $f(R)$ gravity model, and find that the 95\% CL upper limit of ${\rm log_{10}}B_0$ can be reduced by 11\% by our tomographic BAO and RSD measurement.

The observed rapid cooling of the neutron star (NS) located at the center of the supernova remnant Cassiopeia A (Cas A) can be explained in the minimal NS cooling scenario. This consequence may be changed if there exists an extra cooling source, such as axion emission. In this work, we study the Cas A NS cooling in the presence of axion emission, taking account of the temperature evolution in the whole life of the Cas A NS. We obtain a lower limit on the axion decay constant, $f_a\gtrsim (5-7)\times 10^8$ GeV, if the star has an envelope with a thin carbon layer. This is as strong as existing limits imposed by other astrophysical observations such as SN1987A.

We have developed a new automated small-scale magnetic flux rope (SSMFR) detection algorithm based on the Grad-Shafranov (GS) reconstruction technique. We have applied this detection algorithm to the Wind spacecraft in-situ measurements during 1996 - 2016, covering two solar cycles, and successfully detected a total number of 74,241 small-scale magnetic flux rope events with duration from 9 to 361 minutes. This large number of small-scale magnetic flux ropes has not been discovered by any other previous studies through this unique approach. We perform statistical analysis of the small-scale magnetic flux rope events based on our newly developed database, and summarize the main findings as follows. (1) The occurrence of small-scale flux ropes has strong solar cycle dependency with a rate of a few hundreds per month on average. (2) The small-scale magnetic flux ropes in the ecliptic plane tend to align along the Parker spiral. (3) In low speed ($<$ 400 km/s) solar wind, the flux ropes tend to have lower proton temperature and higher proton number density, while in high speed ($\ge$ 400 km/s) solar wind, they tend to have higher proton temperature and lower proton number density. (4) Both the duration and scale size distributions of the small-scale magnetic flux ropes obey a power law. (5) The waiting time distribution of small-scale magnetic flux ropes can be fitted by an exponential function (for shorter waiting times) and a power law function (for longer waiting times). (6) The wall-to-wall time distribution obeys double power laws with the break point at 60 minutes (corresponding to the correlation length). (7) The small-scale magnetic flux ropes tend to accumulate near the heliospheric current sheet (HCS). The entire database is available at \url{http://fluxrope.info} and in machine readable format in this article.

Sky survey is one of the most important motivations to improve the astrophysics development, especially when using new photometric bands. We are performing the SAGE (Stellar Abundance and Galactic Evolution) survey with a self-designed SAGE photometric system, which is composed of eight photometric bands. The project mainly aims to study the stellar atmospheric parameters of $\sim$0.5 billion stars in the $\sim12,000$ deg$^2$ of the northern sky, which mainly focuses on the Galactic sciences, as well as some extragalactic sciences. This work introduces the detailed data reduction process of the testing field NGC\,6791, including the data reduction of single-exposure image and stacking multi-exposure images, and properties of the final catalogue.

To investigate in more details of Stellar Abundance and Galactic Evolution (SAGE) and in a huge sample, we are performing a northern sky photometric survey named SAGES with the SAGE photometric system, which consists of 8 filters: Str\"omgren-$u$, SAGE-$v$, SDSS $g$, $r$, $i$, DDO-$51$, $H\alpha_{wide}$, and $H\alpha _{narrow}$, including three Sloan broadband filters, three intermediate-band filters and two narrow-band filters, and one newly-designed narrow-band filter. SAGES covers $\sim$12,000 square degrees of the northern sky with $\delta > -5 ^{\circ}$, excluding the Galactic disk ($|b|<10^{\circ}$) and the sky area of 12 hr $<$ R.A. $<$ 18\,hr. The photometric detection limit depth at signal-to-noise ratio $5\sigma$ can be as deep as $V\sim$20\,mag. The SAGES will produce a depth-uniformed photometric catalogue for $\sim$500 million stars with atmospheric parameters including effective temperature $T_{\rm eff}$, surface gravity log\,g, and metallicity [Fe/H], as well as interstellar extinction to each individual target. In this work, we will briefly introduce the SAGE photometric system, the SAGE survey, and a preliminary test field of the open cluster NGC\,6791 and around.

Shack-Hartmann wavefront sensors using wound fibre bundles are desired for multi-object adaptive optical systems to provide large multiplex positioned by Starbugs. The use of the large-sized wound fibre bundle provides the exibility to use more sub-apertures wavefront sensor for ELTs. These compact wavefront sensors take advantage of large focal surfaces such as the Giant Magellan Telescope. The focus of this paper is to study the wound fibre image bundle structure defects effect on the centroid measurement accuracy of a Shack-Hartmann wavefront sensor. We use the first moment centroid method to estimate the centroid of a focused Gaussian beam sampled by a simulated bundle. Spot estimation accuracy with wound fibre image bundle and its structure impact on wavefront measurement accuracy statistics are addressed. Our results show that when the measurement signal to noise ratio is high, the centroid measurement accuracy is dominated by the wound fibre image bundle structure, e.g. tile angle and gap spacing. For the measurement with low signal to noise ratio, its accuracy is influenced by the read noise of the detector instead of the wound fibre image bundle structure defects. We demonstrate this both with simulation and experimentally. We provide a statistical model of the centroid and wavefront error of a wound fibre image bundle found through experiment.

We present unique {and additional} observational evidence for the self-generation of small-scale coherent magnetic flux rope structures in the solar wind. Such structures with durations between 9 and 361 minutes are identified from Wind in-situ spacecraft measurements through the Grad-Shafranov (GS) reconstruction approach. The event occurrence counts are on the order of 3,500 per year on average and have a clear solar cycle dependence. We build a database of small-scale magnetic flux ropes from twenty-year worth of Wind spacecraft data. We show a power-law distribution of the wall-to-wall time corresponding well to the inertial range turbulence, which agrees with relevant observations and numerical simulation results. We also provide the axial current density distribution from the GS-based observational analysis, which yields a non-Gaussian probability density function consistent with numerical simulation results.

The Kepler-9 system harbors three known transiting planets. The system holds significant interest for several reasons. First, the outer two planets exhibit a period ratio that is close to a 2:1 orbital commensurability, with attendant dynamical consequences. Second, both planets lie in the planetary mass "desert" that is generally associated with the rapid gas agglomeration phase of the core accretion process. Third, there exist attractive prospects for accurately measuring both the sky-projected stellar spin-orbit angles as well as the mutual orbital inclination between the planets in the system. Following the original \textit{Kepler} detection announcement in 2010, the initially reported orbital ephemerides for Kepler-9~b and c have degraded significantly, due to the limited time base-line of observations on which the discovery of the system rested. Here, we report new ground-based photometric observations and extensive dynamical modeling of the system. These efforts allow us to photometrically recover the transit of Kepler-9~b, and thereby greatly improve the predictions for upcoming transit mid-times. Accurate ephemerides of this system are important in order to confidently schedule follow-up observations of this system, for both in-transit Doppler measurements as well as for atmospheric transmission spectra taken during transit.

We analyze sunspots rotation and magnetic transients in NOAA AR 11429 during two X-class (X5.4 and X1.3) flares using the data from the Helioseismic and Magnetic Imager on board the \emph{Solar Dynamics Observatory}. A large leading sunspot with positive magnetic polarity rotated counterclockwise. As expected, the rotation was significantly affected by the two flares. The magnetic transients induced by the flares were clearly evident in the sunspots with negative polarity. They were moving across the sunspots with speed of order $3-7\ \rm km \ s^{-1}$. Furthermore, the trend of magnetic flux evolution of these sunspots exhibited changes associated with the flares. These results may shed light on the understanding of the evolution of sunspots.

Imaging bundles provide a convenient way to translate a spatially coherent image, yet conventional imaging bundles made from silica fibre optics typically remain expensive with large losses due to poor filling factors (~40%). We present the characterisation of a novel polymer imaging bundle made from poly(methyl methacrylate) (PMMA) that is considerably cheaper and a better alternative to silica imaging bundles over short distances (~1 m; from the middle to the edge of a telescope's focal plane). The large increase in filling factor (92% for the polymer imaging bundle) outweighs the large increase in optical attenuation from using PMMA (1 dB/m) instead of silica (10^{-3} dB/m). We present and discuss current and possible future multi-object applications of the polymer imaging bundle in the context of astronomical instrumentation including: field acquisition, guiding, wavefront sensing, narrow-band imaging, aperture masking, and speckle imaging. The use of PMMA limits its use in low light applications (e.g. imaging of galaxies), however it is possible to fabricate polymer imaging bundles from a range of polymers that are better suited to the desired science.

The high Antarctic plateau provides exceptional conditions for conducting infrared observations of the cosmos on account of the cold, dry and stable atmosphere above the ice surface. This paper describes the scientific goals behind the first program to examine the time-varying universe in the infrared from Antarctica - the Kunlun Infrared Sky Survey (KISS). This will employ a small (50 cm aperture) telescope to monitor the southern skies in the 2.4um Kdark window from China's Kunlun station at Dome A, on the summit of the Antarctic plateau, through the uninterrupted 4-month period of winter darkness. An earlier paper discussed optimisation of the Kdark filter for the best sensitivity (Li et al 2016). This paper examines the scientific program for KISS. We calculate the sensitivity of the camera for the extrema of observing conditions that will be encountered. We present the parameters for sample surveys that could then be carried out for a range of cadences and sensitivities. We then discuss several science programs that could be conducted with these capabilities, involving star formation, brown dwarfs and hot Jupiters, exoplanets around M dwarfs, the terminal phases of stellar evolution, discovering fast transients as part of multi-wavelength campaigns, embedded supernova searches, reverberation mapping of active galactic nuclei, gamma ray bursts and the detection of the cosmic infrared background. Accepted for publication in PASA, 04/08/16.

The rotation of sunspots of 2 yr in two different solar cycles is studied with the data from the Helioseismic and Magnetic Imager on board the \it Solar Dynamics Observatory \rm and the Michelson Doppler Imager instrument on board the \it Solar and Heliospheric Observataory.\rm We choose the $\alpha$ sunspot groups and the relatively large and stable sunspots of complex active regions in our sample. In the year of 2003, the $\alpha$ sunspot groups and the preceding sunspots tend to rotate counterclockwise and have positive magnetic polarity in the northern hemisphere. In the southern hemisphere, the magnetic polarity and rotational tendency of the $\alpha$ sunspot groups and the preceding sunspots are opposite to the northern hemisphere. The average rotational speed of these sunspots in 2003 is about $0^{\circ}.65 \rm \ hr^{-1}$. From 2014 January to 2015 February, the $\alpha$ sunspot groups and the preceding sunspots tend to rotate clockwise and have negative magnetic polarity in the northern hemisphere. The patterns of rotation and magnetic polarity of the southern hemisphere are also opposite to those of the northern hemisphere. The average rotational speed of these sunspots in 2014/2015 is about $1^{\circ}.49 \rm \ hr^{-1}$. The rotation of the relatively large and stable preceding sunspots and that of the $\alpha$ sunspot groups located in the same hemisphere have opposite rotational direction in 2003 and 2014/2015.

The Kunlun Infrared Sky Survey will be the first comprehensive exploration of the time varying Universe in the infrared. A key feature in optimizing the scientific yield of this ambitious research program is the choice of the survey passband. In particular the survey aims to maximally exploit the unique thermal and atmospheric conditions pertaining to the high Antarctic site. By simulating the expected signal-to-noise for varying filter properties within the so-called "K_DARK" 2.4um window, filter performance can be tuned and best-case designs are given covering a range of conditions.

We explore planetary migration scenarios for formation of high inclination Neptune Trojans (NTs) and how they are affected by the planetary migration of Neptune and Uranus. If Neptune and Uranus's eccentricity and inclination were damped during planetary migration, then their eccentricities and inclinations were higher prior and during migration than their current values. Using test particle integrations we study the stability of primordial NTs, objects that were initially Trojans with Neptune prior to migration. We also study Trans-Neptunian objects captured into resonance with Neptune and becoming NTs during planet migration. We find that most primordial NTs were unstable and lost if eccentricity and inclination damping took place during planetary migration. With damping, secular resonances with Neptune can increase a low eccentricity and inclination population of Trans-Neptunian objects increasing the probability that they are captured into 1:1 resonance with Neptune, becoming high inclination NTs. We suggest that the resonant trapping scenario is a promising and more effective mechanism explaining the origin of NTs that is particularly effective if Uranus and Neptune experienced eccentricity and inclination damping during planetary migration.

The High Efficiency and Resolution Multi Element Spectrograph, HERMES, is a facility-class optical spectrograph for the Anglo-Australian Telescope (AAT). It is designed primarily for Galactic Archaeology, the first major attempt to create a detailed understanding of galaxy formation and evolution by studying the history of our own galaxy, the Milky Way. The goal of the GALAH survey is to reconstruct the mass assembly history of the Milky Way through a detailed chemical abundance study of one million stars. The spectrograph is based at the AAT and is fed by the existing 2dF robotic fiber positioning system. The spectrograph uses volume phase holographic gratings to achieve a spectral resolving power of 28,000 in standard mode and also provides a high-resolution mode ranging between 40,000 and 50,000 using a slit mask. The GALAH survey requires an SNR greater than 100 for a star brightness of V ?= 14 in an exposure time of one hour. The total spectral coverage of the four channels is about 100 nm between 370 and 1000 nm for up to 392 simultaneous targets within the 2-degree field of view. HERMES has been commissioned over three runs, during bright time in October, November, and December 2013, in parallel with the beginning of the GALAH pilot survey, which started in November 2013. We present the first-light results from the commissioning run and the beginning of the GALAH survey, including performance results such as throughput and resolution, as well as instrument reliability.

Hydromagnetic waves, especially those of frequencies in the range of a few milli-Hz to a few Hz observed in the Earth's magnetosphere, are categorized as Ultra Low Frequency (ULF) waves or pulsations. They have been extensively studied due to their importance in the interaction with radiation belt particles and in probing the structures of the magnetosphere. We developed an approach in examining the toroidal standing Aflv\'{e}n waves in a background magnetic field by recasting the wave equation into a Klein-Gordon (KG) form along individual field lines. The eigenvalue solutions to the system are characteristic of a propagation type when the corresponding eigen-frequency is greater than a cut-off frequency and an evanescent type otherwise. We apply the approach to a compressed dipole magnetic field model of the inner magnetosphere, and obtain the spatial profiles of relevant parameters and the spatial wave forms of harmonic oscillations. We further extend the approach to poloidal mode standing Alfv\'{e}n waves along field lines. In particular, we present a quantitative comparison with a recent spacecraft observation of a poloidal standing Alfv\'{e}n wave in the Earth's magnetosphere. Our analysis based on KG equation yields consistent results which agree with the spacecraft measurements of the wave period and the amplitude ratio between the magnetic field and electric field perturbations. We also present computational results of eigenvalue solutions to the compressional poloidal mode waves in the compressed dipole magnetic field geometry.

The ability to position multiple miniaturized wavefront sensors precisely over large focal surfaces are advantageous to multi-object adaptive optics. The Australian Astronomical Observatory (AAO) has prototyped a compact and lightweight Shack-Hartmann wavefront-sensor that fits into a standard Starbug parallel fibre positioning robot. Each device makes use of a polymer coherent fibre imaging bundle to relay an image produced by a microlens array placed at the telescope focal plane to a re-imaging camera mounted elsewhere. The advantages of the polymer fibre bundle are its high-fill factor, high-throughput, low weight, and relatively low cost. Multiple devices can also be multiplexed to a single low-noise camera for cost efficiencies per wavefront sensor. The use of fibre bundles also opens the possibility of applications such as telescope field acquisition, guiding, and seeing monitors to be positioned by Starbugs. We present the design aspects, simulations and laboratory test results.

We demonstrate how observations of pulsars can be used to help navigate a spacecraft travelling in the solar system. We make use of archival observations of millisecond pulsars from the Parkes radio telescope in order to demonstrate the effectiveness of the method and highlight issues, such as pulsar spin irregularities, which need to be accounted for. We show that observations of four millisecond pulsars every seven days using a realistic X-ray telescope on the spacecraft throughout a journey from Earth to Mars can lead to position determinations better than approx. 20km and velocity measurements with a precision of approx. 0.1m/s.

GNOSIS is a prototype astrophotonic instrument that utilizes OH suppression fibres consisting of fibre Bragg gratings and photonic lanterns to suppress the 103 brightest atmospheric emission doublets between 1.47-1.7 microns. GNOSIS was commissioned at the 3.9-meter Anglo-Australian Telescope with the IRIS2 spectrograph to demonstrate the potential of OH suppression fibres, but may be potentially used with any telescope and spectrograph combination. Unlike previous atmospheric suppression techniques GNOSIS suppresses the lines before dispersion and in a manner that depends purely on wavelength. We present the instrument design and report the results of laboratory and on-sky tests from commissioning. While these tests demonstrated high throughput and excellent suppression of the skylines by the OH suppression fibres, surprisingly GNOSIS produced no significant reduction in the interline background and the sensitivity of GNOSIS and IRIS2 is about the same as IRIS2. It is unclear whether the lack of reduction in the interline background is due to physical sources or systematic errors as the observations are detector noise-dominated. OH suppression fibres could potentially impact ground-based astronomy at the level of adaptive optics or greater. However, until a clear reduction in the interline background and the corresponding increasing in sensitivity is demonstrated optimized OH suppression fibres paired with a fibre-fed spectrograph will at least provide a real benefits at low resolving powers.

The background noise between 1 and 1.8 microns in ground-based instruments is dominated by atmospheric emission from hydroxyl molecules. We have built and commissioned a new instrument, GNOSIS, which suppresses 103 OH doublets between 1.47 - 1.7 microns by a factor of ~1000 with a resolving power of ~10,000. We present the first results from the commissioning of GNOSIS using the IRIS2 spectrograph at the AAT. The combined throughput of the GNOSIS fore-optics, grating unit and relay optics is ~36 per cent, but this could be improved to ~46 per cent with a more optimal design. We measure strong suppression of the OH lines, confirming that OH suppression with fibre Bragg gratings will be a powerful technology for low resolution spectroscopy. The integrated OH suppressed background between 1.5 and 1.7 microns is reduced by a factor of 9 compared to a control spectrum using the same system without suppression. The potential of low resolution OH suppressed spectroscopy is illustrated with example observations. The GNOSIS background is dominated by detector dark current below 1.67 microns and by thermal emission above 1.67 microns. After subtracting these we detect an unidentified residual interline component of ~ 860 +/ 210 ph/s/m^2/micron/arcsec^2. This component is equally bright in the suppressed and control spectra. We have investigated the possible source of the interline component, but were unable to discriminate between a possible instrumental artifact and intrinsic atmospheric emission. Resolving the source of this emission is crucial for the design of fully optimised OH suppression spectrographs. The next generation OH suppression spectrograph will be focussed on resolving the source of the interline component, taking advantage of better optimisation for a FBG feed. We quantify the necessary improvements for an optimal OH suppressing fibre spectrograph design.

For pulsar projects it is often necessary to predict the pulse phase in advance, for example, when preparing for new observations. Interpolation of the pulse phase between existing measurements is also often required, for example, when folding X-ray or gamma-ray observations according to the radio pulse phase. Until now these procedures have been done using various ad hoc methods. The purpose of this paper is to show how to interpolate or predict the pulse phase optimally using statistical models of the various noise processes and the phase measurement uncertainty.

We investigate the constraints on the scalar, vector and spin-3/2 dark matter interaction with the standard model particles, from the observations of dark matter relic density, the direct detection experiments of CDMS and XENON, and the indirect detection of the antiproton-to-proton ratio by PAMELA. A model independent way is adopted by constructing general 4-particle operators up to dimension 6 for the effective interaction between dark matter and standard model particles. We find that the constraints from different experiments are complementary with each other. Comparison among these constraints may exclude some effective models of dark matter and limit some parameters of others. The spin-independent direct detection gives strong constraints for some operators, while the indirect detection of antiproton-to-proton data can be more sensitive than direct detection or relic density for light dark matter (whose mass less than 70 GeV) in some cases. The constraints on some operators for spin-3/2 dark matter are shown to be similar to those on their analogous operators for Dirac fermionic dark matter. There are still some operators not sensitive to the current dark matter direct and indirect search experiments.

In this work we study the constraints on the dark matter interaction with the standard model particles, from the observations of dark matter relic density, the direct detection experiments of CDMS and XENON, and the indirect detection of the antiproton-to-proton ratio by PAMELA. A model independent way is adopted in the study by constructing the effective interaction operators between dark matter and standard model particles. The most general 4-fermion operators are investigated. We find that the constraints from different observations are complementary with each other. Especially the spin independent scattering gives very strong constraints for corresponding operators. In some cases the indirect detection of antiproton-to-proton data can actually be more sensitive than the direct detection or relic density for light dark matter (less than 70 GeV).

We investigate whether it is possible that viable microbes could have been transported to Earth from the planets in extra-solar systems by means of natural vehicles such as ejecta expelled by comet or asteroid impacts on such planets. The probabilities of close encounters with other solar systems are taken into account as well as the limitations of bacterial survival times inside ejecta in space, caused by radiation and DNA decay. The conclusion is that no potentially DNA/RNA life-carrying ejecta from another solar system in the general Galactic star field landed on Earth before life already existed on Earth, not even if microbial survival time in space is as long as tens of millions of years. However, if the Sun formed initially as a part of a star cluster, as is commonly assumed, we cannot rule out the possibility of transfer of life from one of the sister systems to us. Likewise, there is a possibility that some extra-solar planets carry life that originated in our solar system. It will be of great interest to identify the members of the Sun's birth cluster of stars and study them for evidence for planets and life on the planets. The former step may be accomplished by the GAIA mission, the latter step by the SIM and DARWIN missions. Therefore it may not be too long until we have experimental knowledge on the question whether the natural transfer of life from one solar system to another has actually taken place.

The inner two planets around the 55 Cancri were found to be trapped in the 3:1 mean motion resonance. In this paper, we study the dynamics of this extra-solar planetary system. Our numerical investigation confirms the existence of the 3:1 resonance and implies a complex orbital motion. Different stable motion types, with and without the apsidal corotation, are found. Due to the high eccentricities in this system, we apply a semi-analytical method based on a new expansion of the Hamiltonian of the planar three-body problem in the discussion. We analyse the occurrence of the apsidal corotation in this mean motion resonance and its influence on the stability of the system.