We present the results of our systematic search for the reddest far-infrared (FIR) and submillimeter (sub-mm) galaxies using the data from the Herschel Multi-tiered Extragalactic Survey (HerMES) and the SCUBA2 Cosmological Legacy Survey (S2CLS). The red FIR galaxies are "500~$\mu$m risers," whose spectral energy distributions increase with wavelength across the three FIR passbands of the Spectral and Photometric Imaging REceiver (SPIRE) of Herschel. Within 106.5 deg$^2$ of the HerMES fields, we have selected 629 500 $\mu$m risers. The red sub-mm galaxies are "SPIRE dropouts," which are prominent detections in the S2CLS 850 $\mu$m data but are extremely weak or invisible in the SPIRE bands. Within the 2.98 deg$^2$ common area of HerMES and S2CLS, we have selected 95 such objects. These very red sources could be dusty starbursts at high redshifts ($z\gtrsim 4$-6) because the peak of their cold-dust emission heated by star formation is shifted to the reddest FIR/sub-mm bands. The surface density of 500 $\mu$m risers is $\sim$8.2 deg$^{-2}$ at the $\geq 20$ mJy level in 500 $\mu$m, while that of SPIRE dropouts is $\sim$19.3 deg$^{-2}$ at the $\geq 5$ mJy level in 850 $\mu$m. Both type of objects could span a wide range of redshifts, however. Using deep radio data in these fields to further select the ones likely at the highest redshifts, we find that the surface density of $z>6$ candidates is 5.5 deg$^{-2}$ among 500 $\mu$m risers and is 0.8-13.6 deg$^{-2}$ among SPIRE dropouts. If this is correct, the dust-embedded star formation processes in such objects could contribute comparably as Lyman-break galaxies to the global SFR density at $z>6$.

We study properties of magnetohydrodynamic (MHD) eigenmodes by decomposing the data of MHD simulations into linear MHD modes - namely the Alfven, slow magnetosonic, and fast magnetosonic modes. We drive turbulence with a mixture of solenoidal and compressive driving, while varying the Alfven Mach number (MA), plasma beta, and the sonic Mach number from sub-sonic to trans-sonic. We find that the proportion of fast and slow modes in the mode mixture increases with increasing compressive forcing. This proportion of the magnetosonic modes can also become the dominant fraction in the mode mixture. The anisotropy of the modes is analyzed by means of their structure functions. The Alfven mode anisotropy is consistent with the Goldreich-Sridhar theory. We find a transition from weak to strong Alfvenic turbulence as we go from low to high MA. The slow mode properties are similar to the Alfven mode. On the other hand the isotropic nature of fast modes is verified in the cases where the fast mode is a significant fraction of the mode mixture. The fast mode behavior does not show any transition in going from low to high MA. We find indications that there is some interaction between the different modes and the properties of the dominant mode can affect the properties of the weaker modes. This work identifies the conditions under which magnetosonic modes can be a major fraction of turbulent astrophysical plasmas, including the regime of weak turbulence. Important astrophysical implications for cosmic ray transport and magnetic reconnection are discussed.

Since September 2018, LAMOST starts a new 5-year medium-resolution spectroscopic survey (MRS) using bright/gray nights. We present the scientific goals of LAMOST-MRS and propose a near optimistic strategy of the survey. A complete footprint is also provided. Not only the regular medium-resolution survey, but also a time-domain spectroscopic survey is being conducted since 2018 and will be end in 2023. According to the detailed survey plan, we expect that LAMOST-MRS can observe about 2 million stellar spectra with ~7500 and limiting magnitude of around G=15 mag. Moreover, it will also provide about 200 thousand stars with averagely 60-epoch observations and limiting magnitude of G~14 mag. These high quality spectra will give around 20 elemental abundances, rotational velocities, emission line profiles as well as precise radial velocity with uncertainty less than 1 km/s. With these data, we expect that LAMOST can effectively leverage sciences on stellar physics, e.g. exotic binary stars, detailed observation of many types of variable stars etc., planet host stars, emission nebulae, open clusters, young pre-main-sequence stars etc.

Deep observation of the High Energy Stereoscopic System (HESS) on the most extended pulsar wind nebula HESS J1825-137 reveals an enhanced energy-dependent morphology, providing useful information on the particle transport mechanism in the nebula. We find that the energy-dependent morphology is consistent with a diffusion-dominated transport of electrons/positrons. It provides an alternative possible interpretation for the unusually large spatial extent (i.e., $\gtrsim 100$pc) of the nebula, which could then be attributed to the diffusion of escaping electrons/positrons from a compact plerion. The influence of various model parameters on the energy-dependent extent of the nebula is studied in the diffusion-dominated scenario. We also show that the energy-dependent morphology of the nebula may also be used to study the spin-down history of the pulsar.

We report velocity-delay maps for prominent broad emission lines, Ly_alpha, CIV, HeII and H_beta, in the spectrum of NGC5548. The emission-line responses inhabit the interior of a virial envelope. The velocity-delay maps reveal stratified ionization structure. The HeII response inside 5-10 light-days has a broad single-peaked velocity profile. The Ly_alpha, CIV, and H_beta responses peak inside 10 light-days, extend outside 20 light-days, and exhibit a velocity profile with two peaks separated by 5000 km/s in the 10 to 20 light-day delay range. The velocity-delay maps show that the M-shaped lag vs velocity structure found in previous cross-correlation analysis is the signature of a Keplerian disk with a well-defined outer edge at R=20 light-days. The outer wings of the M arise from the virial envelope, and the U-shaped interior of the M is the lower half of an ellipse in the velocity-delay plane. The far-side response is weaker than that from the near side, so that we see clearly the lower half, but only faintly the upper half, of the velocity--delay ellipse. The delay tau=(R/c)(1-sin(i))=5 light-days at line center is from the near edge of the inclined ring, giving the inclination i=45 deg. A black hole mass of M=7x10^7 Msun is consistent with the velocity-delay structure. A barber-pole pattern with stripes moving from red to blue across the CIV and possibly Ly_alpha line profiles suggests the presence of azimuthal structure rotating around the far side of the broad-line region and may be the signature of precession or orbital motion of structures in the inner disk. Further HST observations of NGC 5548 over a multi-year timespan but with a cadence of perhaps 10 days rather than 1 day could help to clarify the nature of this new AGN phenomenon.

Interesting chemically peculiar field stars may reflect their stellar evolution history and their possible origin in a different environment from where they are found now, which is one of the most important research fields in Galactic archaeology. To explore this further, we have used the CN-CH bands around 4000 A to identify N-rich metal-poor field stars in LAMOST DR3. Here we expand our N-rich metal-poor field star sample to ~100 stars in LAMOST DR5, where 53 of them are newly found in this work. We investigate light elements of the common stars between our sample and APOGEE DR14. While Mg, Al, and Si abundances generally agree with the hypothesis that N-rich metal-poor field stars come from enriched populations in globular clusters, it is still inconclusive for C, N, and O. After integrating the orbits of our N-rich field stars and a control sample of normal metal-poor field stars, we find that N-rich field stars have different orbital parameter distributions compared to the control sample, specifically, apocentric distances, maximum vertical amplitude (Zmax), orbital energy, and z direction angular momentum (Lz). The orbital parameters of N-rich field stars indicate that most of them are inner-halo stars. The kinematics of N-rich field stars support their possible GC origin. The spatial and velocity distributions of our bona fide N-rich field star sample are important observational evidence to constrain simulations of the origin of these interesting objects.

Over the last two decades, the Submillimetre Common-User Bolometer Array (SCUBA) and SCUBA-2 on the James Clerk Maxwell Telescope (JCMT) achieved gread success in discovering the population of dusty starburst galaxies in the early universe. The SCUBA-2 surveys at 450 micron and 850 micron set important constraints on the obscured star formation over cosmic time, and in combination of deep optical and near-IR data, allows the study of protoclusters and structure formation. However, the current submillimeter (submm) surveys by JCMT are still limited by area of sky coverage (confusion limit mapping of only a few deg^2), which prevent a systematic study of large samples of the obscured galaxy population. In this white paper, we review the studies of the submm galaxies with current submillimeter/millimeter (submm/mm) observations, and discuss the important science with the new submm instruments in the next decade. In particular, with a 10 times faster mapping speed of the new camera, we will expect deep 850 micron surveys over 10 to 100 times larger sky area to i) largely increase the sample size of submm detections toward the highest redshift, ii) improve our knowledge of galaxy and structure formation in the early universe.

Standard stellar evolution model predicts a severe depletion of lithium (Li) abundance during the first dredge up process (FDU). Yet a small fraction of giant stars are still found to preserve a considerable amount of Li in their atmospheres after the FDU. Those giants are usually identified as Li-rich by a widely used criterion, A(Li) > 1.5 dex. A large number of works dedicated to searching for investigating this minority of the giant family, and the amount of Li-rich giants, has been largely expanded on, especially in the era of big data. In this paper, we present a catalog of Li-rich giants found from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) survey with Li abundances derived from a template-matching method developed for LAMOST low-resolution spectra.

Based on the medium-high resolution (R~ 20,000), modest signal-to-noise ratio (S/N > 70) FLAMES-GIRAFFE spectra, we investigated the copper abundances of 129 red giant branch stars in the Galactic bulge with [Fe/H] from -1.14 to 0.46 dex. The copper abundances are derived from both local thermodynamic equilibrium (LTE) and nonlocal thermodynamic equilibrium (NLTE) with the spectral synthesis method. We find that the NLTE effects for Cu I lines show a clear dependence on metallicity, and they gradually increase with decreasing [Fe/H] for our sample stars. Our results indicate that the NLTE effects of copper are important not only for metal-poor stars but also for supersolar metal-rich ones and the LTE results underestimate the Cu abundances. We note that the [Cu/Fe] trend of the bulge stars is similar to that of the Galactic disk stars spanning the metallicity range of -1.14 < [Fe/H] < 0.0 dex and the [Cu/Fe] ratios increase with increasing metallicity when [Fe/H] is from~-1.2 to~-0.5 dex, favoring a secondary (metallicity-dependent) production of Cu.

In this letter we propose that the X-ray and the TeV observations in the vicinity of Geminga can be understood in the framework of anisotropic diffusion of injected electrons/positrons. This interpretation only requires the turbulence in the vicinity of Geminga to be sub-Alfv\'enic with the local mean magnetic field direction approximately aligned with our line of sight towards Geminga, without invoking extreme conditions for the environment, such as an extremely small diffusion coefficient and a weak magnetic field of strength $<1\mu$G as suggested in previous literature.

We model the ultraviolet spectra of the Seyfert 1 galaxy NGC~5548 obtained with the Hubble Space Telescope during the 6-month reverberation-mapping campaign in 2014. Our model of the emission from NGC 5548 corrects for overlying absorption and deblends the individual emission lines. Using the modeled spectra, we measure the response to continuum variations for the deblended and absorption-corrected individual broad emission lines, the velocity-dependent profiles of Ly$\alpha$ and C IV, and the narrow and broad intrinsic absorption features. We find that the time lags for the corrected emission lines are comparable to those for the original data. The velocity-binned lag profiles of Ly$\alpha$ and C IV have a double-peaked structure indicative of a truncated Keplerian disk. The narrow absorption lines show delayed response to continuum variations corresponding to recombination in gas with a density of $\sim 10^5~\rm cm^{-3}$. The high-ionization narrow absorption lines decorrelate from continuum variations during the same period as the broad emission lines. Analyzing the response of these absorption lines during this period shows that the ionizing flux is diminished in strength relative to the far-ultraviolet continuum. The broad absorption lines associated with the X-ray obscurer decrease in strength during this same time interval. The appearance of X-ray obscuration in $\sim\,2012$ corresponds with an increase in the luminosity of NGC 5548 following an extended low state. We suggest that the obscurer is a disk wind triggered by the brightening of NGC 5548 following the decrease in size of the broad-line region during the preceding low-luminosity state.

The equipartition of magnetic and thermal energy in the interstellar medium (ISM) indicates the magneto-hydrodynamic nature of the interstellar turbulence, which can be decomposed into three wave modes: Alfv\'en, fast and slow magnetosonic modes\cite{Herlofson50,LG01,CL03}. Even for studies performed in the case of subsonic turbulence, the magnetosonic modes have been revealed to be more effective than the Alfv\'enic modes for processes such as cosmic ray (CR) transport and acceleration \cite{YL04,Lynn14}. The multiphase nature of ISM and diversity of driving mechanisms give rise to spatial variation of turbulence properties. Nevertheless, the employed model of magneto-hydrodynamic turbulence is often oversimplified being assumed to be only Alfv\'enic or even hydrodynamic due to a lack of observational evidence. Here we report the employment of our novel method, the signature from polarization analysis (SPA), on unveiling the plasma modes in interstellar turbulence. Its application leads to the first discovery of magnetosonic modes in the Cygnus X region. It is found that the magnetosonic modes overlap to a high degree with Fermi cocoon, consistent with theoretical expectations. Moreover, through comparison with the spectrum at other wavelengths, the plasma modes of turbulence are unveiled in active star formation zones. The SPA casts light on the plasma modes composition of the Galactic turbulence, and marks the onset of a new era in the study of interstellar turbulence and accordingly our understandings of relevant processes including cosmic ray transport and star formation.

The discovery of Li-rich giant has cast a new challenge for the standard stellar evolution models, and to resolve this issue, the number of this type object has been rapidly increased because of the development of worldwide surveys these days. Taking advantage of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope survey, 44 newly Li-rich giants are reported, which are confirmed with high-resolution observations. Based on the high-resolution and high signal-to-noise spectra, we derived the atmospheric parameters and elemental abundances with the spectral synthesis method. We performed a detailed analysis of their evolutionary stages, infrared excess, projected rotational velocity (v sin i), and the stellar population. We find that (1) The Li-rich giants concentrate at the evolutionary status of the red giant branch bump, red clump, and asymptotic giant branch; (2) Three of them are fast rotators and none exhibit infrared excess. Our results imply that the origins of Li enrichment are most likely to be associated with the extra mixing in the stellar interior, and the external sources maybe only make a minor contribution. Moreover, various Li-rich episodes take place at different evolutionary stages.

Magnetic fields have important or dominant effects in many areas of astrophysics, but have been very difficult to quantify. Spectropolarimetry from Ground State Alignment (GSA) has been suggested as a direct tracer of magnetic field in interstellar diffuse medium. Owing to the long life of the atoms on ground states the Larmor precession in an external magnetic field imprints the direction of the field onto the polarization of absorbing species. This provides a unique tool for studies of sub-gauss magnetic fields using polarimetry of UV, optical and radio lines. Many spectral lines with strong signals from GSA are in the UV band. By discerning magnetic fields in gas with different dynamical properties, high spectral resolution measurement of spectral polarization will allow the study of 3D magnetic field distribution and interstellar turbulence. GSA provides also a unique chance to map 3D direction of magnetic field on small scales, e.g., disks, where grain alignment is unreliable. The range of objects suitable for studies is extremely wide and includes magnetic fields in the interplanetary medium, in the interstellar medium, and in circumstellar regions as well as diffuse media in extragalactic objects.

Sub-Gauss magnetic fields are crucial for many physical processes from Galactic $kpc$ scale to stellar $Au$ scale, but are nontrivial to detect. All present magnetic diagnoses trace only one component of a magnetic field at best. Here we report the first observational results that unveil 3D topology as well as the strength of a sub-milliGauss magnetic field beyond our solar system. We found that two weak neutral iron absorption lines from the ground state in the atmosphere of $89~Her$ produced counterintuitive high-amplitude polarizations and consistent polarization angles, exactly in line with the theoretical prediction from ground state alignment (GSA). Our analysis reveals the first sub-AU scale magnetic field on $89~Her$, that is $1.3kpc$ away, has a 3D orientation aligned to the stellar outflow axis and a strength of $70\mu{G} \lesssim B \lesssim150\mu{G}$, thus substantially improving the accuracy by five orders of magnitude compared to the previous $10G$ upper limit set by non-detection of Zeeman effect. This long-awaited discovery is anticipated to usher in a new era of precise measurement of sub-gauss magnetic field in the Universe.

Studying the formation and evolution of galaxies at the earliest cosmic times, and their role in reionization, requires the deepest imaging possible. Ultra-deep surveys like the HUDF and HFF have pushed to mag \mAB$\,\sim\,$30, revealing galaxies at the faint end of the LF to $z$$\,\sim\,$9$\,-\,$11 and constraining their role in reionization. However, a key limitation of these fields is their size, only a few arcminutes (less than a Mpc at these redshifts), too small to probe large-scale environments or clustering properties of these galaxies, crucial for advancing our understanding of reionization. Achieving HUDF-quality depth over areas $\sim$100 times larger becomes possible with a mission like the Wide Field Infrared Survey Telescope (WFIRST), a 2.4-m telescope with similar optical properties to HST, with a field of view of $\sim$1000 arcmin$^2$, $\sim$100$\times$ the area of the HST/ACS HUDF. This whitepaper motivates an Ultra-Deep Field survey with WFIRST, covering $\sim$100$\,-\,$300$\times$ the area of the HUDF, or up to $\sim$1 deg$^2$, to \mAB$\,\sim\,$30, potentially revealing thousands of galaxies and AGN at the faint end of the LF, at or beyond $z$\,$\sim$\,9$\,-\,$10 in the epoch of reionization, and tracing their LSS environments, dramatically increasing the discovery potential at these redshifts. (Note: This paper is a somewhat expanded version of one that was submitted as input to the Astro2020 Decadal Survey, with this version including an Appendix (which exceeded the Astro2020 page limits), describing how the science drivers for a WFIRST Ultra Deep Field might map into a notional observing program, including the filters used and exposure times needed to achieve these depths.)

Recent near-IR power-spectra and panchromatic Extragalactic Background Light measurements provide upper limits on the near-IR surface brightness (SB>31 mag/arcsec^2) that may come from Pop III stars and accretion disks around resulting stellar-mass black holes (BHs) in the epoch of First Light (z=7-17). Physical parameters for zero metallicity Pop III stars at z>7 can be estimated from MESA stellar evolution models through helium-depletion, and for BH accretion disks from quasar microlensing results and multicolor accretion models. Second-generation stars can form at higher multiplicity, so that BH accretion disks may be fed by Roche-lobe overflow from lower-mass companions in their AGB stage. The near-IR SB constraints can be used to calculate the number of caustic transits behind lensing clusters that JWST and the 25~39 m ground-based telescopes may detect for both Pop III stars and stellar mass BH accretion disks. Because Pop III stars and stellar mass BH accretion disks have sizes of a few x 10^-11 arcsec at z>7, typical caustic magnifications can be mu=10^4~10^5, with rise times of hours and decline times of < 1 year for cluster transverse velocities of v<1000 km/s. Microlensing by intracluster medium objects can modify transit magnifications, and lengthen visibility times. Depending on BH masses, accretion-disk radii and feeding efficiencies, stellar-mass BH accretion-disk caustic transits could outnumber those from Pop III stars. To observe Pop III caustic transits directly may require monitoring 3~30 lensing clusters to AB< 29 mag over a decade or more. Such a program must be started with JWST in Cycle 1, and -- depending on the role of microlensing in the Intra Cluster Light -- should be continued for decades with the GMT and TMT, where JWST and the ground-based telescopes each will play a unique and strongly complementary role.

Although many high-energy neutrinos detected by the IceCube telescope are believed to have anextraterrestrial origin, their astrophysical sources remain a mystery. Recently, an unprecedenteddiscovery of a high-energy muon neutrino event coincident with a multiwavelength flare from ablazar, TXS 0506+056, shed some light on the origin of the neutrinos. It is usually believed that ablazar is produced by a relativistic jet launched from an accreting supermassive black hole (SMBH).Here we show that the high-energy neutrino event can be interpreted by the inelastic hadronuclearinteractions between the accelerated cosmic-ray protons in the relativistic jet and the dense gasclouds in the vicinity of the SMBH. Such a scenario only requires a moderate proton power in thejet, which could be much smaller than that required in the conventional hadronic model whichinstead calls upon the photomeson process. Meanwhile, the flux of the multiwavelength flare fromthe optical to gamma-ray band can be well explained by invoking a second radiation zone in thejet at a larger distance to the SMBH. In our model, the neutrino emission lasts a shorter time thanthe multiwavelength flare so the neutrino event is not necessarily correlated with the flare but it is probably accompanied by a spectrum hardening above a few GeV.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

We present Hubble Space Telescope WFC3-IR imaging in the fields of six apparently bright dusty star-forming galaxies (DSFGs) at $z$ = 2-4 identified by their rest-frame far-infrared colors using the Planck and Herschel space facilities. We detect near-infrared counterparts for all six submillimeter sources, allowing us to undertake strong-lensing analyses. One field in particular stands out for its prominent giant arcs, PLCK G165.7+67.0 (G165). After combining the color and morphological information, we identify 11 sets of image multiplicities in this one field. We construct a strong-lensing model constrained by this lensing evidence, which uncovers a bimodal spatial mass distribution, and from which we measure a mass of $(2.6 \pm 0.11)$ $\times$ $10^{14}$ $M_{\odot}$ within $\sim$250 kpc. The bright ($S_{350}$ $\approx$ 750 mJy) DSFG appears as two images: a giant arc with a spatial extent of 4.5" that is merging with the critical curve, and a lower-magnification counterimage that is detected in our new longer-wavelength ground- and space-based imaging data. Using our ground-based spectroscopy, we calculate a dynamical mass of $1.3^{+0.04}_{-0.70} \times 10^{15}$ $M_{\odot}$ to the same fixed radius, although this value may be inflated relative to the true value if the velocity distribution is enhanced in the line-of-sight direction. We suggest that the bimodal mass taken in combination with the weak X-ray flux and low SZ decrement may be explained as a pre-merger for which the intracluster gas is diluted along the line of sight, while the integrated surface mass density is supercritical to strong-lensing effects.

Radial velocity is one of key measurements in understanding the fundamental properties of stars, stellar clusters and the Galaxy. A plate of stars in the Kepler field were observed in May of 2018 with the medium-resolution spectrographs of LAMOST, aiming to test the performance of this new system which is the upgraded equipment of LAMOST after the first five-year regular survey.We present our analysis on the radial velocity measurements (RVs) derived from these data. The results show that slight and significant systematic errors exist among the RVs obtained from the spectra collected by different spectrographs and exposures, respectively. After correcting the systematic errors with different techniques, the precision of RVs reaches ~1.3, ~1.0, ~0.5 and ~0.3 km/s at S/Nr = 10, 20, 50, and 100, respectively. Comparing with the RVs of the standard stars of the APOGEE survey, our RVs are calibrated with a zero-point shift of ~7 km/s. The results indicate that the LAMOST medium-resolution spectroscopic system may provide RVs in a reasonable accuracy and precision for the selected targets.

We present deep $J$ and $H$-band images in the extended Great Observatories Origins Deep Survey-North (GOODS-N) field covering an area of 0.22 $\rm{deg}^{2}$. The observations were taken using WIRCam on the 3.6-m Canada France Hawaii Telescope (CFHT). Together with the reprocessed $K_{\rm s}$-band image, the $5\sigma$ limiting AB magnitudes (in 2" diameter apertures) are 24.7, 24.2, and 24.4 AB mag in the $J$, $H$, and $K_{\rm s}$ bands, respectively. We also release a multi-band photometry and photometric redshift catalog containing 93598 sources. For non-X-ray sources, we obtained a photometric redshift accuracy $\sigma_{\mathrm{NMAD}}=0.036$ with an outlier fraction $\eta = 7.3\%$. For X-ray sources, which are mainly active galactic nuclei (AGNs), we cross-matched our catalog with the updated 2M-CDFN X-ray catalog from Xue et al. (2016) and found that 658 out of 683 X-ray sources have counterparts. $GALEX$ UV data are included in the photometric redshift computation for the X-ray sources to give $\sigma_{\mathrm{NMAD}} = 0.040$ with $\eta=10.5\%$. Our approach yields more accurate photometric redshift estimates compared to previous works in this field. In particular, by adopting AGN-galaxy hybrid templates, our approach delivers photometric redshifts for the X-ray counterparts with fewer outliers compared to the 3D-HST catalog, which fit these sources with galaxy-only templates.

We study the requirement on the jet power in the conventional p-$\gamma$ models (photopion production and Bethe-Heitler pair production) for TeV BL Lac objects. We select a sample of TeV BL Lac objects whose SEDs are difficult to be explained by the one-zone leptonic model. Based on the relation between the p-$\gamma$ interaction efficiency and the opacity of $\gamma\gamma$ absorption, we find that the detection of TeV emission poses upper limits on the p-$\gamma$ interaction efficiencies in these sources and hence minimum jet powers can be derived accordingly. We find that the obtained minimum jet powers exceed the Eddington luminosity of the supermassive black holes. Implications for the accretion mode of the supermassive black hole in these BL Lac objects and the origin of their TeV emissions are discussed.

Various studies have implied the existence of a gaseous halo around the Galaxy extending out to 100 kpc. Galactic cosmic rays (CRs) that propagate to the halo, either by diffusion or by convection with the possibly existing large-scale Galactic wind, can interact with the gas therein and produce gamma-rays via proton-proton collision. We calculate the cosmic ray distribution in the halo and the gamma-ray flux, and explore the dependence of the result on model parameters such as diffusion coefficient, CR luminosity, CR spectral index. We find that the current measurement of isotropic gamma-ray background at $\lesssim$TeV with Fermi Large Area Telescope already approaches a level that can provide interesting constraints on the properties of Galactic cosmic ray (e.g., with CR luminosity $L_{CR}\leq 10^{41}$erg/s). We also discuss the possibilities of the Fermi bubble and IceCube neutrinos originating from the proton-proton collision between cosmic rays and gas in the halo, as well as the implication of our results for the baryon budget of the hot circumgalactic medium of our Galaxy. Given that the isotropic gamma-ray background is likely to be dominated by unresolved extragalactic sources, future telescopes may extract more individual sources from the IGRB, and hence put even more stringent restriction on the relevant quantities (such as Galactic cosmic ray luminosity and baryon budget in the halo) in the presence of a turbulent halo that we consider.

We study the shadows cast by near-extremal Kerr-MOG black holes for different values of the parameter in modified gravity (MOG). In particular, we consider an isotropic emitter orbiting near such black holes and analytically compute the positions, fluxes and redshift factors of their images. The size of the shadow decreases when the modified parameter is increased. For each shadow, the images of the emitter appear on a special part of the shadow which has a rich structure. The primary image and secondary images are similar to those produced for the near-extremal (high spin) Kerr black hole, but the near-extremal Kerr-MOG black hole can have a spin ($\hat{J}/M^2_{\alpha}$) which is finitely lower than 1. When the modified parameter is varied, the typical positions of the corresponding images do not change, nor does the typical redshift factor associated with the primary image. However, another typical redshift factor associated with the secondary image increases when the modified parameter is increased. We also find that the fluxes increase in that case. These images appear periodically with period greater than that of Kerr. This provides an alternative signature away from the Kerr case which may be tested by the Event Horizon Telescope.

Luminous and ultra-luminous infrared galaxies ((U)LIRGs) are rare today but are increasingly abundant at high redshifts. They are believed to be dusty starbursts, and hence should have high rates of supernovae (multiple events per year). Due to their extremely dusty environment, however, such supernovae could only be detected in restframe infrared and longer wavelengths, where our current facilities lack the capability of finding them individually beyond the local universe. We propose a new technique for higher redshifts, which is to search for the presence of supernovae through the variability of the integrated rest-frame infrared light of the entire hosts. We present a pilot study to assess the feasibility of this technique. We exploit a unique region, the "IRAC Dark Field" (IDF), that the Spitzer Space Telescope has observed for more than 14 years in 3--5 micron. The IDF also has deep far-infrared data (200--550 micron) from the Herschel Space Observatory that allow us to select high-redshift (U)LIRGs. We obtain a sample of (U)LIRGs that have secure optical counterparts, and examine their light curves in 3--5 micron. While the variabilities could also be caused by AGNs, we show that such contaminations can be identified. We present two cases where the distinct features in their light curves are consistent with multiple supernovae overlapping in time. Searching for supernovae this way will be relevant to the James Webb Space Telescope (JWST) to probe high-redshift (U)LIRGs into their nuclear regions where JWST will be limited by its resolution.

Chemical abundances in the X-ray halos (also known as the intracluster medium, ICM) of clusters and groups of galaxies can be measured via prominent emission line features in their X-ray spectra. Elemental abundances are footprints of time-integrated yields of various stellar populations that have left their specific abundance patterns prior to and during the cluster and group evolution. We aim to constrain nitrogen abundances in the CHEmical Evolution RGS Sample (CHEERS), which contains 44 nearby groups and clusters of galaxies, in order to have a better understanding of their chemical enrichment. We examine the high-resolution spectra of the CHEERS sample carefully and take into account various systematic effects in the spectral modelling. We compare the observed abundance ratios with those in the Galactic stellar populations, as well as predictions from stellar yields (low- and intermediate-mass stars, massive stars and degenerate stars). The nitrogen abundance can only be well constrained ($\gtrsim3\sigma$) in one cluster of galaxies and seven groups of galaxies. The [O/Fe] -- [Fe/H] relation of the ICM is comparable to that for the Galaxy, while both [N/Fe] and [N/O] ratios of the ICM are higher than in the Galaxy. Future studies on nitrogen radial distributions are required to tell whether the obtained higher [N/Fe] and [N/O] ratios are biased due to the small extraction region ($r/r_{500}\lesssim0.05$) that we adopt here. Since abundances of odd-$Z$ elements are more sensitive to the initial metallicity of stellar populations, accurate abundance measurements of N, Na and Al are required to better constrain the chemical enrichment in the X-ray halos of clusters and groups of galaxies.

About one percent of giants\upcite{Brown1989} are detected to have anomalously high lithium (Li) abundances in their atmospheres, conflicting directly with the prediction of the standard stellar evolution models\upcite{Iben1967}, and making the production and evolution of Li more intriguing, not only in the sense of the Big Bang nucleosynthesis\upcite{Cyburt2016,Spite1982} or the Galactic medium\upcite{Tajitsu2015}, but also the evolution of stars. Decades of efforts have been put into explaining why such outliers exist\upcite{Sackmann1999, Denissenkov2004, Charbonnel2010}, yet the origins of Li-rich giants are still being debated. Here we report the discovery of the most Li-rich giant known to date, with a super-high Li abundance of 4.51. This rare phenomenon was snapshotted together with another short-term event that the star is experiencing its luminosity bump on the red giant branch. Such high Li abundance indicates that the star might be at the very beginning of its Li-rich phase, which provides a great opportunity to investigate the origin and evolution of Li in the Galaxy. A detailed nuclear simulation is presented with up-to-date reaction rates to recreate the Li enriching process in this star. Our results provide tight constraints on both observational and theoretical points of view, suggesting that low-mass giants can produce Li inside themselves to a super high level via $^{7}$Be transportation during the red giant phase.

We present our deep $Y$-band imaging data of a two square degree field within the F22 region of the VIMOS VLT Deep Survey. The observations were conducted using the WIRCam instrument mounted at the Canada--France--Hawaii Telescope (CFHT). The total on-sky time was 9 hours, distributed uniformly over 18 tiles. The scientific goals of the project are to select faint quasar candidates at redshift $z>2.2$, and constrain the photometric redshifts for quasars and galaxies. In this paper, we present the observation and the image reduction, as well as the photometric redshifts that we derived by combining our $Y$-band data with the CFHTLenS $u^*g'r'i'z'$ optical data and UKIDSS DXS $JHK$ near-infrared data. With $J$-band image as reference total $\sim$80,000 galaxies are detected in the final mosaic down to $Y$-band $5\sigma$ point source limiting depth of 22.86 mag. Compared with the $\sim$3500 spectroscopic redshifts, our photometric redshifts for galaxies with $z<1.5$ and $i'\lesssim24.0$ mag have a small systematic offset of $|\Delta{z}|\lesssim0.2$, 1$\sigma$ scatter $0.03<\sigma_{\Delta z} < 0.06$, and less than 4.0% of catastrophic failures. We also compare to the CFHTLenS photometric redshifts, and find that ours are more reliable at $z\gtrsim0.6$ because of the inclusion of the near-infrared bands. In particular, including the $Y$-band data can improve the accuracy at $z\sim 1.0-2.0$ because the location of the 4000\AA-break is better constrained. The $Y$-band images, the multi-band photometry catalog and the photometric redshifts are released at \url{http://astro.pku.edu.cn/astro/data/DYI.html}.

Bi-site time-series photometric and high-resolution spectroscopic observations were made for the double-mode high-amplitude $\delta$ Scuti star VX\ Hya. The fundamental frequency $f_{0}=4.4763\ \rm{c\ days^{-1}}$, the first overtone $f_{1}=5.7897\ \rm{c\ days^{-1}}$ and 23 harmonics and linear combinations of $f_{0}$ and $f_{1}$ are detected by pulsation analysis. From the spectroscopic data, we get $\rm{[Fe/H] = -0.2\pm0.1\ dex}$. The period change rate of the fundamental mode is obtained by using the Fourier-phase diagram method, which gives the value of $(1/P_{0})(dP_{0}/dt)=(1.81\pm0.09) \times 10^{-7}\ \rm{yr}^{-1}$. With these results from the observations, we perform theoretical explorations with the stellar evolution code MESA, and constrain the models by fitting $f_{0}$, $f_{1}$, and $(1/P_{0})(dP_{0}/dt)$ within $3\sigma$ deviations. The results show that the period change of VX\ Hya could be ascribed to the evolutionary effect. The stellar parameters of VX\ Hya could be derived as: the mass $2.385\pm0.025\ M_{\odot}$, the luminosity $\log(L/L_{\odot})=1.93\pm0.02$ and the age $(4.43\pm0.13)\times 10^8$ years. VX\ Hya is found to locate at the post-main-sequence stage with a helium core and a hydrogen-burning shell on the H${-}$R diagram.

Spectroscopic observations play essential roles in astrophysics. They are crucial for determining physical parameters in the universe, providing information about the chemistry of various astronomical environments. The proper execution of the spectroscopic analysis requires accounting for all the physical effects that are compatible to the signal-to-noise ratio. We find in this paper the influence on spectroscopy from the atomic/ground state alignment owing to anisotropic radiation and modulated by interstellar magnetic field, has significant impact on the study of interstellar gas. In different observational scenarios, we comprehensively demonstrate how atomic alignment influences the spectral analysis and provide the expressions for correcting the effect. The variations are even more pronounced for multiplets and line ratios. We show the variation of the deduced physical parameters caused by the atomic alignment effect, including alpha-to-iron ratio ([X/Fe]) and ionisation fraction. Synthetic observations are performed to illustrate the visibility of such effect with current facilities. A study of PDRs in $\rho$ Ophiuchi cloud is presented to demonstrate how to account for atomic alignment in practice. Our work has shown that due to its potential impact, atomic alignment has to be included in an accurate spectroscopic analysis of the interstellar gas with current observational capability.

Spectroscopic observations play a fundamental role in astrophysics. They are crutial to determine important physical parameters, provide information about the composition of various objects in the universe, as well as depict motions in the universe. However, spectroscopic studies often do not consider the influence of magnetic fields. In this paper, we explore the influence of magnetic fields on the spectroscopic observations using the concept of atomic alignment. Synthetic spectra are generated to show the measurable changes of the spectra due to atomic alignment. The influences of atomic alignment on absorption from DLAs, emission from H\,{\sc ii} Regions, submillimeter fine-structure lines from star forming regions are presented as examples to show this effect in diffuse gas. Furthermore, we demonstrate the influence of atomic alignment on physical parameters derived from atomic line ratios, such as the alpha-to-iron ratio([X/Fe]), interstellar temperature, and ionization rate. We conclude that Ground State Alignment (GSA) should be taken into consideration in the error budget of spectroscopic studies with high signal-to-noise(S/N) ratio.

The copper abundances of 29 metal-poor stars are determined based on the high resolution, high signal-to-noise ratio spectra from the UVES spectragraph at the ESO VLT telescope. Our sample consists of the stars of the Galactic halo, thick- and thin-disk with [Fe/H] ranging from ~ -3.2 to ~ 0.0 dex. The non-local thermodynamic equilibrium (NLTE) effects of Cu I lines are investigated, and line formation calculations are presented for an atomic model of copper including 97 terms and 1089 line transitions. We adopted the recently calculated photo-ionization cross-sections of Cu I, and investigated the hydrogen collision by comparing the theoretical and observed line profiles of our sample stars. The copper abundances are derived for both local thermodynamic equilibrium (LTE) and NLTE based on the spectrum synthesis methods. Our results show that the NLTE effects for Cu I lines are important for metal-poor stars, in particular for very metal-poor stars, and these effects depend on the metallicity. For very metal-poor stars, the NLTE abundance correction reaches as large as ~ +0.5 dex compared to standard LTE calculations. Our results indicate that [Cu/Fe] is under-abundant for metal-poor stars (~ -0.5 dex) when the NLTE effects are included.

The lithium abundances in a few percent of giants exceed the value predicted by the standard stellar evolution models, and the mechanisms of Li enhancement are still under debate. The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) survey has obtained over six million spectra in the past five years, and thus provides a great opportunity to search these rare objects and to more clearly understand the mechanisms of Li enhancement. Based on the high-resolution spectrum we obtained the stellar parameters ($T_\mathrm{eff}$, $\log g$, [Fe/H]), and determined the elemental abundances of Li, C, N, $\alpha$, Fe-peak, r-process, s-process elements, and the projected rotational velocity. For a better understanding of the effect of mixing processes, we also derived the $^{12}\rm{C}$ to $^{13}\rm{C}$ ratio, and constrained the evolutionary status of TYC\,3251-581-1 based on the BaSTI stellar isochrones. The super Li-rich giant TYC\,3251-581-1 has $\rm{A(Li)} = 3.51$, the average abundance of two lithium lines at $\lambda = 6708$ \AA\ and 6104 \AA\ based on the non-local thermodynamic equilibrium (NLTE) analysis. The atmospheric parameters show that our target locates on the luminosity function bump. The low carbon isotopic ratio ($^{12}\rm{C}/^{13}\rm{C} = 9.0 $), a slow rotational velocity $v\sin i = 2.2\ \rm{km\,s}^{-1}$, and no sign of IR excess suggest that additional mixing after first dredge up (FDU) should occur to bring internal synthesized Li to the surface. The low carbon ($[\rm{C}/\rm{Fe}] \sim -0.34$ ) and enhanced nitrogen ($[\rm{N}/\rm{Fe}] \sim 0.33$) are also consistent with the sign of mixing. Given the evolutionary stage of TYC\,3251-581-1 with the relatively low $^{12}\rm{C}/^{13}\rm{C}$, the internal production which replenishes Li in the outer layer is the most likely origin of Li enhancement for this star.

The transport of cosmic rays (CRs) is crucial for the understanding of almost all high-energy phenomena. Both pre-existing large-scale magnetohydrodynamic (MHD) turbulence and locally generated turbulence through plasma instabilities are important for the CR propagation in astrophysical media. The potential role of the resonant instability triggered by CR pressure anisotropy to regulate the parallel spatial diffusion of low-energy CRs ($\lesssim 100$ GeV) in the interstellar and intracluster medium of galaxies (ISM and ICM) has been showed in previous theoretical works. This work aims to study the gyroresonance instability via direct numerical simulations, in order to access quantitatively the wave-particle scattering rates. For this we employ a 1D PIC-MHD code to follow the growth and saturation of the gyroresonance instability. We extract from the simulations the pitch-angle diffusion coefficient $D_{\mu\mu}$ produced by the instability during the linear and saturation phases, and a very good agreement (within a factor of 3) is found with the values predicted by the quasilinear theory (QLT). Our results support the applicability of the QLT for modeling the scattering of low-energy CRs by the gyroresonance instability in the complex interplay between this instability and the large-scale MHD turbulence.

We present the results of a VLT MUSE/FORS2 and Spitzer survey of a unique compact lensing cluster CLIO at z = 0.42, discovered through the GAMA survey using spectroscopic redshifts. Compact and massive clusters such as this are understudied, but provide a unique prospective on dark matter distributions and for finding background lensed high-z galaxies. The CLIO cluster was identified for follow up observations due to its almost unique combination of high mass and dark matter halo concentration, as well as having observed lensing arcs from ground based images. Using dual band optical and infra-red imaging from FORS2 and Spitzer, in combination with MUSE optical spectroscopy we identify 89 cluster members and find background sources out to z = 6.49. We describe the physical state of this cluster, finding a strong correlation between environment and galaxy spectral type. Under the assumption of a NFW profile, we measure the total mass of CLIO to be M$_{200} = (4.49 \pm 0.25) \times 10^{14}$ M$_\odot$. We build and present an initial strong-lensing model for this cluster, and measure a relatively low intracluster light (ICL) fraction of 7.21 $\pm$ 1.53% through galaxy profile fitting. Due to its strong potential for lensing background galaxies and its low ICL, the CLIO cluster will be a target for our 110 hour JWST 'Webb Medium-Deep Field' (WMDF) GTO program.

Magnetic fields play important roles in many astrophysical processes. However, there is no universal diagnostic for the magnetic fields in the interstellar medium (ISM) and each magnetic tracer has its limitation. Any new detection method is thus valuable. Theoretical studies have shown that submillimetre fine-structure lines are polarised due to atomic alignment by Ultraviolet (UV) photon-excitation, which opens up a new avenue to probe interstellar magnetic fields. We will, for the first time, perform synthetic observations on the simulated three-dimensional ISM to demonstrate the measurability of the polarisation of submillimetre atomic lines. The maximum polarisation for different absorption and emission lines expected from various sources, including Star-Forming Regions (SFRs) are provided. Our results demonstrate that the polarisation of submillimetre atomic lines is a powerful magnetic tracer and add great value to the observational studies of the submilimetre astronomy.

Brightest Cluster Galaxies (BCGs) are believed to have assembled most of their stars early in time and, therefore, should be passively evolving at low redshifts and appear "red-and-dead." However, there have been reports that a minority of low-redshift BCGs still have ongoing star formation rates (SFR) of a few to even $\sim$100 $M_\odot/yr$. Such BCGs are found in "cool-core" ("CC") clusters, and their star formation is thought to be fueled by "cooling flow." To further investigate the implications of low-redshift, star-forming BCGs, we perform a systematic search using the 22$\mu$m data ("W4" band) from the Wide-field Infrared Survey Explorer (WISE) on the GMBCG catalog, which contains 55,424 BCGs at $0.1\lesssim z\lesssim 0.55$ identified in the Sloan Digital Sky Survey (SDSS). Our sample consists of 389 BCGs that are bright in W4 ("W4BCGs"), most being brighter than 5 mJy. While some ($\lesssim 20\%$) might host AGN, most W4BCGs should owe their strong mid-IR emissions to dust-enshrouded star formation. Their median total IR luminosity ($L_{IR}$) is $5\times10^{11} L_{\odot}$ (SFR $\sim$50 $M_{\odot}/yr)$, and 27\% of the whole sample has $L_{IR}>10^{12} L_{\odot}$ (SFR $>$100 $M_{\odot}/yr$). Using ten W4BCGs that have Chandra X-ray data, we show that seven of them are possibly in CC clusters. However, in most cases (five out of seven) the mass deposition rate cannot account for the observed SFR. This casts doubt to the idea that cooling flows are the cause of the star formation in non-quiescent BCGs.

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

We examine the fraction of massive ($M_{*}>10^{10} M_{\odot}$), compact star-forming galaxies (cSFGs) that host an active galactic nucleus (AGN) at $z\sim2$. These cSFGs are likely the direct progenitors of the compact quiescent galaxies observed at this epoch, which are the first population of passive galaxies to appear in large numbers in the early Universe. We identify cSFGs that host an AGN using a combination of Hubble WFC3 imaging and Chandra X-ray observations in four fields: the Chandra Deep Fields, the Extended Groth Strip, and the UKIDSS Ultra Deep Survey field. We find that $39.2^{+3.9}_{-3.6}$\% (65/166) of cSFGs at $1.4<z<3.0$ host an X-ray detected AGN. This fraction is 3.2 times higher than the incidence of AGN in extended star-forming galaxies with similar masses at these redshifts. This difference is significant at the $6.2\sigma$ level. Our results are consistent with models in which cSFGs are formed through a dissipative contraction that triggers a compact starburst and concurrent growth of the central black hole. We also discuss our findings in the context of cosmological galaxy evolution simulations that require feedback energy to rapidly quench cSFGs. We show that the AGN fraction peaks precisely where energy injection is needed to reproduce the decline in the number density of cSFGs with redshift. Our results suggest that the first abundant population of massive, quenched galaxies emerged directly following a phase of elevated supermassive black hole growth and further hints at a possible connection between AGN and the rapid quenching of star formation in these galaxies.

The fraction of Lyman-$\alpha$ emitters among the galaxy population has been found to increase from $z\sim0$ to $z\sim6$ and drop dramatically at $z>6$. This drop has been interpreted as an effect of an increasingly neutral intergalactic medium with increasing redshift, while a LyC escape fraction evolving with redshift. We report the result of a large VLT/FORS2 program aiming to confirm spectroscopically a large galaxy sample at $z\geq6$ that has been selected in several independent fields through the Lyman Break technique. Combining those data with archival data, we create a large and homogeneous sample of $z\sim6$ galaxies ($N=127$), complete in terms of Ly$\alpha$ detection at $>95\%$ for EW(Ly$\alpha)\geq25\AA$. We use this sample to derive a new measurement of the LAE fraction at $z\sim6$ and derive the physical properties of these galaxies through spectral energy distribution fitting. We find a median LAE fraction at $z\sim6$ lower than in previous studies. The main difference between LAEs and non-LAEs is that the latter are significantly dustier. Using predictions of our SED fitting code accounting for nebular emission, we find an effective Ly$\alpha$ escape fraction $f^{eff}_{esc}(Ly\alpha)=0.23^{+0.36}_{-0.17}$ remarkably consistent with the value derived by comparing UV luminosity function with Ly$\alpha$ luminosity function. We conclude that the drop in the LAE fraction from $z\sim6$ to $z>6$ is less dramatic than previously found and the effect of an increasing IGM neutral fraction is possibly observed at $5<z<6$. Based on our derived $f^{eff}_{esc}(Ly\alpha)$, we find that the IGM has a relatively small impact on Ly$\alpha$ photon visibility at $z\sim6$, with a lower limit for the IGM transmission to \lya\ photons, $T_{IGM}\gtrsim0.20$, likely due to the presence of outflows. [abdridged]

During the Space Telescope and Optical Reverberation Mapping Project (STORM) observations of NGC 5548, the continuum and emission-line variability became de-correlated during the second half of the 6-month long observing campaign. Here we present Swift and Chandra X-ray spectra of NGC 5548 obtained as a part of the campaign. The Swift spectra show that excess flux (relative to a power-law continuum) in the soft X-ray band appears before the start of the anomalous emission-line behavior, peaks during the period of the anomaly, and then declines. This is a model-independent result suggesting that the soft excess is related to the anomaly. We divide the Swift data into on- and off-anomaly spectra to characterize the soft excess via spectral fitting. The cause of the spectral differences is likely due to a change in the intrinsic spectrum rather than being due to variable obscuration or partial covering. The Chandra spectra have lower signal-to-noise ratios, but are consistent with Swift data. Our preferred model of the soft excess is emission from an optically thick, warm Comptonizing corona, the effective optical depth of which increases during the anomaly. This model simultaneously explains all the three observations: the UV emission line flux decrease, the soft-excess increase, and the emission line anomaly.

We analyze a sample of $z$-dropout galaxies in the CANDELS GOODS South and UDS fields that have been targeted by a dedicated spectroscopic campaign aimed at detecting their Ly$\alpha$ line. Deep IRAC observations at 3.6 and 4.5 $\mu$m are used to determine the strength of optical emission lines affecting these bands at z$\sim$6.5-6.9 in order to i) investigate possible physical differences between Ly$\alpha$ emitting and non-emitting sources; ii) constrain the escape fraction of ionizing photons; iii) provide an estimate of the specific star-formation rate at high redshifts. We find evidence of strong [OIII]+H$\beta$ emission in the average (stacked) SEDs of galaxies both with and without Ly$\alpha$ emission. The blue IRAC [3.6]-[4.5] color of the stack with detected Ly$\alpha$ line can be converted into a rest-frame equivalent width EW([OIII]+H$\beta$)=1500$^{+530}_{-440}\AA$ assuming a flat intrinsic stellar continuum. This strong optical line emission enables a first estimate of f$_{esc}\lesssim$20% on the escape fraction of ionizing photons from Ly$\alpha$ detected objects. The objects with no Ly$\alpha$ line show less extreme EW([OIII]+H$\beta$)=520$^{+170}_{-150}\AA$ suggesting different physical conditions of the HII regions with respect to Ly$\alpha$-emitting ones, or a larger f$_{esc}$. The latter case is consistent with a combined evolution of f$_{esc}$ and the neutral hydrogen fraction as an explanation of the lack of bright Ly$\alpha$ emission at z$>$6. A lower limit on the specific star formation rate, SSFR$>$9.1$Gyr^{-1}$ for $M_{star}=2 \times 10^9 M_{\odot}$ galaxies at these redshifts can be derived from the spectroscopically confirmed sample.

We present a 0.4-8$\mu$m multi-wavelength photometric catalog in the Extended Groth Strip (EGS) field. This catalog is built on the Hubble Space Telescope (HST) WFC3 and ACS data from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS), and it incorporates the existing HST data from the All-wavelength Extended Groth strip International Survey (AEGIS) and the 3D-HST program. The catalog is based on detections in the F160W band reaching a depth of F160W=26.62 AB (90% completeness, point-sources). It includes the photometry for 41457 objects over an area of $\approx 206$ arcmin$^2$ in the following bands: HST ACS F606W and F814W; HST WFC3 F125W, F140W and F160W; CFHT/Megacam $u^*$, $g'$, $r'$, $i'$ and $z'$; CFHT/WIRCAM $J$, $H$ and $K_\mathrm{S}$; Mayall/NEWFIRM $J1$, $J2$, $J3$, $H1$, $H2$, $K$; Spitzer IRAC $3.6\mu$m, $4.5\mu$m, $5.8\mu$m and $8.0\mu$m. We are also releasing value-added catalogs that provide robust photometric redshifts and stellar mass measurements. The catalogs are publicly available through the CANDELS repository.

The reionization of the Universe is one of the most important topics of present day astrophysical research. The most plausible candidates for the reionization process are star-forming galaxies, which according to the predictions of the majority of the theoretical and semi-analytical models should dominate the HI ionizing background at z~3. We aim at measuring the Lyman continuum escape fraction, which is one of the key parameters to compute the contribution of star-forming galaxies to the UV background. We have used ultra-deep U-band imaging (U=30.2mag at 1sigma) by LBC/LBT in the CANDELS/GOODS-North field, as well as deep imaging in COSMOS and EGS fields, in order to estimate the Lyman continuum escape fraction of 69 star-forming galaxies with secure spectroscopic redshifts at 3.27<z<3.40 to faint magnitude limits (L=0.2L*, or equivalently M1500~-19). We have measured through stacks a stringent upper limit (<1.7% at 1sigma) for the relative escape fraction of HI ionizing photons from bright galaxies (L>L*), while for the faint population (L=0.2L*) the limit to the escape fraction is ~10%. We have computed the contribution of star-forming galaxies to the observed UV background at z~3 and we have found that it is not enough to keep the Universe ionized at these redshifts, unless their escape fraction increases significantly (>10%) at low luminosities (M1500>-19). We compare our results on the Lyman continuum escape fraction of high-z galaxies with recent estimates in the literature and discuss future prospects to shed light on the end of the Dark Ages. In the future, strong gravitational lensing will be fundamental to measure the Lyman continuum escape fraction down to faint magnitudes (M1500~-16) which are inaccessible with the present instrumentation on blank fields.

We present the results of an optical spectroscopic monitoring program targeting NGC 5548 as part of a larger multi-wavelength reverberation mapping campaign. The campaign spanned six months and achieved an almost daily cadence with observations from five ground-based telescopes. The H$\beta$ and He II $\lambda$4686 broad emission-line light curves lag that of the 5100 $\AA$ optical continuum by $4.17^{+0.36}_{-0.36}$ days and $0.79^{+0.35}_{-0.34}$ days, respectively. The H$\beta$ lag relative to the 1158 $\AA$ ultraviolet continuum light curve measured by the Hubble Space Telescope is roughly $\sim$50% longer than that measured against the optical continuum, and the lag difference is consistent with the observed lag between the optical and ultraviolet continua. This suggests that the characteristic radius of the broad-line region is $\sim$50% larger than the value inferred from optical data alone. We also measured velocity-resolved emission-line lags for H$\beta$ and found a complex velocity-lag structure with shorter lags in the line wings, indicative of a broad-line region dominated by Keplerian motion. The responses of both the H$\beta$ and He II $\lambda$4686 emission lines to the driving continuum changed significantly halfway through the campaign, a phenomenon also observed for C IV, Ly $\alpha$, He II(+O III]), and Si IV(+O IV]) during the same monitoring period. Finally, given the optical luminosity of NGC 5548 during our campaign, the measured H$\beta$ lag is a factor of five shorter than the expected value implied by the $R_\mathrm{BLR} - L_\mathrm{AGN}$ relation based on the past behavior of NGC 5548.

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 a multi-wavelength photometric catalog in the COSMOS field as part of the observations by the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). The catalog is based on Hubble Space Telescope Wide Field Camera 3 (HST/WFC3) and Advanced Camera for Surveys (ACS) observations of the COSMOS field (centered at RA: $10^h00^m28^s$, Dec:$+02^{\circ}12^{\prime}21^{\prime\prime}$). The final catalog has 38671 sources with photometric data in forty two bands from UV to the infrared ($\rm \sim 0.3-8\,\mu m$). This includes broad-band photometry from the HST, CFHT, Subaru, VISTA and Spitzer Space Telescope in the visible, near infrared and infrared bands along with intermediate and narrow-band photometry from Subaru and medium band data from Mayall NEWFIRM. Source detection was conducted in the WFC3 F160W band (at $\rm 1.6\,\mu m$) and photometry is generated using the Template FITting algorithm. We further present a catalog of the physical properties of sources as identified in the HST F160W band and measured from the multi-band photometry by fitting the observed spectral energy distributions of sources against templates.

Diffusion of cosmic rays (CRs) is the key process of understanding their propagation and acceleration. We employ the description of spatial separation of magnetic field lines in MHD turbulence in Lazarian & Vishniac (1999) to quantify the divergence of magnetic field on scales less than the injection scale of turbulence and show this divergence induces superdiffusion of CR in the direction perpendicular to the mean magnetic field. The perpendicular displacement squared increases, not as distance $x$ along magnetic field, which is the case for a regular diffusion, but as the $x^{3}$ for freely streaming CRs. The dependence changes to $x^{3/2}$ for the CRs propagating diffusively along magnetic field. In the latter case we show that it is important to distinguish the perpendicular displacement in respect to the mean field and to the local magnetic field. We consider how superdiffusion changes the acceleration of CRs in shocks and show how it decreases efficiency of the CRs acceleration in perpendicular shocks. We also demonstrate that in the case when small-scale magnetic field is being generated in the pre-shock region, an efficient acceleration can take place for the CRs streaming without collisions along magnetic loops.

We conduct a multiwavelength continuum variability study of the Seyfert 1 galaxy NGC 5548 to investigate the temperature structure of its accretion disk. The 19 overlapping continuum light curves (1158 to 9157 angstroms) combine simultaneous HST , Swift , and ground-based observations over a 180 day period from 2014 January to July. Light-curve variability is interpreted as the reverberation response of the accretion disk to irradiation by a central time-varying point source. Our model yields the disk inclination, i, temperature T1 at 1 light day from the black hole, and a temperature-radius slope, alpha. We also infer the driving light curve and find that it correlates poorly with both the hard and soft X-ray light curves, suggesting that the X-rays alone may not drive the ultraviolet and optical variability over the observing period. We also decompose the light curves into bright, faint, and mean accretion-disk spectra. These spectra lie below that expected for a standard blackbody accretion disk accreting at L/LEdd = 0.1

[Methods]. We obtained high-precision radial velocities with HARPS on the ESO 3.6 m telescope and determined precise stellar elemental abundances (~0.01 dex) using MIKE spectra on the Magellan 6.5m telescope. [Results]. Our data indicate the presence of a planet with a minimum mass of 26 Earth masses around the solar twin HIP 68468. The planet is a super-Neptune, but unlike the distant Neptune in our solar system (30 AU), HIP 68468c is close-in, with a semi-major axis of 0.66 AU, similar to that of Venus. The data also suggest the presence of a super-Earth with a minimum mass of 2.9 Earth masses at 0.03 AU; if the planet is confirmed, it will be the fifth least massive radial velocity planet discovery to date and the first super-Earth around a solar twin. Both isochrones (5.9 Gyr) and the abundance ratio [Y/Mg] (6.4 Gyr) indicate an age of about 6 billion years. The star is enhanced in refractory elements when compared to the Sun, and the refractory enrichment is even stronger after corrections for Galactic chemical evolution. We determined a NLTE Li abundance of 1.52 dex, which is four times higher than what would be expected for the age of HIP 68468. The older age is also supported by the low log(R'HK) (-5.05) and low jitter. Engulfment of a rocky planet of 6 Earth masses can explain the enhancement in both lithium and the refractory elements. [Conclusions]. The super-Neptune planet candidate is too massive for in situ formation, and therefore its current location is most likely the result of planet migration that could also have driven other planets towards its host star, enhancing thus the abundance of lithium and refractory elements in HIP 68468. The intriguing evidence of planet accretion warrants further observations to verify the existence of the planets that are indicated by our data and to better constrain the nature of the planetary system around this unique star.

We performed numerical calculations to test the suggestion by Desiati & Lazarian (2013) that the anisotropies of TeV cosmic rays may arise from their interactions with the heliosphere. For this purpose, we used a magnetic field model of the heliosphere and performed direct numerical calculations of particle trajectories. Unlike earlier papers testing the idea, we did not employ time-reversible techniques that are based on Liouville's theorem. We showed numerically that for scattering by the heliosphere the conditions of Liouville's theorem are not satisfied and the adiabatic approximation and time-reversibility of the particle trajectories are not valid. Our results indicate sensitivity to the magnetic structure of the heliospheric magnetic field, and we expect that this will be useful for probing this structure in future research.

For the first time, we present an extensive study of stars with individual non-local thermodynamic equilibrium (NLTE) abundances for 17 chemical elements from Li to Eu in a sample of stars uniformly distributed over the $-2.62 \le$ [Fe/H] $\le +0.24$ metallicity range that is suitable for the Galactic chemical evolution research. The star sample has been kinematically selected to trace the Galactic thin and thick disks and halo. We find new and improve earlier results as follows. (i) The element-to-iron ratios for Mg, Si, Ca, and Ti form a MP plateau at a similar height of 0.3~dex, and the knee occurs at common [Fe/H] $\simeq -0.8$. The knee at the same metallicity is observed for [O/Fe], and the MP plateau is formed at [O/Fe] = 0.61. (ii) The upward trend of [C/O] with decreasing metallicity exists at [Fe/H] $< -1.2$, supporting the earlier finding of Akerman et al. (iii) An underabundance of Na relative to Mg in the [Fe/H] $< -1$ stars is nearly constant, with the mean [Na/Mg] $\simeq -0.5$. (iv) The K/Sc, Ca/Sc, and Ti/Sc ratios form well-defined trends, suggesting a common site of the K-Ti production. (v) Sr follows the Fe abundance down to [Fe/H] $\simeq -2.5$, while Zr is enhanced in MP stars. (vi) The comparisons of our results with some widely used Galactic evolution models are given. The use of the NLTE element abundances raises credit to the interpretation of the data in the context of the chemical evolution of the Galaxy.

We present new results on [CII]158$\mu$ m emission from four galaxies in the reionization epoch. These galaxies were previously confirmed to be at redshifts between 6.6 and 7.15 from the presence of the Ly$\alpha$ emission line in their spectra. The Ly$\alpha$ emission line is redshifted by 100-200 km/s compared to the systemic redshift given by the [CII] line. These velocity offsets are smaller than what is observed in z~3 Lyman break galaxies with similar UV luminosities and emission line properties. Smaller velocity shifts reduce the visibility of Ly$\alpha$ and hence somewhat alleviate the need for a very neutral IGM at z~7 to explain the drop in the fraction of Ly$\alpha$ emitters observed at this epoch. The galaxies show [CII] emission with L[CII]=0.6-1.6 x10$^8 L_\odot$: these luminosities place them consistently below the SFR-L[CII] relation observed for low redshift star forming and metal poor galaxies and also below z =5.5 Lyman break galaxies with similar star formation rates. We argue that previous undetections of [CII] in z~7 galaxies with similar or smaller star formation rates are due to selection effects: previous targets were mostly strong Ly$\alpha$ emitters and therefore probably metal poor systems, while our galaxies are more representative of the general high redshift star forming population .

Turbulence in the weakly collisional intracluster medium of galaxies (ICM) is able to generate strong thermal velocity anisotropies in the ions (with respect to the local magnetic field direction), if the magnetic moment of the particles is conserved in the absence of Coulomb collisions. In this scenario, the anisotropic pressure magnetohydrodynamic (AMHD) turbulence shows a very different statistical behaviour from the standard MHD one and is unable to amplify seed magnetic fields, in disagreement with previous cosmological MHD simulations which are successful to explain the observed magnetic fields in the ICM. On the other hand, temperature anisotropies can also drive plasma instabilities which can relax the anisotropy. This work aims to compare the relaxation rate with the growth rate of the anisotropies driven by the turbulence. We employ quasilinear theory to estimate the ions scattering rate due to the parallel firehose, mirror, and ion-cyclotron instabilities, for a set of plasma parameters resulting from AMHD simulations of the turbulent ICM. We show that the ICM turbulence can sustain only anisotropy levels very close to the instabilities thresholds. We argue that the AMHD model which bounds the anisotropies at the marginal stability levels can describe the Alfvenic turbulence cascade in the ICM.

During an intensive Hubble Space Telescope (HST) Cosmic Origins Spectrograph (COS) UV monitoring campaign of the Seyfert~1 galaxy NGC 5548 performed from 2014 February to July, the normally highly correlated far-UV continuum and broad emission-line variations decorrelated for ~60 to 70 days, starting ~75 days after the first HST/COS observation. Following this anomalous state, the flux and variability of the broad emission lines returned to a more normal state. This transient behavior, characterised by significant deficits in flux and equivalent width of the strong broad UV emission lines, is the first of its kind to be unambiguously identified in an active galactic nucleus reverberation mapping campaign. The largest corresponding emission-line flux deficits occurred for the high-ionization collisionally excited lines, C IV and Si IV(+O IV]), and also He II(+O III]), while the anomaly in Ly-alpha was substantially smaller. This pattern of behavior indicates a depletion in the flux of photons with E_{\rm ph} > 54 eV, relative to those near 13.6 eV. We suggest two plausible mechanisms for the observed behavior: (i) temporary obscuration of the ionizing continuum incident upon BLR clouds by a moving veil of material lying between the inner accretion disk and inner BLR, perhaps resulting from an episodic ejection of material from the disk, or (ii) a temporary change in the intrinsic ionizing continuum spectral energy distribution resulting in a deficit of ionizing photons with energies > 54 eV, possibly due to a transient restructuring of the Comptonizing atmosphere above the disk. Current evidence appears to favor the latter explanation.

We present ground-based optical photometric monitoring data for NGC 5548, part of an extended multi-wavelength reverberation mapping campaign. The light curves have nearly daily cadence from 2014 January to July in nine filters (\emph{BVRI} and \emph{ugriz}). Combined with ultraviolet data from the \emph{Hubble Space Telescope} and \emph{Swift}, we confirm significant time delays between the continuum bands as a function of wavelength, extending the wavelength coverage from 1158\,\AA\ to the $z$ band ($\sim\!9160$\,\AA). We find that the lags at wavelengths longer than the {\it V} band are equal to or greater than the lags of high-ionization-state emission lines (such as He\,{\sc ii}\,$\lambda 1640$ and $\lambda 4686$), suggesting that the continuum-emitting source is of a physical size comparable to the inner broad-line region (BLR). The trend of lag with wavelength is broadly consistent with the prediction for continuum reprocessing by an accretion disk with $\tau \propto \lambda^{4/3}$. However, the lags also imply a disk radius that is 3 times larger than the prediction from standard thin-disk theory, assuming that the bolometric luminosity is 10\% of the Eddington luminosity ($L = 0.1L_{\rm Edd}$). Using optical spectra from the Large Binocular Telescope, we estimate the bias of the interband continuum lags due to BLR emission observed in the filters. We find that the bias for filters with high levels of BLR contamination ($\sim\! 20\%$) can be important for the shortest continuum lags, and likely has a significant impact on the {\it u} and {\it U} bands owing to Balmer continuum emission.

Using the far-infrared data obtained by the Herschel Space Observatory, we study the relation between the infrared luminosity (L_IR) and the dust temperature (T) of dusty starbursting galaxies at high redshifts (high-z). We focus on the total infrared luminosity from the cold-dust component (L_IR^(cd)), whose emission can be described by a modified black body (MBB) of a single temperature (T_mbb). An object on the (L_IR^(cd), T_mbb) plane can be explained by the equivalent of the Stefan-Boltzmann law for a MBB with an effective radius of R_eff. We show that R_eff is a good measure of the combined size of the dusty starbursting regions (DSBRs) of the host galaxy. In at least one case where the individual DSBRs are well resolved through strong gravitational lensing, R_eff is consistent with the direct size measurement. We show that the observed L_IR-T relation is simply due to the limited R_eff (<~ 2 kpc). The small R_eff values also agree with the compact sizes of the DSBRs seen in the local universe. However, previous interferometric observations to resolve high-z dusty starbursting galaxies often quote much larger sizes. This inconsistency can be reconciled by the blending effect when considering that the current interferometry might still not be of sufficient resolution. From R_eff we infer the lower limits to the volume densities of the star formation rate ("minSFR3D") in the DSBRs, and find that the $L_{IR}$-$T$ relation outlines a boundary on the (L_IR^(cd), T) plane, below which is the "zone of avoidance" in terms of minSFR3D.

Theoretical studies have revealed that dust grains are usually moving fast through the turbulent interstellar gas, which could have significant effects upon interstellar chemistry by modifying grain accretion. This effect is investigated in this work on the basis of numerical gas-grain chemical modeling. Major features of the grain motion effect in the typical environment of dark clouds (DC) can be summarised as follows: 1) decrease of gas-phase (both neutral and ionic) abundances and increase of surface abundances by up to 2-3 orders of magnitude; 2) shifts of the existing chemical jumps to earlier evolution ages for gas-phase species and to later ages for surface species by factors of about ten; 3) a few exceptional cases in which some species turn out to be insensitive to this effect and some other species can show opposite behaviors too. These effects usually begin to emerge from a typical DC model age of about 10^5 yr. The grain motion in a typical cold neutral medium (CNM) can help overcome the Coulomb repulsive barrier to enable effective accretion of cations onto positively charged grains. As a result, the grain motion greatly enhances the abundances of some gas-phase and surface species by factors up to 2-6 or more orders of magnitude in the CNM model. The grain motion effect in a typical molecular cloud (MC) is intermediate between that of the DC and CNM models, but with weaker strength. The grain motion is found to be important to consider in chemical simulations of typical interstellar medium.

Two distinct halo populations were found in the solar neighborhood by a series of works. They can be clearly separated by [alpha\Fe] and several other elemental abundance ratios including [Cu/Fe]. Very recently, a non-local thermodynamic equilibrium (non-LTE) study revealed that relatively large departures exist between LTE and non-LTE results in copper abundance analysis. We aim to derive the copper abundances for the stars from the sample of Nissen et al (2010) with both LTE and non-LTE calculations. Based on our results, we study the non-LTE effects of copper and investigate whether the high-alpha population can still be distinguished from the low-alpha population in the non-LTE [Cu/Fe] results. Our differential abundance ratios are derived from the high-resolution spectra collected from VLT/UVES and NOT/FIES spectrographs. Applying the MAFAGS opacity sampling atmospheric models and spectrum synthesis method, we derive the non-LTE copper abundances based on the new atomic model with current atomic data obtained from both laboratory and theoretical calculations. The copper abundances determined from non-LTE calculations are increased by 0.01 to 0.2 dex depending on the stellar parameters compared with the LTE results. The non-LTE [Cu/Fe] trend is much flatter than the LTE one in the metallicity range -1.6<[Fe/H]<-0.8. Taking non-LTE effects into consideration, the high- and low-alpha stars still show distinguishable copper abundances, which appear even more clear in a diagram of non-LTE [Cu/Fe] versus [Fe/H]. The non-LTE effects are strong for copper, especially in metal-poor stars. Our results confirmed that there are two distinct halo populations in the solar neighborhood. The dichotomy in copper abundance is a peculiar feature of each population, suggesting that they formed in different environments and evolved obeying diverse scenarios.

We study the evolution of the core (r<1 kpc) and effective (r<r_e) stellar-mass surface densities, in star-forming and quiescent galaxies. Since z=3, both populations occupy distinct, linear relations in log(Sigma_e) and log(Sigma_1) vs. log(M). These structural relations exhibit slopes and scatter that remain almost constant with time while their normalizations decline. For SFGs, the normalization declines by less than a factor of 2 from z=3, in both Sigma_e and Sigma_1. Such mild declines suggest that SFGs build dense cores by growing along these relations. We define this evolution as the structural main sequence (Sigma-MS). Quiescent galaxies follow different relations (Sigma^Q_e, Sigma^Q_1) off the Sigma-MS by having higher densities than SFGs of the same mass and redshift. The normalization of Sigma^Q_e declines by a factor of 10 since z=3, but only a factor of 2 in Sigma^Q_1. Thus, the common denominator for quiescent galaxies at all redshifts is the presence of a dense stellar core, and the formation of such cores in SFGs is the main requirement for quenching. Expressed in 2D as deviations off the SFR-MS and off Sigma^Q_1 at each redshift, the distribution of massive galaxies forms a universal, L-shaped sequence that relates two fundamental physical processes: compaction and quenching. Compaction is a process of substantial core-growth in SFGs relative to that in the Sigma-MS. This process increases the core-to-total mass and Sersic index, thereby, making compact SFGs. Quenching occurs once compact SFGs reach a maximum central density above Sigma^Q_1 > 9.5 M_sun/kpc^2. This threshold provides the most effective selection criterion to identify the star-forming progenitors of quiescent galaxies at all redshifts.

Using the public data from the Herschel wide field surveys, we study the far-infrared properties of optical-selected quasars from the Sloan Digital Sky Survey. Within the common area of $\sim 172~deg^2$, we have identified the far-infrared counterparts for 354 quasars, among which 134 are highly secure detections in the Herschel $250~\mu m$ band (signal-to-noise ratios $\geq5$). This sample is the largest far-infrared quasar sample of its kind, and spans a wide redshift range of $0.14{\leq}z\leq 4.7$. Their far-infrared spectral energy distributions, which are due to the cold dust components within the host galaxies, are consistent with being heated by active star formation. In most cases ($\gtrsim80$\%), their total infrared luminosities as inferred from only their far-infrared emissions ($L_{IR}^{(cd)}$) already exceed $10^{12}~L_{\odot}$, and thus these objects qualify as ultra-luminous infrared galaxies. There is no correlation between $L_{IR}^{(cd)}$ and the absolute magnitudes, the black hole masses or the X-ray luminosities of the quasars, which further support that their far-infrared emissions are not due to their active galactic nuclei. A large fraction of these objects ($\gtrsim50\text{--}60\%$) have star formation rates $\gtrsim 300~M_{\odot}yr^{-1}$. Such extreme starbursts among optical quasars, however, is only a few per cent. This fraction varies with redshift, and peaks at around $z\approx2$. Among the entire sample, 136 objects have secure estimates of their cold-dust temperatures ($T$), and we find that there is a dramatic increasing trend of $T$ with increasing $L_{IR}^{(cd)}$. We interpret this trend as the envelope of the general distribution of infrared galaxies on the ($T$, $L_{IR}^{(cd)}$) plane.

We study the damping from neutral-ion collisions of both incompressible and compressible magnetohydrodynamic (MHD) turbulence in partially ionized medium. We start from the linear analysis of MHD waves applying both single-fluid and two-fluid treatments. The damping rates derived from the linear analysis are then used in determining the damping scales of MHD turbulence. The physical connection between the damping scale of MHD turbulence and cutoff boundary of linear MHD waves is investigated. Our analytical results are shown to be applicable in a variety of partially ionized interstellar medium (ISM) phases and solar chromosphere. As a significant astrophysical utility, we introduce damping effects to propagation of cosmic rays in partially ionized ISM. The important role of turbulence damping in both transit-time damping and gyroresonance is identified.

Tracing magnetic fields is crucial as magnetic fields play an important role in many astrophysical processes. Earlier studies have demonstrated that Ground State Alignment (GSA) is a unique way to detect weak magnetic fields (1G> B> 1exp(-15)G) in diffuse media, they consider the situation when the pumping source is a point source, which applies when the star is very far away from the diffuse media. In this paper, we explore the GSA in the presence of extended radiation fields. For the radiation fields with a clear geometric structure, we consider the alignment in circumstellar medium, binary systems, disc, and Local Interstellar Medium (LISM). For the radiation fields with unidentified pumping sources, we apply the method of multipole expansion and discuss the GSA induced by each component. We demonstrate that for general radiation fields, it is adequate to consider the contribution from dipole and quadrupole radiation components. We find that in general polarization of absorption arizing from GSA coincides with the projection of magnetic field in the 2D sky with 90 degree degeneracy. We conclude that the GSA is a unique tool to detect the direction of weak magnetic field, and it can be applied to diffuse media in any radiation field.

We describe the first results from a six-month long reverberation-mapping experiment in the ultraviolet based on 170 observations of the Seyfert 1 galaxy NGC 5548 with the Cosmic Origins Spectrograph on the Hubble Space Telescope. Significant correlated variability is found in the continuum and broad emission lines, with amplitudes ranging from ~30% to a factor of two in the emission lines and a factor of three in the continuum. The variations of all the strong emission lines lag behind those of the continuum, with He II 1640 lagging behind the continuum by ~2.5 days and Lyman alpha 1215, C IV 1550, and Si IV 1400 lagging by ~5-6 days. The relationship between the continuum and emission lines is complex. In particular, during the second half of the campaign, all emission-line lags increased by a factor of 1.3-2 and differences appear in the detailed structure of the continuum and emission-line light curves. Velocity-resolved cross-correlation analysis shows coherent structure in lag versus line-of-sight velocity for the emission lines; the high-velocity wings of C IV respond to continuum variations more rapidly than the line core, probably indicating higher velocity BLR clouds at smaller distances from the central engine. The velocity-dependent response of Lyman alpha, however, is more complex and will require further analysis.

Recent intensive Swift monitoring of the Seyfert 1 galaxy NGC 5548 yielded 282 usable epochs over 125 days across six UV/optical bands and the X-rays. This is the densest extended AGN UV/optical continuum sampling ever obtained, with a mean sampling rate <0.5 day. Approximately daily HST UV sampling was also obtained. The UV/optical light curves show strong correlations (r_max = 0.57 - 0.90) and the clearest measurement to date of interband lags. These lags are well-fit by a \tau propto \lambda^4/3 wavelength dependence, with a normalization that indicates an unexpectedly large disk radius of 0.35 +/- 0.05 lt-day at 1367 A, assuming a simple face-on model. The U-band shows a marginally larger lag than expected from the fit and surrounding bands, which could be due to Balmer continuum emission from the broad-line region as suggested by Korista and Goad. The UV/X-ray correlation is weaker (r_max < 0.45) and less consistent over time. This indicates that while Swift is beginning to measure UV/optical lags in general agreement with accretion disk theory (although the derived size is larger than predicted), the relationship with X-ray variability is less well understood. Combining this accretion disk size estimate with those from quasar microlensing studies suggests that AGN disk sizes scale approximately linearly with central black hole mass over a wide range of masses.

We address the problem of the difference of line widths of neutrals and ions observed from molecular clouds and explore whether this difference can arise from the effects of magnetohydrodynamic (MHD) turbulence acting on partially ionized gas. Among the three fundamental modes of MHD turbulence, we find fast modes do not contribute to linewidth differences, whereas slow modes can have an effect on different line widths for certain parameters. We focus on Alfv\'{e}nic component because they contain most of the turbulent energy, and consider the damping of this component taking into account both neutral-ion collisions and neutral viscosity. We consider different regimes of turbulence corresponding to different media magnetizations and turbulent drivings. In the case of super-Alfv\'{e}nic turbulence, when the damping scale of Alfv\'{e}nic turbulence is below $l_A$, where $l_A$ is the injection scale of anisotropic GS95-type turbulence, the linewidth difference does not depend on the magnetic field strength. While for other turbulent regimes, the dependence is present. For instance, the difference between the squares of the neutral and ion velocity dispersions in strong sub-Alfv\'{e}nic turbulence allows evaluation of magnetic field. We discuss earlier findings on the neutral-ion linewidth differences in the literature and compare the expressions for magnetic field we obtain with those published earlier.

We investigated the copper abundances for $64$ late-type stars in the Galactic disk and halo with effective temperatures from $5400$ K to $6700$ K and [Fe/H] from $-1.88$ to $-0.17$. For the first time, the copper abundances are derived using both local thermodynamic equilibrium (LTE) and non-local thermodynamic equilibrium (non-LTE) calculations. High resolution ($R > 40,000$), high signal-to-noise ratio ($S/N > 100$) spectra from the FOCES spectrograph are used. The atmospheric models are calculated based on the MAFAGS opacity sampling code. All the abundances are derived using the spectrum synthesis methods. Our results indicate that the non-LTE effects of copper are important for metal-poor stars, showing a departure of $\sim 0.17$ dex at the metallicity $\sim -1.5$. We also find that the copper abundances derived from non-LTE calculations are enhanced compared with those from LTE. The enhancements show clear dependence on the metallicity, which gradually increase with decreasing [Fe/H] for our program stars, leading to a flatter distribution of [Cu/Fe] with [Fe/H] than previous work. There is a hint that the thick- and thin-disk stars have different behaviors in [Cu/Fe], and a bending for disk stars may exist.

We present the results of a new, ultra-deep, near-infrared imaging survey executed with the Hawk-I imager at the ESO VLT, of which we make all the data public. This survey, named HUGS (Hawk-I UDS and GOODS Survey), provides deep, high-quality imaging in the K and Y bands over the CANDELS UDS and GOODS-South fields. We describe here the survey strategy, the data reduction process, and the data quality. HUGS delivers the deepest and highest quality K-band images ever collected over areas of cosmological interest, and ideally complements the CANDELS data set in terms of image quality and depth. The seeing is exceptional and homogeneous, confined to the range 0.38"-0.43". In the deepest region of the GOODS-S field, (which includes most of the HUDF) the K-band exposure time exceeds 80 hours of integration, yielding a 1-sigma magnitude limit of ~28.0 mag/sqarcsec. In the UDS field the survey matches the shallower depth of the CANDELS images reaching a 1-sigma limit per sq.arcsec of ~27.3mag in the K band and ~28.3mag in the Y-band, We show that the HUGS observations are well matched to the depth of the CANDELS WFC3/IR data, since the majority of even the faintest galaxies detected in the CANDELS H-band images are also detected in HUGS. We present the K-band galaxy number counts produced by combining the HUGS data from the two fields. We show that the slope of the number counts depends sensitively on the assumed distribution of galaxy sizes, with potential impact on the estimated extra-galactic background light (abridged).

The Herschel very wide-field surveys have charted hundreds of square degrees in multiple far-IR (FIR) bands. While the Sloan Digital Sky Survey (SDSS) is currently the best resource for optical counterpart identifications over such wide areas, it does not detect a large number of Herschel FIR sources and leaves their nature undetermined. As a test case, we studied seven "SDSS-invisible", very bright 250um sources (S_{250} > 55 mJy) in the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) fields where we have a rich multi-wavelength data set. We took a new approach to decompose the FIR sources, using the near-IR or the optical images directly for position priors. This is an improvement over the previous decomposition efforts where the priors are from mid-IR data that still suffer from the source blending problem in the first place. We found that in most cases the single Herschel sources are made of multiple components that are not necessarily at the same redshifts. Our decomposition succeeded in identifying and extracting their major contributors. We show that these are all ULIRGs at z ~ 1--2 whose high L_IR is mainly due to dust-obscured star formation. Most of them would not be selected as sub-mm galaxies. They all have complicated morphologies indicative of merger or violent instability, and their stellar populations are heterogeneous in terms of stellar masses, ages and formation histories. Their current ULIRG phases are of various degrees of importance in their stellar mass assembly. Our practice provides a promising starting point to develop an automatic routine to reliably study bright Herschel sources.

We derive physical properties of 10 submillimeter galaxies located in the CANDELS coverage of the GOODS-S field. The galaxies were first identified as submillimeter sources with the LABOCA bolometer and subsequently targeted for 870um continuum observation with ALMA. The high angular resolution of the ALMA imaging allows secure counterparts to be identified in the CANDELS multiband dataset. The CANDELS data provide deep photometric data from UV through near-infrared wavelengths. Using synthetic spectral energy distributions, we derive photometric redshifts, stellar masses, extinction, ages, and the star formation history. The redshift range is z=1.65-4.76, with two of the galaxies located at z>4. Two SMG counterparts have stellar masses 2-3 orders of magnitude lower than the rest. The remaining SMG counterparts have stellar masses around 1x10^11 Msun. The stellar population in the SMGs is typically older than the expected duration of the submillimeter phase, suggesting that the star formation history of submillimeter galaxies is more complex than a single burst. Non-parametric morphology indices suggest that the SMG counterparts are among the most asymmetric systems compared with galaxies of the same stellar mass and redshift. The HST images shows that 3 of the SMGs are associated with on-going mergers. The remaining counterparts are isolated. Estimating the dust and molecular gas mass from the submm fluxes, and comparing with our stellar masses shows that the molecular gas mass fraction of SMGs is ~28% and that the final stellar mass is likely to be (1-2)x10^11 Msun.

We revisit propagation of galactic cosmic rays (CRs) in light of recent advances in CR diffusion theory in realistic interstellar turbulence. We use a tested model of turbulence in which it has been shown that fast modes dominate scattering of CRs. As a result, propagation becomes inhomogeneous and environment dependent. By adopting the formalism of the nonlinear theory developed by Yan & Lazarian, we calculate the diffusion of CRs self-consistently from first principles. We assume a two- phase model for the Galaxy to account for different damping mechanisms of the fast modes, and we find that the energy dependence of the diffusion coefficient is mainly affected by medium properties. We show that it gives a correct framework to interpret some of the recent CR puzzles.

Recent advances in both the MHD turbulence theory and cosmic ray observations call for revisions in the paradigm of cosmic ray transport. We use the models of magnetohydrodynamic turbulence that were tested in numerical simulation, in which turbulence is injected at large scale and cascades to to small scales. We shall present the nonlinear results for cosmic ray transport, in particular, the cross field transport of CRs and demonstrate that the concept of cosmic ray subdiffusion in general does not apply and the perpendicular motion is well described by normal diffusion with M_A^4 dependence. Moreover, on scales less than injection scale of turbulence, CRs' transport becomes super-diffusive. Quantitative predictions for both the normal diffusion on large scale and super diffusion are confronted with recent numerical simulations. Implication for shock acceleration is briefly discussed.

A correct description of cosmic-ray (CR) diffusion in turbulent plasma is essential for many astrophysical and heliospheric problems. This paper aims to present physical diffusion behavior of CRs in actual turbulent magnetic fields, model of which has been numerically tested. We perform test particle simulations in compressible magnetohydrodynamic turbulence. We obtain scattering and spatial diffusion coefficients by tracing particle trajectories. We find no resonance gap for pitch-angle scattering at 90$^\circ$. Our result confirms the dominance of mirror interaction with compressible modes for most pitch angles as revealed by the nonlinear theory. For cross-field transport, our results are consistent with normal diffusion predicted earlier for large scales. The diffusion behavior strongly depends on the Alfvenic Mach number and particle's parallel mean free path. We, for the first time, numerically derive the dependence of M_A^4 for perpendicular diffusion coefficient with respect to the mean magnetic field. We conclude that CR diffusion coefficients are anisotropic in sub-Alfvenic turbulence and spatially correlated to the local turbulence properties. On scales smaller than the injection scale, we find that CRs are superdiffusive. We emphasize the importance of our results in a wide range of astrophysical processes, including magnetic reconnection.

Cosmic ray (CR) transport and acceleration is determined by the properties of magnetic turbulence. Recent advances in MHD turbulence call for revisions in the paradigm of cosmic ray transport. We use the models of magnetohydrodynamic turbulence that were tested in numerical simulation, in which turbulence is injected at large scale and cascades to to small scales. We shall address the issue of the transport of CRs, both parallel and perpendicular to the magnetic field. We shall demonstrate compressible fast modes are dominant cosmic ray scatterer from both quasilinear and nonlinear theories. We shall also show that the self-generated wave growth by CRs are constrained by preexisting turbulence and discuss the process in detail in the context of shock acceleration at supernova remnants and their implications. In addition, we shall dwell on the nonlinear growth of kinetic gyroresonance instability of cosmic rays induced by large scale compressible turbulence. This gyroresonance of cosmic rays on turbulence is demonstrated an important scattering mechanism in addition to direct interaction with the compressible turbulence. The feedback of the instability on large scale turbulence should be included in future simulations.

We fit the most updated broadband spectral energy distribution from radio to X-rays for NGC 4258 with a coupled accretion-jet model that surrounding a Kerr black hole (BH), where both the jet and the warped H_2O maser disk are assumed to be triggered by a spinning BH through Blandford-Znajek mechanism and Bardeen-Petterson effect respectively. The accretion flow consists with an inner radiatively inefficient accretion flow (RIAF) and an outer truncated standard thin disk, where the transition radius R_tr~3*10^3Rg for NGC 4258 based on the width and variability of its narrow Fe K$\alpha$ line. The hybrid jet formation model, as a variant of Blandford-Znajek model, is used to model the jet power. Therefore, we can estimate the accretion rate and BH spin through the two observed quantities--X-ray emission and jet power, where the observed jet power is estimated from the low-frequency radio emission. Through this method, we find that the BH of NGC 4258 should be mildly spinning with dimensionless spin parameter a_*=0.7\pm0.2. The outer thin disk mainly radiates at near infrared waveband and the jet contributes predominantly at radio waveband. Using above estimated BH spin and the inferred accretion rate at the region of the maser disk based on the physical existence of the H_2 O maser, we find that the warp radius is ~8.6*10^4 R_g if it is driven by the Bardeen-Petterson effect, which is consistent with the observational result very well.

We present observations of a luminous galaxy at redshift z=6.573 --- the end of the reioinization epoch --- which has been spectroscopically confirmed twice. The first spectroscopic confirmation comes from slitless HST ACS grism spectra from the PEARS survey (Probing Evolution And Reionization Spectroscopically), which show a dramatic continuum break in the spectrum at restframe 1216 A wavelength. The second confirmation is done with Keck + DEIMOS. The continuum is not clearly detected with ground-based spectra, but high wavelength resolution enables the Lyman alpha emission line profile to be determined. We compare the line profile to composite line profiles at redshift z=4.5. The Lyman alpha line profile shows no signature of a damping wing attenuation, confirming that the intergalactic gas is ionized at redshift z=6.57. Spectra of Lyman breaks at yet higher redshifts will be possible using comparably deep observations with IR-sensitive grisms, even at redshifts where Lyman alpha is too attenuated by the neutral IGM to be detectable using traditional spectroscopy from the ground.

We discuss how the model of magnetic reconnection in the presence of turbulence proposed inLazarian & Vishniac 1999 makes the reconnection rate independent either of resistivity or microscopic plasma effects, but determined entirely by the magnetic field line wandering induced by turbulence. We explain that the model accounts for both fast and slow regimes of reconnection and that this property naturally induces flares of reconnection in low beta plasma environments. In addition, we show that the model involves volume reconnection which can convert a substantial part of the energy into energetic particles. It is important that the reconnection induces an efficient acceleration of the first order Fermi type. Finally, we relate the properties of the reconnection with the observed properties of gamma ray bursts and provide evidence supporting the explanation of gamma ray bursts based on energy release via reconnection.

Observational studies of magnetic fields are crucial. We introduce a process "ground state alignment" as a new way to determine the magnetic field direction in diffuse medium. The alignment is due to anisotropic radiation impinging on the atom/ion. The consequence of the process is the polarization of spectral lines resulting from scattering and absorption from aligned atomic/ionic species with fine or hyperfine structure. The magnetic field induces precession and realign the atom/ion and therefore the polarization of the emitted or absorbed radiation reflects the direction of the magnetic field. The atoms get aligned at their low levels and, as the life-time of the atoms/ions we deal with is long, the alignment induced by anisotropic radiation is susceptible to extremely weak magnetic fields ($1{\rm G}\gtrsim B\gtrsim 10^{-15}$G). In fact, the effects of atomic/ionic alignment were studied in the laboratory decades ago, mostly in relation to the maser research. Recently, the atomic effect has been already detected in observations from circumstellar medium and this is a harbinger of future extensive magnetic field studies. A unique feature of the atomic realignment is that they can reveal the 3D orientation of magnetic field. In this article, we shall review the basic physical processes involved in atomic realignment. We shall also discuss its applications to interplanetary, circumstellar and interstellar magnetic fields. In addition, our research reveals that the polarization of the radiation arising from the transitions between fine and hyperfine states of the ground level can provide a unique diagnostics of magnetic fields in the Epoch of Reionization.

We analyze the spectral energy distributions (SEDs) of Lyman break galaxies (LBGs) at z=1-3 selected using the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) UVIS channel filters. These HST/WFC3 observations cover about 50 sq. arcmin in the GOODS-South field as a part of the WFC3 Early Release Science program. These LBGs at z=1-3 are selected using dropout selection criteria similar to high redshift LBGs. The deep multi-band photometry in this field is used to identify best-fit SED models, from which we infer the following results: (1) the photometric redshift estimate of these dropout selected LBGs is accurate to within few percent; (2) the UV spectral slope (beta) is redder than at high redshift (z>3), where LBGs are less dusty; (3) on average, LBGs at z=1-3 are massive, dustier and more highly star-forming, compared to LBGs at higher redshifts with similar luminosities (0.1L*<~L<~2.5L*), though their median values are similar within 1-sigma uncertainties. This could imply that identical dropout selection technique, at all redshifts, find physically similar galaxies; and (4) stellar masses of these LBGs are directly proportional to their UV luminosities with a logarithmic slope of ~0.46, and star-formation rates are proportional to their stellar masses with a logarithmic slope of ~0.90. These relations hold true --- within luminosities probed in this study --- for LBGs from z~1.5 to 5. The star-forming galaxies selected using other color-based techniques show similar correlations at z~2, but to avoid any selection biases, and for direct comparison with LBGs at z>3, a true Lyman break selection at z~2 is essential. The future HST UV surveys, both wider and deeper, covering a large luminosity range are important to better understand LBG properties, and their evolution.

One key goal of the Hubble Space Telescope Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey is to track galaxy evolution back to z ~ 8. Its two-tiered "wide and deep" strategy bridges significant gaps in existing near-infrared surveys. Here we report on z ~ 8 galaxy candidates selected as F105W-band dropouts in one of its deep fields, which covers 50.1 square arcmin to 4 ks depth in each of three near-infrared bands in the Great Observatories Origins Deep Survey southern field. Two of our candidates have J<26.2 mag, and are > 1 mag brighter than any previously known F105W-dropouts. We derive constraints on the bright-end of the rest-frame ultraviolet luminosity function of galaxies at z ~ 8, and show that the number density of such very bright objects is higher than expected from the previous Schechter luminosity function estimates at this redshift. Another two candidates are securely detected in Spitzer Infrared Array Camera images, which are the first such individual detections at z ~ 8. Their derived stellar masses are on the order of a few x 10^9 M_sun, from which we obtain the first measurement of the high-mass end of the galaxy stellar mass function at z ~ 8. The high number density of very luminous and very massive galaxies at z ~ 8, if real, could imply a large stellar-to-halo mass ratio and an efficient conversion of baryons to stars at such an early time.

Turbulence is ubiquitous in astrophysics. It radically changes many astrophysical phenomena, in particular, the propagation and acceleration of cosmic rays. We present the modern understanding of compressible magnetohydrodynamic (MHD) turbulence, in particular its decomposition into Alfv\'en, slow and fast modes, discuss the density structure of turbulent subsonic and supersonic media, as well as other relevant regimes of astrophysical turbulence. All this information is essential for understanding the energetic particle acceleration that we discuss further in the review. For instance, we show how fast and slow modes accelerate energetic particles through the second order Fermi acceleration, while density fluctuations generate magnetic fields in pre-shock regions enabling the first order Fermi acceleration of high energy cosmic rays. Very importantly, however, the first order Fermi cosmic ray acceleration is also possible in sites of magnetic reconnection. In the presence of turbulence this reconnection gets fast and we present numerical evidence supporting the predictions of the Lazarian & Vishniac (1999) model of fast reconnection. The efficiency of this process suggests that magnetic reconnection can release substantial amounts of energy in short periods of time. As the particle tracing numerical simulations show that the particles can be efficiently accelerated during the reconnection, we argue that the process of magnetic reconnection may be much more important for particle acceleration than it is currently accepted. In particular, we discuss the acceleration arising from reconnection as a possible origin of the anomalous cosmic rays measured by Voyagers as well as the origin cosmic ray excess in the direction of Heliotail.

We present ultra-deep J and Ks imaging observations covering a 30' * 30' area of the Extended Chandra Deep Field-South (ECDFS) carried out by our Taiwan ECDFS Near-Infrared Survey (TENIS). The median 5-sigma limiting magnitudes for all detected objects in the ECDFS reach 24.5 and 23.9 mag (AB) for J and Ks, respectively. In the inner 400 arcmin^2 region where the sensitivity is more uniform, objects as faint as 25.6 and 25.0 mag are detected at 5-sigma. So this is by far the deepest J and Ks datasets available for the ECDFS. To combine the TENIS with the Spitzer IRAC data for obtaining better spectral energy distributions of high-redshift objects, we developed a novel deconvolution technique (IRACLEAN) to accurately estimate the IRAC fluxes. IRACLEAN can minimize the effect of blending in the IRAC images caused by the large point-spread functions and reduce the confusion noise. We applied IRACLEAN to the images from the Spitzer IRAC/MUSYC Public Legacy in the ECDFS survey (SIMPLE) and generated a J+Ks selected multi-wavelength catalog including the photometry of both the TENIS near-infrared and the SIMPLE IRAC data. We publicly release the data products derived from this work, including the J and Ks images and the J+Ks selected multiwavelength catalog.

We demonstrate a new way of studying interplanetary magnetic field -- Ground State Alignment (GSA). Instead of sending thousands of space probes, GSA allows magnetic mapping with any ground telescope facilities equipped with spectropolarimeter. The polarization of spectral lines that are pumped by the anisotropic radiation from the Sun is influenced by the magnetic realignment, which happens for magnetic field (<1G). As a result, the linear polarization becomes an excellent tracer of the embedded magnetic field. The method is illustrated by our synthetic observations of the Jupiter's Io and comet Halley. Polarization at each point was constructed according to the local magnetic field detected by spacecrafts. Both spatial and temporal variations of turbulent magnetic field can be traced with this technique as well. The influence of magnetic field on the polarization of scattered light is discussed in detail. For remote regions like the IBEX ribbons discovered at the boundary of interstellar medium, GSA provides a unique diagnostics of magnetic field.

We report on deep near-infrared F125W (J) and F160W (H) Hubble Space Telescope Wide Field Camera 3 images of the z=6.42 quasar J1148+5251 to attempt to detect rest-frame near-ultraviolet emission from the host galaxy. These observations included contemporaneous observations of a nearby star of similar near-infrared colors to measure temporal variations in the telescope and instrument point spread function (PSF). We subtract the quasar point source using both this direct PSF and a model PSF. Using direct subtraction, we measure an upper limit for the quasar host galaxy of m_J>22.8, m_H>23.0 AB mag (2 sigma). After subtracting our best model PSF, we measure a limiting surface brightness from 0.3"-0.5" radius of mu_J > 23.5, mu_H > 23.7 AB magarc (2 sigma). We test the ability of the model subtraction method to recover the host galaxy flux by simulating host galaxies with varying integrated magnitude, effective radius, and S\'ersic index, and conducting the same analysis. These models indicate that the surface brightness limit (mu_J > 23.5 AB magarc) corresponds to an integrated upper limit of m_J > 22 - 23 AB mag, consistent with the direct subtraction method. Combined with existing far-infrared observations, this gives an infrared excess log(IRX) > 1.0 and corresponding ultraviolet spectral slope beta > -1.2\pm0.2. These values match those of most local luminous infrared galaxies, but are redder than those of almost all local star-forming galaxies and z~6 Lyman break galaxies.

Using a sample of BzK-selected galaxies at z~2 identified from the CFHT/WIRCAM near-infrared survey of GOODS-North, we discuss the relation between star formation rate (SFR), specific star formation rate (SSFR), and stellar mass (M_{*}), and the clustering of galaxies as a function of these parameters. For star-forming galaxies (sBzKs), the UV-based SFR, corrected for extinction, scales with the stellar mass as SFR ~ M_{*}^{alpha} with alpha = 0.74+/-0.20 down to M_{*} ~ 10^{9} M_{solar}, indicating a weak dependence on the stellar mass of the SSFR. We also measure the angular correlation function and hence infer the correlation length for sBzK galaxies as a function of M_{*}, SFR, and SSFR, as well as K-band apparent magnitude. We show that passive galaxies (pBzKs) are more strongly clustered than sBzK galaxies at a given stellar mass, mirroring the color-density relation seen at lower redshifts. We also find that the correlation length of sBzK galaxies ranges from 4 to 20 h^{-1}Mpc, being a strong function of M_{K}, M_{*}, and SFR. On the other hand, the clustering dependence on SSFR changes abruptly at 2x10^{-9} yr^{-1}, which is the typical value for "main sequence" star-forming galaxies at z~2. We show that the correlation length reaches a minimum at this characteristic value, and is larger for galaxies with both smaller and larger SSFRs; a dichotomy that is only marginally implied from the predictions of the semi-analytical models. Our results suggest that there are two types of environmental effects at work at z~2. Stronger clustering for relatively quiescent galaxies implies that the environment has started to play a role in quenching star formation. At the same time, stronger clustering for galaxies with elevated SSFRs ("starbursts") might be attributed to an increased efficiency for galaxy interactions and mergers in dense environments.

Observational studies of magnetic fields are vital as magnetic fields play a crucial role in various astrophysical processes, including star formation, accretion of matter, transport processes (e.g., transport of heat), and cosmic rays. We identified a process "ground state alignment" as a new way to determine the magnetic field direction in diffuse medium. The alignment is due to anisotropic radiation impinging on the atom/ion, while the magnetic field induces precession and realign the atom/ion and therefore the polarization of the emitted or absorbed radiation reflects the direction of the magnetic field. The atoms get aligned at their low levels and, as the life-time of the atoms/ions we deal with is long, the alignment induced by anisotropic radiation is susceptible to weak magnetic fields ($1{\rm G}\gtrsim B\gtrsim 10^{-15}$G). Compared to the upper level Hanle effect, atomic realignment is most suitable for the studies of magnetic field in the diffuse medium, where magnetic field is relatively weak. In fact, the effects of atomic/ionic alignment, including the realignment in magnetic field, were studied in the laboratory decades ago, mostly in relation to the maser research. Recently, the atomic effect has been already detected in observations from circumstellar medium and this is a harbinger of future extensive magnetic field studies. A unique feature of the atomic realignment is that they can reveal the 3D orientation of magnetic field. In this article, we shall review the basic physical processes involved in atomic realignment and its applications to interplanetary, circumstellar and interstellar magnetic fields. In addition, our research reveals that the polarization of the radiation arising from the transitions between fine and hyperfine states of the ground level can provide a unique diagnostics of magnetic fields, including those in the Early Universe.

We present a study on Spectral Energy Distributions, Morphologies, and star formation for an IRAC-selected extremely red object sample in the GOODS Chandra Deep Field-South. This work was enabled by new HST/WFC3 near-IR imaging from the CANDELS survey as well as the deepest available X-ray data from Chandra 4 Ms observations. This sample consists of 133 objects with the 3.6um limiting magnitude of [3.6] = 21.5, and is approximately complete for galaxies with M >10^{11}M_sun at 1.5 < z < 2.5. We classify this sample into two types, quiescent and star-forming galaxies, in the observed infrared color-color ([3.6]-[24] vs K-[3.6]) diagram. The further morphological study of this sample show a consistent result with the observed color classification. The classified quiescent galaxies are bulge dominated and star-forming galaxies in the sample have disk or irregular morphologies. Our observed infrared color classification is also consistent with the rest-frame color (U-V vs V-J) classification. We also found that quiescent and star-forming galaxies are well separated in the nonparametric morphology parameter (Gini vs M_{20}) diagram measuring their concentration and clumpiness: quiescent galaxies have Gini coefficient higher than 0.58 and star forming galaxies have Gini coefficient lower that 0.58. We argue that the star formation quenching process must lead to or be accompanied by the increasing galaxy concentration. One prominent morphological feature of this sample is that disks are commonly seen in this massive galaxy sample at 1.5 < z < 2.5: 30% of quiescent galaxies and 70% of star forming galaxies with M >10^{11}M_sun have disks in their rest-frame optical morphologies. The prevalence of these extended, relatively undisturbed disks challenges the merging scenario as the main mode of massive galaxy formation.

Using volume-limited samples drawn from The Sloan Digital Sky Survey Data Release 7 (SDSS DR7), we measure the tidal environment of galaxies, which we characterize by the ellipticity e of the potential field calculated from the smoothed spatial number density 1 + {\delta} of galaxies. We analyze if galaxy properties, including color, Dn4000, concentration and size correlate with e, in addition to depending on 1 + {\delta}. We find that there exists a transition smoothing scale at which correlations/anti-correlations with e reverse. This transition scale is well represented by the distance to the 3rd- nearest-neighbor of a galaxy in a volume limited sample with Mr < -20 which has a distribution peaked at ~ 2 h-1Mpc. We further demonstrate that this scale corresponds to that where the correlation between the color of galaxies and environmental density 1+{\delta} is the strongest. For this optimal smoothing R0 no additional correlations with e are observed. The apparent dependence on tidal ellipticity e at other smoothing scales Rs can be viewed as a geometric effect, arising from the cross correlation be- tween (1+{\delta}o) and e(Rs). We perform the same analysis on numerical simulations with semi-analytical modeling (SAM) of galaxy formation. The e dependence of the galaxy properties shows similar behavior to that in the SDSS, although the color-density cor- relation is significantly stronger in the SAM. The 'optimal adaptive smoothing scale' in the SAM is also closely related to the distance to the 3rd-nearest-neighbor of a galaxy, and its characteristic value is consistent with, albeit slightly smaller than for SDSS.

The primary goal of the Taiwan ECDFS Near-Infrared Survey (TENIS) is to find well screened galaxy candidates at z>7 (z' dropout) in the Extended Chandra Deep Field-South (ECDFS). To this end, TENIS provides relatively deep J and Ks data (~25.3 ABmag, 5-sigma) for an area of 0.5*0.5 degree. Leveraged with existing data at mid-infrared to optical wavelengths, this allows us to screen for the most luminous high-z objects, which are rare and thus require a survey over a large field to be found. We introduce new color selection criteria to select a z>7 sample with minimal contaminations from low-z galaxies and Galactic cool stars; to reduce confusion in the relatively low angular resolution IRAC images, we introduce a novel deconvolution method to measure the IRAC fluxes of individual sources. Illustrating perhaps the effectiveness at which we screen out interlopers, we find only one z>7 candidate, TENIS-ZD1. The candidate has a weighted z_phot of 7.8, and its colors and luminosity indicate a young (45M years old) starburst galaxy with a stellar mass of 3.2*10^10 M_sun. The result matches with the observational luminosity function analysis and the semi-analytic simulation result based on the Millennium Simulations, which may over predict the volume density for high-z massive galaxies. The existence of TENIS-ZD1, if confirmed spectroscopically to be at z>7, therefore poses a challenge to current theoretical models for how so much mass can accumulate in a galaxy at such a high redshift.

In the first of a series of forthcoming publications, we present a panchromatic catalog of 102 visually-selected early-type galaxies (ETGs) from observations in the Early Release Science (ERS) program with the Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST) of the Great Observatories Origins Deep Survey-South (GOODS-S) field. Our ETGs span a large redshift range, 0.35 < z < 1.5, with each redshift spectroscopically-confirmed by previous published surveys of the ERS field. We combine our measured WFC3 ERS and ACS GOODS-S photometry to gain continuous sensitivity from the rest-frame far-UV to near-IR emission for each ETG. The superior spatial resolution of the HST over this panchromatic baseline allows us to classify the ETGs by their small-scale internal structures, as well as their local environment. By fitting stellar population spectral templates to the broad-band photometry of the ETGs, we determine that the average masses of the ETGs are comparable to the characteristic stellar mass of massive galaxies, 11< log(M [Solar]) < 12. By transforming the observed photometry into the GALEX FUV and NUV, Johnson V, and SDSS g' and r' bandpasses we identify a noteworthy diversity in the rest-frame UV-optical colors and find the mean rest-frame (FUV-V)=3.5 and (NUV-V)=3.3, with 1$\sigma$ standard deviations approximately equal to 1.0. The blue rest-frame UV-optical colors observed for most of the ETGs are evidence for star-formation during the preceding gigayear, but no systems exhibit UV-optical photometry consistent with major recent (<~50 Myr) starbursts. Future publications which address the diversity of stellar populations likely to be present in these ETGs, and the potential mechanisms by which recent star-formation episodes are activated, are discussed.

Cosmic ray (CR) transport and acceleration is essential for many astrophysical problems, e.g., CMB foreground, ionization of molecular clouds and all high energy phenomena. Recent advances in MHD turbulence call for revisions in the paradigm of cosmic ray transport. We use the models of magnetohydrodynamic turbulence that were tested in numerical simulation, in which turbulence is injected at large scale and cascades to to small scales. We shall address the issue of the transport of CRs, both parallel and perpendicular to the magnetic field and show that the issue of cosmic ray subdiffusion is only important for restricted cases when the ambient turbulence is far from that suggested by numerical simulations. Moreover, on scales less than injection scale of turbulence, CRs's transport becomes super-diffusive. We also shall discuss the nonlinear growth of kinetic gyroresonance instability of cosmic rays induced by large scale compressible turbulence. This gyroresonance of cosmic rays on turbulence is demonstrated an important scattering mechanism in addition to direct interaction with the compressible turbulence. The feedback of the instability on large scale turbulence cannot be neglected, and should be included in future simulations.

High-energy gamma ray emission has been detected recently from supernovae remnants (SNRs) and their surroundings. The existence of molecular clouds near some of the SNRs suggests that the gamma rays originate predominantly from p-p interactions with cosmic rays accelerated at a closeby SNR shock wave. Here we investigate the acceleration of cosmic rays and the gamma ray production in the cloud self-consistently by taking into account the interactions of the streaming instability and the background turbulence both at the shock front and in the ensuing propagation to the clouds. We focus on the later evolution of SNRs, when the conventional treatment of the streaming instability is valid but the magnetic field is enhanced due to either Bell's current instability and/or due to the dynamo generation of magnetic field in the precursor region. We calculate the time dependence of the maximum energy of the accelerated particles. This result is then used to determine the diffusive flux of the runaway particles escaping the shock region, from which we obtain the gamma spectrum consistent with observations. Finally, we check the self-consistency of our results by comparing the required level of diffusion with the level of the streaming instability attainable in the presence of turbulence damping. The energy range of cosmic rays subject to the streaming instability is able to produce the observed energy spectrum of gamma rays.

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at $z\sim1.5-8$, and to study Type Ia SNe beyond $z>1.5$. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) is designed to document the first third of galactic evolution, over the approximate redshift (z) range 8--1.5. It will image >250,000 distant galaxies using three separate cameras on the Hubble Space Telescope, from the mid-ultraviolet to the near-infrared, and will find and measure Type Ia supernovae at z>1.5 to test their accuracy as standardizable candles for cosmology. Five premier multi-wavelength sky regions are selected, each with extensive ancillary data. The use of five widely separated fields mitigates cosmic variance and yields statistically robust and complete samples of galaxies down to a stellar mass of 10^9 M_\odot to z \approx 2, reaching the knee of the ultraviolet luminosity function (UVLF) of galaxies to z \approx 8. The survey covers approximately 800 arcmin^2 and is divided into two parts. The CANDELS/Deep survey (5\sigma\ point-source limit H=27.7 mag) covers \sim 125 arcmin^2 within GOODS-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (EGS, COSMOS, and UDS) and covers the full area to a 5\sigma\ point-source limit of H \gtrsim 27.0 mag. Together with the Hubble Ultra Deep Fields, the strategy creates a three-tiered "wedding cake" approach that has proven efficient for extragalactic surveys. Data from the survey are nonproprietary and are useful for a wide variety of science investigations. In this paper, we describe the basic motivations for the survey, the CANDELS team science goals and the resulting observational requirements, the field selection and geometry, and the observing design. The Hubble data processing and products are described in a companion paper.

We present an analysis of the size growth seen in early-type galaxies over 10 Gyr of cosmic time. Our analysis is based on a homogeneous synthesis of published data from 17 spectroscopic surveys observed at similar spatial resolution, augmented by new measurements for galaxies in the Gemini Deep Deep Survey. In total, our sample contains structural data for 465 galaxies (mainly early-type) in the redshift range 0.2<z<2.7. The size evolution of passively-evolving galaxies over this redshift range is gradual and continuous, with no evidence for an end or change to the process around z~1, as has been hinted at by some surveys which analyze subsets of the data in isolation. The size growth appears to be independent of stellar mass, with the mass-normalized half-light radius scaling with redshift as R_e (1+z)^(-1.62 +/- 0.34). Surprisingly, this power law seems to be in good agreement with the recently reported continuous size evolution of UV-bright galaxies in the redshift range z~0.5-3.5. It is also in accordance with the predictions from recent theoretical models.

We present a sample of 17 newly discovered ultracool dwarf candidates later than ~M8, drawn from 231.90 arcmin2 of {\it Hubble Space Telescope} Wide Field Camera 3 infrared imaging. By comparing the observed number counts for 17.5<J_125<25.5 AB mag to an exponential disk model, we estimate a vertical scale height of z_scl=290 +- 25 (random) +- 30 (systematic) pc for a binarity fraction of f_b=0. While our estimate is roughly consistent with published results, we suggest that the differences can be attributed to sample properties, with the present sample containing far more substellar objects than previous work. We predict the object counts should peak at J_{125}~24 AB mag due to the exponentially-declining number density at the edge of the disc. We conclude by arguing that trend in scale height with spectral type may breakdown for brown dwarfs since they do not settle onto the main sequence.

We study the nonlinear growth of kinetic gyroresonance instability of cosmic rays (CRs) induced by large scale compressible turbulence. This feedback of cosmic rays on turbulence was shown to induce an important scattering mechanism in addition to direct interaction with the compressible turbulence. The linear growth is bound to saturate due to the wave-particle interactions. By balancing increase of CR anisotropy via the large scale compression and its decrease via the wave-particle scattering, we find the steady state solutions. The nonlinear suppression due to the wave-particle scattering limit the energy range of CRs that can excite the instabilities and be scattered by the induced slab waves. The direct interaction with large scale compressible modes still appears to be the dominant mechanism for isotropization of high energy cosmic rays (> 100 GeV).

We describe the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) Early Release Science (ERS) observations in the Great Observatories Origins Deep Survey (GOODS) South field. The new WFC3 ERS data provide calibrated, drizzled mosaics in the UV filters F225W, F275W, and F336W, as well as in the near-IR filters F098M (Ys), F125W (J), and F160W (H) with 1-2 HST orbits per filter. Together with the existing HST Advanced Camera for Surveys (ACS) GOODS-South mosaics in the BViz filters, these panchromatic 10-band ERS data cover 40-50 square arcmin at 0.2-1.7 {\mu}m in wavelength at 0.07-0.15" FWHM resolution and 0.090" Multidrizzled pixels to depths of AB\simeq 26.0-27.0 mag (5-{\sigma}) for point sources, and AB\simeq 25.5-26.5 mag for compact galaxies. In this paper, we describe: a) the scientific rationale, and the data taking plus reduction procedures of the panchromatic 10-band ERS mosaics; b) the procedure of generating object catalogs across the 10 different ERS filters, and the specific star-galaxy separation techniques used; and c) the reliability and completeness of the object catalogs from the WFC3 ERS mosaics. The excellent 0.07-0.15" FWHM resolution of HST/WFC3 and ACS makes star- galaxy separation straightforward over a factor of 10 in wavelength to AB\simeq 25-26 mag from the UV to the near-IR, respectively.

We present ASAS data starting 25 days before the discovery of the recent type IIn SN 2010jl, and we compare its light curve to other luminous IIn SNe, showing that it is a luminous (M_I ~ -20.5) event. Its host galaxy, UGC 5189, has a low gas-phase oxygen abundance (12 + log(O/H) = 8.2), which reinforces the emerging trend that over-luminous core-collapse supernovae are found in the low-metallicity tail of the galaxy distribution, similar to the known trend for the hosts of long GRBs. We compile oxygen abundances from the literature and from our own observations of UGC 5189, and we present an unpublished spectrum of the luminous type Ic SN 2010gx that we use to estimate its host metallicity. We discuss these in the context of host metallicity trends for different classes of core-collapse objects. The earliest generations of stars are known to be enhanced in [O/Fe] relative to the Solar mixture; it is therefore likely that the stellar progenitors of these overluminous supernovae are even more iron-poor than they are oxygen-poor. A number of mechanisms and massive star progenitor systems have been proposed to explain the most luminous core-collapse supernovae; any successful theory will need to include the emerging trend that points towards low-metallicity for the massive progenitor stars. This trend for very luminous supernovae to strongly prefer low-metallicity galaxies should be taken into account when considering various aspects of the evolution of the metal-poor early universe. (abridged)

We study large-scale structures from numerical simulations, paying particular attention to supercluster-like structures. A grid-density-contour based algorithm is adopted to locate connected groups. With the increase of the linking density threshold, the foam- like cosmic web is subsequently broken into individual supercluster-like groups and further halos which are in accordance to groups with the linking density threshold {\rho}/{\rho}= 1 + {\delta} = 80. By analyzing sets of cosmological simulations with varying cosmological parameters, we find that an universal mass function exists not only for halos but also for low-density supercluster-like groups until the linking density threshold decreases to a density where the global percolation of large-scale structures occurs. We further show that the mass functions of different groups can be well described by the Jenkins form with the parameters being dependent on the linking density threshold. On the other hand, these low- density supercluster-like groups cannot be directly associated with the predictions from the excursion set theory with effective barriers obtained from dynamical collapse models, and the peak exclusion effect must be taken into account. Including such an effect, the mass function of groups with the linking density threshold 1 + {\delta} = 16 is in good agreements with that from the excursion set theory with a nearly flat effective barrier. A simplified analysis of the ellipsoidal collapse model indicates that the barrier for collapses along two axes to form filaments is approximately flat in scales, thus we define groups identified with 1 + {\delta} = 16 as filaments. We then further study the halo-filament conditional mass function and the filament-halo conditional mass function, and compare them with the predictions from the two-barrier excursion set theory. The shape statistics for filaments are also presented.

The observed number counts of high-redshift galaxy candidates have been used to build up a statistical description of star-forming activity at redshift z >~ 7, when galaxies reionized the Universe. Standard models predict that a high incidence of gravitational lensing will probably distort measurements of flux and number of these earliest galaxies. The raw probability of this happening has been estimated to be ~ 0.5 percent, but can be larger owing to observational biases. Here we report that gravitational lensing is likely to dominate the observed properties of galaxies with redshifts of z >~ 12, when the instrumental limiting magnitude is expected to be brighter than the characteristic magnitude of the galaxy sample. The number counts could be modified by an order of magnitude, with most galaxies being part of multiply imaged systems, located less than 1 arcsec from brighter foreground galaxies at z ~ 2. This lens-induced association of high-redshift and foreground galaxies has perhaps already been observed among a sample of galaxy candidates identified at z ~ 10.6. Future surveys will need to be designed to account for a significant gravitational lensing bias in high-redshift galaxy samples.

We report the first results from the Hubble Infrared Pure Parallel Imaging Extragalactic Survey, which utilizes the pure parallel orbits of the Hubble Space Telescope to do deep imaging along a large number of random sightlines. To date, our analysis includes 26 widely separated fields observed by the Wide Field Camera 3, which amounts to 122.8 sq.arcmin in total area. We have found three bright Y098-dropouts, which are candidate galaxies at z >~ 7.4. One of these objects shows an indication of peculiar variability and its nature is uncertain. The other two objects are among the brightest candidate galaxies at these redshifts known to date L>2L*. Such very luminous objects could be the progenitors of the high-mass Lyman break galaxis (LBGs) observed at lower redshifts (up to z~5). While our sample is still limited in size, it is much less subject to the uncertainty caused by "cosmic variance" than other samples because it is derived using fields along many random sightlines. We find that the existence of the brightest candidate at z~7.4 is not well explained by the current luminosity function (LF) estimates at z~8. However, its inferred surface density could be explained by the prediction from the LFs at z~7 if it belongs to the high-redshift tail of the galaxy population at z~7.

The recent Fermi observation of GRB 080916C shows that the bright photosphere emission associated with a putative fireball is missing, which suggests a Poynting-flux-dominated outflow. We propose a model of gamma-ray burst (GRB) prompt emission in the Poynting-flux-dominated regime, namely, the Internal-Collision-induced MAgnetic Reconnection and Turbulence (ICMART) model. It is envisaged that the GRB central engine launches an intermittent, magnetically-dominated wind, and that in the GRB emission region, the ejecta is still moderately magnetized. Similar to the internal shock (IS) model, the mini-shells interact internally at the traditional internal shock radius. Most of these early collision have little energy dissipation, but serve to distort the ordered magnetic field lines. At a certain point, the distortion of magnetic field configuration reaches the critical condition to allow fast reconnection seeds to occur, which induce relativistic MHD turbulence in the interaction regions. The turbulence further distorts field lines easing additional magnetic reconnections, resulting in a runway release of the stored magnetic field energy (an ICMART event). Particles accelerated in the ICMART region radiate synchrotron photons that power the observed gamma-rays. Each ICMART event corresponds to a broad pulse in the GRB lightcurve, and a GRB is composed of multiple ICMART events. This model retains the merits of the IS and other models, but may overcome several difficulties/issues faced by the IS model (e.g. low efficiency, fast cooling, electron number excess, Amati/Yonetoku relation inconsistency, and missing bright photosphere). It predicts two-component variability time scales, and a decreasing Ep and polarization degree during each ICMART event. The model may be applied to most Fermi LAT GRBs that have time-resolved, featureless Band-function spectra (abridged).

We present the details and early results from a deep near-infrared survey utilising the NICMOS instrument on the Hubble Space Telescope centred around massive M_* > 10^11 M_0 galaxies at 1.7 < z < 2.9 found within the Great Observatories Origins Deep Survey (GOODS) fields. The GOODS NICMOS Survey (GNS) was designed to obtain deep F160W (H-band) imaging of 80 of these massive galaxies, as well as other colour selected objects such as Lyman-break drop-outs, BzK objects, Distant Red Galaxies, EROs, Spitzer Selected EROs, BX/BM galaxies, as well as sub-mm galaxies. We present in this paper details of the observations, our sample selection, as well as a description of features of the massive galaxies found within our survey fields. This includes: photometric redshifts, rest-frame colours, and stellar masses. We furthermore provide an analysis of the selection methods for finding massive galaxies at high redshifts, including colour selection, and how galaxy populations selected through different methods overlap. We find that a single colour selection method cannot locate all of the massive galaxies, with no one method finding more than 70 percent. We however find that the combination of these colour methods finds nearly all the massive galaxies, as selected by photometric redshifts with the exception of apparently rare blue massive galaxies. By investigating the rest-frame (U-B) vs. M_B diagram for these galaxies we furthermore show that there exists a bimodality in colour-magnitude space at z < 2, driven by stellar mass, such that the most massive galaxies are systematically red up to z~2.5, while lower mass galaxies tend to be blue. We also discuss the number densities for galaxies with stellar masses M_* > 10^11 M_0, whereby we find an increase of a factor of eight between z = 3 and z = 1.5, demonstrating that this is an epoch when massive galaxies establish most of their mass.

We combine new high sensitivity ultraviolet (UV) imaging from the Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST) with existing deep HST/Advanced Camera for Surveys (ACS) optical images from the Great Observatories Origins Deep Survey (GOODS) program to identify UV-dropouts, which are Lyman break galaxy (LBG) candidates at z~1-3. These new HST/WFC3 observations were taken over 50 sq.arcmin in the GOODS-South field as a part of the Early Release Science program. The uniqueness of these new UV data is that they are observed in 3 UV/optical (WFC3 UVIS) channel filters (F225W, F275W and F336W), which allows us to identify three different sets of UV-dropout samples. We apply Lyman break dropout selection criteria to identify F225W-, F275W- and F336W-dropouts, which are z~1.7, 2.1 and 2.7 LBG candidates, respectively. Our results are as follows: (1) these WFC3 UVIS filters are very reliable in selecting LBGs with z~2.0, which helps to reduce the gap between the well studied z~>3 and z~0 regimes, (2) the combined number counts agrees very well with the observed change in the surface densities as a function of redshift when compared with the higher redshift LBG samples; and (3) the best-fit Schechter function parameters from the rest-frame UV luminosity functions at three different redshifts fit very well with the evolutionary trend of the characteristic absolute magnitude, and the faint-end slope, as a function of redshift. This is the first study to illustrate the usefulness of the WFC3 UVIS channel observations to select z<3 LBGs. The addition of the new WFC3 on the HST has made it possible to uniformly select LBGs from z~1 to z~9, and significantly enhance our understanding of these galaxies using HST sensitivity and resolution.

We present a large sample of candidate galaxies at z~7--10, selected in the HUDF using the new observations made by the HST/WFC3. Our sample is composed of 20 z-dropouts, 15 Y-dropouts, and 20 J-dropouts. The surface densities of the z-dropouts are close to what predicted by earlier studies, however, those of the Y- and J-dropouts are quite unexpected. While no Y- or J-dropouts have been found at AB < 28.0 mag, their surface densities seem to increase sharply at fainter levels. While some of these candidates seem to be close to foreground galaxies and thus could possibly be gravitationally lensed, the overall surface densities after excluding such cases are still much higher than what would be expected if the luminosity function does not evolve from z~7 to 10. Motivated by such steep increases, we tentatively propose a set of Schechter function parameters to describe the LFs at z~8 and 10. As compared to their counterpart at z~7, here L* decreases by ~ 6.5x and Phi* increases by 17--90x. Although such parameters are not yet demanded by the existing observations, they are allowed and seem to agree with the data better than other alternatives. If these LFs are still valid beyond our current detection limit, this would imply a sudden emergence of a large number of low-luminosity galaxies when looking back in time to z~10, which, while seemingly exotic, would naturally fit in the picture of the cosmic hydrogen reionization. These early galaxies could easily account for the ionizing photon budget required by the reionization, and they would imply that the global star formation rate density might start from a very high value at z~10, rapidly reach the minimum at z~7, and start to rise again towards z~6. In this scenario, the majority of the stellar mass that the universe assembled through the reionization epoch seems still undetected by current observations at z~6. [Abridged]

We present results on the size evolution of passively evolving galaxies at 1<z<2 drawn from the Wide Field Camera 3 Early Release Science program. Our sample was constructed using an analog to the passive BzK selection criterion, which isolates galaxies with little or no on-going star formation at z>1.5. We identify 30 galaxies in ~40 square arcmin to H<25 mag. We supplement spectroscopic redshifts from the literature with photometric redshifts determined from the 15-band photometry from 0.22-8 micron. We determine effective radii from Sersic profile fits to the H-band image using an empirical PSF. We find that size evolution is a strong function of stellar mass, with the most massive (M* ~ 10^11 Msol) galaxies undergoing the most rapid evolution from z~2 to the present. Parameterizing the size evolution as (1+z)^{-alpha}, we find a tentative scaling between alpha and stellar mass of alpha ~ -1.8+1.4 log(M*/10^9 Msol). We briefly discuss the implications of this result for our understanding of the dynamical evolution of the red galaxies.

We report the discovery of eight new quasars in one extragalactic field (five degree centered at RA=$08^h58^m08.2^s$, Dec=$01^o32'29.7''$) with the LAMOST commissioning observations on December 18, 2009. These quasars, with $i$ magnitudes from 16.44 to 19.34 and redshifts from 0.898 to 2.773, were not identified in the SDSS spectroscopic survey, though six of them with redshifts less than 2.5 were selected as quasar targets in SDSS. Except one source without near-IR $Y$-band data, seven of these eight new quasars meet a newly proposed quasar selection criterion involving both near-IR and optical colors. Two of them were found in the 'redshift desert' for quasars ($z$ from 2.2 to 3), indicating that the new criterion is efficient for recovering the missing quasars with similar optical colors as stars. Although LAMOST met some problems during the commissioning observations, we were still able to identify other 38 known SDSS quasars in this field, with $i$ magnitudes from 16.24 to 19.10 and redshifts from 0.297 to 4.512. Our identifications imply that a substantial fraction of quasars may be missing in the previous quasar surveys. The implication of our results to the future LAMOST quasar survey is discussed.

The redshift range from 2.2 to 3, is known as the 'redshift desert' of quasars because quasars with redshift in this range have similar optical colors as normal stars and are thus difficult to be found in optical sky surveys. A quasar candidate, SDSS J085543.40-001517.7, which was selected by a recently proposed criterion involving near-IR $Y-K$ and optical $g-z$ colors, was identified spectroscopically as a new quasar with redshift of 2.427 by the LAMOST commissioning observation in December 2009 and confirmed by the observation made with the NAOC/Xinglong 2.16m telescope in March 2010. This quasar was not targeted in the SDSS spectroscopic survey because it locates in the stellar locus of the optical color-color diagrams, while it is clearly separated from stars in the $Y-K$ vs. $g-z$ diagram. Comparing with other SDSS quasars we found this new quasar with $i$ magnitude of 16.44 is apparently the brightest one in the redshift range from 2.3 to 2.7. From the spectral properties we derived its central black hole mass as $(1.4\sim3.9) \times 10^{10} M_\odot$ and the bolometric luminosity as $3.7\times 10^{48}$ \ergs, which indicates that this new quasar is intrinsically very bright and belongs to the most luminous quasars in the universe. Our identification supports that quasars in the redshift desert can be found by the quasar selection criterion involving the near-IR colors. More missing quasars are expected to be recovered by the future LAMOST spectroscopic surveys, which is important to the study of the cosmological evolution of quasars at redshift higher than 2.2.

We study diffusion of Cosmic Rays (CRs) in turbulent magnetic fields using test particle simulations. Electromagnetic fields are produced in direct numerical MHD simulations of turbulence and used as an input for particle tracing, particle feedback on turbulence being ignored. Statistical transport coefficients from the test particle runs are compared with earlier analytical predictions. We find qualitative correspondence between them in various aspects of CR diffusion. In the incompressible case, that we consider in this paper, the dominant scattering mechanism occurs to be the non-resonant mirror interactions with the slow-mode perturbations. Perpendicular transport roughly agrees with being produced by magnetic field wandering.

Dust grains can be efficiently accelerated and shattered in warm ionized medium (WIM) because of the turbulent motion. This effect is enhanced in starburst galaxies, where gas is ionized and turbulence is sustained by massive stars. Moreover, dust production by Type II supernovae (SNe II) can be efficient in starburst galaxies. In this paper, we examine the effect of shattering in WIM on the dust grains produced by SNe II. We find that although the grains ejected from SNe II are expected to be biased to large sizes ($a\ga 0.1 \micron$, where $a$ is the grain radius) because of the shock destruction in supernova remnants, the shattering in WIM is efficient enough in $\sim 5$ Myr to produce small grains if the metallicity is nearly solar or more. The production of small grains by shattering steepens the extinction curve. Thus, steepening of extinction curves by shattering should always be taken into account for the system where the metallicity is solar and the starburst age is typically larger than 5 Myr. These conditions may be satisfied not only in nearby starbursts but also in high redshift ($z>5$) quasars.

A broad continuum excess in the near-infrared, peaking in the rest-frame at 2-5 micron, is detected in a spectroscopic sample of 88 galaxies at 0.5<z<2.0 taken from the Gemini Deep Deep Survey. Line emission from polycyclic aromatic hydrocarbons (PAHs) at 3.3 micron alone cannot explain the excess, which can be fit by a spectral component consisting of a template of PAH emission lines superposed on a modified blackbody of temperature T~850 K. The luminosity of this near-infrared excess emission at 3 micron is found to be correlated with the star formation rate of the galaxy. The origin of the near-infrared excess is explored by examining similar excesses observed locally in massive star forming regions, reflection and planetary nebulae, post-asymptotic giant branch stars and in the galactic cirrus. We also consider the potential contribution from dust heated around low-luminosity active galactic nuclei. We conclude that the most likely explanation for the 2-5 micron excess is the contribution from circumstellar disks around massive young stellar objects seen in the integrated light of high-redshift galaxies. Assuming circumstellar disks extend down to lower masses, as they do in our own Galaxy, the excess emission presents us with an exciting opportunity to measure the formation rate of planetary systems at cosmic epochs before our own solar system formed.

We discuss a new type of dust acceleration mechanism that acts in a turbulent magnetized medium. The magnetohydrodynamic (MHD) turbulence can accelerate grains through resonant as well as nonresonant interactions. We show that the magnetic compression provides higher velocities for super-Alfvenic turbulence and can accelerate an extended range of grains in warm media compared to gyroresonance. While fast modes dominate the acceleration for the large grains, slow modes can be important for sub-micron grains. We provide comprehensive discussion of all the possible grain acceleration mechanisms in interstellar medium. We show that supersonic velocities are attainable for Galactic dust grains. We discuss the consequence of the acceleration. The implications for extinction curve, grain alignment, chemical abundance, etc, are provided.

We present the results of NICMOS imaging of a sample of 16 high mass passively evolving galaxies with 1.3<z<2, taken primarily from the Gemini Deep Deep Survey. Around 80% of galaxies in our sample have spectra dominated by stars with ages >1 Gyr. Our rest-frame R-band images show that most of these objects have compact regular morphologies which follow the classical R^1/4 law. These galaxies scatter along a tight sequence in the Kormendy relation. Around one-third of the massive red objects are extraordinarily compact, with effective radii under one kiloparsec. Our NICMOS observations allow the detection of such systems more robustly than is possible with optical (rest-frame UV) data, and while similar systems have been seen at z>2, this is the first time such systems have been detected in a rest-frame optical survey at 1.3<z<2. We refer to these compact galaxies as "red nuggets". Similarly compact massive galaxies are completely absent in the nearby Universe. We introduce a new "stellar mass Kormendy relation" (stellar mass density vs size) which isolates the effects of size evolution from those of luminosity and color evolution. The 1.1 < z < 2 passive galaxies have mass densities that are an order of magnitude larger then early type galaxies today and are comparable to the compact distant red galaxies at 2 < z < 3. We briefly consider mechanisms for size evolution in contemporary models focusing on equal-mass mergers and adiabatic expansion driven by stellar mass loss. Neither of these mechanisms appears able to transform the high-redshift Kormendy relation into its local counterpart. <ABRIDGED>

We investigate shattering and coagulation of dust grains in turbulent interstellar medium (ISM). The typical velocity of dust grain as a function of grain size has been calculated for various ISM phases based on a theory of grain dynamics in compressible magnetohydrodynamic turbulence. In this paper, we develop a scheme of grain shattering and coagulation and apply it to turbulent ISM by using the grain velocities predicted by the above turbulence theory. Since large grains tend to acquire large velocity dispersions as shown by earlier studies, large grains tend to be shattered. Large shattering effects are indeed seen in warm ionized medium (WIM) within a few Myr for grains with radius $a\ga 10^{-6}$ cm. We also show that shattering in warm neutral medium (WNM) can limit the largest grain size in ISM ($a\sim 2\times 10^{-5} \mathrm{cm}$). On the other hand, coagulation tends to modify small grains since it only occurs when the grain velocity is small enough. Coagulation significantly modifies the grain size distribution in dense clouds (DC), where a large fraction of the grains with $a<10^{-6}$ cm coagulate in 10 Myr. In fact, the correlation among $R_V$, the carbon bump strength, and the ultraviolet slope in the observed Milky Way extinction curves can be explained by the coagulation in DC. It is possible that the grain size distribution in the Milky Way is determined by a combination of all the above effects of shattering and coagulation. Considering that shattering and coagulation in turbulence are effective if dust-to-gas ratio is typically more than $\sim 1/10$ of the Galactic value, the regulation mechanism of grain size distribution should be different between metal-poor and metal-rich environments.

We calculate the major pair fraction and derive the major merger fraction and rate for 82 massive ($M_{*}>10^{11}M_{\odot}$) galaxies at $1.7 < z < 3.0$ utilising deep HST NICMOS data taken in the GOODS North and South fields. For the first time, our NICMOS data provides imaging with sufficient angular resolution and depth to collate a sufficiently large sample of massive galaxies at z $>$ 1.5 to reliably measure their pair fraction history. We find strong evidence that the pair fraction of massive galaxies evolves with redshift. We calculate a pair fraction of $f_{m}$ = 0.29 +/- 0.06 for our whole sample at $1.7 < z < 3.0$. Specifically, we fit a power law function of the form $f_{m}=f_{0}(1+z)^{m}$ to a combined sample of low redshift data from Conselice et al. (2007) and recently acquired high redshift data from the GOODS NICMOS Survey. We find a best fit to the free parameters of $f_{0}$ = 0.008 +/- 0.003 and $m$ = 3.0 +/- 0.4. We go on to fit a theoretically motivated Press-Schechter curve to this data. This Press-Schechter fit, and the data, show no sign of levelling off or turning over, implying that the merger fraction of massive galaxies continues to rise with redshift out to z $\sim$ 3. Since previous work has established that the merger fraction for lower mass galaxies turns over at z $\sim$ 1.5 - 2.0, this is evidence that higher mass galaxies experience more mergers earlier than their lower mass counterparts, i.e. a galaxy assembly downsizing. Finally, we calculate a merger rate at z = 2.6 of $\Re$ $<$ 5 $\times$ 10$^{5}$ Gpc$^{-3}$ Gyr$^{-1}$, which experiences no significant change to $\Re$ $<$ 1.2 $\times$ 10$^{5}$ Gpc$^{-3}$ Gyr$^{-1}$ at z = 0.5. This corresponds to an average $M_{*}>10^{11}M_{\odot}$ galaxy experiencing 1.7 +/- 0.5 mergers between z = 3 and z = 0.

Astrophysical fluids, including interstellar and interplanetary medium, are magnetized and turbulent. Their appearance, evolution, and overall properties are determined by the magnetic turbulence that stirs it. We argue that examining magnetic turbulence at a fundamental level is vital to understanding many processes. A point that frequently escapes the attention of researchers is that magnetic turbulence cannot be confidently understood only using "brute force" numerical approaches. In this review we illustrate this point on a number of examples, including intermittent heating of plasma by turbulence, interactions of turbulence with cosmic rays and effects of turbulence on the rate of magnetic reconnection. We show that the properties of magnetic turbulence may vary considerably in various environments, e.g. imbalanced turbulence in solar wind differs from balanced turbulence and both of these differ from turbulence in partially ionized gas. Appealing for the necessity of more observational data on magnetic fields, we discuss a possibility of studying interplanetary turbulence using alignment of Sodium atoms in the tail of comets.

We measure the sizes of 82 massive (M>10^11 M_sun) galaxies at 1.7<z<3 utilizing deep HST NICMOS data taken in the GOODS North and South fields. Our sample is almost an order of magnitude larger than previous studies at these redshifts, providing the first statistical study of massive galaxy sizes at z>2, confirming the extreme compactness of these galaxies. We split our sample into disk-like (n<2) and spheroid-like (n>2) galaxies based on their Sersic indices, and find that at a given stellar mass disk-like galaxies at z~2.3 are a factor of 2.6+/-0.3 smaller than present day equal mass systems, and spheroid-like galaxies at the same redshifts are 4.3+/-0.7 smaller than comparatively massive elliptical galaxies today. At z>2 our results are compatible with both a leveling off, or a mild evolution in size. Furthermore, the high density (~2x10^10 M_sun kpc^-3) of massive galaxies at these redshifts, which are similar to present day globular clusters, possibly makes any further evolution in sizes beyond z=3 unlikely.

We discuss a new technique of studying magnetic fields in diffuse astrophysical media, e.g. interstellar and intergalactic gas/plasma. This technique is based on the angular momentum alignment of atoms and ions in their ground or metastable states. As the life-time of atoms in such states is long, the alignment induced by anisotropic radiation is susceptible to weak magnetic fields ($1{\rm G}\gtrsim B\gtrsim0.1\mu$G). The alignment reveals itself in terms of the polarization of the absorbed and emitted light. A variety of atoms with fine or hyperfine splitting of the ground or metastable states exhibit the alignment and the resulting polarization degree in some cases exceeds 20%. We show that in the case of absorption the polarization direction is either parallel or perpendicular to magnetic field, while more complex dependencies emerge for the case of emission of aligned atoms. We show that the corresponding studies of magnetic fields can be performed with optical and UV polarimetry. A unique feature of these studies is that they can reveal the 3D orientation of magnetic field. In addition, we point out that the polarization of the radiation arising from the transitions between fine and hyperfine states of the ground level can provide yet another promising diagnostics of magnetic fields, including the magnetic fields in the Early Universe. We mention several cases of interplanetary, circumstellar and interstellar magnetic fields for which the studies of magnetic fields using ground state atomic alignment effect are promising.

We address the problem of particle acceleration in solar flares by fast modes which may be excited during the reconnection and undergo cascade and are subjected to damping. We extend the calculations beyond quasilinear approximation and compare the acceleration and scattering by transit time damping and gyroresonance interactions. We find that the acceleration is dominated by the so called transit time damping mechanism. We estimate the total energy transferred into particles, and show that our approach provides sufficiently accurate results We compare this rate with energy loss rate. Scattering by fast modes appears to be sufficient to prevent the protons from escaping the system during the acceleration. Confinement of electrons, on the other hand, requires the existence of plasma waves. Electrons can be accelerated to GeV energies through the process described here for solar flare conditions.

We continue our studies of atomic alignment in diffuse media, in particularly, in interstellar and circumstellar media, with the goal of developing new diagnostics of magnetic fields in these environments. We understand atomic alignment as alignment of atoms or ions in their ground state. Such atoms are sensitive to weak magnetic fields. In particular, we provide predictions of the polarization that arises from astrophysically important aligned atoms (ions) with fine structure of the ground level, namely, OI and SII and Ti II. Unlike our earlier papers which dealt with weak fields only, a substantial part of our current paper is devoted to the studies of atomic alignment when magnetic fields get strong enough to affect the emission from the excited level, i.e. with the regime when the magnetic splitting is comparable to the line-width. This is a regime of Hanle effect modified by the atomic alignment. Using an example of emission and absorption lines of SII ion we demonstrate how polarimetric studies can probe magnetic fields in circumstellar regions and accretion disks. In addition, we show that atomic alignment induced by anisotropic radiation can induce substantial variations of magnetic dipole transitions within the ground state, thus affecting abundance studies based on this emission. Moreover, the radio emission is polarized, provides a new way to study magnetic fields, e.g. at the epoch of Universe reionization.

Numerical simulations shed light onto earlier not trackable problem of magnetohydrodynamic (MHD) turbulence. They allowed to test the predictions of different models and choose the correct ones. Inevitably, this progress calls for revisions in the picture of cosmic ray (CR) transport. It also shed light on the problems with the present day numerical modeling of CR. In this paper we focus on the analytical way of describing CR propagation and scattering, which should be used in synergy with the numerical studies. In particular, we use recently established scaling laws for MHD modes to obtain the transport properties for CRs. We include nonlinear effects arising from large scale trapping, to remove the 90 degree divergence. We determine how the efficiency of the scattering and CR mean free path depend on the characteristics of ionized media, e.g. plasma $\beta$, Coulomb collisional mean free path. Implications for particle transport in interstellar medium and solar corona are discussed. We also examine the perpendicular transport of CRs. Perpendicular transport depends on the comparison of parallel mean free path and the injection scale of the turbulence, as well as the Alfv\'enic Mach number. Normal turbulence does not allow subdiffusion unless there are slab waves. The critical scale below which subdiffusion applies is provided. These results can be used to compare with the numerical calculations, provided that these calculations use the structure of magnetic field which is consistent with the numerical studies of MHD turbulence.

Using a set of deep imaging data obtained by the ACS on the HST shortly after its deployment, Yan, Windhorst & Cohen (2003) found a large number of F775W-band dropouts (i-dropouts), which are consistent with being galaxies at z~6. The surface density of i-dropouts thus derived, however, is an order of magnitude higher than those subsequent studies found in other deep ACS fields, including the HUDF. Here we revisit this problem, using both the existing and the new data. We confirm that the large overdensity of i-dropouts does exist in this field, and that their optical-to-IR colors are similar to those in the HUDF. However, we have discovered that the i-dropout overdensity is accompanied with an even larger excess of faint field objects in this region and its vicinity. This large excess of field objects is most likely caused by the tidal ``debris'' and/or halo stars related to an interacting galaxy pair in the Virgo Cluster, M60/NGC4647, which lies several arcminutes away from the region where the excess is found. This excess population is also red in color, and the red wing of its color distribution continuously extends to the regime where the i-dropouts reside. While we still cannot completely rule out the possibility that the overdensity of i-dropouts might be a genuine large-scale structure of galaxies at z~6, we prefer the interpretation that most of them are part of the excess stellar population related to M60/NGC4647. Future spectroscopic work will be needed to identify the nature of this i-dropout overdensity. (Abridged)

We consider the propagation of cosmic rays in turbulent magnetic fields. We use the models of magnetohydrodynamic turbulence that were tested in numerical simulations, in which the turbulence is injected on large scale and cascades to small scales. Our attention is focused on the models of the strong turbulence, but we also briefly discuss the effects that the weak turbulence and the slab Alfv\'enic perturbations can have. The latter are likely to emerge as a result of instabilities with in the cosmic ray fluid itself, e.g., beaming and gyroresonance instabilities of cosmic rays. To describe the interaction of cosmic rays with magnetic perturbations we develop a non-linear formalism that extends the ordinary Quasi-Linear Theory (QLT) that is routinely used for the purpose. This allows us to avoid the usual problem of 90 degree scattering and enable our computation of the mean free path of cosmic rays. We apply the formalism to the cosmic ray propagation in the galactic halo and in the Warm Ionized medium (WIM). In addition, we address the issue of the transport of cosmic rays perpendicular to the mean magnetic field and show that the issue of cosmic ray subdiffusion (i.e., propagation with retracing the trajectories backwards, which slows down the diffusion) is only important for restricted cases when the ambient turbulence is far from what numerical simulations suggest to us. As a result, this work provides formalism that can be applied for calculating cosmic ray propagation in a wide variety of circumstances.

We present the rest-frame ultraviolet through near infrared spectral energy distribution for an interacting Lyman break galaxy at a redshift z=4.42, the highest redshift merging system known with clearly resolved tidal features. The two objects in this system - HDF-G4 and its previously unidentified companion - are both B_{435} band dropouts, have similar V_{606}-i_{775} and i_{775}-z_{850} colors, and are separated by 1", which at z=4.42 corresponds to 7 kpc projected nuclear separation; all indicative of an interacting system. Fits to stellar population models indicate a stellar mass of M_\star = 2.6\times 10^{10} M_\odot, age of \tau_\star = 720 My, and exponential star formation history with an e-folding time \tau_0 = 440 My. Using these derived stellar populations as constraints, we model the HDF-G4 system using hydrodynamical simulations, and find that it will likely evolve into a quasar by z\sim3.5, and a quiescent, compact spheroid by z\sim 2.5 similar to those observed at z > 2. And, the existence of such an object supports galaxy formation models in which major mergers drive the high redshift buildup of spheroids and black holes.

We have used the Hubble Space Telescope's Advanced Camera for Surveys (Ford et al. 2003) to measure the cumulative mass density in morphologically-selected early-type galaxies over the redshift range 0.8 < z < 1.7. Our imaging data set covers four well-separated sight-lines, and is roughly intermediate (in terms of both depth and area) between the GOODS/GEMS imaging data, and the images obtained in the Hubble Deep Field campaigns. Our images contain 144 galaxies with ultra-deep spectroscopy obtained as part of the Gemini Deep Deep Survey. These images have been analyzed using a new purpose-written morphological analysis code which improves the reliability of morphological classifications by adopting a 'quasi-Petrosian' image thresholding technique. We find that at z \~ 1 about 80% of the stars living in the most massive galaxies reside in early-type systems. This fraction is similar to that seen in the local Universe. However, we detect very rapid evolution in this fraction over the range 0.8 < z < 1.7, suggesting that over this redshift range the strong morphology-mass relationship seen in the nearby Universe is beginning to fall into place. By comparing our images to published spectroscopic classifications, we show that little ambiguity exists in connecting spectral classes to morphological classes for spectroscopically quiescent systems. However, the mass density function of early-type galaxies is evolving more rapidly than that of spectroscopically quiescent systems, which we take as further evidence that we are witnessing the formation of massive early-type galaxies over the 0.8 < z < 1.7 redshift range.

We present a catalog of spectro-photometric redshifts for 1308 galaxies from the GRism ACS Program for Extragalactic Science (GRAPES) observations with the Hubble Space Telescope. These low-resolution spectra between 6000 A and 9500 A are supplemented with U, J, H, and Ks from various facilities, resulting in redshifts computed with ~40 spectral bins per galaxy. For 81 galaxies between 0.5<z<1.5 with spectroscopic redshifts, the standard deviation in the fractional error in (1+z) is 0.046. With this catalog, we compute the B-band luminosity function in this redshift range from 72 galaxies. Owing to the depth of the GRAPES survey, we are able to accurately constrain the faint-end slope by going to M_B~-18 mag at 0.8<z<1.2, nearly two magnitudes fainter than previous studies. The faint-end slope is alpha=-1.32+-0.07. When compared to numerous published values at various redshifts, we find strong evidence for a steepening of the faint-end slope with redshift which is expected in the hierarchical formation scenario of galaxies.

We present redshifts for 115 emission line objects in the Hubble Ultra Deep Field (HUDF) identified through the GRism ACS Program for Extragalactic Science (GRAPES) project using the slitless grism spectroscopy mode of the ACS Camera on the Hubble Space Telescope (HST). The sample was selected by an emission line search on all extracted 1-dimensional GRAPES spectra. We identify the emission lines using line wavelength ratios where multiple lines are detected in the grism wavelength range (5800A < lambda < 9600A), and using photometric redshift information where multiple lines are unavailable. We then derive redshifts using the identified lines. Our redshifts are accurate to delta(z) = 0.009, based on both statistical uncertainty estimates and comparison with published ground-based spectra. Over 40% of our sample is fainter than typical magnitude limits for ground-based spectroscopy (with i_{AB}>25 mag). Such emission lines would likely remain undiscovered without our deep survey. The emission line objects fall into 3 categories: 1) Most are low to moderate redshift galaxies (0 < z < 2), including many actively star forming galaxies with strong HII regions; 2) 9 are high redshift (4 < z < 7) Lyman-alpha emitters; and 3) at least 3 are candidate AGNs.

We describe a compact cluster of massive red galaxies at z=1.51 discovered in one of the Gemini Deep Deep Survey (GDDS) fields. Deep imaging with the Near Infrared Camera and Multi Object Spectrometer (NICMOS) on the Hubble Space Telescope reveals a high density of galaxies with red optical to near-IR colors surrounding a galaxy with a spectroscopic redshift of 1.51. Mid-IR imaging with Infrared Array Camera (IRAC) on the Spitzer Space telescope shows that these galaxies have spectral energy distributions that peak between 3.6 and 4.5 microns. Fits to 12-band photometry reveal 12 or more galaxies with spectral shapes consistent with z = 1.51. Most are within ~170 co-moving kpc of the GDDS galaxy. Deep F814W images with the Advanced Camera for Surveys (ACS) on HST reveal that these galaxies are a mix of early-type galaxies, disk galaxies and close pairs. The total stellar mass enclosed within a sphere of 170 kpc in radius is > 8E+11 solar masses. The colors of the most massive galaxies are close to those expected from passive evolution of simple stellar populations (SSP) formed at much higher redshifts. We suggest that several of these galaxies will merge to form a single, very massive galaxy by the present day. This system may represent an example of a short-lived dense group or cluster core typical of the progenitors of massive clusters in the present day and suggests the red sequence was in place in over-dense regions at early times.

We show that atomic alignment presents a reliable way to study topology of astrophysical magnetic fields. The effect of atomic alignment arises from modulation of the relative population of the sublevels of atomic ground state pumped by anisotropic radiation flux. As such aligned atoms precess in the external magnetic field and this affects the properties of the polarized radiation arising from both scattering and absorption by the atoms. As the result the polarizations of emission and absorption lines depend on the 3D geometry of the magnetic field as well as the direction and anisotropy of incident radiation. We consider a subset of astrophysically important atoms with hyperfine structure. For emission lines we obtain the dependencies of the direction of linear polarization on the directions of magnetic field and the incident pumping radiation. For absorption lines we establish when the polarization is perpendicular and parallel to magnetic field. For both emission and absorption lines we find the dependence on the degree of polarization on the 3D geometry of magnetic field. We claim that atomic alignment provides a unique tool to study magnetic fields in circumstellar regions, AGN, interplanetary and interstellar medium. This tool allows studying of 3D topology of magnetic fields and establish other important astrophysical parameters. We consider polarization arising from both atoms in the steady state and also as they undergo individual scattering of photons. We exemplify the utility of atomic alignment for studies of astrophysical magnetic fields by considering a case of Na alignment in a comet wake.

The relative population of the fine structure sublevels of an atom's ground state is affected by radiative transitions induced by an anisotropic radiation flux. This causes the alignment of atomic angular momentum. In terms of observational consequences for the interstellar and intergalactic medium, this results in the polarization of the absorption lines. In the paper we consider the conditions necessary for this effect and provide calculations of polarization from a few astrophysically important atoms and ions with multiple upper and lower levels for an arbitrary orientation of magnetic fields to the a) source of optical pumping, b) direction of observation, c) absorbed source. We also consider an astrophysically important ``degenerate'' case when the source of optical pumping coincides with the source of the absorbed radiation. We present analytical expressions that relate the degree of linear polarization and the intensity of absorption to the 3D orientation of the magnetic field with respect to the pumping source, the source of the absorbed radiation, and the direction of observations. We discuss how all these parameters can be determined via simultaneous observations of several absorption lines and suggest graphical means that are helpful in practical data interpretation. We prove that studies of absorption line polarization provide a unique tool to study 3D magnetic field topology in various astrophysical conditions.

Using the deep Spitzer IRAC observations of the Great Observatories Origins Deep Survey (GOODS), we study the stellar masses and star formation histories of galaxies at z ~ 6. Our study is based on the i-band dropout sample selected from the GOODS southern and northern fields (~ 330 arcmin^2 in total), several of which already have spectroscopic confirmations. In total, we derive stellar masses for 53 i-band dropouts that have robust IRAC detections. These galaxies have typical stellar masses of ~ 1e10 M_sun and typical ages of a couple of hundred million years, consistent with earlier results based on a smaller sample of z ~ 6 galaxies in the HUDF. We also study 79 i-band dropouts that are invisible in the IRAC data and find that they are typically less massive by a factor of ten. Based on our mass estimates, we derive a lower limit to the global stellar mass density at z ~ 6. Considering the range of systematic uncertainties in the derived stellar masses, this lower limit is 1.1--6.7e6 M_sun/Mpc^3 (co-moving), which is 0.2--1.1% of the present-day value. The prospect of detecting the progenitors of the most massive galaxies at yet higher redshifts is explored: a deep, wide-field near-IR survey using our current technology could possibly result in positive detections at z > 7. We also investigate the implication of our results for reionzation, and find that the progenitors of the galaxies comparable to those in our sample, even in the most optimized (probably unrealistic) scenario, cannot sustain the reionization for a period longer than ~ 2 million years. Thus most of the photons required for reionization must have been provided by other sources, such as the progenitors of the dwarf galaxies that are far below our current detection capability. (Abridged)

We propose collisionless damping of fast MHD waves as an important mechanism for the heating and acceleration of winds from rotating stars. Stellar rotation causes magnetic field lines anchored at the surface to form a spiral pattern and magneto-rotational winds can be driven. If the structure is a magnetically dominated, fast MHD waves generated at the surface can propagate almost radially outward and cross the field lines. The propagating waves undergo collisionless damping owing to interactions with particles surfing on magnetic mirrors that are formed by the waves themselves. The damping is especially effective where the angle between the wave propagation and the field lines becomes moderately large ($\sim 20$ to $80^{\circ}$). The angle tends naturally to increase into this range because the field in magneto-rotational winds develops an increasingly large azimuthal component. The dissipation of the wave energy produces heating and acceleration of the outflow. We show using specified wind structures that this damping process can be important in both solar-type stars and massive stars that have moderately large rotation rates. This mechanism can play a role in coronae of young solar-type stars which are rapidly rotating and show X-ray luminosities much larger than the sun. The mechanism could also be important for producing the extended X-ray emitting regions inferred to exist in massive stars of spectral type middle B and later.

We investigate the properties of massive galaxies at z=1-3.5 using HST observations, ground-based near-IR imaging, and Spitzer Space Telescope observations at 3-24 micron. We identify 153 distant red galaxies (DRGs) with J-K > 2.3 mag (Vega) in the southern GOODS field. This sample is approximately complete in stellar mass for passively evolving galaxies above 10^11 solar masses and z < 3. The galaxies identified by this selection are roughly split between objects whose optical and near-IR rest-frame light is dominated by evolved stars combined with ongoing star formation, and galaxies whose light is dominated by heavily reddened starbursts. Very few of the galaxies (< 10%) have no indication of current star formation. Using SFR estimates that include the reradiated IR emission, the DRGs at z=1.5-3 with stellar masses > 10^11 solar masses have specific SFRs (SFRs per unit stellar mass) ranging from 0.2 to 10 Gyr^-1, with a mean value of ~2.4 Gyr^-1. The DRGs with stellar masses > 10^11 solar masses and 1.5 < z < 3 have integrated specific SFRs greater the global value over all galaxies. In contrast, we find that galaxies at z = 0.3-0.75 with these stellar masses have integrated specific SFRs less than the global value, and more than an order of magnitude lower than that for massive DRGs at z = 1.5-3. At z < 1, lower-mass galaxies dominate the overall cosmic mass assembly. This suggests that the bulk of star formation in massive galaxies occurs at early cosmic epochs and is largely complete by z~1.5. [Abridged]

(Abridged) We describe the cascade of plasma waves or turbulence injected, presumably by reconnection, at scales comparable to the size of a solar flare loop to scales comparable to particle gyroradii, and evaluate their damping by various mechanisms. We show that the classical viscous damping is unimportant for magnetically dominated or low beta plasmas and the primary damping mechanism is the collisionless damping by the background particles. We show that the damping rate is proportional to the total random momentum density of the particles. For solar flare conditions this means that in most flares, except the very large ones, the damping is dominated by thermal background electrons. For large flares one requires acceleration of essentially all background electrons into a nonthermal distribution so that the accelerated electrons can be important in the damping of the waves. In general, damping by thermal or nonthermal protons is negligible compared to that of electrons except for quasi-perpendicular propagating waves or for rare proton dominated flares with strong nuclear gamma-ray line emission. Using the rate for damping we determine the critical scale below which the damping becomes important and the spectrum of the turbulence steepens. This critical scale, however, has strong dependence on the angle of propagation with respect to the magnetic field direction. The waves can cascade down to very small scales, such as the gyroradii of the particles at small angles (quasi-parallel propagation) and possibly near 90 degree (quasi-perpendicular propagation) giving rise to a highly anisotropic spectral distribution.

We use data from the first epoch of observations with the IRAC/Spitzer for the GOODS to detect and study a collection of LBGs at z ~ 6 to 5 in the HUDF, six of which have spectroscopic confirmation. At these redshifts, IRAC samples rest-frame optical light in the range 0.5 to 0.8 um, where the effects of dust extinction are smaller and the sensitivity to light from evolved stars is greater than at shorter, rest-frame UV observable from the ground or with the HST. As such, it provides useful constraints on the ages and masses of these galaxies' stellar populations. We find that the SEDs for many of these galaxies are best fitted by models of stellar populations with masses of a few x 1e10 M_sun, and with ages of a few hundred Myrs, values quite similar to those derived for typical LBGs at z ~ 3. When the universe was only 1 Gyr old, some galaxies had already formed a mass of stars approaching that of the present-day Milky Way, and that they started forming those stars at z > 7, and in some cases much earlier. We find that the lower limits to the space density for galaxies in this mass range are consistent with predictions from recent hydrodynamic simulations of structure formation in a LCDM universe. All objects in our samples are consistent with having solar metallicity, suggesting that they might have already been significantly polluted by metals. The values for dust reddening derived from the model fitting are low or zero, and we find that some of the galaxies have rest-frame UV colors that are even bluer than those predicted by the stellar population models to which we compare them. These colors might be attributed to the presence of very massive stars (> 100 M_sun), or by weaker intergalactic HI absorption than what is commonly assumed.

We summarize the design and performance of the James Webb Space Telescope that is to be launched to an L2 orbit in 2011, and how it is designed, in particular, to study the epochs of First Light, Reionization and Galaxy Assembly.

Most astrophysical fluids are turbulent and magnetized. Fluctuations of polarization provide a promising way to study astrophysical magnetic turbulence. We discuss polarization that arises from grains and atoms aligned in respect to magnetic field, describe the processes of alignment, explain when the alignment reflects the direction of magnetic field. We show that temporal fluctuations of interplanetary magnetic field can be studied using aligned sodium atoms.

Recent advances in understanding of magnetohydrodynamic (MHD) turbulence call for revisions in the picture of cosmic ray transport. In this paper we use recently obtained scaling laws for MHD modes to obtain the scattering frequency for cosmic rays. We account for the turbulence cutoff arising from both collisional and collisionless damping. We obtain the scattering rate and show that fast modes provide the dominant contribution to cosmic ray scattering for the typical interstellar conditions in spite of the fact that fast modes are subjected to damping. We determine how the efficiency of the scattering depends on the characteristics of ionized media, e.g. plasma $\beta$. We show that streaming instability is suppressed by the ambient MHD turbulence.

We present grism spectra taken with the Advanced Camera for Surveys to identify 29 red sources with (i'-z') > 0.9 in the Hubble Ultra Deep Field (HUDF). Of these 23 are found to be galaxies at redshifts between z=5.4 and 6.7, identified by the break at 1216 angstrom due to IGM absorption; two are late type dwarf stars with red colors; and four are galaxies with colors and spectral shape similar to dust reddened or old galaxies at redshifts z~1-2. This constitutes the largest uniform, flux-limited sample of spectroscopically confirmed galaxies at such faint fluxes z'<27.5. Many are also among the most distant spectroscopically confirmed galaxies (at redshifts up to z=6.7). We find a significant overdensity of galaxies at redshifts z=5.9 +/- 0.2. Nearly two thirds of the galaxies in our sample (15/23) belong to this peak. Taking into account the selection function and the redshift sensitivity of the survey, we get a conservative overdensity of at least a factor of two along the line-of-sight. The galaxies found in this redshift peak are also localized in the plane of the sky in a non-random manner, occupying about half of the ACS chip. Thus the volume overdensity is a factor of four. The star-formation rate derived from detected sources in this overdense region is sufficient to reionize the local IGM.

We discuss the damping of strong magnetic turbulence downstream of the forward shock of young supernova remnants (SNR). We find that strong magnetic fields, that have been produced by the streaming instability in the upstream region of the shock, or by other kinetic instabilities at the shock, will be efficiently reduced, so the region of enhanced magnetic field strength would typically have a thickness of the order l_d=(10^{16}-10^{17}) cm. The non-thermal X-ray filaments observed in young SNR are thus likely limited by the magnetic field and not by the energy losses of the radiating electrons. Consequently the thickness of the filaments would not be a measure of the magnetic field strength and claims of efficient cosmic-ray acceleration on account of a run-away streaming instability appear premature.

Population of levels of the hyperfine and fine split ground state of an atom is affected by radiative transitions induced by anisotropic radiation flux. Such aligned atoms precess in the external magnetic field and this affects properties of polarized radiation arising from both scattering and absorption by atoms. As the result the degree of light polarization depends on the direction of the magnetic field. This provides a perspective tool for studies of astrophysical magnetic fields using optical and UV polarimetry. We provide calculations for several atoms and ions that can be used to study magnetic fields in interplanetary medium, interstellar medium, circumstellar regions and quasars.

We have discovered an extended Lyman alpha plume associated with a compact source at redshift 5.4 in slitless spectroscopic data from the Grism ACS Program for Extragalactic Science (GRAPES) project. The spatial extent of the emission is about 6 x 1.5 kpc (1 x 0.25 arcsec). Combining our grism data and the broadband images from the Hubble UltraDeep Field (UDF) images, we find a Lyman alpha line flux of 2e-17 erg/cm2/s and surface brightness 7e-17 erg/cm2/s/arcsec2. The UDF images show diffuse continuum emission associated with the Lyman alpha plume (hereafter UDF 5225), with three embedded knots. The morphology of UDF 5225 is highly suggestive of a galaxy in assembly. It is moreover possible that the prominent Lyman alpha emission from this object is due to an active nucleus, and that we are seeing the simultaneous growth through accretion of a galaxy and its central black hole. Followup observations at higher spectral resolution could test this hypothesis.

The quasilinear theory for cosmic ray propagation is a well known and widely accepted theory. In this paper, we discuss the different contributions to the pitch-angle Fokker-Planck coefficient from large and small scales for slab geometry using the damping model of dynamical turbulence. These examinations will give us a hint on the limitation range where quasilinear approximation is a good approximation.

We calculate the relative grain-grain motions arising from interstellar magnetohydrodynamic (MHD) turbulence. The MHD turbulence includes both fluid motions and magnetic fluctuations. While the fluid motions accelerate grains through hydro-drag, the electromagnetic fluctuations accelerate grains through resonant interactions. We consider both incompressive (Alfv\'{e}n) and compressive (fast and slow) MHD modes and use descriptions of MHD turbulence obtained in Cho & Lazarian (2002). Calculations of grain relative motion are made for realistic grain charging and interstellar turbulence that is consistent with the velocity dispersions observed in diffuse gas, including cutoff of the turbulence from various damping processes. We show that fast modes dominate grain acceleration, and can drive grains to supersonic velocities. Grains are also scattered by gyroresonance interactions, but the scattering is less important than acceleration for grains moving with sub-Alfv\'{e}nic velocities. Since the grains are preferentially accelerated with large pitch angles, the supersonic grains will be aligned with long axes perpendicular to the magnetic field. We compare grain velocities arising from MHD turbulence with those arising from photoelectric emission, radiation pressure and H$_{2}$ thrust. We show that for typical interstellar conditions turbulence should prevent these mechanisms from segregating small and large grains. Finally, gyroresonant acceleration is bound to preaccelerate grains that are further accelerated in shocks. Grain-grain collisions in the shock may then contribute to the overabundance of refractory elements in the composition of galactic cosmic rays.

Recent advances in understanding of magnetohydrodynamic (MHD) turbulence call for revisions in the picture of cosmic ray transport. In this paper we use recently obtained scaling laws for MHD modes to obtain the scattering frequency for cosmic rays. Using quasilinear theory we calculate gyroresonance with MHD modes (Alfv\'{e}nic, slow and fast) and transit-time damping (TTD) by fast modes. We provide calculations of cosmic ray scattering for various phases of interstellar medium with realistic interstellar turbulence driving that is consistent with the velocity dispersions observed in diffuse gas. We account for the turbulence cutoff arising from both collisional and collisionless damping. We obtain analytical expressions for diffusion coefficients that enter Fokker-Planck equation describing cosmic ray evolution. We obtain the scattering rate and show that fast modes provide the dominant contribution to cosmic ray scattering for the typical interstellar conditions in spite of the fact that fast modes are subjected to damping. We determine how the efficiency of the scattering depends on the characteristics of ionized media, e.g. plasma $\beta$. We calculate the range of energies for which the streaming instability is suppressed by the ambient MHD turbulence.

Using early data from the Infrared Array Camera (IRAC) on the Spitzer Space Telescope, taken for the Great Observatories Origins Deep Survey (GOODS), we identify and study objects that are well-detected at 3.6um, but are very faint (and in some cases, invisible) in the Hubble Ultra Deep Field (HUDF) ACS and NICMOS images and in very deep VLT Ks-band imaging. We select a sample of 17 objects with f_nu(3.6um)/f_nu(z_{850})>20. The analysis of their spectral energy-distributions (SEDs) from 0.4 to 8.0um shows that the majority of these objects cannot be satisfactorily explained without a well-evolved stellar population. We find that most of them can be well fitted by a simple two-component model, where the primary component represents a massive, old population that dominates the strong IR emission, while the secondary component represents a low amplitude, on-going star-formation process that accounts for the weak optical fluxes. Their estimated photometric redshifts (z_p) range from 1.6 to 2.9, with the median at z_p=2.4. For the simple star formation histories considered here, their corresponding stellar masses range from 0.1 to 1.6e11 M_sun for a Chabrier Initial Mass Function (IMF). Their median rest-frame Ks-band absolute magnitude is -22.9 mag in AB system, or 1.5xL*(K) for present--day elliptical galaxies. In the scenario of pure luminosity evolution, such objects may be direct progenitors for at least 14 to 51% of the local population of early--type galaxies. Due to the small cosmic volume of the HUDF, however, this simple estimate could be affected by other effects such as cosmic variance and the strong clustering of massive galaxies. A full analysis of the entire GOODS area is now underway to assess such effects.

We report results from our $z\simeq 5.5$--7 galaxy search in the HST Ultra Deep Field (UDF). Using the 400-orbit of ACS data, we found 108 plausible $5.5\leq z\leq 6.5$ (or $z\simeq 6$ for short) candidates to $m_{AB}(z_{850})=30.0$ mag. The contamination to the sample, either due to image artifacts or known types of astronomical objects, is likely negligible. The inferred surface densities of $z\simeq 6$ galaxies are consistent with our earlier predictions from $m_{AB}(z_{850})=26.5$ to 28.5 mag. After correcting for detection incompleteness, the counts of $z\simeq 6$ candidates to $m_{AB}(z_{850})=29.2$ mag suggests that the faint-end slope of the galaxy luminosity function (LF) at this redshift is likely between $\alpha=-1.8$ and -1.9, which is sufficient to account for the entire Lyman photon budget necessary to complete the reionization of the universe at $z\simeq 6$. We also searched for $z\simeq 6.5$--7 candidates using the UDF NICMOS data, and have found four candidates to $J_{110}=27.2$ mag. However, the infrared colors of three candidates cannot be easily explained by galaxies in this redshift range. We tentatively derive an upper limit to the cumulative surface density of galaxies at $z\simeq 7$ of 0.36 per arcmin$^2$ to $J_{110}=26.6$ mag, which suggest a noticeable drop in the LF amplitude from $z\simeq 6$ to $z\simeq 7$.

The Spitzer Space Telescope has undertaken the deepest ever observations of the 24 micron sky in the ELAIS-N1 field as part of GOODS Science Verification observations. We present the shape of the 24 micron source counts in the flux range 20-1000 microJy, discuss the redshift distribution and nature of these sources with particular emphasis on their near-infrared counterparts.

We present deep unbiased spectroscopy of the Hubble Ultra Deep Field (UDF) carried out using the slitless grism spectroscopy mode of the Advance Camera for Surveys on the Hubble Space Telescope (HST). The GRIsm ACS Program for Extragalactic Science (GRAPES) achieves continuum detection as faint as $z_{AB}=27.2$ using 40 orbits ($9.2 \times 10^4$ seconds) on HST. The data were taken at four orientation angles to correct for the overlap of spectra. GRAPES data provide a unique, uninterrupted, low resolution (R=100) spectral coverage for $5500\AA < \lambda < 10500\AA$, and allow us to detect high redshift galaxies at $4 < z < 7$ whether they have $\lya$ lines or just show the Lyman Break, as well as find low luminosity AGNs in an unbiased fashion. This paper describes in detail the observations and the data reduction, and examines the quality of the extracted spectra. Subsequent papers will deal with the analysis of the data. The extracted and calibrated GRAPES spectra will be available from MAST at STScI.

In this paper, we address which sources contributed most of the reionizing photons. Our argument assumes that the reionization ended around z ~ 6 and that it was a relatively quick process, i.e., that there was a non-negligible fraction of neutral hydrogen in the Universe at somewhat earlier epochs. Starting from our earlier estimate of the luminosity function (LF) of galaxies at z ~ 6, we quantitatively show that the major sources of reionization are most likely galaxies with L < L*. Our approach allows us to put stronger constraints to the LF of galaxies at z ~ 6. To have the Universe completely ionized at this redshift, the faint-end slope of the LF should be steeper than $\alpha=-1.6$, which is the value measured at lower redshifts (z ~ 3), unless either the normalization (Phi*) of the LF or the clumping factor of the ionized hydrogen has been significantly underestimated. If Phi* is actually lower than what we assumed by a factor of two, a steep slope close to $\alpha=-2.0$ is required. Our LF predicts a total of 50 -- 80 z ~ 6 galaxies in the HST Ultra Deep Field (UDF) to a depth of AB=28.4 mag, which can be used to constraint both Phi* and $\alpha$. We conclude that the least luminous galaxies existing at this redshift should reach as low as some critical luminosity in order to accumulate the entire reionizing photon budget. On the other hand, the existence of significant amounts of neutral hydrogen at slightly earlier epochs, e.g. z ~ 7, requires that the least luminous galaxies should not be fainter than another critical value (i.e., the LF should cut-off at this point).

Interstellar dust grains exhibit complex dynamics which is essential for understanding many key interstellar processes that involve dust, including grain alignment, grain growth, grain shattering etc. Grain rotational and translational motions are affected not only by gaseous collisions, but also by interactions with ions, photons, magnetic fields etc. Some of those interactions, e.g. interactions of ions with the dipole electric moment of dust grains, require the quantum nature of the process to be accounted for. Similarly, coupling of rotational and vibrational degrees of freedom in a grain happens due to relaxation processes, among which the process related to nuclear spins frequently is the dominant one. This coupling modifies substantially both the dynamics of rotational and translational motions by inducing grain flips. The flipping averages our certain systematic torques and forces that act on grains. As the rate of flipping is larger for smaller grains, these grains can be ``thermally trapped'', i.e. rotate at a thermal rate in spite of the presence of systematic torques. Moreover, a subset of small grains with high dipole moments may rotate at subthermal rates due to high damping arising from grain emission in microwave range of frequencies. Translational dynamics of grains is mostly dominated by grain interactions with magnetohydrodynamic turbulence. Efficient turbulent mixing of dust grains limits the degree to which grains of different sizes may be segregated in space.

Recent advances in understanding of magnetohydrodynamic (MHD) turbulence call for substantial revisions in the picture of cosmic ray transport. In this paper we use recently obtained scaling laws for MHD modes to calculate the scattering frequency for cosmic rays in the ISM. We consider gyroresonance with MHD modes (Alfvenic, slow and fast) and transit-time damping (TTD) by fast modes. We provide calculations of cosmic ray scattering for various phases of interstellar medium with realistic interstellar turbulence driving that is consistent with the velocity dispersions observed in diffuse gas. We account for the turbulence cutoff arising from both collisional and collisionless damping. We obtain analytical expressions for diffusion coefficients that enter Fokker-Planck equation describing cosmic ray evolution. We calculate the scattering rate and parallel spatial diffusion coefficients of cosmic rays for both Alfvenic and fast modes. We conclude that fast modes provides the dominant contribution to cosmic ray scattering for the typical interstellar conditions in spite of the fact that fast modes are subjected to damping. We show that the efficiency of the scattering depends on the plasma beta since it determines the damping of the fast modes. We also show that the streaming instability is modified in the presence of turbulence.

We present new Hubble Space Telescope STIS, high-resolution optical imaging of a sample of 13 submillimeter (submm) luminous galaxies, for which the optical emission has been pinpointed either through radio-1.4 GHz or millimeter interferometry. We find a predominance of irregular and complex morphologies in the sample, suggesting that mergers are likely common for submm galaxies. The component separation in these objects are on average a factor two larger than local galaxies with similarly high bolometric luminosities. The sizes and star formation rates of the submm galaxies are consistent with the maximal star formation rate densities of 20 Msun kpc^{-2} in local starburst galaxies (Lehnert & Heckman 1996). We derive quantitative morphological information for the optical galaxies hosting the submm emission; total and isophotal magnitudes, Petrosian radius, effective radius, concentration, aspect ratio, surface brightness, and asymmetry. We compare these morphological indices with those of other galaxies lying within the same STIS images. Most strikingly, we find ~70% of the submm galaxies to be extraordinarily large and elongated relative to the field population, regardless of optical magnitude. Comparison of the submm galaxy morphologies with those of optically selected galaxies at z~2-3 reveal the submm galaxies to be a morphologically distinct population, with generally larger sizes, higher concentrations and more prevalent major-merger configurations.

We discuss a new mechanism of dust acceleration that acts in a turbulent magnetized medium. The magnetohydrodynamic (MHD) turbulence includes both fluid motions and magnetic fluctuations. We show that while the fluid motions bring about grain motions through the drag, the electromagnetic fluctuations, can accelerate grains through resonant interactions. In this paper we calculate grain acceleration by gyroresonance in Cold Neutral Medium. We consider both incompressible and compressible MHD modes. We show that fast modes dominate grain acceleration. For the parameters chosen, fast modes render supersonic velocities to grains that may shatter grains and enable efficient absorption of heavy elements. Since the grains are preferentially accelerated with large pitch angles, the supersonic grains get aligned with long axes perpendicular to the magnetic field.

Recent results suggest that z ~ 6 marks the end of the reionization era. A large sample of objects at z ~ 6, therefore, will be of enormous importance, as it will enable us to observationally determine the exact epoch of the reionization and the sources that are responsible for it. With the HST Advanced Camera for Surveys (ACS) coming on line, we now have an unique opportunity to discover a significant number of objects at z ~ 6. The pure parallel mode implemented for the Wide Field Camera (WFC) has greatly enhanced this ability. We present our preliminary analysis of a deep ACS/WFC parallel field at |b|=74.4^o. We find 30 plausible z ~ 6 candidates, all of which have S/N > 7 in the F850LP-band. The major source of contamination could be faint Galactic cool dwarfs, and we estimated that they would contribute at most 4 objects to our candidate list. We derived the cumulative number density of galaxies at 6.0 <= z <= 6.5 as 2.3 arcmin^{-2} to a limit of 28.0 mag in the F850LP-band, which is slightly higher than our prediction. If this is not due to an underestimated contamination rate, it could possibly imply that the faint-end slope of the z ~ 6 luminosity function is steeper than alpha=-1.6. At the very least, our result suggests that galaxies with L<L^* do exist in significant number at z ~ 6 and they could be the major sources that contributed to the reionizing photons.

Understanding of the nature of gamma-ray bursts (GRBs) is one of the challenging problem facing the astrophysics community. The field of gamma ray bursts is a rapidly developing one and we expect that new missions, like SWIFT and GLAST will bring the field to a new quantitative level. The detail of the explosion, the formation of the fireball, its propagation, the generation of the shocks, the source of the magnetic fields and soft photons, the particle acceleration process, and the details of the radiation process are all still outstanding questions. The goal of this review is to attract the attention of the community to a number of physical processes (MHD turbulence and particle acceleration, magnetic reconnection, etc.) that seem to be very relevant for these sources.

Recent advances in understanding of magnetohydrodynamic (MHD) turbulence call for substantial revisions in our understanding of cosmic ray transport. In this paper we use gyroresonance recently obtained scaling laws for MHD modes to calculate the scattering frequency for cosmic rays in the ISM. We consider gyroresonance with MHD modes (Alfvenic, slow and fast) and transit-time damping (TTD) by fast modes. We conclude that the gyroresonance with fast modes is the dominant contribution to cosmic ray scattering for the typical interstellar conditions. In contrast to earlier studies, we find that Alfvenic and slow modes are inefficient because they are far from isotropy usually assumed.

Turbulence is the most common state of astrophysical flows. In typical astrophysical fluids, turbulence is accompanied by strong magnetic fields, which has a large impact on the dynamics of the turbulent cascade. Recently, there has been a significant breakthrough on the theory of magnetohydrodynamic (MHD) turbulence. For the first time we have a scaling model that is supported by both observations and numerical simulations. We review recent progress in studies of both incompressible and compressible turbulence. We compare Iroshnikov-Kraichnan and Goldreich-Sridhar models, and discuss scalings of Alfven, slow, and fast waves. We discuss the implications of this new insight into MHD turbulence for cosmic ray transport.

The Chandra Deep Field South (CDF-S) 1Ms exposure produced a catalog of 346 X-ray sources, of which 59 were not visible on the VLT/FORS1 and the ESO-MPI/WFI deep R-band images to a limit of R_{vega}=26.1--26.7 mag. Using the first release of the ESO VLT/ISAAC JHK_s data on the CDF-S, we identified six of the twelve such objects that were within the coverage of these IR observations. The VLT/FORS1 I-band data further confirms that five of these six objects are undetected in the optical. The photometric properties of these six counterparts are compared against those of the optically brighter counterparts of Chandra sources in the same field. We found that the location of these optically brighter Chandra sources in the near-IR color space was bifurcated, with the color of one branch being consistent with that of E/S0 galaxies at 0 <= z <= 1.5 and the other branch being consistent with that of unreddened AGN/QSOs at 0 <= z <= 3.5. The six counterparts that we identified seemed to lie on the E/S0 branch and its extension, suggesting that these X-ray source hosts are mostly luminous E/S0 galaxies (M_V ~ -20 mag in AB system) at 1 <= z <= 2.5. On the other hand, some of them can also be explained by AGN/QSOs over a wide redshift range (0 <= z <= 5) if a range of internal extinction (A_V=0--1 mag) is allowed. However, the later interpretation requires fine-tuning extinction together with redshift for these objects individually. If they are indeed AGN/QSOs, the most luminous of them is just barely qualified for being a QSO. Finally, we point out that neither high-redshift (z > 5) star-forming galaxies nor irregular galaxies at lower redshift can be a viable explanation to the nature of these six counterparts.

Systematic surveys are being proposed to discover a significant number of galaxies at z ~ 6, which is now suggested as the epoch when the reionization era of the Universe ends. To plan such surveys, we need a reasonable expectation of the surface density of high redshift galaxies at different flux limits. Here we present a simple prediction of the surface density of 5.5 < z < 6.5 galaxies in the optical regime, extrapolating from what is already known about galaxies at z ~ 3. This prediction is consistent with the results of nearly all known searches for objects at z ~ 6, giving confidence that we may use it to plan optimal combination of survey depth and sky coverage in searching for such objects. We suggest that the most efficient strategy with existing ground-based facilities is to do medium-depth (m_{AB}= 24.0 -- 24.5 mag), wide-field (a couple of square degrees) survey using a wide-field camera at a 4m-class telescope. As the predicted surface density at this brightness level is very sensitive to the value of L*, the result of such a survey can be easily used to constrain the luminosity evolution from z ~ 3 to 6.

We discuss the relative grain motions due to MHD turbulence in interstellar medium. It has been known for decades that turbulent drag is an efficient way to induce grain relative motions. However, earlier treatments disregarded magnetic field and used Kolmogorov turbulence. Unlike hydro turbulence, MHD turbulence is anisotropic on small scales. Moreover, compressible modes are important for MHD and magnetic perturbations can directly interact with grains. We provide calculations of grain relative motion for realistic interstellar turbulence driving that is consistent with the velocity dispersions observed in diffuse gas and for realistic grain charging. We account for the turbulence cutoff arising from abmipolar drag. Our results on grain shattering are consistent with the customary accepted cutoff size. We obtain grain velocities for turbulence with parameters consistent with those in HI and dark clouds. These velocities are smaller than those in earlier papers, where MHD effects were disregarded. Finally, we consider grain velocities arising from photoelectric emission, radiation pressure and the thrust due to molecular hydrogen formation. These are lower than relative velocities induced by turbulence. We conclude that turbulence should prevent these mechanisms from segregating grains by size.

Interstellar medium is turbulent and this induces relative motions of dust grains. We calculate relative velocities of charged grains in a partially ionized magnetized gas. We account for anisotropy of magnetohydrodynamic (MHD) turbulence, grain coupling with magnetic field, and the turbulence cutoff arising from the ambipolar drag. We obtain grain velocities for turbulence with parameters consistent with those in HI and dark clouds. Those velocities are smaller than those in earlier papers, where MHD effects were disregarded. Finally, we compare grain velocities arising from photoelectric emission, radiation pressure and thrusts from molecular hydrogen formation. We conclude that turbulence should prevent segregation of grains of different sizes.

Galactic HI is a gas that is coupled to magnetic field because of its fractional ionization. Many properties of HI are affected by turbulence. Recently, there has been a significant breakthrough on the theory of magnetohydrodynamic (MHD) turbulence. For the first time in the history of the subject we have a scaling model that is supported by numerical simulations. We review recent progress in studies of both incompressible and compressible turbulence. We also discuss the new regime of MHD turbulence that happens below the scale at which conventional turbulent motions get damped by viscosity. The viscosity in the case of HI is being produced by neutrals and truncates the turbulent cascade at parsec scales. We show that below this scale magnetic fluctuations with a shallow spectrum persist and point out to a possibility of the resumption of the MHD cascade after ions and neutrals decouple. We discuss the implications of the new insight into MHD turbulence for cosmic ray transport and grain dynamics.

The edge-on, nearby spiral galaxy NGC5907 has long been used as the prototype of a ``non-interacting'' warped galaxy. We report here the discovery of two interactions with companion dwarf galaxies that substantially change this picture. First, a faint ring structure is discovered around this galaxy that is likely due to the tidal disruption of a companion dwarf spheroidal galaxy. The ring is elliptical in shape with the center of NGC5907 close to one of the ring's foci. This suggests the ring material is in orbit around NGC5907. No gaseous component to the ring has been detected either with deep H$\alpha$ images or in Very Large Array (VLA) HI 21-cm line maps. The visible material in the ring has an integrated luminosity <= 10^8 L_sun and its brightest part has a color R-I ~ 0.9. All of these properties are consistent with the ring being a tidally-disrupted dwarf spheroidal galaxy. Second, we find that NGC5907 has a dwarf companion galaxy, PGC54419, projected to be only 36.9 kpc from the center of NGC5907, close in radial velocity (\Delta V = 45 km s^{-1}) to the giant spiral galaxy. This dwarf is seen at the tip of the HI warp and in the direction of the warp. Hence, NGC5907 can no longer be considered ``non-interacting,'' but is obviously interacting with its dwarf companions much as the Milky Way interacts with its dwarf galaxies. These results, coupled with the finding by others that dwarf galaxies tend to be found around giant galaxies, suggest that tidal interaction with companions, even if containing a mere 1% of the mass of the parent galaxy, might be sufficient to excite the warps found in the disks of many large spiral galaxies.

We present nine color CCD intermediate-band spectrophotometry of a two square degree field centered on the old open cluster M67, from 3890$\rm \AA$ to nearly 1$\mu$. These observations are taken as a part of the BATC (Beijing-Arizona-Taipei-Connecticut) Color Survey of the Sky, for both scientific and calibration reasons. With these data we show that the BATC survey can reach its goal of obtaining spectrophotometry to a zero point accuracy of 0.01 mag, and down to V = 21 with 0.3 mag random error. We fit the color-magnitude diagrams (CMDs) with Worthey's theoretical models. The net result is the excellent fit of the 4.0 Gyr, [Fe/H] = $-0.10$ model to our data, including a good fit to the main sequence (MS) turn-off. Our data are consistent with a toy model with 50\% of the stars in M67 being binaries and a random distribution of binary mass-ratios, although other models with different mass-ratio distributions cannot be ruled out. The spatial distribution and mass function (MF) of stars in M67 show marked effects of dynamical evolution and evaporation of stars from the cluster. Blue stragglers and binary stars are the most condensed within the cluster, with degree of condensation depending on mass.We find M67 to have an elongated shape, oriented at an angle of $15^{\circ}$ relative to the galactic plane. Within its tidal radius, the observed MF of M67 between 1.2 $\rm M_\odot$ and $\rm 0.8 M_\odot$ has a Salpeter slope $\rm \eta = -1.93 \pm 0.66$. For stars of mass below 0.8 $\rm M_\odot$, $\rm \eta \sim 0$. It is plausible that the leveling-off of the MF at lower masses is a result of evaporation of lower mass stars in this mass range at a rate of one every $\sim 10^7$ years. If so, it is plausible that the IMF of M67 has the canonical field value of $\rm \eta = -2.0$.