The discoveries of LIGO/Virgo black holes in recent years have revitalized the study of primordial black holes. In this work, we investigate a mechanism where primordial black holes are formed by vacuum bubbles that randomly nucleate during inflation through quantum tunneling. After inflation, these bubbles typically run into the ambient radiation fluid with a large Lorentz factor. In our previous work, we assumed the bubble fields are strongly coupled to the standard model particles so that the bubble wall is impermeable. Here we complete this picture by considering bubbles interacting with the fluid only through gravity. By studying the scenario in several limits, we found that black holes could form in either subcritical or supercritical regime. Depending on the model parameters, the resulting mass spectrum of the black holes could be wide or narrow, and may develop two peaks separated by a large mass range. With different spectra, these black holes may account for the LIGO/Virgo black holes, supermassive black holes, and may play an important role in dark matter.

Modern models of terrestrial planet formation require solids depletion interior to 0.5-0.7 au in the planetesimal disk to explain the small mass of Mercury. Earth and Venus analogues emerge after ~100 Myr collisional growth while Mercury form in the diffusive tails of the planetesimal disk. We carried out 250 N-body simulations of planetesimal disks with mass confined to 0.7-1.0 au to study the statistics of close encounters which were recently proposed as an explanation for the high iron mass fraction in Mercury by Deng (2020). We formed 39 Mercury analogues in total and all proto-Mercury analogues were scattered inward by proto-Venus. Proto-Mercury typically experiences 6 extreme close encounters (closest approach smaller than 6 Venus radii) with Proto-Venus after Proto-Venus acquires 0.7 Venus Mass. At such close separation, the tidal interaction can already affect the orbital motion significantly such that the N-body treatment itself is invalid. More and closer encounters are expected should tidal dissipation of orbital angular momentum accounted. Hybrid N-body hydrodynamic simulations, treating orbital and encounter dynamics self-consistently, are desirable to evaluate the degree of tidal mantle stripping of proto-Mercury.

The primary task of the 1.26-m telescope jointly operated by the National Astronomical Observatory and Guangzhou University is photometric observations of the g, r, and i bands. A data processing pipeline system was set up with mature software packages, such as IRAF, SExtractor, and SCAMP, to process approximately 5 GB of observational data automatically every day. However, the success ratio was significantly reduced when processing blurred images owing to telescope tracking error; this, in turn, significantly constrained the output of the telescope. We propose a robust automated photometric pipeline (RAPP) software that can correctly process blurred images. Two key techniques are presented in detail: blurred star enhancement and robust image matching. A series of tests proved that RAPP not only achieves a photometric success ratio and precision comparable to those of IRAF but also significantly reduces the data processing load and improves the efficiency.

Significant progress in the classification of Fermi unassociated sources , has led to an increasing number of blazars are being found. The optical spectrum is effectively used to classify the blazars into two groups such as BL Lacs and flat spectrum radio quasars (FSRQs). However, the accurate classification of the blazars without optical spectrum information, i.e., blazars of uncertain type (BCUs), remains a significant challenge. In this paper, we present a principal component analysis (PCA) and machine learning hybrid blazars classification method. The method, based on the data from Fermi LAT 3FGL Catalog, first used the PCA to extract the primary features of the BCUs and then used a machine learning algorithm to further classify the BCUs. Experimental results indicate that the that the use of PCA algorithms significantly improved the classification. More importantly, comparison with the Fermi LAT 4FGL Catalog, which contains the spectral classification of those BCUs in the Fermi-LAT 3FGL Catalog, reveals that the proposed classification method in the study exhibits higher accuracy than currently established methods; specifically, 151 out of 171 BL Lacs and 19 out of 24 FSRQs are correctly classified.

Quasi-biennial oscillations (QBOs) are considered as a fundamental mode of solar magnetic activity at low latitudes ($\leq50^\circ$). However, the evolutionary aspect and the hemispheric distribution of solar QBOs at high latitudes ($\geq60^\circ$) are rarely studied. Here, a relatively novel time-frequency analysis technique, named the synchrosqueezed wavelet transform, is applied to extract the main components of the polar faculae in the northern and southern hemispheres for the time interval from August 1951 to December 1998. It is found as the following: (1) Apart from the 22-year Hale cycle, the 17-year extended activity cycle, and the 11-year Schwabe cycle, the QBOs have been estimated as a prominent timescale of solar magnetic activity at high latitudes; (2) the QBOs of the polar faculae are coherent in the two hemispheres, but the temporal (phase) and the spatial (amplitude) variations of solar QBOs occur unevenly on both hemispheres; and (3) for the 11-year period mode, the northern hemisphere begins three months earlier than that in the southern one. Moreover, the spatial and temporal distributions of the hemispheric QBOs differ from those of the 11-year Schwabe cycle mode in the two hemispheres. Our findings could be helpful to improve our knowledge on the physical origin of the spatial distribution of solar QBOs at high latitudes, and could also provide more constraints on solar dynamo models introduced to characterize the different components of the solar magnetic activity cycle.

Primordial bubbles that possibly nucleate through quantum tunneling during inflation in a multidimensional potential might have left some relic detectable at the present time. These bubbles turn into black holes during the radiation era, which may account for the LIGO black holes, supermassive black holes, and may play an important role in dark matter. Typically, these black holes are surrounded by an energy deficit in the form of a spherical sound wave packet propagating outwards. In this work we study how this perturbation of the cosmic plasma dissipates before the time of recombination, leading to spectral distortions in CMB. We find that there may exist some rare regions on the last scattering surface containing huge black holes, which have produced potentially detectable point-like signals of $\mu$-type distortions.

In the recent years, the field of fast radio bursts (FRBs) is thriving and growing rapidly. It is of interest to study cosmology by using FRBs with known redshifts. In the present work, we try to test the possible cosmic anisotropy with the simulated FRBs. In particular, we only consider the possible dipole in FRBs, rather than the cosmic anisotropy in general, while the analysis is only concerned with finding the rough number of necessary data points to distinguish a dipole from a monopole structure through simulations. Noting that there is no a large sample of actual data of FRBs with known redshifts by now, simulations are necessary to this end. We find that at least 2800, 190, 100 FRBs are competent to find the cosmic dipole with amplitude 0.01, 0.03, 0.05, respectively. Unfortunately, even 10000 FRBs are not competent to find the tiny cosmic dipole with amplitude of ${\cal O}(10^{-3})$. On the other hand, at least 20 FRBs with known redshifts are competent to find the cosmic dipole with amplitude 0.1. We expect that such a big cosmic dipole could be ruled out by using only a few tens of FRBs with known redshifts in the near future.

In the early stages of a protoplanetary disk, when its mass is a significant fraction of its star's, turbulence generated by gravitational instability (GI) should feature significantly in the disk's evolution. At the same time, the disk may be sufficiently ionised for magnetic fields to play some role in the dynamics. Though usually neglected, the impact of magnetism on the GI may be critical, with consequences for several processes: the efficiency of accretion, spiral structure formation, fragmentation, and the dynamics of solids. In this paper, we report on global three-dimensional magnetohydrodynamical simulations of a self-gravitating protoplanetary disk using the meshless finite mass (MFM) Lagrangian technique. We confirm that GI spiral waves trigger a dynamo that amplifies an initial magnetic field to nearly thermal amplitudes (plasma beta < 10), an order of magnitude greater than that generated by the magneto-rotational instability alone. We also determine the dynamo's nonlinear back reaction on the gravitoturbulent flow: the saturated state is substantially hotter, with an associated larger Toomre parameter and weaker, more 'flocculent' spirals. But perhaps of greater import is the dynamo's boosting of accretion via a significant Maxwell stress; mass accretion is enhanced by factors of several relative to either pure GI or pure MRI. Our simulations use ideal MHD, an admittedly poor approximation in protoplanetary disks, and thus future studies should explore the full gamut of non-ideal MHD. In preparation for that, we exhibit a small number of Ohmic runs that reveal that the dynamo, if anything, is stronger in a non-ideal environment. This work confirms that magnetic fields are a potentially critical ingredient in gravitoturbulent young disks, possibly controlling their evolution, especially via their enhancement of (potentially episodic) accretion.

Mercury has an unusually large metal core comprising ~70% of its mass comparing to all other terrestrial planets in the solar system. Giant impacts can remove a significant fraction of the silicate mantle of a chondritic proto-Mercury and form the iron rich present-day Mercury. However, such high-temperature giant impacts seem at odds with the retainment of moderately volatile elements on present-day Mercury (Peplowski et al. 2011). We simulated a series of hyperbolic encounters between proto-Mercury and proto-Venus, which may occur in the chaotic early solar system. Tidal disruption of proto-Mercury always removes part of its silicate mantle while its iron core remains intact. We find, in favourable cases, four close encounters with fast spinning projectiles (resulting from previous encounters) can lead to present-day Mercury iron fraction. More encounters are needed when the spin and orbital angular momentum are not always aligned. These hyperbolic encounters have various outcomes, such as orbital decay, binary planets and change of spin rates. These results suggest the importance of proper treatment of close encounters in N-body simulations of planetary accretion.

Gravitational wave (GW) has become one of the most active fields in physics and astronomy since the first direct detection of GW event in 2015. As is well known, multiple images of GW events are possible through the gravitational lenses. Here, we propose a novel mirror imaging mechanism for GW events different from the gravitational lens. In this work, we predict that there are two types of GW mirror imaging phenomena caused by the neutron star located in Milky Way or the same host galaxy of GW source, which might be detected within a life period of man (namely the time delay $\Delta t$ can be a few years to a few tens of years). It is expected to witness this predicted GW mirror imaging phenomenon in the near future. In the long term, the observations of this novel GW mirror imaging phenomenon might help us to find numerous neutron stars unseen by other means, and learn more about the complicated internal structures of neutron stars, as well as their equations of state.

Temporal variation of the solar coronal rotation appears to be very complex and its relevances to the eleven-year solar activity cycle are still unclear. Using the modified coronal index for the time interval from 1939 January 1 to 2019 May 31, the systematic regularities of the solar coronal rotation are investigated. Our main findings are as follows: (1) from a global point of view, the synodic coronal rotation period with a value of 27.5 days is the only significant period at the periodic scales shorter than 64 days; (2) the coronal rotation period exhibit an obviously decreasing trend during the considered time interval, implying the solar corona accelerates its global rotation rate in the long run; (3) there exist significant periods of 3.25, 6.13, 9.53, and 11.13 years in the period length of the coronal rotation, providing an evidence that the coronal rotation should be connected with the quasi-biennial oscillation, the eleven-year solar cycle, and the 22-year Hale cycle (or the magnetic activity reversal); and (4) the phase relationship between the coronal rotation period and the solar magnetic activity is not only time-dependent but also frequency-dependent. For a small range around the 11- year cycle band, there is a systematic trend in the phase, and the small mismatch in this band brings out the phase to drift. The possible mechanism for the above analysis results is discussed.

Quasi-biennial oscillation (QBO) of solar magnetic activities is intrinsic to dynamo mechanism, but still far from fully understood. In this work, the phase and amplitude asymmetry of solar QBO of Halpha flare activity in the northern and southern hemispheres is studied by the ensemble empirical mode decomposition, the cross-correlation analysis, and the wavelet transform technique. The following results are found: (1) solar QBO of Halpha flare index in the two hemispheres has a complicated phase relationship, but does not show any systematic regularity; (2) the solar cycle mode of solar Halpha flare index in the northern hemisphere generally leads that in the southern one by 9 months for the time interval from 1966 to 2014. The possible origin of these results is discussed.

Studies on the periodic variation and the phase relationship between different solar activity indicators are useful for understanding the long-term evolution of solar activity cycle. Here we report the statistical analysis of grouped solar flare (GSF) and sunspot number (SN) during the time interval from January 1965 to March 2009. We find that, 1) the significant periodicities of both GSF and SN are related to the differential rotation periodicity, the quasi-biennial oscillation (QBO), and the eleven-year Schwabe cycle (ESC), but the specific values are not absolutely identical; 2) the ESC signal of GSF lags behind that of SN with an average of 7.8 months during the considered time interval, implying that the systematic phase delays between GSF and SN originate from the inter-solar-cycle signal. Our results may provide evidence about the storage of magnetic energy in the corona.

The giant impact hypothesis for Moon formation successfully explains the dynamic properties of the Earth-Moon system but remains challenged by the similarity of isotopic fingerprints of the terrestrial and lunar mantles. Moreover, recent geochemical evidence suggests that the Earth's mantle preserves ancient (or "primordial") heterogeneity that predates the Moon-forming giant impact. Using a new hydrodynamical method, we here show that Moon-forming giant impacts lead to a stratified starting condition for the evolution of the terrestrial mantle. The upper layer of the Earth is compositionally similar to the disk, out of which the Moon evolves, whereas the lower layer preserves proto-Earth characteristics. As long as this predicted compositional stratification can at least partially be preserved over the subsequent billions of years of Earth mantle convection, the compositional similarity between the Moon and the accessible Earth's mantle is a natural outcome of realistic and high-probability Moon-forming impact scenarios. The preservation of primordial heterogeneity in the modern Earth not only reconciles geochemical constraints but is also consistent with recent geophysical observations. Furthermore, for significant preservation of a proto-Earth reservoir, the bulk composition of the Earth-Moon system may be systematically shifted towards chondritic values.

The magneto-rotational instability (MRI) is one of the most important processes in sufficiently ionized astrophysical disks. Grid-based simulations, especially those using the local shearing box approximation, provide a powerful tool to study the ensuing nonlinear turbulence. On the other hand, while meshless methods have been widely used in both cosmology, galactic dynamics, and planet formation they have not been fully deployed on the MRI problem. We present local unstratified and vertically stratified MRI simulations with two meshless MHD schemes: a recent implementation of SPH MHD (Price2012), and a MFM MHD scheme with a constrained gradient divergence cleaning scheme, as implemented in the GIZMO code \citep{Hopkins2017}. Concerning variants of the SPH hydro force formulation we consider both the "vanilla" SPH and the PSPH variant included in GIZMO. We find, as expected, that the numerical noise inherent in these schemes affects turbulence significantly. A high order kernel, free of the pairing instability, is necessary. Both schemes can adequately simulate MRI turbulence in unstratified shearing boxes with net vertical flux. The turbulence, however, dies out in zero-net-flux unstratified boxes, probably due to excessive and numerical dissipation. In zero-net-flux vertically stratified simulations, MFM can reproduce the MRI dynamo and its characteristic butterfly diagram for several tens of orbits before ultimately decaying. In contrast, extremely strong toroidal fields, as opposed to sustained turbulence, develop in equivalent simulations using SPH MHD. This unphysical state in SPH MHD is likely caused by a combination of excessive artificial viscosity, numerical resistivity, and the relatively large residual errors in the divergence of the magnetic field remaining even after cleaning procedures are applied.

Gas clumps formed within massive gravitationally unstable circumstellar discs are potential seeds of gas giant planets, brown dwarfs and companion stars. Simulations show that competition between three processes -- migration, gas accretion and tidal disruption -- establishes what grows from a given seed. Here we investigate the robustness of numerical modelling of clump migration and accretion with the codes PHANTOM, GADGET, SPHINX, SEREN, GIZMO-MFM, SPHNG and FARGO. The test problem comprises a clump embedded in a massive disc at an initial separation of 120 AU. There is a general qualitative agreement between the codes, but the quantitative agreement in the planet migration rate ranges from $\sim 10$% to $\sim 50$%, depending on the numerical setup. We find that the artificial viscosity treatment and the sink particle prescription may account for much of the differences between the codes. In order to understand the wider implications of our work, we also attempt to reproduce the planet evolution tracks from our hydrodynamical simulations with prescriptions from three previous population synthesis studies. We find that the disagreement amongst the population synthesis models is far greater than that between our hydrodynamical simulations. The results of our code comparison project are therefore encouraging in that uncertainties in the given problem are probably dominated by the physics not yet included in the codes rather than by how hydrodynamics is modelled in them.

The Earths core formation process has decisive effect in the chemical differentiation between the Earths core and its mantle. Here, we propose a new core formation model which is caused by a special giant impact. This model suggests that the impactors core can be kept intact by its own sticky mantle under appropriate impacting conditions and let it merge into the targets core without contact with the targets mantle. We call this special giant impact that caused the new core formation mode as glue ball impact model (GBI). By simulating hundreds of giant impacts with the sizes from planetesimals to planets, the conditions that can lead to GBI have been found out. If with small impact angle (i.e., less than 20 degree), small impact velocity and small impactors mass but larger than 0.07 Mearth, there is a good chance to produce a GBI at the final stage of the Earths accretion. We find that it will be much easier to have GBIs at the late stage of the Earths accretion rather than at the early stage of it. The GBI model will pose a great challenge to many problems between the equilibrium of Earths core and mantle. It provides an additional source for the excess of highly siderophile elements in the Earths mantle and also brings excessive lithophile elements to the Earths core. The GBI model may shed light on the study of Moon-formation and chemical differentiations of the pro-Earth.

Giant impacts (GIs) are common in the late stage of planet formation. The Smoothed Particle Hydrodynamics (SPH) method is widely used for simulating the outcome of such violent collisions, one prominent example being the formation of the Moon. However, a decade of numerical studies in various areas of computational astrophysics has shown that the standard formulation of SPH suffers from several shortcomings such as artificial surface tension and its tendency to promptly damp turbulent motions on scales much larger than the physical dissipation scale, both resulting in the suppression of mixing. In order to quantify how severe these limitations are when modeling GIs we carried out a comparison of simulations with identical initial conditions performed with the standard SPH as well as with the novel Lagrangian Meshless Finite Mass (MFM) method in the GIZMO code. We confirm the lack of mixing between the impactor and target when SPH is employed, while MFM is capable of driving vigorous sub-sonic turbulence and leads to significant mixing between the two bodies. Modern SPH variants with artificial conductivity, a different formulation of the hydro force or reduced artificial viscosity, do not improve mixing as significantly. Angular momentum is conserved similarly well in both methods, but MFM does not suffer from spurious transport induced by artificial viscosity, resulting in a slightly higher angular momentum of the proto-lunar disk. Furthermore, SPH initial conditions exhibit an unphysical density discontinuity at the core-mantle boundary which is easily removed in MFM.

Using the Hilbert-Huang transform technique, we investigate the midrange periodicities in solar radio flux at 2800 MHz (F10.7) and sunspot areas (SAs) from February 1, 1947 to September 30, 2016. The following prominent results are found: (1) The quasi-periodic oscillations of both data sets are not identical, such as the rotational cycle, the midrange periodicities, and the Schwabe cycle. In particular, the midrange periodicities ranging from 37.9 days to 297.3 days are related to the magnetic Rossby-type waves; 2) The 1.3-year and 1.7-year fluctuations in solar activity indicators are surface manifestations (from photosphere to corona) of magnetic flux changes generated deep inside the Sun; 3) At the timescale of the Schwabe cycle, \textbf{the complicated phase relationships} in the three intervals (1947-1958, 1959-1988, and 1989-2016) agree with the produced periodicities of the magnetic Rossby-type waves. \textbf{The findings indicate that the magnetic Rossby-type waves are the possible physical mechanism behind the midrange periodicities of solar activity indicators. Moreover, the significant change in the relationship between photospheric and coronal activity took place after the maximum of solar cycle 22 could be interpreted by the magnetic Rossby-type waves.

The cosmological principle assumes that the universe is homogeneous and isotropic on cosmic scales. There exist many works testing the cosmic homogeneity and/or the cosmic isotropy of the universe in the literature. In fact, some observational hints of the cosmic anisotropy have been claimed. However, we note that the paucity of the data considered in the literature might be responsible for the `found' cosmic anisotropy. So, it might disappear in a large enough sample. Very recently, the Pantheon sample consisting of 1048 type Ia supernovae (SNIa) has been released, which is the largest spectroscopically confirmed SNIa sample to date. In the present work, we test the cosmic anisotropy in the Pantheon SNIa sample by using three methods, and hence the results from different methods can be cross-checked. All the results obtained by using the hemisphere comparison (HC) method, the dipole fitting (DF) method and HEALPix suggest that no evidence for the cosmic anisotropy is found in the Pantheon SNIa sample.

Recently there has been much interest in light dark matter, especially ultra-light axions, as they may provide a solution to the core-cusp problem at the center of galaxies. Since very light bosons can have a de Broglie wavelength that is of astrophysical size, they can smooth out the centers of galaxies to produce a core, as opposed to vanilla dark matter models, and so it has been suggested that this solves the core-cusp problem. In this work, we critically examine this claim. While an ultra-light particle will indeed lead to a core, we examine whether the relationship between the density of the core and its radius matches the data over a range of galaxies. We first review data that shows the core density of a galaxy $\rho_c$ varies as a function of the core radius $R_c$ as $\rho_c\propto1/R_c^\beta$ with $\beta\approx1$. We then compare this to theoretical models. We examine a large class of light scalar dark matter models, governed by some potential $V$. For simplicity, we take the scalar to be complex with a global $U(1)$ symmetry in order to readily organize solutions by a conserved particle number. However, we expect our central conclusions to persist even for a real scalar, and furthermore, a complex scalar matches the behavior of a real scalar in the non-relativistic limit, which is the standard regime of interest. For any potential $V$, we find the relationship between $\rho_c$ and $R_c$ for ground state solutions is always in one of the following regimes: (i) $\beta\gg1$, or (ii) $\beta\ll1$, or (iii) unstable, and so it never matches the data. We also find similar conclusions for virialized dark matter, more general scalar field theories, degenerate fermion dark matter, superfluid dark matter, and general polytropes. We conclude that the solution to the core-cusp problem is more likely due to either complicated baryonic effects or some other type of dark matter interactions.

The cosmological principle is one of the cornerstones in modern cosmology. It assumes that the universe is homogeneous and isotropic on cosmic scales. Both the homogeneity and the isotropy of the universe should be tested carefully. In the present work, we are interested in probing the possible preferred direction in the distribution of type Ia supernovae (SNIa). To our best knowledge, two main methods have been used in almost all of the relevant works in the literature, namely the hemisphere comparison (HC) method and the dipole fitting (DF) method. However, the results from these two methods are not always approximately coincident with each other. In this work, we test the cosmic anisotropy by using these two methods with the Joint Light-Curve Analysis (JLA) and simulated SNIa datasets. In many cases, both methods work well, and their results are consistent with each other. However, in the cases with two (or even more) preferred directions, the DF method fails while the HC method still works well. This might shed new light on our understanding of these two methods.

Vacuum bubbles may nucleate and expand during the cosmic inflation. When inflation ends, the bubbles run into the ambient plasma, producing strong shocks followed by underdensity waves, which propagate outwards. The bubbles themselves eventually form black holes with a wide distribution of masses. It has been recently suggested that such black holes may account for LIGO observations and may provide seeds for supermassive black holes observed at galactic centers. They may also provide a significant part or even the whole of the dark matter. We estimate the spectral $\mu$-distortion of the CMB induced by expanding shocks and underdensities. The predicted distortions averaged over the sky are well below the current bounds, but localized peaks due to the largest black holes impose constraints on the model parameters.

This paper presents a fast algorithm for obtaining high-accuracy subpixel translation of low PSNR images. Instead of locating the maximum point on the upsampled images or fitting the peak of correlation surface, the proposed algorithm is based on the measurement of centroid on the cross correlation surface by Modified Moment method. Synthetic images, real solar images and standard testing images with white Gaussian noise added were tested, and the results show that the accuracies of our algorithm are comparable with other subpixel registration techniques and the processing speed is higher. The drawback is also discussed at the end of this paper.

We use high spatial and temporal resolution observations, simultaneously obtained with the New Vacuum Solar Telescope and Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory, to investigate the high-frequency oscillations above a sunspot umbra. A novel time--frequency analysis method, namely the synchrosqueezing transform (SST), is employed to represent their power spectra and to reconstruct the high-frequency signals at different solar atmospheric layers. A validation study with synthetic signals demonstrates that SST is capable to resolving weak signals even when their strength is comparable with the high-frequency noise. The power spectra, obtained from both SST and the Fourier transform, of the entire umbral region indicate that there are significant enhancements between 10 and 14 mHz (labeled as 12 mHz) at different atmospheric layers. Analyzing the spectrum of a photospheric region far away from the umbra demonstrates that this 12~mHz component exists only inside the umbra. The animation based on the reconstructed 12 mHz component in AIA 171 \AA\ illustrates that an intermittently propagating wave first emerges near the footpoints of coronal fan structures, and then propagates outward along the structures. A time--distance diagram, coupled with a subsonic wave speed ($\sim$ 49 km s$^{-1}$), highlights the fact that these coronal perturbations are best described as upwardly propagating magnetoacoustic slow waves. Thus, we first reveal the high-frequency oscillations with a period around one minute in imaging observations at different height above an umbra, and these oscillations seem to be related to the umbral perturbations in the photosphere.

A novel time--frequency technique, called the synchrosqueezing transform (SST), is used to investigate the midterm periodic variations of magnetic fields on the solar surface. The Magnetic Plage Strength Index (MPSI) and the Mount Wilson Sunspot Index (MWSI), measured daily by the Mount Wilson Observatory between 1970 January 19 and 2012 January 22, are selected. The short-, mid, and longer-term periodicities are represented and decomposed by the SST with hardly any mode mixing. This demonstrates that the SST is a useful time--frequency analysis technique to characterize the periodic modes of helioseismic data. Apart from the fundamental modes of the annual periodicity, $\sim$27 day rotational cycle and $\sim$11 year solar cycle, the SST reveals several midterm periodicities in the two magnetic activity indices, specifically, $\sim$157 day (i.e., Rieger-type periodicity), and $\sim$1.3 and 1.7 years. The periodic modes, with 116.4 and 276.2 day periodicity in the MPSI, with 108.5 and 251.6 day periodicity in the MWSI, and the 157.7 day periodicity in the two indices, are in better accord with those significant periodicities derived from the Rossby waves theoretical model. This study suggests that the modes are caused by the Rossby waves. For the 1.30 and 1.71 year periodicity of the MPSI, and the 1.33 and 1.67 year periodicity of the MWSI, our analysis infers that they are related to those periodicity with the same timescale in the interior of the Sun and in the high atmospheric layers.

Vacuum bubbles may nucleate during the inflationary epoch and expand, reaching relativistic speeds. After inflation ends, the bubbles are quickly slowed down, transferring their momentum to a shock wave that propagates outwards in the radiation background. The ultimate fate of the bubble depends on its size. Bubbles smaller than certain critical size collapse to ordinary black holes, while in the supercritical case the bubble interior inflates, forming a baby universe, which is connected to the exterior region by a wormhole. The wormhole then closes up, turning into two black holes at its two mouths. We use numerical simulations to find the masses of black holes formed in this scenario, both in subcritical and supercritical regime. The resulting mass spectrum is extremely broad, ranging over many orders of magnitude. For some parameter values, these black holes can serve as seeds for supermassive black holes and may account for LIGO observations.

As a dedicated solar radio interferometer, the MingantU SpEctral RadioHeliograph (MUSER) generates massive observational data in the frequency range of 400 MHz -- 15 GHz. High-performance imaging forms a significantly important aspect of MUSER's massive data processing requirements. In this study, we implement a practical high-performance imaging pipeline for MUSER data processing. At first, the specifications of the MUSER are introduced and its imaging requirements are analyzed. Referring to the most commonly used radio astronomy software such as CASA and MIRIAD, we then implement a high-performance imaging pipeline based on the Graphics Processing Unit (GPU) technology with respect to the current operational status of the MUSER. A series of critical algorithms and their pseudo codes, i.e., detection of the solar disk and sky brightness, automatic centering of the solar disk and estimation of the number of iterations for clean algorithms, are proposed in detail. The preliminary experimental results indicate that the proposed imaging approach significantly increases the processing performance of MUSER and generates images with high-quality, which can meet the requirements of the MUSER data processing.

Since Lorentz invariance plays an important role in modern physics, it is of interest to test the possible Lorentz invariance violation (LIV). The time-lag (the arrival time delay between light curves in different energy bands) of Gamma-ray bursts (GRBs) has been extensively used to this end. However, to our best knowledge, one or more particular cosmological models were assumed {\it a priori} in (almost) all of the relevant works in the literature. So, this makes the results on LIV in those works model-dependent and hence not so robust in fact. In the present work, we try to avoid this problem by using a model-independent approach. We calculate the time delay induced by LIV with the cosmic expansion history given in terms of cosmography, without assuming any particular cosmological model. Then, we constrain the possible LIV with the observational data, and find weak hints for LIV.

We carry out simulations of gravitationally unstable disks using smoothed particle hydrodynamics(SPH) and the novel Lagrangian meshless finite mass (MFM) scheme in the GIZMO code (Hopkins 2015). Our aim is to understand the cause of the non-convergence of the cooling boundary for fragmentation reported in the literature. We run SPH simulations with two different artificial viscosity implementations, and compare them with MFM, which does not employ any artificial viscosity. With MFM we demonstrate convergence of the critical cooling time scale for fragmentation at \beta_{crit} =3.. Non-convergence persists in SPH codes, although it is significantly mitigated with schemes having reduced artificial viscosity such as inviscid SPH (ISPH) (Cullen & Dehnen 2010). We show how the non-convergence problem is caused by artificial fragmentation triggered by excessive dissipation of angular momentum in domains with large velocity derivatives. With increased resolution such domains become more prominent. Vorticity lags behind density due to numerical viscous dissipation in these regions, promoting collapse with longer cooling times. Such effect is shown to be dominant over the competing tendency of artificial viscosity to diminish with increasing resolution. When the initial conditions are first relaxed for several orbits, the flow is more regular, with lower shear and vorticity in non-axisymmetric regions, aiding convergence. Yet MFM is the only method that converges exactly. Our findings are of general interest as numerical dissipation via artificial viscosity or advection errors can also occur in grid-based codes. Indeed for the FARGO code values of \beta_{crit} significantly higher than our converged estimate have been reported in the literature. Finally, we discuss implications for giant planet formation via disk instability.

As a dedicated synthetic aperture radio interferometer, the MingantU SpEctral Radioheliograph (MUSER), initially known as the Chinese Spectral RadioHeliograph (CSRH), has entered the stage of routine observation. More than 23 million data records per day need to be effectively managed to provide high performance data query and retrieval for scientific data reduction. In light of these massive amounts of data generated by the MUSER, in this paper, a novel data management technique called the negative database (ND) is proposed and used to implement a data management system for the MUSER. Based on the key-value database, the ND technique makes complete utilization of the complement set of observational data to derive the requisite information. Experimental results showed that the proposed ND can significantly reduce storage volume in comparison with a relational database management system (RDBMS). Even when considering the time needed to derive records that were absent, its overall performance, including querying and deriving the data of the ND, is comparable with that of an RDBMS. The ND technique effectively solves the problem of massive data storage for the MUSER, and is a valuable reference for the massive data management required in next-generation telescopes.

All next generation ground-based and space-based solar telescopes require a good quality assessment metric in order to evaluate their imaging performance. In this paper, a new image quality metric, the median filter gradient similarity (MFGS) is proposed for photospheric images. MFGS is a no-reference/blind objective image quality metric (IQM) by a measurement result between 0 and 1 and has been performed on short-exposure photospheric images captured by the New Vacuum Solar Telescope (NVST) of the Fuxian Solar Observatory and by the Solar Optical Telescope (SOT) onboard the Hinode satellite, respectively. The results show that: (1)the measured value of MFGS changes monotonically from 1 to 0 with degradation of image quality; (2)there exists a linear correlation between the measured values of MFGS and root-mean-square-contrast (RMS-contrast) of granulation; (3)MFGS is less affected by the image contents than the granular RMS-contrast. Overall, MFGS is a good alternative for the quality assessment of photospheric images.

The volume of data generated by modern astronomical telescopes is extremely large and rapidly growing. However, current high-performance data processing architectures/frameworks are not well suited for astronomers because of their limitations and programming difficulties. In this paper, we therefore present OpenCluster, an open-source distributed computing framework to support rapidly developing high-performance processing pipelines of astronomical big data. We first detail the OpenCluster design principles and implementations and present the APIs facilitated by the framework. We then demonstrate a case in which OpenCluster is used to resolve complex data processing problems for developing a pipeline for the Mingantu Ultrawide Spectral Radioheliograph. Finally, we present our OpenCluster performance evaluation. Overall, OpenCluster provides not only high fault tolerance and simple programming interfaces, but also a flexible means of scaling up the number of interacting entities. OpenCluster thereby provides an easily integrated distributed computing framework for quickly developing a high-performance data processing system of astronomical telescopes and for significantly reducing software development expenses.

The New Vacuum Solar Telescope (NVST) is a 1-m solar telescope that aims to observe the fine structures in both the photosphere and the chromosphere of the Sun. The observational data acquired simultaneously from one channel for the chromosphere and two channels for the photosphere bring great challenges to the data storage of NVST. The multi-channel instruments of NVST, including scientific cameras and multi-band spectrometers, generate at least 3 terabytes data per day and require high access performance while storing massive short-exposure images. It is worth studying and implementing a storage system for NVST which would balance the data availability, access performance and the cost of development. In this paper, we build a distributed data storage system (DDSS) for NVST and then deeply evaluate the availability of real-time data storage on a distributed computing environment. The experimental results show that two factors, i.e., the number of concurrent read/write and the file size, are critically important for improving the performance of data access on a distributed environment. Referring to these two factors, three strategies for storing FITS files are presented and implemented to ensure the access performance of the DDSS under conditions of multi-host write and read simultaneously. The real applications of the DDSS proves that the system is capable of meeting the requirements of NVST real-time high performance observational data storage. Our study on the DDSS is the first attempt for modern astronomical telescope systems to store real-time observational data on a low-cost distributed system. The research results and corresponding techniques of the DDSS provide a new option for designing real-time massive astronomical data storage system and will be a reference for future astronomical data storage.

The New Vacuum Solar Telescope (NVST) is a 1-meter vacuum solar telescope that aims to observe the fine structures of active regions on the Sun. The main tasks of the NVST are high resolution imaging and spectral observations, including the measurements of the solar magnetic field. The NVST has been collecting more than 20 million FITS files since it began routine observations in 2012 and produces a maximum observational records of 120 thousand files in a day. Given the large amount of files, the effective archiving and retrieval of files becomes a critical and urgent problem. In this study, we implement a new data archiving system for the NVST based on the Fastbit Not Only Structured Query Language (NoSQL) database. Comparing to the relational database (i.e., MySQL; My Structured Query Language), the Fastbit database manifests distinctive advantages on indexing and querying performance. In a large scale database of 40 million records, the multi-field combined query response time of Fastbit database is about 15 times faster and fully meets the requirements of the NVST. Our study brings a new idea for massive astronomical data archiving and would contribute to the design of data management systems for other astronomical telescopes.

The Chinese Spectral RadioHeliograph (CSRH) is a synthetic aperture radio interferometer built in Inner Mongolia, China. As a solar-dedicated interferometric array, CSRH is capable of producing high quality radio images at frequency range from 400 MHz to 15 GHz with high temporal, spatial, and spectral resolution.To implement high cadence imaging at wide-band and obtain more than 2 order higher multiple frequencies, the implementation of the data processing system for CSRH is a great challenge. It is urgent to build a pipeline for processing massive data of CSRH generated every day. In this paper, we develop a high performance distributed data processing pipeline (DDPP) built on the OpenCluster infrastructure for processing CSRH observational data including data storage, archiving, preprocessing, image reconstruction, deconvolution, and real-time monitoring. We comprehensively elaborate the system architecture of the pipeline and the implementation of each subsystem. The DDPP is automatic, robust, scalable and manageable. The processing performance under multi computers parallel and GPU hybrid system meets the requirements of CSRH data processing. The study presents an valuable reference for other radio telescopes especially aperture synthesis telescopes, and also gives an valuable contribution to the current and/or future data intensive astronomical observations.

G-band bright points (GBPs) are regarded as good manifestations of magnetic flux concentrations. We aim to investigate the relationship between the dynamic properties of GBPs and their longitudinal magnetic field strengths. High spatial and temporal resolution observations were recorded simultaneously with G-band filtergrams and Narrow-band Filter Imager (NFI) Stokes I and V images with Hinode /Solar Optical Telescope. The GBPs are identified and tracked in the G-band images automatically, and the corresponding longitudinal magnetic field strength of each GBP is extracted from the calibrated NFI magnetograms by a point-to-point method. After categorizing the GBPs into five groups by their longitudinal magnetic field strengths, we analyze the dynamics of GBPs of each group. The results suggest that with increasing longitudinal magnetic field strengths of GBPs correspond to a decrease in their horizontal velocities and motion ranges as well as by showing more complicated motion paths. This suggests that magnetic elements showing weaker magnetic field strengths prefer to move faster and farther along straighter paths, while stronger ones move more slowly in more erratic paths within a smaller region. The dynamic behaviors of GBPs with different longitudinal magnetic field strengths can be explained by that the stronger flux concentrations withstand the convective flows much better than weaker ones.

Umbral dots (UDs) are small isolated brightenings observed in sunspot umbrae. They are convective phenomena existing inside umbrae. UDs are usually divided into central UDs (CUDs) and peripheral UDs (PUDs) according to their positions inside an umbra. Our purpose is to investigate UD properties and analyze their relationships, and further to find whether or not the properties depend on umbral magnetic field strengths. Thus, we selected high-resolution TiO images of four active regions (ARs) taken under the best seeing conditions with the \textit{New Vacuum Solar Telescope} in the Fuxian Solar Observatory of the Yunnan Astronomical Observatory, China. The four ARs (NOAA 11598, 11801, 12158, and 12178) include six sunspots. A total of 1220 CUDs and 603 PUDs were identified. Meanwhile, the radial component of the vector magnetic field of the sunspots taken with the \textit{Helioseismic and Magnetic Imager} on-board the \textit{Solar Dynamics Observatory} was used to analyze relationships between UD properties and umbral magnetic field strengths. We find that diameters and lifetimes of UDs exhibit an increasing trend with the brightness, but velocities do not. Moreover, diameters, intensities, lifetimes and velocities depend on the surrounding magnetic field. A CUD diameter was found larger, the CUD brighter, its lifetime longer, and its motion slower in a weak umbral magnetic field environment than in a strong one.

Six high-resolution TiO-band image sequences from the New Vacuum Solar Telescope (NVST) are used to investigate the properties of intergranular bright points (igBPs). We detect the igBPs using a Laplacian and morphological dilation algorithm (LMD) and track them using a three-dimensional segmentation algorithm automatically, and then investigate the morphologic, photometric and dynamic properties of igBPs, in terms of equivalent diameter, the intensity contrast, lifetime, horizontal velocity, diffusion index, motion range and motion type. The statistical results confirm the previous studies based on G-band or TiO-band igBPs from the other telescopes. It illustrates that the TiO data from the NVST have a stable and reliable quality, which are suitable for studying the igBPs. In addition, our method is feasible to detect and track the igBPs in the TiO data from the NVST. With the aid of the vector magnetograms obtained from the Solar Dynamics Observatory /Helioseismic and Magnetic Imager, the properties of igBPs are found to be influenced by their embedded magnetic environments strongly. The area coverage, the size and the intensity contrast values of igBPs are generally larger in the regions with higher magnetic flux. However, the dynamics of igBPs, including the horizontal velocity, the diffusion index, the ratio of motion range and the index of motion type are generally larger in the regions with lower magnetic flux. It suggests that the absence of strong magnetic fields in the medium makes it possible for the igBPs to look smaller and weaker, diffuse faster, move faster and further in a straighter path.

We study the motions of G band bright points (GBPs) in the quiet Sun to obtain the characteristics of different motion types. A high resolution image sequence taken with the Hinode/Solar Optical Telescope (SOT) is used, and GBPs are automatically tracked by segmenting 3D evolutional structures in a space time cube. After putting the GBPs that do not move during their lifetimes aside, the non stationary GBPs are categorized into three types based on an index of motion type. Most GBPs that move in straight or nearly straight lines are categorized into a straight motion type, a few moving in rotary paths into a rotary motion, and the others fall into a motion type we called erratic. The mean horizontal velocity is 2.18 km/s, 1.63 km/s and 1.33 km/s for straight, erratic and rotary motion type, respectively. We find that a GBP drifts at a higher and constant velocity during its whole life if it moves in a straight line. However, it has a lower and variational velocity if it moves in a rotary path. The diffusive process is ballistic, super and sub diffusion for straight, erratic and rotary motion type, respectively. The corresponding diffusion index and coefficients are 2.13 and 850 km2/s, 1.82 and 331 km2/s, 0.73 and 13 km2/s. In terms of direction of motion, it is homogeneous and isotropical, and usually persists between neighbouring frames, no matter what motion type a GBP belongs to.

The Sudbury Neutrino Observatory (SNO) has confirmed the standard solar model and neutrino oscillations through the observation of neutrinos from the solar core. In this paper we present a search for neutrinos associated with sources other than the solar core, such as gamma-ray bursters and solar flares. We present a new method for looking for temporal coincidences between neutrino events and astrophysical bursts of widely varying intensity. No correlations were found between neutrinos detected in SNO and such astrophysical sources.

We report results from a combined analysis of solar neutrino data from all phases of the Sudbury Neutrino Observatory. By exploiting particle identification information obtained from the proportional counters installed during the third phase, this analysis improved background rejection in that phase of the experiment. The combined analysis resulted in a total flux of active neutrino flavors from 8B decays in the Sun of (5.25 \pm 0.16(stat.)+0.11-0.13(syst.))\times10^6 cm^{-2}s^{-1}. A two-flavor neutrino oscillation analysis yielded \Deltam^2_{21} = (5.6^{+1.9}_{-1.4})\times10^{-5} eV^2 and tan^2{\theta}_{12}= 0.427^{+0.033}_{-0.029}. A three-flavor neutrino oscillation analysis combining this result with results of all other solar neutrino experiments and the KamLAND experiment yielded \Deltam^2_{21} = (7.41^{+0.21}_{-0.19})\times10^{-5} eV^2, tan^2{\theta}_{12} = 0.446^{+0.030}_{-0.029}, and sin^2{\theta}_{13} = (2.5^{+1.8}_{-1.5})\times10^{-2}. This implied an upper bound of sin^2{\theta}_{13} < 0.053 at the 95% confidence level (C.L.).

This paper details the solar neutrino analysis of the 385.17-day Phase-III data set acquired by the Sudbury Neutrino Observatory (SNO). An array of $^3$He proportional counters was installed in the heavy-water target to measure precisely the rate of neutrino-deuteron neutral-current interactions. This technique to determine the total active $^8$B solar neutrino flux was largely independent of the methods employed in previous phases. The total flux of active neutrinos was measured to be $5.54^{+0.33}_{-0.31}(stat.)^{+0.36}_{-0.34}(syst.)\times 10^{6}$ cm$^{-2}$ s$^{-1}$, consistent with previous measurements and standard solar models. A global analysis of solar and reactor neutrino mixing parameters yielded the best-fit values of $\Delta m^2 = 7.59^{+0.19}_{-0.21}\times 10^{-5}{eV}^2$ and $\theta = 34.4^{+1.3}_{-1.2}$ degrees.