Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of $\mathrm{\pm 10~s}$, $\mathrm{\pm 500~s}$, and $\mathrm{\pm 1000~s}$ relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of $(1.13~-~2.44) \times 10^{11}~\rm{cm^{-2}}$ at 5 MeV to $8.0 \times 10^{7}~\rm{cm^{-2}}$ at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be $(5.4~-~7.0)\times 10^{9}~\rm{cm^{-2}}$ for the three time windows.

M31N 2017-09a is a classical nova and was observed for some 160 days following its initial eruption, during which time it underwent a number of bright secondary outbursts. The light-curve is characterized by continual variation with excursions of at least 0.5 magnitudes on a daily time-scale. The lower envelope of the eruption suggests that a single power-law can describe the decline rate. The eruption is relatively long with $t_2 = 111$, and $t_3 = 153$ days.

General relativity is a fully conservative theory, but there exist other possible metric theories of gravity. We consider non-conservative ones with a parameterized post-Newtonian (PPN) parameter, $\zeta_2$. A non-zero $\zeta_2$ induces a self-acceleration for the center of mass of an eccentric binary pulsar system, which contributes to the second time derivative of the pulsar spin frequency, $\ddot{\nu}$. In our work, using the method in Will (1992), we provide an improved analysis with four well-timed, carefully-chosen binary pulsars. In addition, we extend Will's method and derive $\zeta_2$'s effect on the third time derivative of the spin frequency, $\dddot{\nu}$. For PSR B1913+16, the constraint from $\dddot{\nu}$ is even tighter than that from $\ddot{\nu}$. We combine multiple pulsars with Bayesian inference, and obtain an upper limit, $\left|\zeta_{2}\right|<1.3\times10^{-5}$ at 95% confidence level, assuming a flat prior in $\log_{10} \left| \zeta_{2}\right|$. It improves the existing bound by a factor of three. Moreover, we propose an analytical timing formalism for $\zeta_2$. Our simulated times of arrival with simplified assumptions show binary pulsars' capability in limiting $\zeta_{2}$, and useful clues are extracted for real data analysis in future. In particular, we discover that for PSRs B1913+16 and J0737$-$3039A, $\dddot{\nu}$ can yield more constraining limits than $\ddot{\nu}$.

We combine asteroseismology, optical high-resolution spectroscopy, and kinematic analysis for 26 halo red giant branch stars in the \textit{Kepler} field in the range of $-2.5<[\mathrm{{Fe}/{H}}]<-0.6$. After applying theoretically motivated corrections to the seismic scaling relations, we obtain an average mass of $0.96\pm 0.02\,\mathrm{M_{\odot}}$ for our sample of halo stars. This maps into an age of $\sim 8\,\mathrm{Gyr}$, significantly younger than independent age estimates of the Milky Way stellar halo. Although our results suggest that further improvements are needed to accurately determine ages of metal-poor stars, asteroseismology still provides robust relative ages. The lack of a mass dispersion among our sample constrains a relative age dispersion to $<18\%$, corresponding to $<2\,\mathrm{Gyr}$. The precise chemical abundances allow us to separate most of the stars with [{Fe}/{H}]$>-1.7$ into two [{Mg}/{Fe}] groups. While [$\alpha$/{Fe}] and [{Eu}/{Mg}] ratios are different between the two subsamples, [$s$/Eu], where $s$ stands for Ba, La, Ce, and Nd, does not show a significant difference. These abundance ratios suggest that the chemical evolution of the low-Mg population is contributed by type~Ia supernovae, but not by low-to-intermediate mass asymptotic giant branch stars, providing a constraint on its star formation timescale as $100\,\mathrm{Myr}<\tau<300\,\mathrm{Myr}$. We also do not detect any significant mass difference between the two [{Mg}/{Fe}] groups, thus suggesting that their formation epochs are not separated by more than 1.5 Gyr.

The Advanced LIGO and Virgo detectors opened a new era to study black holes (BHs) in our Universe. A population of stellar-mass binary BHs (BBHs) are discovered to be heavier than previously expected. These heavy BBHs provide us an opportunity to achieve multiband observation with ground-based and space-based gravitational-wave (GW) detectors. In this work, we use BBHs discovered by the LIGO/Virgo Collaboration as indubitable examples, and study in great detail the prospects for multiband observation with GW detectors in the near future. We apply the Fisher matrix to spinning, non-precessing inspiral-merger-ringdown waveforms, while taking the motion of space-based GW detectors fully into account. Our analysis shows that, detectors with decihertz sensitivity are expected to log stellar-mass BBH signals with very large signal-to-noise ratio, and provide accurate parameter estimation, including the sky location and time to coalescence. Furthermore, the combination of multiple detectors will achieve unprecedented measurement of BBH properties. As an explicit example, we present the multiband sensitivity to the generic dipole radiation for BHs, which is vastly important for the equivalence principle in the foundation of gravitation, in particular for those theories that predict curvature-induced scalarization of BHs.

We present observations of CO(3-2) and $^{13}$CO(3-2) emission near the supernebula in the dwarf galaxy NGC 5253, which contains one of the best examples of a potential globular cluster in formation. The 0.3" resolution images reveal an unusual molecular cloud, "Cloud D1", coincident with the radio-infrared supernebula. The ~6-pc diameter cloud has a linewidth, $\Delta$ v = 21.7 km/s, that reflects only the gravitational potential of the star cluster residing within it. The corresponding virial mass is 2.5 x 10$^5$ M$_\odot$. The cluster appears to have a top-heavy initial mass function, with $M_{low}$~1-2 M$_\odot$. Cloud D1 is optically thin in CO(3-2) probably because the gas is hot. Molecular gas mass is very uncertain but constitutes < 35% of the dynamical mass within the cloud boundaries. In spite of the presence of an estimated ~1500-2000 O stars within the small cloud, the CO appears relatively undisturbed. We propose that Cloud D1 consists of molecular clumps or cores, possibly star-forming, orbiting with more evolved stars in the core of the giant cluster.

We present a new moving group clustered in kinematics, spatial position and elemental abundances. Its spatial position is around the center of the Local Arm of the Milky Way. A convergent point method was taken to select candidate member stars.\textbf{ Among 206 candidate member stars, 74 are pre-main-sequence stars and some of them have stellar disks.} We presume those pre-main sequence stars belong to Orion nebula. We suggest this moving group is caused by density wave of the Local Arm passing by.

A large number of new members ($\sim$150) of the Cetus Stream (CS) were identified from their clustering features in dynamical space using 6D kinematic data by combining LAMOST DR5 and Gaia DR2 surveys. They map a diffuse structure that extends over at least 100 degrees in the northern and southern Galactic hemispheres, at heliocentric distances between 20 to 50 kpc. Taking advantage of this expanded dataset, we model the stream with a suite of tailored N-body simulations. Our findings exclude the possibility that the NGC 5824 globular cluster is the core of the progenitor of the stream, as postulated by previous studies. Our best models, which successfully reproduce the features of the CS indicate that the progenitor is likely a dwarf galaxy of $\sim$ 2$\times$10$^9$M$_{\odot}$, with a diffuse disc morphology. The merger occured $\sim$ 5 Gyr ago and since then it has experienced approximately eight apo-center passages. Our results suggest that NGC 5824 was either a globular cluster situated off-centre in the dwarf progenitor or, alternatively, it was the nuclear star cluster of another dwarf galaxy that has very similar orbit as the progenitor of the CS. In both scenarios, the progenitor systems would leave streams around NGC 5824, but with distinct distance distributions. To discriminate between these scenarios, the detection and accurate distance measurements of the predicted stream around the GC are crucial, which will be possible in the upcoming LSST era. Our simulations also predict that part of the Southern Cetus stream is very likely the newly discovered Palca stream, and possibly related to another, more diffuse Southern substructure, the Eridanus-Pheonix overdensity.

Tiny violation of Lorentz invariance has been the subject of theoretic study and experimental test for a long time. We use the Standard-Model Extension (SME) framework to investigate the effect of the minimal Lorentz violation on the structure of a neutron star. A set of hydrostatic equations with modifications from Lorentz violation are derived, and then the modifications are isolated and added to the Tolman-Oppenheimer-Volkoff (TOV) equation as the leading-order Lorentz-violation corrections in relativistic systems. A perturbation solution to the leading-order modified TOV equations is found. The quadrupole moments due to the anisotropy in the structure of neutron stars are calculated and used to estimate the quadrupole radiation of a spinning neutron star with the same deformation. The calculation puts forward a new test for Lorentz invariance in the strong-field regime when continuous gravitational waves are observed in the future.

182 single-lined hot subdwarf stars are identified by using spectra from the sixth and seventh data release (DR6 and DR7) of the Large Sky Area Multi-Object Fibre Spectroscopic Telescope (LAMOST) survey. We classified all the hot subdwarf stars using a canonical classification scheme, and got 89 sdB, 37 sdOB, 26 sdO, 24 He-sdOB, 3 He-sdO and 3 He-sdB stars, respectively. Among these stars, 108 hot subdwarfs are newly discovered, while 74 stars were reported by previous catalogs. The atmospheric parameters of these stars were obtained by fitting the hydrogen (H) and helium (He) lines with non-local thermodynamic equilibrium (non-LTE) model atmospheres. The atmospheric parameters confirm the two He sequences and the two subgroups of He-sdOB stars in our samples, which were found by previous studies in the T eff -log(nHe/nH) diagram. Our results demonstrate different origins of field hot subdwarf stars and extreme horizontal branch (EHB) stars in globular clusters (GCs), and provide strict observational limits on the formation and evolution models of the different sub-types of these evolved objects. Based on the results, we evaluated the completeness of the Geier et al. (2019) catalog. We found the fraction of hot subdwarf stars is between 10% and 60%, depending on the brightness of the sample. A more accurate estimation for the hot subdwarf fraction can be obtained when similar results from composite spectra will become available.

Being able to image active regions on the Sun's far side is useful for modeling the global-scale magnetic field around the Sun, and for predicting the arrival of major active regions that rotate around the limb onto the near side. Helioseismic methods have already been developed to image the Sun's far-side active regions using near-side high-cadence Doppler-velocity observations; however, the existing methods primarily explore the 3-, 4-, and 5-skip helioseismic waves, leaving room for further improvement in the imaging quality by including waves with more multi-skip waves. Taking advantage of the facts that 6-skip waves have the same target-annuli geometry as 3- and 4-skip waves, and that 8-skip waves have the same target-annuli geometry as 4-skip waves, we further develop a time--distance helioseismic code to include a total of 14 sets of measurement schemes. We then apply the new code on the SDO/HMI-observed Dopplergrams, and find that the new code provides substantial improvements over the existing codes in mapping newly-emerged active regions and active regions near both far-side limbs. Comparing 3 months of far-side helioseismic images with the STEREO/EUVI-observed 304A images, we find that 97.3% of the helioseismically detected far-side active regions that are larger than a certain size correspond to an observed region with strong EUV brightening. The high reliability of the new imaging tool will potentially allow us to further calibrate the far-side helioseismic images into maps of magnetic flux.

Since Gaia DR2 was released, many velocity structures in the disk have been revealed such as large scale ridge-like patterns in the phase space. Both kinematic information and stellar elemental abundances are needed to reveal their evolution history. We have used labels from the APOGEE survey to predict elemental abundances for a huge amount of low resolution spectra from the LAMOST survey. Deep learning with artificial neural networks can automatically draw on physically sensible features in the spectrum for their predictions. Abundances of 12 individual elements: [C/Fe], [N/Fe], [O/Fe], [Mg/Fe], [Al/Fe], [Si/Fe], [S/Fe], [Cl/Fe], [Ca/Fe], [Ti/Fe], [Mn/Fe] and [Ni/Fe] along with basic stellar labels $T_{\textrm{eff}}$, log $g$, metallicity ([M/H] and [Fe/H]) and [$\alpha$/M] for 1 063 386 stars have been estimated. Then those stars were cross matched with Gaia DR2 data to obtain kinematic parameters. We presented distributions of chemical abundances in the $V_{\phi}$ versus \textit{R} coordinate. Our results extend the chemical characterization of the ridges in the (R, $V_{\phi}$) plane to about $R = $13 kpc toward the anticenter direction. In addition, radial elemental abundance gradients for disk stars with abs(\textit{z})$ < 0.5$ kpc are investigated and we fitted a line for median abundance values of bins of stars with galactocentric distance between $R > 7.84$ kpc and $R < 15.84$ kpc. The radial metallicity gradients for disk stars are respectively -0.0475 $\pm $0.0015 for $R < 13.09$ kpc and -0.0173 $\pm $0.0028 for $R < 13.09$ kpc. Gradients for other elemental abundances are also obtained, for example, [$\alpha$/M] gradients for disk stars is 0.0030 $\pm $0.0002.

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

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

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

We identify member stars of more than 90 open clusters in the LAMOST survey. With the method of Fang et al.(2018), the chromospheric activity (CA) indices logR'CaK for 1091 member stars in 82 open clusters and logR'H{\alpha} for 1118 member stars in 83 open clusters are calculated. The relations between the average logR'CaK, logR'H{\alpha} in each open cluster and its age are investigated in different Teff and [Fe/H] ranges. We find that CA starts to decrease slowly from logt = 6.70 to logt = 8.50, and then decreases rapidly until logt = 9.53. The trend becomes clearer for cooler stars. The quadratic functions between logR' and logt with 4000K < Teff < 5500K are constructed, which can be used to roughly estimate ages of field stars with accuracy about 40% for logR'CaK and 60% for logR'H{\alpha}.

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

We investigate the scalar-tensor gravity of Damour and Esposito-Far\`ese (DEF), which predicts non-trivial phenomena in the nonperturbative strong-field regime for neutron stars (NSs). Instead of solving the modified Tolman-Oppenheimer-Volkoff equations, we construct reduced-order surrogate models, coded in the pySTGROM package, to predict the relations of a NS radius, mass, and effective scalar coupling to its central density. Our models are accurate at $\sim1\%$ level and speed up large-scale calculations by two orders of magnitude. As an application, we use pySTGROM and Markov-chain Monte Carlo techniques to constrain parameters in the DEF theory, with five well-timed binary pulsars, the binary NS (BNS) inspiral GW170817, and a hypothetical BNS inspiral in the Advanced LIGO and next-generation GW detectors. In the future, as more binary pulsars and BNS mergers are detected, our surrogate models will be helpful in constraining strong-field gravity with essential speed and accuracy.

We compared the results from the JSOC pipeline for horizontal flow components and the perturbations of the speed of sound at set of depths with equivalent results from an independently implemented pipeline using a different time-distance inversion scheme. Our inversion pipeline allows inversion for all quantities at once while allowing minimisation of the crosstalk between them. This gives us an opportunity to discuss the possible biases present in the JSOC data products. For the tests we used the subtractive optimally localised averaging (SOLA) method with a minimisation of the cross-talk. We compared three test inversions for each quantity at each target depth. At first, we used the JSOC setup to reproduce the JSOC results. Subsequently, we used the extended pipeline to improve these results by incorporating more independent travel-time measurements but keeping the JSOC-indicated localisation in the Sun. Finally, we inverted for flow components and sound-speed perturbations using a localisation kernel with properties advertised in the JSOC metadata. We successfully reproduced the horizontal flow components. The sound-speed perturbations are strongly affected by the high level of the cross-talk in JSOC products. This leads to larger amplitudes in the inversions for the sound-speed perturbations. Different results were obtained when a target function localised around the target depth was used. This is a consequence of non-localised JSOC averaging kernels. We add that our methodology also allows inversion for the vertical flow.

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

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

388 hot subdwarf stars have been identified by using the Hertzsprung-Russell (HR) diagram built from the second data release (DR2) of the Gaia mission. By analyzing their observed LAMOST spectra, we characterized 186 sdB, 73 He-sdOB, 65 sdOB, 45 sdO, 12 He-sdO and 7 He-sdB stars. The atmospheric parameters of these stars (e.g., T eff , log g, log(nHe/nH)) are obtained by fitting the hydrogen (H) and helium (He) line profiles with synthetic spectra calculated from non-Local Thermodynamic Equilibrium (non-LTE) model atmospheres. Among these stars, we have 135 new identified hot subdwarfs which have not been cataloged before. Although 253 stars appear in the catalog by Geier et al. (2017) , but only 91 of them have atmospheric parameters. Together with the 294 hot subdwarf stars found by Lei et al. (2018), we identified 682 hot subdwarf stars in total by using the Gaia HR-diagram and LAMOST spectra. These results demonstrate the efficiency of our method to combine large surveys to search for hot subdwarf stars. We found a distinct gap in our He-sdOB stars based on their He abundance, which is also presented in extreme horizontal branch (EHB) stars of the globular cluster (GC) {\omega} Cen. The number fraction of the sample size for the two sub-groups is very different between the two counterparts. However, the distinct gap between the H-sdB stars and He-sdOB stars in {\omega} Cen is not visible in our sample. More interestingly, the He-sdB population with the highest He abundance in our sample is completely missing in {\omega} Cen. The discrepancy between our field hot subdwarf stars and the EHB stas in {\omega} Cen indicate different origins for the two counterparts.

Employing a new machine learning method, named hierarchical extreme learning machine (HELM) algorithm, we identified 56 hot subdwarf stars in the first data release (DR1) of the Large Sky Area Multi-Object Fibre Spectroscopic Telescope (LAMOST) survey. The atmospheric parameters of the stars are obtained by fitting the profiles of hydrogen (H) Balmer lines and helium (He) lines with synthetic spectra calculated from non-Local Thermodynamic Equilibrium (NLTE) model atmospheres. Five He-rich hot subdwarf stars were found in our sample with their log(nHe/nH) > -1 , while 51 stars are He-poor sdB, sdO and sdOB stars. We also confirmed the two He sequences of hot subdwarf stars found by Edelmann et al. (2003) in Teff - log(nHe/nH) diagram. The HELM algorithm works directly on the observed spectroscopy and is able to filter out spectral properties without supplementary photometric data. The results presented in this study demonstrate that the HELM algorithm is a reliable method to search for hot subdwarf stars after a suitable training is performed, and it is also suitable to search for other objects which have obvious features in their spectra or images.

We selected 4593 hot subdwarf candidates from the Gaia DR2 Hertzsprung-Russell (HR) diagram. By combining the sample with LAMOST DR5, we identified 294 hot subdwarf stars, including 169 sdB, 63 sdOB, 31 He-sdOB, 22 sdO, 7 He-sdO and 2 He-sdB stars. The atmospheric parameters (e.g., T_{eff} , log g, log(nHe/nH)) are obtained by fitting the hydrogen (H) and helium (He) line profiles with synthetic spectra. Two distinct He sequences of hot subdwarf stars are clearly presented in the T_{eff} - log g diagram. We found that the He-rich sequence consists of the bulk of sdB and sdOB stars as well as all of the He-sdB, He-sdO and He-sdOB stars in our samples, while all the stars in the He-weak sequence belong to the sdO spectral type, combined with a few sdB and sdOB stars. We demonstrated that the combination of Gaia DR2 and LAMOST DR5 allows one to uncover a huge number of new hot subdwarf stars in our Galaxy.

In this paper we report on the characterization of the Hamamatsu VUV4 (S/N: S13370-6152) Vacuum Ultra-Violet (VUV) sensitive Silicon Photo-Multipliers (SiPMs) as part of the development of a solution for the detection of liquid xenon scintillation light for the nEXO experiment. Various SiPM features, such as: dark noise, gain, correlated avalanches, direct crosstalk and Photon Detection Efficiency (PDE) were measured in a dedicated setup at TRIUMF. SiPMs were characterized in the range $163 \text{ } \text{K} \leq \text{T}\leq 233 \text{ } \text{K}$. At an over voltage of $3.1\pm0.2$ V and at $\text{T}=163 \text{ }\text{K}$ we report a number of Correlated Avalanches (CAs) per pulse in the $1 \upmu\text{s}$ interval following the trigger pulse of $0.161\pm0.005$. At the same settings the Dark-Noise (DN) rate is $0.137\pm0.002 \text{ Hz/mm}^{2}$. Both the number of CAs and the DN rate are within nEXO specifications. The PDE of the Hamamatsu VUV4 was measured for two different devices at $\text{T}=233 \text{ }\text{K}$ for a mean wavelength of $189\pm7\text{ nm}$. At $3.6\pm0.2$ V and $3.5\pm0.2$ V of over voltage we report a PDE of $13.4\pm2.6\text{ }\%$ and $11\pm2\%$, corresponding to a saturation PDE of $14.8\pm2.8\text{ }\%$ and $12.2\pm2.3\%$, respectively. Both values are well below the $24\text{ }\%$ saturation PDE advertised by Hamamatsu. More generally, the second device tested at $3.5\pm0.2$ V of over voltage is below the nEXO PDE requirement. The first one instead yields a PDE that is marginally close to meeting the nEXO specifications. This suggests that with modest improvements the Hamamatsu VUV4 MPPCs could be considered as an alternative to the FBK-LF SiPMs for the final design of the nEXO detector.

Using JVLA data of Sgr A* and NRAO 530 (as a calibrator) at 5.5, 9 and 33 GHz, we developed a procedure for the reduction of wideband data. We have demonstrated that, correcting for residual interferometer errors such as residual delays, astronomers can now achieve high-fidelity radio images with a dynamic range exceeding 1,000,000:1. We outline the detailed procedure, noting that it can have broad application to the analysis of broadband continuum observations. We apply this procedure to observations of a distant blazar, NRAO 530, revealing its structures in unprecedented detail. Our 5.5-GHz image shows that the structure of NRAO 530 is prominently characterized by a moderately curved western jet terminating at a hot spot. Close to the radio core, an abrupt bending of the jet is revealed in the high-resolution (<100 mas) images at 33 GHz, showing an evolution of the position angle of the jet from the north at the VLBI scale (50 mas, or ~400 pc projected), increasing toward the west at the larger VLA scales (1 arcsec, or ~10 kpc). The continuation of the jet axis drift forms the curved western jet extending out to 20 arcsec. In contrast, a faint and broad counter-jet is present on the eastern side with a curvature antisymmetric to the western jet. The eastern jet terminates at a bright hotspot, forming an edge-brightened diffuse lobe. The observed contrast in brightness between the western and eastern jets suggests that the jets on the VLA scale are mildly relativistic. The radiation from the western jet is boosted while the radiation from the receding eastern jet is plausibly suppressed owing to the relativistic Doppler effect.

We write in response to the call from the 2020 Decadal Survey to submit white papers illustrating the most pressing scientific questions in astrophysics for the coming decade. We propose exploration as the central question for the Decadal Committee's discussions.The history of astronomy shows that paradigm changing discoveries are not driven by well formulated scientific questions, based on the knowledge of the time. They were instead the result of the increase in discovery space fostered by new telescopes and instruments. An additional tool for increasing the discovery space is provided by the analysis and mining of the increasingly larger amount of archival data available to astronomers. Revolutionary observing facilities, and the state of the art astronomy archives needed to support these facilities, will open up the universe to new discovery. Here we focus on exploration for compact objects and multi messenger science. This white paper includes science examples of the power of the discovery approach, encompassing all the areas of astrophysics covered by the 2020 Decadal Survey.

We have observed the active star $\xi$ Boo A (HD 131156A) with high precision broadband linear polarimetry contemporaneously with circular spectropolarimetry. We find both signals are modulated by the 6.43 day rotation period of $\xi$ Boo A. The signals from the two techniques are 0.25 out of phase, consistent with the broadband linear polarization resulting from differential saturation of spectral lines in the global transverse magnetic field. The mean magnitude of the linear polarization signal is ~4 ppm/G but its structure is complex and the amplitude of the variations suppressed relative to the longitudinal magnetic field. The result has important implications for current attempts to detect polarized light from hot Jupiters orbiting active stars in the combined light of the star and planet. In such work stellar activity will manifest as noise, both on the time scale of stellar rotation, and on longer time scales - where changes in activity level will manifest as a baseline shift between observing runs.

Based on the [Fe/H] versus [Mg/Fe] diagram and distances from APOGEE data release 14, we compare the spatial distributions, the l-Vlos diagram and the abundance gradients between high-[Mg/Fe] and low-[Mg/Fe] sequences. The two sequences are clearly shown at 5<|Z|<10 kpc in the metallicity range of -1.6 <[Fe/H] <-0.7, where the halo at |Z| > 10 kpc consists of low-[Mg/Fe] stars only. In the intermediate-metallicity range of -1.1 <[Fe/H]<-0.7, a [Mg/Fe] gradient is detected for stars at |Z|=10-30 kpc and it flattens out at |Z|>30 kpc. The l-Vlos diagram is adopted to separate halo stars from the disk by defining the transition metallicity, which is of [Fe/H]~ -1.1 dex for the high-[Mg/Fe] sequence and of [Fe/H]~-0.7 dex for the low-[Mg/Fe] sequence. The R and |Z| distributions for the high-[Mg/Fe] sequence, the thick disk at -1.1<[Fe/H]<-0.7 and the in situ halo at -1.6<[Fe/H]<-1.1, have a cutoff at R~15 kpc and |Z|~10 kpc, beyond which low-[Mg/Fe] halo stars are the main contributions. In the metallicity range of -1.6<[Fe/H]<-0.7, there is a negative metallicity gradient for the high-[Mg/Fe] halo at |Z|<8-10 kpc, while only a marginal or no slope in the [Fe/H] versus |Z| diagram for the low-[Mg/Fe] halo at |Z|<8-10 kpc, beyond which both the high-[Mg/Fe] halo and low-[Mg/Fe] halo flatten out toward |Z| > 20 kpc. These results indicate a complicated formation history of the Galaxy and we may see a hint of a three-section halo, i.e. the inner in situ halo within $|Z|~8-10$ kpc, the intermediately outer dual-mode halo at |Z|~10-30 kpc, and the extremely outer accreted halo with |Z|>30 kpc.

We examine eight known single-eccentric planetary systems in light of recently released large data archives and new analysis techniques. For four of these systems (HD 7449, HD 65216, HD 89744, HD 92788) we find evidence for additional long-period companions. HD 65216c is a Jupiter analog, with a period of 14.7 yr, $e=0.18$, and m sin $i$ of 2M_Jup, whilst the remaining candidate companions move on as-yet-incomplete orbits. Our results highlight the importance of revisiting the analysis of known exoplanetary systems when new data become available, particularly given the possibility that poorly-sampled data might previously have led to the detection of a 'false-positive' single eccentric planet, when the system in question actually contains two (or more) planets on near-circular orbits.

The hadronic interaction of cosmic rays with solar atmosphere can produce high energy gamma rays. The gamma-ray luminosity is correlated both with the flux of primary cosmic rays and the intensity of the solar magnetic field. The gamma rays below 200 GeV have been observed by $Fermi$ without any evident energy cutoff. The bright gamma-ray flux above 100 GeV has been detected only during solar minimum. The only available data in TeV range come from the HAWC observations, however outside the solar minimum. The ARGO-YBJ dataset has been used to search for sub-TeV/TeV gamma rays from the Sun during the solar minimum from 2008 to 2010, the same time period covered by the Fermi data. A suitable model containing the Sun shadow, solar disk emission and inverse-Compton emission has been developed, and the chi-square minimization method was used to quantitatively estimate the disk gamma-ray signal. The result shows that no significant gamma-ray signal is detected and upper limits to the gamma-ray flux at 0.3$-$7 TeV are set at 95\% confidence level. In the low energy range these limits are consistent with the extrapolation of the Fermi-LAT measurements taken during solar minimum and are compatible with a softening of the gamma-ray spectrum below 1 TeV. They provide also an experimental upper bound to any solar disk emission at TeV energies. Models of dark matter annihilation via long-lived mediators predicting gamma-ray fluxes > $10^{-7}$ GeV $cm^{-2}$ $s^{-1}$ below 1 TeV are ruled out by the ARGO-YBJ limits.

We study the radial metallicity gradient $\Delta[M/H]/\Delta R_g$ as a function of [Mg/Fe] and $|Z|$ with the help of a guiding radius based on the Apache Point Observatory Galactic Evolution Experiment and Gaia and then analyze the radial migration effect on the radial metallicity gradient and metallicity-rotation gradient between the Galactic thin and thick disks. The derived trend of gradient $\Delta[M/H]/\Delta R_g$ versus [Mg/Fe] shows a transition at [Mg/Fe] $\sim 0.18$ dex, below which the gradient is negative and varies a little as [Mg/Fe] increases; however, it changes sharply in [Mg/Fe] ranges of 0.16-0.18, above which the gradient increases linearly with increasing [Mg/Fe], being a positive value at [Mg/Fe]$\gtrsim 0.22$ dex. These positive gradients in the high-[Mg/Fe] populations are found at $|Z| < 0.8$ kpc, and there are nearly no gradients toward higher $|Z|$. By comparing the metallicity distributions, the radial metallicity gradients $\Delta[M/H]/\Delta R$ and the metallicity-rotation gradients between the total sample and $|R-R_g|<2$ kpc subsample (or $|R-R_g|>2$ kpc subsample), we find that, for the thick disk, blurring flattens the gradient $\Delta[M/H]/\Delta R$ and favors metal-poor high-eccentricity stars. These stars are responsible for the measured positive metallicity-rotation gradient of the thick disk.

We present the X-ray timing results of the new black hole candidate (BHC) MAXI J1535-571 during its 2017 outburst from Hard X-ray Modulation Telescope (\emph{Insight}-HXMT) observations taken from 2017 September 6 to 23. Following the definitions given by \citet{Belloni2010}, we find that the source exhibits state transitions from Low/Hard state (LHS) to Hard Intermediate state (HIMS) and eventually to Soft Intermediate state (SIMS). Quasi-periodic oscillations (QPOs) are found in the intermediate states, which suggest different types of QPOs. With the large effective area of \emph{Insight}-HXMT at high energies, we are able to present the energy dependence of the QPO amplitude and centroid frequency up to 100 keV which is rarely explored by previous satellites. We also find that the phase lag at the type-C QPOs centroid frequency is negative (soft lags) and strongly correlated with the centroid frequency. By assuming a geometrical origin of type-C QPOs, the source is consistent with being a high inclination system.

We present chemical abundances of 15 stars in the {\gamma} Leo moving group based on high-resolution spectra with the Subaru High Dispersion Spectrograph. The sample was picked up by applying wavelet transform to UVW velocity compo- nents of stars in the solar neighbourhood. Both photometric and spectroscopic method have been used to determine the stellar parameters of stars. Abun- dances of 11 elements including Na, Mg, Al, Si, Ca, Ti, Cr, Fe, Ni, Y and Ba are measured. Our results show that the member stars display a wide metallicity distribution with abundance ratios similar to Milky way disk stars. We presume that the {\gamma} Leo moving group is originated from dynamical effects probably related to the Galactic spiral arms.

In this paper, we report the discovery of a new late-B type unbound hyper-runaway star (LAMOST-HVS4) from the LAMOST spectroscopic surveys. According to its atmospheric parameters, it is either a B-type main sequence (MS) star or a blue horizontal branch (BHB) star. Its Galactocentric distance and velocity are 30.3 +/- 1.6 kpc and 586 +/- 7 km/s if it is an MS star, and they are 13.2 +/- 3.7 kpc and 590 +/- 7 km/s if a BHB star. We track its trajectories back, and find that the trajectories intersect with the Galactic disk and the Galactic center lies outside of the intersection region at the 3 sigma confidence level. We investigate a number of mechanisms that could be responsible for the ejection of the star, and find that it is probably ejected from the Galactic disk by supernova explosion or multiple-body interactions in dense young stellar clusters.

Using a machine learning (ML) method, we mine DB white dwarfs (DBWDs) from the Sloan Digital Sky Survey (SDSS) Data Release (DR) 12 and DR14. The ML method consists of two parts: feature extraction and classification. The least absolute shrinkage and selection operator (LASSO) is used for the spectral feature extraction by comparing high quality data of a positive sample group with negative sample groups. In both the training and testing sets, the positive sample group is composed of a selection of 300 known DBWDs, while the negative sample groups are obtained from all types of SDSS spectra. In the space of the LASSO detected features, a support vector machine is then employed to build classifiers that are used to separate the DBWDs from the non DBWDs for each individual type. Depending on the classifiers, the DBWD candidates are selected from the entire SDSS dataset. After visual inspection, 2808 spectra (2029 objects) are spectroscopically confirmed. By checking the samples with the literature, there are 58 objects with 60 spectra that are newly identified, including a newly discovered AM CVn. Finally, we measure their effective temperatures (Teff), surface gravities (log g), and radial velocities, before compiling them into a catalog.

We present the first results from a joint XMM-Newton and NuSTAR observation of the gamma-ray emitting millisecond pulsar compact binary PSR J2129-0429. X-ray emission up to about 40 keV is detected and the joint spectrum can be modeled with a power-law plus a neutron star atmosphere model. At a distance of 1.4 kpc, the 0.3-79 keV luminosity is $3.56\times10^{32}$ erg s$^{-1}$. We also detected the 0.64-day binary orbital period with a double-peaked structure across the wavebands. By combining the updated Fermi GeV data, we modeled the broadband spectral energy distribution as well as the X-ray modulation with an intrabinary model involving shock interaction between pulsar wind and outflow from the companion star. Lastly, we reported a high-resolution X-ray image provided by Chandra to rule out the proposed pulsar wind nebula associated with PSR J2129-0429.

In the minimal supersymmetric standard model (MSSM) the lightest superparticle (LSP) can be a TeV neutralino (mainly dominated by higgsino or wino) which serves as a dark matter candidate with correct thermal relic density. In this work we confront the 1-2 TeV neutralino dark matter with the latest direct and indirect detections from PandaX and AMS-02/DAMPE. Considering various scenarios with decoupled sfermions, with \textit{A}-mediated annihilation, with squark or stop coannihilation, we find that the parameter space is stringently constrained by the direct detection limits. In the allowed parameter space, the TeV neutralino dark matter annihilation contribution to the anti-proton flux is found to agree with the AMS-02 data while its contribution to eletron/positron flux is too small to cause any visible excess. The current survived parameter space can be mostly covered by the future direct detection experiment LZ7.2T.

This paper reports on the observation of the sidereal large-scale anisotropy of cosmic rays using data collected by the ARGO-YBJ experiment over 5 years (2008$-$2012). This analysis extends previous work limited to the period from 2008 January to 2009 December,near the minimum of solar activity between cycles 23 and 24.With the new data sample the period of solar cycle 24 from near minimum to maximum is investigated. A new method is used to improve the energy reconstruction, allowing us to cover a much wider energy range, from 4 to 520 TeV. Below 100 TeV, the anisotropy is dominated by two wide regions, the so-called "tail-in" and "loss-cone" features. At higher energies, a dramatic change of the morphology is confirmed. The yearly time dependence of the anisotropy is investigated. Finally, no noticeable variation of cosmic-ray anisotropy with solar activity is observed for a median energy of 7 TeV.

PSR~J2021+4026 showed a sudden decrease in the gamma-ray emission at the glitch that occurred around 2011, October 16, and a relaxation of the flux to the pre-glitch state at around 2014 December. We report X-ray analysis results of the data observed by XMM-Newton on 2015 December 20 in the post-relaxation state. To examine any change in the X-ray emission, we compare the properties of the pulse profiles and spectra at the low gamma-ray flux state and at the post-relaxation state. The phase-averaged spectra for both states can be well described by a power-law component plus a blackbody component. The former is dominated by unpulsed emission and is probably originated from the pulsar wind nebula as reported by Hui et al (2015). The emission property of the blackbody component is consistent with the emission from the polar cap heated by the back-flow bombardment of the high-energy electrons or positrons that were accelerated in the magnetosphere. We found no significant change in the X-ray emission properties between two states. We suggest that the change of the X-ray luminosity is at an order of ~4%, which is difficult to measure with the current observations. We model the observed X-ray light curve with the heated polar cap emission and we speculate that the observed large pulsed fraction is owing to asymmetric magnetospheric structure.

A correlation between the secondary cosmic ray flux and the near-earth electric field intensity, measured during thunderstorms, has been found by analyzing the data of the ARGO-YBJ experiment, a full coverage air shower array located at the Yangbajing Cosmic Ray Laboratory (4300 m a. s. l., Tibet, China). The counting rates of showers with different particle multiplicities, have been found to be strongly dependent upon the intensity and polarity of the electric field measured during the course of 15 thunderstorms. In negative electric fields (i.e. accelerating negative charges downwards), the counting rates increase with increasing electric field strength. In positive fields, the rates decrease with field intensity until a certain value of the field EFmin (whose value depends on the event multiplicity), above which the rates begin increasing. By using Monte Carlo simulations, we found that this peculiar behavior can be well described by the presence of an electric field in a layer of thickness of a few hundred meters in the atmosphere above the detector, which accelerates/decelerates the secondary shower particles of opposite charge, modifying the number of particles with energy exceeding the detector threshold. These results, for the first time, give a consistent explanation for the origin of the variation of the electron/positron flux observed for decades by high altitude cosmic ray detectors during thunderstorms.

Study of primary stars lying in Sirius-like systems with various masses of WD companions and orbital separations is one of the key aspects to understand the origin and nature of Barium (Ba) stars. In this paper, based on high resolution and high S/N spectra, we present systematic analysis of photospheric abundances for 18 FGK primary stars of Sirius-like systems including six giants and 12 dwarfs. Atmospheric parameters, stellar masses, and abundances of 24 elements (C, Na, Mg, Al, Si, S, K, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Sr, Y, Zr, Ba, La, Ce and Nd) are determined homogeneously. The abundance patterns in these sample stars show that most of the elements in our sample follow the behavior of field stars with similar metallicity. As expected, s-process elements in four known Ba giants show overabundance. A weak correlation was found between anomalies of s-process elemental abundance and orbital separation, suggesting the orbital separation of the binaries could not be the main constraint to differentiate strong Ba stars from mild Ba stars. Our study shows that the large mass (>0.51 M ) of a WD companion in a binary system is not a sufficient condition to form a Ba star, even if the separation between the two components is small. Although not sufficient it seems to be a necessary condition since Ba stars with lower mass WDs in the observed sample were not found. Our results support that [s/Fe] and [hs/ls] ratios of Ba stars are anti-correlated with the metallicity. However, the different levels of s-process overabundance among Ba stars may not to be dominated mainly by the metallicity.

We use the LAMOST spectra of member stars in Pleiades, M34, Praesepe, and Hyades to study how chromospheric activity vary as a function of mass and rotation at different age. We measured excess equivalent widths of H$\alpha$, H$\beta$, and Ca~{\sc ii} K based on estimated chromospheric contributions from old and inactive field dwarfs, and excess luminosities are obtained by normalizing bolometric luminosity, for more than 700 late-type stars in these open clusters. Results indicate two activity sequences in cool spot coverage and H$\alpha$ excess emission among GK dwarfs in Pleiades and M dwarfs in Praesepe and Hyades, paralleling with well known rotation sequences. A weak dependence of chromospheric emission on rotation exists among ultra fast rotators in saturated regime with Rossby number Ro$\lesssim0.1$. In the unsaturated regime, chromospheric and coronal emission show similar dependence on Ro, but with a shift toward larger Ro, indicating chromospheric emission gets easily saturated than coronal emission, and/or convective turnover time-scales based on X-ray data do not work well with chromospheric emission. More interestingly, our analysis show fully convective slow rotators obey the rotation-chromospheric activity relation similar to hotter stars, confirming the previous finding. We found correlations among H$\alpha$, H$\beta$, and Ca~{\sc ii} K emissions, in which H$\alpha$ losses are more important than Ca~{\sc ii} K for cooler and more active stars. In addition, a weak correlation is seen between chromospheric emission and photospheric activity that shows dependency on stellar spectral type and activity level, which provides some clues on how spot configuration vary as a function of mass and activity level.

On 2016 September 20, the Interface Region Imaging Spectrograph observed an active region during its earliest emerging phase for almost 7 hours. The Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory observed continuous emergence of small-scale magnetic bipoles with a rate of $\sim$10$^{16}$ Mx~s$^{-1}$. The emergence of magnetic fluxes and interactions between different polarities lead to frequent occurrence of ultraviolet (UV) bursts, which exhibit as intense transient brightenings in the 1400 \AA{} images. In the meantime, discrete small patches with the same magnetic polarity tend to move together and merge, leading to enhancement of the magnetic fields and thus formation of pores (small sunspots) at some locations. The spectra of these UV bursts are characterized by the superposition of several chromospheric absorption lines on the greatly broadened profiles of some emission lines formed at typical transition region temperatures, suggesting heating of the local materials to a few tens of thousands of kelvin in the lower atmosphere by magnetic reconnection. Some bursts reveal blue and red shifts of $\sim$100~km~s$^{-1}$ at neighboring pixels, indicating the spatially resolved bidirectional reconnection outflows. Many such bursts appear to be associated with the cancellation of magnetic fluxes with a rate of the order of $\sim$10$^{15}$ Mx~s$^{-1}$. We also investigate the three-dimensional magnetic field topology through a magneto-hydrostatic model and find that a small fraction of the bursts are associated with bald patches (magnetic dips). Finally, we find that almost all bursts are located in regions of large squashing factor at the height of $\sim$1 Mm, reinforcing our conclusion that these bursts are produced through reconnection in the lower atmosphere.

Given the recent constraints from the dark matter (DM) direct detections, we examine a light GeV-scale (2-30 GeV) neutralino DM in the alignment limit of the Minimal Supersymmetric Standard Model (MSSM). In this limit without decoupling, the heavy CP-even scalar $H$ plays the role of the Standard Model (SM) Higgs boson while the other scalar $h$ can be rather light so that the DM can annihilate through the $h$ resonance or into a pair of $h$ to achieve the observed relic density. With the current collider and cosmological constraints, we find that such a light neutralino DM above 6 GeV can be excluded by the XENON-1T (2017) limits while the survivied parameter space below 6 GeV can be fully covered by the future germanium-based light dark matter detections (such as CDEX), by the Higgs coupling precison measurements or by the production process $e^+e^- \to hA$ at an electron-positron collider (Higgs factory).

We scrutinize the parameter space of the SM-like Higgs boson in the minimal supersymmetric standard model (MSSM) under current experimental constraints. The constraints are from (i) the precison electroweak data and various flavor observables; (ii) the direct 22 separate ATLAS searches in Run-1; (iii) the latest LHC Run-2 Higgs data and tri-lepton search of electroweakinos. We perform a scan over the parameter space and find that the Run-2 data can further exclude a part of parameter space. For the property of the SM-like Higgs boson, its gauge couplings further approach to the SM values with a deviation below 0.1\%, while its Yukawa couplings $hb\bar{b}$ and $h\tau\bar\tau$ can still sizably differ from the SM predictions by several tens percent.

Radial velocity observations from three instruments reveal the presence of a 4 M_jup planet candidate orbiting the K giant HD 76920. HD 76920b has an orbital eccentricity of 0.856$\pm$0.009, making it the most eccentric planet known to orbit an evolved star. There is no indication that HD 76920 has an unseen binary companion, suggesting a scattering event rather than Kozai oscillations as a probable culprit for the observed eccentricity. The candidate planet currently approaches to about four stellar radii from its host star, and is predicted to be engulfed on a $\sim$100 Myr timescale due to the combined effects of stellar evolution and tidal interactions.

Using the Very Large Array, we have investigated a non-thermal radio filament (NTF) recently found very near the Galactic black hole and its radio counterpart, SgrA*. While this NTF -- the Sgr A West Filament (SgrAWF) -- shares many characteristics with the population of NTFs occupying the central few hundred parsecs of the Galaxy, the SgrAWF has the distinction of having an orientation and sky location that suggest an intimate physical connection to SgrA*. We present 3.3 and 5.5 cm images constructed using an innovative methodology that yields a very high dynamic range, providing an unprecedentedly clear picture of the SgrAWF. While the physical association of the SgrAWF with SgrA* is not unambiguous, the images decidedly evoke this interesting possibility. Assuming that the SgrAWF bears a physical relationship to SgrA*, we examine the potential implications. One is that SgrA* is a source of relativistic particles constrained to diffuse along ordered local field lines. The relativistic particles could also be fed into the local field by a collimated outflow from SgrA*, perhaps driven by the Poynting flux accompanying the black hole spin in the presence of a magnetic field threading the event horizon. Second, we consider the possibility that the SgrAWF is the manifestation of a low-mass-density cosmic string that has become anchored to the black hole. The simplest form of these hypotheses would predict that the filament be bi-directional, whereas the SgrAWF is only seen on one side of SgrA*, perhaps because of the dynamics of the local medium.

We report three new barium (Ba) dwarfs lying in Sirius-like systems, which provides direct evidence that Ba dwarfs are companions to white dwarfs (WDs). Atmospheric parameters, stellar masses, and chemical abundances of 25 elements, including light, $\alpha$, Fe-peak and s-process elements, are derived from high resolution and high S/N spectra. Enhancement of s-process elements with [s/Fe] ratios between 0.4 and 0.6 confirm them as mild barium stars. The estimated metallicities ($-$0.31, $-$0.06, 0.13) of BD+68$^\circ$1027, RE~J0702+129 and BD+80$^\circ$670 are in the range of known Ba dwarfs and giants. As expected, observed indices of [hs/ls], [s/Fe] and [C/Fe] show anticorrelation with metallicity. AGB progenitor masses are estimated for the WD companions of RE~J0702+129 (1.47 $M_{\odot}$) and BD+80$^\circ$670 (3.59 $M_{\odot}$), which confirms the predicted range of progenitor AGB masses (1.5 $\sim$ 4 $M_{\odot}$) for unseen WDs around Ba dwarfs. Surface abundances of s-process elements in RE~J0702+129 and BD+80$^\circ$670 are compared with AGB models and they are in close agreement, within predicted accretion efficiencies and pollution factors for Ba stars. These results support that the origin of s-process overabundances in Ba dwarfs is similar to Ba giants via McClure hypothesis in which Ba stars accumulate s-process elements through mass transfer from their host companions during AGB phase.

We report ALMA observations of $^{12}$CO(3-2) and $^{13}$CO(3-2) in the gas-poor dwarf galaxy NGC 5253. These 0.3"(5.5 pc) resolution images reveal small, dense molecular gas clouds that are located in kinematically distinct, extended filaments. Some of the filaments appear to be falling into the galaxy and may be fueling its current star formation. The most intense CO(3-2) emission comes from the central $\sim$100 pc region centered on the luminous radio-infrared HII region known as the supernebula. The CO(3-2) clumps within the starburst region are anti-correlated with H$\alpha$ on $\sim$5 pc scales, but are well-correlated with radio free-free emission. Cloud D1, which enshrouds the supernebula, has a high $^{12}$CO/$^{13}$CO ratio, as does another cloud within the central 100 pc starburst region, possibly because the clouds are hot. CO(3-2) emission alone does not allow determination of cloud masses as molecular gas temperature and column density are degenerate at the observed brightness, unless combined with other lines such as $^{13}$CO.

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

Meridional circulation is a crucial component of the Sun's internal dynamics, but its inference in the deep interior is complicated by a systematic center-to-limb effect in helioseismic measurement techniques. Previously, an empirical method, removing travel-time shifts measured for east-west traveling waves in the equatorial area from those measured for north-south traveling waves in the central meridian area, was used, but its validity and accuracy need to be assessed. Here we develop a new method to separate the center-to-limb effect, \delta\tau_CtoL, and meridional-flow-induced travel-time shifts, \delta\tau_MF, in a more robust way. Using 7-yr observations of the SDO/HMI, we exhaustively measure travel-time shifts between two surface locations along the solar disk's radial direction for all azimuthal angles and all skip distances. The measured travel-time shifts are a linear combination of \delta\tau_CtoL and \delta\tau_MF, which can be disentangled through solving the linear equation set. The \delta\tau_CtoL is found isotropic relative to the azimuthal angle, and the \delta\tau_MF are then inverted for the meridional circulation. Our inversion results show a three-layer flow structure, with equatorward flow found between about 0.82 and 0.91 R_sun for low latitude areas and between about 0.85 and 0.91 R_sun for higher latitude areas. Poleward flows are found below and above the equatorward flow zones, indicating a double-cell circulation in each hemisphere.

Based on the effective field theory (EFT) of nonsingular cosmologies, we build a stable model, without the ghost and gradient instabilities, of bounce inflation (inflation is preceded by a cosmological bounce). We perform a full simulation for the evolution of scalar perturbation, and find that the perturbation spectrum has a large-scale suppression (as expected), which is consistent with the power deficit of the cosmic microwave background (CMB) TT-spectrum at low multipoles, but unexpectedly, it also shows itself one marked lower valley, which actually provides a better fit to the dip at multipole $l\sim 20$. The depth of valley is relevant with the physics around the bounce scale, which is model-dependent.

We have cross-matched the LAMOST DR3 with the Gaia DR1 TGAS catalogs and obtained a sample of 166,827 stars with reliable kinematics. A technique based on the wavelet transform was applied to detect significant overdensities in velocity space among five subsamples divided by spatial position. In total, 16 significant overdensities of stars with very similar kinematics were identified. Among these, four are new stream candidates and the rest are previously known groups. Both the U-V velocity and metallicity distributions of the local sample show a clear gap between the Hercules structure and the Hyades-Pleiades structure. The U-V positions of these peaks shift with the spatial position. Following a description of our analysis, we speculate on possible origins of our stream candidates.

The ARGO-YBJ experiment, a full coverage extensive air shower (EAS) detector located at high altitude (4300 m a.s.l.) in Tibet, China, has smoothly taken data, with very high stability, since November 2007 to the beginning of 2013. The array consisted of a carpet of about 7000 m$^2$ Resistive Plate Chambers (RPCs) operated in streamer mode and equipped with both digital and analog readout, providing the measurement of particle densities up to few particles per cm$^2$. The unique detector features (full coverage, readout granularity, wide dynamic range, etc) and location (very high altitude) allowed a detailed study of the lateral density profile of charged particles at ground very close to the shower axis and its description by a proper lateral distribution function (LDF). In particular, the information collected in the first 10 m from the shower axis have been shown to provide a very effective tool for the determination of the shower development stage ("age") in the energy range 50 TeV - 10 PeV. The sensitivity of the age parameter to the mass composition of primary Cosmic Rays is also discussed.

A glitch of a pulsar is known as a sudden increase in the spin frequency and spin-down rate (frequency time derivative), and it can be caused by a sudden rel\ ease of the stress built up in the solid crust of the star or pinned vortices in the superfluid interior. PSR J2021+4026 is the first pulsar that shows a significant change in the gamma-ray flux and pulse profile at the glitch that occurred around 2011 October 16. We report the results of timing and spectral analysis of PSR~J2021+4026 using $\sim$ 8~yr Fermi-LAT data. We find that the pulsar stayed at a high spin-down rate ($\sim 4\%$ higher than the pre-glitch value) and a low gamma-ray state ($\sim 18\%$ lower) for about 3~yr after the glitch. Around 2014 December, the spin-down rate and gamma-ray flux gradually returned to pre-glitch values within a time scale of a few months. The phase-resolved spectra and pulse profiles after the relaxation are also consistent with those before the glitch. The observed long-term evolution of the spin-down rate and the gamma-ray flux indicates that the glitch triggered a mode change in the global magnetosphere. We speculate that the glitch changed the local magnetic field structure around the polar cap and/or the inclination angle of the dipole axis, leading to a change in the electric current circulating in the magnetosphere.

Red clump stars are fundamental distance indicators in astrophysics, although theoretical stellar models predict a dependence of absolute magnitudes with ages. This effect is particularly strong below 2 Gyr, but even above this limit a mild age dependence is still expected. We use seismically identified red clump stars in the Kepler field for which we have reliable distances, masses and ages from the SAGA survey to first explore this effect. By excluding red clump stars with masses larger than 1.6 Msun (corresponding to ages younger than 2 Gyr), we derive robust calibrations linking intrinsic colors to absolute magnitudes in the following photometric systems: Str\"omgren $by$, Johnson $BV$, Sloan $griz$, 2MASS $JHK_s$ and WISE $W1W2W3$. With the precision achieved we also detect a slope of absolute magnitudes 0.020(0.003) mag per Gyrin the infrared, implying that distance calibrations of clump stars can be off by up to 0.2 mag in the infrared (over the range from 2 Gyr to 12 Gyr) if their ages are unknown. Even larger uncertainties affect optical bands, because of the stronger interdependency of absolute magnitudes on colors and age. Our distance calibrations are ultimately based on asteroseismology, and we show how the distance scale can be used to test the accuracy of seismic scaling relations. Within the uncertainties our calibrations are in agreement with those built upon local red clump with Hipparcos} parallaxes, although we find a tension which if confirmed would imply that scaling relations overestimate radii of red clump stars by 2(+-20%. Data-releases post Gaia DR1 will provide an important testbed for our results.

The ARGO-YBJ detector, located at the Yangbajing Cosmic Ray Laboratory (4300 m a. s. l., Tibet, China), was a full coverage air shower array dedicated to gamma ray astronomy and cosmic ray studies. The wide field of view (~ 2 sr) and high duty cycle (> 86%), made ARGO-YBJ suitable to search for short and unexpected gamma ray emissions like gamma ray bursts (GRBs). Between 2007 November 6 and 2013 February 7, 156 satellite-triggered GRBs (24 of them with known redshift) occurred within the ARGO-YBJ field of view. A search for possible emission associated to these GRBs has been made in the two energy ranges 10-100 GeV and 10-1000 GeV. No significant excess has been found in time coincidence with the satellite detections nor in a time window of one hour after the bursts. Taking into account the EBL absorption, upper limits to the energy fluence at 99% of confidence level have been evaluated,with values ranging from ~ 10-5 erg cm-2 to ~10-1 erg cm-2.

By using a new method of forced-field extrapolation, we study the emerging flux region AR 11850 observed by the Interface Region Imaging Spectrograph (IRIS) and Solar Dynamical Observatory (SDO). Our results suggest that the bright points (BPs) in this emerging region have responses in lines formed from the upper photosphere to the transition region, with a relatively similar morphology. They have an oscillation of several minutes according to the Atmospheric Imaging Assembly (AIA) data at 1600 and 1700 A . The ratio between the BP intensities measured in 1600 A and 1700 A filtergrams reveals that these BPs are heated differently. Our analysis of the Helioseismic and Magnetic Imager (HMI) vector magnetic field and the corresponding topology in AR11850 indicates that the BPs are located at the polarity inversion line (PIL) and most of them related with magnetic reconnection or cancelation. The heating of the BPs might be different due to different magnetic topology. We find that the heating due to the magnetic cancelation would be stronger than the case of bald patch reconnection. The plasma density rather than the magnetic field strength could play a dominant role in this process. Based on physical conditions in the lower atmosphere, our forced-field extrapolation shows consistent results between the bright arcades visible in slit-jaw image (SJI) 1400 A and the extrapolated field lines that pass through the bald patches. It provides a reliable observational evidence for testing the mechanism of magnetic reconnection for the BPs and arcades in emerging flux region, as proposed in simulation works.

We present the second release of value-added catalogues of the LAMOST Spectroscopic Survey of the Galactic Anticentre (LSS-GAC DR2). The catalogues present values of radial velocity $V_{\rm r}$, atmospheric parameters --- effective temperature $T_{\rm eff}$, surface gravity log$g$, metallicity [Fe/H], $\alpha$-element to iron (metal) abundance ratio [$\alpha$/Fe] ([$\alpha$/M]), elemental abundances [C/H] and [N/H], and absolute magnitudes ${\rm M}_V$ and ${\rm M}_{K_{\rm s}}$ deduced from 1.8 million spectra of 1.4 million unique stars targeted by the LSS-GAC since September 2011 until June 2014. The catalogues also give values of interstellar reddening, distance and orbital parameters determined with a variety of techniques, as well as proper motions and multi-band photometry from the far-UV to the mid-IR collected from the literature and various surveys. Accuracies of radial velocities reach 5kms$^{-1}$ for late-type stars, and those of distance estimates range between 10 -- 30 per cent, depending on the spectral signal-to-noise ratios. Precisions of [Fe/H], [C/H] and [N/H] estimates reach 0.1dex, and those of [$\alpha$/Fe] and [$\alpha$/M] reach 0.05dex. The large number of stars, the contiguous sky coverage, the simple yet non-trivial target selection function and the robust estimates of stellar radial velocities and atmospheric parameters, distances and elemental abundances, make the catalogues a valuable data set to study the structure and evolution of the Galaxy, especially the solar-neighbourhood and the outer disk.

Silicon Photo-Multipliers (SiPM) are regarded as novel photo-detectors to replace conventional Photo-Multiplier Tubes (PMTs). However, the breakdown voltage dependence on the ambient temperature results in a gain variation of $\sim$3$\% /^{\circ} \mathrm C$. This severely limits the application of this device in experiments with wide range of operating temperature, especially in space missions. An experimental setup was established to investigate the temperature and bias voltage dependence of gain for the Multi-Pixel Photon Counter (MPPC). The gain and breakdown voltage dependence on operating temperature of an MPPC can be approximated by a linear function, which is similar to the behavior of a zener diode. The measured temperature coefficient of the breakdown voltage is $(59.4 \pm 0.4$ mV)$/^{\circ} \mathrm C$. According to this fact, an analog power supply based on two zener diodes and an operational amplifier was designed with a positive temperature coefficient. The measured temperature dependence for the designed power supply is between 63.65 to 64.61~mV/$^{\circ} \mathrm C$ at different output voltages. The designed power supply can bias the MPPC at an over-voltage with a temperature variation of $\sim$ 5~mV$/^{\circ} \mathrm C$. The gain variation of the MPPC biased at over-voltage of 2~V was reduced from 2.8$\% /^{\circ} \mathrm C$ to 0.3$\% /^{\circ} \mathrm C$ when biased the MPPC with the designed power supply for gain control. Detailed design and performance of the analog power supply in the temperature range from -42.7$^{\circ}\mathrm{C}$ to 20.9$^{\circ}\mathrm{C}$ will be discussed in this paper.

We report the observations of an electron vortex magnetic hole corresponding to a new type of coherent structures in the magnetosheath turbulent plasma using the Magnetospheric Multiscale (MMS) mission data. The magnetic hole is characterized by a magnetic depression, a density peak, a total electron temperature increase (with a parallel temperature decrease but a perpendicular temperature increase), and strong currents carried by the electrons. The current has a dip in the center of the magnetic hole and a peak in the outer region of the magnetic hole. The estimated size of the magnetic hole is about 0.23 \r{ho}i (~ 30 \r{ho}e) in the circular cross-section perpendicular to its axis, where \r{ho}i and \r{ho}e are respectively the proton and electron gyroradius. There are no clear enhancement seen in high energy electron fluxes, but an enhancement in the perpendicular electron fluxes at ~ 90{\deg} pitch angles inside the magnetic hole is seen, implying that the electron are trapped within it. The variations of the electron velocity components Vem and Ven suggest that an electron vortex is formed by trapping electrons inside the magnetic hole in the circular cross-section (in the M-N plane). These observations demonstrate the existence of a new type of coherent structures behaving as an electron vortex magnetic hole in turbulent space plasmas as predicted by recent kinetic simulations.

Oxygen and carbon are important elements in stellar populations. Their behavior refers to the formation history of the stellar populations. C and O abundances would also obviously influence stellar opacities and the overall metal abundance $Z$. With observed high-quality spectroscopic properties, we construct stellar models with C and O elements, to give more accurate ages for 70 metal-poor dwarfs, which have been determined to be high-$\alpha$ halo, low-$\alpha$ halo and thick-disk stars. Our results show that high-$\alpha$ halo stars are relatively older than low-$\alpha$ halo stars by around 2.0 Gyr. The thick-disk population has an age range between the two halo populations. The age distribution profiles indicate that high-$\alpha$ halo and low-$\alpha$ halo stars match the in situ accretion simulation by Zolotov et al., and the thick-disk stars might be formed in a relatively quiescent and long-lasting process. We also note that stellar ages are very sensitive to O abundance, since the ages clearly increase with increasing [O/Fe] values. Additionally, we obtain several stars with peculiar ages, including 2 young thick-disk stars and 12 stars older than the universe age.

Full functional and performance tests were performed many times before the Hard X-ray Modulation Telescope (HXMT) launch. During one of the tests, the count rate curves of the 18 High Energy Detectors (HED) have been found increased consistently within an interval of time. A further study on the correlation between the count rate and rainfall was carried out,and the increased net spectrum was also analyzed. The analysis results indicate that the short-lived 222Rn decay products (214Pb and 214Bi) in rainwater were responsible for the transient changes of the background radiation spectra in HEDs. The results show that the HXMT/HEDs have a good detection sensitivity on X/gamma rays, and the detector calibration results are effective.

We report the detection of two new planets orbiting the K giants HD 86950 and HD 222076, based on precise radial velocities obtained with three instruments: AAT/UCLES, FEROS, and CHIRON. HD 86950b has a period of 1270$\pm$57 days at $a=2.72\pm$0.08 AU, and m sin $i=3.6\pm$0.7 Mjup. HD 222076b has $P=871\pm$19 days at $a=1.83\pm$0.03 AU, and m sin $i=1.56\pm$0.11 Mjup. These two giant planets are typical of the population of planets known to orbit evolved stars. In addition, we find a high-amplitude periodic velocity signal ($K\sim$50 m/s) in HD 29399, and show that it is due to stellar variability rather than Keplerian reflex motion. We also investigate the relation between planet occurrence and host-star metallicity for the 164-star Pan-Pacific Planet Search sample of evolved stars. In spite of the small sample of PPPS detections, we confirm the trend of increasing planet occurrence as a function of metallicity found by other studies of planets orbiting evolved stars.

Nuclear fusion reactions are the most important processes in nature to power stars and produce new elements, and lie at the center of the understanding of nucleosynthesis in the universe. It is critically important to study the reactions in full plasma environments that are close to true astrophysical conditions. By using laser-driven counter-streaming collisionless plasmas, we studied the fusion D$+$D$\rightarrow n +^3$He in a Gamow-like window around 27 keV. The results show that astrophysical nuclear reaction yield can be modulated significantly by the self-generated electromagnetic fields and the collective motion of the plasma. This plasma-version mini-collider may provide a novel tool for studies of astrophysics-interested nuclear reactions in plasma with tunable energies in earth-based laboratories.

We use the spectra of Pleiades and field stars from LAMOST DR2 archive to study how spottedness and activity vary as a function of mass at young ages. We obtained standard TiO band strength by measuring TiO bands near 7050 $\AA$ from LAMOST spectra (R$\approx$1800) for large sample of field GKM dwarfs with solar metallicity. Analysis show that active dwarfs, including late G- and early K-type, have extra TiO absorption compare to inactive counterparts, indicating the presence of cool spots on their surface. Active late K- and M-dwarfs show deeper TiO2 and shallower TiO4 compare to inactive stars at a given TiO5, which could be partly explained through cool spots. We estimated cool spot fractional coverage for 304 Pleiades candidates by modelling their TiO2 (&TiO5) band strength with respect to standard value. Results show that surface of large fraction of K- and M-type members have very large spot coverage ($\sim$50%). We analysed a correlation between spot coverage, rotation and the amplitude of light variation, and found spot coverage on slow rotators ($R_{o}$ > 0.1) increases with decreasing Rossby Number $R_{o}$. Interestingly, we detected a saturation-like feature for spot coverage in fast rotators with a saturation level of 40%-50%. In addition, spot distribution in hotter fast rotators show more symmetrical compare to slow rotators. More interestingly, we detected large spot coverage in many M type members with no or little light variation. In bigger picture, these results provide important constraints for stellar dynamo on these cool active stars.

We selected 82 u-band variable objects based on the u-band photometry data from SCUSS and SDSS, in the field of LAMOST Complete Spectroscopic Survey of Pointing Area at Southern Galactic Cap. The magnitude variation of the targets is restricted to larger than 0.2 mag and limiting magnitude down to 19.0 mag in u-band. According to the spectra from LAMOST, there are 11 quasars with red-shift between 0.4 and 1.8, 60 variable stars and 11 unidentified targets. The variable stars include one active M-dwarf with a series of Balmer emission lines, seven Horizontal Branch stars containing six RR Lyrae stars matching with SIMBAD, and one giant, one AGB star and two RR Lyrae candidates by different colour selections. All these variable stars mainly locate near the main sequence in the g-r verse u-g diagram. The quasars are well distinguished from stars by both u-g colour and variation in u-band.

Local helioseismology has opened new frontiers in our quest for understanding of the internal dynamics and dynamo on the Sun. Local helioseismology reconstructs subsurface structures and flows by extracting coherent signals of acoustic waves traveling through the interior and carrying information about subsurface perturbations and flows, from stochastic oscillations observed on the surface. The initial analysis of the subsurface flow maps reconstructed from the 5 years of SDO/HMI data by time-distance helioseismology reveals the great potential for studying and understanding of the dynamics of the quiet Sun and active regions, and the evolution with the solar cycle. In particular, our results show that the emergence and evolution of active regions are accompanied by multi-scale flow patterns, and that the meridional flows display the North-South asymmetry closely correlating with the magnetic activity. The latitudinal variations of the meridional circulation speed, which are probably related to the large-scale converging flows, are mostly confined in shallow subsurface layers. Therefore, these variations do not necessarily affect the magnetic flux transport. The North-South asymmetry is also pronounced in the variations of the differential rotation ("torsional oscillations"). The calculations of a proxy of the subsurface kinetic helicity density show that the helicity does not vary during the solar cycle, and that supergranulation is a likely source of the near-surface helicity.

Local helioseismology provides a unique opportunity to investigate the subsurface structure and dynamics of active regions and their effect on the large-scale flows and global circulation of the Sun. We use measurements of plasma flows in the upper convection zone, provided by the Time-Distance Helioseismology Pipeline developed for analysis of solar oscillation data obtained by Helioseismic and Magnetic Imager (HMI) on Solar Dynamics Observatory (SDO), to investigate the subsurface dynamics of emerging active region NOAA 11726. The active region emergence was detected in deep layers of the convection zone about 12 hours before the first bipolar magnetic structure appeared on the surface, and 2 days before the emergence of most of the magnetic flux. The speed of emergence determined by tracking the flow divergence with depth is about 1.4 km/s, very close to the emergence speed in the deep layers. As the emerging magnetic flux becomes concentrated in sunspots local converging flows are observed beneath the forming sunspots. These flows are most prominent in the depth range 1-3 Mm, and remain converging after the formation process is completed. On the larger scale converging flows around active region appear as a diversion of the zonal shearing flows towards the active region, accompanied by formation of a large-scale vortex structure. This process occurs when a substantial amount of the magnetic flux emerged on the surface, and the converging flow pattern remains stable during the following evolution of the active region. The Carrington synoptic flow maps show that the large-scale subsurface inflows are typical for active regions. In the deeper layers (10-13 Mm) the flows become diverging, and surprisingly strong beneath some active regions. In addition, the synoptic maps reveal a complex evolving pattern of large-scale flows on the scale much larger than supergranulation.

Our aim is to characterize the polarized continuum emission properties including intensity, polarization position angle, and polarization percentage of Sgr A* at $\sim$100 (3.0 mm), $\sim$230 (1.3 mm), $\sim$345 (0.87 mm), $\sim$500 (0.6 mm), and $\sim$700 GHz (0.43 mm). We report continuum emission properties of Sgr A* at the above frequency bands, based on the Atacama Large Millimeter Array (ALMA) observations. We measured flux densities of Sgr A* from ALMA single pointing and mosaic observations. We performed sinusoidal fittings to the observed (XX-YY)/I intensity ratios, to derive the polarization position angles and polarization percentages. We successfully detect polarized continuum emission from all observed frequency bands. We observed lower Stokes I intensity at $\sim$700 GHz than that at $\sim$500 GHz, which suggests that emission at $\gtrsim$500 GHz is from optically thin part of a synchrotron emission spectrum. Both the Stokes I intensity and the polarization position angle at our highest observing frequency of $\sim$700 GHz, may be varying with time. However, we do not yet detect variation in the polarization percentage at $>$500 GHz. The polarization percentage at $\sim$700 GHz is likely lower than that at $\sim$500 GHz. By comparing the $\sim$500 GHz and $\sim$700 GHz observations with the observations at lower frequency bands, we suggest that the intrinsic polarization position angle of Sgr A* is varying with time. This paper also reports the measurable polarization properties from the observed calibration quasars. The future simultaneous multi-frequency polarization observations are required for clarifying the time and frequency variation of polarization position angle and polarization percentage.

Air radon concentration measurement is useful for research on geophysical effects, but it is strongly sensitive to site geology and many geophysical and microclimatic processes such as wind, ventilation, air humidity and so on that induce very big fluctuations on the concentration of radon in air. On the contrary, monitoring the radon concentration in soil by measuring the thermal neutron flux reduces environmental effects. In this paper we report some experimental results on the natural thermal neutron flux as well as the concentration of air radon and its variations at 4300 m a.s.l. These results were obtained with unshielded thermal neutron scintillation detectors (en-detectors) and radon monitors located inside the ARGO-YBJ experimental hall. The correlation of these variations with the lunar month and 1-year period is undoubtedly confirmed. A method for earthquakes prediction provided by a global net of the en-detectors is currently under study.

We report on a measurement of thermal neutrons, generated by the hadronic component of extensive air showers (EAS), by means of a small array of EN-detectors developed for the PRISMA project (PRImary Spectrum Measurement Array), novel devices based on a compound alloy of ZnS(Ag) and $^{6}$LiF. This array has been operated within the ARGO-YBJ experiment at the high altitude Cosmic Ray Observatory in Yangbajing (Tibet, 4300 m a.s.l.). Due to the tight correlation between the air shower hadrons and thermal neutrons, this technique can be envisaged as a simple way to estimate the number of high energy hadrons in EAS. Coincident events generated by primary cosmic rays of energies greater than 100 TeV have been selected and analyzed. The EN-detectors have been used to record simultaneously thermal neutrons and the air shower electromagnetic component. The density distributions of both components and the total number of thermal neutrons have been measured. The correlation of these data with the measurements carried out by ARGO-YBJ confirms the excellent performance of the EN-detector.

We report linearly polarized continuum emission properties of Sgr A* at $\sim$492 GHz, based on the Atacama Large Millimeter Array (ALMA) observations. We used the observations of the likely unpolarized continuum emission of Titan, and the observations of C\textsc{i} line emission, to gauge the degree of spurious polarization. The Stokes I flux of 3.6$\pm$0.72 Jy during our run is consistent with extrapolations from the previous, lower frequency observations. We found that the continuum emission of Sgr A* at $\sim$492 GHz shows large amplitude differences between the XX and the YY correlations. The observed intensity ratio between the XX and YY correlations as a function of parallactic angle may be explained by a constant polarization position angle of $\sim$158$^{\circ}$$\pm$3$^{\circ}$. The fitted polarization percentage of Sgr A* during our observational period is 14\%$\pm$1.2\%. The calibrator quasar J1744-3116 we observed at the same night can be fitted to Stokes I = 252 mJy, with 7.9\%$\pm$0.9\% polarization in position angle P.A. = 4.1$^{\circ}$$\pm$4.2$^{\circ}$. The observed polarization percentage and polarization position angle in the present work appear consistent with those expected from longer wavelength observations in the period of 1999-2005. In particular, the polarization position angle at 492 GHz, expected from the previously fitted 167$^{\circ}$$\pm$7$^{\circ}$ intrinsic polarization position angle and (-5.6$\pm$0.7)$\times$10$^{5}$ rotation measure, is 155$^{+9}_{-8}$, which is consistent with our new measurement of polarization position angle within 1$\sigma$. The polarization percentage and the polarization position angle may be varying over the period of our ALMA 12m Array observations, which demands further investigation with future polarization observations.

Under current LHC and dark matter constraints, the general NMSSM can have self-interacting dark matter to explain the cosmological small structure. In this scenario, the dark matter is the light singlino-like neutralino ($\chi$) which self-interacts through exchanging the light singlet-like scalars ($h_1, a_1$). These light scalars and neutralinos inevitably interact with the 125 GeV SM-like Higgs boson ($h_{SM}$), which cause the Higgs exotic decays $h_{SM} \to h_1 h_1, a_1 a_1, \chi\chi$. We first demonstrate the parameter space required by the explanation of the cosmological small structure and then display the Higgs exotic decays. We find that in such a parameter space the Higgs exotic decays can have branching ratios of a few percent, which should be accessible in the future $e^+e^-$ colliders.

We studied the magnetic topology of active region 12158 on 2014 September 10 and compared it with the observations before and early in the flare which begins at 17:21 UT (SOL2014-09-10T17:45:00). Our results show that the sigmoidal structure and flare ribbons of this active region observed by SDO/AIA can be well reproduced from a Grad-Rubin non linear force free field extrapolation method. Various inverse-S and -J shaped magnetic field lines, that surround a coronal flux rope, coincide with the sigmoid as observed in different extreme ultraviolet wavelengths, including its multi-threaded curved ends. Also, the observed distribution of surface currents in the magnetic polarity where it was not prescribed is well reproduced. This validates our numerical implementation and set-up of the Grad-Rubin method. The modeled double inverse-J shaped Quasi-Separatrix Layer (QSL) footprints match the observed flare ribbons during the rising phase of the flare, including their hooked parts. The spiral-like shape of the latter may be related to a complex pre-eruptive flux rope with more than one turn of twist, as obtained in the model. These ribbon-associated flux-rope QSL-footprints are consistent with the new standard flare model in 3D, with the presence of a hyperbolic flux tube located below an inverse tear drop shaped coronal QSL. This is a new step forward forecasting the locations of reconnection and ribbons in solar flares, and the geometrical properties of eruptive flux ropes.

We present observations of the $^{12}$CO(6-5) line and 686GHz continuum emission in NGC253 with the Submillimeter Array at an angular resolution of ~4arcsec. The $^{12}$CO(6-5) emission is clearly detected along the disk and follows the distribution of the lower $^{12}$CO line transitions with little variations of the line ratios in it. A large-velocity gradient analysis suggests a two-temperature model of the molecular gas in the disk, likely dominated by a combination of low-velocity shocks and the disk wide PDRs. Only marginal $^{12}$CO(6-5) emission is detected in the vicinity of the expanding shells at the eastern and western edges of the disk. While the eastern shell contains gas even warmer (T$_{\rm kin}$>300~K) than the hot gas component (T$_{\rm kin}$=300K) of the disk, the western shell is surrounded by gas much cooler (T$_{\rm kin}$=60K) than the eastern shell but somewhat hotter than the cold gas component of the disk (for similar H$_2$ and CO column densities), indicative of different (or differently efficient) heating mechansisms. The continuum emission at 686GHz in the disk agrees well in shape and size with that at lower (sub-)millimeter frequencies, exhibiting a spectral index consistent with thermal dust emission. We find dust temperatures of ~10-30K and largely optically thin emission. However, our fits suggest a second (more optically thick) dust component at higher temperatures (T$_{\rm d}$>60K), similar to the molecular gas. We estimate a global dust mass of ~10$^6$Msun for the disk translating into a gas-to-dust mass ratio of a few hundred consistent with other nearby active galaxies.

New observations of Sgr A have been carried out with the VLA using the broadband (2 GHz) continuum mode at 5.5 GHz, covering the central 30 pc region of the RBZ at the Galactic center. Using the MS-MFS algorithms in CASA, we have imaged Sgr A with a resolution of 1", achieving an rms 8 $\mu$Jy/beam, and a dynamic range 100,000:1.The radio image is compared with X-ray, CN emission-line and Paschen-$\alpha$ images obtained using Chandra, SMA and HST/NICMOS, respectively. We discuss several prominent radio features. The "Sgr A West Wings" extend 5 pc from the NW and SE tips of the ionized "Mini-spiral" in Sgr A West to positions located 2.9 and 2.4 arc min to the NW and SE of Sgr A*, respectively. The NW wing, along with several other prominent features, including the "NW Streamers", form an elongated radio lobe (NW lobe), oriented nearly perpendicular to the Galactic plane. This radio lobe, with a size of 14.4 pc x 7.3 pc, has a known X-ray counterpart. A row of three thermally emitting rings is observed in the NW lobe. A field containing numerous amorphous radio blobs extends for a distance of ~2 arc min beyond the tip of the SE wing; these features coincide with the SE X-ray lobe. Most of the amorphous radio blobs in the NW and SE lobes have Paschen-$\alpha$ counterparts, suggesting that a shock interaction of ambient gas concentrations with a collimated nuclear wind (outflow) that may be driven by radiation force from the central star cluster within the CND. Finally, we remark on a prominent radio feature located within the shell of the Sgr A East SNR. Because this feature -- the "Sigma Front" -- correlates well in shape and orientation with the nearby edge of the CND, we propose that it is a reflected shock wave resulting from the impact of the Sgr A East blast wave on the CND.

We use archival HARPS spectra to detect three planets orbiting the M3 dwarf Wolf1061 (GJ 628). We detect a 1.36 Mearth minimum-mass planet with an orbital period P = 4.888d (Wolf1061b), a 4.25 Mearth minimum-mass planet with orbital period P = 17.867d (Wolf1061c), and a likely 5.21 Mearth minimum-mass planet with orbital period P = 67.274d (Wolf1061d). All of the planets are of sufficiently low mass that they may be rocky in nature. The 17.867d planet falls within the habitable zone for Wolf 1061 and the 67.274d planet falls just outside the outer boundary of the habitable zone. There are no signs of activity observed in the bisector spans, cross-correlation full-width-half-maxima, Calcium H & K indices, NaD indices, or H-alpha indices near the planetary periods. We use custom methods to generate a cross-correlation template tailored to the star. The resulting velocities do not suffer the strong annual variation observed in the HARPS DRS velocities. This differential technique should deliver better exploitation of the archival HARPS data for the detection of planets at extremely low amplitudes.

We report on the extensive multi-wavelength observations of the blazar Markarian 421 (Mrk 421) covering radio to gamma-rays, during the 4.5 year period of ARGO-YBJ and Fermi common operation time, from August 2008 to February 2013. In particular, thanks to the ARGO-YBJ and Fermi data, the whole energy range from 100 MeV to 10 TeV is covered without any gap. In the observation period, Mrk 421 showed both low and high activity states at all wavebands. The correlations among flux variations in different wavebands were analyzed. Seven large flares, including five X-ray flares and two GeV gamma-ray flares with variable durations (3-58 days), and one X-ray outburst phase were identified and used to investigate the variation of the spectral energy distribution with respect to a relative quiescent phase. During the outburst phase and the seven flaring episodes, the peak energy in X-rays is observed to increase from sub-keV to few keV. The TeV gamma-ray flux increases up to 0.9-7.2 times the flux of the Crab Nebula. The behavior of GeV gamma-rays is found to vary depending on the flare, a feature that leads us to classify flares into three groups according to the GeV flux variation. Finally, the one-zone synchrotron self-Compton model was adopted to describe the emission spectra. Two out of three groups can be satisfactorily described using injected electrons with a power-law spectral index around 2.2, as expected from relativistic diffuse shock acceleration, whereas the remaining group requires a harder injected spectrum. The underlying physical mechanisms responsible for different groups may be related to the acceleration process or to the environment properties.

Alfv\'{e}nic turbulent cascade perpendicular and parallel to the background magnetic field is studied accounting for anisotropic dispersive effects and turbulent intermittency. The perpendicular dispersion and intermittency make the perpendicular-wavenumber magnetic spectra steeper and speed up production of high ion-cyclotron frequencies by the turbulent cascade. On the contrary, the parallel dispersion makes the spectra flatter and decelerate the frequency cascade above the ion-cyclotron frequency. Competition of the above factors results in spectral indices distributed in the interval [-2,-3], where -2 is the index of high-frequency space-filling turbulence, and -3 is the index of low-frequency intermittent turbulence formed by tube-like fluctuations. Spectra of fully intermittent turbulence fill a narrower range of spectral indices [-7/3,-3], which almost coincides with the range of indexes measured in the solar wind. This suggests that the kinetic-scale turbulent spectra are shaped mainly by dispersion and intermittency. A small mismatch with measured indexes of about 0.1 can be associated with damping effects not studied here.

The events recorded by ARGO-YBJ in more than five years of data collection have been analyzed to determine the diffuse gamma-ray emission in the Galactic plane at Galactic longitudes 25{\deg} < l < 100{\deg} and Galactic latitudes . The energy range covered by this analysis, from ~350 GeV to ~2 TeV, allows the connection of the region explored by Fermi with the multi-TeV measurements carried out by Milagro. Our analysis has been focused on two selected regions of the Galactic plane, i.e., 40{\deg} < l < 100{\deg} and 65{\deg} < l < 85{\deg} (the Cygnus region), where Milagro observed an excess with respect to the predictions of current models. Great care has been taken in order to mask the most intense gamma-ray sources, including the TeV counterpart of the Cygnus cocoon recently identified by ARGO-YBJ, and to remove residual contributions. The ARGO-YBJ results do not show any excess at sub-TeV energies corresponding to the excess found by Milagro, and are consistent with the predictions of the Fermi model for the diffuse Galactic emission. From the measured energy distribution we derive spectral indices and the differential flux at 1 TeV of the diffuse gamma-ray emission in the sky regions investigated.

Helioseismic and magnetohydrodynamic waves are abundant in and above sunspots. Through cross-correlating oscillation signals in the photosphere observed by the SDO/HMI, we reconstruct how waves propagate away from virtual wave sources located inside a sunspot. In addition to the usual helioseismic wave, a fast-moving wave is detected traveling along the sunspot's radial direction from the umbra to about 15 Mm beyond the sunspot boundary. The wave has a frequency range of 2.5 - 4.0 mHz with a phase velocity of 45.3 km/s, substantially faster than the typical speeds of Alfven and magnetoacoustic waves in the photosphere. The observed phenomenon is consistent with a scenario of that a magnetoacoustic wave is excited at approximately 5 Mm beneath the sunspot, and its wavefront travels to and sweeps across the photosphere with a speed higher than the local magnetoacoustic speed. The fast-moving wave, if truly excited beneath the sunspot's surface, will help open a new window to study the internal structure and dynamics of sunspots.

We have compared the stellar parameters, temperature, gravity and metallicity, between the LAMOST-DR2 and SDSS-DR12/APOGEE database for stars in common. It is found that the LAMOST database provides a better red-clump feature than the APOGEE database in the Teff versus logg diagram. With this advantage, we have separated red clump stars from red giant stars, and attempt to establish the calibrations between the two datasets for the two groups of stars respectively. It shows that there is a good consistency in temperature with a calibration close to the one-to-one line, and we can establish a satisfied metallicity calibration of [Fe/H]_APO=1.18[Fe/H]_APO+0.11 with a scatter of 0.08 dex for both red clump and red giant branch samples. For gravity, there is no any correlation for red clump stars between the two databases, and scatters around the calibrations of red giant stars are substantial. We found two main sources of the scatters of logg for red giant stars. One is a group of stars with $0.00253*Teff-8.67<logg<2.6 locating at the forbidden region, and the other is the contaminated red clump stars, which could be picked out from the unmatched region where stellar metallicity is not consistent with the position in the Teff versus logg diagram. After excluding stars at the two regions, we have established two calibrations for red giant stars, logg_APO=0.000615*Teff_LAMO+0.697*logg_LAMO-2.208$ (sigma=0.150) for [Fe/H]>-1 and logg_APO=0.000874*Teff_LAMO+0.588*logg_LAMO-3.117 (sigma=0.167) for [Fe/H]<-1. The calibrations are valid for stars with Teff=3800-5400 K and logg=0-3.8 dex, and are useful in a work aiming to combine the LAMOST and APOGEE databases in the future study. In addition, we find that an SVM method based on seismic logg is a good way to greatly improve the accuracy of gravity for these two regions at least in the LAMOST database.

The ARGO-YBJ experiment has been in stable data taking from November 2007 till February 2013 at the YangBaJing Cosmic Ray Observatory (4300 m a.s.l.). The detector consists of a single layer of Resistive Plate Chambers (RPCs) ( about 6700 m^2}) operated in streamer mode. The signal pick-up is obtained by means of strips facing one side of the gas volume. The digital readout of the signals, while allows a high space-time resolution in the shower front reconstruction, limits the measurable energy to a few hundred TeV. In order to fully investigate the 1-10 PeV region, an analog readout has been implemented by instrumenting each RPC with two large size electrodes facing the other side of the gas volume. Since December 2009 the RPC charge readout has been in operation on the entire central carpet (about 5800 m^2). In this configuration the detector is able to measure the particle density at the core position where it ranges from tens to many thousands of particles per m^2. Thus ARGO-YBJ provides a highly detailed image of the charge component at the core of air showers. In this paper we describe the analog readout of RPCs in ARGO-YBJ and discuss both the performance of the system and the physical impact on the EAS measurements.

The ARGO-YBJ experiment is a full-coverage air shower detector located at the Yangbajing Cosmic Ray Observatory (Tibet, People's Republic of China, 4300 m a.s.l.). The high altitude, combined with the full-coverage technique, allows the detection of extensive air showers in a wide energy range and offer the possibility of measuring the cosmic ray proton plus helium spectrum down to the TeV region, where direct balloon/space-borne measurements are available. The detector has been in stable data taking in its full configuration from November 2007 to February 2013. In this paper the measurement of the cosmic ray proton plus helium energy spectrum is presented in the region 3-300 TeV by analyzing the full collected data sample. The resulting spectral index is $\gamma = -2.64 \pm 0.01$. These results demonstrate the possibility of performing an accurate measurement of the spectrum of light elements with a ground based air shower detector.

A local dwarf galaxy, NGC 5253, has a young super star cluster that may provide an example of highly efficient star formation. Here we report the detection and imaging, with the Submillimeter Array, of the J= 3-2 rotational transition of CO at the location of the massive cluster associated with the supernebula. The gas cloud is hot, dense, quiescent, and extremely dusty. Its gas-to-dust ratio is lower than the Galactic value, which we attribute to dust enrichment by Wolf-Rayet stars within the embedded star cluster. Its star formation efficiency exceeds 50%, ten times higher than clouds in the Milky Way: this cloud is a factory of stars and soot. We suggest that high efficiency results from the force-feeding of star formation by a streamer of gas falling into the galaxy.

The search for Gamma Ray Burst (GRB) emission in the energy range 1-100 GeV in coincidence with the satellite detection has been carried out using the Astrophysical Radiation with Ground-based Observatory at YangBaJing (ARGO-YBJ) experiment. The high altitude location (4300 m a.s.l.), the large active surface ($\sim$ 6700 m$^2$ of Resistive Plate Chambers), the wide field of view ($\sim 2~$sr, limited only by the atmospheric absorption) and the high duty cycle ($>$ 86 %) make the ARGO-YBJ experiment particularly suitable to detect short and unexpected events like GRBs. With the scaler mode technique, i.e., counting all the particles hitting the detector with no measurement of the primary energy and arrival direction, the minimum threshold of $\sim$ 1 GeV can be reached, overlapping the direct measurements carried out by satellites. During the experiment lifetime, from December 17, 2004 to February 7, 2013, a total of 206 GRBs occurring within the ARGO-YBJ field of view (zenith angle $\theta$ $\le$ 45$^{\circ}$) have been analyzed. This is the largest sample of GRBs investigated with a ground-based detector. Two lightcurve models have been assumed and since in both cases no significant excess has been found, the corresponding fluence upper limits in the 1-100 GeV energy region have been derived, with values as low as 10$^{-5}~$erg cm$^{-2}$. The analysis of a subset of 24 GRBs with known redshift has been used to constrain the fluence extrapolation to the GeV region together with possible cutoffs under different assumptions on the spectrum.

The measurement of cosmic ray energy spectra, in particular for individual species, is an essential approach in finding their origin. Locating the "knees" of the spectra is an important part of the approach and has yet to be achieved. Here we report a measurement of the mixed Hydrogen and Helium spectrum using the combination of the ARGO-YBJ experiment and of a prototype Cherenkov telescope for the LHAASO experiment. A knee feature at 640+/-87 TeV, with a clear steepening of the spectrum, is observed. This gives fundamental inputs to galactic cosmic ray acceleration models.

Photochemical escape is an important process for oxygen escape from present Mars. In this work, a 1-D Monte-Carlo Model is developed to calculate escape rates of energetic oxygen atoms produced from O2+ dissociative recombination reactions (DR) under 1, 3, 10, and 20 times present solar XUV fluxes. We found that although the overall DR rates increase with solar XUV flux almost linearly, oxygen escape rate increases from 1 to 10 times present solar XUV conditions but decreases when increasing solar XUV flux further. Analysis shows that atomic species in the upper thermosphere of early Mars increases more rapidly than O2+ when increasing XUV fluxes. While the latter is the source of energetic O atoms, the former increases the collision probability and thus decreases the escape probability of energetic O. Our results suggest that photochemical escape be a less important escape mechanism than previously thought for the loss of water and/or CO2 from early Mars.

It is suggested by Chen {\it et al.} that the Earth's surface Ultraviolet irradiance ($280-400$ nm) could influence the Earth's surface temperature variation by "Highly Excited Water Vapor" (HEWV) effect. In this manuscript, we reconstruct the developing history of the ozone layer depth variation from 1860 to 2011 based on the HEWV effect. It is shown that the reconstructed ozone layer depth variation correlates with the observational variation from 1958 to 2005 very well ($R=0.8422$, $P>99.9\%$). From this reconstruction, we may limit the spectra band of the surface Ultraviolet irradiance referred in HEWV effect to Ultraviolet B ($280-320$ nm).

As each solar cycle progresses, remnant magnetic flux from active regions (ARs) migrates poleward to cancel the old-cycle polar field. We describe this polarity reversal process during Cycle 24 using four years (2010.33--2014.33) of line-of-sight magnetic field measurements from the Helioseismic and Magnetic Imager. The total flux associated with ARs reached maximum in the north in 2011, more than two years earlier than the south; the maximum is significantly weaker than Cycle 23. The process of polar field reversal is relatively slow, north-south asymmetric, and episodic. We estimate that the global axial dipole changed sign in October 2013; the northern and southern polar fields (mean above 60$^\circ$ latitude) reversed in November 2012 and March 2014, respectively, about 16 months apart. Notably, the poleward surges of flux in each hemisphere alternated in polarity, giving rise to multiple reversals in the north. We show that the surges of the trailing sunspot polarity tend to correspond to normal mean AR tilt, higher total AR flux, or slower mid-latitude near-surface meridional flow, while exceptions occur during low magnetic activity. In particular, the AR flux and the mid-latitude poleward flow speed exhibit a clear anti-correlation. We discuss how these features can be explained in a surface flux transport process that includes a field-dependent converging flow toward the ARs, a characteristic that may contribute to solar cycle variability.

Whether natural factors could interpret the rise of the Earth's surface temperature is still controversial. Though numerous recent researches have reported apparent correlations between solar activity and the Earth's climate, solar activity has encountered a big problem when describing the rapid global warming after 1970s. Our investigation shows the good positive correlations between the Earth's surface Ultraviolet irradiance (280-400 nm) and the Earth's surface temperature both in temporal and spatial variations by analyzing the global surface Ultraviolet irradiance (280-400 nm) and global surface temperature data from 1980-1999. The rise of CO$_2$ cannot interpret the good positive correlations, and we could even get an opposite result to the good correlations when employing the rise of CO$_2$ to describe the relation between them. Based on the good positive correlations, we suggest a new effect, named "Highly Excited Water Vapor" (HEWV) effect, which can interpret how the Sun influences the Earth's surface temperature reasonably, including the rapid warming after 1970s.

We present dust polarization and CO molecular line images of NGC 7538 IRS1. We combined data from the SMA, CARMA and JCMT telescopes to make images with 2.5 arcsec resolution at 230 and 345 GHz. The images show a remarkable spiral pattern in both the dust polarization and molecular outflow. These data dramatically illustrate the interplay between a high infall rate onto IRS1 and a powerful outflow disrupting the dense, clumpy medium surrounding the star. The images of the dust polarization and the CO outflow presented here provide observational evidence for the exchange of energy and angular momentum between the infall and the outflow. The spiral dust pattern, which rotates through over 180 degrees from IRS1, may be a clumpy filament wound up by conservation of angular momentum in the infalling material. The redshifted CO emission ridge traces the dust spiral closely through the MM dust cores, several of which may contain protostars. We propose that the CO maps the boundary layer where the outflow is ablating gas from the dense gas in the spiral.

Linear properties of kinetic Alfv\'{e}n waves (KAWs) and kinetic slow waves (KSWs) are studied in the framework of two-fluid magnetohydrodynamics. We obtain the wave dispersion relations that are valid in a wide range of the wave frequency $\omega $ and plasma-to-magnetic pressure ratio $\beta $. The KAW frequency can reach and exceed the ion cyclotron frequency at ion kinetic scales, whereas the KSW frequency remains sub-cyclotron. At $\beta \sim 1$, the plasma and magnetic pressure perturbations of both modes are in anti-phase, so that there is nearly no total pressure perturbations. However, these modes exhibit also several opposite properties. At high $% \beta $, the electric polarization ratios of KAWs and KSWs are opposite at the ion gyroradius scale, where KAWs are polarized in sense of electron gyration (right-hand polarized) and KSWs are left-hand polarized. The magnetic helicity $\sigma \sim 1$ for KAWs and $\sigma \sim -1$ for KSWs, and the ion Alfv\'{e}n ratio $R_{Ai}\ll 1$ for KAWs and $R_{Ai}\gg 1$ for KSWs. We also found transition wavenumbers where KAWs change their polarization from left- to right-hand. These new properties can be used to discriminate KAWs and KSWs when interpreting kinetic-scale electromagnetic fluctuations observed in various solar-terrestrial plasmas. This concerns, in particular, identification of modes responsible for kinetic-scale pressure-balanced fluctuations and turbulence in the solar wind.

The extended TeV gamma-ray source ARGO J2031+4157 (or MGRO J2031+41) is positionally consistent with the Cygnus Cocoon discovered by $Fermi$-LAT at GeV energies in the Cygnus superbubble. Reanalyzing the ARGO-YBJ data collected from November 2007 to January 2013, the angular extension and energy spectrum of ARGO J2031+4157 are evaluated. After subtracting the contribution of the overlapping TeV sources, the ARGO-YBJ excess map is fitted with a two-dimensional Gaussian function in a square region of $10^{\circ}\times 10^{\circ}$, finding a source extension $\sigma_{ext}$= 1$^{\circ}$.8$\pm$0$^{\circ}$.5. The observed differential energy spectrum is $dN/dE =(2.5\pm0.4) \times 10^{-11}(E/1 TeV)^{-2.6\pm0.3}$ photons cm$^{-2}$ s$^{-1}$ TeV$^{-1}$, in the energy range 0.2-10 TeV. The angular extension is consistent with that of the Cygnus Cocoon as measured by $Fermi$-LAT, and the spectrum also shows a good connection with the one measured in the 1-100 GeV energy range. These features suggest to identify ARGO J2031+4157 as the counterpart of the Cygnus Cocoon at TeV energies. The Cygnus Cocoon, located in the star-forming region of Cygnus X, is interpreted as a cocoon of freshly accelerated cosmic rays related to the Cygnus superbubble. The spectral similarity with Supernova Remnants indicates that the particle acceleration inside a superbubble is similar to that in a SNR. The spectral measurements from 1 GeV to 10 TeV allows for the first time to determine the possible spectrum slope of the underlying particle distribution. A hadronic model is adopted to explain the spectral energy distribution.

Solar subsurface zonal and meridional-flow profiles during the rising phase of solar cycle 24 are studied using time-distance helioseismology technique. The faster zonal bands in the torsional-oscillation pattern show strong hemispheric asymmetries and temporal variations in both width and speed. The faster band in the northern hemisphere is located closer to the equator than the band in the southern hemisphere, and migrates past the equator when the magnetic activity in the southern hemisphere is reaching maximum. The meridional-flow speed decreases substantially with the increase of magnetic activity, and the flow profile shows two zonal structures in each hemisphere. The residual meridional flow, after subtracting a mean meridional-flow profile, converges toward the activity belts and shows faster and slower bands like the torsional-oscillation pattern. More interestingly, the meridional-flow speed above latitude 30 degree shows an anti-correlation with the poleward-transporting magnetic flux, slower when the following-polarity flux is transported and faster when the leading-polarity flux is transported. It is expected that this phenomenon slows the process of magnetic cancellation and polarity reversal in the high-latitude areas.

We analyze the kinematics of thick disk and halo stars observed by the Large sky Area Multi-Object Fiber Spectroscopic Telescope. We have constructed a sample of 7,993 F, G and K nearby main-sequence stars (\textit{d} $<$ 2 kpc) with estimates of position (x, y, z) and space velocity ($U$, $V$, $W$) based on color and proper motion from the SDSS DR9 catalog. Three `phase-space overdensities' are identified in [\textit{V}, $\sqrt{U^{2}+2V^{2}}$] with significance levels of $\sigma$ $>$ 3. %[L$_{Z}$, eccentricity], [L$_{Z}$, L$_{\bot}$], and [V$_{az}$, V$_{\triangle}E$]. Two of them (Hyades-Pleiades stream, Arcturus-AF06 stream) have been identified previously. We also find evidence for a new stream (centered at \textit{V} $\sim$ -180 km s$^{-1}$) in the halo. The formation mechanisms of these three streams are analyzed. Our results support the hypothesis the Arcturus-AF06 stream and the new stream originated from the debris of a disrupted satellite, while Hyades-Pleiades stream has a dynamical origin.

We studied the temporal evolution of the magnetic topology of the active region (AR) 11158 based on the reconstructed three-dimensional magnetic fields in the corona. The \nlfff\ extrapolation method was applied to the 12 minutes cadence data obtained with the \hmi\ (HMI) onboard the \sdo\ (SDO) during five days. By calculating the squashing degree factor Q in the volume, the derived quasi-separatrix layers (QSLs) show that this AR has an overall topology, resulting from a magnetic quadrupole, including an hyperbolic flux tube (HFT) configuration which is relatively stable at the time scale of the flare ($\sim 1-2$ hours). A strong QSL, which corresponds to some highly sheared arcades that might be related to the formation of a flux rope, is prominent just before the M6.6 and X2.2 flares, respectively. These facts indicate the close relationship between the strong QSL and the high flare productivity of AR 11158. In addition, with a close inspection of the topology, we found a small-scale HFT which has an inverse tear-drop structure above the aforementioned QSL before the X2.2 flare. It indicates the existence of magnetic flux rope at this place. Even though a global configuration (HFT) is recognized in this AR, it turns out that the large-scale HFT only plays a secondary role during the eruption. In final, we dismiss a trigger based on the breakout model and highlight the central role of the flux rope in the related eruption.

Helioseismology provides important constraints for the solar dynamo problem. However, the basic properties and even the depth of the dynamo process, which operates also in other stars, are unknown. Most of the dynamo models suggest that the toroidal magnetic field that emerges on the surface and forms sunspots is generated near the bottom of the convection zone, in the tachocline. However, there is a number of theoretical and observational problems with justifying the deep-seated dynamo models. This leads to the idea that the subsurface angular velocity shear may play an important role in the solar dynamo. Using helioseismology measurements of the internal rotation and meridional circulation, we investigate a mean-field MHD model of dynamo distributed in the bulk of the convection zone but shaped in a near-surface layer. We show that if the boundary conditions at the top of the dynamo region allow the large-scale toroidal magnetic fields to penetrate into the surface, then the dynamo wave propagates along the isosurface of angular velocity in the subsurface shear layer, forming the butterfly diagram in agreement with the Parker-Yoshimura rule and solar-cycle observations. Unlike the flux-transport dynamo models, this model does not depend on the transport of magnetic field by meridional circulation at the bottom of the convection zone, and works well when the meridional circulation forms two cells in radius, as recently indicated by deep-focus time-distance helioseismology analysis of the SDO/HMI and SOHO/MDI data. We compare the new dynamo model with various characteristics if the solar magnetic cycles, including the cycle asymmetry (Waldmeier's relations) and magnetic `butterfly' diagrams.

The energy spectrum of cosmic Hydrogen and Helium nuclei has been measured, below the so-called "knee", by using a hybrid experiment with a wide field-of-view Cherenkov telescope and the Resistive Plate Chamber (RPC) array of the ARGO-YBJ experiment at 4300 m above sea level. The Hydrogen and Helium nuclei have been well separated from other cosmic ray components by using a multi-parameter technique. A highly uniform energy resolution of about 25% is achieved throughout the whole energy range (100 TeV - 700 TeV). The observed energy spectrum is compatible with a single power law with index gamma=-2.63+/-0.06.

At a distance of 8 kpc, the center of our Galaxy is the nearest galactic nucleus, and has been the subject of numerous key projects undertaken by great observatories such as Chandra, Spitzer, and Herschel. However, there are still no surveys of molecular gas properties in the Galactic center with less than 30" (1 pc) resolution. There is also no sensitive polarization survey of this region, despite numerous nonthermal magnetic features apparently unique to the central 300 parsecs. In this paper, we outline the potential the VLASS has to fill this gap. We assess multiple considerations in observing the Galactic center, and recommend a C-band survey with 10 micro-Jy continuum RMS and sensitive to molecular gas with densities greater than 10^4 cm^{-3}, covering 17 square degrees in both DnC and CnB configurations ( resolution ~5"), totaling 750 hours of observing time. Ultimately, we wish to note that the upgraded VLA is not just optimized for fast continuum surveys, but has a powerful correlator capable of simultaneously observing continuum emission and dozens of molecular and recombination lines. This is an enormous strength that should be fully exploited and highlighted by the VLASS, and which is ideally suited for surveying the center of our Galaxy.

Deep imaging of the Sgr A complex at 6 cm wavelength with the B and C configurations of the Karl G. Jansky VLA has revealed a new population of faint radio filaments. Like their brighter counterparts that have been observed throughout the Galactic Center on larger scales, these filaments can extend up to ~10 parsecs, and in most cases are strikingly uniform in brightness and curvature. Comparison with a survey of Paschen-alpha emission reveals that some of the filaments are emitting thermally, but most of these structures are nonthermal: local magnetic flux tubes illuminated by synchrotron emission. The new image reveals considerable filamentary substructure in previously known nonthermal filaments (NTFs). Unlike NTFs previously observed on larger scales, which tend to show a predominant orientation roughly perpendicular to the Galactic plane, the NTFs in the vicinity of the Sgr A complex are relatively randomly oriented. Two well-known radio sources to the south of Sgr A - sources E and F - consist of numerous quasi-parallel filaments that now appear to be particularly bright portions of a much larger, strongly curved, continuous, nonthermal radio structure that we refer to as the "Southern Curl". It is therefore unlikely that sources E and F are HII regions or pulsar wind nebulae. The Southern Curl has a smaller counterpart on the opposite side of the Galactic Center - the Northern Curl - that, except for its smaller scale and smaller distance from the center, is roughly point-reflection symmetric with respect to the Southern Curl. The curl features indicate that some field lines are strongly distorted, presumably by mass flows. The point symmetry about the center then suggests that the flows originate near the center and are somewhat collimated.

Based on our deep image of Sgr A using broadband data observed with the Jansky VLA at 6 cm, we present a new perspective of the radio bright zone at the Galactic center. We further show the radio detection of the X-ray Cannonball, a candidate neutron star associated with the Galactic center SNR Sgr A East. The radio image is compared with the Chandra X-ray image to show the detailed structure of the radio counterparts of the bipolar X-ray lobes. The bipolar lobes are likely produced by the winds from the activities within Sgr A West, which could be collimated by the inertia of gas in the CND, or by the momentum driving of Sgr A*; and the poloidal magnetic fields likely play an important role in the collimation. The less-collimated SE lobe, in comparison to the NW one, is perhaps due to the fact that the Sgr A East SN might have locally reconfigured the magnetic field toward negative galactic latitudes. In agreement with the X-ray observations, the time-scale of ~ $1\times10^4$ yr estimated for the outermost radio ring appears to be comparable to the inferred age of the Sgr A East SNR.

The ARGO-YBJ detector is an extensive air shower array that has been used to monitor the northern $\gamma$-ray sky at energies above 0.3 TeV from 2007 November to 2013 January. In this paper, we present the results of a sky survey in the declination band from $-10^{\circ}$ to $70^{\circ}$, using data recorded over the past five years. With an integrated sensitivity ranging from 0.24 to $\sim$1 Crab units depending on the declination, six sources have been detected with a statistical significance greater than 5 standard deviations. Several excesses are also reported as potential $\gamma$-ray emitters. The features of each source are presented and discussed. Additionally, $95\%$ confidence level upper limits of the flux from the investigated sky region are shown. Specific upper limits for 663 GeV $\gamma$-ray AGNs inside the ARGO-YBJ field of view are reported. The effect of the absorption of $\gamma$-rays due to the interaction with extragalactic background light is estimated.

To study the build up of a magnetic flux rope before a major flare and coronal mass ejection (CME), we compute the magnetic helicity injection, twist accumulation, and the topology structure of the three dimensional magnetic field, which is derived by the nonlinear force-free field model. The Extreme-ultraviolet Imaging Telescope on board the Solar and Heliospheric Observatory observed a series of confined flares without any CME before a major flare with a CME at 23:02 UT on 2005 January 15 in active region NOAA 10720. We derive the vector velocity at eight time points from 18:27 UT to 22:20 UT with the differential affine velocity estimator for vector magnetic fields, which were observed by the Digital Vector Magnetograph at Big Bear Solar Observatory. The injected magnetic helicity is computed with the vector magnetic and velocity fields. The helicity injection rate was (-16.47 \pm 3.52) \times 10^{40} Mx^2/hr. We find that only about 1.8% of the injected magnetic helicity became finally the internal helicity of the magnetic flux rope, whose twist increasing rate was -0.18 \pm 0.08 Turns/hr. The quasi-separatrix layers (QSLs) of the three dimensional magnetic field are computed by evaluating the squashing degree, Q. We find that the flux rope was wrapped by QSLs with large Q values, where the magnetic reconnection induced by the continuously injected magnetic helicity further produced the confined flares. We suggest that the flux rope was built up and heated by the magnetic reconnection in the QSLs.

We report new results based on the analysis of the SMA and CARMA observations of NGC 7538\,IRS\,1 at 1.3 and 3.4 mm with sub-arcsec resolutions. With angular resolutions $\sim$ 0\farcs7, the SMA and CARMA observations show that the continuum emission at 1.3 and 3.4 mm from the hyper-compact \ion{H}{2} region IRS\,1 is dominated by a compact source with a tail-like extended structure to the southwest of IRS\,1. With a CARMA B-array image at 1.3 mm convolved to 0\farcs1, we resolve the hyper-compact \ion{H}{2} region into two components: an unresolved hyper-compact core, and a north-south extension with linear sizes of $<270$ AU and $\sim$2000 AU, respectively. The fine structure observed with CARMA is in good agreement with the previous VLA results at centimeter wavelengths, suggesting that the hyper-compact \ion{H}{2} region at the center of IRS\,1 is associated with an ionized bipolar outflow. We image the molecular lines OCS(19-18) and CH$_3$CN(12-11) as well as $^{13}$CO(2-1) surrounding IRS\,1, showing a velocity gradient along the southwest-northeast direction. The spectral line profiles in $^{13}$CO(2-1), CO(2-1), and HCN(1-0) observed toward IRS\,1 show broad redshifted absorption, providing evidence for gas infall with rates in the range of $3-10\times10^{-3}$ M$_\odot$ yr$^{-1}$ inferred from our observations.

We report the VLA detection of the radio counterpart of the X-ray object referred to as the "Cannonball", which has been proposed to be the remnant neutron star resulting from the creation of the Galactic Center supernova remnant, Sagittarius A East. The radio object was detected both in our new VLA image from observations in 2012 at 5.5 GHz and in archival VLA images from observations in 1987 at 4.75 GHz and in the period from 1990 to 2002 at 8.31 GHz. The radio morphology of this object is characterized as a compact, partially resolved point source located at the northern tip of a radio "tongue" similar to the X-ray structure observed by Chandra. Behind the Cannonball, a radio counterpart to the X-ray plume is observed. This object consists of a broad radio plume with a size of 30\arcsec$\times$15\arcsec, followed by a linear tail having a length of 30\arcsec. The compact head and broad plume sources appear to have relatively flat spectra ($\propto\nu^\alpha$) with mean values of $\alpha=-0.44\pm0.08$ and $-0.10\pm0.02$, respectively; and the linear tail shows a steep spectrum with the mean value of $-1.94\pm0.05$. The total radio luminosity integrated from these components is $\sim8\times10^{33}$ erg s$^{-1}$, while the emission from the head and tongue amounts for only $\sim1.5\times10^{31}$ erg s$^{-1}$. Based on the images obtained from the two epochs' observations at 5 GHz, we infer the proper motion of the object: $\mu_\alpha = 0.001 \pm0.003$ arcsec yr$^{-1}$ and $\mu_\delta = 0.013 \pm0.003$ arcsec yr$^{-1}$. With an implied velocity of 500 km s$^{-1}$, a plausible model can be constructed in which a runaway neutron star surrounded by a pulsar wind nebula was created in the event that produced Sgr A East. The inferred age of this object, assuming that its origin coincides with the center of Sgr A East, is approximately 9000 years.

Some experiments have found multi-core events in cosmic rays extensive air showers which should be interpreted by hadronic interaction theory. In this paper, the multi-core events are reproduced by Monte Carlo simulation with CORSIKA. The origin of each sub-cores is tracked back from the observation level. The interaction mechanism and original particles of sub-core are studied in this paper.

We present a set of white dwarf-main sequence (WDMS) binaries identified spectroscopically from the Large sky Area Multi-Object fiber Spectroscopic Telescope (LAMOST, also called the Guo Shou Jing Telescope) pilot survey. We develop a color selection criteria based on what is so far the largest and most complete Sloan Digital Sky Survey (SDSS) DR7 WDMS binary catalog and identify 28 WDMS binaries within the LAMOST pilot survey. The primaries in our binary sample are mostly DA white dwarfs except for one DB white dwarf. We derive the stellar atmospheric parameters, masses, and radii for the two components of 10 of our binaries. We also provide cooling ages for the white dwarf primaries as well as the spectral types for the companion stars of these 10 WDMS binaries. These binaries tend to contain hot white dwarfs and early-type companions. Through cross-identification, we note that nine binaries in our sample have been published in the SDSS DR7 WDMS binary catalog. Nineteen spectroscopic WDMS binaries identified by the LAMOST pilot survey are new. Using the 3$\sigma$ radial velocity variation as a criterion, we find two post-common-envelope binary candidates from our WDMS binary sample.

Meridional flow in the solar interior plays an important role in redistributing angular momentum and transporting magnetic flux inside the Sun. Although it has long been recognized that the meridional flow is predominantly poleward at the Sun's surface and in its shallow interior, the location of the equatorward return flow and the meridional flow profile in the deeper interior remain unclear. Using the first two years of continuous helioseismology observations from the Solar Dynamics Observatory / Helioseismic Magnetic Imager, we analyze travel times of acoustic waves that propagate through different depths of the solar interior carrying information about the solar interior dynamics. After removing a systematic center-to-limb effect in the helioseismic measurements and performing inversions for flow speed, we find that the poleward meridional flow of a speed of 15 m/s extends in depth from the photosphere to about 0.91 R_sun. An equatorward flow of a speed of 10 m/s is found between 0.82 to 0.91 R_sun in the middle of the convection zone. Our analysis also shows evidence of that the meridional flow turns poleward again below 0.82 R_sun, indicating an existence of a second meridional circulation cell below the shallower one. This double-cell meridional circulation profile with an equatorward flow shallower than previously thought suggests a rethinking of how magnetic field is generated and redistributed inside the Sun.

For the M active catalogue of Guoshoujing Telescope (LAMOST), 933 sources are presented in at least two exposures. We found that many M active stars show chromospheric variabilities in the Ca H, H$_{\alpha}$, H$_{\beta}$, and H$_{\gamma}$ lines on short or long timescales.

We present a spectroscopic catalog of 67082 M dwarfs from the LAMOST pilot survey. For each spectrum of the catalog, spectral subtype, radial velocity, equivalent width of H${\alpha}$, a number of prominent molecular band indices and the metal sensitive parameter $\zeta$ are provided . Spectral subtype have been estimated by a remedied Hammer program (Original Hammer: Covey et al. 2007), in which indices are reselected to obtain more accurate auto-classified spectral subtypes. All spectra in this catalog have been visually inspected to confirm the spectral subtypes. Radial velocities have been well measured by our developed program which uses cross-correlation method and estimates uncertainty of radial velocity as well. We also examine the magnetic activity properties of M dwarfs traced by H${\alpha}$ emission line. The molecular band indices included in this catalog are temperature or metallicity sensitive and can be used for future analysis of the physical properties of M dwarfs. The catalog is available on the website \url{http://sciwiki.lamost.org/MCatalogPilot/}.

We present a spectroscopically identified catalogue of 72 DA white dwarfs from the LAMOST pilot survey. 35 are found to be new identifications after cross-correlation with the Eisenstein et al. and Villanova catalogues. The effective temperature and gravity of these white dwarfs are estimated by Balmer lines fitting. Most of them are hot white dwarfs. The cooling times and masses of these white dwarfs are estimated by interpolation in theoretical evolution tracks. The peak of mass distribution is found to be $\sim$ 0.6 $M_\odot$ which is consistent with prior work in the literature. The distances of these white dwarfs are estimated using the method of Synthetic Spectral Distances. All of these WDs are found to be in the Galactic disk from our analysis of space motions. Our sample supports the expectation white dwarfs with high mass are concentrated near the plane of Galactic disk

We present chromospheric activity index $S\rm_{HK}$ measurements for over 13,000 F, G and K disk stars with high signal-to-noise ratio ($>$ 60) spectra in the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) spectroscopic sample. A parameter $\delta$S is defined as the difference between $S\rm_{HK}$ and a `zero' emission line fitted by several of the most inactive stars. The $S\rm_{HK}$ indices of subgiant stars tend to be much lower than dwarfs, which provide a way to distinguish dwarfs and giants with relatively low resolution spectra. Cooler stars are generally more active and display a larger scatter than hotter stars. Stars associated with the thick disk are in general less active than those of the thin disk. The fraction of K dwarfs that are active drops with vertical distance from the Galactic plane. Metallicity affects $S\rm_{HK}$ measurements differently among F, G and K dwarfs in this sample. Using the open clusters NGC 2420, M67 and NGC6791 as calibrations, ages of most field stars in this SDSS sample range from 3-8 Gyr.

We have measured the CH G band (CH(G)) index for evolved stars in the globular cluster M3 based on the Sloan Digital Sky Survey (SDSS) spectroscopic survey. It is found that there is a useful way to select red giant branch (RGB) stars from the contamination of other evolved stars such as asymptotic giant branch (AGB) and red horizontal branch (RHB) stars by using the CH(G) index versus (g-r) diagram if the metallicity is known from the spectra. When this diagram is applied to field giant stars with similar metallicity, we establish a calibration of CH(G)=1.625(g-r)-1.174(g-r)^2-0.934. This method is confirmed by stars with [Fe/H]~2.3 where spectra of member stars in globular clusters M15 and M92 are available in the SDSS database. We thus extend this kind of calibration to every individual metallicity bin ranging from [Fe/H] ~ -3.0 to [Fe/H] ~ 0.0 by using field red giant stars with 0.4 < (g-r) < 1.0. The metallicity-dependent calibrations give CH(G) =1.625(g-r)-1.174(g-r)^2+ 0.060[Fe/H]-0.830 for -3.0 < [Fe/H] <-1.2 and CH(G)=0.953(g-r)-0.655(g-r)^2+0.060[Fe/H]-0.650 for -1.2 < [Fe/H] < 0.0. The calibrations are valid for the SDSS spectroscopic data set, and they cannot be applied blindly to other data sets. With the two calibrations, a significant number of the contaminating stars (AGB and RHB stars) were excluded and thus a clear sample of red giant stars is obtained by selecting stars within 0.05 mag of the calibration. The sample is published online and it is expected that this large and clean sample of RGB stars will provide new information on the formation and evolution of the Galaxy.

We report the observation of a very high energy \gamma-ray source, whose position is coincident with HESS J1841-055. This source has been observed for 4.5 years by the ARGO-YBJ experiment from November 2007 to July 2012. Its emission is detected with a statistical significance of 5.3 standard deviations. Parameterizing the source shape with a two-dimensional Gaussian function we estimate an extension \sigma=(0.40(+0.32,-0.22}) degree, consistent with the HESS measurement. The observed energy spectrum is dN/dE =(9.0-+1.6) x 10^{-13}(E/5 TeV)^{-2.32-+0.23} photons cm^{-2} s^{-1} TeV^{-1}, in the energy range 0.9-50 TeV. The integral \gamma-ray flux above 1 TeV is 1.3-+0.4 Crab units, which is 3.2-+1.0 times the flux derived by HESS. The differences in the flux determination between HESS and ARGO-YBJ, and possible counterparts at other wavelengths are discussed.

An investigation on correlations between photospheric current helicity and subsur- face kinetic helicity is carried out by analyzing vector magnetograms and subsurface velocities for two rapidly developing active regions. The vector magnetograms are from the SDO/HMI (Solar Dynamics Observatory / Helioseismic and Magnetic Im- ager) observed Stokes parameters, and the subsurface velocity is from time-distance data-analysis pipeline using HMI Dopplergrams. Over a span of several days, the evo- lution of the weighted current helicity shows a tendency similar to that of the weighted subsurface kinetic helicity, attaining a correlation coefficient above 0.60 for both ac- tive regions. Additionally, there seems to be a phase lag between the evolutions of the unweighted current and subsurface kinetic helicities for one of the active regions. The good correlation between these two helicities indicate that there is some intrinsic con- nection between the interior dynamics and photospheric magnetic twistedness inside active regions, which may help to interpret the well-known hemispheric preponder- ance of current-helicity distribution.

The continuous high spatial-resolution Doppler observation of the Sun by Solar Dynamics Observatory / Helioseismic and Magnetic Imager allows us to compute helioseismic k-omega power-spectrum diagram using only oscillations inside a sunspot. Individual modal ridges can be clearly seen with reduced power in the k-omega diagram constructed by use of 40-hour observation of a stable and round sunspot. Comparing with the k-omega diagram obtained from a quiet-Sun region, inside the sunspot the f-mode ridge gets more power reduction than p-mode ridges, especially at high wavenumber. The p-mode ridges all shift toward lower-wavenumber (or higher-frequency) areas for a given frequency (or wavenumber), implying an increase of phase velocity beneath the sunspot. This probably results from acoustic waves' travel across the inclined magnetic field of the sunspot penumbra. Line-profile asymmetries exhibited in the p-mode ridges are more significant in the sunspot than in quiet Sun. Convection inside the sunspot is also highly suppressed, and its characteristic spatial scale is substantially larger than the typical convection scale of quiet Sun. These observational facts demand a better understanding of magnetoconvection and interactions of helioseismic waves with magnetic field.

We performed 3D numerical simulations of the solar surface wave field for the quiet Sun and for three models with different localized sound-speed variations in the interior with: (i) deep, (ii) shallow, and (iii) two-layer structures. We used simulated data generated by two different codes which use the same standard solar model as a background model, but utilize two different integration techniques and use different models of stochastic wave excitation. Acoustic travel times were measured from all data sets using the time-distance helioseismology technique and compared with the ray theory predictions, frequently used for helioseismic travel-time inversions. It is found that the measured travel-time shifts agree well with the ray theory in both cases with and without phase-speed filtering for the shallow and deep perturbations. This testing verifies the whole measuring-filtering-inversion procedure for sound-speed anomalies inside the Sun. It is shown, that the phase-speed filtering, frequently used to improve the signal-to-noise ratio does not introduce significant systematic errors. Results of the sound-speed inversion procedure show good agreement with the background sound-speed profiles in all cases. Due to its smoothing nature, the inversion procedure overestimates sound speed variations in areas with sharp gradients of the sound-speed profile.

The meridional flow in the Sun is an axisymmetric flow that is generally poleward directed at the surface, and is presumed to be of fundamental importance in the generation and transport of magnetic fields. Its true shape and strength, however, is debated. We present a numerical simulation of helioseismic wave propagation in the whole solar interior in the presence of a prescribed, stationary, single-cell, deep meridional circulation serving as a test-bed for helioseismic measurement techniques. A deep-focusing time-distance helioseismology technique is applied to the artificial data showing that it can in fact be used to measure the effects of the meridional flow very deep in the solar convection zone. It is shown that the ray-approximation which is commonly used for interpretation of helioseismology measurements remains a reasonable approximation even for the very long distances between 12 and 42 degrees corresponding to depths between 52 and 195 Mm considered here. From the measurement noise we extrapolate that on the order of a full solar cycle may be needed to probe the flow all the way to the base of the convection zone.

As one of the brightest active blazars in both X-ray and very high energy $\gamma$-ray bands, Mrk 501 is very useful for physics associated with jets from AGNs. The ARGO-YBJ experiment is monitoring it for $\gamma$-rays above 0.3 TeV since November 2007. Starting from October 2011 the largest flare since 2005 is observed, which lasts to about April 2012. In this paper, a detailed analysis is reported. During the brightest $\gamma$-rays flaring episodes from October 17 to November 22, 2011, an excess of the event rate over 6 $\sigma$ is detected by ARGO-YBJ in the direction of Mrk 501, corresponding to an increase of the $\gamma$-ray flux above 1 TeV by a factor of 6.6$\pm$2.2 from its steady emission. In particular, the $\gamma$-ray flux above 8 TeV is detected with a significance better than 4 $\sigma$. Based on time-dependent synchrotron self-Compton (SSC) processes, the broad-band energy spectrum is interpreted as the emission from an electron energy distribution parameterized with a single power-law function with an exponential cutoff at its high energy end. The average spectral energy distribution for the steady emission is well described by this simple one-zone SSC model. However, the detection of $\gamma$-rays above 8 TeV during the flare challenges this model due to the hardness of the spectra. Correlations between X-rays and $\gamma$-rays are also investigated.

We statistically study the property of emerging flux regions (EFRs) and the upper solar atmosphere response to the flux emergence by using data from the Helioseismic and Magnetic Imager (HMI) and the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO). Parameters including the total emerged flux, the flux growth rate, the maximum area, the duration of the emergence and the separation speed of the opposite polarities are adopted to delineate the property of the EFRs. The response of the upper atmosphere is addressed by the response of the atmosphere at different wavelengths (and thus at different temperatures). According to our results, the total emerged fluxes are in the range of (0.44 -- 11.2)$\times10^{19}$ Mx while the maximum area ranges from 17 to 182 arcsec$^2$. The durations of the emergence are between 1 and 12 hours, which are positively correlated to both the total emerged flux and the maximum area. The maximum distances between the opposite polarities are 7 -- 25 arcsec and are also correlated to the duration positively. The separation speeds are from 0.05 to 1.08 km s$^{-1}$, negatively correlated to the duration. The derived flux growth rates are (0.1 -- 1.3)$\times10^{19}$ Mx hr$^{-1}$, which are positively correlated to the total emerging flux. The upper atmosphere responds to the flux emergence in the 1600\AA\ chromospheric line first, and then tens and hundreds of seconds later, in coronal lines, such as the 171\AA\ (T=10$^{5.8}$ K) and 211\AA\ (T=10$^{6.3}$ K) lines almost simultaneously, suggesting the successively heating of atmosphere from the chromosphere to the corona.

The results of Monte-Carlo simulations of Extensive Air Shower are presented to show the difference of hadronic component content at various altitudes with the aim to choose an optimal altitude for the PRISMA-like experiment. CORSIKA program for EAS simulations with QGSJET and GHEISHA models was used to calculate the number of hadrons reaching the observational level inside a ring of 50 m radius around the EAS axis. Then the number of neutrons produced by the hadronic component was calculated using an empirical relationship between the two components. We have tested the results with the ProtoPRISMA array at sea level, and recorded neutrons are close to the simulation results.

Investigating the properties of magnetic flux emergence is one of the most important problems of solar physics. In this study we present a newly developed deep-focus time-distance measurement scheme which is able to detect strong emerging flux events in the deep solar interior, before the flux becomes visible on the surface. We discuss in detail the differences between our method and previous methods, and demonstrate step-by-step how the signal-to-noise (S/N) ratio is increased. The method is based on detection of perturbations in acoustic phase travel times determined from cross-covariances of solar oscillations observed on the surface. We detect strong acoustic travel-time reductions of an order of 12 - 16 seconds at a depth of 42 - 75 Mm. These acoustic anomalies are detected 1 - 2 days before high peaks in the photospheric magnetic flux rate implying that the average emerging speed is 0.3 - 0.6 km/s. The results of this work contribute to our understanding of solar magnetism and benefit space weather forecasting.

We report on a systematic center-to-limb variation in measured helioseismic travel times, which must be taken into account for an accurate determination of solar interior meridional flows. The systematic variation, found in time-distance helioseismology analysis using SDO/HMI and SDO/AIA observations, is different in both travel-time magnitude and variation trend for different observables. It is not clear what causes this systematic effect. Subtracting the longitude-dependent east-west travel times, obtained along the equatorial area, from the latitude-dependent north-south travel times, obtained along the central meridian area, gives remarkably similar results for different observables. We suggest this as an effective procedure for removing the systematic center-to-limb variation. The subsurface meridional flows obtained from inversion of the corrected travel times are approximately 10 m/s slower than those obtained without removing the systematic effect. The detected center-to-limb variation may have important implications in the derivation of meridional flows in the deep interior, and needs a better understanding.

We report the observation of TeV gamma-rays from the Cygnus region using the ARGO-YBJ data collected from 2007 November to 2011 August. Several TeV sources are located in this region including the two bright extended MGRO J2019+37 and MGRO J2031+41. According to the Milagro data set, at 20 TeV MGRO J2019+37 is the most significant source apart from the Crab Nebula. No signal from MGRO J2019+37 is detected by the ARGO-YBJ experiment, and the derived flux upper limits at 90% confidence level for all the events above 600 GeV with medium energy of 3 TeV are lower than the Milagro flux, implying that the source might be variable and hard to be identified as a pulsar wind nebula. The only statistically significant (6.4 standard deviations) gamma-ray signal is found from MGRO J2031+41, with a flux consistent with the measurement by Milagro.

We present high angular resolution continuum observations of the high-mass protostar NGC 7538S with BIMA and CARMA at 3 and 1.4 mm, VLA observations at 1.3, 2, 3.5 and 6 cm, and archive IRAC observations from the Spitzer Space Observatory, which detect the star at 4.5, 5.8, and 8 $\mu$m. The star looks rather unremarkable in the mid-IR. The excellent positional agreement of the IRAC source with the VLA free-free emission, the OH, CH$_3$OH, H$_2$O masers, and the dust continuum confirms that this is the most luminous object in the NGC 7538S core. The continuum emission at millimeter wavelengths is dominated by dust emission from the dense cold cloud core surrounding the protostar. Including all array configurations, the emission is dominated by an elliptical source with a size of ~ 8" x 3". If we filter out the extended emission we find three compact mm-sources inside the elliptical core. The strongest one, $S_A$, coincides with the VLA/IRAC source and resolves into a double source at 1.4 mm, where we have sub-arcsecond resolution. The measured spectral index, $\alpha$, between 3 and 1.4 mm is ~ 2.3, and steeper at longer wavelengths, suggesting a low dust emissivity or that the dust is optically thick. We argue that the dust in these accretion disks is optically thick and estimate a mass of an accretion disk or infalling envelope surrounding S$_A$ to be ~ 60 solar masses.

A wide field of view Cherenkov/fluorescence telescope array is one of the main components of the Large High Altitude Air Shower Observatory project. To serve as Cherenkov and fluorescence detectors, a flexible and mobile design is adopted for easy reconfiguring of the telescope array. Two prototype telescopes have been constructed and successfully run at the site of the ARGO-YBJ experiment in Tibet. The features and performance of the telescopes are presented.

We present the initial-final mass relation derived from 10 white dwarfs in wide binaries that consist of a main sequence star and a white dwarf. The temperature and gravity of each white dwarf was measured by fitting theoretical model atmospheres to the observed spectrum using a $\chi^{2}$ fitting algorithm. The cooling time and mass was obtained using theoretical cooling tracks. The total age of each binary was estimated from the chromospheric activity of its main sequence component to an uncertainty of about 0.17 dex in log \textit{t} The difference between the total age and white dwarf cooling time is taken as the main sequence lifetime of each white dwarf. The initial mass of each white dwarf was then determined using stellar evolution tracks with a corresponding metallicity derived from spectra of their main sequence companions, thus yielding the initial-final mass relation. Most of the initial masses of the white dwarf components are between 1 - 2 M$_{\odot}$. Our results suggest a correlation between the metallicity of a white dwarf's progenitor and the amount of post-main-sequence mass loss it experiences - at least among progenitors with masses in the range of 1 - 2 M$_{\odot}$. A comparison of our observations to theoretical models suggests that low mass stars preferentially lose mass on the red giant branch.

We present a catalog of 80 very wide fragile binary candidates (projected separations > 10000 AU) from the SDSS DR8 spectral archive. The pairs were selected based on proper motion, radial velocity, metallicity and photometric parallax criteria. The angular separations of these pairs range from 3" to 250". The peak in the metallicity distribution of these pairs is about -0.5 dex of solar metallicity. Space motions and reduced proper motion diagrams indicate all these pairs are members of the disk. The chromospheric activity index SHK of each component in 38 binary candidates having spectra of high signal-to-noise ratio and member stars of three open clusters (NGC2420, M67 and NGC6791) were measured. The SHK vs. color relation for these binary candidates is consistent with the trend seen in these open clusters. The ages implied by this relation suggest that fragile wide pairs can survive longer than 8 Gyr.

Light transport in graded index media follows a curved trajectory determined by the Fermat's principle. Besides the effect of variation of the refractive index on the transport of radiative intensity, the curved ray trajectory will induce geometrical effects on the transport of polarization ellipse. This paper presents a complete derivation of vector radiative transfer equation for polarized radiation transport in absorption, emission and scattering graded index media. The derivation is based on the analysis of the conserved quantities for polarized light transport along curved trajectory and a novel approach. The obtained transfer equation can be considered as a generalization of the classic vector radiative transfer equation that is only valid for uniform refractive index media. Several variant forms of the transport equation are also presented, which include the form for Stokes parameters defined with a fixed reference and the Eulerian forms in the ray coordinate and in several common orthogonal coordinate systems.

On the basis of a recently derived color-metallicity relation and stellar parameters from the Sloan Digtal Sky Survey Data Release 7 spectroscopic survey, a large sample of red horizontal branch candidates have been selected to serve as standard candles. The metallicity and kinematic distributions of these stars indicate that they mainly originate from the thick-disk and the halo populations. The typical thick disk is characterized by the first group peaking at [Fe/H] ~ -0.6, Vrot ~ 170 km/s with a vertical scale height around |Z| ~ 1.2 kpc, while stars with [Fe/H] < -0.9 are dominated by the halo population. Two sub-populations of the halo are suggested by the RHB stars peaking at [Fe/H] ~-1.3: one component with Vrot > 0 km/s (Halo I) shows a sign of metallicity gradient in the [Fe/H] versus |Z| diagram, while the other with Vrot < 0 km/s (Halo II) does not. The Halo I mainly clumps at the inner halo with R < 10 kpc and the Halo II comes both from the inner halo with R < 10 kpc and the outer halo with R > 10 kpc based on the star distribution in the R versus |Z| diagram.

Based on SDSS DR8, we have selected a sample of 1728 red horizontal branch stars with |Z|<3 kpc by using a color-metallicity relation and stellar parameters. The sample stars clearly trace a typical thick disk population with peaks at |Z|=1.26 kpc and [Fe/H]=-0.54. The vertical metallicity gradient of the thick disk is estimated in two ways. One is a fit to the Gaussian peaks of the metallicity histograms of the thick disk by subtracting minor contributions from the thin disk and the inner halo based on the Besancon Galaxy model. The resulting gradient is -0.12+-0.01 dex/kpc for 0.5 <|Z| <3 kpc. The other method is to linearly fit the data based on stars with 1 <|Z| <3 kpc being the main component of the thick disk. Five subgroups are then selected in different directions in the X-|Z| plane to investigate the difference in the vertical metallicity gradient between the Galactocenter and anti-Galactocenter directions. We found that a vertical gradient o-0.22+-0.07 dex/kpc is detected for five directions except for one involving the pollution of stars from the bulge. The results indicate that the vertical gradient is dominant, but a radial gradient has a minor contribution for the thick disk population represented by RHB stars with 1 <|Z| <3 kpc. The present work strongly suggests the existence of a metallicity gradient in the thick disk, which is thought to be negligible in most previous works in the literature.

The ionized core in the Sgr B2 Main star-forming region was imaged using the Submillimeter Array archival data observed for the H26$\alpha$ line and continuum emission at 0.86 millimeter with an angular resolution 0.3\arcsec. Eight hyper-compact H26$\alpha$ emission sources were detected with a typical size in the range of 1.6--20$\times10^2$ AU and electron density of 0.3--3$\times10^7$ cm$^{-3}$, corresponding to the emission measure 0.4--8.4$\times10^{10}$ cm$^{-6}$ pc. The H26$\alpha$ line fluxes from the eight hyper-compact HII sources imply that the ionization for each of the sources must be powered by a Lyman continuum flux from an O star or a cluster of B stars. The most luminous H26$\alpha$ source among the eight detected requires an O6 star that appears to be embedded in the ultra-compact HII region F3. In addition, $\sim$ 23 compact continuum emission sources were also detected within the central 5\arcsec$\times$3\arcsec\,($\sim0.2$ pc) region. In the assumption of a power-law distribution for the dust temperature, with the observed brightness temperature of the dust emission we determined the physical properties of the submillimeter emission sources showing that the molecular densities are in the range of 1--10$\times10^8$ cm$^{-3}$, surface densities between 13 to 150 $g$ cm$^{-2}$, and total gas masses in the range from 5 to $\gtrsim$ 200 $M_\odot$ which are 1 or 2 orders of magnitude greater than the corresponding values of the Bonnor-Ebert mass. With a mean free-fall time scale of 2$\times10^3$ y, each of the massive protostellar cores are undergoing gravitational collapse to form new massive stars in the Sgr B2 Main core.

We report new results on outflow and infall in the star forming cores W3-SE SMA-1 and SMA-2 based on analysis of $\sim2.5\arcsec$ resolution observations of the molecular lines HCN(3-2), HCO$^+$(3-2), N$_2$H$^+$(3-2) and CH$_3$OH(5$_{2,3}-4_{1,3}$) with the Submillimeter Array. A high-velocity bipolar outflow originating from the proto-stellar core SMA-1 was observed in the HCN(3-2) line, with a projected outflow axis in position angle 48$\arcdeg$. The detection of the outflow is confirmed from other molecular lines. An inverse P-Cygni profile in the HCN(3-2) line toward SMA-1 suggests that at least one of the double cores accretes matters from the molecular core. A filamentary structure in the molecular gas surrounds SMA-1 and SMA-2. Based on the SMA observations, our analysis suggests that the double pre-stellar cores SMA-1 and SMA-2 result from fragmentation in the collapsing massive molecular core W3-SE, and it is likely that they are forming intermediate to high-mass stars which will be new members of a star cluster in the W3-SE region.

Quasi-periodic, propagating fast mode magnetosonic waves in the corona were difficult to observe in the past due to relatively low instrument cadences. We report here evidence of such waves directly imaged in EUV by the new SDO AIA instrument. In the 2010 August 1 C3.2 flare/CME event, we find arc-shaped wave trains of 1-5% intensity variations (lifetime ~200 s) that emanate near the flare kernel and propagate outward up to ~400 Mm along a funnel of coronal loops. Sinusoidal fits to a typical wave train indicate a phase velocity of 2200 +/- 130 km/s. Similar waves propagating in opposite directions are observed in closed loops between two flare ribbons. In the k-$\omega$ diagram of the Fourier wave power, we find a bright ridge that represents the dispersion relation and can be well fitted with a straight line passing through the origin. This k-$\omega$ ridge shows a broad frequency distribution with indicative power at 5.5, 14.5, and 25.1 mHz. The strongest signal at 5.5 mHz (period 181 s) temporally coincides with quasi-periodic pulsations of the flare, suggesting a common origin. The instantaneous wave energy flux of $(0.1-2.6) \times 10^7 ergs/cm^2/s$ estimated at the coronal base is comparable to the steady-state heating requirement of active region loops.

1156+295 is a flat-spectrum quasar which is loud at radio and gamma-ray. Previous observations of the source revealed a radio morphology on pc to kpc scales consistent with a helical jet model. In our present research, this source was observed with the VLBA at 86, 43 and 15 GHz on four epochs from 10 May 2003 to 13 March 2005 aiming at studying the structure of the innermost jet in order to understand the relation between the helical structure and the astrophysical processes in the central engine. A core-jet structure with six jet components is identified. The apparent transverse velocities of the six jet components derived from proper motion measurements are in the range between 3.6 c and 11.6 c. The overall jet shape shows oscillatory morphology with multiple curvatures on pc scales which might be indicative of a helical pattern. Models of helical jet are discussed on the basis of both Kelvin-Helmholtz (K-H) instability and jet precession. The K-H instability model shows better agreement with the observed data. The overall radio structure on the scale from sub-pc to kpc appears to be fitted with a hydrodynamic model with the fundamental helical mode in Kelvin-Helmholtz (K-H) instability. This helical mode with an initial characteristic wavelength of 0.2 pc is excited at the base of the jet on the scale of 0.005 pc (or 1000R_s, the typical size of the broad line region for a super massive black hole of $4.3\times10^8M_{\odot}$). A presessing jet model can also fit the observed jet structure on the scale between 10 pc and 300 pc. However, additional astrophysical processes may be required for the presessing jet model in order to explain the bendings on the inner jet structure (1 to 10 pc) and re-collimation of the large scale jet outflow (>300 pc).

The Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory (SDO/HMI) provides continuous full-disk observations of solar oscillations. We develop a data-analysis pipeline based on the time-distance helioseismology method to measure acoustic travel times using HMI Doppler-shift observations, and infer solar interior properties by inverting these measurements. The pipeline is used for routine production of near-real-time full-disk maps of subsurface wave-speed perturbations and horizontal flow velocities for depths ranging from 0 to 20 Mm, every eight hours. In addition, Carrington synoptic maps for the subsurface properties are made from these full-disk maps. The pipeline can also be used for selected target areas and time periods. We explain details of the pipeline organization and procedures, including processing of the HMI Doppler observations, measurements of the travel times, inversions, and constructions of the full-disk and synoptic maps. Some initial results from the pipeline, including full-disk flow maps, sunspot subsurface flow fields, and the interior rotation and meridional flow speeds, are presented.

Based on high resolution and high signal-to-noise ratio (S/N) spectra analysis of 90 solar type stars, we have established several new metallicity calibrations in Teff range [5600, 6500] K based on red spectra with the wavelength range of 560-880 nm. The new metallicity calibrations are applied to determine the metallicity of solar analogs selected from SDSS spectra. There is a good consistent result with the adopted value presented in SDSS-DR7 and a small scatter of 0.26 dex for stars with S/N > 50 is obtained. This study provides a new reliable way to derive the metallicity for solar-like stars with low resolution spectra. In particular, our calibrations are useful for finding metal-rich stars, which are missing in SSPP.

We present the chromospheric activity (CA) levels, metallicities and full space motions for 41 F, G, K and M dwarf stars in 36 wide binary systems. Thirty-one of the binaries, contain a white dwarf component. In such binaries the total age can be estimated by adding the cooling age of the white dwarf to an estimate of the progenitor's main sequence lifetime. To better understand how CA correlates to stellar age, 14 cluster member stars were also observed. Our observations demonstrate for the first time that in general CA decays with age from 50 Myr to at least 8 Gyr for stars with 1.0 < V-I < 2.4. However, little change occurs in CA level for stars with V-I < 1.0 between 1 Gyr and 5 Gyr, consistent with the results of Pace et al. (2009). Our sample also exhibits a negative correlation between stellar age and metallicity, a positive correlation between stellar age and W space velocity component and the W velocity dispersion increases with age. Finally, the population membership of these wide binaries is examined based upon their U, V, W kinematics, metallicity and CA. We conclude that wide binaries are similar to field and cluster stars in these respects. More importantly, they span a much more continuous range in age and metallicity than is afforded by nearby clusters.

We present the first aperture synthesis unbiased spectral line survey toward an extragalactic object. The survey covered the 40 GHz frequency range between 202 and 242 GHz of the 1.3 mm atmospheric window. We find that 80% of the observed band shows molecular emission, with 73 features identified from 15 molecular species and 6 isotopologues. The 13C isotopic substitutions of HC3N and transitions from H2(18)O, 29SiO, and CH2CO are detected for the first time outside the Galaxy. Within the broad observed band, we estimate that 28% of the total measured flux is due to the molecular line contribution, with CO only contributing 9% to the overall flux. We present maps of the CO emission at a resolution of 2.9"x1.9" which, though not enough to resolve the two nuclei, recover all the single-dish flux. The 40 GHz spectral scan has been modelled assuming LTE conditions and abundances are derived for all identified species. The chemical composition of Arp 220 shows no clear evidence of an AGN impact on the molecular emission but seems indicative of a purely starburst-heated ISM. The overabundance of H2S and the low isotopic ratios observed suggest a chemically enriched environment by consecutive bursts of star formation, with an ongoing burst at an early evolutionary stage. The large abundance of water (~10^-5), derived from the isotopologue H2(18)O, as well as the vibrationally excited emission from HC3N and CH3CN are claimed to be evidence of massive star forming regions within Arp 220. Moreover, the observations put strong constraints on the compactness of the starburst event in Arp 220. We estimate that such emission would require ~2-8x10^6 hot cores, similar to those found in the Sgr B2 region in the Galactic center, concentrated within the central 700 pc of Arp 220.

We report on a new phenomenon of `alignment' of supergranulation cells in the polar regions of the Sun. Recent high-resolution datasets obtained by the Solar Optical Telescope onboard the Hinode satellite enabled us to investigate supergranular structures in high-latitude regions of the Sun. We have carried out a local helioseismology time-distance analysis of the data, and detected acoustic travel-time variations due to the supergranular flows. The supergranulation cells in both the north and south polar regions show systematic alignment patterns in the north-south direction. The south-pole datasets obtained in a month-long Hinode campaign indicate that the supergranulation alignment property may be quite common in the polar regions. We also discuss the latitudinal dependence of the supergranulation cell sizes; the data show that the east-west cell size decreases towards higher latitudes.

We investigate how helioseismic waves that originate from effective point sources interact with a sunspot. These waves are reconstructed from observed stochastic wavefields on the Sun by cross-correlating photospheric Doppler-velocity signals. We select the wave sources at different locations relative to the sunspot, and investigate the p- and f-mode waves separately. The results reveal a complicated picture of waveform perturbations caused by the wave interaction with the sunspot. In particular, it is found that for waves originating from outside of the sunspot, p-mode waves travel across the sunspot with a small amplitude reduction and slightly higher speed, and wave amplitude and phase get mostly restored to the quiet-Sun values after passing the sunspot. The f-mode wave experiences some amplitude reduction passing through the sunspot, and the reduced amplitude is not recovered after that. The wave-propagation speed does not change before encountering the sunspot and inside the sunspot, but the wavefront becomes faster than the reference wave after passing through the sunspot. For waves originating from inside the sunspot umbra, both f- and p-mode waves show significant amplitude reductions and faster speed during all courses of propagation. A comparison of positive and negative time lags of cross-correlation functions shows an apparent asymmetry in the waveform changes for both the f- and p-mode waves. We suggest that the waveform variations of the helioseismic waves interacting with a sunspot found in this article can be used for developing a method of waveform heliotomography, similar to the waveform tomography of the Earth.

Local helioseismology provides unique information about the subsurface structure and dynamics of sunspots and active regions. However, because of complexity of sunspot regions local helioseismology diagnostics require careful analysis of systematic uncertainties and physical interpretation of the inversion results. We present new results of comparison of the ring-diagram analysis and time-distance helioseismology for active region NOAA 9787, for which a previous comparison showed significant differences in the subsurface sound-speed structure, and discuss systematic uncertainties of the measurements and inversions. Our results show that both the ring-diagram and time-distance techniques give qualitatively similar results, revealing a characteristic two-layer seismic sound-speed structure consistent with the results for other active regions. However, a quantitative comparison of the inversion results is not straightforward. It must take into account differences in the sensitivity, spatial resolution and the averaging kernels. In particular, because of the acoustic power suppression, the contribution of the sunspot seismic structure to the ring-diagram signal can be substantially reduced. We show that taking into account this effect reduces the difference in the depth of transition between the negative and positive sound-speed variations inferred by these methods. Further detailed analysis of the sensitivity, resolution and averaging properties of the local helioseismology methods is necessary for consolidation of the inversion results. It seems to be important that both methods indicate that the seismic structure of sunspots is rather deep and extends to at least 20 Mm below the surface, putting constraints on theoretical models of sunspots.

The power of solar acoustic waves is reduced inside sunspots mainly due to absorption, emissivity reduction, and local suppression. The coefficients of these power-reduction mechanisms can be determined by comparing time-distance cross-covariances obtained from sunspots and from the quiet Sun. By analyzing 47 active regions observed by SOHO/MDI without using signal filters, we have determined the coefficients of surface absorption, deep absorption, emissivity reduction, and local suppression. The dissipation in the quiet Sun is derived as well. All of the cross-covariances are width corrected to offset the effect of dispersion. We find that absorption is the dominant mechanism of the power deficit in sunspots for short travel distances, but gradually drops to zero at travel distances longer than about 6 degrees. The absorption in sunspot interiors is also significant. The emissivity-reduction coefficient ranges from about 0.44 to 1.00 within the umbra and 0.29 to 0.72 in the sunspot, and accounts for only about 21.5% of the umbra's and 16.5% of the sunspot's total power reduction. Local suppression is nearly constant as a function of travel distance with values of 0.80 and 0.665 for umbrae and whole sunspots respectively, and is the major cause of the power deficit at large travel distances.

We report a study of the H30$\alpha$ line emission at 1.3 mm from the region around Sgr A* made with the Submillimeter Array at a resolution of 2\arcsec\ over a field of 60\arcsec\ (2 parsec) and a velocity range of -360 to +345 \kms. This field encompasses most of the Galactic center's "minispiral". With an isothermal homogeneous HII model, we determined the physical conditions of the ionized gas at specific locations in the Northern and Eastern Arms from the H30$\alpha$ line data along with Very Large Array data from the H92$\alpha$ line at 3.6 cm and from the radio continuum emission at 1.3 cm. The typical electron density and kinetic temperature in the minispiral arms are 3-21$\times10^4$ cm$^{-3}$ and 5,000-13,000 K, respectively. The H30$\alpha$ and H92$\alpha$ line profiles are broadened due to the large velocity shear within and along the beam produced by dynamical motions in the strong gravitational field near Sgr A*. We constructed a 3D model of the minispiral using the orbital parameters derived under the assumptions that the gas flows are in Keplerian motion. The gas in the Eastern Arm appears to collide with the Northern Arm flow in the "Bar" region, which is located 0.1-0.2 parsec south of and behind Sgr A*. Finally, a total Lyman continuum flux of $3\times10^{50}$ photons s$^{-1}$ is inferred from the assumption that the gas is photoionized and the ionizing photons for the high-density gas in the minispiral arms are from external sources, which is equivalent to $\sim250$ O9-type zero-age-main-sequence stars.

The Large High Altitude Air Shower Observatory project is proposed to study high energy gamma ray astronomy ( 40 GeV-1 PeV ) and cosmic ray physics ( 20 TeV-1 EeV ). The wide field of view Cherenkov telescope array, as a component of the LHAASO project, will be used to study energy spectrum and compositions of cosmic ray by measuring the total Cherenkov light generated by air showers and shower maximum depth. Two prototype telescopes have been in operation since 2008. The pointing accuracy of each telescope is crucial to the direction reconstruction of the primary particles. On the other hand the primary energy reconstruction relies on the shape of the Cherenkov image on the camera and the unrecorded photons due to the imperfect connections between photomultiplier tubes. UV bright stars are used as point-like objects to calibrate the pointing and to study the optical properties of the camera, the spot size and the fractions of unrecorded photons in the insensitive areas of the camera.

We present high resolution images of the bipolar outflow from W51e2, which are produced from the Submillimeter Array archival data observed for CO(3-2) and HCN(4-3) lines with angular resolutions of 0.8" x 0.6" and 0.3" x 0.2", respectively. The images show that the powerful outflow originates from the protostellar core W51e2-E rather than from the ultracompact HII region W51e2-W. The kinematic timescale of the outflow from W51e2-E is about 1000 yr, younger than the age (~5000 yr) of the ultracompact HII region W51e2-W. A large mass loss rate of ~1 x 10^{-3} M_sun yr^{-1} and a high mechanical power of 120 L_sun are inferred, suggesting that an O star or a cluster of B stars are forming in W51e2-E. The observed outflow activity along with the inferred large accretion rate indicates that at present W51e2-E is in a rapid phase of star formation.

We show that under certain conditions, subsurface structures in the solar interior can alter the average acoustic power observed at the photosphere above them. By using numerical simulations of wave propagation, we show that this effect is large enough for it to be potentially used for detecting emerging active regions before they appear on the surface. In our simulations, simplified subsurface structures are modeled as regions with enhanced or reduced acoustic wave speed. We investigate the dependence of the acoustic power above a subsurface region on the sign, depth, and strength of the wave speed perturbation. Observations from the Solar and Heliospheric Observatory/Michelson Doppler Imager (SOHO/MDI) prior and during the emergence of NOAA active region 10488 are used to test the use of acoustic power as a potential precursor of magnetic flux emergence.

We report the results from observations of H30$\alpha$ line emission in Sgr A West with the Submillimeter Array at a resolution of 2\arcsec and a field of view of about 40\arcsec. The H30$\alpha$ line is sensitive to the high-density ionized gas in the minispiral structure. We compare the velocity field obtained from H30$\alpha$ line emission to a Keplerian model, and our results suggest that the supermassive black hole at Sgr A* dominates the dynamics of the ionized gas. However, we also detect significant deviations from the Keplerian motion, which show that the impact of strong stellar winds from the massive stars along the ionized flows and the interaction between Northern and Eastern arms play significant roles in the local gas dynamics.

We present high-resolution continuum images of the W51e2 complex processed from archival data of the Submillimeter Array (SMA) at 0.85 and 1.3 mm and the Very Large Array (VLA) at 7 and 13 mm. We also made line images and profiles of W51e2 for three hydrogen radio recombination lines (H26\alpha, H53\alpha, and H66\alpha) and absorption of two molecular lines of HCN(4-3) and CO(2-1). At least four distinct continuum components have been detected in the 3" region of W51e2 from the SMA continuum images at 0.85 and 1.3 mm with resolutions of 0.3"x0.2" and 1.4"x0.7", respectively. The west component, W51e2-W, coincides with the UC HII region reported from previous radio observations. The H26\alpha line observation reveals an unresolved hyper-compact ionized core (<0.06" or <310 AU) with a high electron temperature of 1.2x10^4 K, with corresponding emission measure EM>7x10^{10} pc cm^{-6} and electron density N_e>7x10^6 cm^{-3}. The inferred Lyman continuum flux implies that the HII region W51e2-W requires a newly formed massive star, an O8 star or a cluster of B-type stars, to maintain the ionization. The east component, W51e2-E, has a total mass of ~140 M_{\sun} according to our SED analysis and a large infall rate of > 1.3x10^{-3} M_{\sun}yr^{-1} inferred from the absorption of HCN. W51e2-E appears to be the accretion center in W51e2 and to host one or more growing massive proto-stars. Located 2" northwest from W51e2-E, W51e2-NW is not detected in the continuum emission at \lambda>=7 mm. Along with the maser activities previously observed, our analysis suggests that W51e2-NW is at an earlier phase of star formation. W51e2-N is located 2" north of W51e2-E and has only been detected at 1.3 mm with a lower angular resolution (~1"), suggesting that it is a primordial, massive gas clump in the W51e2 complex.

We report new results from CARMA observations of both continuum and HCO+(1-0) emission at 3.4 mm from W3-SE, a molecular core of intermediate mass, together with the continuum observations at 1.1 and 0.85/0.45 mm with the SMA and JCMT. A continuum emission core elongated from SE to NW (~10"), has been observed at the and further resolved into a double source with the SMA at 1.1 mm, with a separation of ~4". Together with the measurements from the Spitzer and MSX at mid-IR, we determined the SED of W3-SE and fit it with a thermal dust emission model, suggesting the presence of two dust components with different temperatures. The emission at mm/submm wavelengths is dominated by a major cold (~41 K) with a mass of ~65 Msun. In addition, there is a weaker hot component (~400 K) which accounts for emission in the mid-IR, suggesting that a small fraction of dust has been heated by newly formed stars. We also imaged the molecular core in the HCO+(1-0) line using CARMA at an angular resolution ~6". With the CARMA observations, we have verified the presence of a blue-dominated double peak profile toward this core. The line profile cannot be explained by infall alone. The broad velocity wings of the line profile suggest that other kinematics such as outflows within the central 6" of the core likely dominate the resulting spectrum. The kinematics of the sub-structures of this core suggest that the molecular gas outside the main component appears to be dominated by the bipolar outflow originated from the dust core with a dynamical age of >30000 yr. Our analysis, based on the observations at wavelengths from mm/submm to mid-IR suggests that the molecular core W3-SE hosts a group of newly formed young stars and protostars.

We analyze a solar active region observed by the Hinode CaII H line using the time-distance helioseismology technique, and infer wave-speed perturbation structures and flow fields beneath the active region with a high spatial resolution. The general subsurface wave-speed structure is similar to the previous results obtained from SOHO/MDI observations. The general subsurface flow structure is also similar, and the downward flows beneath the sunspot and the mass circulations around the sunspot are clearly resolved. Below the sunspot, some organized divergent flow cells are observed, and these structures may indicate the existence of mesoscale convective motions. Near the light bridge inside the sunspot, hotter plasma is found beneath, and flows divergent from this area are observed. The Hinode data also allow us to investigate potential uncertainties caused by the use of phase-speed filter for short travel distances. Comparing the measurements with and without the phase-speed filtering, we find out that inside the sunspot, mean acoustic travel times are in basic agreement, but the values are underestimated by a factor of 20-40% inside the sunspot umbra for measurements with the filtering. The initial acoustic tomography results from Hinode show a great potential of using high-resolution observations for probing the internal structure and dynamics of sunspots.

The purpose of this work is to image solar far-side active regions using acoustic signals with three skips and improve the quality of existing images. The mapping of far-side active regions was first made possible using the helioseismic holography technique by use of four-skip acoustic signals. The quality of far-side images was later improved with the combination of four- and five-skip signals using the time-distance helioseismology technique. In this work, we explore the possibility of making three-skip far-side images of active regions, and improving the image quality by combining the three-skip images with the images obtained from existing techniques. A new method of combining images is proposed that increases the signal-to-noise ratio and reduces the appearance of spurious features.

Based on photometric data of the central parts of eight globular clusters and one open cluster presented by An and his collaborators, we select red horizontal branch (RHB) stars in the (g-r)0-g0 diagram and make a statistical study of the distributions of their colors and absolute magnitudes in the SDSS ugriz system. Meanwhile, absolute magnitudes in the Johnson VRI system are calculated through the translation formulae between gri and VRI in the literature. The calibrations of absolute magnitude as functions of metallicity and age are established by linear regressions of the data. It is found that metallicity coefficients in these calibrations decrease, while age coefficients increase, from the blue $u$ filter to the red $z$ filter. The calibration of Mi= 0.06[Fe/H]+0.040t+0.03 has the smallest scatter of 0.04 mag, and thus i is the best filter in the $ugriz$ system when RHB stars are used for distance indicators. The comparison of the MI calibration from our data with that from red clump stars indicates that the previous suggestion that the $I$ filter is better than the V filter in distance determination may not be true because of its significant dependence on age.

The solar tachocline at the bottom of the convection zone is an important region for the dynamics of the Sun and the solar dynamo. In this region, the sound speed inferred by global helioseismology exhibits a bump of approximately 0.4% relative to the standard solar model. Global helioseismology does not provide any information on possible latitudinal variations or asymmetries between the Northern and Southern hemisphere. Here, we develop a time-distance helioseismology technique, including surface- and deep-focusing measurement schemes and a combination of both, for two-dimensional tomographic imaging of the solar tachocline that infers radial and latitudinal variations in the sound speed. We test the technique using artificial solar oscillation data obtained from numerical simulations. The technique successfully recovers major features of the simplified tachocline models. The technique is then applied to SOHO/MDI medium-l data and provides for the first time a full two-dimensional sound-speed perturbation image of the solar tachocline. The one-dimensional radial profile obtained by latitudinal averaging of the image is in good agreement with the previous global helioseismology result. It is found that the amplitude of the sound-speed perturbation at the tachocline varies with latitude, but it is not clear whether this is in part or fully an effect of instrumental distortion. Our initial results demonstrate that time-distance helioseismology can be used to probe the deep interior structure of the Sun, including the solar tachocline.

We present new results based on high-resolution observations of Sgr A West at the Galactic center with the VLA at 1.3 cm. We measured proper motions for 71 compact HII components. We also investigated radial velocities in the LSR velocity using the H92a line data. Combining proper motion and radial velocity measurements, we have determined the 3D velocity distribution in Sgr A West. We find that the three ionized streams (Northern Arm, Eastern Arm, and Western Arc) can be modeled with three bundles of Keplerian orbits around Sgr A*. We determined the five orbital parameters for each of them using LSQ fitting to the locii of the streams. Our results confirm earlier results on the streams in the Western Arc and the Northern Arm to be in Keplerian orbits, suggesting that the stream in the Eastern Arm is also consistent with an elliptical orbit. Both the Northern and Eastern Arm streams have high eccentricities, while the Western Arc stream is nearly circular. All three streams orbit around Sgr A* in a counterclockwise sense (viewed from the Earth). We also report an ionized nebula associated with IRS 8, including a bow shock in radio continuum emission which shows excellent agreement with near IR observations. From the H92a line data, we find evidence for interaction between the IRS 8 nebula and the Northern Arm stream. Other new morphological features revealed in our high-resolution image include: 1) a helical structure in the Northern Arm, suggesting that MHD plays an important role in the motion of the ionized gas, in addition to the dynamics determined by the central gravitational field and 2) a linear feature in the IRS 16 region, suggesting the compressed edge of the Northern Arm may result from the collective winds and radiation pressure from the high mass stars in the IRS16 cluster.

We report on a signature of chromospheric downflows in two emerging-flux regions detected by time-distance helioseismology analysis. We use both chromospheric intensity oscillation data in the Ca II H line and photospheric Dopplergrams in the Fe I 557.6nm line obtained by Hinode/SOT for our analyses. By cross-correlating the Ca II oscillation signals, we have detected a travel-time anomaly in the plage regions; outward travel times are shorter than inward travel times by 0.5-1 minute. However, such an anomaly is absent in the Fe I data. These results can be interpreted as evidence of downflows in the lower chromosphere. The downflow speed is estimated to be below 10 km/s. This result demonstrates a new possibility of studying chromospheric flows by time-distance analysis.

Based on the kernel estimator and wavelet technique, we have identified 22 moving group candidates in the solar neighborhood from a sample which includes around 14,000 dwarfs and 6000 giants. Six of them were previously known as the Hercules stream, the Sirus-UMa stream, the Hyades stream, the Caster group, the Pleiades stream, and the IC 2391; five of them have also been reported by other authors. 11 moving group candidates, not previously reported in the literature, showprominent structures in dwarf or giant samples.Acatalog of moving group candidates in the solar neighborhood is presented in this work.

The complex dynamics that lead to the emergence of active regions on the sun are poorly understood. One possibility is that magnetic structures (flux tubes, etc.) rise from below the surface by self induction and convection that lead to the formation of active regions and sunspots on the solar surface. For space weather forecasting, one would like to detect the subsurface structures before they reach the surface. The goal of this study is to investigate whether sound speed perturbations associated with subsurface structures could affect the acoustic power observed at the solar surface above them. Possible mechanisms for this effect are wave reflection, scattering or diffraction. By using numerical simulations of wave propagation in the solar interior, we investigate whether observations of the acoustic power can be used to detect emerging active regions before they appear on the surface. In the simulations, subsurface structures are modeled as regions with enhanced or reduced acoustic wavespeed. We show how the acoustic power above a subsurface region depends on the sign, depth and strength of the wavespeed perturbation. For comparison, we analyze observations from SOHO/MDI of the emergence of solar active region NOAA 10488.

Energetic flares are observed in the Galactic supermassive black hole Sagittarius A* from radio to X-ray wavelengths. On a few occasions, simultaneous flares have been detected in IR and X-ray observations, but clear counterparts at longer wavelengths have not been seen. We present a flare observed over several hours on 2006 July 17 with the Chandra X-Ray Observatory, the Keck II telescope, the Caltech Submillimeter Observatory, and the Submillimeter Array. All telescopes observed strong flare events, but the submillimeter peak is found to occur nearly 100 minutes after the X-ray peak. Submillimeter polarization data show linear polarization in the excess flare emission, increasing from 9% to 17% as the flare passes through its peak, consistent with a transition from optically thick to thin synchrotron emission. The temporal and spectral behavior of the flare require that the energetic electrons responsible for the emission cool faster than expected from their radiative output. This is consistent with adiabatic cooling in an expanding emission region, with X-rays produced through self-Compton scattering, although not consistent with the simplest model of such expansion. We also present a submillimeter flare that followed a bright IR flare on 2005 July 31. Compared to 2006, this event had a larger peak IR flux and similar submillimeter flux, but it lacked measurable X-ray emission. It also showed a shorter delay between the IR and submillimeter peaks. Based on these events we propose a synchrotron and self-Compton model to relate the submillimeter lag and the variable IR/X-ray luminosity ratio.

The aim of this paper is to extend our previous study of the solar-cycle variations of the meridional flows and to investigate their latitudinal and longitudinal structure in the subphotospheric layer, especially their variations in magnetic regions. Helioseismology observations indicate that mass flows around active regions are dominated by inflows into those regions. On average, those local flows are more important around leading magnetic polarities of active regions than around the following polarities, and depend on the evolutionary stage of particular active regions. We present a statistical study based on MDI/SOHO observations of 1996-2002 and show that this effect explains a significant part of the cyclic change of meridional flows in near-equatorial regions, but not at higher latitudes. A different mechanism driving solar-cycle variations of the meridional flow probably operates.

Far-side images of solar active regions have become one of the routine products of helioseismic observations, and are of importance for space weather forecasting by allowing the detection of sunspot regions before they become visible on the Earth side of the Sun. An accurate assessment of the quality of the far-side maps is difficult, because there are no direct observations of the solar far side to verify the detections. In this paper we assess far-side imaging based on the time-distance helioseismology method, by using numerical simulations of solar oscillations in a spherical solar model. Localized variations in the speed of sound in the surface and subsurface layers are used to model the perturbations associated with sunspots and active regions. We examine how the accuracy of the resulting far-side maps of acoustic travel times depends on the size and location of active regions. We investigate potential artifacts in the far-side imaging procedure, such as those caused by the presence of active regions on the solar near side, and suggest how these artifacts can be identified in the real Sun far-side images obtained from SOHO/MDI and GONG data.

Aims: We searched for the molecular outflows from fifteen molecular clouds associated with ultra-compact and compact HII (UCHII/HII) regions and discussed possible gas heating mechanism. Methods: Mapping observations of CO $J=3-2$ and $J=2-1$ lines were carried out with the KOSMA 3m-Telescope towards the 15 HII regions/molecular cloud complexes. Results: Ten molecular outflows were identified out of the fifteen HII region/molecular cloud complexes. The higher outflow detection rate (67%) suggested that such outflows are as common in high mass star forming regions as those in low mass star forming regions, which is consistent with the results of other authors. The observations also showed that the outflow might occur in the HII region. The integrated CO line intensity ratios (${\rm R_{I_{CO(3-2)}/I_{CO(2-1)}}}$) were determined from the core component of the spectra as well as from both the blue and red wings. Maximum line intensity ratios from the wings and core components appeared to be related to the mid-infrared sources imaged by Midcourse Space Experiment (MSX). The relationship between the maximum line intensity ratios and MSX sources indicates that the molecular gas could be heated by the emission of dust associated with massive stars. Based on maser observations reported in the literature, we found that ${\rm H_{2}O}$ masers were only detected in seven regions. The ${\rm H_{2}O}$ masers in these regions are located near the MSX sources and within the maximum line intensity ratio regions, suggesting that ${\rm H_{2}O}$ masers occur in relatively warm environments.

By analyzing numerically simulated solar oscillation data, we study the influence of non-uniform distribution of acoustic wave amplitude, acoustic source strength, and perturbations of the sound speed on the shifts of acoustic travel times measured by the time-distance helioseismology method. It is found that for short distances, the contribution to the mean travel time shift caused by non-uniform distribution of acoustic sources in sunspots may be comparable to (but smaller than) the contribution from the sound speed perturbation in sunspots, and that it has the opposite sign to the sound-speed effect. This effect may cause some underestimation of the negative sound-speed perturbations in sunspots just below the surface, that was found in previous time-distance helioseismology inferences. This effect cannot be corrected by artificially increasing the amplitude of oscillations in sunspots. For large time-distance annuli, the non-uniform distribution of wavefields does not have significant effects on the mean travel times, and thus the sound-speed inversion results. The measured travel time differences, which are used to determine the mass flows beneath sunspots, can also be systematically shifted by this effect, but only in an insignificant magnitude.

Observations of two H$_2$CO ($3_{03}-2_{02}$ and $3_{21}-2_{20}$) lines and continuum emission at 1.3 mm towards Sgr B2(N) and Sgr B2(M) have been carried out with the SMA. The mosaic maps of Sgr B2(N) and Sgr B2(M) in both continuum and lines show a complex distribution of dust and molecular gas in both clumps and filaments surrounding the compact star formation cores. We have observed a decelerating outflow originated from the Sgr B2(M) core, showing that both the red-shifted and blue-shifted outflow components have a common terminal velocity. This terminal velocity is 58$\pm$2 km s$^{-1}$. It provides an excellent method in determination of the systematic velocity of the molecular cloud. The SMA observations have also shown that a large fraction of absorption against the two continuum cores is red-shifted with respect to the systematic velocities of Sgr B2(N) and Sgr B2(M), respectively, suggesting that the majority of the dense molecular gas is flowing into the two major cores where massive stars have been formed. We have solved the radiative transfer in a multi-level system with LVG approximation. The observed H$_2$CO line intensities and their ratios can be adequately fitted with this model for the most of the gas components. However, the line intensities between the higher energy level transition H$_2$CO ($3_{21}-2_{20}$) and the lower energy level transition H$_2$CO ($3_{03}-2_{02}$) is reversed in the red-shifted outflow region of Sgr B2(M), suggesting the presence of inversion in population between the ground levels in the two K ladders (K$_{-1}$= 0 and 2). The possibility of weak maser processes for the H$_2$CO emission in Sgr B2(M) is discussed.

We review our work on Galactic open clusters in recent years, and introduce our proposed large program for the LOCS (LAMOST Open Cluster Survey). First, based on the most complete open clusters sample with metallicity, age and distance data as well as kinematic information, some preliminary statistical analysis regarding the spatial and metallicity distributions is presented. In particular, a radial abundance gradient of - 0.058$\pm$ 0.006 dex kpc$^{-1}$ was derived, and by dividing clusters into age groups we show that the disk abundance gradient was steeper in the past. Secondly, proper motions, membership probabilities, and velocity dispersions of stars in the regions of two very young open clusters are derived. Both clusters show clear evidence of mass segregation, which provides support for the ``primordial'' mass segregation scenarios. Based on the great advantages of the forthcoming LAMOST facility, we have proposed a detailed open cluster survey with LAMOST (the LOCS). The aim, feasibility, and the present development of the LOCS are briefly summarized.

We use the magnetic butterfly diagram to determine the speed of the magnetic flux transport on the solar surface towards the poles. The manifestation of the flux transport is clearly visible as elongated structures extended from the sunspot belt to the polar regions. The slopes of these structures are measured and interpreted as meridional magnetic flux transport speed. Comparison with the time-distance helioseismology measurements of the mean speed of the meridional flows at the depth of 3.5--12 Mm shows a generally good agreement, but the speeds of the flux transport and the meridional flow are significantly different in areas occupied by the magnetic field. The local circulation flows around active regions, especially the strong equatorward flows on the equatorial side of active regions affect the mean velocity profile derived by helioseismology, but do not influence the magnetic flux transport. The results show that the mean longitudinally averaged meridional flow measurements by helioseismology may not be used directly in solar dynamo models for describing the magnetic flux transport, and that it is necessary to take into account the longitudinal structure of these flows.

Results from initial helioseismic observations by Solar Optical Telescope onboard Hinode are reported. It has been demonstrated that intensity oscillation data from Broadband Filter Imager can be used for various helioseismic analyses. The k-omega power spectra, as well as corresponding time-distance cross-correlation function that promises high-resolution time-distance analysis below 6-Mm travelling distance, were obtained for G-band and CaII-H data. Subsurface supergranular patterns have been observed from our first time-distance analysis. The results show that the solar oscillation spectrum is extended to much higher frequencies and wavenumbers, and the time-distance diagram is extended to much shorter travel distances and times than they were observed before, thus revealing great potential for high-resolution helioseismic observations from Hinode.

We present high angular resolution (0.7'') observations made with the Very Large Array (VLA) of the radio recombination line (RRL) H53alpha and radio continuum emission at 43 GHz from the ultraluminous infrared galaxy (ULIRG) Arp 220. The 43 GHz continuum emission shows a compact structure (~2'') with two peaks separated by ~1'', the East (E) and West (W) components, that correspond to each galactic nucleus of the merger. The spectral indices for both the E and W components, using radio continuum images at 8.3 and 43 GHz are typical of synchrotron emission (alpha ~ -1.0). Our 43 GHz continuum and H53alpha line observations confirm the flux densities predicted by the models proposed by Anantharamaiah et al. This agreement with the models implies the presence of high-density (~ 100,000 cm^-3) compact HII regions (~ 0.1 pc) in Arp 220. The integrated H53alpha line emission is stronger toward the non-thermal radio continuum peaks, which are also coincident with the peaks of molecular emission of the H2CO. The coincidence between the integrated H53alpha and the H2CO maser line emission suggests that the recent star forming regions, traced by the high density gas, are located mainly in regions that are close to the two radio continuum peaks. A velocity gradient of ~ 0.30 km/s/pc in the H53alpha RRL is observed toward the E component and a second velocity gradient of ~ 0.15 km/s/pc is detected toward the W component. The orientations of these velocity gradients are in agreement with previous CO, HI and OH observations. The kinematics of the high-density ionized gas traced by the H53alpha line are consistent with two counter rotating disks as suggested by the CO and HI observations.

It is of great importance to monitor large solar active regions in the far-side of the Sun for space weather forecast, in particular, to predict their appearance before they rotate into our view from the solar east limb. Local helioseismology techniques, including helioseismic holography and time-distance, have successfully imaged solar far-side active regions. In this Letter, we further explore the possibility of imaging and improving the image quality of solar far-side active regions by use of time-distance helioseismology. In addition to the previously used scheme with four acoustic signal skips, a five-skip scheme is also included in this newly developed technique. The combination of both four- and five-skip far-side images significantly enhances the signal-to-noise ratio in the far-side images, and reduces spurious signals. The accuracy of the far-side active region imaging is also assessed using one whole year solar observation.

We derive the proper motions, membership probabilities, and velocity dispersions of stars in the regions of the young (about 2-4 Myr-old) open clusters NGC 2244 (the central cluster in the Monoceros R2 association) and NGC 6530 (the dominant cluster in the Sgr OB1 association) from photographic plate material obtained at Shanghai Astronomical Observatory, with time baselines of 34 and 87 years, respectively. Both clusters show clear evidence of mass segregation, but they do not exhibit any significant velocity-mass (or, equivalently, a velocity-luminosity) dependence. This provides strong support for the suggestion that the observed mass segregation is -- at least partially -- due to the way in which star formation has proceeded in these complex star-forming regions (``primordial'' mass segregation). Based on arguments related to the clusters' published initial mass functions, in conjunction with our new measurements of their internal velocity dispersions (35 and 8 km/s for NGC 2244 and NGC 6530, respectively), we provide strong arguments in favor of the dissolution of NGC 2244 on very short time-scales, while we speculate that NGC 6530 may be more stable against the effects of internal two-body relaxation. However, this latter object may well be destroyed by the strong tidal field prevalent at its location in the Galactic plane in the direction of the Galactic Center.

We apply time-distance helioseismology, local correlation tracking and Fourier spatial-temporal filtering methods to realistic supergranule scale simulations of solar convection and compare the results with high-resolution observations from the SOHO Michelson Doppler Imager (MDI). Our objective is to investigate the surface and sub-surface convective structures and test helioseismic measurements. The size and grid of the computational domain are sufficient to resolve various convective scales from granulation to supergranulation. The spatial velocity spectrum is approximately a power law for scales larger than granules, with a continuous decrease in velocity amplitude with increasing size. Aside from granulation no special scales exist, although a small enhancement in power at supergranulation scales can be seen. We calculate the time-distance diagram for f- and p-modes and show that it is consistent with the SOHO/MDI observations. From the simulation data we calculate travel time maps for surface gravity waves (f-mode). We also apply correlation tracking to the simulated vertical velocity in the photosphere to calculate the corresponding horizontal flows. We compare both of these to the actual large-scale (filtered) simulation velocities. All three methods reveal similar large scale convective patterns and provide an initial test of time-distance methods.

We present a direct comparison between two different techniques time-distance helioseismology and a local correlation tracking method for measuring mass flows in the solar photosphere and in a near-surface layer: We applied both methods to the same dataset (MDI high-cadence Dopplergrams covering almost the entire Carrington rotation 1974) and compared the results. We found that after necessary corrections, the vector flow fields obtained by these techniques are very similar. The median difference between directions of corresponding vectors is 24 degrees, and the correlation coefficients of the results for mean zonal and meridional flows are 0.98 and 0.88 respectively. The largest discrepancies are found in areas of small velocities where the inaccuracies of the computed vectors play a significant role. The good agreement of these two methods increases confidence in the reliability of large-scale synoptic maps obtained by them.

Recent progress in realistic simulations of solar convection have given us an unprecedented opportunity to evaluate the robustness of solar interior structures and dynamics obtained by methods of local helioseismology. We present results of testing the time-distance method using realistic simulations. By computing acoustic wave propagation time and distance relations for different depths of the simulated data, we confirm that acoustic waves propagate into the interior and then turn back to the photosphere. This demonstrates that in the numerical simulations properties of acoustic waves (p-modes) are similar to the solar conditions, and that these properties can be analyzed by the time-distance technique. For the surface gravity waves (f-mode), we calculate perturbations of their travel times, caused by localized downdrafts, and demonstrate that the spatial pattern of these perturbations (representing so-called sensitivity kernels) is similar to the patterns obtained from the real Sun, displaying characteristic hyperbolic structures. We then test the time-distance measurements and inversions by calculating acoustic travel times from a sequence of vertical velocities at the photosphere of the simulated data, and inferring a mean 3D flow fields by performing inversion based on the ray approximation. The inverted horizontal flow fields agree very well with the simulated data in subsurface areas up to 3 Mm deep, but differ in deeper areas. Due to the cross-talk effects between the horizontal divergence and downward flows, the inverted vertical velocities are significantly different from the mean convection velocities of the simulation dataset.

The millimeter/submillimeter wavelength polarization of Sgr A* is known to be variable in both magnitude and position angle on time scales down to a few hours. The unstable polarization has prevented measurements made at different frequencies and different epochs from yielding convincing measurements of Faraday rotation in this source. Here we present observations made with the Submillimeter Array polarimeter at 227 and 343 GHz with sufficient sensitivity to determine the rotation measure at each band without comparing position angles measured at separate epochs. We find the 10-epoch mean rotation measure to be (-5.6+/-0.7)x10^5 rad/m^2; the measurements are consistent with a constant value. We conservatively assign a 3sigma upper limit of 2x10^5 rad/m^2 to rotation measure changes, which limits accretion rate fluctuations to 25%. This rotation measure detection limits the accretion rate to less than 2x10^-7 M_sun/yr if the magnetic field is near equipartition, ordered, and largely radial, while a lower limit of 2x10^-9 M_sun/yr holds even for a sub-equipartition, disordered, or toroidal field. The mean intrinsic position angle is 167+/-7 degrees and we detect variations of 31(+18/-9) degrees. These variations must originate in the submillimeter photosphere, rather than arising from rotation measure changes.

Observations of H$_{2}$CO lines and continuum at 1.3 mm towards Sgr B2(N) and Sgr B2(M) cores were carried out with the SMA. We imaged H$_{2}$CO line absorption against the continuum cores and the surrounding line emission clumps. The results show that the majority of the dense gas is falling into the major cores where massive stars have been formed. The filaments and clumps of the continuum and gas are detected outside of Sgr B2(N) and Sgr B2(M) cores. Both the spectra and moment analysis show the presence of outflows from Sgr B2(M) cores. The H$_{2}$CO gas in the red-shifted outflow of Sgr B2(M) appears to be excited by a non-LTE process which might be related to the shocks in the outflow.

We report on the submillimeter properties of Sagittarius A* derived from observations with the Submillimeter Array and its polarimeter. We find that the spectrum of Sgr A* between 230 and 690 GHz is slightly decreasing when measured simultaneously, indicating a transition to optically thin emission around 300-400 GHz. We also present very sensitive and well calibrated measurements of the polarization of Sgr A* at 230 and 345 GHz. With these data we are able to show for the first time that the polarization of Sgr A* varies on hour timescales, as has been observed for the total intensity. On one night we find variability that may arise from a polarized "blob" orbiting the black hole. Finally, we use the ensemble of observations to determine the rotation measure. This represents the first statistically significant rotation measure determination and the only one made without resorting to comparing position angles measured at separate epochs. We find a rotation measure of (-5.6+/-0.7)x10^5 rad/m^2, with no evidence for variability on inter-day timescales at the level of the measurement error. The stability constrains interday fluctuations in the accretion rate to 8%. The mean intrinsic polarization position angle is 167+/-7 degrees and we detect variations of 31+18/-9 degrees. This separation of intrinsic polarization changes and possible rotation measure fluctuations is now possible because of the frequency coverage and sensitivity of our data. The observable rotation measure restricts the accretion rate to the range 2x10^{-7} Msun/yr to 2x10^{-9} Msun/yr, if the magnetic field is near equipartition and ordered.

It is well known that the observed amplitude of solar oscillations is lower in sunspots than in quiet regions of the Sun. We show that this local reduction in oscillation amplitudes combined with the phase-speed filtering procedure in time-distance helioseismic analyses could be a source of systematic errors in the range of 5 - 40% in the measured travel-time anomalies of acoustic waves around sunspots. Removing these travel time artifacts is important for correctly inferring the subsurface structure of sunspots. We suggest an empirical correction procedure and illustrate its usage for a small sunspot. This work utilizes data from MDI/SOHO.

We present new Very Large Array (VLA) observations toward NGC 253 of the recombination line H53alpha (43 GHz) at an angular resolution of 1.''5 X 1.''0. The free-free emission at 43 GHz is estimated to be ~100 mJy, implying a star formation rate of ~1.3 Msun/yr in the nuclear region of this starburst galaxy. A reanalysis is made for previously reported H92alpha observations carried out with angular resolution of 1.''5 X 1.''0 and 0.''36 X 0.''21. Based on the line and continuum emission models used for the 1.''5 X 1.''0 angular resolution observations, the RRLs H53alpha and H92alpha are tracers of the high-density (~10^5 cm-3) and low-density (~10^3 cm^-3) thermally ionized gas components in NGC 253, respectively. The velocity fields observed in the H53alpha and H92alpha lines 1.''5 X 1.''0) are consistent. The velocity gradient in the central ~18 pc of the NE component, as observed in both the H53alpha and H92alpha lines, is in the opposite direction to the velocity gradient determined from the CO observations. The enclosed virial mass, as deduced from the H53alpha velocity gradient over the NE component, is ~5 X 10^6 Msun in the central ~18 pc region. The H92alpha line observations at high angular resolution (0.''36 X 0.''21) reveal a larger velocity gradient, along a P.A. ~ -45 deg on the NE component, of ~110 km/s/arcsec. The dynamical mass estimated using the high angular resolution H92alpha data (~7 X 10^6 Msun) supports the existence of an accreted massive object in the nuclear region of NGC 253.

Recent observations of helioseismic holography revealed that magnetic fields inclined to the line-of-sight direction could cause systematic variations in measured acoustic phase shifts, and that the surface magnetic field may shift the phases and impair the coherence of acoustic waves (known as "showerglass effect"). We examine how these affect time-distance helioseismology measurements. It is confirmed that the inclined magnetic field could cause variations in acoustic travel times inside sunspot penumbra, however, inversions of the measured times show that this effect only slightly shifts the location of negative sound-speed variations near the solar surface, but does not change the inverted deeper interior structures. Regarding to the showerglass effect, we find that outgoing and incoming travel time perturbations through sunspots are significantly smaller than those reported from helioseismic holography. In addition, our second-skip cross-correlation experiments demonstrate that inside sunspots, the half of the double-skip travel times are very similar to the mean single-skip travel times. We conclude that these surface magnetism effects do not cause considerable systematic errors in time-distance helioseismology of active regions.

Although it is widely accepted that most galaxies have supermassive black holes (SMBHs) at their centers^{1-3}, concrete proof has proved elusive. Sagittarius A* (Sgr A*)^4, an extremely compact radio source at the center of our Galaxy, is the best candidate for proof^{5-7}, because it is the closest. Previous Very Long Baseline Interferometry (VLBI) observations (at 7mm) have detected that Sgr A* is ~2 astronomical unit (AU) in size^8, but this is still larger than the "shadow" (a remarkably dim inner region encircled by a bright ring) arising from general relativistic effects near the event horizon^9. Moreover, the measured size is wavelength dependent^{10}. Here we report a radio image of Sgr A* at a wavelength of 3.5mm, demonstrating that its size is \~1 AU. When combined with the lower limit on its mass^{11}, the lower limit on the mass density is 6.5x10^{21} Msun pc^{-3}, which provides the most stringent evidence to date that Sgr A* is an SMBH. The power-law relationship between wavelength and intrinsic size (The size is proportional to wavelength^{1.09}), explicitly rules out explanations other than those emission models with stratified structure, which predict a smaller emitting region observed at a shorter radio wavelength.

We report new simultaneous near-infrared/sub-millimeter/X-ray observations of the SgrA* counterpart associated with the massive 3-4x10**6 solar mass black hole at the Galactic Center. The main aim is to investigate the physical processes responsible for the variable emission from SgrA*. The observations have been carried out using the NACO adaptive optics (AO) instrument at the European Southern Observatory's Very Large Telescope and the ACIS-I instrument aboard the Chandra X-ray Observatory as well as the Submillimeter Array SMA on Mauna Kea, Hawaii, and the Very Large Array in New Mexico. We detected one moderately bright flare event in the X-ray domain and 5 events at infrared wavelengths.

Using the Submillimeter Array, we have made the first high angular resolution measurements of the linear polarization of Sagittarius A* at submillimeter wavelengths, and the first detection of intra-day variability in its linear polarization. We detected linear polarization at 340 GHz (880um) at several epochs. At the typical resolution of 1.4"x2.2", the expected contamination from the surrounding (partially polarized) dust emission is negligible. We found that both the polarization fraction and position angle are variable, with the polarization fraction dropping from 8.5% to 2.3% over three days. This is the first significant measurement of variability in the linear polarization fraction in this source. We also found variability in the polarization and total intensity within single nights, although the relationship between the two is not clear from these data. The simultaneous 332 and 342 GHz position angles are the same, setting a one-sigma rotation measure (RM) upper limit of 7x10^5 rad/m^2. From position angle variations and comparison of "quiescent" position angles observed here and at 230 GHz we infer that the RM is a few times 10^5 rad/m^2, a factor of a few below our direct detection limit. A generalized model of the RM produced in the accretion flow suggests that the accretion rate at small radii must be low, below 10^{-6} to 10^{-7} M_{sun}/year depending on the radial density and temperature profiles, but in all cases below the gas capture rate inferred from X-ray observations.

We observed Sgr A* using the Very Large Array (VLA) and the Giant Metrewave Radio Telescope (GMRT) at multiple cm and mm wavelengths on 17 June 2003. The measured flux densities of Sgr A*, together with those obtained from the Submillimeter Array (SMA) and the Keck II 10 m telescope on the same date, are used to construct a simultaneous spectrum of Sgr A* from 90 cm to 3.8 $\mu$m. The simultaneous spectrum shows a spectral break at about 3.6 cm, a possible signature of synchrotron self-absorption of the strong radio outburst which occurred near epoch 17 July 2003. At 90 cm, the flux density of Sgr A* is $0.22\pm0.06$ Jy, suggesting a sharp decrease in flux density at wavelengths longer than 47 cm. The spectrum at long cm wavelengths appears to be consistent with free-free absorption by a screen of ionized gas with a cutoff $\sim$100 cm. This cutoff wavelength appears to be three times longer than that of $\sim$30 cm suggested by Davies, Walsh & Booth based on observations in 1974 and 1975. Our analysis suggests that the flux densities of Sgr A* at wavelengths longer than 30 cm could be attenuated and modulated by stellar winds from massive stars close to Sgr A*.

Submillimeter Array observations of Orion-KL at 1'' resolution in the 440mu/690GHz band reveal new insights about the continuum and line emission of the region. The 440mu continuum flux density measurement from source I allows us to differentiate among the various proposed physical models: Source I can be well modeled by a ``normal'' protostellar SED consisting of a proton-electron free-free emission component at low frequencies and a strong dust component in the submillimeter bands. Furthermore, we find that the protostellar object SMA1 is clearly distinct from the hot core. The non-detection of SMA1 at cm and infrared wavelengths suggests that it may be one of the youngest sources in the entire Orion-KL region. The molecular line maps show emission mainly from the sources I, SMA1 and the hot core peak position. An analysis of the CH$3CN(37_K-36_K) K-ladder (K=0...3) indicates a warm gas component of the order 600+-200K. In addition, we detect a large fraction (~58%) of unidentified lines and discuss the difficulties of line identifications at these frequencies.

The central 2x1 kpc of the starburst galaxy NGC 253 has been imaged using the Submillimeter Array at a 60 pc resolution in the J=2-1 transitions of 12CO, 13CO, and C18O as well as in the 1.3 mm continuum. Molecular gas and dust are mainly in the circumnuclear disk of ~500 pc radius, with warm (~40 K) and high area-filling factor gas in its central part. Two gas shells or cavities have been discovered in the circumnuclear disk. They have ~100 pc diameters and have large velocity widths of 80-100 km/s, suggestive of expansion at ~50 km/s. Modeled as an expanding bubble, each shell has an age of ~0.5 Myr and needed kinetic energy of ~1E46 J as well as mean mechanical luminosity of ~1E33 W for its formation. The large energy allows each to be called a superbubble. A ~10^6 Msun super star cluster can provide the luminosity, and could be a building block of the nuclear starburst in NGC 253. Alternatively, a hypernova can also be the main source of energy for each superbubble, not only because it can provide the mechanical energy and luminosity but also because the estimated rate of superbubble formation and that of hypernova explosions are comparable. Our observations indicate that the circumnuclear molecular disk harboring the starburst is highly disturbed on 100 pc or smaller scales, presumably by individual young clusters and stellar explosions, in addition to globally disturbed in the form of the well-known superwind.

The Submillimeter Array (SMA), a collaboration between the Smithsonian Astrophysical Observatory and the Academica Sinica Institute for Astronomy and Astrophysics of Taiwan, is an eight-element radio-interferometer designed to operate throughout the major atmospheric windows from about 180 to 900 GHz. In an effort to mitigate the effects of atmospheric instabilities which limit the phase coherence of the array especially in the higher frequency bands, the array was designed to allow simultaneous operation of a low frequency receiver (<350 GHz) with a high frequency receiver (>330 GHz). The overlap region of 330-350 GHz was included to facilitate dual polarization measurements in the frequency range considered to offer the highest sensitivity for continuum observations with the array. So far, the array is equipped with working SIS receivers covering the frequency ranges 176-256 GHz, 260-350 GHz, and 600-700 GHz, and single frequency operation has been routine in the lower two frequency bands for the past year. More recently, with the completion of IF hardware required to make full use of the SMA cross-correlator, dual receiver operation became possible. We have since made a number of Galactic and extra-galactic astronomical observations in dual-band mode with the hopes of using the 230 GHz receiver as a phase reference to enable improved interferometry in the 650 GHz band. We will present the current antenna and receiver performance, some of the first interferometric images in the 650 GHz receiver band, and our initial attempts at phase transfer.

We present a spectrum of Sgr A* observed simultaneously on June 17, 2003 at wavelengths from 90 to 0.7 cm with the VLA. In the spectrum, we also include the measurements of Sgr A* observed on the same day with the GMRT at 49 cm, the SMA at 0.89 mm and the Keck II at 3.8 $\mu$m. The spectrum at the centimeter wavelengths suggests the presence of a break wavelength at 3.8 cm (8 GHz). The spectral index is alpha=0.43+-0.03 (propto nu(alpha}) at 3.8 cm and shorter wavelengths. The spectrum between 3.8 cm and 49 cm can be described by a power law with spectral index of alpha=0.10+-0.03$. We detected Sgr A* with 0.22+-0.06 Jy at 90 cm, suggesting a sharp decrease in flux density at the wavelengths longer than 49 cm. The best fit to the spectrum at the wavelengths longer than lambda_b appears to be consistent with free-free absorption by a screen of ionized gas with a turnover wavelength at nu(tau_ff}=1) = 100 cm (300 MHz). This turnover wavelength appears to be three times longer than that of 30 cm (1 GHz) as suggested by Davies et al. (1976) based on the observations in 1994 and 1995. Our analysis suggests that stellar winds from the massive stars near Sgr A* could modulate the flux density at the wavelengths longer than 30 cm (or frequencies below 1 GHz).

We present multi--frequency radio observational results of the quasar 3C~48. The observations were carried out with the Very Large Array (VLA) at five frequencies of 0.33, 1.5, 4.8, 8.4, and 22.5 GHz, and with the Multi--Element Radio Linked Interferometer Network (MERLIN) at the two frequencies of 1.6 and 5 GHz. The source shows a one--sided jet to the north within 1\arcsec, which then extends to the northeast and becomes diffuse. Two bright components (N2 and N3), containing most of the flux density are present in the northern jet. The spectral index of the two components is $\alpha_{N2}\sim-0.99\pm0.12$ and $\alpha_{N3}\sim-0.84\pm0.23$ ($S\propto\nu^{\alpha}$). Our images show the presence of an extended structure surrounding component N2, suggestive of strong interaction between the jet and the interstellar medium (ISM) of the host galaxy. A steep--spectrum component, labelled as S, located 0.25 ardsec southwest to the flat--spectrum component which could be the core of 3C 48, is detected at a significance of $>15\sigma$. Both the location and the steepness of the spectrum of component S suggest the presence of a counter--jet in 3C 48.

We present images of NRAO530 observed with the EVN (VLBI) at 5 GHz, the MERLIN at 1 .6 and 5 GHz, and the VLA at 5 and 8 GHz showing the complex morphology on scales from pc to kpc. The VLBI image shows a core-jet structure indicating a somehow oscillation trajectory on a scale of 30 mas, north to the strongest compact component (core). A core-jet structure extended to several hundreds mas at about P.A. -50 deg and a distant component located 11 arcsec west to the core are detected in both the MERLIN and the VLA observations. An arched structure of significant emission between the core and the distant component is also revealed in both the MERLIN image at 1.6 cm and the VLA images at 8.4 and 5 GHz. The core component shows a flat spectrum with alpha = -0.02 (S proportional to the frequency power -alpha) while alpha = 0.8 for the distant component. The steep spectrum of the distant component and the detection of the arched emission suggests that the western distant component is a lobe or a hot-spot powered by the nucleus of NRAO530. A patch of diffuse emission, 12 arcsec nearly east (P.A. = 70 deg) to the core component, is also observed with the VLA at 5 GHz, suggesting a presence of a counter lobe in the source.

We have imaged a region of about 5' extent surrounding Sgr A* in the HI 21 cm-line absorption using the Very Large Array. A Gaussian decomposition of the optical depth spectra at positions within about 2' (approx. 5 pc at 8.5 kpc) of Sgr A* detects a wide line underlying the many narrow absorption lines. The wide line has a mean peak optical depth of 0.32 +/- 0.12 centered at a mean velocity of V(lsr) = -4 +/- 15 km/s. The mean full width at half maximum is 119 +/- 42 km/s. Such a wide line is absent in the spectra at positions beyond about 2' from Sgr A*. The position-velocity diagrams in optical depth reveal that the wide line originates in various components of the circumnuclear disk (radius approx. 1.3') surrounding Sgr A*. These components contribute to the optical depth of the wide line in different velocity ranges. The position-velocity diagrams do not reveal any diffuse feature which could be attributed to a large number of HI clouds along the line of sight to Sgr A*. Consequently, the wide line has no implications either to a global population of shocked HI clouds in the Galaxy or to the energetics of the interstellar medium as was earlier thought.

We present high angular resolution (0.6'') observations made with the VLA of the radio continuum at 8.3 and 43 GHz as well as H92 Alpha and H53 Alpha radio recombination lines from the nearby (3 Mpc) starburst galaxy M82. In the continuum we report 19 newly identified sources at 8.3 GHz and 5 at 43 GHz that were unknown previously. The spatial distribution of the H92 Alpha line is inhomogeneous; we identify 27 features. The line and continuum emission are modeled using a collection of HII regions at different distances from the nucleus assuming a single-density component and two-density components. The high-density component has a density of 4 X 10^{4} cm^{-3}. However, the bulk of the ionization is in regions with densities which are typically a factor 10 lower. The gas kinematics, using the H92 Alpha line, confirms the presence of steep velocity gradient (26 km s^{-1} arcsec^{-1}) in the nuclear region. As this steep gradient is observed not only on the major axis but also at large distances along a band of PA of 150 degrees, the interpretation in terms of x2 orbits elongated along the minor axis of the bar, which would be observed at an angle close to the inclination of the main disk, seems inadequate. Ad-hoc radial motions must be introduced to reproduce the pattern of the velocity field. Different families of orbits are indicated as we detect a signature in the kinematics at the transition between the two plateaus observed in the NIR light distribution. The H92 Alpha line also reveals the base of the outflow where the injection towards the halo on the Northern side occurs. The kinematical pattern suggests a connection between the gas flowing in the plane of M82 towards the center; this behavior most likely originates due to the presence of a bar and the outflow out of the plane.

We have calculated proper motions and re-evaluated the membership probabilities of 810 stars in the area of two NGC objects, NGC 1817 and NGC 1807. We have obtained absolute proper motions from 25 plates in the reference system of the Tycho-2 Catalogue. The plates have a maximum epoch difference of 81 years; and they were taken with the double astrograph at Zo-Se station of Shanghai Observatory, which has an aperture of 40 cm and a plate scale of 30 arcsec/mm. The average proper motion precision is 1.55 mas/yr. These proper motions are used to determine the membership probabilities of stars in the region, based on there being only one very extended physical cluster: NGC 1817. With that aim, we have applied and compared parametric and non-parametric approaches to cluster/field segregation. We have obtained a list of 169 probable member stars.

Using the partially completed Submillimeter Array with five antennas, we have observed the CO J=2-1 and 3-2 emission from the envelope surrounding the carbon star V Hya. The high-angular resolution (2"-4") maps show that V Hya is powering a bipolar molecular jet having an extreme velocity of 70-185 km/s. The axis of this high velocity jet is perpendicular to the major axis of the flattened disk-like envelope, which is expanding with a velocity of ~16 km/s.There is a third kinematic component, a medium-velocity wind having a deprojected velocity of 40-120 km/s moving along the disk plane. Both the high velocity jet and the medium velocity wind have a dynamical time scale of a few hundred years. The flattened structure and the collimated jet observed in V Hya suggests that the formation of asymmetrical structure proceeds while the central star is still in the AGB phase

Submillimeter Array (SMA) observations of the high-mass star-forming region IRAS 18089-1732 in the 1 mm and 850 $\mu$m band with 1 GHz bandwidth reveal a wealth of information. We present the observations of 34 lines from 16 different molecular species. Most molecular line maps show significant contributions from the outflow, and only few molecules are confined to the inner core. We present and discuss the molecular line observations and outline the unique capabilities of the SMA for future imaging line surveys at high spatial resolution.

We present initial results of an ongoing search for interferometric calibrators at submillimeter (sub-mm) wavelengths with the Submillimeter Array (SMA). Powerful radio galaxies are commonly used as calibrators at centimeter and millimeter wavelengths, but many are not strong enough to serve as calibrators at sub-mm wavelengths because of their rapidly declining flux densities toward shorter wavelengths. The inability to find a calibrator close to the target source may limit or even prevent us from imaging many interesting sources at sub-mm wavelengths. Here, we investigate whether high-mass protostellar objects and ultracompact HII regions can serve as useful calibrators for the SMA. The dust emission associated with these objects makes them among the brightest sub-mm sources in the sky. Our observations at 0.85 mm (345 GHz) with an angular resolution of ~3" reveal that although a large fraction of the dust emission originates from an extended ``halo'' component, a compact unresolved component often remains that when sufficiently strong may serve as a useful calibrator. These observations also provide a first glimpse at the small-scale distribution of dust around ultracompact HII regions and high-mass protostellar objects at sub-mm wavelengths. We discuss the origin of the core-halo structure seen in many sources, and conclude with suggestions for future searches for calibrators with the SMA.

We have detected the intrinsic size of Sagittarius A*, the Galactic Center radio source associated with a supermassive black hole, showing that the short-wavelength radio emission arises from very near the event horizon of the black hole. Radio observations with the Very Long Baseline Array show that the source has a size of $24 \pm 2$ Schwarzschild radii at 7 mm wavelength. In one of eight 7-mm epochs we also detect an increase in the intrinsic size of $60^{+25}_{-17}$%. These observations place a lower limit to the mass density of Sgr A* of $1.4\times10^{4}$ solar masses per cubic astronomical unit.

We have mapped the ultracompact HII region, G5.89-0.39, and its molecular surroundings with the Submillimeter Array at 2".8 x 1".8 angular resolution in 1.3 mm continuum, SiO(5-4), and eight other molecular lines. We have resolved for the first time the highly energetic molecular outflow in this region. At this resolution, the outflow is definitely bipolar and appears to originate in a 1.3 mm continuum source. The continuum source peaks in the center of the HII region. The axis of the outflow lines up with a recently discovered O5V star.

We present the results of a 3.3-year project to monitor the flux density of Sagittarius A* at 2.0, 1.3, and 0.7 cm with the VLA. The fully calibrated light curves for Sgr A* at all three wavelengths are presented. Typical errors in the flux density are 6.1%, 6.2%, and 9.2% at 2.0, 1.3, and 0.7 cm, respectively. There is preliminary evidence for a bimodal distribution of flux densities, which may indicate the existence of two distinct states of accretion onto the supermassive black hole. At 1.3 and 0.7 cm, there is a tail in the distribution towards high flux densities. Significant variability is detected at all three wavelengths, with the largest amplitude variations occurring at 0.7 cm. The rms deviation of the flux density of Sgr A* is 0.13, 0.16, and 0.21 Jy at 2.0, 1.3, and 0.7 cm, respectively. During much of this monitoring campaign, Sgr A* appeared to be relatively quiescent compared to results from previous campaigns. At no point during the monitoring campaign did the flux density of Sgr A* more than double its mean value. The mean spectral index of Sgr A* is alpha=0.20+/-0.01, with a standard deviation of 0.14. The spectral index appears to depend linearly on the observed flux density at 0.7 cm with a steeper index observed during outbursts. This correlation is consistent with the expectation for outbursts that are self-absorbed at wavelengths of 0.7 cm or longer and inconsistent with the effects of simple models for interstellar scintillation. Much of the variability of Sgr A*, including possible time lags between flux density changes at the different wavelengths, appears to occur on time scales less than the time resolution of our observations (8 days). Future observations should focus on the evolution of the flux density on these time scales.

SMA observations of the massive star-forming region IRAS 18089-1732 in the 1mm and 850mu band reveal outflow and disk signatures in different molecular lines. The SiO(5--4) data show a collimated outflow in the northern direction. In contrast, the HCOOCH3(20--19) line, which traces high-density gas, is confined to the very center of the region and shows a velocity gradient across the core. The HCOOCH3 velocity gradient is not exactly perpendicular to the outflow axis but between an assumed disk plane and the outflow axis. We interpret these HCOOCH3 features as originating from a rotating disk that is influenced by the outflow and infall. Based on the (sub-)mm continuum emission, the mass of the central core is estimated to be around 38M_sun. The dynamical mass derived from the HCOOCH3 data is 22Msun, of about the same order as the core mass. Thus, the mass of the protostar/disk/envelope system is dominated by its disk and envelope. The two frequency continuum data of the core indicate a low dust opacity index beta ~ 1.2 in the outer part, decreasing to beta ~ 0.5 on shorter spatial scales.

A large radio outburst from Sgr A* was observed during the VLA weekly monitoring program at 2 cm, 1.3 cm and 7 mm, nearly coincident with the brightest X-ray flare detected to date with the XMM-Newton X-ray Observatory on 2002-10-03. The flux density of 1.9$\pm$0.2 Jy measured at 7 mm exceeds the mean value (1.00$\pm$0.01 Jy) by a factor of $\sim 2$, one of the two highest increases observed during the past three years (June 2000-October 2003), while less significant increases in flux densities were observed at 1.3 cm and 2 cm. The radio observation started 13.5 hrs after the onset of the X-ray flare (which had occurred over a 45 min duration) and continued for 1.3 hrs. During the observation, there was no significant ($<3\sigma$) change in the radio flux densities at all the three wavelengths, indicating that the radio outburst varied on a timescale of $>1$hr. A spectral index of $\alpha=2.4^{+0.3}_{-0.6}$ (${\rm S \propto\nu^\alpha}$) was derived for the outburst component, consistent with an optically thick nonthermal synchrotron source. These results suggest that energetic electrons responsible for the radio outburst might be produced via a process associated with the X-ray flare, then transported to large radii, producing the observed radio outburst. The observation is the first evidence for a correlated variation in the radio and X-ray emissions from Sgr A*.

We report the results from recent observations of Sgr A* at short-/sub-millimeter wavelengths made with the partially finished Sub-Millimeter Array (SMA) on Mauna Kea. A total of 25 epochs of observations were carried out over the past 15 months in 2001 March to 2002 May. Noticeable variations in flux density at 1.3 mm were observed showing three ``flares''. The SMA observations suggest that Sgr A* highly increases towards submillimeter wavelengths during a flare suggesting the presence of a break wavelength in spectral index around 3 mm. A cross-correlation of the SMA data at 1 mm with the VLA data at 1 cm show a global delay of $t_{delay} > 3d$, suggesting that sub-millimeter wavelengths tend to peak first. Only marginal day-to-day variations in flux density (2-3 $\sigma$) have been detected at 1.3 mm. No significant flares on a short time scale ($\sim1$ hr) have been observed at 1.3 mm. We also failed to detect significant periodic signals at a level of 5% (3$\sigma$) from Sgr A* in a periodic searching window ranging from 10 min to 2.5 hr. The flares observed at the wavelengths between short-centimeter and sub-millimeter might be a result of collective mass ejections associated with X-ray flares that originate from the inner region of the accretion disk near the supermassive black hole.

The Galactic center provides a unique astrophysical laboratory for us to study various astrophysical processes. In this paper, we review and outline the latest results from observations of Sgr~A$^*$ in terms of source structure and variations in flux density. Sgr~A$^*$ phenomenon represents a typical case of low radiative efficiency accretion flow surrounding a supermassive black hole in low luminosity AGNs. Many pending astrophysical problems found from observations of Sgr A$^*$ have challenged the existing astrophysical theories. Current theoretical models of Sgr A$^*$ are also reviewed.

Absolute proper motions, their corresponding errors and membership probabilities of 501 stars in the intermediate-age open cluster NGC2548 region are determined from MAMA measurements of 10 photographic plates. The plates have the maximum epoch difference of 82 years and they were taken with the double astrograph at Zo-Se station of Shanghai Observatory, which has an aperture of 40 cm and a plate scale of 30"/mm. The average proper motion precision is 1.18 mas/yr. These proper motions are used to determine the membership probabilities of stars in the region. The number of stars with membership probabilities higher than 0.7 is 165.

At radio wavelengths, images of the compact radio source Sagittarius A* (Sgr A*) in the Galactic Center are scatter broadened with a lambda^2 dependence due to an intervening ionized medium. We present VLBI observations of Sgr A* at 86 GHz using a six station array including the VLBA antennas at Pie Town, Fort Davis and Los Alamos, the 12m antenna at Kitt Peak and the millimeter arrays at Hat Creek and Owens Valley. To avoid systematic errors due to imperfect antenna calibration, the data were modeled using interferometric closure information. The data are best modeled by a circular Gaussian brightness distribution of FWHM 0.18 +- 0.02 mas. The data are also shown to be consistent with an elliptical model corresponding to the scattering of a point source. The source structure in the N-S direction, which is less well determined than in the E-W direction due to the limited N-S (u,v) coverage of the array, is constrained to be less than 0.27 mas by these measurements. These results are consistent with extrapolations of intrinsic structure estimates obtained with VLBI at 7mm wavelength assuming the intrinsic size of Sgr A* has a greater dependence than lambda^0.9 with wavelength.

We report the presence of a 106-day cycle in the radio variability of Sgr A* based on an analysis of data observed with the Very Large Array (VLA) over the past 20 years. The pulsed signal is most clearly seen at 1.3 cm with a ratio of cycle frequency to frequency width f/Delta_f= 2.2+/-0.3. The periodic signal is also clearly observed at 2 cm. At 3.6 cm the detection of a periodic signal is marginal. No significant periodicity is detected at both 6 and 20 cm. Since the sampling function is irregular we performed a number of tests to insure that the observed periodicity is not the result of noise. Similar results were found for a maximum entropy method and periodogram with CLEAN method. The probability of false detection for several different noise distributions is less than 5% based on Monte Carlo tests. The radio properties of the pulsed component at 1.3 cm are spectral index alpha ~ 1.0+/- 0.1 (for S nu^alpha), amplitude Delta S=0.42 +/- 0.04 Jy and characteristic time scale Delta t_FWHM ~ 25 +/- 5 days. The lack of VLBI detection of a secondary component suggests that the variability occurs within Sgr A* on a scale of ~5 AU, suggesting an instability of the accretion disk.

We present observations of radio recombination lines from the starburst galaxy Arp 220 at 1.4, 8.1, 84, 96 and 207 GHz (sensitive upper limit for the 1.4 GHz line and firm detections at the other frequencies), and the radio continuum spectrum between 330 MHz and 207 GHz. We show that a model with three components of ionized gas with different densities and area covering factors can consistently explain both RRL and continuum data. The total mass of ionized gas in the three components is 3 x 10^7 M_sun requiring 3 x 10^5 O5 stars with a total Lyman continuum production rate (NLyc) of 1.3 x 10^{55} photons /s. These values imply a dust extinction A_V ~ 45 magnitudes and an SFR of ~240 M_sun/yr. The NLyc of ~3% associated with the high density HII regions implies similar SFR at recent epochs. The data is also consistent with multiple starbursts of very high SFR and short durations. The derived value of 24 for the IR-excess favours a starburst rather than an AGN as the origin of the observed FIR luminosity. (the abstract has been abridged)

This paper presents some statistical correlations of 72 northern spiral galaxies. The results show that early-type spirals that are brighter, and thicker, and the axis ratios ($H_{z}/H_r$) of the disk tend to be smaller along the Hubble sequence. We also find that $H_{z}/H_r$ correlates strongly with the galaxy's color, and early-type spirals have larger values of $H_{z}/H_r$. The inclinations obtained by fitting the pattern of a spiral structure with a logarithmic spiral form are nearly the same as those obtained by using the formulas of Aaronson et al. (1980). Finally, the mean measured pitch angles for the different Hubble sequences in the Third Catalogue of Bright Galaxies by de Vaucouleurs et al. (1991) are derived.

Photographic photometries of 789 stars in the region of the open clusters NGC1750 and NGC1758 are derived from a set of 8 photometric plates. According to the astronomical data including proper motions, positions and membership probabilities of individual stars, several astrophysical parameters for these two clusters are determined, such as HR diagrams, ages, distances, luminosity functions, masses and kinematics, etc. It is found that the distances are ($525\pm 48$pc) and ($794\pm 73$pc) with ages of $1.5\times10^8$yr and $6.3\times10^8$yr for NGC1750 and NGC1758 respectively. Both clusters show no significant mass segregation effects because of their relatively young ages. The analysis of the proper motions imply that their velocity distributions are isotropy.

Synchrotron radiation from active galactic nuclei (AGN) is often highly polarized. We present a search for linear polarization with the Very Large Array (VLA) at 4.8 GHz and 8.4 GHz from the nearest AGN, Sagittarius A*. As a part of this study we used spectro-polarimetric data that were sensitive to a rotation measure (RM) as large as 3.5 x 10^6 rad m^-2 at 4.8 GHz and 1.5 x 10^7 rad m^-2 at 8.4 GHz. The upper limit to the linear polarization of Sgr A* over a broad range of RM is 0.2% at both frequencies. We also present continuum observations with the VLA at 4.8 GHz which give an upper limit of 0.1% for RMs less than 10^4 rad m^-2. We conclude that depolarization is unlikely to occur in the Galacter Center scattering medium. However, it is possible for depolarization to occur in the accretion region of Sgr A* if the outer scale of turbulence is small enough. We also consider the implications of a very low intrinsic polarization for Sgr A*.

By means of near-simultaneous multi-wavelength VLBA measurements, we determine for the first time the intrinsic size of Sgr A* to be 3.6 AU by <1 AU with the major axis oriented essentially north-south. Contrary to previous expectation that the intrinsic structure of Sgr A* is observable only at wavelength shorter than 1 mm, we can discern the intrinsic source size at 7~mm because (1) the scattering size along the minor axis is half that along the major axis, and (2) the near simultaneous multi-wavelength mapping of Sgr A* makes it possible to extrapolate precisely the minor axis scattering angle at 7~mm. The intrinsic size and shape place direct constraints on the various theoretical models for Sgr A*.

Recent proper motion studies of stars at the very center of the Galaxy strongly suggest that Sagittarius (Sgr) A*, the compact nonthermal radio source at the Galactic Center, is a 2.5 million solar mass black hole. By means of near-simultaneous multi-wavelength Very Long Baseline Array measurements, we determine for the first time the intrinsic size and shape of Sgr A* to be 72 Rsc by < 20 Rsc, with the major axis oriented essentially north-south, where Rsc (= 7.5 x 10^{11} cm) is the Schwarzschild radius for a 2.5 million solar mass black hole. Contrary to previous expectation that the intrinsic structure of Sgr A* is observable only at wavelengths shorter than 1 mm, we can discern the intrinsic source size at 7 mm because (1) the scattering size along the minor axis is half that along the major axis, and (2) the near simultaneous multi-wavelength mapping of Sgr A* with the same interferometer makes it possible to extrapolate precisely the minor axis scattering angle at 7 mm. The intrinsic size and shape place direct constraints on the various emission models for Sgr A*. In particular, the advection dominated accretion flow model may have to incorporate a radio jet in order to account for the structure of Sgr A*.

We report results from a multiwavelength campaign to measure the simultaneous spectrum of the super-massive black hole candidate Sgr A* in the Galactic Center from cm to mm-wavelengths using the VLA, BIMA, the Nobeyama 45m, and the IRAM 30m telescopes. The observations confirm that the previously detected mm-excess is an intrinsic feature of the spectrum of Sgr A*. The excess can be interpreted as due to the presence of an ultra-compact component of relativistic plasma with a size of a few Schwarzschild radii near the black hole. If so, Sgr A* might offer the unique possibility to image the putative black hole against the background of this component with future mm-VLBI experiments.

A simplified frequency-integrated radiative transfer equation is solved to study Compton scatterings in the corona of the disk by using numerical iterating method. We find that the vertical thickness of the corona cannot be used as the typical length to measure the optical depth of the corona. A semi-analytical approach is proposed to calculate the energy dissipations in the corona of the disk. We demonstrate that our approach can reproduce the numerical solutions to an accuracy of <2 %.

We report first results from a multiwavelength campaign to measure the simultaneous spectrum of Sgr A* from cm to mm wavelengths. The observations confirm that the previously detected submm-excess is not due to variability; the presence of an ultracompact component with a size of a few Schwarzschild radii is inferred. In a VLA survey of LINER galaxies, we found Sgr A*-like nuclei in one quarter of the galaxies searched, suggesting a link between those low-power AGN and the Galactic Center.