Time series data of celestial objects are commonly used to study valuable and unexpected objects such as extrasolar planets and supernova in time domain astronomy. Due to the rapid growth of data volume, traditional manual methods are becoming extremely hard and infeasible for continuously analyzing accumulated observation data. To meet such demands, we designed and implemented a special tool named AstroCatR that can efficiently and flexibly reconstruct time series data from large-scale astronomical catalogues. AstroCatR can load original catalogue data from Flexible Image Transport System (FITS) files or databases, match each item to determine which object it belongs to, and finally produce time series datasets. To support the high-performance parallel processing of large-scale datasets, AstroCatR uses the extract-transform-load (ETL) preprocessing module to create sky zone files and balance the workload. The matching module uses the overlapped indexing method and an in-memory reference table to improve accuracy and performance. The output of AstroCatR can be stored in CSV files or be transformed other into formats as needed. Simultaneously, the module-based software architecture ensures the flexibility and scalability of AstroCatR. We evaluated AstroCatR with actual observation data from The three Antarctic Survey Telescopes (AST3). The experiments demonstrate that AstroCatR can efficiently and flexibly reconstruct all time series data by setting relevant parameters and configuration files. Furthermore, the tool is approximately 3X faster than methods using relational database management systems at matching massive catalogues.

Astronomical observation data require long-term preservation, and the rapid accumulation of observation data makes it necessary to consider the cost of long-term archive storage. In addition to low-speed disk-based online storage, optical disk or tape-based offline storage can be used to save costs. However, for astronomical research that requires historical data (particularly time-domain astronomy), the performance and energy consumption of data-accessing techniques cause problems because the requested data (which are organized according to observation time) may be located across multiple storage devices. In this study, we design and develop a tool referred to as AstroLayout to redistribute the observation data using spatial aggregation. The core algorithm uses graph partitioning to generate an optimized data placement according to the original observation data statistics and the target storage system. For the given observation data, AstroLayout can copy the long-term archive in the target storage system in accordance with this placement. An efficiency evaluation shows that AstroLayout can reduce the number of devices activated when responding to data-access requests in time-domain astronomy research. In addition to improving the performance of data-accessing techniques, AstroLayout can also reduce the storage systems power consumption. For enhanced adaptability, it supports storage systems of any media, including optical disks, tapes, and hard disks.

In the era of big data astronomy, next generation telescopes and large sky surveys produce data sets at the TB or even PB level. Due to their large data volumes, these astronomical data sets are extremely difficult to transfer and analyze using personal computers or small clusters. In order to offer better access to data, data centers now generally provide online science platforms that enable analysis close to the data. The Chinese Virtual Observatory (China-VO) is one of the member projects in the International Virtual Observatory Alliance and it is dedicated to providing a research and education environment where globally distributed astronomy archives are simple to find, access, and interoperate. In this study, we summarize highlights of the work conducted at the China-VO, as well the experiences and lessons learned during the full life-cycle management of astronomical data. Finally, We discuss the challenges and future trends for astronomical science platforms.

In this paper, we investigate the magnetohydrodynamical structure of a jet powered by a spinning black hole, where electromagnetic fields and fluid motion are governed by the Grad-Shafranov equation and the Bernoulli equation, respectively. Assuming steady and axisymmetric jet structure, the global solution is uniquely determined with prescribed plasma loading into the jet and the poloidal shape of the outmost magnetic field line. We apply this model to the jet in the center of nearby radio galaxy M87, and we find it can naturally explain the slow flow acceleration and the flow velocity stratification within $10^5$ gravitational radii from the central black hole. In particular, we find the extremal black hole spin is disfavored by the flow velocity measurements, if the plasma loading to the jet is dominated by the electron/positron pair production at the jet base.

We study the angular momentum transport inside the hot Jupiters under the the influences of gravitational and thermal forcing. Due to the strong stellar irradiation, radiative region develops on top of the convective region. Internal gravity waves are launched at the radiative-convective boundaries (RCBs). The thermal response is dynamical and plays an important role in the angular momentum transport. By separating the gravitational and thermal forcing terms, we identify the thermal effects for increasing the angular momentum transport. For the low frequency (in the co-rotating frame with planets) prograde (retrograde) tidal frequency, the angular momentum flux is positive (negative). The tidal interactions tends to drive the planet to the synchronous state. We find that the angular momentum transport associated with the internal gravity wave is very sensitive to relative position between the RCB and the penetration depth of the thermal forcing. If the RCB is in the vicinity of the thermal forcing penetration depth, even with small amplitude thermal forcing, the thermally driven angular momentum flux could be much larger than the flux induced by gravitational forcing. The thermally enhanced torque could drive the planet to the synchronous state in as short as a few $10^4$ years.

Branches of cosmic inflationary models, such as slow-roll inflation, predict a background of primordial gravitational waves that imprints a unique odd-parity B-mode pattern in the Cosmic Microwave Background (CMB) at amplitudes that are within experimental reach. The BICEP/Keck (BK) experiment targets this primordial signature, the amplitude of which is parameterized by the tensor-to-scalar ratio r, by observing the polarized microwave sky through the exceptionally clean and stable atmosphere at the South Pole. B-mode measurements require an instrument with exquisite sensitivity, tight control of systematics, and wide frequency coverage to disentangle the primordial signal from the Galactic foregrounds. BICEP Array represents the most recent stage of the BK program, and comprises four BICEP3-class receivers observing at 30/40, 95, 150 and 220/270 GHz. The 30/40 GHz receiver will be deployed at the South Pole during the 2019/2020 austral summer. After 3 full years of observations with 30,000+ detectors, BICEP Array will measure primordial gravitational waves to a precision $\sigma(r)$ between 0.002 and 0.004, depending on foreground complexity and the degree of lensing removal. In this paper we give an overview of the instrument, highlighting the design features in terms of cryogenics, magnetic shielding, detectors and readout architecture as well as reporting on the integration and tests that are ongoing with the first receiver at 30/40 GHz.

Families of cosmic inflation models predict a primordial gravitational-wave background that imprints B-mode polarization pattern in the Cosmic Microwave Background (CMB). High sensitivity instruments with wide frequency coverage and well-controlled systematic errors are needed to constrain the faint B-mode amplitude. We have developed antenna-coupled Transition Edge Sensor (TES) arrays for high-sensitivity polarized CMB observations over a wide range of millimeter-wave bands. BICEP Array, the latest phase of the BICEP/Keck experiment series, is a multi-receiver experiment designed to search for inflationary B-mode polarization to a precision $\sigma$(r) between 0.002 and 0.004 after 3 full years of observations, depending on foreground complexity and the degree of lensing removal. We describe the electromagnetic design and measured performance of BICEP Array low-frequency 40-GHz detector, their packaging in focal plane modules, and optical characterization including efficiency and beam matching between polarization pairs. We summarize the design and simulated optical performance, including an approach to improve the optical efficiency due to mismatch losses. We report the measured beam maps for a new broad-band corrugation design to minimize beam differential ellipticity between polarization pairs caused by interactions with the module housing frame, which helps minimize polarized beam mismatch that converts CMB temperature to polarization ($T \rightarrow P$) anisotropy in CMB maps.

The BICEP/Keck experiment (BK) is a series of small-aperture refracting telescopes observing degree-scale Cosmic Microwave Background (CMB) polarization from the South Pole in search of a primordial $B$-mode signature. This $B$-mode signal arises from primordial gravitational waves interacting with the CMB, and has amplitude parametrized by the tensor-to-scalar ratio $r$. Since 2016, BICEP3 and the Keck Array have been observing with 4800 total antenna-coupled transition-edge sensor detectors, with frequency bands spanning 95, 150, 220, and 270 GHz. Here we present the optical performance of these receivers from 2016 to 2019, including far-field beams measured in situ with an improved chopped thermal source and instrument spectral response measured with a field-deployable Fourier Transform Spectrometer. As a pair differencing experiment, an important systematic that must be controlled is the differential beam response between the co-located, orthogonally polarized detectors. We generate per-detector far-field beam maps and the corresponding differential beam mismatch that is used to estimate the temperature-to-polarization leakage in our CMB maps and to give feedback on detector and optics fabrication. The differential beam parameters presented here were estimated using improved low-level beam map analysis techniques, including efficient removal of non-Gaussian noise as well as improved spatial masking. These techniques help minimize systematic uncertainty in the beam analysis, with the goal of constraining the bias on $r$ induced by temperature-to-polarization leakage to be subdominant to the statistical uncertainty. This is essential as we progress to higher detector counts in the next generation of CMB experiments.

The BICEP/Keck (BK) experiment aims to detect the imprint of primordial gravitational waves in the Cosmic Microwave Background polarization, which would be direct evidence of the inflation theory. While the tensor-to-scalar ratio has been constrained to be r_0.05 < 0.06 at 95% c.l., further improvements on this upper limit are hindered by polarized Galactic foreground emissions and removal of gravitational lensing polarization. The 30/40 GHz receiver of the BICEP Array (BA) will deploy at the end of 2019 and will constrain the synchrotron foreground with unprecedented accuracy within the BK sky patch. We will show the design of the 30/40 GHz detectors and test results summarizing its performance. The low optical and atmospheric loading at these frequencies requires our TES detectors to have low saturation power in order to be photon-noise dominated. To realize the low thermal conductivity required from a 250 mK base temperature, we developed new bolometer leg designs. We will present the relevant measured detector parameters: G, Tc, Rn, Psat , and spectral bands, and noise spectra. We achieved a per bolometer NEP including all noise components of 2.07E-17 W/sqrt(Hz), including an anticipated photon noise level 1.54E-17 W/sqrt(Hz).

Radio frequency interference (RFI) detection and excision are key steps in the data-processing pipeline of the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Because of its high sensitivity and large data rate, FAST requires more accurate and efficient RFI flagging methods than its counterparts. In the last decades, approaches based upon artificial intelligence (AI), such as codes using convolutional neural networks (CNNs), have been proposed to identify RFI more reliably and efficiently. However, RFI flagging of FAST data with such methods has often proved to be erroneous, with further manual inspections required. In addition, network construction as well as preparation of training data sets for effective RFI flagging has imposed significant additional workloads. Therefore, rapid deployment and adjustment of AI approaches for different observations is impractical to implement with existing algorithms. To overcome such problems, we propose a model called RFI-Net. With the input of raw data without any processing, RFI-Net can detect RFI automatically, producing corresponding masks without any alteration of the original data. Experiments with RFI-Net using simulated astronomical data show that our model has outperformed existing methods in terms of both precision and recall. Besides, compared with other models, our method can obtain the same relative accuracy with fewer training data, thus reducing the effort and time required to prepare the training data set. Further, the training process of RFI-Net can be accelerated, with overfittings being minimized, compared with other CNN codes. The performance of RFI-Net has also been evaluated with observing data obtained by FAST and the Bleien Observatory. Our results demonstrate the ability of RFI-Net to accurately identify RFI with fine-grained, high-precision masks that required no further modification.

We describe an on-sky demonstration of a microwave-multiplexing readout system in one of the receivers of the Keck Array, a polarimetry experiment observing the cosmic microwave background at the South Pole. During the austral summer of 2018-2019, we replaced the time-division multiplexing readout system with microwave-multiplexing components including superconducting microwave resonators coupled to radio-frequency superconducting quantum interference devices at the sub-Kelvin focal plane, coaxial-cable plumbing and amplification between room temperature and the cold stages, and a SLAC Microresonator Radio Frequency system for the warm electronics. In the range 5-6 GHz, a single coaxial cable reads out 528 channels. The readout system is coupled to transition-edge sensors, which are in turn coupled to 150-GHz slot-dipole phased-array antennas. Observations began in April 2019, and we report here on an initial characterization of the system performance.

Large arrays of cryogenic detectors, including transition-edge sensors (TESs) or magnetic micro-calorimeters (MMCs), are needed for future experiments across a wide range of applications. Complexities in integration and cryogenic wiring have driven efforts to develop cryogenic readout technologies with large multiplexing factors while maintaining minimal readout noise. One such example is the microwave SQUID multiplexer ($\mu$mux), which couples an incoming TES or magnetic calorimeter signal to a unique GHz-frequency resonance that is modulated in frequency. Here, we present a hybrid scheme combining the microwave SQUID multiplexer with code division multiplexing: the impedance-modulated code-division multiplexer (Z-CDM), which may enable an order of magnitude increase in multiplexing factor particularly for low-bandwidth signal applications.

The MAJORANA DEMONSTRATOR is sensitive to rare events near its energy threshold, including bosonic dark matter, solar axions, and lightly ionizing particles. In this analysis, a novel training set of low energy small-angle Compton scatter events is used to determine the efficiency of pulse shape analysis cuts, and we present updated bosonic dark matter and solar axion results from an 11.17 kg-y dataset using a 5 keV analysis threshold.

The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4.

In this paper, we construct a framework for investigating magnetohydrodynamical jet structure of spinning black holes (BHs), where electromagnetic fields and fluid motion are governed by the Grad-Shafranov equation and the Bernoulli equation, respectively. Assuming steady and axisymmetric jet structure, we can self-consistently obtain electromagnetic fields, fluid energy density and velocity within the jet, given proper plasma loading and boundary conditions. Specifically, we structure the two coupled governing equations as two eigenvalue problems, and develop full numerical techniques for solving them. As an example, we explicitly solve the governing equations for the split monopole magnetic field configuration and simplified plasma loading on the stagnation surface where the poloidal fluid velocity vanishes. As expected, we find the rotation of magnetic field lines is dragged down by fluid inertia, and the fluid as a whole does not contribute to energy extraction from the central BH, i.e., the magnetic Penrose process is not working. However, if we decompose the charged fluid as two oppositely charged components, we find the magnetic Penrose process does work for one of the two components when the plasma loading is low enough.

The Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands: 27, 39, 93, 145, 225 and 280 GHz. The initial configuration of SO will have three small-aperture 0.5-m telescopes (SATs) and one large-aperture 6-m telescope (LAT), with a total of 60,000 cryogenic bolometers. Our key science goals are to characterize the primordial perturbations, measure the number of relativistic species and the mass of neutrinos, test for deviations from a cosmological constant, improve our understanding of galaxy evolution, and constrain the duration of reionization. The SATs will target the largest angular scales observable from Chile, mapping ~10% of the sky to a white noise level of 2 $\mu$K-arcmin in combined 93 and 145 GHz bands, to measure the primordial tensor-to-scalar ratio, $r$, at a target level of $\sigma(r)=0.003$. The LAT will map ~40% of the sky at arcminute angular resolution to an expected white noise level of 6 $\mu$K-arcmin in combined 93 and 145 GHz bands, overlapping with the majority of the LSST sky region and partially with DESI. With up to an order of magnitude lower polarization noise than maps from the Planck satellite, the high-resolution sky maps will constrain cosmological parameters derived from the damping tail, gravitational lensing of the microwave background, the primordial bispectrum, and the thermal and kinematic Sunyaev-Zel'dovich effects, and will aid in delensing the large-angle polarization signal to measure the tensor-to-scalar ratio. The survey will also provide a legacy catalog of 16,000 galaxy clusters and more than 20,000 extragalactic sources.

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 present an analysis of meteorological data from the second generation of the Kunlun Automated Weather Station (KLAWS-2G) at Dome A, Antarctica during 2015 and 2016. We find that a strong temperature inversion exists for all the elevations up to 14 m that KLAWS-2G can reach, and lasts for more than 10 hours for 50% or more of the time when temperature inversion occurs. The average wind speeds at 4 m elevation are 4.2 m/s and 3.8 m/s during 2015 and 2016, respectively. The strong temperature inversion and moderate wind speed lead to a shallow turbulent boundary layer height at Dome A. By analyzing the temperature and wind shear profiles, we note telescopes should be elevated by at least 8 m above the ice. We also find that the duration of temperature inversions, and the wind speed, vary considerably from year to year. Therefore, long-term and continuous data are still needed for the site survey at Dome A.

The next generation of cryogenic CMB and submillimeter cameras under development require densely instrumented sensor arrays to meet their science goals. The readout of large numbers ($\sim$10,000--100,000 per camera) of sub-Kelvin sensors, for instance as proposed for the CMB-S4 experiment, will require substantial improvements in cold and warm readout techniques. To reduce the readout cost per sensor and integration complexity, efforts are presently focused on achieving higher multiplexing density while maintaining readout noise subdominant to intrinsic detector noise. Highly-multiplexed cold readout technologies in active development include Microwave Kinetic Inductance Sensors (MKIDs) and microwave rf-SQUIDs. Both exploit the high quality factors of superconducting microwave resonators to densely channelize sub-Kelvin sensors into the bandwidth of a microwave transmission line. We present advancements in the development of a new warm readout system for microwave SQUID multiplexing, the SLAC Superconducting Microresonator RF electronics, or SMuRF. The SMuRF system is unique in its ability to track each tone, minimizing the total RF power required to read out each resonator, thereby significantly reducing the linearity requirements on the cold and warm readout. Here, we present measurements of the readout noise and linearity of the first full SMuRF system, including a demonstration of closed-loop tone tracking on a 528 channel cryogenic microwave SQUID multiplexer. SMuRF is being explored as a potential readout solution for several future CMB projects including Simons Observatory, BICEP Array, CCAT-prime, Ali-CPT, and CMB-S4. Parallel development of the platform is underway to adapt SMuRF to read out both MKID and fast X-ray TES calorimeter arrays.

The structure of force-free, steady and axisymmetric magnetosphere of a neutron star (NS) is governed by the Grad-Shafranov (GS) equation, which is a second-order differential equation but degrades to first-order on the light surface (LS). The key to numerically solving the GS equation is to enable magnetic field lines smoothly cross the LS, and crossing a straight LS in flat spacetime has been a well-studied problem. But the numerical algorithm implementation becomes complicate in the presence of a bent LS, e.g. in curved spacetime, since there is no suitable computation grid adapted to it. We propose to circumvent this grid-LS mismatch problem by introducing a new coordinate frame designed such that the LS in it is a straight line. As an application, we investigate the general relativistic (GR) effect in magnetosphere structure of rapidly rotating pulsars in detail, where the LS is bent towards the central NS. We split the GR effect into two parts, curvature and frame-dragging; measure each of them and examine their dependence on the NS mass and the angular velocity for pulsars embedded in aligned dipole and multipole magnetic fields. Qualitatively speaking, we find that the curvature effect compactifies the magnetic field lines near the NS, therefore reduces the open magnetic flux and the Poynting luminosity, while the frame-dragging effect contributes a minor part in shaping the magnetosphere structure but plays a role in enhancing the spacelike current generation.

The Greenland Telescope project has recently participated in an experiment to image the supermassive black hole shadow at the center of M87 using Very Long Baseline Interferometry technique in April of 2018. The antenna consists of the 12-m ALMA North American prototype antenna that was modified to support two auxiliary side containers and to withstand an extremely cold environment. The telescope is currently at Thule Air Base in Greenland with the long-term goal to move the telescope over the Greenland ice sheet to Summit Station. The GLT currently has a single cryostat which houses three dual polarization receivers that cover 84-96 GHz, 213-243 GHz and 271-377 GHz bands. A hydrogen maser frequency source in conjunction with high frequency synthesizers are used to generate the local oscillator references for the receivers. The intermediate frequency outputs of each receiver cover 4-8 GHz and are heterodyned to baseband for digitization within a set of ROACH-2 units then formatted for recording onto Mark-6 data recorders. A separate set of ROACH-2 units operating in parallel provides the function of auto-correlation for real-time spectral analysis. Due to the stringent instrumental stability requirements for interferometry a diagnostic test system was incorporated into the design. Tying all of the above equipment together is the fiber optic system designed to operate in a low temperature environment and scalable to accommodate a larger distance between the control module and telescope for Summit Station. A report on the progress of the above electronics instrumentation system will be provided.

The \textsc{Majorana Demonstrator} is an ultra low-background experiment searching for neutrinoless double-beta decay in $^{76}$Ge. The heavily shielded array of germanium detectors, placed nearly a mile underground at the Sanford Underground Research Facility in Lead, South Dakota, also allows searches for new exotic physics. Free, relativistic, lightly-ionizing particles with electrical charges less than $e$ are forbidden by the standard model but predicted by some of its extensions. If such particles exist, they might be detected in the \textsc{Majorana Demonstrator} by searching for multiple- detector events with individual-detector energy depositions down to 1 keV. This search is background free and no candidate events have been found in 285 days of data taking. New direct-detection limits are set for the flux of lightly ionizing particles for charges as low as $e$/1000.

The three Antarctic Survey Telescopes (AST3) aim to carry out time domain imaging survey at Dome A, Antarctica. The first of the three telescopes (AST3-1) was successfully deployed on January 2012. AST3-1 is a 500\,mm aperture modified Schmidt telescope with a 680\,mm diameter primary mirror. AST3-1 is equipped with a SDSS $i$ filter and a 10k $\times$ 10k frame transfer CCD camera, reduced to 5k $\times$ 10k by electronic shuttering, resulting in a 4.3 deg$^2$ field-of-view. To verify the capability of AST3-1 for a variety of science goals, extensive commissioning was carried out between March and May 2012. The commissioning included a survey covering 2000 deg$^2$ as well as the entire Large and Small Magellanic Clouds. Frequent repeated images were made of the center of the Large Magellanic Cloud, a selected exoplanet transit field, and fields including some Wolf-Rayet stars. Here we present the data reduction and photometric measurements of the point sources observed by AST3-1. We have achieved a survey depth of 19.3\,mag in 60 s exposures with 5\,mmag precision in the light curves of bright stars. The facility achieves sub-mmag photometric precision under stable survey conditions, approaching its photon noise limit. These results demonstrate that AST3-1 at Dome A is extraordinarily competitive in time-domain astronomy, including both quick searches for faint transients and the detection of tiny transit signals.

The gamma-ray burst GRB 170817A associated with GW170817 is subluminous and subenergetic compared with other typical short GRBs. It may be due to a relativistic jet viewed off-axis, or a structured jet, or cocoon emission. Giant flares from magnetars may possibly be ruled out. However, the luminosity and energetics of GRB 170817A is coincident with that of magnetar giant flares. After the coalescence of the binary neutron star, a hypermassive neutron star may be formed. The hypermassive neutron star may have magnetar-strength magnetic field. During the collapse of the hypermassive neutron star, the magnetic field energy will also be released. This giant-flare-like event may explain the the luminosity and energetics of GRB 170817A. Bursts with similar luminosity and energetics are expected in future neutron star-neutron star or neutron star-black hole mergers.

With many large science equipment constructing and putting into use, astronomy has stepped into the big data era. The new method and infrastructure of big data processing has become a new requirement of many astronomers. Cloud computing, Map/Reduce, Hadoop, Spark, etc. many new technology has sprung up in recent years. Comparing to the high performance computing(HPC), Data is the center of these new technology. So, a new computing architecture infrastructure is necessary, which can be shared by both HPC and big data processing. Based on Astronomy Cloud project of Chinese Virtual Observatory (China-VO), we have made much efforts to optimize the designation of the hybrid computing platform. which include the hardware architecture, cluster management, Job and Resource scheduling.

We have analyzed the multipolar magnetic field structure variation at neutron star surface by means of the catastrophic eruption model, and find that the variation of the geometry of multipolar fields on the magnetar surface could result in the catastrophic rearrangement of the magnetosphere, which provides certain physical mechanism for the outburst of giant flares. The magnetospheric model we adopted consists of two assumptions: a helically twisted flux rope is suspended in an ideal force-free magnetosphere around the magnetar, and a current sheet emerges during the flux rope evolution. Magnetic energy accumulates during the flux rope's gradual evolution along with the variation of magnetar surface magnetic structure before the eruption. The two typical behaviors, either state transition or catastrophic escape, would take place once the flux rope loses equilibrium, thus tremendous accumulated energy is radiated. We have investigated the equilibrium state of the flux rope and the energy release affected by different multipolar structures, and find structures that could trigger violent eruption and provide the radiation approximately 0.5$\%$ of the total magnetic energy during the giant flare outburst. Our results provide certain multipolar structures of the neutron star's magnetic field with an energy release percentage 0.42$\%$ in the state transition and 0.51$\%$ in the catastrophic escape case, which are sufficient for the previously reported energy release from SGR 1806-20 giant flares.

The Kepler observations indicate that many exoplanets are super-Earths, which brings about a puzzle for the core-accretion scenario. Since observed super-Earths are in the range of critical mass, they would accrete gas efficiently and become gas giants. Theoretically, super-Earths are predicted to be rare in the core-accretion framework. To resolve this contradiction, we propose that the tidally-forced turbulent diffusion may affect the heat transport inside the planet. Thermal feedback induced by turbulent diffusion is investigated. We find that the tidally-forced turbulence would generate pseudo-adiabatic regions within radiative zones, which pushes the radiative-convective boundaries (RCBs) inwards. This would decrease the cooling luminosity and enhance the Kelvin-Helmholtz (KH) timescale. For a given lifetime of protoplanetary disks (PPDs), there exists a critical threshold for the turbulent diffusivity, $\nu_{\rm critical}$. If $\nu_{\rm turb}>\nu_{\rm critical} $, the KH timescale is longer than the disk lifetime and the planet would become a super-Earth rather than a gas giant. We find that even a small value of turbulent diffusion has influential effects on evolutions of super-Earths. $\nu_{\rm critical}$ increases with the core mass. We further ascertain that, within the minimum mass extrasolar nebula (MMEN), $\nu_{\rm critical}$ increases with the semi-major axis. This may explain the feature that super-Earths are common in inner PPD regions, while gas giants are common in the outer PPD regions. The predicted envelope mass fraction (EMF) is not fully consistent with observations. We discuss physical processes, such as late core assembly and mass loss mechanisms, that may be operating during super-Earth formation.

We present a fully integrated InGaAs/InP negative feedback avalanche diode (NFAD) based free-running single-photon detector (SPD) designed for accurate lidar applications. A free-piston Stirling cooler is used to cool down the NFAD with a large temperature range, and an active hold-off circuit implemented in a field programmable gate array is applied to further suppress the afterpulsing contribution. The key parameters of the free-running SPD including photon detection efficiency (PDE), dark count rate (DCR), afterpulse probability, and maximum count rate (MCR) are dedicatedly optimized for lidar application in practice. We then perform a field experiment using a Mie lidar system with 20 kHz pulse repetition frequency to compare the performance between the free-running InGaAs/InP SPD and a commercial superconducting nanowire single-photon detector (SNSPD). Our detector exhibits good performance with 1.6 Mcps MCR (0.6 {\mu}s hold-off time), 10% PDE, 950 cps DCR, and 18% afterpulse probability over 50 {\mu}s period. Such performance is worse than the SNSPD with 60% PDE and 300 cps DCR. However, after performing a specific algorithm that we have developed for afterpulse and count rate corrections, the lidar system performance in terms of range-corrected signal (Pr2) distribution using our SPD agrees very well with the result using the SNSPD, with only a relative error of ~2%. Due to the advantages of low-cost and small size of InGaAs/InP NFADs, such detector provides a practical solution for accurate lidar applications.

The structure of steady axisymmetric force-free magnetosphere of a Kerr black hole (BH) is governed by a second-order partial differential equation of $A_\phi$ depending on two "free" functions $\Omega(A_\phi)$ and $I(A_\phi)$, where $A_\phi$ is the $\phi$ component of the vector potential of the electromagnetic field, $\Omega$ is the angular velocity of the magnetic field lines and $I$ is the poloidal electric current. In this paper, we investigate the solution uniqueness. Taking asymptotically uniform field as an example, analytic studies imply that there are infinitely many solutions approaching uniform field at infinity, while only a unique one is found in general relativistic magnetohydrodynamic simulations. To settle down the disagreement, we reinvestigate the structure of the governing equation and numerically solve it with given constraint condition and boundary condition. We find that the constraint condition (field lines smoothly crossing the light surface (LS)) and boundary conditions at horizon and at infinity are connected via radiation conditions at horizon and at infinity, rather than being independent. With appropriate constraint condition and boundary condition, we numerically solve the governing equation and find a unique solution. Contrary to naive expectation, our numerical solution yields a discontinuity in the angular velocity of the field lines and a current sheet along the last field line crossing the event horizon. We also briefly discuss the applicability of the perturbation approach to solving the governing equation.

Virtual Observatory (VO) is a data-intensively online astronomical research and education environment, which takes advantages of advanced information technologies to achieve seamless and global access to astronomical information. AstroCloud is a cyber-infrastructure for astronomy research initiated by Chinese Virtual Observatory (China-VO) project, and also a kind of physical distributed platform which integrates lots of tasks such as telescope access proposal management, data archiving, data quality control, data release and open access, cloud based data processing and analysis. It consists of five application channels, i.e. observation, data, tools, cloud and public and is acting as a full lifecycle management system and gateway for astronomical data and telescopes. Physically, the platform is hosted in six cities currently, i.e. Beijing, Nanjing, Shanghai, Kunming, Lijiang and Urumqi, and serving more than 17 thousand users. Achievements from international Virtual Observatories and Cloud Computing are adopted heavily. In the paper, backgrounds of the project, architecture, Cloud Computing environment, key features of the system, current status and future plans are introduced.

We propose that the magnetosphere reconfiguration induced by magnetic twists in the closed field line region can account for the mode-switching of intermittent pulsars. We carefully investigate the properties of axisymmetric force-free pulsar magnetospheres with magnetic twists in closed field line region around the polar caps. The magnetosphere with twisted closed lines leads to enhanced spin-down rates. The enhancement in spin-down rate depends on the size of region with twisted closed lines. Typically, it is increased by a factor of $\sim2$, which is consistent with the intermittent pulsars' spin down behavior during the `off' and `on' states. We find there is a threshold of maximal twist angle $\Delta\phi_{\rm thres}\sim1$. The magnetosphere is stable only if the closed line twist angle is less than $\Delta\phi_{\rm thres}$. Beyond this value, the magnetosphere becomes unstable and gets untwisted. The spin-down rate would reduce to its off-state value. The quasi-periodicity in spin-down rate change can be explained by long-term activities in star's crust and the untwisting induced by MHD instability. The estimated duration time of on-state is about one week, consistent with observations. Due to the MHD instability, there exists an upper limit for the spin down ratio ($f\sim3$) between the on-state and the off-state, if the Y-point remains at the light cylinder.

As the cyber-infrastructure for Astronomical research from Chinese Virtual Observatory (China-VO) project, AstroCloud has been archived solid progresses during the last one year. Proposal management system and data access system are re-designed. Several new sub-systems are developed, including China-VO PaperData, AstroCloud Statics and Public channel. More data sets and application environments are integrated into the platform. LAMOST DR1, the largest astronomical spectrum archive was released to the public using the platform. The latest progresses will be introduced.

We reinvestigate the structure of a steady axisymmetic force-free magnetosphere around a Kerr black hole (BH). The BH magnetosphere structure is governed by a second-order differential equation of $A_\phi$ depending on two `free' functions $\Omega$ and $I$, where $A_\phi$ is the $\phi$ component of the vector potential of the electromagnetic field, $\Omega$ is the angular velocity of the magnetic field lines and $I$ is the poloidal electric current. While the two functions $\Omega$ and $I$ are not arbitrarily given, which need to be self-consistently determined along with the differential equation. Based on the perturbation approach we proposed in paper I \citep{Pan2015a}, in this paper, we self-consistently sort out two boundary conditions governing $\Omega$ and $I$, and interpret these conditions mathematically and physically. Making use of the boundary conditions, we prove that all magnetic field lines crossing the infinite-redshift surface also penetrate the event horizon. Furthermore, we argue that the BH Meissner effect does not work in force-free magnetosphere due to the perfect conductivity.

Blandford-Znajek (BZ) mechanism describes a process extracting rotation energy from a spinning black hole (BH) via magnetic field lines penetrating the event horizon of central BH. In this paper, we present a perturbation approach to study force-free jets launched by the BZ mechanism, and its two immediate applications: (1) we present a high-order split monopole perturbation solution to the BZ mechanism, which accurately pins down the energy extraction rate $\dot E$ and well describes the structure of BH magnetosphere for all range of BH spins ($0\leq a\leq 1$); (2) the approach yields an exact constraint for the monopole field configuration in the Kerr spacetime, $I = \Omega (1-A_\phi^2)$, where $A_\phi$ is the $\phi-$component of the vector potential of electromagnetic field, $\Omega$ is the angular velocity of magnetic field lines and $I$ is the poloidal electric current. The constraint is of particular importance to benchmark the accuracy of numerical simulations.

In this paper, we propose a general method for perturtative solutions to Blandford-Znajek mechanism. Instead of solving the nonlinear Grad-Shafranov equation directly, we introduce an alternative way to determine relevant physical quantities based on the horizon boundary condition and the convergence requirement. Both the angular velocity $\Omega$ of magnetic field lines, the toroidal magnetic field $B^\phi$ and the total electric current $I$ are self-consistently specified according to our method. As an example, stationary axisymmetric and force-free jet models around rotating black holes are self-consistently constructed according to the method we proposed. This jet solution distinguishes itself from prior known analytic solutions in that it is highly collimated and asymptotically approaches a magnetic cylinder. This jet solution is helically twisted, since toroidal magnetic field is generated when the black hole spin is taken into account. For a given magnetic flux threading the black hole, the jet power and energy extraction rate of the collimated jet are compared with previous solutions. We find that our new solution agrees better with current state-of-the-art numerical simulation results. Some interesting properties of the collimated jet and effects of field line rotation on the jet stability are also briefly discussed.

We study the effect of large-scale magnetic fields on the non-axisymmetric inertial-acoustic modes (also called p-modes) trapped in the innermost regions of accretion discs around black holes (BHs). These global modes could provide an explanation for the high-frequency quasi-periodic oscillations (HFQPOs) observed in BH X-ray binaries. There may be observational evidence for the presence of such large-scale magnetic fields in the disks since episodic jets are observed in the same spectral state when HFQPOs are detected. We find that a large-scale poloidal magnetic field can enhance the corotational instability and increase the growth rate of the purely hydrodynamic overstable p-modes. In addition, we show that the frequencies of these overstable p-modes could be further reduced by such magnetic fields, making them agree better with observations.

We study force-free magnetospheres in the Blandford-Znajek process from rapidly rotating black holes by adopting the near-horizon geometry of near-extreme Kerr black holes (near-NHEK). It is shown that the Znajek regularity condition on the horizon can be directly derived from the resulting stream equation. In terms of the condition, we split the full stream equation into two separate equations. Approximate solutions around the rotation axis are derived. They are found to be consistent with previous solutions obtained in the asymptotic region. The solutions indicate energy and angular-momentum extraction from the hole.

The Blandford-Znajek (BZ) mechanism describes a physical process for the energy extraction from a spinning black hole (BH), which is believed to power a great variety of astrophysical sources, such as active galactic nuclei (AGNs) and Gamma ray bursts (GRBs). The only known analytic solution to the BZ mechanism is a split monopole perturbation solution up to $O(a^2)$, where $a$ is the spin parameter of a Kerr black hole. In this paper, we extend the monopole solution to higher order $\sim O(a^4)$. We carefully investigate the structure of the BH magnetosphere, including the angular velocity of magnetic field lines $\Omega$, the toroidal magnetic field $B^\phi$ as well as the poloidal electric current $I$. In addition, the relevant energy extraction rate $\dot E$ and the stability of this high-order monopole perturbation solution are also examined.

AstroCloud is a cyber-Infrastructure for Astronomy Research initiated by Chinese Virtual Observatory (China-VO) under funding support from NDRC (National Development and Reform commission) and CAS (Chinese Academy of Sciences). The ultimate goal of this project is to provide a comprehensive end-to-end astronomy research environment where several independent systems seamlessly collaborate to support the full lifecycle of the modern observational astronomy based on big data, from proposal submission, to data archiving, data release, and to in-situ data analysis and processing. In this paper, the architecture and key designs of the AstroCloud platform are introduced, including data access middleware, access control and security framework, extendible proposal workflow, and system integration mechanism.

AstroCloud is a cyber-Infrastructure for Astronomy Research initiated by Chinese Virtual Observatory (China-VO) under funding support from NDRC (National Development and Reform commission) and CAS (Chinese Academy of Sciences){\url{http://astrocloud.china-vo.org}}\citep{O8-5_Cui_adassxxiv}. To archive the astronomical data in China, we present the implementation of the astronomical data archiving system (ADAS). Data archiving and quality control are the infrastructure for the AstroCloud. Throughout the data of the entire life cycle, data archiving system standardized data, transferring data, logging observational data, archiving ambient data, And storing these data and metadata in database. Quality control covers the whole process and all aspects of data archiving.

Data access and interoperability module connects the observation proposals, data, virtual machines and software. According to the unique identifier of PI (principal investigator), an email address or an internal ID, data can be collected by PI's proposals, or by the search interfaces, e.g. conesearch. Files associated with the searched results could be easily transported to cloud storages, including the storage with virtual machines, or several commercial platforms like Dropbox. Benefitted from the standards of IVOA (International Observatories Alliance), VOTable formatted searching result could be sent to kinds of VO software. Latter endeavor will try to integrate more data and connect archives and some other astronomical resources.

We propose a physical mechanism to explain giant flares and radio afterglows in terms of a magnetospheric model containing both a helically twisted flux rope and a current sheet (CS). With the appearance of CS, we solve a mixed boundary value problem to get the magnetospheric field based on a domain decomposition method. We investigate properties of the equilibrium curve of the flux rope when the CS is present in background multipolar fields. In response to the variations at the magnetar surface, it quasi-statically evolves in stable equilibrium states. The loss of equilibrium occurs at a critical point and, beyond that point, it erupts catastrophically. New features show up when the CS is considered. Especially, we find two kinds of physical behaviors, i.e., catastrophic state transition and catastrophic escape. Magnetic energy would be released during state transitions. The released magnetic energy is sufficient to drive giant flares. The flux rope would go away from the magnetar quasi-statically, which is inconsistent with the radio afterglow. Fortunately, in the latter case, i.e., the catastrophic escape, the flux rope could escape the magnetar and go to infinity in a dynamical way. This is more consistent with radio afterglow observations of giant flares. We find that the minor radius of flux rope has important implications for its eruption. Flux ropes with larger minor radius are more prone to erupt. We stress that the CS provides an ideal place for magnetic reconnection, which would further enhance the energy release during eruptions.

The China Dark Matter Experiment reports results on light WIMP dark matter searches at the China Jinping Underground Laboratory with a germanium detector array with a total mass of 20 g. The physics threshold achieved is 177 eVee ("ee" represents electron equivalent energy) at 50% signal efficiency. With 0.784 kg-days of data, exclusion region on spin-independent coupling with the nucleon is derived, improving over our earlier bounds at WIMP mass less than 4.6 GeV.

A low-background, high-purity germanium detector has been used to search for evidence of low-energy, bremsstrahlung-generated solar axions. An upper bound of $1.36\times 10^{-11}$ $(95% CL)$ is placed on the direct coupling of DFSZ model axions to electrons. The prospects for the sensitivity of the Majorana Demonstrator array of point-contact germanium detectors to solar axions are discussed in the context of the model-independent annual modulation due to the seasonal variation of the earth-sun distance.

We propose a catastrophic eruption model for magnetar's enormous energy release during giant flares, in which a toroidal and helically twisted flux rope is embedded within a force-free magnetosphere. The flux rope stays in stable equilibrium states initially and evolves quasi-statically. Upon the loss of equilibrium point is reached, the flux rope cannot sustain the stable equilibrium states and erupts catastrophically. During the process, the magnetic energy stored in the magnetosphere is rapidly released as the result of destabilization of global magnetic topology. The magnetospheric energy that could be accumulated is of vital importance for the outbursts of magnetars. We carefully establish the fully open fields and partially open fields for various boundary conditions at the magnetar surface and study the relevant energy thresholds. By investigating the magnetic energy accumulated at the critical catastrophic point, we find that it is possible to drive fully open eruptions for dipole dominated background fields. Nevertheless, it is hard to generate fully open magnetic eruptions for multipolar background fields. Given the observational importance of the multipolar magnetic fields in the vicinity of the magnetar surface, it would be worthwhile to explore the possibility of the alternative eruption approach in multipolar background fields. Fortunately, we find that flux ropes may give rise to partially open eruptions in the multipolar fields, which involve only partial opening up of background fields. The energy release fractions are greater for cases with central-arcaded multipoles than those with central-caved multipoles emerged in background fields. Eruptions would fail only when the centrally-caved multipoles become extremely strong.

The China Dark matter Experiment collaboration reports the first experimental limit on WIMP dark matter from 14.6 kg-day of data taken with a 994 g p-type point-contact germanium detector at the China Jinping underground Laboratory where the rock overburden is more than 2400 m. The energy threshold achieved was 400 eVee. According to the 14.6 kg-day live data, we placed the limit of N= 1.75 * 10^{-40} cm^{2} at 90% confidence level on the spin-independent cross-section at WIMP mass of 7 GeV before differentiating bulk signals from the surface backgrounds.

We propose a catastrophic magnetospheric model for magnetar precursors and their successive giant flares. Axisymmetric models of the magnetosphere, which contain both a helically twisted flux rope and a current sheet, are established based on force-free field configurations. In this model, the helically twisted flux rope would lose its equilibrium and erupt abruptly in response to the slow and quasi-static variations at the ultra-strongly magnetized neutron star's surface. In a previous model without current sheets, only one critical point exists in the flux rope equilibrium curve. New features show up in the equilibrium curves for the flux rope when current sheets appear in the magnetosphere. The causal connection between the precursor and the giant flare, as well as the temporary re-entry of the quiescent state between the precursor and the giant flare, can be naturally explained. Magnetic energy would be released during the catastrophic state transitions. The detailed energetics of the model are also discussed. The current sheet created by the catastrophic loss of equilibrium of the flux rope provides an ideal place for magnetic reconnection. We point out the importance of magnetic reconnection for further enhancement of the energy release during eruptions.

The CDEX Collaboration has been established for direct detection of light dark matter particles, using ultra-low energy threshold p-type point-contact germanium detectors, in China JinPing underground Laboratory (CJPL). The first 1 kg point-contact germanium detector with a sub-keV energy threshold has been tested in a passive shielding system located in CJPL. The outputs from both the point-contact p+ electrode and the outside n+ electrode make it possible to scan the lower energy range of less than 1 keV and at the same time to detect the higher energy range up to 3 MeV. The outputs from both p+ and n+ electrode may also provide a more powerful method for signal discrimination for dark matter experiment. Some key parameters, including energy resolution, dead time, decay times of internal X-rays, and system stability, have been tested and measured. The results show that the 1 kg point-contact germanium detector, together with its shielding system and electronics, can run smoothly with good performances. This detector system will be deployed for dark matter search experiments.

Weakly Interacting Massive Particles (WIMPs) are the candidates of dark matter in our universe. Up to now any direct interaction of WIMP with nuclei has not been observed yet. The exclusion limits of the spin-independent cross section of WIMP-nucleon which have been experimentally obtained is about 10^{-7}pb at high mass region and only 10^{-5}pb} at low mass region. China Jin-Ping underground laboratory CJPL is the deepest underground lab in the world and provides a very promising environment for direct observation of dark matter. The China Dark Matter Experiment (CDEX) experiment is going to directly detect the WIMP flux with high sensitivity in the low mass region. Both CJPL and CDEX have achieved a remarkable progress in recent two years. The CDEX employs a point-contact germanium semi-conductor detector PCGe whose detection threshold is less than 300 eV. We report the measurement results of Muon flux, monitoring of radioactivity and Radon concentration carried out in CJPL, as well describe the structure and performance of the 1 kg PCGe detector CDEX-1 and 10kg detector array CDEX-10 including the detectors, electronics, shielding and cooling systems. Finally we discuss the physics goals of the CDEX-1, CDEX-10 and the future CDEX-1T detectors.

The physical origin of high-frequency QPOs (HFQPOs) in black-hole X-ray binaries remains an enigma despite many years of detailed observational studies. Although there exists a number of models for HFQPOs, many of these are simply "notions" or "concepts" without actual calculation derived from fluid or disk physics. Future progress requires a combination of numerical simulations and semi-analytic studies to extract physical insights. We review recent works on global oscillation modes in black-hole accretion disks, and explain how, with the help of general relativistic effects, the energy stored in the disk differential rotation can be pumped into global spiral density modes in the disk, making these modes grow to large amplitudes under certain conditions ("corotational instability"). These modes are robust in the presence of disk magnetic fields and turbulence. The computed oscillation mode frequencies are largely consistent with the observed values for HFQPOs in BH X-ray binaries. The approximate 2:3 frequency ratio is also expected from this model. The connection of HFQPOs with other disk properties (such as production of episodic jets) is also discussed.

We study the effect of large-scale magnetic fields on the non-axisymmetric Rossby wave instability (RWI) in accretion discs. The instability develops around a density bump, which is likely present in the transition region between the active zone and dead zone of protoplanetary discs. Previous works suggest that the vortices resulting from the RWI may facilitate planetesimal formation and angular momentum transport. We consider discs threaded by a large-scale poloidal magnetic field, with a radial field component at the disc surface. Such field configurations may lead to the production of magnetic winds or jets. In general, the magnetic field can affect the RWI even when it is sub-thermal (plasma $\beta\sim 10$). For infinitely thin discs, the instability can be enhanced by about 10 percent. For discs with finite thickness, with a radial gradient of the magnetic field strength, the RWI growth rate can increase significantly (by a factor of $\sim 2$) as the field approaches equipartition ($\beta \sim 1$). Our result suggests that the RWI can continue to operate in discs that produce magnetic winds.

We address a primary question regarding the physical mechanism that triggers the energy release and initiates the onset of eruptions in the magnetar magnetosphere. A self-consistent stationary, axisymmetric model of the magnetar magnetosphere is constructed based on a force-free magnetic field configuration which contains a helically twisted force-free flux rope. Given the complex multipolar magnetic fields at the magnetar surface, we also develop a convenient numerical scheme to solve the GS equation. Depending on the surface magnetic field polarity, there exist two kinds of magnetic field configurations, inverse and normal. For these two kinds of configurations, variations of the flux rope equilibrium height in response to gradual surface physical processes, such as flux injections and crust motions, are carefully examined. We find that equilibrium curves contain two branches, one represents a stable equilibrium branch, the other an unstable equilibrium branch. As a result, the evolution of the system shows a catastrophic behavior: when the magnetar surface magnetic field evolves slowly, the height of flux rope would gradually reach a critical value beyond which stable equilibriums can no longer be maintained. Subsequently the flux rope would lose equilibrium and the gradual quasi-static evolution of the magnetar magnetosphere will be replaced by a fast dynamical evolution. In addition to flux injections, the relative motion of active regions would give rise to the catastrophic behavior and lead to magnetic eruptions as well. We propose that a gradual process could lead to a sudden release of magnetosphere energy on a very short dynamical timescale, without being initiated by a sudden fracture in the crust of the magnetar. Some implications of our model are also discussed.

The Rossby wave instability (RWI) is a promising mechanism for producing large-scale vortices in protoplanetary discs. The instability operates around a density bump in the disc, and the resulting vortices may facilitate planetesimal formation and angular momentum transfer in the disc dead zone. Most previous works on the RWI deal with two-dimensional (height-integrated) discs. However, vortices may have different dynamical behaviours in 3D than in 2D. Recent numerical simulations of the RWI in 3D global discs by Meheut et al. have revealed intriguing vertical structure of the vortices, including appreciable vertical velocities. In this paper we present a linear analysis of the RWI in 3D global models of isothermal discs. We calculate the growth rates of the Rossby modes (of various azimuthal wave numbers m = 2 - 6) trapped around the fiducial density bump and carry out 3D numerical simulations to compare with our linear results. We show that the 3D RWI growth rates are only slightly smaller than the 2D growth rates, and the velocity structures seen in the numerical simulations during the linear phase are in agreement with the velocity eigenfunctions obtained in our linear calculations. This numerical benchmark shows that numerical simulations can accurately describe the instability. The angular momentum transfer rate associated with Rossby vortices is also studied.

Motivated by coronal mass ejection studies, we construct general relativistic models of a magnetar magnetosphere endowed with strong magnetic fields. The equilibrium states of the stationary, axisymmetric magnetic fields in the magnetar magnetosphere are obtained as solutions of the Grad-Shafranov equation in a Schwarzschild spacetime. To understand the magnetic energy buildup in the magnetar magnetosphere, a generalized magnetic virial theorem in the Schwarzschild metric is newly derived. We carefully address the question whether the magnetar magnetospheric magnetic field can build up sufficient magnetic energy to account for the work required to open up the magnetic field during magnetar giant flares. We point out the importance of the Aly-Sturrock constraint, which has been widely studied in solar corona mass ejections, as a reference state in understanding magnetar energy storage processes. We examine how the magnetic field can possess enough energy to overcome the Aly-Sturrock energy constraint and open up. In particular, general relativistic (GR) effects on the Aly-Sturrock energy constraint in the Schwarzschild spacetime are carefully investigated. It is found that, for magnetar outbursts, the Aly-Sturrock constraint is more stringent, i.e., the Aly-Sturrock energy threshold is enhanced due to the GR effects. In addition, neutron stars with greater mass have a higher Aly-Sturrock energy threshold and are more difficult to erupt. This indicates that magnetars are probably not neutron stars with extreme mass. For a typical neutron star with mass of $1-2 M_{\odot}$, we further explore the effects of cross-field current effects, caused by the mass loading, on the possibility of stored magnetic field energy exceeding the Aly-Sturrock threshold.

The force-free (or low inertia) limit of magnetohydrodynamics (MHD) can be applied to many astrophysical objects, including black holes, neutron stars, and accretion disks, where the electromagnetic field is so strong that the inertia and pressure of the plasma can be ignored. This is difficult to achieve with the standard MHD numerical methods because they still have to deal with plasma inertial terms even when these terms are much smaller than the electromagnetic terms. Under the force free approximation, the plasma dynamics is entirely determined by the magnetic field. The plasma provides the currents and charge densities required by the dynamics of electromagnetic fields, but these currents carry no inertia. We present a high order Godunov scheme to study such force-free electrodynamics. We have implemented weighted essentially non-oscillatory (WENO) spatial interpolations in our scheme. An exact Riemann solver is implemented, which requires spectral decomposition into characteristic waves. We advance the magnetic field with the constrained transport (CT) scheme to preserve the divergence free condition to machine round-off error. We apply the third order total variation diminishing (TVD) Runge-Kutta scheme for the temporal integration. The mapping from face-centered variables to volume-centered variables is carefully considered. Extensive testing are performed to demonstrate the ability of our scheme to address force-free electrodynamics correctly. We finally apply the scheme to study relativistic magnetically dominated tearing instabilities and neutron star magnetospheres.

We carry out 2-D high resolution numerical simulations of type I planet migration with different disk viscosities. We find that the planet migration is strongly dependent on disk viscosities. Two kinds of density wave damping mechanisms are discussed. Accordingly, the angular momentum transport can be either viscosity dominated or shock dominated, depending on the disk viscosities. The long term migration behavior is different as well. Influences of the Rossby vortex instability on planet migration are also discussed. In addition, we investigate very weak shock generation in inviscid disks by small mass planets and compare the results with prior analytic results.

We investigate the global nonaxisymmetric Rossby vortex instability in a differentially rotating, compressible magnetized accretion disk with radial density structures. Equilibrium magnetic fields are assumed to have only the toroidal component. Using linear theory analysis, we show that the density structure can be unstable to nonaxisymmetric modes. We find that, for the magnetic field profiles we have studied, magnetic fields always provide a stabilizing effect to the unstable Rossby vortex instability modes. We discuss the physical mechanism of this stabilizing effect. The threshold and properties of the unstable modes are also discussed in detail. In addition, we present linear stability results for the global magnetorotational instability when the disk is compressible.

Sagittarius A* is a compact radio source at the Galactic center, powered by accretion of fully ionized plasmas into a supermassive black hole. However, the radio emission cannot be produced through the thermal synchrotron process by a gravitationally bounded flow. General relativistic magneto-hydrodynamical(GRMHD) simulations of black hole accretion show that there are strong unbounded outflows along the accretion. With the flow structure around the black hole given by GRMHD simulations, we investigate whether thermal synchrotron emission from these outflows may account for the observed radio emission. We find that simulations producing relatively high values of plasma beta cannot produce the radio flux level without exceeding the X-ray upper limit set by Chandra observations through the bremsstrahlung process. The predicted radio spectrum is also harder than the observed spectrum both for the one temperature thermal model and a simple nonthermal model with a single power-law electron distribution. The electron temperature needs to be lower than the gas temperature near the black hole to reproduce the observed radio spectrum. A more complete modeling of the radiation processes, including the general relativistic effects and transfer of polarized radiation, will give more quantitative constraints on physical processes in Sgr A* with the current multi-wavelength, multi-epoch, and polarimetric observations of this source.

We explore MHD solutions for envelope expansions with core collapse (EECC) with isothermal MHD shocks in a quasi-spherical symmetry and outline potential astrophysical applications of such magnetized shock flows. MHD shock solutions are classified into three classes according to the downstream characteristics near the core. Class I solutions are those characterized by free-fall collapses towards the core downstream of an MHD shock, while Class II solutions are those characterized by Larson-Penston (LP) type near the core downstream of an MHD shock. Class III solutions are novel, sharing both features of Class I and II solutions with the presence of a sufficiently strong magnetic field as a prerequisite. Various MHD processes may occur within the regime of these isothermal MHD shock similarity solutions, such as sub-magnetosonic oscillations, free-fall core collapses, radial contractions and expansions. We can also construct families of twin MHD shock solutions as well as an `isothermal MHD shock' separating two magnetofluid regions of two different yet constant temperatures. The versatile behaviours of such MHD shock solutions may be utilized to model a wide range of astrophysical problems, including star formation in magnetized molecular clouds, MHD link between the asymptotic giant branch phase to the proto-planetary nebula phase with a hot central magnetized white dwarf, relativistic MHD pulsar winds in supernova remnants, radio afterglows of soft gamma-ray repeaters and so forth.

We investigate self-similar magnetohydrodynamic (MHD) processes in an isothermal self-gravitating fluid with a quasi-spherical symmetry and extend the envelope expansion with core collapse (EECC) solutions of Lou & Shen by incorporating a random magnetic field. Stagnation surfaces of EECC solutions that seperate core collapse and envelope expansion propagate at constant speeds either sub-magnetosonically or super-magnetosonically. Crossing the magnetosonic line twice analytically, there exists an infinite number of discrete magnetized EECC and ECCC solutions. In addition to the EECC shock solution which could change the central accretion rate, the magnetic field can also affect the core accretion rate. As the magnetic parameter $\lambda$ increases, the core accretion rate appropriate for the MHD EWCS becomes larger. Under the frozen-in approximation, magnetic fields in the envelope expansion portion would scale as $B\propto r^{-1}$, while in the core collapse portion they would scale as $B\propto r^{-1/2}$. We discuss several astrophysical applications of EECC similarity solutions to the formation process of proto-planetary nebulae connecting the AGB phase and the planetary nebula, to supernova remnants, collapse of magnetized molecular cloud, to $\textrm{H}\ {\textrm{\small{II}}}$ clouds surrounding massive OB stars and to a certain evolution phase of galaxy clusters.