The authors of this white paper met on 16-17 January 2020 at the New Jersey Institute of Technology, Newark, NJ, for a 2-day workshop that brought together a group of heliophysicists, data providers, expert modelers, and computer/data scientists. Their objective was to discuss critical developments and prospects of the application of machine and/or deep learning techniques for data analysis, modeling and forecasting in Heliophysics, and to shape a strategy for further developments in the field. The workshop combined a set of plenary sessions featuring invited introductory talks interleaved with a set of open discussion sessions. The outcome of the discussion is encapsulated in this white paper that also features a top-level list of recommendations agreed by participants.

We consider the high spin expansion for the null geodesics in the Kerr spacetime. We expand the null geodesic equation successively to higher orders in deviation from extremity. Via the method of matched asymptotic expansion, the radial integrals are obtained analytically. It turns out that the analytic expressions are very sensitive to the value of the shifted Carter constant $q$. We show that for a large $q$, the analytic expressions can be used to study observational electromagnetic signatures for astrophysical black holes like M87*. However, for a small $q$, the high spin expansion method can only be applied to (near-) extreme black holes.

We present a sample of $\sim$ 140,000 primary red clump (RC) stars of spectral signal-to-noise ratios higher than 20 from the LAMOST Galactic spectroscopic surveys, selected based on their positions in the metallicity-dependent effective temperature--surface gravity and color--metallicity diagrams, supervised by high-quality $Kepler$ asteroseismology data. The stellar masses and ages of those stars are further determined from the LAMOST spectra, using the Kernel Principal Component Analysis method, trained with thousands of RCs in the LAMOST-$Kepler$ fields with accurate asteroseismic mass measurements. The purity and completeness of our primary RC sample are generally higher than 80 per cent. For the mass and age, a variety of tests show typical uncertainties of 15 and 30 per cent, respectively. Using over ten thousand primary RCs with accurate distance measurements from the parallaxes of Gaia DR2, we re-calibrate the $K_{\rm s}$ absolute magnitudes of primary RCs by, for the first time, considering both the metallicity and age dependencies. With the the new calibration, distances are derived for all the primary RCs, with a typical uncertainty of 5--10 per cent, even better than the values yielded by the Gaia parallax measurements for stars beyond 3--4 kpc. The sample covers a significant volume of the Galactic disk of $4 \leq R \leq 16$ kpc, $|Z| \leq 5$ kpc, and $-20 \leq \phi \leq 50^{\circ}$. Stellar atmospheric parameters, line-of-sight velocities and elemental abundances derived from the LAMOST spectra and proper motions of Gaia DR2 are also provided for the sample stars. Finally, the selection function of the sample is carefully evaluated in the color-magnitude plane for different sky areas. The sample is publicly available.

RV variable stars are important in astrophysics. The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) spectroscopic survey has provided ~ 6.5 million stellar spectra in its Data Release 4 (DR4). During the survey, ~ 4.7 million unique sources were targeted and ~ 1 million stars observed repeatedly. The probabilities of stars being RV variables are estimated by comparing the observed radial velocity variations with the simulated ones. We build a catalog of 80,702 RV variable candidates with probability greater than 0.60 by analyzing the duplicate-observed multi-epoch sources covered by the LAMOST DR4. Simulations and cross-identifications show that the purity of the catalog is higher than 80%. The catalog consists of 77% binary systems and 7% pulsating stars as well as 16% pollution by single stars. 3,138 RV variables are classified through cross-identifications with published results in literatures. By using the 3,138 sources common to both LAMOST and a collection of published RV variable catalogs we are able to analyze LAMOST's RV variable detection rate. The efficiency of the method adopted in this work relies not only on the sampling frequency of observations but also periods and amplitudes of RV variables. With the progress of LAMOST, Gaia and other surveys, more and more RV variables would will be confirmed and classified. This catalog is valuable for other large-scale surveys, especially for RV variable searches. The catalog will be released according to the LAMOST Data Policy via http://dr4.lamost.org.

We revisit the diagonal ridge feature (diagonal distributions in the $R, v_{\phi}$ plane) found in $Gaia$ and present timing analysis for it between Galactocentric distances of $R=7.5$ and 12 \,kpc, using Main-Sequence-Turn-Off and OB stars selected from the LAMOST Galactic spectroscopic surveys. We recover the ridge pattern in the $R$-$v_{\phi}$ plane color coded by stellar number density and mean radial velocity and find this feature is presented from very young (OB stars, few hundred \,Myr) to very old populations ($\tau > $ 9 \,Gyr), accompanying with the north-south asymmetrical ridges. One of the ridge patterns, with constant angular momentum per unit mass, shows clear variations with different age populations compared. However, the remaining two ones are very stable, clearly implying there might have two kinds of ridge patterns with different dynamical origins and evolution. For the vertical velocity and height, there are no ridge features corresponding to the in-plane Galactoseismogy features, but the metallicity and [$\alpha/Fe$] have some weak signals. This first temporal chemo-kinematical analysis for ridge strongly supports that the ridge pattern might be originated from the phase mixing of disrupting spiral arms implying different ridges have different physical processes and no coupling phenomenon, although other possible scenarios such as Sgr perturbations could not be ruled out.

We report microwave spectral imaging observations of an erupting magnetic flux rope during the early impulsive phase of the X8.2-class limb flare on 2017 September 10, obtained by the Expanded Owens Valley Solar Array. A few days prior to the eruption, when viewed against the disk, the flux rope appeared as a reverse S-shaped dark filament along the magnetic polarity inversion line. During the eruption, the rope exhibited a "hot channel" structure in extreme ultraviolet and soft X-ray passbands sensitive to ~10 MK plasma. The central portion of the flux rope was nearly aligned with the line of sight, which quickly developed into a teardrop-shaped dark cavity during the early phase of the eruption. A long and thin plasma sheet formed below the cavity, interpreted as the reconnection current sheet viewed edge-on. A nonthermal microwave source was present at the location of the central current sheet, which extended upward encompassing the dark cavity. A pair of nonthermal microwave sources were observed for several minutes on both sides of the main flaring region. They shared a similar temporal behavior and spectral property to the central microwave source below the cavity, interpreted as the conjugate footpoints of the erupting flux rope. These observations are broadly consistent with the magnetic topology and the associated energy release scenario suggested in the three-dimensional standard model for eruptive solar flares. In particular, our detection of nonthermal emission at conjugate flux rope footpoints provides solid evidence of particle transport along an erupting magnetic flux rope.

In the standard model of solar flares, a large-scale reconnection current sheet is postulated as the central engine for powering the flare energy release and accelerating particles. However, where and how the energy release and particle acceleration occur remain unclear due to the lack of measurements for the magnetic properties of the current sheet. Here we report the measurement of spatially-resolved magnetic field and flare-accelerated relativistic electrons along a current-sheet feature in a solar flare. The measured magnetic field profile shows a local maximum where the reconnecting field lines of opposite polarities closely approach each other, known as the reconnection X point. The measurements also reveal a local minimum near the bottom of the current sheet above the flare loop-top, referred to as a "magnetic bottle". This spatial structure agrees with theoretical predictions and numerical modeling results. A strong reconnection electric field of ~4000 V/m is inferred near the X point. This location, however, shows a local depletion of microwave-emitting relativistic electrons. These electrons concentrate instead at or near the magnetic bottle structure, where more than 99% of them reside at each instant. Our observations suggest that the loop-top magnetic bottle is likely the primary site for accelerating and/or confining the relativistic electrons.

Distance measurements of supernova remnants (SNRs) are essential and important. Accurate estimates of physical size, dust masses, and some other properties of SNRs depend critically on accurate distance measurements. However, the determination of SNR distances is still a tough task. Red clump stars (RCs) have a long history been used as standard candles. In this work, we take RCs as tracers to determine the distances to a large group of SNRs in the inner disk. We first select RC stars based on the near-infrared (IR) color-magnitude diagram (CMD). Then, the distance to and extinction of RC stars are calculated. To extend the measurable range of distance, we combine near-IR photometric data from the 2MASS survey with the deeper UKIDSS and VVV surveys. With the help of the Gaia parallaxes, we also remove contaminants including dwarfs and giants. Because an SN explosion compresses the surrounding interstellar medium, the SNR region would become denser and exhibit higher extinction than the surroundings. The distance of a SNR is then recognized by the position where the extinction and its gradient is higher than that of the ambient medium. A total of 63 SNRs' distances in the Galactic inner disk are determined and divided into three Levels A, B, and C with decreasing reliability. The distances to 43 SNRs are well determined with reliability A or B. The diameters and dust masses of SNRs are estimated with the obtained distance and extinction.

By combining LAMOST DR4 and Gaia DR2 common red clump stars with age and proper motion, we analyze the amplitude evolution of the stellar warp independently of any assumption with a simple model. The greatest height of the warp disk increases with Galactocentric distance in different populations and it is dependent on the age: the younger stellar populations exhibit stronger warp features than the old ones, accompanied with the warp amplitude $\gamma ({\rm age})$ decreasing with age and its first derivative $\dot{\gamma} ({\rm age})$ is different from zero. The azimuth of line of nodes $\phi _w$ is stable at $-$5 degree without clear time evolution, which perfectly confirms some of previous works. All these self-consistent evidences support that our Galactic warp should most likely be a long lived, but non-steady structure and not a transient one, which is supporting the warp is originated from gas infall onto the disk or other hypotheses that suppose that the warp mainly affects to the gas, and consequently younger populations tracing the gas are stronger than the older ones. In other words, the Galactic warp is induced by the non-gravitational interaction over the disk models.

Nonthermal electrons accelerated in solar flares produce electromagnetic emission in two distinct, highly complementary domains - hard X-rays (HXRs) and microwaves (MWs). This paper reports MW imaging spectroscopy observations from the Expanded Owens Valley Solar Array of an M1.2 flare that occurred on 2017 September 9, from which we deduce evolving coronal parameter maps. We analyze these data jointly with the complementary Reuven Ramaty High-Energy Solar Spectroscopic Imager HXR data to reveal the spatially-resolved evolution of the nonthermal electrons in the flaring volume. We find that the high-energy portion of the nonthermal electron distribution, responsible for the MW emission, displays a much more prominent evolution (in the form of strong spectral hardening) than the low-energy portion, responsible for the HXR emission. We show that the revealed trends are consistent with a single electron population evolving according to a simplified trap-plus-precipitation model with sustained injection/acceleration of nonthermal electrons, which produces a double-powerlaw with steadily increasing break energy.

We installed two sets of Astronomical Site Monitoring System(ASMS) at Lijiang Observatory(GMG), for the running of the 2.4-meter Lijiang optical telescope(LJT) and the 1.6-meter Multi-channel Photometric Survey Telescope (Mephisto). The Mephistro is under construction. ASMS has been running on robotic mode since 2017. The core instruments: Cloud Sensor, All-Sky Camera and Autonomous-DIMM that are developed by our group, together with the commercial Meteorological Station and Sky Quality Meter, are combined into the astronomical optical site monitoring system. The new Cloud Sensor's Cloud-Clear Relationship is presented for the first time, which is used to calculate the All-Sky cloud cover. We designed the Autonomous-DIMM located on a tower, with the same height as LJT. The seeing data have been observed for a full year. ASMS's data for the year 2019 are also analysed in detail, which are valuable to observers.

Strong quasar-galaxy lensing provides a powerful tool to probe the inter-stellar medium (ISM) of the lens galaxy using radiation from the background quasar. Using the Cosmic Origins Spectrograph (COS) on board the Hubble Space Telescope, we study the cold ISM properties of the lens galaxy in B1152+199 at a redshift of z=0.4377. Since existing optical extinction and X-ray absorption measurements of the lens have revealed a large amount of cold ISM, we expected to detect a damped Lya absorption (DLA) system in the near ultraviolet spectrum; however, our upper limit on the HI column density is several orders of magnitude below the expectation. We also marginally detect OI and CII absorption lines associated with the lens galaxy in the COS spectrum. Thus, the lens galaxy is identified as a ghostly DLA system, and further investigations of these ghostly DLA systems would be important to characterize the biases of using DLAs to probe the matter density of the universe. Although preliminary, the most likely explanation of the non-detection of the DLA is because of the Lya emission of the lens galaxy that fills in the absorption trough, with a Lya luminosity of 4e42 erg/s.

We explore the fundamental relations governing the radial and vertical velocity dispersions of stars in the Milky Way, from combined studies of complementary surveys including GALAH, LAMOST, APOGEE, the NASA $Kepler$ and K2 missions, and $Gaia$ DR2. We find that different stellar samples, even though they target different tracer populations and employ a variety of age estimation techniques, follow the same set of fundamental relations. We provide the clearest evidence to date that, in addition to the well-known dependence on stellar age, the velocity dispersions of stars depend on orbital angular momentum $L_z$, metallicity and height above the plane $|z|$, and are well described by a multiplicatively separable functional form. The dispersions have a power-law dependence on age with exponents of 0.441$\pm 0.007$ and 0.251$\pm 0.006$ for $\sigma_z$ and $\sigma_R$ respectively, and the power law is valid even for the oldest stars. For the solar neighborhood stars, the apparent break in the power law for older stars, as seen in previous studies, is due to the anti-correlation of $L_z$ with age. The dispersions decrease with increasing $L_z$ until we reach the Sun's orbital angular momentum, after which $\sigma_z$ increases (implying flaring in the outer disc) while $\sigma_R$ flattens. The dispersions increase with decreasing metallicity, suggesting that the dispersions increase with birth radius. The dispersions also increase linearly with $|z|$. The same set of relations that work in the solar neighborhood also work for stars between $3<R/{\rm kpc}<20$. Finally, the high-[$\alpha$/Fe] stars follow the same relations as the low-[$\alpha$/Fe] stars.

We report the detection of emission from a non-thermal electron distribution in a small solar microflare (GOES class A5.7) observed by the Nuclear Spectroscopic Telescope Array (NuSTAR), with supporting observation by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). The flaring plasma is well accounted for by a thick-target model of accelerated electrons collisionally thermalizing within the loop, akin to the "coronal thick target" behavior occasionally observed in larger flares. This is the first positive detection of non-thermal hard X-rays from the Sun using a direct imager (as opposed to indirectly imaging instruments). The accelerated electron distribution has a spectral index of 6.3 +/- 0.7, extends down to at least 6.5 keV, and deposits energy at a rate of ~2x1027 erg/s, heating the flare loop to at least 10 MK. The existence of dominant non-thermal emission in X-rays down to <5 keV means that RHESSI emission is almost entirely non-thermal, contrary to what is usually assumed in RHESSI spectroscopy. The ratio of non-thermal to thermal energies is similar to that of large flares, in contrast to what has been found in previous studies of small RHESSI flares. We suggest that a coronal thick target may be a common property of many small microflares based on the average electron energy and collisional mean free path. Future observations of this kind will enable understanding of how flare particle acceleration changes across energy scales, and will aid the push toward the observational regime of nanoflares, which are a possible source of significant coronal heating.

The SunPy Project developed a 13-question survey to understand the software and hardware usage of the solar physics community. 364 members of the solar physics community, across 35 countries, responded to our survey. We found that 99$\pm$0.5% of respondents use software in their research and 66% use the Python scientific software stack. Students are twice as likely as faculty, staff scientists, and researchers to use Python rather than Interactive Data Language (IDL). In this respect, the astrophysics and solar physics communities differ widely: 78% of solar physics faculty, staff scientists, and researchers in our sample uses IDL, compared with 44% of astrophysics faculty and scientists sampled by Momcheva and Tollerud (2015). 63$\pm$4% of respondents have not taken any computer-science courses at an undergraduate or graduate level. We also found that most respondents utilize consumer hardware to run software for solar-physics research. Although 82% of respondents work with data from space-based or ground-based missions, some of which (e.g. the Solar Dynamics Observatory and Daniel K. Inouye Solar Telescope) produce terabytes of data a day, 14% use a regional or national cluster, 5% use a commercial cloud provider, and 29% use exclusively a laptop or desktop. Finally, we found that 73$\pm$4% of respondents cite scientific software in their research, although only 42$\pm$3% do so routinely.

Drifting pulsation structures (DPSs) are important radio fine structures usually observed at the beginning of eruptive solar flares. It has been suggested that DPSs carry important information on the energy release processes in solar flares. We study DPS observed in an X8.2-class flare on 2017 September 10 in the context of spatial and spectral diagnostics provided by microwave, EUV, and X-ray observations. We describe DPS and its sub-structures that were observed for the first time. We use a new wavelet technique to reveal characteristic periods in DPS and their frequency bands. Comparing the periods of pulsations found in this DPS with those in previous DPSs we found new very short periods in the 0.09-0.15 s range. We present EOVSA images and spectra of microwave sources observed during the DPS. This DPS at its very beginning has pulsations in two frequency bands (1000-1300 MHz and 1600--1800 MHz) which are interconnected by fast drifting bursts. We show that these double-band pulsations started just at the moment when the ejected filament splits apart in a tearing motion at the location where a signature of the flare current sheetlater appeared. Using the standard flare model and previous observations of DPSs, we interpret these double-band pulsations as a radio signature of superthermal electrons trapped in the rising magnetic rope and flare arcade at the moment when the flare magnetic reconnection starts. The results are discussed in a scenario with the plasmoid in the rising magnetic rope.

We present a new investigation of the dependence of H2O maser detection rates and properties on the mid-IR AGN luminosity, L_AGN, and the obscuring column density, N_H, based on mid-IR and hard X-ray photometry. Based on spectral energy distribution fitting that allows for decomposition of the black hole accretion and star-formation components in the mid-infrared, we show that the megamaser (disk maser) detection rate increases sharply for galaxies with 12 micron AGN luminosity L^{AGN}_{12 micron} greater than 10^42 erg/s, from ~<3%(~<2%) to ~12%(~5%). By using the ratio of the observed X-ray to mid-IR AGN luminosity as an indicator of N_H, we also find that the megamaser (disk maser) detection rates are boosted to 15%(7%) and 20%(9%) for galaxies with N_H >= 10^23 cm^{-2} and N_H >= 10^{24} cm^{-2}, respectively. Combining these column density cuts with a constraint for high L^{AGN}_{12 micron} (>=10^42 erg/s) predicts further increases in the megamaser (disk maser) detection rates to 19%(8%) and 27%(14%), revealing unprecedented potential boosts of the megamaser and disk maser detection rates by a factor of 7-15 relative to the current rates, depending on the chosen sample selection criteria. A noteworthy aspect of these new predictions is that the completeness rates are only compromised mildly, with the rates remaining at the level of ~95%(~50%) for sources with N_H >= 10^{23} cm^{-2} (N_H >= 10^24 cm^-2). Applying these selection methods to the current X-ray AGN surveys predicts the detection of >~15 new megamaser disks.

By combining the photometric, spectroscopic, and astrometric information of the stars in the sightline of SNRs, the distances to and the extinctions of 32 Galactic supernova remnants (SNRs) are investigated. The stellar atmospheric parameters are from the SDSS$-$DR14$/$APOGEE and LAMOST$-$DR5$/$LEGUE spectroscopic surveys. The multi-band photometry, from optical to infrared, are collected from the {\it Gaia}, APASS, Pan--STARRS1, 2MASS, and {\it WISE} surveys. With the calibrated {\it Gaia} distances of individual stars, the distances to 15 of 32 SNRs are well determined from their produced extinction and association with molecular clouds. The upper limits of distance are derived for 3 SNRs. The color excess ratios $E(g_{\rm P1}-\lambda) / E(g_{\rm P1}-r_{\rm P1})$ of 32 SNRs are calculated, and their variation with wavebands is fitted by a simple dust model. The inferred dust grain size distribution bifurcates: while the graphite grains have comparable size to the average ISM dust, the silicate grains are generally larger. Along the way, the average extinction law from optical to near-infrared of the Milky Way is derived from the 1.3 million star sample and found to agree with the CCM89 law with $R_{\rm V}=3.15$.

We present a catalog of 5,290 RR Lyrae stars (RRLs) with metallicities estimated from spectra of the LAMOST Experiment for Galactic Understanding and Exploration (LEGUE) and the Sloan Extension for Galactic Understanding and Exploration (SEGUE) surveys. Nearly 70 per cent of them (3,642 objects) also have systemic radial velocities measured. Given the pulsating nature of RRLs, metallicity estimates are based on spectra of individual exposures, by matching them with the synthetic templates. The systemic radial velocities are measured by fitting the observed velocity as a function of phase assuming an empirical pulsating velocity template curve. Various tests show that our analyses yield metallicities with a typical precision of 0.20\,dex and systemic radial velocities with uncertainties ranging from 5 to 21\,km\,s$^{-1}$ (depending on the number of radial velocity measurements available for a given star). Based on the well calibrated near-infrared $PM_{W1}Z$ or $PM_{K_{\rm s}}Z$, and $M_{V}$-[Fe/H] relations, precise distances are derived for these RRLs. Finally, we include Gaia DR2 proper motions in our catalog. The catalog should be very useful for various Galactic studies, especially of the Galactic halo.

We present a large, homogeneous catalogue of molecular clouds within 4 kpc from the Sun at low Galactic latitudes ($|b|$ $<$ 10\degr) with unprecedented accurate distance determinations. Based on the three-dimensional dust reddening map and estimates of colour excesses and distances of over 32 million stars presented in Chen et al, we have identified 567 dust/molecular clouds with a hierarchical structure identification method and obtained their distance estimates by a dust model fitting algorithm. The typical distance uncertainty is less than 5 per cent. As far as we know, this is the first large catalogue of molecular clouds in the Galactic plane with distances derived in a direct manner. The clouds are seen to lie along the Sagittarius, Local and Perseus Arms. In addition to the known structures, we propose the existence of a possible {\it spur}, with a pitch angle of about 34\degr, connecting the Local and the Sagittarius Arms in the fourth quadrant. We have also derived the physical properties of those molecular clouds. The distribution of cloud properties in different parameter spaces agrees grossly with the previous results. Our cloud sample is an ideal starting point to study the concentration of dust and gas in the solar vicinity and their star formation activities.

We present a catalogue of 3,305 H$\alpha$ emission-line point sources observed with the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) in the vicinity fields of M31 and M33 during September 2011 and January 2016. The catalogue contains 1,487 emission-line stars, 532 emission-line nebulae including 377 likely planetary nebulae (PNe), 83 H~{\textsc{ii}} regions candidates and 20 possible supernovae remnants (SNRs) and 1,286 unknown objects. Among them, 24 PN candidates, 19 H~{\sc ii} region candidates, 10 SNR candidates and 1 symbiotic star candidate are new discoveries. Radial velocities and fluxes estimated from the H$\alpha$ line and those quantities of seven other major emission lines including H$\beta$, [O~{\textsc{iii}}]~$\lambda$4959, [O~{\textsc{iii}}]~$\lambda$5007, [N~{\textsc{ii}}]~$\lambda$6548, [N~{\textsc{ii}}]~$\lambda$6583, [S~{\textsc{ii}}]~$\lambda$6717 and [S~{\textsc{ii}}]~$\lambda$6731 lines of all the catalogued sources yielded from the LAMOST spectra are also presented in our catalogue. Our catalogue is an ideal starting point to study the chemistry properties and kinematics of M31 and M33.

We show in detail that the recently derived expression for evaluating the integrated Sachs-Wolfe (ISW) temperature shift in the cosmic microwave background (CMB) caused by individual embedded (compensated) lenses is equivalent to the conventional approach for flat background cosmologies. The conventional approach requires evaluating an integral of the time derivative of the lensing potential, whereas the new Fermat potential approach is simpler and only requires taking a derivative of the potential part of the time delay.

Nonthermal loop-top sources in solar flares are the most prominent observational signature that suggests energy release and particle acceleration in the solar corona. Although several scenarios for particle acceleration have been proposed, the origin of the loop-top sources remains unclear. Here we present a model that combines a large-scale magnetohydrodynamic simulation of a two-ribbon flare with a particle acceleration and transport model for investigating electron acceleration by a fast-mode termination shock at the looptop. Our model provides spatially resolved electron distribution that evolves in response to the dynamic flare geometry. We find a concave-downward magnetic structure located below the flare termination shock, induced by the fast reconnection downflows. It acts as a magnetic trap to confine the electrons at the looptop for an extended period of time. The electrons are energized significantly as they cross the shock front, and eventually build up a power-law energy spectrum extending to hundreds of keV. We suggest that this particle acceleration and transport scenario driven by a flare termination shock is a viable interpretation for the observed nonthermal loop-top sources.

We perform analysis of the three-dimensional kinematics of Milky Way disk stars in mono-age populations. We focus on stars between Galactocentric distances of $R=6$ and 14 \,kpc, selected from the combined LAMOST DR4 red clump giant stars and Gaia DR2 proper motion catalogue. We confirm the 3D asymmetrical motions of recent works and provide time tagging of the Galactic outer disk asymmetrical motions near the anticenter direction out to Galactocentric distances of 14\,kpc. Radial Galactocentric motions reach values up to 10 km s$^{-1}$, depending on the age of the population, and present a north-south asymmetry in the region corresponding to density and velocity substructures that were sensitive to the perturbations in the early 6 \,Gyr. After that time, the disk stars in this asymmetrical structure have become kinematically hotter, and are thus not sensitive to perturbations, and we find the structure is a relatively younger population. With quantitative analysis, we find stars both above and below the plane at $R\gtrsim 9$ kpc that exhibit bending mode motions of which the sensitive duration is around 8 \,Gyr. We speculate that the in-plane asymmetries might not be mainly caused by a fast rotating bar, intrinsically elliptical outer disk, secular expansion of the disk, or streams. Spiral arm dynamics, out-of-equilibrium models, minor mergers or others are important contributors. Vertical motions might be dominated by bending and breathing modes induced by complicated inner or external perturbers. It is likely that many of these mechanisms are coupled together.

Solar flare termination shocks have been suggested as one of the promising drivers for particle acceleration in solar flares, yet observational evidence remains rare. By utilizing radio dynamic spectroscopic imaging of decimetric stochastic spike bursts in an eruptive flare, Chen et al. found that the bursts form a dynamic surface-like feature located at the ending points of fast plasma downflows above the looptop, interpreted as a flare termination shock. One piece of observational evidence that strongly supports the termination shock interpretation is the occasional split of the emission band into two finer lanes in frequency, similar to the split-band feature seen in fast-coronal-shock-driven type II radio bursts. Here we perform spatially, spectrally, and temporally resolved analysis of the split-band feature of the flare termination shock event. We find that the ensemble of the radio centroids from the two split-band lanes each outlines a nearly co-spatial surface. The high-frequency lane is located slightly below its low frequency counterpart by ~0.8 Mm, which strongly supports the shock upstream-downstream interpretation. Under this scenario, the density compression ratio across the shock front can be inferred from the frequency split, which implies a shock with a Mach number of up to 2.0. Further, the spatiotemporal evolution of the density compression along the shock front agrees favorably with results from magnetohydrodynamics simulations. We conclude that the detailed variations of the shock compression ratio may be due to the impact of dynamic plasma structures in the reconnection outflows, which results in distortion of the shock front.

The hot massive luminous blue variables (LBVs) represent an important evolutionary phase of massive stars. Here, we report the discovery of a new LBV -- LAMOST J0037+4016 in the distant outskirt of the Andromeda galaxy. It is located in the south-western corner (a possible faint spiral arm) of M31 with an unexpectedly large projection distance of $\sim$ 22 kpc from the center. The optical light curve shows a 1.2 mag variation in $V$ band and its outburst and quiescence phases both last over several years. The observed spectra indicate an A-type supergiant at epoch close to the outburst phase and a hot B-type supergiant with weak [Fe II] emission lines at epoch of much dimmer brightness. The near-infrared color-color diagram further shows it follows the distribution of Galactic and M31 LBVs rather than B[e] supergiants. All the existing data strongly show that LAMOST J0037+4016 is an LBV. By spectral energy distribution fitting, we find it has a luminosity ($4.42 \pm 1.64$)$\times 10^5$ $L_{\odot}$ and an initial mass $\sim 30$ $M_{\odot}$, indicating its nature of less luminosity class of LBV.

We investigate the three-dimensional asymmetrical kinematics and present time stamps of the Milky Way disk between Galactocentric distances of $R=12$ and 15 \,kpc, using red clump stars selected from the LAMOST Galactic survey, also with proper motion measurements provided by the Gaia DR2. We discover velocity substructure above the Galactic plane corresponding to a density dip found recently ("south-middle opposite" density structure[R $\sim$ 12-15 \,kpc, Z $\sim$ 1.5 \,kpc] discovered in \citet{Wang2018b, Wang2018c}) in the radial and azimuthal velocity. For the vertical velocity, we detect clear vertical bulk motions or bending mode motions, which has no clear north-south asymmetry corresponding to the in-plane asymmetrical features. In the subsample of stars with different ages, we find that there is little temporal evolution of the in-plane asymmetry from 0$-$14 \,Gyr, which means the structure is sensitive to the perturbations in almost cosmic time possibly. We propose that the possible scenario of this asymmetric velocity structure is caused by the mechanisms generated in-plane, rather than vertical perturbations.

Gaia DR2 has revealed new small-scale and large-scale patterns in the phase-space distribution of stars in the Milky Way. In cylindrical Galactic coordinates $(R,\phi,z)$, ridge-like structures can be seen in the \vphiR{} plane and asymmetric arch-like structures in the \vphivR{} plane. We show that the ridges are also clearly present when the third dimension of the \vphiR{} plane is represented by $\langle z \rangle$, $\langle V_z \rangle$, $\langle V_R \rangle$, $\langle$[Fe/H]$\rangle$ and $\langle[\alpha/{\rm Fe}]\rangle$. The maps suggest that stars along the ridges lie preferentially close to the Galactic midplane ($|z|<0.2$ kpc), and have metallicity and $\alpha$ elemental abundance similar to that of the Sun. We show that phase mixing of disrupting spiral arms can generate both the ridges and the arches. It also generates discrete groupings in orbital energy $-$ the ridges and arches are simply surfaces of constant energy. We identify 8 distinct ridges in the \gaia{} data: six of them have constant energy while two have constant angular momentum. Given that the signature is strongest for stars close to the plane, the presence of ridges in $\langle z \rangle$ and $\langle V_z \rangle$ suggests a coupling between planar and vertical directions. We demonstrate, using N-body simulations that such coupling can be generated both in isolated discs and in discs perturbed by an orbiting satellite like the Sagittarius dwarf galaxy.

We present three dimensional (3D) dust mapping of 12 supernova remnants (SNRs) in the Galactic anti-center (Galactic longitude $l$ between 150\degr\ and 210\degr) based on a recent 3D interstellar extinction map. The dust distribution of the regions which cover the full extents in the radio continuum for the individual SNRs are discussed. Four SNRs show significant spatial coincidences between molecular clouds (MCs) revealed from the 3D extinction mapping and the corresponding radio features. The results confirm the interactions between these SNRs and their surrounding MCs. Based on these correlations, we provide new distance estimates of the four SNRs, G189.1+3.0 (IC443, $d=1729^{+116}_{-94} \rm \,pc$), G190.9-2.2 ($d=1036^{+17}_{-81} \rm \,pc$), G205.5+0.5 ($d=941^{+96}_{-94}$ or $1257^{+92}_{-101} \rm \,pc$) and G213.0-0.6 ($d=1146^{+79}_{-80} \rm \,pc$). In addition, we find indirect evidences of potential interactions between SNRs and MCs for three other SNRs. New distance constraints are also given for these three SNRs.

We use two unsupervised machine learning algorithms, Shared Nearest Neighbor (SNN) and \textit{EnLink}, as a single approach to identify stellar groupings in the Orion star-forming complex as an application to the 5-dimensional astrometric data from Gaia DR2. The algorithms present two different ways to limit user bias when evaluating the relative weights among the astrometric parameters, automatically determined by the machine and through a standard procedure by monitoring several outcome measures. Both algorithms complement each other and produce similar stellar groups. Because SNN groups have a much smaller spread in proper motions compared to \textit{EnLink}, we use \textit{EnLink}, which requires no input, as a first pass tool for group identification and validation. We then used the SNN algorithm to dissect the Orion star-forming complex. We identify 21 spatially- and kinematically-coherent groups in the Orion complex, nine of which previously unknown. The groups show a wide distribution of distances extending as far as about 150 pc in front of the star-forming Orion clouds, to about 50 pc beyond them where we find, unexpectedly, three groups. Our results expose to view the wealth of sub-structure in the OB association, within and beyond the classical Blaauw Orion OBI sub-groups. A full characterization of the new groups is of the essence as it offers the potential to unveil how star formation proceeds globally in large complexes such as Orion. The data and code that generated the groups is provided in this Letter.

Accurate determinations of atmospheric parameters (effective temperature $T_{\rm eff}$, surface gravity log $g$ and metallicity [Fe/H]) and distances for large complete samples are of vital importance for various Galactic studies. We have developed a photometric method to select red giant stars and estimate their atmospheric parameters from the photometric colors provided by the SkyMapper Southern Survey (SMSS) data release (DR) 1.1, using stars in common with the LAMOST Galactic spectroscopic surveys as a training set. Distances are estimated with two different approaches: one based on the Gaia DR2 parallaxes for nearby ($d \leq 4.5$ kpc) bright stars and another based on the absolute magnitudes predicted by intrinsic color $(g-i)_0$ and photometric metallicity [Fe/H] for distant ($d > 4.5$ kpc) faint stars. Various tests show that our method is capable of delivering atmospheric parameters with a precision of $\sim$80 K for $T_{\rm eff}$, $\sim$0.18 dex for [Fe/H] and $\sim$0.35 dex for log $g$, but with a significant systematic error at log $g \sim$ 2.3. For distances delivered from $(g-i)_0$ and photometric [Fe/H], our test with the member stars of globular clusters show a median uncertainty of 16 per cent with a negligible zero-point offset. Using this method, atmospheric parameters and distances of nearly one million red giant stars are derived from SMSS DR1.1. Proper motion measurements from Gaia DR2 are available for almost all of the red giant stars, and radial velocity measurements from several large spectroscopic surveys are available for 44 per cent of these. This sample will be accessible online at https://yanghuang0.wixsite.com/yangh/research .

We present an analysis of Chandra spectra of five gravitationally lensed active galactic nuclei. We confirm the previous detections of FeK$\alpha$ emission lines in most images of these objects with high significance. The line energies range from 5.8 to 6.8 keV with widths from unresolved to 0.6 keV, consistent with emission close to spinning black holes viewed at different inclination angles. We also confirm the positive offset from the Iwasawa-Taniguchi effect, the inverse correlation between the FeK$\alpha$ equivalent width and the X-ray luminosity in AGN, where our measured equivalent widths are larger in lensed quasars. We attribute this effect to microlensing, and perform a microlensing likelihood analysis to constrain the emission size of the relativistic reflection region and the spin of supermassive black holes, assuming that the X-ray corona and the reflection region, responsible for the iron emission line, both follow power-law emissivity profiles. The microlensing analysis yields strong constraints on the spin and emissivity index of the reflection component for Q2237+0305, with $a > 0.92$ and $n > 5.4$. For the remaining four targets, we jointly constrain the two parameters, yielding $a=0.8\pm0.16$ and an emissivity index of $n=4.0\pm 0.8$, suggesting that the relativistic X-ray reflection region is ultra-compact and very close to the innermost stable circular orbits of black holes, which are spinning at close to the maximal value. We successfully constrain the half light radius of the emission region to $< 2.4$ $r_g$ ($r_g = GM/c^2$) for Q2237+0305 and in the range 5.9--7.4 $r_g$ for the joint sample.

Asteroseismology is a promising tool to study Galactic structure and evolution because it can probe the ages of stars. Earlier attempts comparing seismic data from the {\it Kepler} satellite with predictions from Galaxy models found that the models predicted more low-mass stars compared to the observed distribution of masses. It was unclear if the mismatch was due to inaccuracies in the Galactic models, or the unknown aspects of the selection function of the stars. Using new data from the K2 mission, which has a well-defined selection function, we find that an old metal-poor thick disc, as used in previous Galactic models, is incompatible with the asteroseismic information. We show that spectroscopic measurements of [Fe/H] and [$\alpha$/Fe] elemental abundances from the GALAH survey indicate a mean metallicity of $\log (Z/Z_{\odot})=-0.16$ for the thick disc. Here $Z$ is the effective solar-scaled metallicity, which is a function of [Fe/H] and [$\alpha$/Fe]. With the revised disc metallicities, for the first time, the theoretically predicted distribution of seismic masses show excellent agreement with the observed distribution of masses. This provides an indirect verification of the asteroseismic mass scaling relation is good to within five percent. Using an importance-sampling framework that takes the selection function into account, we fit a population synthesis model of the Galaxy to the observed seismic and spectroscopic data. Assuming the asteroseismic scaling relations are correct, we estimate the mean age of the thick disc to be about 10 Gyr, in agreement with the traditional idea of an old $\alpha$-enhanced thick disc.

The Sun is an active star that can have a direct impact on the Earth, its magnetosphere, and the technological infrastructure on which modern society depends. Among the phenomena that drive "space weather" are fast solar wind streams and co-rotating interaction regions, solar flares, coronal mass ejections, the shocks they produce, and the energetic particles they accelerate. Radio emission from these and associated phenomena offer unique diagnostic possibilities that complement those available at other wavelengths. Here, the relevant space weather drivers are briefly described, the potential role of radio observations is outlined, and the requirements of an instrument to provide them are provided: specifically, ultrabroadband imaging spectropolarimetry. The insights provided by radio observations of space weather drivers will not only inform the science of space weather, they will pave the way for new tools for forecasting and "nowcasting" space weather. They will also serve as an important touchstone against which local environment of exoplanets and the impact of "exo-space weather" can be evaluated.

Identification of the mechanisms responsible for heating the solar chromosphere and corona remains an outstanding problem, one of great relevance to late-type stars as well. There has been tremendous progress in the past decade, largely driven by new instruments, new observations, and sophisticated modeling efforts. Despite this progress, gaps remain. We briefly discuss the need for radio coverage of the 3D solar atmosphere and discuss the requirements.

With the photometric data from the SDSS survey, the spectroscopic data from the SDSS/SEGUE and the LAMOST surveys, and the astrometric data from the Gaia DR2, we have identified 71 highly-probable member stars of the GD-1 cold stellar stream spread along almost its entire length (i.e. from 126 to 203 degree in Right Ascension). With the accurate spectroscopic (i.e. metallicity and line-of-sight velocity) and astrometric (i.e. proper motions) information, the position-velocity diagrams, i.e. $\phi_{1}$-$\mu_{\alpha}$, $\phi_{1}$-$\mu_{\delta}$ and $\phi_{1}$-${v_{\rm gsr}}$, of the GD-1 stream are well mapped. The stream has an average metallicity [Fe/H] $= -1.95 \pm 0.23$. The rich information of member stars of the stream now available allow one not only to model its origin, but also to place strong constraints on the mass distribution and the gravitational potential of the Milky Way.

Magnetic reconnection is a fundamental physical process in many laboratory, space, and astrophysical plasma contexts. Solar flares serve as an outstanding laboratory to study the magnetic reconnection and the associated energy release and conversion processes under plasma conditions difficult to reproduce in the laboratory, and with considerable spatiotemporal details not possible elsewhere in astrophysics. Here we emphasize the unique power of remote-sensing observations of solar flares at radio wavelengths. In particular, we discuss the transformative technique of broadband radio dynamic imaging spectroscopy in making significant contributions to addressing several outstanding challenges in magnetic reconnection, including the capability of pinpointing magnetic reconnection sites, measuring the time-evolving reconnecting magnetic fields, and deriving the spatially and temporally resolved distribution function of flare-accelerated electrons.

We have investigated the distributions of stellar azimuthal and radial velocity components $V_{\Phi}$ and $V_{R}$ in the vertical position-velocity plane $Z$-$V_{Z}$ across the Galactic disc of $6.34 \lesssim R \lesssim 12.34$\,kpc and $|\Phi| \lesssim 7.5^{\circ}$ using a Gaia and Gaia-LAMOST sample of stars. As found in previous works, the distributions exhibit significant spiral patterns. The $V_{R}$ distributions also show clear quadrupole patterns, which are the consequence of the well-known tilt of the velocity ellipsoid. The observed spiral and quadrupole patterns in the phase space plane vary strongly with radial and azimuthal positions. The phase spirals of $V_{\Phi}$ become more and more relaxed as $R$ increases. The spiral patterns of $V_{\Phi}$ and $V_{R}$ and the quadrupole patterns of $V_{R}$ are strongest at $-2^{\circ} < \Phi < 2^{\circ}$ but negligible at $4^{\circ} < \Phi < 6^{\circ}$ and $-6^{\circ} < \Phi < -4^{\circ}$. Our results suggest an external origin of the phase spirals. In this scenario, the intruder, most likely the previously well-known Sagittarius dwarf galaxy, passed through the Galactic plane in the direction towards either Galactic center or anti-center. The azimuthal variations of the phase spirals also help us constrain the passage duration of the intruder. A detailed model is required to reproduce the observed radial and azimuthal variations of the phase spirals of $V_{\Phi}$ and $V_{R}$.

Exotic compact objects (ECOs) have recently become an exciting research subject, since they are speculated to have a special response to the incident gravitational waves (GWs) that leads to GW echoes. We show that energy carried by GWs can easily cause the event horizon to form out of a static ECO --- leaving no echo signals towards spatial infinity. To show this, we use the ingoing Vaidya spacetime and take into account the back reaction due to incoming GWs. Demanding that an ECO does not collapse into a black hole puts an upper bound on the compactness of the ECO, at the cost of less distinct echo signals for smaller compactness. The trade-off between echoes' detectability and distinguishability leads to a fine tuning of ECO parameters for LIGO to find distinct echoes. We also show that an extremely compact ECO that can survive the gravitational collapse and give rise to GW echoes might have to expand its surface in a non-causal way.

Solar flares involve the sudden release of magnetic energy in the solar corona. Accelerated nonthermal electrons have often been invoked as the primary means for transporting the bulk of the released energy to the lower solar atmosphere. However, significant challenges remain for this scenario, especially in accounting for the large number of accelerated electrons inferred from observations. Propagating magnetohydrodynamics (MHD) waves, particularly those with subsecond/second-scale periods, have been proposed as an alternative means for transporting the released flare energy, likely alongside the electron beams, while observational evidence remains elusive. Here we report a possible detection of such waves in the late impulsive phase of a two-ribbon flare. This is based on ultrahigh cadence dynamic imaging spectroscopic observations of a peculiar type of decimetric radio bursts obtained by the Karl G. Jansky Very Large Array. Radio imaging at each time and frequency pixel allows us to trace the spatiotemporal motion of the source, which agrees with the implications of the frequency drift pattern in the dynamic spectrum. The radio source, propagating at 1000--2000 km s$^{-1}$ in projection, shows close spatial and temporal association with transient brightenings on the flare ribbon. In addition, multitudes of subsecond-period oscillations are present in the radio emission. We interpret the observed radio bursts as short-period MHD wave packets propagating along newly reconnected magnetic flux tubes linking to the flare ribbon. The estimated energy flux carried by the waves is comparable to that needed to account for the plasma heating during the late impulsive phase of this flare.

We present a catalog of stellar age and mass estimates for a sample of 640\,986 red giant branch (RGB) stars of the Galactic disk from the LAMOST Galactic Spectroscopic Survey (DR4). The RGB stars are distinguished from the red clump stars utilizing period spacing derived from the spectra with a machine learning method based on kernel principal component analysis (KPCA). Cross-validation suggests our method is capable of distinguishing RC from RGB stars with only 2 per cent contamination rate for stars with signal-to-noise ratio (SNR) higher than 50. The age and mass of these RGB stars are determined from their LAMOST spectra with KPCA method by taking the LAMOST - $Kepler$ giant stars having asteroseismic parameters and the LAMOST-TGAS sub-giant stars based on isochrones as training sets. Examinations suggest that the age and mass estimates of our RGB sample stars with SNR $>$ 30 have a median error of 30 per cent and 10 per cent, respectively. Stellar ages are found to exhibit positive vertical and negative radial gradients across the disk, and the age structure of the disk is strongly flared across the whole disk of $6<R<13$\,kpc. The data set demonstrates good correlations among stellar age, [Fe/H] and [$\alpha$/Fe]. There are two separate sequences in the [Fe/H] -- [$\alpha$/Fe] plane: a high--$\alpha$ sequence with stars older than $\sim$\,8\,Gyr and a low--$\alpha$ sequence composed of stars with ages covering the whole range of possible ages of stars. We also examine relations between age and kinematic parameters derived from the Gaia DR2 parallax and proper motions. Both the median value and dispersion of the orbital eccentricity are found to increase with age. The vertical angular momentum is found to fairly smoothly decrease with age from 2 to 12\,Gyr, with a rate of about $-$50\,kpc\,km\,s$^{-1}$\,Gyr$^{-1}$. A full table of the catalog is public available online.

We present the chemodynamic structure of the solar neighbourhood using 62 814 stars within a 500 pc sphere of the Sun observed by GALAH and with astrometric parameters from Gaia DR2. We measure the velocity dispersion for all three components (vertical, radial, and tangential) and find that it varies smoothly with [Fe/H] and [$\alpha$/Fe] for each component. The vertical component is especially clean, with $\sigma_{v_z}$ increasing from a low of $8$ km s$^{-1}$ at solar-[$\alpha$/Fe] and [Fe/H] to a high of more than 50 km s$^{-1}$ for more metal-poor and [$\alpha$/Fe] enhanced populations. We find no evidence of a large decrease in the velocity dispersion of the higher-[$\alpha$/Fe] populations as claimed in analysis prior to Gaia DR2, although the trend of increasing velocity dispersion with [$\alpha$/Fe] for the same metallicity does significantly flatten at high-[$\alpha$/Fe]. The eccentricity distribution for local stars varies most strongly as a function of [$\alpha$/Fe], where stars with [$\alpha$/Fe]$<0.1$ dex having generally circular orbits ($e<0.15$), while the median eccentricity increases rapidly for more [$\alpha$/Fe]-enhanced stellar populations up to $e\sim0.35$. These [$\alpha$/Fe]-enhanced populations have guiding radii consistent with origins in the inner Galaxy. Of the stars with metallicities much higher than the local ISM ([Fe/H]>0.1 dex), we find that more than 70\% have $e<0.2$ and are likely observed in the solar neighbourhood through churning/migration rather than blurring effects, as the epicyclic motion for these stars is not large enough to reach the radii at which they were likely born based on their metallicity.

We present new three-dimensional (3D) interstellar dust reddening maps of the Galactic plane in three colours, E(G-Ks), E(Bp-Rp) and E(H-Ks). The maps have a spatial angular resolution of 6 arcmin and covers over 7000 deg$^2$ of the Galactic plane for Galactic longitude 0 deg $<$ $l$ $<$ 360 deg and latitude $|b|$ $<$ $10$ deg. The maps are constructed from robust parallax estimates from the Gaia Data Release 2 (Gaia DR2) combined with the high-quality optical photometry from the Gaia DR2 and the infrared photometry from the 2MASS and WISE surveys. We estimate the colour excesses, E(G-Ks), E(Bp-Rp) and E(H-Ks), of over 56 million stars with the machine learning algorithm Random Forest regression, using a training data set constructed from the large-scale spectroscopic surveys LAMOST, SEGUE and APOGEE. The results reveal the large-scale dust distribution in the Galactic disk, showing a number of features consistent with the earlier studies. The Galactic dust disk is clearly warped and show complex structures possibly spatially associated with the Sagittarius, Local and Perseus arms. We also provide the empirical extinction coefficients for the Gaia photometry that can be used to convert the colour excesses presented here to the line-of-sight extinction values in the Gaia photometric bands.

In eruptive solar flares, termination shocks (TSs), formed when high-speed reconnection outflows collide with closed dense flaring loops, are believed to be one of the possible candidates for plasma heating and particle acceleration. In this work, we perform resistive magnetohydrodynamic simulations in a classic Kopp-Pneuman flare configuration to study the formation and evolution of TSs, and analyze in detail the dynamic features of TSs and variations of the shock strength in space and time. This research focuses on the fast reconnection phase when plasmoids form and produce small-scale structures inside the flare current sheet. It is found that the TS emerges once the downward outflow colliding with closed magnetic loops becomes super-magnetosonic, and immediately becomes highly dynamical. The morphology of a TS can be flat, oblique, or curved depending on the detailed interactions between the outflows/plasmoids and the highly dynamic plasma in the looptop region. The TS becomes weaker when a plasmoid is crossing through, or may even be destroyed by well developed plasmoids and then re-constructed above the plasmoids. We also perform detailed statistical analysis on important physical quantities along and across the shock front. The density and temperature ratios range from 1 to 3 across the TS front, and the pressure ratio typically has larger values up to 10. We show that weak guide fields do not strongly affect the Mach number and compression ratios, and the TS length becomes slightly larger in the case with thermal conduction.

The Advanced Space-based Solar Observatory (ASO-S) mission aims to explore two most spectacular eruptions in the Sun: solar flares and coronal mass ejections (CMEs), and their magnetism. For the studies of CMEs, the payload Lyman-alpha Solar Telescope (LST) has been proposed. It includes a traditional white-light coronagraph and a Lyman-alpha coronagraph which opens a new window to CME observations. Polarization measurements taken by white-light coronagraphs are crucial to derive fundamental physical parameters of CMEs. To make such measurements, there are two options of Stokes polarimeter which have been used by existing white-light coronagraphs for space missions. One uses a single or triple linear polarizers, the other involves both a half-wave plate and a linear polarizer. We find that the former option subjects to less uncertainty in the derived Stokes vector propagated from detector noise. The latter option involves two plates which are prone to internal reflections and may have a reduced transmission factor. Therefore, the former option is adopted as our Stokes polarimeter scheme for LST. Based on the parameters of the intended linear polarizer(s) colorPol provided by CODIXX and the half-wave plate 2-APW-L2-012C by Altechna, it is further shown that the imperfect maximum transmittance of the polarizer significantly increases the variance amplification of Stokes vector by at least about 50% when compared with the ideal case. The relative errors of Stokes vector caused by the imperfection of colorPol polarizer and the uncertainty due to the polarizer assembling in the telescope are estimated to be about 5%. Among the considered parameters, we find that the dominant error comes from the uncertainty in the maximum transmittance of the polarizer.

Solar flares are due to the catastrophic release of magnetic energy in the Sun's corona, resulting in plasma heating, mass motions, particle acceleration, and radiation emitted from radio to $\gamma$-ray wavelengths. They are associated with global coronal eruptions of plasma into the interplanetary medium---coronal mass ejections---that can result in a variety of "space weather" phenomena. Flares release energy over a vast range of energies, from $\sim\!10^{23}$ ergs (nanoflares) to more than $10^{32}$ ergs. Solar flares are a phenomenon of general astrophysical interest, allowing detailed study of magnetic energy release, eruptive processes, shock formation and propagation, particle acceleration and transport, and radiative processes. Observations at radio wavelengths offer unique diagnostics of the physics of flares. To fully exploit these diagnostics requires the means of performing time-resolved imaging spectropolarimetry. Recent observations with the Jansky Very Large Array (JVLA) and the Expanded Owens Valley Solar Array (EOVSA), supported by extensive development in forward modeling, have demonstrated the power of the approach. The ngVLA has the potential to bring our understanding of flare processes to a new level through its combination of high spatial resolution, broad frequency range, and imaging dynamic range---especially when used in concert with multi-wavelength observations and data at hard X-ray energies.

We have investigated the metallicity distributions of mono-age stellar populations across the disc of $6 \lesssim R \lesssim 12$\,kpc and $|Z| \lesssim 2$\,kpc using samples selected from the main-sequence turn-off and sub-giant (MSTO-SG) stars targeted by the LAMOST Galactic Spectroscopic surveys. Both the mean values and the profiles of the distributions exhibit significant variations with age and position. We confirm that the oldest ($>11$\,Gyr) stars have nearly flat radial [Fe/H] gradients at all heights above the disc but show negative vertical [Fe/H] gradients. For stars younger than 11\,Gyr, the radial [Fe/H] gradients flatten with $|Z|$, while the vertical [Fe/H] gradients flatten with $R$. Stars of 4--6\,Gyr exhibit steeper negative radial [Fe/H] gradients than those of either younger or older ages. Values of [$\alpha$/Fe] of mono-age stellar populations also show significant radial and vertical gradients, with patterns varying with age. The [Fe/H] distribution profiles of old ($>8$\,Gyr) stars vary little with $R$, while those of younger stars exhibit strong radial variations, probably a consequence of significant radial migration. The [$\alpha$/Fe] radial distribution profiles show opposite patterns of variations with age compared to those of [Fe/H]. We have also explored the impacts of stellar mixing by epicycle motions (blurring) on the [Fe/H] and [$\alpha$/Fe] distributions, and found that blurring mainly change the widths of the distribution profiles. Our results suggest that the disc may have experienced a complex assemblage history, in which both the "inside-out" and "upside-down" formation processes may have played an important role.

Magnetic reconnection, the central engine that powers explosive phenomena throughout the Universe, is also perceived as one of the principal mechanisms for accelerating particles to high energies. Although various signatures of magnetic reconnection have been frequently reported, observational evidence that links particle acceleration directly to the reconnection site has been rare, especially for space plasma environments currently inaccessible to $\textit{in situ}$ measurements. Here we utilize broadband radio dynamic imaging spectroscopy available from the Karl G. Jansky Very Large Array to observe decimetric type III radio bursts in a solar jet with high angular ($\sim$20$''$), spectral ($\sim$1 %), and temporal resolution (50 milliseconds). These observations allow us to derive detailed trajectories of semi-relativistic (tens of keV) electron beams in the low solar corona with unprecedentedly high angular precision ($<0''.65$). We found that each group of electron beams, which corresponds to a cluster of type III bursts with 1-2-second duration, diverges from an extremely compact region ($\sim$600 km$^2$) in the low solar corona. The beam-diverging sites are located behind the erupting jet spire and above the closed arcades, coinciding with the presumed location of magnetic reconnection in the jet eruption picture supported by extreme ultraviolet/X-ray data and magnetic modeling. We interpret each beam-diverging site as a reconnection null point where multitudes of magnetic flux tubes join and reconnect. Our data suggest that the null points likely consist of a high level of density inhomogeneities possibly down to 10-km scales. These results, at least in the present case, strongly favor a reconnection-driven electron acceleration scenario.

This study expands the coverage and improves the homogeneity of the distribution of MILES template stars in the parameter space, as well as extends the wavelength coverage of the template spectra to the far red beyond the Ca\,{\sc{ii}} triplet. To achieve this we have carried out a major observational campaign using the OMR long-slit spectrograph mounted on the NAOC 2.16\,m telescope and the YFOSC long-slit spectrograph mounted on the YNAO 2.4\,m telescope. The original sample is based on the MILES library, supplemented by 918 stars selected from PASTEL database. In total, 822 OMR and 1,324 YFOSC spectra, covering respectively the wavelength ranges $\lambda\lambda$3800--5180 and $\lambda\lambda$5150--9000, have been collected and reduced. The spectra have a mean resolution FWHM (full-width at half-maximum) of $\sim 3.3$\,{\AA} and are wavelength- and flux-calibrated to an accuracy of $\sim 20$\,km\,s$^{-1}$ and $\sim 5$\,per\,cent, respectively. The spectra are further corrected for systematic errors in the wavelength calibration to an accuracy of $\sim 4$\,km\,s$^{-1}$ by cross-correlating with the theoretical spectra. Almost all the spectra have an average signal to noise ratio (SNR) better than 100 per pixel. Combined with the MILES spectra, there are now 1,731, 1,542, 1,324 and 1,273 stars with spectra covering respectively $\lambda\lambda$3800--5180, $\lambda\lambda$3800--7500, $\lambda\lambda$5150--9000 and $\lambda\lambda$3800--9000. This paper describes our template star selection, the observation and data reduction, and presents the reduced spectra collected hitherto.

We analyze the spectra of high temperature \fexxiv~lines observed by \emph{Hinode}/EIS during the impulsive phase of the X8.3--class flare on September 10, 2017. The line profiles are broad, show pronounced wings, and clearly depart from a single Gaussian shape. The lines can be well fitted with $\kappa$ distributions, with values of $\kappa$ varying between~$\approx$1.7 to 3. The regions where we observe the non-Gaussian profiles coincide with the location of high-energy ($\approx$100--300 keV) HXR sources observed by \emph{RHESSI}, suggesting the presence of particle acceleration or turbulence, also confirmed by the observations of a non-thermal microwave sources with the \emph{Expanded Owens Valley Solar Array} (EOVSA) at and above the HXR looptop source. We also investigate the effect of taking into account $\kappa$ distributions in the temperature diagnostics based on the ratio of the \fexxiii~263.76~\AA~and \fexxiv~255.1~\AA~EIS lines. We found that these lines can be formed at much higher temperatures than expected (up to log($T$\,[K])\,$\approx$\,7.8), if departures from Maxwellian distributions are taken into account. Although larger line widths are expected because of these higher formation temperatures, the observed line widths still imply non-thermal broadening in excess of ~200\,\kps. The non-thermal broadening related to HXR emission is better interpreted by turbulence rather than chromospheric evaporation.

We present a detailed determination and analysis of 3D stellar mass distribution of the Galactic disk for mono-age populations using a sample of 0.93 million main-sequence turn-off and subgiant stars from the LAMOST Galactic Surveys. Our results show (1) all stellar populations younger than 10\,Gyr exhibit strong disk flaring, which is accompanied with a dumpy vertical density profile that is best described by a $sech^n$ function with index depending on both radius and age; (2) Asymmetries and wave-like oscillations are presented in both the radial and vertical direction, with strength varying with stellar populations; (3) As a contribution by the Local spiral arm, the mid-plane stellar mass density at solar radius but 400--800\,pc (3--6$^\circ$) away from the Sun in the azimuthal direction has a value of $0.0594\pm0.0008$\,$M_\odot$/pc$^3$, which is 0.0164\,$M_\odot$/pc$^3$ higher than previous estimates at the solar neighborhood. The result causes doubts on the current estimate of local dark matter density; (4) The radial distribution of surface mass density yields a disk scale length evolving from $\sim$4\,kpc for the young to $\sim$2\,kpc for the old populations. The overall population exhibits a disk scale length of $2.48\pm0.05$\,kpc, and a total stellar mass of $3.6(\pm0.1)\times10^{10}$\,$M_\odot$ assuming $R_{\odot}=8.0$\,kpc, and the value becomes $4.1(\pm0.1)\times10^{10}$\,$M_\odot$ if $R_{\odot}=8.3$\,kpc; (5) The disk has a peak star formation rate ({\rm SFR}) changing from 6--8\,Gyr at the inner to 4--6\,Gyr ago at the outer part, indicating an inside-out assemblage history. The 0--1\,Gyr population yields a recent disk total {\rm SFR} of $1.96\pm0.12$\,$M_\odot$/yr.

Using a sample of about 123,000 stars with accurate 3D velocity measurements from the LAMOST-TGAS data, we confirm the kinematic signature of the Galactic warp found by Sch\"onrich & Dehnen recently. The data reveal a clear trend of increasing mean vertical velocity $\overline{V_{z}}$ as a function of absolute vertical angular momentum $L_{z}$ and azimuthal velocity $V_{\phi}$ for guiding center radius $R_{g}$ between 6.0 and 10.5 kpc. The trend is consistent with a large-scale Galactic warp. Similar to Sch\"onrich & Dehnen, we also find a wave-like pattern of $\overline{V_{z}}$ versus $L_{z}$ with an amplitude of $\sim 0.9$ km s$^{-1}$ on a scale of $\sim 2.0$ kpc, which could arise from bending waves or a winding warp. Finally, we confirm a prominent, localized peak in $\overline{V_z}$ near $L_z \sim 2150$ kpc km s$^{-1}$ (corresponding to $R_{g} \sim 9$ kpc and $V_{\phi} \sim 255$ km s$^{-1}$). The additional line-of-sight velocity information from LAMOST reveals that stars in this feature have a large, inward radial velocity of $V_{R} \sim -13.33 \pm 0.59$ km s$^{-1}$ and a small radial velocity dispersion $\sigma_{R} \sim 25.27 \pm 0.89$ km s$^{-1}$, suggesting that a stellar stream gives rise to this feature.

We report the first science results from the newly completed Expanded Owens Valley Solar Array (EOVSA), which obtained excellent microwave imaging spectroscopy observations of SOL2017-09-10, a classic partially-occulted solar limb flare associated with an erupting flux rope. This event is also well-covered by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) in hard X-rays (HXRs). We present an overview of this event focusing on microwave and HXR data, both associated with high-energy nonthermal electrons, and discuss them within the context of the flare geometry and evolution revealed by extreme ultraviolet (EUV) observations from the Atmospheric Imaging Assembly aboard the Solar Dynamics Observatory (SDO/AIA). The EOVSA and RHESSI data reveal the evolving spatial and energy distribution of high-energy electrons throughout the entire flaring region. The results suggest that the microwave and HXR sources largely arise from a common nonthermal electron population, although the microwave imaging spectroscopy provides information over a much larger volume of the corona.

We present a new catalogue of 18 080 radial velocity standard stars selected from the APOGEE data. These RV standard stars are observed at least three times and have a median stability ($3\sigma_{\rm RV}$) around 240 m s$^{-1}$ over a time baseline longer than 200 days. They are largely distributed in the northern sky and could be extended to the southern sky by the future APOGEE-2 survey. Most of the stars are red giants ($J - K_{\rm s} \ge 0.5$) owing to the APOGEE target selection criteria. Only about ten per cent of them are main-sequence stars. The $H$ band magnitude range of the stars is 7-12.5 mag with the faint limit much fainter than the magnitudes of previous RV standard stars. As an application, we show the new set of standard stars to determine the radial velocity zero points of the RAVE, the LAMOST {and the Gaia-RVS} Galactic spectroscopic surveys.

We have developed a novel technique based on a clustering algorithm which searches for kinematically- and chemically-clustered stars in the APOGEE DR12 Cannon data. As compared to classical chemical tagging, the kinematic information included in our methodology allows us to identify stars that are members of known globular clusters with greater confidence. We apply our algorithm to the entire APOGEE catalog of 150,615 stars whose chemical abundances are derived by the Cannon. Our methodology found anti-correlations between the elements Al and Mg, Na and O, and C and N previously identified in the optical spectra in globular clusters, even though we omit these elements in our algorithm. Our algorithm identifies globular clusters without a priori knowledge of their locations in the sky. Thus, not only does this technique promise to discover new globular clusters, but it also allows us to identify candidate streams of kinematically- and chemically-clustered stars in the Milky Way.

We present a new millimeter CO-line observation towards supernova remnant (SNR) CTB 87, which was regarded purely as a pulsar wind nebula (PWN), and an optical investigation of a coincident surrounding superbubble. The CO observation shows that the SNR delineated by the radio emission is projectively covered by a molecular cloud (MC) complex at V$_{\rm {LSR}}$ = $-60$ to $-54$ km s$^{-1}$. Both the symmetric axis of the radio emission and the trailing X-ray PWN appear projectively to be along a gap between two molecular gas patches at $-58$ to $-57$ km s$^{-1}$. Asymmetric broad profiles of $^{12}$CO lines peaked at $-58$ km s$^{-1}$ are found at the eastern and southwestern edges of the radio emission. This represents a kinematic signature consistent with an SNR-MC interaction. We also find that a superbubble, $\sim 37'$ in radius, appears to surround the SNR from HI 21cm (V$_{\rm {LSR}} \sim -61$ to $-68$ km s$^{-1}$), WISE mid-IR, and optical extinction data. We build a multi-band photometric stellar sample of stars within the superbubble region and find 82 OB star candidates. The likely peak distance in the stars' distribution seems consistent with the distance previously suggested for CTB 87. We suggest the arc-like radio emission is mainly the relic of the part of blastwave that propagates into the MC complex and is now in a radiative stage while the other part of blastwave has been expanding into the low-density region in the superbubble. This scenario naturally explains the lack of the X-ray emission related to the ejecta and blastwave. The SNR-MC interaction also favours a hadronic contribution to the {\gamma}-ray emission from the CTB 87 region.

We present extensive spectroscopic observations of supernova remnant (SNR) S147 collected with the Large sky Area Multi-Object fiber Spectroscopic Telescope (LAMOST). The spectra were carefully sky-subtracted taking into account the complex filamentary structure of S147. We have utilized all available LAMOST spectra toward S147, including sky and stellar spectra. By measuring the prominent optical emission lines including H{\alpha}, [NII]{\lambda}6584, and [SII]{\lambda}{\lambda}6717,6731, we present maps of radial velocity and line intensity ratio covering the whole nebula of S147 with unprecedented detail. The maps spatially correlated well with the complex filamentary structure of S147. For the central 2 deg of S147, the radial velocity varies from -100 to 100 km/s and peaks between ~ 0 and 10 km/s. The intensity ratios of H{\alpha}/[SII]{\lambda}{\lambda}6717,6731, [SII]{\lambda}6717/{\lambda}6731 and H{\alpha}/[NII]{\lambda}6584 peak at about 0.77, 1.35 and 1.48, respectively, with a scatter of 0.17, 0.19 and 0.37, respectively. The intensity ratios are consistent with the literature values. However, the range of variations of line intensity ratios estimated here and representative of the whole nebula, are larger than previously estimated.

We present a three-dimensional modeling of the Milky Way dust distribution by fitting the value-added star catalog of LAMOST spectral survey. The global dust distribution can be described by an exponential disk with scale-length of 3,192 pc and scale height of 103 pc. In this modeling, the Sun is located above the dust disk with a vertical distance of 23 pc. Besides the global smooth structure, two substructures around the solar position are also identified. The one located at $150^{\circ}<l<200^{\circ}$ and $-5^{\circ}<b<-30^{\circ}$ is consistent with the Gould Belt model of \citet{Gontcharov2009}, and the other one located at $140^{\circ}<l<165^{\circ}$ and $0^{\circ}<b<15^{\circ}$ is associated with the Camelopardalis molecular clouds.

Stellar systems composed of single, double, triple or high-order systems are rightfully regarded as the fundamental building blocks of the Milky Way. Binary stars play an important role in formation and evolution of the Galaxy. Through comparing the radial velocity variations from multi-epoch observations, we analyze the binary fraction of dwarf stars observed with the LAMOST. Effects of different model assumptions such as orbital period distributions on the estimate of binary fractions, are investigated. The results based on log-normal distribution of orbital periods reproduce the previous complete analyses better than the power-law distribution. We find that the binary fraction increases with $T_{\rm eff}$ and decreases with [Fe/H]. We first investigate the relation between $\alpha$-elements and binary fraction in such a large sample as the LAMOST. The old stars with high [$\alpha$/Fe] dominate higher binary fraction than young stars with low [$\alpha$/Fe]. At the same mass, former forming stars possess a higher binary fraction than newly forming ones, which may be related with the evolution of the Galaxy.

We present a detailed analysis of the selection function of the LAMOST Spectroscopic Survey of the Galactic Anti-centre (LSS-GAC). LSS-GAC was designed to obtain low resolution optical spectra for a sample of more than 3 million stars in the Galactic anti-centre. The second release of value-added catalogues of the LSS-GAC (LSS-GAC DR2) contains stellar parameters, including radial velocity, atmospheric parameters, elemental abundances and absolute magnitudes deduced from 1.8 million spectra of 1.4 million unique stars targeted by the LSS-GAC between 2011 and 2014. For many studies using this database, such as those investigating the chemodynamical structure of the Milky Way, a detailed understanding of the selection function of the survey is indispensable. In this paper, we describe how the selection function of the LSS-GAC can be evaluated to sufficient detail and provide selection function corrections for all spectroscopic measurements with reliable parameters released in LSS-GAC DR2. The results, to be released as new entries in the LSS-GAC value-added catalogues, can be used to correct the selection effects of the catalogue for scientific studies of various purposes.

In this work, we present a catalog of 2651 carbon stars from the fourth Data Release (DR4) of the Large Sky Area Multi-Object Fiber Spectroscopy Telescope (LAMOST). Using an efficient machine-learning algorithm, we find out these stars from more than seven million spectra. As a by-product, 17 carbon-enhanced metal-poor (CEMP) turnoff star candidates are also reported in this paper, and they are preliminarily identified by their atmospheric parameters. Except for 176 stars that could not be given spectral types, we classify the other 2475 carbon stars into five subtypes including 864 C-H, 226 C-R, 400 C-J, 266 C-N, and 719 barium stars based on a series of spectral features. Furthermore, we divide the C-J stars into three subtypes of CJ( H), C-J(R), C-J(N), and about 90% of them are cool N-type stars as expected from previous literature. Beside spectroscopic classification, we also match these carbon stars to multiple broadband photometries. Using ultraviolet photometry data, we find that 25 carbon stars have FUV detections and they are likely to be in binary systems with compact white dwarf companions.

The Xuyi Schmidt Telescope Photometric Survey of the Galactic Anti-center (XSTPS-GAC) is a photometric sky survey that covers nearly 6 000 deg^2 towards Galactic anti-center in g r i bands. Half of its survey field locates on the Galactic Anti-center disk, which makes XSTPS-GAC highly suitable for searching new open clusters in the GAC region. In this paper, we report new open cluster candidates discovered in this survey, as well as properties of these open cluster candidates, such as age, distance and reddening, derived by isochrone fitting in the color-magnitude diagram (CMD). These open cluster candidates are stellar density peaks detected in the star density maps by applying the method from Koposov et al. (2008). Each candidate is inspected on its true color image composed from XSTPS-GAC three band images. Then its CMD is checked, in order to identify whether the central region stars have a clear isochrone-like trend differing from the background stars. The parameters derived from isochrone fitting for these candidates are mainly based on three band photometry of XSTPS-GAC. Meanwhile, when these new candidates are able to be seen clearly on 2MASS, their parameters are also derived based on 2MASS (J-H, J) CMD. Finally, there are 320 known open clusters rediscovered and 24 new open cluster candidates discovered in this work. Further more, the parameters of these new candidates, as well as another 11 known recovered open clusters, are properly determined for the first time.

Fiber bursts are a type of fine structure that is often superposed on type IV radio continuum emission during solar flares. Although studied for many decades, its physical exciter, emission mechanism, and association with the flare energy release remain unclear, partly due to the lack of simultaneous imaging observations. We report the first dynamic spectroscopic imaging observations of decimetric fiber bursts, which occurred during the rise phase of a long-duration eruptive flare on 2012 March 3, as obtained by the Karl G. Jansky Very Large Array in 1--2 GHz. Our results show that the fiber sources are located near and above one footpoint of the flare loops. The fiber source and the background continuum source are found to be cospatial and share the same morphology. It is likely that they are associated with nonthermal electrons trapped in the converging magnetic fields near the footpoint, as supported by a persistent coronal hard X-ray source present during the flare rise phase. We analyze three groups of fiber bursts in detail with dynamic imaging spectroscopy and obtain their mean frequency-dependent centroid trajectories in projection. By using a barometric density model and magnetic field based on a potential-field extrapolation, we further reconstruct the 3-D source trajectories of fiber bursts, for comparison with expectations from the whistler wave model and two MHD-based models. We conclude that the observed fiber burst properties are consistent with an exciter moving at the propagation velocity expected for whistler waves, or models that posit similar exciter velocities.

Linear perturbation theory is a powerful toolkit for studying black hole spacetimes. However, the perturbation equations are hard to solve unless we can use separation of variables. In the Kerr spacetime, metric perturbations do not separate, but curvature perturbations do. The cost of curvature perturbations is a very complicated metric-reconstruction procedure. This procedure can be avoided using a symmetry-adapted choice of basis functions in highly symmetric spacetimes, such as near-horizon extremal Kerr. In this paper, we focus on this spacetime, and (i) construct the symmetry-adapted basis functions; (ii) show their orthogonality; and (iii) show that they lead to separation of variables of the scalar, Maxwell, and metric perturbation equations. This separation turns the system of partial differential equations into one of ordinary differential equations over a compact domain, the polar angle.

We present a three dimensional (3D) extinction analysis in the region toward the supernova remnant (SNR) S147 (G180.0-1.7) using multi-band photometric data from the Xuyi Schmidt Telescope Photometric Survey of the Galactic Anticentre (XSTPS-GAC), 2MASS and WISE. We isolate a previously unrecognised dust structure likely to be associated with SNR S147. The structure, which we term as "S147 dust cloud", is estimated to have a distance $d$ = 1.22 $\pm$ 0.21 kpc, consistent with the conjecture that S147 is associated with pulsar PSR J0538 + 2817. The cloud includes several dense clumps of relatively high extinction that locate on the radio shell of S147 and coincide spatially with the CO and gamma-ray emission features. We conclude that the usage of CO measurements to trace the SNR associated MCs is unavoidably limited by the detection threshold, dust depletion, and the difficulty of distance estimates in the outer Galaxy. 3D dust extinction mapping may provide a better way to identify and study SNR-MC interactions.

We report the discovery of two new unbound hypervelocity stars (HVSs) from the LAMOST spectroscopic surveys. They are respectively a B2V type star of ~ 7 M$_{\rm \odot}$ with a Galactic rest-frame radial velocity of 502 km/s at a Galactocentric radius of ~ 21 kpc and a B7V type star of ~ 4 M$_{\rm \odot}$ with a Galactic rest-frame radial velocity of 408 km/s at a Galactocentric radius of ~ 30 kpc. The origins of the two HVSs are not clear given their currently poorly measured proper motions. However, the future data releases of Gaia should provide proper motion measurements accurate enough to solve this problem. The ongoing LAMOST spectroscopic surveys are expected to yield more HVSs to form a statistical sample, providing vital constraint on understanding the nature of HVSs and their ejection mechanisms.

We present estimates of stellar age and mass for 0.93 million Galactic disk main sequence turn-off and sub-giant stars from the LAMOST Galactic Spectroscopic Surveys. The ages and masses are determined by matching with stellar isochrones using Bayesian algorithm, utilizing effective temperature $T_{\rm eff}$, absolute magnitude ${\rm M}_V$, metallicity [Fe/H] and $\alpha$-element to iron abundance ratio [$\alpha$/Fe] deduced from the LAMOST spectra. Extensive examinations suggest the age and mass estimates are robust. The overall sample stars have a median error of 34 per cent for the age estimates, and half of the stars older than 2\,Gyr have age uncertainties of only 20--30 per cent. Median error for the mass estimates of the whole sample stars is $\sim8$ per cent. The huge dataset demonstrates good correlations among stellar age, [Fe/H] ([$\alpha$/H]) and [$\alpha$/Fe]. Particularly, double sequence features are revealed in the both the age--[$\alpha$/Fe] and age--[Fe/H]([$\alpha$/H]) spaces. In the [Fe/H]--[$\alpha$/Fe] space, stars of 8--10\,Gyr exhibit both the thin and thick disk sequences, while younger (older) stars show only the thin (thick) disk sequence, indicating that the thin disk became prominent 8--10\,Gyr ago, while the thick disk formed earlier and almost quenched 8\,Gyr ago. Stellar ages exhibit positive vertical and negative radial gradients across the disk, and the outer disk of $R\gtrsim$\,9\,kpc exhibits a strong flare in stellar age distribution.

The Atacama Large Millimeter-submillimeter Array (ALMA) radio telescope has commenced science observations of the Sun starting in late 2016. Since the Sun is much larger than the field of view of individual ALMA dishes, the ALMA interferometer is unable to measure the background level of solar emission when observing the solar disk. The absolute temperature scale is a critical measurement for much of ALMA solar science, including the understanding of energy transfer through the solar atmosphere, the properties of prominences, and the study of shock heating in the chromosphere. In order to provide an absolute temperature scale, ALMA solar observing will take advantage of the remarkable fast-scanning capabilities of the ALMA 12m dishes to make single-dish maps of the full Sun. This article reports on the results of an extensive commissioning effort to optimize the mapping procedure, and it describes the nature of the resulting data. Amplitude calibration is discussed in detail: a path that utilizes the two loads in the ALMA calibration system as well as sky measurements is described and applied to commissioning data. Inspection of a large number of single-dish datasets shows significant variation in the resulting temperatures, and based on the temperature distributions we derive quiet-Sun values at disk center of 7300 K at lambda=3 mm and 5900 K at lambda=1.3 mm. These values have statistical uncertainties of order 100 K, but systematic uncertainties in the temperature scale that may be significantly larger. Example images are presented from two periods with very different levels of solar activity. At a resolution of order 25 arcsec, the 1.3 mm wavelength images show temperatures on the disk that vary over about a 2000 K range.

Observations of the Sun at millimeter and submillimeter wavelengths offer a unique probe into the structure, dynamics, and heating of the chromosphere; the structure of sunspots; the formation and eruption of prominences and filaments; and energetic phenomena such as jets and flares. High-resolution observations of the Sun at millimeter and submillimeter wavelengths are challenging due to the intense, extended, low- contrast, and dynamic nature of emission from the quiet Sun, and the extremely intense and variable nature of emissions associated with energetic phenomena. The Atacama Large Millimeter/submillimeter Array (ALMA) was designed with solar observations in mind. The requirements for solar observations are significantly different from observations of sidereal sources and special measures are necessary to successfully carry out this type of observations. We describe the commissioning efforts that enable the use of two frequency bands, the 3 mm band (Band 3) and the 1.25 mm band (Band 6), for continuum interferometric-imaging observations of the Sun with ALMA. Examples of high-resolution synthesized images obtained using the newly commissioned modes during the solar commissioning campaign held in December 2015 are presented. Although only 30 of the eventual 66 ALMA antennas were used for the campaign, the solar images synthesized from the ALMA commissioning data reveal new features of the solar atmosphere that demonstrate the potential power of ALMA solar observations. The ongoing expansion of ALMA and solar-commissioning efforts will continue to enable new and unique solar observing capabilities.

Recently Graphics Processing Units (GPUs) have been used to speed up very CPU-intensive gravitational microlensing simulations. In this work, we use the Xeon Phi coprocessor to accelerate such simulations and compare its performance on a microlensing code with that of NVIDIA's GPUs. For the selected set of parameters evaluated in our experiment, we find that the speedup by Intel's Knights Corner coprocessor is comparable to that by NVIDIA's Fermi family of GPUs with compute capability 2.0, but less significant than GPUs with higher compute capabilities such as the Kepler. However, the very recently released second generation Xeon Phi, Knights Landing, is about 5.8 times faster than the Knights Corner, and about 2.9 times faster than the Kepler GPU used in our simulations. We conclude that the Xeon Phi is a very promising alternative to GPUs for modern high performance microlensing simulations.

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.

We present an X-ray photometric analysis of six gravitationally lensed quasars, with observation campaigns spanning from 5 to 14 years, measuring the total (0.83 - 21.8 keV restframe), soft (0.83 - 3.6 keV), and hard (3.6 - 21.8 keV) band image flux ratios or each epoch. Using the ratios of the model-predicted macro-magnifications as baselines, we build differential microlensing light curves and obtain joint likelihood functions for the average X-ray emission region sizes. Our analysis yields a Probability Distribution Function for the average half-light radius of the X-Ray emission region in the sample that peaks slightly above 1 gravitational radius and with nearly indistinguishable 68% confidence (one-sided) upper limits of 17.8 and 18.9 gravitational radii for the soft and hard X-ray emitting regions, assuming a mean stellar mass of 0.3 solar masses. We see hints of energy dependent microlensing between the soft and hard bands in two of the objects. In a separate analysis on the root-mean-square of the microlensing variability, we find significant differences between the soft and hard bands but the sign of the difference is not consistent across the sample. This suggests the existence of some kind of spatial structure to the X-ray emission in an otherwise extremely compact source. We also discover a correlation between the RMS microlensing variability and the average microlensing amplitude.

How impulsive magnetic energy release leads to solar eruptions and how those eruptions are energized and evolve are vital unsolved problems in Heliophysics. The standard model for solar eruptions summarizes our current understanding of these events. Magnetic energy in the corona is released through drastic restructuring of the magnetic field via reconnection. Electrons and ions are then accelerated by poorly understood processes. Theories include contracting loops, merging magnetic islands, stochastic acceleration, and turbulence at shocks, among others. Although this basic model is well established, the fundamental physics is poorly understood. HXR observations using grazing-incidence focusing optics can now probe all of the key regions of the standard model. These include two above-the-looptop (ALT) sources which bookend the reconnection region and are likely the sites of particle acceleration and direct heating. The science achievable by a direct HXR imaging instrument can be summarized by the following science questions and objectives which are some of the most outstanding issues in solar physics (1) How are particles accelerated at the Sun? (1a) Where are electrons accelerated and on what time scales? (1b) What fraction of electrons is accelerated out of the ambient medium? (2) How does magnetic energy release on the Sun lead to flares and eruptions? A Focusing Optics X-ray Solar Imager (FOXSI) instrument, which can be built now using proven technology and at modest cost, would enable revolutionary advancements in our understanding of impulsive magnetic energy release and particle acceleration, a process which is known to occur at the Sun but also throughout the Universe.

Application of fitting techniques to obtain physical parameters---such as ages, metallicities, and $\alpha$-element to iron ratios---of stellar populations is an important approach to understand the nature of both galaxies and globular clusters (GCs). In fact, fitting methods based on different underlying models may yield different results, and with varying precision. In this paper, we have selected 22 confirmed M31 GCs for which we do not have access to previously known spectroscopic metallicities. Most are located at approximately one degree (in projection) from the galactic center. We performed spectroscopic observations with the 6.5 m MMT telescope, equipped with its Red Channel Spectrograph. Lick/IDS absorption-line indices, radial velocities, ages, and metallicities were derived based on the $\rm EZ\_Ages$ stellar population parameter calculator. We also applied full spectral fitting with the ULySS code to constrain the parameters of our sample star clusters. In addition, we performed $\chi^2_{\rm min}$ fitting of the clusters' Lick/IDS indices with different models, including the Bruzual & Charlot models (adopting Chabrier or Salpeter stellar initial mass functions and 1994 or 2000 Padova stellar evolutionary tracks), the GALEV, and the Thomas et al. models. For comparison, we collected their $UVBRIJK$ photometry from the Revised Bologna Catalogue (v.5) to obtain and fit the GCs' spectral energy distributions (SEDs). Finally, we performed fits using a combination of Lick/IDS indices and SEDs. The latter results are more reliable and the associated error bars become significantly smaller than those resulting from either our Lick/IDS indices-only or our SED-only fits.

Accurate determination of stellar atmospheric parameters and elemental abundances is crucial for Galactic archeology via large-scale spectroscopic surveys. In this paper, we estimate stellar atmospheric parameters -- effective temperature T_{\rm eff}, surface gravity log g and metallicity [Fe/H], absolute magnitudes M_V and M_{Ks}, {\alpha}-element to metal (and iron) abundance ratio [{\alpha}/M] (and [{\alpha}/Fe]), as well as carbon and nitrogen abundances [C/H] and [N/H] from the LAMOST spectra with amultivariate regressionmethod based on kernel-based principal component analysis, using stars in common with other surveys (Hipparcos, Kepler, APOGEE) as training data sets. Both internal and external examinations indicate that given a spectral signal-to-noise ratio (SNR) better than 50, our method is capable of delivering stellar parameters with a precision of ~100K for Teff, ~0.1 dex for log g, 0.3 -- 0.4mag for M_V and M_{Ks}, 0.1 dex for [Fe/H], [C/H] and [N/H], and better than 0.05 dex for [{\alpha}/M] ([{\alpha}/Fe]). The results are satisfactory even for a spectral SNR of 20. The work presents first determinations of [C/H] and [N/H] abundances from a vast data set of LAMOST, and, to our knowledge, the first reported implementation of absolute magnitude estimation directly based on the observed spectra. The derived stellar parameters for millions of stars from the LAMOST surveys will be publicly available in the form of value-added catalogues.

Photometric data from the Xuyi Schmidt Telescope Photometric Survey of the Galactic Anticentre (XSTPS-GAC) and the Sloan Digital Sky Survey (SDSS) are used to derive the global structure parameters of the smooth components of the Milky Way. The data, which cover nearly 11,000 deg$^2$ sky area and the full range of Galactic latitude, allow us to construct a globally representative Galactic model. The number density distribution of Galactic halo stars is fitted with an oblate spheroid that decays by power law. The best-fit yields an axis ratio and a power law index $\kappa=0.65$ and $p=2.79$, respectively. The $r$-band differential star counts of three dwarf samples are then fitted with a Galactic model. The best-fit model yielded by a Markov Chain Monte Carlo analysis has thin and thick disk scale heights and lengths of $H_{1}=$ 322\,pc and $L_{1}=$2343\,pc, $H_{2}=$794\,pc and $L_{2}=$3638\,pc, a local thick-to-thin disk density ratio of $f_2=$11\,per\,cent, and a local density ratio of the oblate halo to the thin disk of $f_h=$0.16\,per\,cent. The measured star count distribution, which is in good agreement with the above model for most of the sky area, shows a number of statistically significant large scale overdensities, including some of the previously known substructures, such as the Virgo overdensity and the so-called "north near structure", and a new feature between 150\degr $< l < $ 240\degr~and $-1$5\degr $< b < $ $-$5\degr, at an estimated distance between 1.0 and 1.5\,kpc. The Galactic North-South asymmetry in the anticentre is even stronger than previously thought.

The rotation curve (RC) of the Milky Way out to $\sim$ 100 kpc has been constructed using $\sim$ 16,000 primary red clump giants (PRCGs) in the outer disk selected from the LSS-GAC and the SDSS-III/APOGEE survey, combined with $\sim$ 5700 halo K giants (HKGs) selected from the SDSS/SEGUE survey. To derive the RC, the PRCG sample of the warm disc population and the HKG sample of halo stellar population are respectively analyzed using a kinematical model allowing for the asymmetric drift corrections and re-analyzed using the spherical Jeans equation along with measurements of the anisotropic parameter $\beta$ currently available. The typical uncertainties of RC derived from the PRCG and HKG samples are respectively 5-7 km/s and several tens km/s. We determine a circular velocity at the solar position, $V_c (R_0)$ = 240 $\pm$ 6 km/s and an azimuthal peculiar speed of the Sun, $V_{\odot}$ = 12.1 $\pm$ 7.6 km/s, both in good agreement with the previous determinations. The newly constructed RC has a generally flat value of 240 km/s within a Galactocentric distance $r$ of 25 kpc and then decreases steadily to 150 km/s at $r$ $\sim$ 100 kpc. On top of this overall trend, the RC exhibits two prominent localized dips, one at $r$ $\sim$ 11 kpc and another at $r$ $\sim$ 19 kpc. From the newly constructed RC, combined with other constraints, we have built a parametrized mass model for the Galaxy, yielding a virial mass of the Milky Way's dark matter halo of $0.90^{+0.07}_{-0.08} \times 10^{12}$ ${\rm M}_{\odot}$ and a local dark matter density, $\rho_{\rm \odot, dm} = 0.32^{+0.02}_{-0.02}$ GeV cm$^{-3}$.

We present results from the the first campaign of dedicated solar observations undertaken by the \textit{Nuclear Spectroscopic Telescope ARray} ({\em NuSTAR}) hard X-ray telescope. Designed as an astrophysics mission, {\em NuSTAR} nonetheless has the capability of directly imaging the Sun at hard X-ray energies ($>$3~keV) with an increase in sensitivity of at least two magnitude compared to current non-focusing telescopes. In this paper we describe the scientific areas where \textit{NuSTAR} will make major improvements on existing solar measurements. We report on the techniques used to observe the Sun with \textit{NuSTAR}, their limitations and complications, and the procedures developed to optimize solar data quality derived from our experience with the initial solar observations. These first observations are briefly described, including the measurement of the Fe K-shell lines in a decaying X-class flare, hard X-ray emission from high in the solar corona, and full-disk hard X-ray images of the Sun.

Using a Bayesian technology we derived distances and extinctions for over 100,000 red giant stars observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey by taking into account spectroscopic constraints from the APOGEE stellar parameters and photometric constraints from 2MASS, as well as a prior knowledge on the Milky Way. Derived distances are compared with those from four other independent methods, the Hipparcos parallaxes, star clusters, APOGEE red clump stars, and asteroseismic distances from APOKASC (Rodrigues et al. 2014) and SAGA Catalogues (Casagrande et al. 2014). These comparisons covers four orders of magnitude in the distance scale from 0.02 kpc to 20 kpc. The results show that our distances agree very well with those from other methods: the mean relative difference between our Bayesian distances and those derived from other methods ranges from -4.2% to +3.6%, and the dispersion ranges from 15% to 25%. The extinctions toward all stars are also derived and compared with those from several other independent methods: the Rayleigh-Jeans Color Excess (RJCE) method, Gonzalez's two-dimensional extinction map, as well as three-dimensional extinction maps and models. The comparisons reveal that, overall, estimated extinctions agree very well, but RJCE tends to overestimate extinctions for cool stars and objects with low logg.

An observation from the Interface Region Imaging Spectrograph reveals coherent oscillations in the loops of an M1.6 flare on 2015 March 12. Both the intensity and Doppler shift of Fe~{\sc{xxi}}~1354.08\AA{}~show clear oscillations with a period of $\sim$25 seconds. Remarkably similar oscillations were also detected in the soft X-ray flux recorded by the Geostationary Operational Environmental Satellites (GOES). With an estimated phase speed of $\sim$2420~km~s$^{-1}$ and a derived electron density of at least 5.4$\times$10$^{10}$ cm$^{-3}$, the observed short-period oscillation is most likely the global fast sausage mode of a hot flare loop. We find a phase shift of $\sim$$\pi$/2 (1/4 period) between the Doppler shift oscillation and the intensity/GOES oscillations, which is consistent with a recent forward modeling study of the sausage mode. The observed oscillation requires a density contrast between the flare loop and coronal background of a factor $\geqslant$42. The estimated phase speed of the global mode provides an lower limit of the Alfv\'en speed outside the flare loop. We also find an increase of the oscillation period, which might be caused by the separation of the loop footpoints with time.

Jets present ubiquitously in both quiet and active regions on the Sun. They are widely believed to be driven by magnetic reconnection. A fan-spine structure has been frequently reported in some coronal jets and flares, regarded as a signature of ongoing magnetic reconnection in a topology consisting of a magnetic null connected by a fan-like separatrix surface and a spine. However, for small-scale chromospheric jets, clear evidence of such structures is rather rare, although they are implied in earlier works that show an inverted-Y-shaped feature. Here we report high-resolution (0$"$.16) observations of a small-scale chromospheric jet obtained by the New Solar Telescope (NST) using 10830~{\AA} filtergrams. Bi-directional flows were observed across the separatrix regions in the 10830~\AA{} images, suggesting that the jet was produced due to magnetic reconnection. At the base of the jet, a fan-spine structure was clearly resolved by the NST, including the spine and the fan-like surface, as well as the loops before and after the reconnection. A major part of this fan-spine structure, with the exception of its bright footpoints and part of the base arc, was invisible in the extreme ultraviolet and soft X-ray images (observed by the Atmosphere Imaging Assembly and the X-Ray Telescope, respectively), indicating that the reconnection occurred in the upper chromosphere. Our observations suggest that the evolution of this chromospheric jet is consistent with a two-step reconnection scenario proposed by \citet{Torok2009}

Solar flares - the most powerful explosions in the solar system - are also efficient particle accelerators, capable of energizing a large number of charged particles to relativistic speeds. A termination shock is often invoked in the standard model of solar flares as a possible driver for particle acceleration, yet its existence and role have remained controversial. We present observations of a solar flare termination shock and trace its morphology and dynamics using high-cadence radio imaging spectroscopy. We show that a disruption of the shock coincides with an abrupt reduction of the energetic electron population. The observed properties of the shock are well-reproduced by simulations. These results strongly suggest that a termination shock is responsible, at least in part, for accelerating energetic electrons in solar flares.

The effect of gravitational Faraday rotation was predicted in the 1950s, but there is currently no practical method for measuring this effect. Measuring this effect is important because it will provide new evidence for correctness of general relativity, in particular, in the strong field limit. We predict that the observed degree and angle of the X-ray polarization of a cosmologically distant quasar microlensed by the random star field in a foreground galaxy or cluster lens vary rapidly and concurrently with flux during caustic-crossing events using the first simulation of quasar X-ray microlensing polarization light curves. Therefore, it is possible to detect gravitational Faraday rotation by monitoring the X-ray polarization of gravitationally microlensed quasars. Detecting this effect will also confirm the strong gravity nature of quasar X-ray emission.

Asteroseismology allows for deriving precise values of surface gravity of stars. The accurate asteroseismic determinations now available for large number of stars in the Kepler fields can be used to check and calibrate surface gravities that are currently being obtained spectroscopically for a huge numbers of stars targeted by large-scale spectroscopic surveys, such as the on-going Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) Galactic survey. The LAMOST spectral surveys have obtained a large number of stellar spectra in the Kepler fields. Stellar atmospheric parameters of those stars have been determined with the LAMOST Stellar Parameter Pipeline at Peking University (LSP3), by template matching with the MILES empirical spectral library. In the current work, we compare surface gravities yielded by LSP3 with those of two asteroseismic samples - the largest sample from Huber et al. (2014) and the most accurate sample from Hekker et al. (2012, 2013). We find that LSP3 surface gravities are in good agreement with asteroseismic values of Hekker et al. (2012, 2013), with a dispersion of about 0.2 dex. Except for a few cases, asteroseismic surface gravities of Huber et al. (2014) and LSP3 spectroscopic values agree for a wide range of surface gravities. However, some patterns of differences can be identified upon close inspection. Potential ways to further improve the LSP3 spectroscopic estimation of stellar atmospheric parameters in the near future are briefly discussed. The effects of effective temperature and metallicity on asteroseismic determinations of surface gravities for giant stars are also discussed.

We review results from our monitoring observations of several lensed quasars performed in the optical, UV, and X-ray bands. Modeling of the multi-wavelength light curves provides constraints on the extent of the optical, UV, and X-ray emission regions. One of the important results of our analysis is that the optical sizes as inferred from the microlensing analysis are significantly larger than those predicted by the theoretical-thin-disk estimate. In a few cases we also constrain the slope of the size-wavelength relation. Our size constraints of the soft and hard X-ray emission regions of quasars indicate that in some objects of our sample the hard X-ray emission region is more compact than the soft and in others the soft emission region is smaller. This difference may be the result of the relative strengths of the disk-reflected (harder and extended) versus corona-direct (softer and compact) components in the quasars of our sample. Finally, we present the analysis of several strong microlensing events where we detect an evolution of the relativistic Fe line profile as the magnification caustic traverses the accretion disk. These caustic crossings are used to provide constraints on the innermost stable circular orbit (ISCO) radius and the accretion disk inclination angle of the black hole in quasar RX J1131-1231.

We present empirical metallicity-dependent calibrations of effective temperature against colours for dwarfs of luminosity classes IV and V and for giants of luminosity classes II and III, based on a collection from the literature of about two hundred nearby stars with direct effective temperature measurements of better than 2.5 per cent. The calibrations are valid for an effective temperature range 3,100 - 10,000 K for dwarfs of spectral types M5 to A0 and 3,100 - 5,700 K for giants of spectral types K5 to G5. A total of twenty-one colours for dwarfs and eighteen colours for giants of bands of four photometric systems, i.e. the Johnson ($UBVR_{\rm J}I_{\rm J}JHK$), the Cousins ($R_{\rm C}I_{\rm C}$), the Sloan Digital Sky Survey (SDSS, $gr$) and the Two Micron All Sky Survey (2MASS, $JHK_{\rm s}$), have been calibrated. Restricted by the metallicity range of the current sample, the calibrations are mainly applicable for disk stars ([Fe/H]$\,\gtrsim\,-1.0$). The normalized percentage residuals of the calibrations are typically 2.0 and 1.5 per cent for dwarfs and giants, respectively. Some systematic discrepancies at various levels are found between the current scales and those available in the literature (e.g. those based on the infrared flux method IRFM or spectroscopy). Based on the current calibrations, we have re-determined the colours of the Sun. We have also investigated the systematic errors in effective temperatures yielded by the current on-going large scale low- to intermediate-resolution stellar spectroscopic surveys. We show that the calibration of colour ($g-K_{\rm s}$) presented in the current work provides an invaluable tool for the estimation of stellar effective temperature for those on-going or upcoming surveys.

With observations from the Interface Region Imaging Spectrograph (IRIS), we track the complete evolution of $\sim$11 MK evaporation flows in an M1.1 flare on 2014 September 6 and an X1.6 flare on 2014 September 10. These hot flows, as indicated by the blueshifted Fe~{\sc{xxi}}~1354.08\AA{}~line, evolve smoothly with a velocity decreasing exponentially from $\sim$200~km~s$^{-1}$ to almost stationary within a few minutes. We find a good correlation between the flow velocity and energy deposition rate as represented by the hard X-Ray flux observed with the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI), or time derivative of the soft X-Ray flux observed with the Geostationary Operational Environmental Satellites (GOES) and the HINODE X-ray Telescope (XRT), which is in general agreement with models of nonthermal electron heating. The maximum blue shift of Fe~{\sc{xxi}}~appears approximately at the same time as or slightly after the impulsive enhancement of the ultraviolet continuum and the Mg~{\sc{ii}}~2798.8\AA{}~line emission, demonstrating that the evaporation flow is closely related to heating of the lower chromosphere. Finally, while the hot Fe~{\sc{xxi}}~1354.08\AA{} line is entirely blueshifted with no obvious rest component, cool chromospheric and transition region lines like Si~{\sc{iv}}~1402.77\AA{} are often not entirely redshifted but just reveal an obvious red wing enhancement at the ribbons, suggesting that the speed of chromospheric condensation might be larger than previously thought.

Using a sample of over 70, 000 red clump (RC) stars with $5$-$10$% distance accuracy selected from the LAMOST Spectroscopic Survey of the Galactic Anti-center (LSS-GAC), we study the radial and vertical gradients of the Galactic disk(s) mainly in the anti-center direction, covering a significant disk volume of projected Galactocentric radius $7 \leq R_{\rm GC} \leq 14$ kpc and height from the Galactic midplane $0 \leq |Z| \leq 3$ kpc. Our analysis shows that both the radial and vertical metallicity gradients are negative across much of the disk volume probed, and exhibit significant spatial variations. Near the solar circle ($7 \leq R_{\rm GC} \leq 11.5$ kpc), the radial gradient has a moderately steep, negative slope of $-0.08$ dex kpc$^{-1}$ near the midplane ($|Z| < 0.1$ kpc), and the slope flattens with increasing $|Z|$. In the outer disk ($11.5 < R_{\rm GC} \leq 14$ kpc), the radial gradients have an essentially constant, much less steep slope of $-0.01$ dex kpc $^{-1}$ at all heights above the plane, suggesting that the outer disk may have experienced an evolution path different from that of the inner disk. The vertical gradients are found to flatten largely with increasing $R_{\rm GC}$. However, the vertical gradient of the lower disk ($0 \leq |Z| \leq 1$ kpc) is found to flatten with $R_{\rm GC}$ quicker than that of the upper disk ($1 < |Z| \leq 3$ kpc). Our results should provide strong constraints on the theory of disk formation and evolution, as well as the underlying physical processes that shape the disk (e.g. gas flows, radial migration, internal and external perturbations).

We report a detailed investigation of the bulk motions of the nearby Galactic stellar disk, based on three samples selected from the LSS-GAC DR2: a global sample containing 0.57 million FGK dwarfs out to $\sim$ 2 kpc, a local subset of the global sample consisting $\sim$ 5,400 stars within 150 pc, and an anti-center sample containing $\sim$ 4,400 AFGK dwarfs and red clump stars within windows of a few degree wide centered on the Galactic anti-center. The global sample is used to construct a three-dimensional map of bulk motions of the Galactic disk from the solar vicinity out to $\sim$ 2 kpc with a spatial resolution of $\sim$ 250 pc. Typical values of the radial and vertical components of bulk motion range from $-$15 km s$^{-1}$ to 15 km s$^{-1}$, while the lag behind the circular speed dominates the azimuthal component by up to $\sim$ 15 km s$^{-1}$. The map reveals spatially coherent, kpc-scale stellar flows in the disk, with typical velocities of a few tens km s$^{-1}$. Bending- and breathing-mode perturbations are clearly visible, and vary smoothly across the disk plane. Our data also reveal higher-order perturbations, such as breaks and ripples, in the profiles of vertical motion versus height. From the local sample, we find that stars of different populations exhibit very different patterns of bulk motion. Finally, the anti-center sample reveals a number of peaks in stellar number density in the line-of-sight velocity versus distance distribution, with the nearer ones apparently related to the known moving groups. The "velocity bifurcation" reported by Liu et al. (2012) at Galactocentric radii 10--11 kpc is confirmed. However, just beyond this distance, our data also reveal a new triple-peaked structure.

We present a catalog of 908 objects observed with the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) in the vicinity fields of M31 and M33, targeted as globular clusters (GCs) and candidates. The targets include known GCs and candidates selected from the literature, as well as new candidates selected from the Sloan Digital Sky Survey (SDSS). Analysis shows that 356 of them are likely GCs of various degree of confidence, while the remaining ones turn out to be background galaxies and quasars, stars and HII regions in M31 or foreground Galactic stars. The 356 likely GCs include 298 bona fide GCs and 26 candidates known in the literature. Three candidates selected from the Revised Bologna Catalog of M31 GCs and candidates (RBC) and one possible cluster from Johnson et al. are confirmed to be bona fide clusters. We search for new GCs in the halo of the M31 amongst the new candidates selected from the SDSS photometry. Based on radial velocities yielded by LAMSOT spectra and visual examination of the SDSS images, we find 28 objects, 5 bona fide and 23 likely GCs. Amongst the five bona fide GCs, three have been recently discovered independently by others, the remaining 25 are our new identifications, including two bona fide ones. The new identified objects fall at projected distances ranging from 13 to 265 kpc from M31. Of the two newly discovered bona fide GCs, one is located near M33, probably a GC belonging to M33. The other bona fide GC falls on the Giant Stream with a projected distance of 78 kpc from M31. Of the 23 newly identified likely GCs, one has a projected distance of about 265 kpc from M31 and could be an intergalactic cluster.

We use 297 042 main sequence turn-off stars selected from the LSS-GAC to determine the radial and vertical gradients of stellar metallicity of the Galactic disk in the anti-center direction. We determine ages of those turn-off stars by isochrone fitting and measure the temporal variations of metallicity gradients. Our results show that the gradients, both in the radial and vertical directions, exhibit significant spatial and temporal variations. The radial gradients yielded by stars of oldest ages (>11 Gyr) are essentially zero at all heights from the disk midplane, while those given by younger stars are always negative. The vertical gradients deduced from stars of oldest ages (>11Gyr) are negative and show only very weak variations with the Galactocentric distance in the disk plane, $R$, while those yielded by younger stars show strong variations with $R$. After being essentially flat at the earliest epochs of disk formation, the radial gradients steepen as age decreases, reaching a maxima (steepest) at age 7-8 Gyr, and then they flatten again. Similar temporal trends are also found for the vertical gradients. We infer that the assemblage of the Milky Way disk may have experienced at least two distinct phases. The earlier phase is probably related to a slow, pressure-supported collapse of gas, when the gas settles down to the disk mainly in the vertical direction. In the later phase, there are significant radial flows of gas in the disk, and the rate of gas inflow near the solar neighborhood reaches a maximum around a lookback time of 7-8 Gyr. The transition of the two phases occurs around a lookback time between 8 and 11 Gyr. The two phases may be responsible for the formation of the Milky Way thick and thin disks, respectively. And, as a consequence, we recommend that stellar age is a natural, physical criterion to distinguish thin and thick disk stars. ... (abridged)

We analyze the optical, UV, and X-ray microlensing variability of the lensed quasar SDSS J0924+0219 using six epochs of Chandra data in two energy bands (spanning 0.4-8.0 keV, or 1-20 keV in the quasar rest frame), 10 epochs of F275W (rest-frame 1089A) Hubble Space Telescope data, and high-cadence R-band (rest-frame 2770A) monitoring spanning eleven years. Our joint analysis provides robust constraints on the extent of the X-ray continuum emission region and the projected area of the accretion disk. The best-fit half-light radius of the soft X-ray continuum emission region is between 5x10^13 and 10^15 cm, and we find an upper limit of 10^15 cm for the hard X-rays. The best-fit soft-band size is about 13 times smaller than the optical size, and roughly 7 GM_BH/c^2 for a 2.8x10^8 M_sol black hole, similar to the results for other systems. We find that the UV emitting region falls in between the optical and X-ray emitting regions at 10^14 cm < r_1/2,UV < 3x10^15 cm. Finally, the optical size is significantly larger, by 1.5*sigma, than the theoretical thin-disk estimate based on the observed, magnification-corrected I-band flux, suggesting a shallower temperature profile than expected for a standard disk.

Active galactic nuclei (AGNs) and quasars are important astrophysical objects to understand. Recently, microlensing observations have constrained the size of the quasar X-ray emission region to be of the order of 10 gravitational radii of the central supermassive black hole. For distances within a few gravitational radii, light paths are strongly bent by the strong gravity field of the central black hole. If the central black hole has nonzero angular momentum (spin), a photon's polarization plane will be rotated by the gravitational Faraday effect. The observed X-ray flux and polarization will then be influenced significantly by the strong gravity field near the source. Consequently, linear gravitational lensing theory is inadequate for such extreme circumstances. We present simple algorithms computing strong lensing effects of Kerr black holes, including effects on polarization. Our algorithms are realized in a program "KERTAP" in two versions: MATLAB and Python. The key ingredients of KERTAP are: a graphic user interface, a {\it backward} ray-tracing algorithm, a polarization propagator dealing with gravitational Faraday rotation, and algorithms computing observables such as flux magnification and polarization angles. Our algorithms can be easily realized in other programming languages such as FORTRAN, C, and C++. The MATLAB version of KERTAP is parallelized using the MATLAB Parallel Computing Toolbox and the Distributed Computing Server. The Python code was sped up using Cython and supports full implementation of MPI using 'mpi4py' package. As an example, we investigate the inclination angle dependence of the observed polarization and the strong lensing magnification of AGN X-ray emission. We conclude that it is possible to perform complex numerical-relativity-related computations using interpreted languages such as MATLAB and Python.

The photometry profile of the integrated Sachs-Wolfe (ISW) effect recently obtained by the Planck consortium by stacking patches of Cosmic Microwave Background (CMB) sky maps around a large number of cosmic voids, contains a cold ring at about half the void's effective radius surrounded by a hot ring near the void's boundary. The source of the temperature structure is assumed to be the ISW effect but the exact cause of the ringed structure is not currently well understood, particularly the outer hot ring. Numerical simulations have suggested that hot/cold ring structures can be produced by motions associated with nonlinear growths of cosmic structures whose gravitational potentials produce the ISW effect. We have recently developed the embedded lens theory and the Fermat potential formalism which can be used to model the ISW effect caused by intervening individual lens inhomogeneities evolving arbitrarily. This theory only requires knowledge of the void's projected mass profile as a function of the passing CMB photons' impact radius and the rate of change of that mass distribution at passage. We present two simple embedded void lens models with evolving mass densities and investigate the ISW effect caused by these lenses. Both models posses expanding mass shells which produce hot rings around central cold regions, consistent with the recent observations. By adding a small over-density at the void's center we can produce the slight positive temperature excess hinted at in Planck's photometric results. We conclude that the embedded lens theory and the Fermat potential formalism is well suited for modeling the ISW effect.

We investigate the evolution of the power law k-essence field in FRWL spacetime. The autonomous dynamical system and critical points are obtained. The corresponding cosmological parameters, such as $\Omega _{\phi }$ and $w_{\phi }$, are calculated at these critical points. We find it is possible to achieve an equation of state crossing through $-1$ for k-essence field. The results we obtained indicate that the power law k-essence dark energy model can be compatible with observations.

We present a simple gravitational lens model to illustrate the ease of using the embedded lensing theory when studying cosmic voids. It confirms the previously used repulsive lensing models for deep voids. We start by estimating magnitude fluctuations and weak lensing shears of background sources lensed by large voids. We find that sources behind large ($\sim$$90\,\rm Mpc$) and deep voids (density contrast about $-0.9$) can be magnified or demagnified with magnitude fluctuations of up to $\sim$$0.05\,\rm mag$ and that the weak-lensing shear can be up to the $\sim$$10^{-2}$ level in the outer regions of large voids. Smaller or shallower voids produce proportionally smaller effects. We investigate the "wiggling" of the primary cosmic microwave background (CMB) temperature anisotropies caused by intervening cosmic voids. The void-wiggling of primary CMB temperature gradients is of the opposite sign to that caused by galaxy clusters. Only extremely large and deep voids can produce wiggling amplitudes similar to galaxy clusters, $\sim$$\rm 15\,\mu K$ by a large void of radius $\sim$$4^\circ$ and central density contrast $-0.9$ at redshift 0.5 assuming a CMB background gradient of $\sim$$\rm10\,\mu K\, arcmin^{-1}$. The dipole signal is spread over the entire void area, and not concentrated at the lens' center as it is for clusters. Finally we use our model to simulate CMB sky maps lensed by large cosmic voids. Our embedded theory can easily be applied to more complicated void models and used to study gravitational lensing of the CMB, to probe dark-matter profiles, to reduce the lensing-induced systematics in supernova Hubble diagrams, as well as study the integrated Sachs-Wolfe effect.

We describe a method using the integrated Sachs-Wolfe (ISW) effect caused by individual inhomogeneities to determine the cosmological parameters, $H_0$, $\Omega_{\rm m}$, and $\Omega_\Lambda$, etc. This ISW-redshift test requires detailed knowledge of the internal kinematics of a set of individual density perturbations, e.g., galaxy clusters and/or cosmic voids, in particular their density and velocity profiles, and their mass accretion rates. It assumes the density perturbations are isolated and imbedded (equivalently compensated) and makes use of the newly found relation between the ISW temperature perturbation of the CMB and the Fermat potential of the lens. Given measurements of the amplitudes of the temperature variations in the CMB caused by such clusters or voids at various redshifts and estimates of their angular sizes or masses, one can constrain the cosmological parameters. More realistically, the converse is more likely, i.e., if the background cosmology is sufficiently constrained, measurement of ISW profiles of clusters and voids (e.g., hot and cold spots and rings) can constrain dynamical properties of the dark matter, including accretion, associated with such lenses and thus constrain the evolution of these objects with redshift.

We derive the Fermat potential for a spherically symmetric lens embedded in an FLRW cosmology and use it to investigate the late-time integrated Sachs-Wolfe (ISW) effect, i.e., secondary temperature fluctuations in the cosmic microwave background (CMB) caused by individual large scale clusters and voids. We present a simple analytical expression for the temperature fluctuation in the CMB across such a lens as a derivative of the lens' Fermat potential. This formalism is applicable to both linear and nonlinear density evolution scenarios, to arbitrarily large density contrasts, and to all open and closed background cosmologies. It is much simpler to use and makes the same predictions as conventional approaches. In this approach the total temperature fluctuation can be split into a time-delay part and an evolutionary part. Both parts must be included for cosmic structures that evolve and both can be equally important. We present very simple ISW models for cosmic voids and galaxy clusters to illustrate the ease of use of our formalism. We use the Fermat potentials of simple cosmic void models to compare predicted ISW effects with those recently extracted from WMAP and \emph{Planck} data by stacking large cosmic voids using the aperture photometry method. If voids in the local universe with large density contrasts are no longer evolving we find that the time delay contribution alone predicts values consistent with the measurements. However, we find that for voids still evolving linearly, the evolutionary contribution cancels a significant part of the time delay contribution and results in predicted signals that are much smaller than recently observed.

We develop a method to estimate photometric metallicities by simultaneously fitting the dereddened colors u-g, g-r, r-i and i-z from the SDSS with those predicted by the metallicity-dependent stellar loci. The method is tested with a spectroscopic sample of main-sequence stars in Stripe 82 selected from the SDSS DR9 and three open clusters. With 1 per cent photometry, the method is capable of delivering photometric metallicities precise to about 0.05, 0.12, and 0.18 dex at metallicities of 0.0, -1.0, and -2.0, respectively, comparable to the precision achievable with low-resolution spectroscopy at a signal-to-noise ratio of 10. We apply this method to the re-calibrated Stripe 82 catalog and derive metallicities for about 0.5 million stars of colors 0.3 < g-i < 1.6 mag and distances between 0.3 -- 18 kpc. Potential systematics in the metallicities thus derived, due to the contamination of giants and binaries, are investigated. Photometric distances are also calculated. About 91, 72, and 53 per cent of the sample stars are brighter than r = 20.5, 19.5, and 18.5 mag, respectively. The median metallicity errors are around 0.19, 0.16, 0.11, and 0.085 dex for the whole sample, and for stars brighter than r = 20.5, 19.5, and 18.5 mag, respectively. The median distance errors are 8.8, 8.4, 7.7, and 7.3 per cent for the aforementioned four groups of stars, respectively. The data are publicly available. Potential applications of the data in studies of the distribution, (sub)structure, and chemistry of the Galactic stellar populations, are briefly discussed. The results will be presented in future papers.

We re-estimate the peculiar velocity of the Sun with respect to the local standard of rest using a sample of local stars within 600 pc of the Sun, selected from the LAMOST Spectroscopic Survey of the Galactic Anti-centre (LSS-GAC). The sample consists of 94332 FGK main-sequence stars with well-determined radial velocities and atmospheric parameters. To derive the LSR, two independent analyses are applied to the data. Firstly, we determine the solar motion by comparing the observed velocity distribution to that generated with the analytic formulism of Schonrich & Binney that has been demonstrated to show excellent agreement with rigorous torus-based dynamics modelling by Binney & McMillan. Secondly, we propose that cold populations of thin disc stars, selected by applying an orbital eccentricity cut, can be directly used to determine the LSR without the need of asymmetric drift corrections. Both approaches yield consistent results of solar motion in the direction of Galactic rotation, V_sun, that are much higher than the standard value adopted hitherto, derived from Stromgren's equation. The newly deduced values of V_sun are 1-2 km/s smaller than the more recent estimates derived from the Geneva-Copenhagen Survey sample of stars in the solar neighbourhood (within 100 pc). We attribute the small difference to the presence of several well-known moving groups in the GCS sample that, fortunately, hardly affect the LSS-GAC sample. The newly derived radial and vertical components of the solar motion agree well with the previous studies. In addition, for all components of the solar motion, the values yielded by stars of different spectral types in the LSS-GAC sample are consistent with each other, suggesting that the local disk is well relaxed and that the LSR reported in the current work is robust. Our final recommended LSR is, (U,V,W)_sun = (7.01+/-0.20, 10.13+/-0.12, 4.95+/-0.09) km/s.

We investigate the correlation between extinction and H~{\sc i} and CO emission at intermediate and high Galactic latitudes ($|b|>10\degr$) within the footprint of the Xuyi Schmidt Telescope Photometric Survey of the Galactic anticentre (XSTPS-GAC) on small and large scales. In Paper I (Chen et al. 2014), we present a three-dimensional dust extinction map within the footprint of XSTPS-GAC, covering a sky area of over 6,000\,deg$^2$ at a spatial angular resolution of 6\,arcmin. In the current work, the map is combined with data from gas tracers, including H~{\sc i} data from the Galactic Arecibo L-band Feed Array H~{\sc i} survey and CO data from the Planck mission, to constrain the values of dust-to-gas ratio $DGR=A_V/N({\rm H})$ and CO-to-$\rm H_2$ conversion factor $X_{\rm CO}=N({\rm H_2})/W_{\rm CO}$ for the entire GAC footprint excluding the Galactic plane, as well as for selected star-forming regions (such as the Orion, Taurus and Perseus clouds) and a region of diffuse gas in the northern Galactic hemisphere. For the whole GAC footprint, we find $DGR=(4.15\pm0.01) \times 10^{-22}$\,$\rm mag\,cm^{2}$ and $X_{\rm CO}=(1.72 \pm 0.03) \times 10^{20}$\,$\rm cm^{-2}\,(K\,km\,s^{-1})^{-1}$. We have also investigated the distribution of "CO-dark" gas (DG) within the footprint of GAC and found a linear correlation between the DG column density and the $V$-band extinction: $N({\rm DG}) \simeq 2.2 \times 10^{21} (A_V - A^{c}_{V})\,\rm cm^{-2}$. The mass fraction of DG is found to be $f_{\rm DG}\sim 0.55$ toward the Galactic anticentre, which is respectively about 23 and 124 per cent of the atomic and CO-traced molecular gas in the same region. This result is consistent with the theoretical work of Papadopoulos et al. but much larger than that expected in the $\rm H_2$ cloud models by Wolfire et al.

We apply the Statefinder hierarchy and the growth rate of matter perturbations to discriminate modified Chaplygin gas (MCG), generalized Chaplygin gas (GCG), superfluid Chaplygin gas (SCG), purely kinetic k-essence (PKK), and $\Lambda$CDM model. We plot the evolutional trajectories of these models in the statefinder plane and in the composite diagnostic plane. We find that GCG, MCG, SCG, PKK, and $\Lambda$CDM can be distinguished well from each other at the present epoch by using the composite diagnostic $\{\epsilon(z), S^{(1)}_{5}\}$. Using other combinations, such as $\{S^{(1)}_{3}, S^{(1)}_4\}$, $\{S^{(1)}_{3}, S_{5}\}$, $\{\epsilon(z), S^{(1)}_{3}\}$, and $\{\epsilon(z), S_4 \}$, some of these five dark energy models cannot be distinguished.

We have developed and implemented an iterative algorithm of flux calibration for the LAMOST Spectroscopic Survey of the Galactic anti-center (LSS-GAC). For a given LSS-GAC plate, the spectra are first processed with a set of nominal spectral response curves (SRCs) and used to derive initial stellar atmospheric parameters (effective temperature $T_{\rm eff}$, surface gravity log\,$g$ and metallicity [Fe/H]) as well as dust reddening $E(B-V)$ of all targeted stars. For each of the sixteen spectrographs, several F-type stars of good signal-to-noise ratios (SNRs) are then selected as flux standard stars for further, iterative spectral flux calibration. Comparison of spectrophotometric colours, deduced from the flux-calibrated spectra, with the photometric measurements yields average differences of 0.02$\pm$0.07 and $-$0.04$\pm$0.09\,mag for the $(g-r)$ and $(g-i)$, respectively. The relatively large negative offset in $(g-i)$ is due to the fact that we have opted not to correct for the telluric bands, most notably the atmospheric A-band in the wavelength range of $i$-band. Comparison of LSS-GAC multi-epoch observations of duplicate targets indicates that the algorithm has achieved an accuracy of about 10 per cent in relative flux calibration for the wavelength range 4000 -- 9000\,\AA. The shapes of SRC deduced for the individual LAMOST spectrographs are found to vary by up to 30 per cent for a given night, and larger for different nights, indicating that the derivation of SRCs for the individual plates is essential in order to achieve accurate flux calibration for the LAMOST spectra.

We introduce the LAMOST Stellar Parameter Pipeline at Peking University --- LSP3, developed and implemented for the determinations of radial velocity $V_{\rm r}$ and stellar atmospheric parameters (effective temperature $T_{\rm eff}$, surface gravity log\,$g$, metallicity [Fe/H]) for the LAMOST Spectroscopic Survey of the Galactic Anti-center (LSS-GAC). We describe the algorithms of LSP3 and examine the accuracy of parameters yielded by it. The precision and accuracy of parameters yielded are investigated by comparing results of multi-epoch observations and of candidate members of open and globular clusters, with photometric calibration, as well as with independent determinations available from a number of external databases, including the PASTEL archive, the APOGEE, SDSS and RAVE surveys, as well as those released in the LAMOST DR1. The uncertainties of LSP3 parameters are characterized and quantified as a function of the spectral signal-to-noise ratio (SNR) and stellar atmospheric parameters. We conclude that the current implementation of LSP3 has achieved an accuracy of 5.0\,km\,s$^{-1}$, 150\,K, 0.25\,dex, 0.15\,dex for the radial velocity, effective temperature, surface gravity and metallicity, respectively, for LSS-GAC spectra of FGK stars of SNRs per pixel higher than 10. The LSP3 has been applied to over a million LSS-GAC spectra collected hitherto. Stellar parameters yielded by the LSP3 will be released to the general public following the data policy of LAMOST, together with estimates of the interstellar extinction $E(B-V)$ and stellar distances, deduced by combining spectroscopic and multi-band photometric measurements using a variety of techniques.

As a major component of the LAMOST Galactic surveys, the LAMOST Spectroscopic Survey of the Galactic Anti-centre (LSS-GAC) aims to survey a significant volume of the Galactic thin/thick discs and halo for a contiguous sky area of over 3,400 deg$^2$ centred on the Galactic anti-centre ($|b| \leq 30^{\circ}$, $150 \leq l \leq 210^{\circ}$), and obtain $\lambda\lambda$3700 -- 9000 low resolution ($R \sim 1,800$) spectra for a statistically complete sample of $\sim 3$\,M stars of all colours down to a limiting magnitued of $r$ $\sim$ 17.8\,mag (to 18.5\,mag for limited fields). Together with Gaia, the LSS-GAC will yield a unique dataset to advance our understanding of the structure and assemblage history of the Galaxy, in particular its disk(s). In addition to the main survey, the LSS-GAC will also target hundreds of thousands objects in the vicinity fields of M\,31 and M\,33 and survey a significant fraction (over a million) of randomly selected very bright stars (VB; $r \le 14$ mag) in the northern hemisphere. During the Pilot and the first year Regular Surveys of LAMOST, a total of 1,042,586 [750,867] spectra of a signal to noise ratio S/N(7450\AA) $\ge$ 10 [S/N(4650\AA) $\ge$ 10] have been collected. In this paper, we present a detailed description of the target selection algorithm, survey design, observations and the first data release of value-added catalogues (including radial velocities, effective temperatures, surface gravities, metallicities, values of interstellar extinction, distances, proper motions and orbital parameters) of the LSS-GAC.

We propose a Stellar Locus OuTlier (SLOT) method to determine the binary fraction of main-sequence stars statistically. The method is sensitive to neither the period nor mass-ratio distributions of binaries, and able to provide model-free estimates of binary fraction for large numbers of stars of different populations in large survey volumes. We have applied the SLOT method to two samples of stars from the SDSS Stripe 82, constructed by combining the re-calibrated SDSS photometric data with respectively the spectroscopic information from the SDSS and LAMOST surveys. For the SDSS spectroscopic sample, we find an average binary fraction for field FGK stars of $41%\pm2%$. The fractions decrease toward late spectral types, and are respectively $44%\pm5%$, $43%\pm3%$, $35%\pm5%$, and $28%\pm6%$ for stars of $g-i$ colors between 0.3 -- 0.6, 0.6 -- 0.9, 0.9 -- 1.2, and 1.2 - 1.6\,mag. A modest metallicity dependence is also found. The fraction decreases with increasing metallicity. For stars of [Fe/H] between $-0.5$ - 0.0, $-1.0$ - $-0.5$, $-1.5$ - $-1.0$, and $-2.0$ - $-1.5$\,dex, the inferred binary fractions are $37%\pm3%$, $39%\pm3%$, $50%\pm9%$, and $53%\pm20%$, respectively. We have further divided the sample into stars from the thin disk, the thick disk, the transition zone between them, and the halo. The results suggest that the Galactic thin and thick disks have comparable binary fractions, whereas the Galactic halo contains a significantly larger fraction of binaries. Applying the method to the LAMOST spectroscopic sample yields consistent results. Finally, other potential applications and future work with the method are discussed.

Stellar loci are widely used for selection of interesting outliers, reddening determinations, and calibrations. However, hitherto the dependence of stellar loci on metallicity has not been fully explored and their intrinsic widths are unclear. In this paper, by combining the spectroscopic and re-calibrated imaging data of the SDSS Stripe 82, we have built a large, clean sample of dwarf stars with accurate colors and well determined metallicities to investigate the metallicity dependence and intrinsic widths of the SDSS stellar loci. Typically, one dex decrease in metallicity causes 0.20 and 0.02 mag decrease in colors u-g and g-r, and 0.02 and 0.02 mag increase in colors r-i and i-z, respectively. The variations are larger for metal-rich stars than for metal-poor ones, and for F/G/K stars than for A/M ones. Using the sample, we have performed two dimensional polynomial fitting to the u-g, g-r, r-i, and i-z colors as a function of color g-i and metallicity [Fe/H]. The residuals, at the level of 0.029, 0.008, 0.008 and 0.011 mag for the u-g, g-r, r-i, and i-z colors, respectively can be fully accounted for by the photometric errors and metallicity uncertainties, suggesting that the intrinsic widths of the loci are at maximum a few mmag. The residual distributions are asymmetric, revealing that a significant fraction of stars are binaries. In a companion paper, we will present an unbiased estimate of the binary fraction for field stars. Other potential applications of the metallicity dependent stellar loci are briefly discussed.

In this paper, we propose a spectroscopy based Stellar Color Regression (SCR) method to perform accurate color calibration for modern imaging surveys, taking advantage of millions of stellar spectra now available. The method is straightforward, insensitive to systematic errors in the spectroscopically determined stellar atmospheric parameters, applicable to regions that are effectively covered by spectroscopic surveys, and capable of delivering an accuracy of a few millimagnitudes for color calibration. As an illustration, we have applied the method to the SDSS Stripe 82 data (Ivezic et al; I07 hereafter). With a total number of 23,759 spectroscopically targeted stars, we have mapped out the small but strongly correlated color zero point errors present in the photometric catalog of Stripe 82, and improve the color calibration by a factor of 2 -- 3. Our study also reveals some small but significant magnitude dependence errors in z-band for some CCDs. Such errors are likely to be present in all the SDSS photometric data. Our results are compared with those from a completely independent test based on the intrinsic colors of red galaxies presented by I07. The comparison as well as other tests shows that the SCR method has achieved a color calibration internally consistent at a level of about 5 mmag in u-g, 3 mmag in g-r, and 2 mmag in r-i and i-z, respectively. Given the power of the SCR method, we discuss briefly the potential benefits by applying the method to existing, on-going, and up-coming imaging surveys.

Magnetic reconnection is believed to be the dominant energy release mechanism in solar flares. The standard flare model predicts both downward and upward outflow plasmas with speeds close to the coronal Alfv\'{e}n speed. Yet, spectroscopic observations of such outflows, especially the downflows, are extremely rare. With observations of the newly launched Interface Region Imaging Spectrograph (IRIS), we report the detection of greatly redshifted ($\sim$125 km s$^{-1}$ along line of sight) Fe {\sc{xxi}} 1354.08\AA{} emission line with a $\sim$100 km s$^{-1}$ nonthermal width at the reconnection site of a flare. The redshifted Fe {\sc{xxi}} feature coincides spatially with the loop-top X-Ray source observed by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). We interpret this large redshift as the signature of downward-moving reconnection outflow/hot retracting loops. Imaging observations from both IRIS and the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) also reveal the eruption and reconnection processes. Fast downward-propagating blobs along these loops are also found from cool emission lines (e.g., Si {\sc{iv}}, O {\sc{iv}}, C {\sc{ii}}, Mg {\sc{ii}}) and images of AIA and IRIS. Furthermore, the entire Fe {\sc{xxi}} line is blueshifted by $\sim$260 km s$^{-1}$ at the loop footpoints, where the cool lines mentioned above all exhibit obvious redshift, a result that is consistent with the scenario of chromospheric evaporation induced by downward-propagating nonthermal electrons from the reconnection site.

Magnetic flux ropes (MFRs) are believed to be at the heart of solar coronal mass ejections (CMEs). A well-known example is the prominence cavity in the low corona that sometimes makes up a three-part white-light (WL) CME upon its eruption. Such a system, which is usually observed in quiet-Sun regions, has long been suggested to be the manifestation of an MFR with relatively cool filament material collecting near its bottom. However, observational evidence of eruptive, filament-hosting MFR systems has been elusive for those originating in active regions. By utilizing multi-passband extreme-ultra-violet (EUV) observations from SDO/AIA, we present direct evidence of an eruptive MFR in the low corona that exhibits a hot envelope and a cooler core; the latter is likely the upper part of a filament that undergoes a partial eruption, which is later observed in the upper corona as the coiled kernel of a fast, WL CME. This MFR-like structure exists more than one hour prior to its eruption, and displays successive stages of dynamical evolution, in which both ideal and non-ideal physical processes may be involved. The timing of the MFR kinematics is found to be well correlated with the energy release of the associated long-duration C1.9 flare. We suggest that the long-duration flare is the result of prolonged energy release associated with the vertical current sheet induced by the erupting MFR.

We present a three dimensional extinction map in $r$ band. The map has a spatial angular resolution, depending on latitude, between 3 -- 9\,arcmin and covers the entire XSTPS-GAC survey area of over 6,000\,$\rm deg^2$ for Galactic longitude $\rm 140 \leq$ $l$ $ \leq 220\deg$ and latitude $\rm -40\leq$ $b$ $ \leq 40\deg$. By cross-matching the photometric catalog of the Xuyi Schmidt Telescope Photometric Survey of the Galactic Anticentre (XSTPS-GAC) with those of 2MASS and WISE, we have built a multi-band photometric stellar sample of about 30 million stars and applied spectral energy distribution (SED) fitting to the sample. By combining photometric data from the optical to the near-infrared, we are able to break the degeneracy between the intrinsic stellar colours and the amounts of extinction by dust grains for stars with high photometric accuracy, and trace the extinction as a function of distance for low Galactic latitude and thus highly extincted regions. This has allowed us to derive the best-fit extinction and distance information of more than 13 million stars, which are used to construct the three dimensional extinction map. We have also applied a Rayleigh-Jeans colour excess (RJCE) method to the data using the 2MASS and WISE colour $(H-W2)$. The resulting RJCE extinction map is consistent with the integrated two dimensional map deduced using the best-fit SED algorithm. However for individual stars, the amounts of extinction yielded by the RJCE method suffer from larger errors than those given by the best-fit SED algorithm.

Microlensing has proven an effective probe of the structure of the innermost regions of quasars, and an important test of accretion disk models. We present light curves of the lensed quasar HE 0435-1223 in the R band and in the ultraviolet, and consider them together with X-ray light curves in two energy bands that are presented in a companion paper. Using a Bayesian Monte Carlo method, we constrain the size of the accretion disk in the rest-frame near- and far-UV, and constrain for the first time the size of the X-ray emission regions in two X-ray energy bands. The R-band scale size of the accretion disk is about 10^15.23 cm (~23 r_g), slightly smaller than previous estimates, but larger than would be predicted from the quasar flux. In the UV, the source size is weakly constrained, with a strong prior dependence. The UV to R-band size ratio is consistent with the thin disk model prediction, with large error bars. In soft and hard X-rays, the source size is smaller than ~10^14.8 cm (~10 r_g) at 95% confidence. We do not find evidence of structure in the X-ray emission region, as the most likely value for the ratio of the hard X-ray size to the soft X-ray size is unity. Finally, we find that the most likely value for the mean mass of stars in the lens galaxy is ~0.3 M_sun, consistent with other studies.

Three dimensional interstellar extinction maps provide a powerful tool for stellar population analysis. We use data from the VISTA Variables in the Via Lactea survey together with the Besan\c{c}on stellar population synthesis model of the Galaxy to determine interstellar extinction as a function of distance in the Galactic bulge covering $ -10 < l < 10$ and $-10 < b <5$. We adopted a recently developed method to calculate the colour excess. First we constructed the H-Ks vs. Ks and J-Ks vs. Ks colour-magnitude diagrams based on the VVV catalogues that matched 2MASS. Then, based on the temperature-colour relation for M giants and the distance-colour relations, we derived the extinction as a function of distance. The observed colours were shifted to match the intrinsic colours in the Besan\c{c}on model as a function of distance iteratively. This created an extinction map with three dimensions: two spatial and one distance dimension along each line of sight towards the bulge. We present a 3D extinction map that covers the whole VVV area with a resolution of 6' x 6', using distance bins of 0.5 kpc. The high resolution and depth of the photometry allows us to derive extinction maps for a range of distances up to 10 kpc and up to 30 magnitudes of extinction in $A_{V}$. Integrated maps show the same dust features and consistent values as other 2D maps. We discuss the spatial distribution of dust features in the line of sight, which suggests that there is much material in front of the Galactic bar, specifically between 5-7 kpc. We compare our dust extinction map with high-resolution $\rm ^{12}CO$ maps towards the Galactic bulge, where we find a good correlation between $\rm ^{12}CO$ and $\rm A_{V}$. We determine the X factor by combining the CO map and our dust extinction map. Our derived average value is consistent with the canonical value of the Milky Way.

Since it can describe the rotation curves of galaxies without dark matter and can give rise to accelerated expansion, conformal gravity attracts much attention recently. As a theory of modified gravity, it is important to test conformal gravity with astrophysical observations. Here we constrain conformal gravity with SNIa and Hubble parameter data and investigate whether it suffers from an age problem with the age of APM~08279+5255. We find conformal gravity can accommodate the age of APM~08279+5255 at 3 $\sigma$ deviation, unlike most of dark energy models which suffer from an age problem.

The GLT project is deploying a new submillimeter (submm) VLBI station in Greenland. Our primary scientific goal is to image a shadow of the supermassive black hole (SMBH) of six billion solar masses in M87 at the center of the Virgo cluster of galaxies. The expected SMBH shadow size of 40-50 $\mu$as requires superbly high angular resolution, suggesting that the submm VLBI would be the only way to obtain the shadow image. The Summit station in Greenland enables us to establish baselines longer than 9,000 km with ALMA in Chile and SMA in Hawaii as well as providing a unique $u$--$v$ coverage for imaging M87. Our VLBI network will achieve a superior angular resolution of about 20 $\mu$as at 350 GHz, corresponding to $\sim2.5$ times of the Schwarzschild radius of the supermassive black hole in M87. We have been monitoring the atmospheric opacity at 230 GHz since August. 2011; we have confirmed the value on site during the winter season is comparable to the ALMA site thanks to high altitude of 3,200 m and low temperature of $-50\degr$C. We will report current status and future plan of the GLT project towards our expected first light on 2015--2016.

We present a simplified version of the lowest-order embedded point mass gravitational lens theory and then make the extension of this theory to any embedded transparent lens. Embedding a lens effectively reduces the gravitational potential's range, i.e., partially shields the lensing potential because the lens mass is made a contributor to the mean mass density of the universe and not simply superimposed upon it. We give the time-delay function for the embedded point mass lens from which we can derive the simplified lens equation by applying Fermat's least-time principle. Even though rigorous derivations are only made for the point mass in a flat background, the generalization of the lens equation to lowest-order for any distributed lens in any homogeneous background is obvious. We find from this simplified theory that embedding can introduce corrections above the few percent level in weak lensing shears caused by large clusters but only at large impacts. The potential part of the time delay is also affected in strong lensing at the few percent level. Additionally we again confirm that the presence of a cosmological constant alters the gravitational deflection of passing photons.

We homogeneously analyzed the \chandra\ X-ray observations of 10 gravitational lenses, HE 0047-1756, QJ 0158-4325, SDSS 0246-0805, HE 0435-1223, SDSS 0924+0219, SDSS 1004+4112, HE 1104-1805, PG 1115+080, Q 1355-2257, and Q 2237+0305, to measure the differential X-ray absorption between images, the metallicity, and the dust-to-gas ratio of the lens galaxies. We detected differential absorption in all lenses except SDSS 0924+0219 and HE 1104-1805. This doubles the sample of dust-to-gas ratio measurements in cosmologically distant lens galaxies. We successfully measured the gas phase metallicity of three lenses, Q 2237+0305, SDSS 1004+4112, and B 1152+199 from the X-ray spectra. Our results suggest a linear correlation between metallicity and dust-to-gas ratio (i.e., a constant metal-to-dust ratio), consistent with what is found for nearby galaxies. We obtain an average dust-to-gas ratio $E(B-V)/N_H=1.17^{+0.41}_{-0.31} \times 10^{-22}\rm mag\,cm^2\,atom^{-1}$ in the lens galaxies, with an intrinsic scatter of $\rm0.3\,dex$. Combining these results with data from GRB afterglows and quasar foreground absorbers, we found a mean dust-to-gas ratio $\mdtg,$ now significantly lower than the average Galactic value, $1.7\,\times 10^{-22}\,\rm mag\, cm^{2}\, atoms^{-1}.$ This suggests evolution of dust-to-gas ratios with redshift and lower average metallicities for the higher redshift galaxies, consistent with current metal and dust evolution models of interstellar medium. The slow evolution in the metal-to-dust ratio with redshift implies very rapid dust formation in high redshift ($z>2$) galaxies.

The gravitationally lensed quasar HE 1104-1805 has been observed at a variety of wavelengths ranging from the mid-infrared to X-ray for nearly 20 years. We combine flux ratios from the literature, including recent Chandra data, with new observations from the SMARTS telescope and HST, and use them to investigate the spatial structure of the central regions using a Bayesian Monte Carlo analysis of the microlensing variability. The wide wavelength coverage allows us to constrain not only the accretion disk half-light radius r_1/2, but the power-law slope \xi\ of the size-wavelength relation r_1/2 ~ \lambda^\xi. With a logarithmic prior on the source size, the (observed-frame) R-band half-light radius log(r_1/2/cm) is 16.0+0.3-0.4, and the slope \xi\ is 1.0+0.30-0.56. We put upper limits on the source size in soft (0.4-1.2 keV) and hard (1.2-8 keV) X-ray bands, finding 95% upper limits on log (r_1/2/cm) of 15.33 in both bands. A linear prior yields somewhat larger sizes, particularly in the X-ray bands. For comparison, the gravitational radius, using a black hole mass estimated using the H\beta\ line, is log(r_g/cm) = 13.94. We find that the accretion disk is probably close to face-on, with cos i = 1.0 being four times more likely than cos i = 0.5. We also find probability distributions for the mean mass of the stars in the foreground lensing galaxy, the direction of the transverse peculiar velocity of the lens, and the position angle of the projected accretion disk's major axis (if not face-on).

We use gravitational microlensing to determine the size of the X-ray and optical emission regions of the quadruple lens system Q 2237+0305. The optical half-light radius, log(R_{1/2,V}/cm)=16.41\pm0.18 (at lambda_{rest}=2018 \AA), is significantly larger than the observed soft, log(R_{1/2,soft}/cm)=15.76^{+0.41}_{-0.34} (1.1-3.5 keV in the rest frame), and hard, log(R_{1/2,hard}/cm)=15.46^{+0.34}_{-0.29} (3.5-21.5 keV in the rest frame), band X-ray emission. There is a weak evidence that the hard component is more compact than the soft, with log(R_{1/2,soft}/R_{1/2,hard}) \sim 0.30^{+0.53}_{-0.45}. This wavelength-dependent structure agrees with recent results found in other lens systems using microlensing techniques, and favors geometries in which the corona is concentrated near the inner edge of the accretion disk. While the available measurements are limited, the size of the X-ray emission region appears to be roughly proportional to the mass of the central black hole.

We report observations of type III radio bursts at decimeter wavelengths (type IIIdm bursts) -- signatures of suprathermal electron beams propagating in the low corona -- using the new technique of radio dynamic imaging spectroscopy provided by the recently upgraded Karl G. Jansky Very Large Array (VLA). For the first time, type IIIdm bursts were imaged with high time and frequency resolution over a broad frequency band, allowing electron beam trajectories in the corona to be deduced. Together with simultaneous hard X-ray (HXR) and extreme ultraviolet (EUV) observations, we show these beams emanate from an energy release site located in the low corona at a height below ~15 Mm, and propagate along a bundle of discrete magnetic loops upward into the corona. Our observations enable direct measurements of the plasma density along the magnetic loops, and allow us to constrain the diameter of these loops to be less than 100 km. These over-dense and ultra-thin loops reveal the fundamentally fibrous structure of the Sun's corona. The impulsive nature of the electron beams, their accessibility to different magnetic field lines, and the detailed structure of the magnetic release site revealed by the radio observations indicate that the localized energy release is highly fragmentary in time and space, supporting a bursty reconnection model that involves secondary magnetic structures for magnetic energy release and particle acceleration.

Recent quasar microlensing observations have constrained the sizes of X-ray emission regions to be within about 10 gravitational radii of the central supermassive black hole. Therefore, the X-ray emission from lensed quasars is first strongly lensed by the black hole before it is lensed by the foreground galaxy and star fields. We present a scheme that combines the initial strong lensing of a Kerr black hole with standard linearized microlensing by intervening stars. We find that X-ray microlensed light curves incorporating Kerr strong gravity can differ significantly from standard curves. The amplitude of the fluctuations in the light curves can increase or decrease by ~0.65-0.75 mag by including Kerr strong gravity. Larger inclination angles give larger amplitude fluctuations in the microlensing light curves. Consequently, current X-ray microlensing observations might have under or overestimated the sizes of the X-ray emission regions. We estimate this bias using a simple metric based on the amplitude of magnitude fluctuations. The half light radius of the X-ray emission region can be underestimated up to ~50% or overestimated up to ~20%. Underestimates are found in most situations we have investigated. The only exception is for a disk with large spin, radially flat emission profile, and observed nearly face on, where an overestimate is found. Thus, more accurate microlensing size constraints should be obtainable by including Kerr lensing. The caustic crossing time can differ by months after including Kerr strong gravity. A simultaneous monitoring of gravitational lensed quasars in both X-ray and optical bands with densely sampled X-ray light curves might reveal this feature. We conclude that it should be possible to constrain important parameters such as inclination angles and black hole spins from combined Kerr and microlensing effects.

Studies of the properties of the inner Galactic Bulge depend strongly on the assumptions about the interstellar extinction. Most of the extinction maps available in the literature lack the information about the distance. We combine the observations with the Besancon model of the Galaxy to investigate the variations of extinction along different lines of sight towards the inner Galactic bulge as a function of distance. In addition we study the variations in the extinction law in the Bulge. We construct color-magnitude diagrams with the following sets of colors: H-Ks and J-Ks from the VVV catalogue as well as Ks-[3.6], Ks-[4.5], Ks-[5.8] and Ks-[8.0] from GLIMPSE-II catalogue matched with 2MASS. Using the newly derived temperature-color relation for M giants that match better the observed color-magnitude diagrams we then use the distance-color relations to derive the extinction as a function of distance. The observed colors are shifted to match the intrinsic colors in the Besan\c{c}on model, as a function of distance, iteratively thereby creating an extinction map with three dimensions: two spatial and one distance dimension along each line of sight towards the bulge. Colour excess maps are presented at a resolution of 15' x 15' for 6 different combinations of colors, in distance bins of 1 kpc. The high resolution and depth of the photometry allows to derive extinction maps to 10 kpc distance and up to 35 magnitudes of extinction in Av (3.5 mag in Aks). Integrated maps show the same dust features and consistent values with the other 2D maps. Starting from the color excess between the observations and the model we investigate the extinction law in near-infrared and its variation along different lines of sight.

Recent quasar microlensing observations have constrained the X-ray emission sizes of quasars to be about 10 gravitational radii, one order of magnitude smaller than the optical emission sizes. Using a new ray-tracing code for the Kerr space-time, we find that the observed X-ray flux is strongly influenced by the gravity field of the central black hole, even for observers at moderate inclination angles. We calculate inclination-dependent flux profiles of active galactic nuclei in the optical and X-ray bands by combining the Kerr lensing and projection effects for future references. We further study the dependence of the X-ray-to-optical flux ratio on the inclination angle caused by differential lensing distortion of the X-ray and optical emission, assuming several corona geometries. The strong lensing X-ray-to-optical magnification ratio can change by a factor of ~10 for normal quasars in some cases, and another factor of ~10 for broad absorption line quasars (BALs) and obscured quasars. Comparing our results with the observed distributions in normal and broad absorption line quasars, we find that the inclination angle dependence of the magnification ratios can change the X-ray-to-optical flux ratio distributions significantly. In particular, the mean value of the spectrum slope parameter $\alpha_{ox},$ $0.3838\log F_{2 keV}/F_{2500 {\AA}}$, can differ by ~0.1-0.2 between normal and broad absorption line quasars, depending on corona geometries, suggesting larger intrinsic absorptions in BALs.

3-D astrophysical atmospheres will have random velocity fields. We seek to combine the methods we have developed for solving the 1-D problem with arbitrary flows to those that we have developed for solving the fully 3-D relativistic radiative transfer problem in the case of monotonic flows. The methods developed in the case of 3-D atmospheres with monotonic flows, solving the fully relativistic problem along curves defined by an affine parameter, are very flexible and can be extended to the case of arbitrary velocity fields in 3-D. Simultaneously, the techniques we developed for treating the 1-D problem with arbitrary velocity fields are easily adapted to the 3-D problem. The algorithm we present allows the solution of 3-D radiative transfer problems that include arbitrary wavelength couplings. We use a quasi-analytic formal solution of the radiative transfer equation that significantly improves the overall computation speed. We show that the approximate lambda operator developed in previous work gives good convergence, even neglecting wavelength coupling. Ng acceleration also gives good results. We present tests that are of similar resolution to what has been presented using Monte-Carlo techniques, thus our methods will be applicable to problems outside of our test setup. Additional domain decomposition parallelization strategies will be explored in future work.

In supernova spectroscopy relatively little attention has been given to the properties of optically thick spectral lines in epochs following the photosphere's recession. Most treatments and analyses of post-photospheric optical spectra of supernovae assume that forbidden-line emission comprises most if not all spectral features. However, evidence exists which suggests that some spectra exhibit line profiles formed via optically thick resonance-scattering even months or years after the supernova explosion. To explore this possibility we present a geometrical approach to supernova spectrum formation based on the "Elementary Supernova" model, wherein we investigate the characteristics of resonance-scattering in optically thick lines while replacing the photosphere with a transparent central core emitting non-blackbody continuum radiation, akin to the optical continuum provided by decaying 56Co formed during the explosion. We develop the mathematical framework necessary for solving the radiative transfer equation under these conditions, and calculate spectra for both isolated and blended lines. Our comparisons with analogous results from the Elementary Supernova code SYNOW reveal several marked differences in line formation. Most notably, resonance lines in these conditions form P Cygni-like profiles, but the emission peaks and absorption troughs shift redward and blueward, respectively, from the line's rest wavelength by a significant amount, despite the spherically symmetric distribution of the line optical depth in the ejecta. These properties and others that we find in this work could lead to misidentification of lines or misattribution of properties of line-forming material at post-photospheric times in supernova optical spectra.

A systematic study of RR Lyrae stars is performed based on a selected sample of 655 objects in the Large Magellanic Cloud with observation of long span and numerous measurements by the Optical Gravitational Lensing Experiment III project. The Phase Dispersion Method and linear superposition of the harmonic oscillations are used to derive the pulsation frequency and variation properties. It is found that there exists an Oo I and Oo II dichotomy in the LMC RR Lyrae stars. Due to our strict criteria to identify a frequency, a lower limit of the incidence rate of Blazhko modulation in LMC is estimated in various subclasses of RR Lyrae stars. For fundamental-mode RR Lyrae stars, the rate 7.5% is smaller than previous result. In the case of the first-overtone RR Lyr variables, the rate 9.1% is relatively high. In addition to the Blazhko variables, fifteen objects are identified to pulsate in the fundamental/first-overtone double mode. Furthermore, four objects show a period ratio around 0.6 which makes them very likely the rare pulsators in the fundamental/second-overtone double-mode.

The presence of two stellar populations in the Milky Way bulge has been reported recently. We aim at studying the abundances and kinematics of stars in the outer bulge, thereby providing additional constraints on models of its formation. Spectra of 401 red giant stars in a field at (l,b)=(0{\deg},-10{\deg}) were obtained with FLAMES at the VLT. Stars of luminosities down to below the two bulge red clumps (RCs) are included. From these spectra we measure general metallicities, abundances of Fe and the alpha-elements, and radial velocities (RV) of the stars. These measurements as well as photometric data are compared to simulations with the Besancon and TRILEGAL models of the Galaxy. We confirm the presence of two populations among our sample stars: i) a metal-rich one at [M/H] ~+0.3, comprising about 30% of the sample, with low RV dispersion and low alpha-abundance, and ii) a metal-poor population at [M/H] ~-0.6 with high RV dispersion and high alpha-abundance. The metal-rich population could be connected to the Galactic bar. We identify this population as the carrier of the double RC feature. We do not find a significant difference in metallicity or RV between the two RCs, a small difference in metallicity being probably due to a selection effect. The RV dispersion agrees well with predictions of the Besancon Galaxy model, but the metallicity of the "thick bulge" model component should be shifted to lower metallicity by 0.2 to 0.3dex to well reproduce the observations. We present evidence that the metallicity distribution function depends on the evolutionary state of the sample stars, suggesting that enhanced mass loss preferentially removes metal-rich stars. We also confirm the decrease of \alpha-element over-abundance with increasing metallicity.

We give analytic expressions for image properties of objects seen around point mass lenses embedded in a flat $\Lambda$CDM universe. An embedded lens in an otherwise homogeneous universe offers a more realistic representation of the lens's gravity field and its associated deflection properties than does the conventional linear superposition theory. Embedding reduces the range of the gravitational force acting on passing light beams thus altering all quantities such as deflection angles, amplifications, shears and Einstein ring sizes. Embedding also exhibits the explicit effect of the cosmological constant on these same lensing quantities. In this paper we present these new results and demonstrate how they can be used. The effects of embedding on image properties, although small i.e., usually less than a fraction of a percent, have a more pronounced effect on image distortions in weak lensing where the effects can be larger than 10%. Embedding also introduces a negative surface mass density for both weak and strong lensing, a quantity altogether absent in conventional Schwarzschild lensing. In strong lensing we find only one additional quantity, the potential part of the time delay, which differs from conventional lensing by as much as 4%, in agreement with our previous numerical estimates.

We present four new seasons of optical monitoring data and six epochs of X-ray photometry for the doubly-imaged lensed quasar Q J0158-4325. The high-amplitude, short-period microlensing variability for which this system is known has historically precluded a time delay measurement by conventional methods. We attempt to circumvent this limitation by application of a Monte Carlo microlensing analysis technique, but we are only able to prove that the delay must have the expected sign (image A leads image B). Despite our failure to robustly measure the time delay, we successfully model the microlensing at optical and X-ray wavelengths to find a half light radius for soft X-ray emission log(r_{1/2,X,soft}/cm) = 14.3^{+0.4}_{-0.5}, an upper limit on the half-light radius for hard X-ray emission log(r_{1/2,X,hard}/cm) <= 14.6 and a refined estimate of the inclination-corrected scale radius of the optical R-band (rest frame 3100 Angstrom) continuum emission region of log(r_s/cm) = 15.6+-0.3.

We use the Vista Variables in the Via Lactea (VVV) ESO public survey data to measure extinction values in the complete area of the Galactic bulge covered by the survey at high resolution. We derive reddening values using the method described in Paper I. This is based on measuring the mean (J-Ks) color of red clump giants in small subfields of 2' to 6' in the following bulge area: -10.3<b<+5.1 and -10<l<+10.4. To determine the reddening values E(J-Ks) for each region, we measure the RC color and compare it to the (J-Ks) color of RC stars measured in Baade's window, for which we adopt E(B-V)=0.55. This allows us to construct a reddening map sensitive to small scale variations minimizing the problems arising from differential extinction. The significant reddening variations are clearly observed on spatial scales as small as 2'. We find a good agreement between our extinction measurements and Schlegel maps in the outer bulge, but, as already stated in the literature the Schlegel maps are not reliable for regions within |b| < 6. In the inner regions we compare our results with maps derived from DENIS and Spitzer surveys. While we find good agreement with other studies in the corresponding overlapping regions, our extinction map has better quality due to both higher resolution and a more complete spatial coverage in the Bulge. We investigate the importance of differential reddening and demonstrate the need for high resolution extinction maps for detailed studies of Bulge stellar populations and structure. The extinction variations on scales of up to 2'-6', must be taken into account when analysing the stellar populations of the Bulge.

We propose a new exponential f(R) gravity model with f(R)=(R-\lambda c)e^{\lambda(c/R)^n} and n>3, \lambda\geq 1, c>0 to explain late-time acceleration of the universe. At the high curvature region, the model behaves like the \LambdaCDM model. In the asymptotic future, it reaches a stable de-Sitter spacetime. It is a cosmologically viable model and can evade the local gravity constraints easily. This model share many features with other f(R) dark energy models like Hu-Sawicki model and Exponential gravity model. In it the dark energy equation of state is of an oscillating form and can cross phantom divide line \omega_{de}=-1. In particular, in the parameter range 3< n\leq 4, \lambda \sim 1, the model is most distinguishable from other models. For instance, when n=4, \lambda=1, the dark energy equation of state will cross -1 in the earlier future and has a stronger oscillating form than the other models, the dark energy density in asymptotical future is smaller than the one in the high curvature region. This new model can evade the local gravity tests easily when n>3 and \lambda>1.

Coronal hard X-ray (HXR) and continuum gamma-ray sources associated with the impulsive phase of solar flares have been the subject of renewed interest in recent years. They have been interpreted in terms of thin-target, nonthermal bremsstrahlung emission. This interpretation has led to rather extreme physical requirements in some cases. For example, in one case, essentially all of the electrons in the source must be accelerated to nonthermal energies to account for the coronal HXR source. In other cases, the extremely hard photon spectra of the coronal continuum gamma-ray emission suggest that the low energy cutoff of the electron energy distribution lies in the MeV energy range. Here we consider the role of inverse Compton scattering (ICS) as an alternate emission mechanism in both the ultra- and mildly relativistic regimes. It is known that relativistic electrons are produced during powerful flares; these are capable of up-scattering soft photospheric photons to HXR and gamma-ray energies. Previously overlooked is the fact that mildly relativistic electrons, generally produced in much greater numbers in flares of all sizes, can up-scatter EUV/SXR photons to HXR energies. We also explore ICS on anisotropic electron distributions and show that the resulting emission can be significantly enhanced over an isotropic electron distribution for favorable viewing geometries. We briefly review results from bremsstrahlung emission and reconsider circumstances under which nonthermal bremsstrahlung or ICS would be favored. Finally, we consider a selection of coronal HXR and gamma-ray events and find that in some cases the ICS is a viable alternative emission mechanism.

We present \emph{Chandra} monitoring data for six gravitationally lensed quasars: QJ 0158$-$4325, HE 0435$-$1223, HE 1104$-$1805, SDSS 0924+0219, SDSS 1004+4112, and Q 2237+0305. We detect X-ray microlensing variability in all six lenses with high confidence. We detect energy dependent microlensing in HE 0435$-$1223, SDSS 1004+4112, SDSS 0924+0219 and Q 2237+0305. We present a detailed spectral analysis for each lens, and find that simple power-law models plus Gaussian emission lines give good fits to the spectra. We detect intrinsic spectral variability in two epochs of Q 2237+0305. We detect differential absorption between images in four lenses. We also detect the \feka\ emission line in all six lenses, and the Ni XXVII K$\alpha$ line in two images of Q 2237+0305. The rest frame equivalent widths of the \feka\ lines are measured to be 0.4--1.2 keV, significantly higher than those measured in typical active galactic nuclei of similar X-ray luminosities. This suggests that the \feka\ emission region is more compact or centrally concentrated than the continuum emission region.

We present our long term Chandra X-ray monitoring data for the gravitationally lensed quasar Q2237+0305 with 20 epochs spanning 10 years. We easily detect microlensing variability between the images in the full (0.2--8 keV), soft (0.2--2 keV), and hard (2--8 keV) bands at very high confidence. We also detect, for the first time, chromatic microlensing differences between the soft and hard X-ray bands. The hard X-ray band is more strongly microlensed than the soft band, suggesting that the corona above the accretion disk thought to generate the X-rays has a non-uniform electron distribution, in which the hotter and more energetic electrons occupy more compact regions surrounding the black holes. Both the hard and soft X-ray bands are more strongly microlensed than the optical (restframe UV) emission, indicating that the X-ray emission is more compact than the optical, confirming the microlensing results from other lenses.

We give the lens equation for light deflections caused by point mass condensations in an otherwise spatially homogeneous and flat universe. We assume the signal from a distant source is deflected by a single condensation before it reaches the observer. We call this deflector an embedded lens because the deflecting mass is part of the mean density. The embedded lens equation differs from the conventional lens equation because the deflector mass is not simply an addition to the cosmic mean. We prescribe an iteration scheme to solve this new lens equation and use it to compare our results with standard linear lensing theory. We also compute analytic expressions for the lowest order corrections to image amplifications and distortions caused by incorporating the lensing mass into the mean.

We present the first interferometric observation of a zebra-pattern radio burst with simultaneous high spectral (~ 1 MHz) and high time (20 ms) resolution. The Frequency-Agile Solar Radiotelescope (FASR) Subsystem Testbed (FST) and the Owens Valley Solar Array (OVSA) were used in parallel to observe the X1.5 flare on 14 December 2006. By using OVSA to calibrate the FST the source position of the zebra pattern can be located on the solar disk. With the help of multi-wavelength observations and a nonlinear force-free field (NLFFF) extrapolation, the zebra source is explored in relation to the magnetic field configuration. New constraints are placed on the source size and position as a function of frequency and time. We conclude that the zebra burst is consistent with a double-plasma resonance (DPR) model in which the radio emission occurs in resonance layers where the upper hybrid frequency is harmonically related to the electron cyclotron frequency in a coronal magnetic loop.

We compute time delays for gravitational lensing in a flat LambdaCDM Swiss cheese universe. We assume a primary and secondary pair of light rays are deflected by a single point mass condensation described by a Kottler metric (Schwarzschild with Lambda) embedded in an otherwise homogeneous cosmology. We find that the cosmological constant's effect on the difference in arrival times is non-linear and at most around 0.002% for a large cluster lens; however, we find differences from time delays predicted by conventional linear lensing theory that can reach ~4% for these large lenses. The differences in predicted delay times are due to the failure of conventional lensing to incorporate the lensing mass into the mean mass density of the universe.

We compute the deflection angle to order (m/r_0)^2 and m/r_0*Lambda r_0^2 for a light ray traveling in a flat LambdaCDM cosmology which encounters a completely condensed mass region. We use a Swiss cheese model for the inhomogeneities and find that the most significant correction to the Einstein angle occurs not because of the non-linear terms but instead occurs because the condensed mass is embedded in a background cosmology.The Swiss cheese model predicts a decrease in the deflection angle of ~2% for weakly lensed galaxies behind the rich cluster A1689, and that the reduction can be as large as ~5% for similar rich clusters at z\approx 1. Weak lensing deflection angles caused by galaxies can likewise be reduced by as much as ~4%. We show that the lowest order correction in which Lambda appears is proportional to m/r_0*\sqrt{Lambda r_0^2} and could cause as much as a ~0.02% increase in the deflection angle for light that passes through a rich cluster. The lowest order non-linear correction in the mass is proportional to m/r_0*\sqrt{m/r_0} and can increase the deflection angle by ~0.005% for weak lensing by galaxies.

Using the methods of general relativity Lindquist derived the radiative transfer equation that is correct to all orders in v/c. Mihalas developed a method of solution for the important case of monotonic velocity fields with spherically symmetry. We have developed the generalized atmosphere code PHOENIX, which in 1-D has used the framework of Mihalas to solve the radiative transfer equation (RTE) in 1-D moving flows. We describe our recent work including 3-D radiation transfer in PHOENIX and particularly including moving flows exactly using a novel affine method. We briefly discuss quantitative spectroscopy in supernovae.

In this paper we investigate the physical implications of the dynamical scalar mode in pure Horava-Lifshitz gravity on cosmology. We find that it can produce a scale-invariant power spectrum in UV era if the detailed balance condition on the action is relaxed. This indicates that the physical scalar mode may seed the large scale structure and Horava-Lifshitz cosmology could be a qualified alternative to inflationary scenario.

We show that it is possible to equate the intensity reduction of a light wave caused by weak absorption with a geometrical reduction in intensity caused by a "transverse" conformal transformation of the spacetime metric in which the wave travels. We are consequently able to modify Gordon's optical metric to account for electromagnetic properties of ponderable material whose properties include both refraction and absorption. Unlike refraction alone however, including absorption requires a modification of the optical metric that depends on the eikonal of the wave itself. We derive the distance-redshift relation from the modified optical metric for Friedman-Lema\^itre-Robertson-Walker spacetimes whose cosmic fluid has associated refraction and absorption coefficients. We then fit the current supernovae data and provide an alternate explanation (other than dark energy) of the apparent acceleration of the universe.

We show that Gordon's optical metric on a curved spacetime can be generalized to include absorption by allowing the metric to become complex. We demonstrate its use in the realm of geometrical optics by giving three simple examples. We use one of these examples to compute corrected distance-redshift relations for Friedman-Lema\^itre-Robertson-Walker models in which the cosmic fluid has an associated complex index of refraction that represents grey extinction. We then fit this corrected Hubble curve to the type Ia supernovae data and provide a possible explanation (other than dark energy) of the deviation of these observations from dark matter predictions.

Observations and theoretical calculations have shown the importance of non-spherically symmetric structures in supernovae. Thus, the interpretation of observed supernova spectra requires the ability to solve the transfer equation in 3-D moving atmospheres. We present an implementation of the solution of the radiative transfer equation in 3-D homologously expanding atmospheres in spherical coordinates. The implementation is exact to all orders in v/c. We use a new affine method that makes use of the fact that photons travel on straight lines. We compare our results in 3-D for spherically symmetric test problems with high velocity fields and find excellent agreement. Our well-tested 1-D results are based on methods where the momentum directions vary along the characteristic (co-moving momentum directions). Thus, we are able to verify both the analytic framework and its numerical implementation. Additionally, we have been able to test the parallelization over characteristics. Using 512^2 momentum angles we ran the code on 16,384 Opteron processors and achieved excellent scaling. It is now possible to calculate synthetic spectra from realistic 3D hydro simulations. This should open an era of progress in hydro modeling, similar to that that occurred in the 1980s when 1-D models were confronted with synthetic spectra.

We review Gordon's optical metric and the transport equations for the amplitude and polarization of a geometrical optics wave traveling in a gravity field. We apply the theory to the FLRW cosmologies by associating a refraction index with the cosmic fluid. We then derive an expression for the accumulated effect of a refraction index on the distance redshift relations and fit the Hubble curve of current supernova observations with a non-accelerating cosmological model. We also show that some observational effects caused by inhomogeneities, e.g. the Sachs-Wolfe effect, can be interpreted as being caused by an effective index of refraction, and hence this theory could extend to other speed of light communications such as gravitational radiation and neutrino fluxes.

In this paper, we discuss the non-Guassianity originated from the $\alpha$-vacuum on the CMB anisotropy. For $\alpha$-vacuum, there exist correlation between points in the acausal two patches of de Sitter spactime. Such kind of correlation can lead to large local form non-Guassianity in $\alpha$-vacuum. For the single field slow-roll inflationary scenario, the spacetime is in a quasi-de Sitter phase during the inflation. We will show that the $\alpha$-vacuum in this case will lead to non-Gaussianity with distinguished feature, of a large local form and a very different shape.

In this paper, we develop a simple method for the reconstruction of the string-inspired dark energy model with the lagrangian $ {\cal L} = -V (\phi) \sqrt{1 - \alpha' {\nabla}_{\mu} \phi \nabla^{\mu} \phi + \beta' \phi \Box \phi} $ given by Cai et. al., which may allow the equation-of-state parameter cross the cosmological constant boundary $ (w=-1) $. We reconstruct this model in the light of three forms of parametrization for dynamical dark energy.

We investigate the noncommutative effect on the non-Gaussianities of primordial cosmological perturbation. In the lowest order of string length and slow-roll parameter, we find that in the models with small speed of sound the noncommutative modifications could be observable if assuming a relatively low string scale. In particular, the dominant modification of non-Gaussianity estimator f_{NL} could reach O(1) in DBI inflation and K-inflation. The corrections are sensitive to the speed of sound and the choice of string length scale. Moreover the shapes of the corrected non-Gaussianities are distinct from that of ordinary ones.

We calculate the contribution of the fluctuations with the thermal origin to the inflationary nonGaussianity. We find that even a small component of radiation can lead to a large nonGaussianity. We show that this thermal nonGaussianity always has positive $f_{\rm NL}$. We illustrate our result in the chain inflation model and the very weakly dissipative warm inflation model. We show that $f_{NL}\sim {\cal O}(1)$ is general in such models. If we allow modified equation of state, or some decoupling effects, the large thermal nonGaussianity of order $f_{\rm NL}>5$ or even $f_{\rm NL}\sim 100$ can be produced. We also show that the power spectrum of chain inflation should have a thermal origin. In the Appendix A, we made a clarification on the different conventions used in the literature related to the calculation of $f_{\rm NL}$.

We investigate the corrections of the holographic dark energy to inflation paradigm. We study the evolution of the holographic dark energy in the inflationary universe in detail, and carry out a model-independent analysis on the holographic dark energy correction to the primordial scalar power spectrum. It turns out that the corrections generically make the spectrum redder. To be consistent with the experimental data, there must be a upper bound on the reheating temperature. We also discuss the corrections due to different choices of the infrared cutoff.

We derive the radiative transfer equation for arbitrary stationary relativistic flows in stationary spacetimes, i.e. for steady-state transfer problems. We show how the standard characteristics method of solution developed by Mihalas and used throughout the radiative transfer community can be adapted to multi-dimensional applications with isotropic sources. Because the characteristics always coincide with geodesics and can always be specified by constants, direct integration of the characteristics derived from the transfer equation as commonly done in 1-D applications is not required. The characteristics are known for a specified metric from the geodesics. We give details in both flat and static spherically symmetric spacetimes. This work has direct application in 3-dimensional simulations of supernovae, gamma-ray bursts, and active galactic nuclei, as well as in modeling neutron star atmospheres.

We study a class of generalized inflation models in which the inflaton is coupled to the Ricci scalar by a general $f(\phi, R)$ term. The scalar power spectrum, the spectral index, the running of the spectral index, the tensor mode spectrum and a new consistency relation of the model are calculated. We discuss in detail the issues of how to diagonize the coupled perturbation equations, how to deal with an entropy-like source, and how to determine the initial condition by quantization. By studying some explicit models, we find that rich phenomena such as a blue scalar power spectrum, a large running of the spectral index, and a blue tensor mode spectrum can be obtained.

In the latest in our series of papers on XMM-Newton and ground-based optical follow-up of new candidate magnetic cataclysmic variables (mCVs) found in the Sloan Digital Sky Survey, we report classifications of three systems: SDSS J144659.95+025330.3, SDSS J205017.84-053626.8, and SDSS J210131.26+105251.5. Both the X-ray and optical fluxes of SDSS J1446+02 are modulated on a period of 48.7+/-0.5 min, with the X-ray modulation showing the characteristic energy dependence of photo-electric absorption seen in many intermediate polars (IP). A longer period modulation and radial velocity variation is also seen at a period around 4 hrs, though neither dataset set is long enough to constrain this longer, likely orbital, period well. SDSS J2050-05 appears to be an example of the most highly magnetized class of mCV, a disk-less, stream-fed polar. Its 1.57 hr orbital period is well-constrained via optical eclipse timings; in the X-ray it shows both eclipses and an underlying strong, smooth modulation. In this case, broadly phase-resolved spectral fits indicate that this change in flux is the result of a varying normalization of the dominant component (a 41 keV MEKAL), plus the addition of a partial covering absorber during the lower flux interval. SDSS J2101+10 is a more perplexing system to categorize: its X-ray and optical fluxes exhibit no large periodic modulations; there are only barely detectable changes in the velocity structure of its optical emission lines; the X-ray spectra require only absorption by the interstellar medium; and the temperatures of the MEKAL fits are low, with maximum temperature components of either 10 or 25 keV. We conclude that SDSS J2101+10 can not be an IP, nor likely a polar, but is rather most likely a disc accretor-- a low inclination SW Sex star.

A description of the light curves of microlensed Type Ia supernovae (SNe Ia) as extended and expanding sources is presented. We give examples of what microlensing by stellar-mass Chang-Refsdal lenses can do to a small percentage of supernova light curves. We find that in addition to overall brightening, significant changes in light--curve shapes can also occur. Peaks can be distorted, plateaus can appear, and even secondary peaks can be formed. The effects of both the relative motion of the lens and the supernova and the expansion of the supernova are given and compared. The effects of relative motion are more pronounced when a distant supernova ($z_s \sim 1$) impacts well within the Einstein ring of a nearby microlens ($z_d \sim 0.05$) and are less important for more distant deflectors. We also find that the increase in shear that comes with increased deflector distance tends to reduce the time variability of microlensing. We briefly discuss the probability of observing these effects.

We report follow-up XMM-Newton and ground-based optical observations of the unusual X-ray binary SDSS J102347.67+003841.2 (=FIRST J102347.6+003841), and a new candidate intermediate polar found in the Sloan Digital Sky Survey: SDSS J093249.57+472523.0. SDSS J1023 was observed in its low-state, with similar magnitude/color (V=17.4 and B=17.9), and smooth orbital modulation as seen in most previous observations. We further refine the ephemeris (for photometric minimum) to: HJD(TT)_min= 2453081.8546(3) + E* 0.198094(1) d. It is easily detected in X-rays at an unabsorbed flux (0.01-10.0 keV) of 5x10e-13 erg/cm^2/s. Fitting a variety of models we find that: (i) either a hot (kT>~15 keV) optically thin plasma emission model (bremsstrahlung or MEKAL) or a simple power law can provide adequate fits to the data; (ii) these models prefer a low column density ~10e19 cm^-2; (iii) a neutron star atmosphere plus power law model (as found for quiescent low-mass X-ray binaries) can also produce a good fit (for plausible distances), though only for a much higher column of about 4x10e20 cm^-2 and a very cool atmosphere kT<~50eV. These results support the case that SDSS J1023 is a transient LMXB, and indeed places it in the subclass of such systems whose quiescent X-ray emission is dominated by a hard power law component. Our optical photometry of SDSS J0932 reveals that it is an high inclination eclipsing system. Combined with its optical characteristics -- high excitation emission lines, and brightness, yielding a large F_X/F_opt ratio -- its highly absorbed X-ray spectrum argues that SDSS J0932 is a strong IP candidate. However, only more extensive optical photometry and a detection of its spin or spin-orbit beat frequency can confirm this classification. (abridged)

We analyze the existing solar neutrino experiment data and show the allowed regions. The result from SNO's salt phase itself restricts quite a lot the allowed region's area. Reactor neutrinos play an important role in determining oscillation parameters. KamLAND gives decisive conclusion on the solution to the solar neutrino puzzle, in particular, the spectral distortion in the 766.3 Ty KamLAND data gives another new improvement in the constraint of solar MSW-LMA solutions. We confirm that at 99.73% C.L. the high-LMA solution is excluded.

We report follow-up XMM-Newton and optical observations of three new polars found in the Sloan Digital Sky Survey. Simple modeling of the X-ray spectra, and consideration of the details of the X-ray and optical lightcurves corroborate the polar nature of these three systems and provide further insights into their accretion characteristics. During the XMM-Newton observation of SDSS J072910.68+365838.3, X-rays are undetected apart from a probable flare event, during which we find both the typical hard X-ray bremsstrahlung component and a very strong line O VII (E=0.57 keV), but no evidence of a soft blackbody contribution. In SDSS J075240.45+362823.2 we identify an X-ray eclipse at the beginning of the observation, roughly in phase with the primary minimum of the optical broad band curve. The X-ray spectra require the presence of both hard and soft X-ray components, with their luminosity ratio consistent with that found in other recent XMM-Newton results on polars. Lastly, SDSS J170053.30+400357.6 appears optically as a very typical polar, however its large amplitude optical modulation is 180 degrees out of phase with the variation in our short X-ray lightcurve.

The cataclysmic variable SDSSJ 015543.40+002807.2 is confirmed to be a magnetic system of the AM Herculis type. With an orbital period of 87.13 min, it is also the shortest-period eclipsing Polar known. Monitoring with XMM-Newton finds a high-state light curve dominated by a single X-ray emitting accretion pole located slightly prograde of the secondary star. The hard X-ray spectrum is typical of radial shocks on magnetic white dwarfs (kT ~ 10keV), and there is evidence for a soft X-ray component consistent with reprocessing from the stellar surface. The optical circular polarization is weak (v < 3%) when the accretion rate is high (mV ~ 15.5), due to optically-thick cyclotron emission and the apparent competition between two accreting poles. However, in low states (mV ~ 18), the polarization increases smoothly to the blue, reaching 20% at 4200A, and the flux spectrum displays a rich set of thermally-broadened cyclotron harmonics that indicate a polar field of 29MG. The phase interval preceding the 6.5 min eclipse depicts the development of P-Cygni components followed by complete absorption reversals in the emission lines. This phenomenon is not unexpected for a strongly accreting magnetic system viewed through the cool base of the funnel, and high-quality spectroscopy through this interval will likely lead to new insights into the dynamics of magnetic coupling and gas flow onto the white dwarf.

XMM-Newton observations of the polar SDSSJ155331.12+551614.5 reveal that all the X-ray flux emerges at energies less than 2 keV. The best fit to the spectrum is with a thermal plasma with kT=0.8 keV plus a 20-90 ev black body, yielding a thermal X-ray luminosity of 8-9.5E28 ergs/s. The low temperature and X-ray luminosity, together with the lack of variation of the X-ray flux during the observations, are all consistent with an extremely low accretion rate that puts the system in the bombardment regime of accretion, rather than accretion involving a standoff shock. It is likely that the observed X-rays originate from the M dwarf secondary star, thus providing a base activity level for late main sequence stars in close binaries. SDSSJ132411.57+032050.5 is detected by XMM-Newton at the faint EPIC pn count rate of 0.0012 c/s, giving an upper limit to the X-ray luminosity of 7E28 ergs/s for a distance of 300 pc, which is also consistent with the above scenario.

ESA's large X-ray space observatory XMM-Newton is in its fifth year of operations. We give a general overview of the status of calibration of the five X-ray instruments and the Optical Monitor. A main point of interest in the last year became the cross-calibration between the instruments. A cross-calibration campaign started at the XMM-Newton Science Operation Centre at the European Space Astronomy Centre in collaboration with the Instrument Principle Investigators provides a first systematic comparison of the X-ray instruments EPIC and RGS for various kind of sources making also an initial assessment in cross calibration with other X-ray observatories.

We report on exploratory Chandra observations of the three highest redshift quasars known (z = 5.82, 5.99, and 6.28), all found in the Sloan Digital Sky Survey. These data, combined with a previous XMM-Newton observation of a z = 5.74 quasar, form a complete set of color-selected, z > 5.7 quasars. X-ray emission is detected from all of the quasars at levels that indicate that the X-ray to optical flux ratios of z ~ 6 optically selected quasars are similar to those of lower redshift quasars. The observations demonstrate that it will be feasible to obtain quality X-ray spectra of z ~ 6 quasars with current and future X-ray missions.

We present the 158 standard stars that define the u'g'r'i'z' photometric system. These stars form the basis for the photometric calibration of the Sloan Digital Sky Survey (SDSS). The defining instrument system and filters, the observing process, the reduction techniques, and the software used to create the stellar network are all described. We briefly discuss the history of the star selection process, the derivation of a set of transformation equations for the UBVRcIc system, and plans for future work.

We present multi-instrument optical observations of the High Energy Transient Explorer (HETE-2)/Interplanetary Network (IPN) error box of GRB 010921. This event was the first gamma ray burst (GRB) localized by HETE-2 which has resulted in the detection of an optical afterglow. In this paper we report the earliest known observations of the GRB010921 field, taken with the 0.11-m Livermore Optical Transient Imaging System (LOTIS) telescope, and the earliest known detection of the GRB010921 optical afterglow, using the 0.5-m Sloan Digital Sky Survey Photometric Telescope (SDSS PT). Observations with the LOTIS telescope began during a routine sky patrol 52 minutes after the burst. Observations were made with the SDSS PT, the 0.6-m Super-LOTIS telescope, and the 1.34-m Tautenburg Schmidt telescope at 21.3, 21.8, and 37.5 hours after the GRB, respectively. In addition, the host galaxy was observed with the USNOFS 1.0-m telescope 56 days after the burst. We find that at later times (t > 1 day after the burst), the optical afterglow exhibited a power-law decline with a slope of $\alpha = 1.75 \pm 0.28$. However, our earliest observations show that this power-law decline can not have extended to early times (t < 0.035 day).

The discovery of an optical counterpart to GRB010222 (detected by BeppoSAX; Piro 2001) was announced 4.4 hrs after the burst by Henden (2001a). The Sloan Digital Sky Survey's 0.5m photometric telescope (PT) and 2.5m survey telescope were used to observe the afterglow of GRB010222 starting 4.8 hours after the GRB. The 0.5m PT observed the afterglow in five, 300 sec g' band exposures over the course of half an hour, measuring a temporal decay rate in this short period of F_nu \propto t^{-1.0+/-0.5}. The 2.5m camera imaged the counterpart nearly simultaneously in five filters (u' g' r' i' z'), with r' = 18.74+/-0.02 at 12:10 UT. These multicolor observations, corrected for reddening and the afterglow's temporal decay, are well fit by the power-law F_nu \propto nu^{-0.90+/-0.03} with the exception of the u' band UV flux which is 20% below this slope. We examine possible interpretations of this spectral shape, including source extinction in a star forming region.

Early data taken during commissioning of the SDSS have resulted in the discovery of a very cool white dwarf. It appears to have stronger collision induced absorption from molecular hydrogen than any other known white dwarf, suggesting it has a cooler temperature than any other. While its distance is presently unknown, it has a surprisingly small proper motion, making it unlikely to be a halo star. An analysis of white dwarf cooling times suggests that this object may be a low-mass star with a helium core. The SDSS imaging and spectroscopy also recovered LHS 3250, the coolest previously known white dwarf, indicating that the SDSS will be an effective tool for identifying these extreme objects.

A sample of 4208 objects with magnitude 15 < g* < 22 and colors of main sequence A stars has been selected from 370 square degrees of Sloan Digital Sky Survey (SDSS) commissioning observations. The data is from two long, narrow stripes, each with an opening angle of greater than 60 deg, at Galactic latitudes 36 < abs(b) < 63 on the celestial equator. An examination of the sample's distribution shows that these stars trace considerable substructure in the halo. Large overdensities of A-colored stars in the North at (l,b,R) = (350, 50, 46 kpc) and in the South at (157, -58, 33 kpc) and extending over tens of degrees are present in the halo of the Milky Way. Using photometry to separate the stars by surface gravity, both structures are shown to contain a sequence of low surface gravity stars consistent with identification as a blue horizontal branch (BHB). Both structures also contain a population of high surface gravity stars two magnitudes fainter than the BHB stars, consistent with their identification as blue stragglers (BSs). From the numbers of detected BHB stars, lower limits to the implied mass of the structures are 6x10^6 M_sun and 2x10^6 M_sun. The fact that two such large clumps have been detected in a survey of only 1% of the sky indicates that such structures are not uncommon in the halo. Simple spheroidal parameters are fit to a complete sample of the remaining unclumped BHB stars and yield (at r < 40 kpc) a fit to a halo distribution with flattening (c/a = 0.65+/-0.2) and a density falloff exponent of alpha = -3.2+/-0.3.

Interstellar absorption in the galactic plane is highly variable from one direction to another. In this paper colour excesses and distances from a new open cluster sample are used to investigate the spatial distribution of the interstellar extinction. An inverse method (Tarantola & Valette, 1982) is used to construct the extinction map in the galactic plane below $|b| < 10^{o}$. The $A_{v} (r,l) $ diagrams are compared with those derived from individual stars (Arenou et al. 1992, Neckel & Klare 1980). An analytic expression for the interstellar extinction as a function of galactic longitude and distance in the solar neighborhood is given. The comparison of the model predictions with Hipparcos observations in the 4-dimensional space of ($V$, $B-V$, $H_v$, $r$) shows that our extinction model provides a better fit to the data. However, a new and more detailed extinction model is still lacking.

NGC 4815 is a distant and populous open cluster, which lies in the galactic plane in a region of strong absorption. As a consequence, its membership, spatial distribution and luminosity function are not well determined. In this paper, we present an algorithm which uses both positional and photometric data simultaneously to search for open cluster members. The contribution of the field stars is estimated by our Galaxy model software. The method has been developed for deep CCD photometric observations, but it can also be used to investigate the membership in an n-dimensional space, including proper motions, radial velocity and metallicity. A clean and well defined colour-magnitude diagram (CMD) for NGC 4815 has been obtained after decontaminating from field stars. We have investigated the spatial distribution of cluster members and the luminosity function. We have shown that the projected surface distribution of stars is well-represented by an exponentially decreasing function with a scale length of 1.67 $\pm$ 0.06 arcmin. We have found evidence for mass segregation. The luminosity function for NGC 4815 is similar to that of the Hyades cluster, and shows a gap at $V$ = 15 mag.