In an effort to probe the origin of surface brightness profile (SBP) breaks widely observed in nearby disk galaxies, we carry out a comparative study of stellar population profiles of 635 disk galaxies selected from the MaNGA spectroscopic survey. We classify our galaxies into single exponential (TI), down-bending (TII) and up-bending (TIII) SBP types, and derive their spin parameters and radial profiles of age/metallicity-sensitive spectral features. Most TII (TIII) galaxies have down-bending (up-bending) star formation rate (SFR) radial profiles, implying that abrupt radial changes of SFR intensities contribute to the formation of both TII and TIII breaks. Nevertheless, a comparison between our galaxies and simulations suggests that stellar migration plays a significant role in weakening down-bending $\Sigma_{\star}$ profile breaks. While there is a correlation between the break strengths of SBPs and age/metallicity-sensitive spectral features for TII galaxies, no such correlation is found for TIII galaxies, indicating that stellar migration may not play a major role in shaping TIII breaks, as is evidenced by a good correspondence between break strengths of $\Sigma_{\star}$ and surface brightness profiles of TIII galaxies. We do not find evidence for galaxy spin being a relevant parameter for forming different SBP types, nor do we find significant differences between the asymmetries of galaxies with different SBP types, suggesting that environmental disturbances or satellite accretion in the recent past do not significantly influence the break formation. By dividing our sample into early and late morphological types, we find that galaxies with different SBP types follow nearly the same tight stellar mass-$R_{25}$ relation, which makes the hypothesis that stellar migration alone can transform SBP types from TII to TI and then to TIII highly unlikely.

The recently reported result from XENON1T experiment indicates that there is an excess with $3.5\sigma$ significance in the electronic recoil events. Interpreted as new physics, new sources are needed to account for the electronic scattering. Here, we show that a dark fermion $\psi$ from semi-annihilation of $Z_3$ dark matter $X$ and subsequent decay may be responsible for the excess. The relevant semi-annihilation process is $X+X\rightarrow \overline{X}+V_\mu (\rightarrow \psi + \overline{\psi})$, in which the final $\psi$ has a box-shape energy spectrum. The fast-moving $\psi$ can scatter with electron through an additional gauge boson that mixes with photon kinetically. The shape of the signals in this model can be well consistent with the observed excess.

Inflation in the early universe can generate the nearly conformal invariant fluctuation that leads to the structures we observe at the present. The simple viable Starobinsky $R^2$ inflation has an approximate global scale symmetry. We study the conformal symmetric Weyl $\hat{R}^2$ and general $F(\hat{R})$ theories and demonstrate their equivalence to Einstein gravity coupled with a scalar and a Weyl gauge field. The scalar field in Weyl $\hat{R}^2$ gravity can be responsible for inflation with Starobinsky model as the attractor, potentially distinguishable from the latter by future experiments. The intrinsic Weyl gauge boson becomes massive once the Einstein frame is fixed, and constitutes as a dark matter candidate with mass up to $\sim 5\times 10^{16}$GeV.

We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.

Context. Planet formation is expected to take place in the first million years of a planetary system through various processes, which remain to be tested through observations. Aims. With the recent discovery, using ALMA, of two gaseous spiral arms inside the 120 au cavity and connected to dusty spirals, the famous protoplanetary disk around AB Aurigae presents a strong incentive for investigating the mechanisms that lead to giant planet formation. A candidate protoplanet located inside a spiral arm has already been claimed in an earlier study based on the same ALMA data. Methods. We used SPHERE at the Very Large Telescope (VLT) to perform near-infrared (IR) high-contrast imaging of AB Aur in polarized and unpolarized light in order to study the morphology of the disk and search for signs of planet formation. Results. SPHERE has delivered the deepest images ever obtained for AB Aur in scattered light. Among the many structures that are yet to be understood, we identified not only the inner spiral arms, but we also resolved a feature in the form of a twist in the eastern spiral at a separation of about 30 au. The twist of the spiral is perfectly reproduced with a planet-driven density wave model when projection effects are accounted for. We measured an azimuthal displacement with respect to the counterpart of this feature in the ALMA data, which is consistent with Keplerian motion on a 4-yr baseline. Another point sxce is detected near the edge of the inner ring, which is likely the result of scattering as opposed to the direct emission from a planet photosphere. We tentatively derived mass constraints for these two features. Conclusions. The twist and its apparent orbital motion could well be the first direct evidence of a connection between a protoplanet candidate and its manifestation as a spiral imprinted in the gas and dust distributions.

Conformal scaling invariance should play an important role for understanding the origin and evolution of universe. During inflation period, it appears to be an approximate symmetry, but how it is broken remains uncertain. The appealing $\alpha$-attractor inflation implements the spontaneous breaking of conformal symmetry and a mysterious SO(1,1) global symmetry. To better understand the SO(1,1) symmetry, here we present a systematic treatment of the inflation models with local conformal symmetry in a more general formalism. We find SO(2) is the other possible symmetry in the presence of Weyl gauge field. We also obtain all the analytic solutions that relate the inflation fields between Jordan frame and Einstein frame. We illustrate a class of inflation models with the approximate SO(2) symmetry and trigonometric potential, and find that it can fit the current observations and will be probed by future CMB experiments.

We present our observational results of the 0.87 mm polarized dust emission in the Class 0 protostar B335 obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) at a 0.2" (20 au) resolution. We compared our data at 0.87 mm with those at 1.3 mm from the ALMA archive. The observed polarization orientations at the two wavelengths are consistent within the uncertainty, and the polarization percentages are systematically higher at 1.3 mm than 0.87 mm by a factor of ~1.7, suggesting that the polarized emission originates from magnetically aligned dust grains. We inferred the magnetic field orientations from the observed polarization orientations. We found that the magnetic field changes from ordered and highly pinched to more complicated and asymmetric structures within the inner 100 au scale of B335, and the magnetic field connects to the center along the equatorial plane as well as along the directions which are ~40-60 degrees from the equatorial plane. We performed non-ideal MHD simulations of collapsing dense cores. We found that similar magnetic field structures appear in our simulations of dense cores with the magnetic field and rotational axis slightly misaligned by 15 degrees but not in those with the aligned magnetic field and rotational axis. Our results suggest that the midplane of the inner envelope within the inner 100 au scale of B335 could be warped because of the misaligned magnetic field and rotational axis, and the magnetic field could be dragged by the warped accretion flows.

Motivated by the Weyl scaling gauge symmetry, we present a theoretical framework to explain cosmic inflation and dark matter simultaneously. This symmetry has been resurrected in recent attempts to formulate the gauge theory of gravity. We show the inspired inflation model is well consistent with current observations and will be probed further by future experiments. Furthermore, we clarify and prove the stability of Weyl gauge boson in the general theory with multiple scalars. We show the massive Weyl gauge boson can be a dark matter candidate and give the correct relic abundance.

We used new ALMA $^{13}$CO and C$^{18}$O(3-2) observations obtained at high angular resolution ($\sim$0.2") together with previous CO(3-2) and (6-5) ALMA data and continuum maps at 1.3 and 0.8 mm in order to determine the gas properties (temperature, density, and kinematics) in the cavity and to a lesser extent in the outer disk of GG Tau A, the prototype of a young triple T Tauri star that is surrounded by a massive and extended Keplerian outer disk. By deprojecting, we studied the radial and azimuthal gas distribution and its kinematics. We also applied a new method to improve the deconvolution of the CO data and in particular better quantify the emission from gas inside the cavity. We perform local and nonlocal thermodynamic equilibrium studies in order to determine the excitation conditions and relevant physical parameters inside the ring and in the central cavity. Residual emission after removing a smooth-disk model indicates unresolved structures at our angular resolution, probably in the form of irregular rings or spirals. The outer disk is cold, with a temperature $<20$ K beyond 250 au that drops quickly (r$^{-1}$). The kinematics of the gas inside the cavity reveals infall motions at about 10% of the Keplerian speed. We derive the amount of gas in the cavity, and find that the brightest clumps, which contain about 10% of this mass, have kinetic temperatures 40$-$80 K, CO column densities of a few 10$^{17}$ cm$^{-2}$, and H$_2$ densities around 10$^7$ cm$^{-3}$. Although the gas in the cavity is only a small fraction of the disk mass, the mass accretion rate throughout the cavity is comparable to or higher than the stellar accretion rate. It is accordingly sufficient to sustain the circumstellar disks on a long timescale.

We present CS(7-6) line maps toward the central parsec of the Galactic Center (GC), conducted with the Atacama Large Millimeter/submillimeter Array (ALMA). The primary goal is to find and characterize the gas structure in the inner cavity of the circumnuclear disk (CND) in high resolution (1.3"=0.05 pc). Our large field-of-view mosaic maps -- combining interferometric and single-dish data that recover extended emission - provide a first homogeneous look to resolve and link the molecular streamers in the CND with the neutral nuclear filaments newly detected within the central cavity of the CND. We find that the nuclear filaments are rotating with Keplerian velocities in a nearly face-on orbit with an inclination angle of ~10-20 degree (radius <= 0.5 pc). This is in contrast to the CND which is highly inclined at ~65-80 degree (radius ~2-5 pc). Our analysis suggests a highly warped structure from the CND to the nuclear filaments. This result may hint that the nuclear filaments and the CND were created by different external clouds passing by Sgr A*.

Rosetta observations of 67P/Churyumov-Gerasimenko (67P) reveal that most changes occur in the fallback-generated smooth terrains, vast deposits of granular material blanketing the comet's northern hemisphere. These changes express themselves both morphologically and spectrally across the nucleus, yet we lack a model that describes their formation and evolution. Here we present a self-consistent model that thoroughly explains the activity and mass loss from Hapi's smooth terrains. Our model predicts the removal of dust via re-radiated solar insolation localized within depression scarps that are substantially more ice-rich than previously expected. We couple our model with numerous Rosetta observations to thoroughly capture the seasonal erosion of Hapi's smooth terrains, where local scarp retreat gradually removes the uppermost dusty mantle. As sublimation-regolith interactions occur on rocky planets, comets, icy moons and KBOs, our coupled model and observations provide a foundation for future understanding of the myriad of sublimation-carved worlds.

We present the POL-2 850 $\mu$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $\mu$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $\mu$G and a criticality criterion $\lambda_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component.

We consider the simplest possibility for a model of particle dark matter in which dark matter has only gravitational interaction with the standard model sector. Even in such a case, it is known that the gravitational particle production in an expanding universe may lead to a correct relic abundance depending on the inflation scale and the mass of dark matter particle. We provide a comprehensive and systematic analysis of the gravitational particle production of fermionic and vectorial dark matter, and emphasize that particles which are much heavier than the Hubble parameter but lighter than inflaton can also be produced abundantly.

(Abbreviated) We investigate the interplay between magnetic (B) field, gravity, and turbulence in the fragmentation process of cores within the filamentary infrared dark cloud G34.43+00.24. We observe the magnetic field (B) morphology across G34.43 and compare with the kinematics obtained from N2H+ across the filament. We derive local velocity gradients from N2H+, tracing motion in the plane of sky, and compare with the observed local B field orientations in the plane of sky. Besides a large-scale east-west velocity gradient, we find a close alignment between local B field orientations and local velocity gradients toward the MM1/MM2 ridge. This local correlation in alignment suggests that gas motions are influenced by the magnetic field morphology or vice versa. Additionally, this alignment seems to be getting even closer with increasing integrated emission in N2H+, possibly indicating that a growing gravitational pull is more and more aligning B field and gas motion. We analyze and quantify B field, gravity, turbulence, and their relative importance toward the MM1, MM2 and MM3 regions with various techniques over two scales, a larger clump area at 2 pc scale and the smaller core area at 0.6 pc scale. While gravitational energy, B field, and turbulent pressure all grow systematically from large to small scale, the ratios among the three constituents develop clearly differently over scale. We propose that this varying relative importance between B field, gravity, and turbulence over scale drives and explains the different fragmentation types seen at sub-pc scale (no fragmentation in MM1; aligned fragmentation in MM2; clustered fragmentation in MM3). We discuss uncertainties, subtleties, and the robustness of our conclusion, and we stress the need of a multi-scale joint analysis to understand the dynamics in these systems.

We present the 850 $\mu$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfv\'{e}nic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important.

We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

We report 850~$\mu$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $\rho$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $\mu$G, 136 $\pm$ 69 $\mu$G, and 213 $\pm$ 115 $\mu$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties.

We present our analysis of the magnetic field structures from 6000 au to 100 au scales in the Class 0 protostar B335 inferred from our JCMT POL-2 observations and the ALMA archival polarimetric data. To interpret the observational results, we perform a series of (non-)ideal MHD simulations of the collapse of a rotating non-turbulent dense core, whose initial conditions are adopted to be the same as observed in B335, and generate synthetic polarization maps. The comparison of our JCMT and simulation results suggests that the magnetic field on a 6000 au scale in B335 is pinched and well aligned with the bipolar outflow along the east-west direction. Among all our simulations, the ALMA polarimetric results are best explained with weak magnetic field models having an initial mass-to-flux ratio of 9.6. However, we find that with the weak magnetic field, the rotational velocity on a 100 au scale and the disk size in our simulations are larger than the observational estimates by a factor of several. An independent comparison of our simulations and the gas kinematics in B335 observed with the SMA and ALMA favors strong magnetic field models with an initial mass-to-flux ratio smaller than 4.8. We discuss two possibilities resulting in the different magnetic field strengths inferred from the polarimetric and molecular-line observations, (1) overestimated rotational-to-gravitational energy in B335 and (2) additional contributions in the polarized intensity due to scattering on a 100 au scale.

We propose that there may be a substantial stochastic gravitational wave background from particle origin, mainly from the gravitational three-body decay of the inflaton. The emitted gravitons could constitute a sizable contribution to dark radiation if the mass of inflaton is close to the Planck scale, which can be probed by future CMB experiments that have a sensitivity on the deviation of the effective number of neutrinos in the standard cosmology, $\delta N_{\textrm{eff}}\sim 0.02 - 0.03$. We have also illustrated the spectrum of the radiated gravitational waves, in comparison to the current and future experiments, and found that gravitational waves from particle origin could be the dominant contribution to the energy density at high-frequency domain, but beyond the sensitivity regions of various experiments.

We consider the thermal effects into the evaluation of the dark matter production process. With the assistance of the right handed neutrinos, the freeze-in massive particle dark matter production history can be modified by the two-step phase transitions. The kinematic of decay/inverse decay or annihilation processes can be affected by the finite temperature effects as the Universe cools down. The history of the symmetry respected by the model can be revealed by the DM relic abundance evolution processes. The strong first order electroweak phase transition generated gravitational waves can be probed. The number of extra scalars for the Hierarchy problem can be probed through the Higgs off-shell searches at the LHC.

We present an asteroseismic analysis of the Kepler light curve of KIC 8263801, a red-giant star in the open cluster NGC 6866 that has previously been reported to be a helium-burning red-clump star. We extracted the frequencies of the radial and quadrupole modes from its frequency power spectrum and determined its properties using a grid of evolutionary models constructed with MESA. The oscillation frequencies were calculated using the GYRE code and the surface term was corrected using the Ball & Gizon(2014) prescription. We find that the star has a mass of $M/M_{\odot} = 1.793\pm 0.072$, age $t=1.48\pm 0.21$ Gyr and radius $R/R_{\odot} = 10.53\pm 0.28$. By analyzing the internal structure of the best-fitting model, we infer the evolutionary status of the star KIC 8263801 as being on the ascending part of the red giant branch, and not on the red clump. This result is verified using a purely asteroseismic diagnostic, the $\epsilon_{c}-\Delta\nu_{c}$ diagram which can distinguish red giant branch stars from red clump stars. Finally, by comparing its age with NGC 6866 ($t=0.65 \pm 0.1$ Gyr) we conclude that KIC 8263801 is not a member of this open cluster.

Motivated by the gauge theory of gravity with local scaling symmetry proposed recently in 1712.04537 and 1506.01807, we investigate whether the scalar field therein can be responsible for the inflation. We show that the classical theory would suffer from the difficulty that inflation can start but will never stop. We explore possible solutions by invoking the symmetry breaking through quantum effects. The effective potential of the scalar field is shown to have phenomenologically interesting forms to give viable inflation models. The predictions of physical observables agree well with current cosmological measurements and can be further tested in future experiments searching for primordial gravitational waves.

Studying molecular species in protoplanetary disks is very useful to characterize the properties of these objects, which are the site of planet formation. We attempt to constrain the chemistry of S-bearing molecules in the cold parts of circumstellar disk of GG Tau A. We searched for H$_2$S, CS, SO, and SO$_2$ in the dense disk around GG Tau A with the NOrthem Extended Millimeter Array (NOEMA) interferometer. We detected H$_2$S emission from the dense and cold ring orbiting around GG Tau A. This is the first detection of H$_2$S in a protoplanetary disk. We also detected HCO$^+$, H$^{13}$CO$^+$, and DCO$^+$ in the disk. Upper limits for other molecules, CCS, SO$_2$, SO, HC$_3$N, and $c$-C$_3$H$_2$ are also obtained. The observed DCO$^+$/HCO$^+$ ratio is similar to those in other disks. The observed column densities, derived using our radiative transfer code DiskFit, are then compared with those from our chemical code Nautilus. The column densities are in reasonable agreement for DCO$^{+}$, CS, CCS, and SO$_2$. For H$_2$S and SO, our predicted vertical integrated column densities are more than a factor of 10 higher than the measured values. Our results reinforce the hypothesis that only a strong sulfur depletion may explain the low observed H$_2$S column density in the disk. The H$_2$S detection in GG Tau A is most likely linked to the much larger mass of this disk compared to that in other T Tauri systems.

Recently the EDGES experiment reported an enhanced 21cm absorption signal in the radio wave observation, which may be interpreted as either anomalous cooling of baryons or heating of cosmic microwave background photons. In this paper, we pursue the latter possibility. We point out that dark radiation consisting of axion-like particles can resonantly convert into photons under the intergalactic magnetic field, which can effectively heat up the radiation in the frequency range relevant for the EDGES experiment. This may explain the EDGES anomaly.

Magnetic field is one of the key agents that play a crucial role in shaping molecular clouds and regulating star formation, yet the complete information on the magnetic field is not well constrained due to the limitations in observations. We study the magnetic field in the massive infrared dark cloud G035.39-00.33 from dust continuum polarization observations at 850 $\micron$ with SCUBA-2/POL-2 at JCMT. The magnetic field tends to be perpendicular to the densest part of the main filament (F$_{M}$), whereas it has a less defined relative orientation in the rest of the structure, where it tends to be parallel to some diffuse regions. A mean plane-of-the-sky magnetic field strength of $\sim$50 $\mu$G for F$_{M}$ is obtained using Davis-Chandrasekhar-Fermi method. Based on $^{13}$CO (1-0) line observations, we suggest a formation scenario of F$_{M}$ due to large-scale ($\sim$10 pc) cloud-cloud collision. Using additional NH$_3$ line data, we estimate that F$_{M}$ will be gravitationally unstable if it is only supported by thermal pressure and turbulence. The northern part of F$_{M}$, however, can be stabilized by a modest additional support from the local magnetic field. The middle and southern parts of F$_{M}$ are likely unstable even if the magnetic field support is taken into account. We claim that the clumps in F$_{M}$ may be supported by turbulence and magnetic fields against gravitational collapse. Finally, we identified for the first time a massive ($\sim$200 M$_{\sun}$), collapsing starless clump candidate, "c8", in G035.39-00.33. The magnetic field surrounding "c8" is likely pinched, hinting at an accretion flow along the filament.

Utilizing James Clark Maxwell Telescope (JCMT) 850 $\mu$m SCUPOL dust polarization data, we investigate the configuration of the magnetic ($B$) field in the circumnuclear disk (CND) of the Galactic Center (GC). The SCUPOL data show a highly improved polarization coverage and resolution compared to earlier 100 $\mu$m observations. The 850 $\mu$m data have a resolution and coverage similar to previous 350 $\mu$m polarimetry data. However, with a proper sampling on a 10$"$ grid, we find the 850 $\mu$m data trace the morphological structures of the CND substantially better. Furthermore, because the 850 $\mu$m trace the field deeper into the material near Sgr A*, they represent the highest resolution submillimeter probe to date of the CND magnetic field. The observed $B$-field morphology is well described by a self-similar axisymmetric disk model where the radial infall velocity is one quarter of the rotational velocity. A detailed comparison with higher-resolution interferometric maps from the Submillimeter Array further reveals that the $B$-field aligns with the neutral gas streamers connecting to the CND. Moreover, the innermost observed $B$-field structure also appears to trace and align with the mini-spiral located inside the CND. This suggests that there is one underlying $B$-field structure that is connecting the CND with its streamers and the inner mini-spiral. An estimate of $\beta_{\rm Plasma} \lesssim 1-$based on the global B-field morphology that constrains the azimuthal-to-vertical field strength ratio of around 40 combined with a measurement of the azimuthal velocity--indicates that the B-field appears dynamically significant towards the CND and also onwards to the inner mini-spiral.

We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $\mu$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $\lambda$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical.

We present 850 $\mu$m imaging polarimetry data of the $\rho$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We identify 10 subregions with distinct polarization fractions and angles in the 0.2 pc $\rho$ Oph A core; some of them can be part of a coherent magnetic field structure in the $\rho$ Oph region. The results are consistent with previous observations of the brightest regions of $\rho$ Oph-A, where the degrees of polarization are at a level of a few percents, but our data reveal for the first time the magnetic field structures in the fainter regions surrounding the core where the degree of polarization is much higher ($> 5 \%$). A comparison with previous near-infrared polarimetric data shows that there are several magnetic field components which are consistent at near-infrared and submillimeter wavelengths. Using the Davis-Chandrasekhar-Fermi method, we also derive magnetic field strengths in several sub-core regions, which range from approximately 0.2 to 5 mG. We also find a correlation between the magnetic field orientations projected on the sky with the core centroid velocity components.

In the purely gravitational dark matter scenario, the dark matter particle does not have any interaction except for gravitational one. We study the gravitational particle production of dark matter particle in such a minimal setup and show that correct amount of dark matter can be produced depending on the inflation model and the dark matter mass. In particular, we carefully evaluate the particle production rate from the transition epoch to the inflaton oscillation epoch in a realistic inflation model and point out that the gravitational particle production is efficient even if dark matter mass is much larger than the Hubble scale during inflation as long as it is smaller than the inflaton mass.

HD 176465 is a binary system whose both components are solar-like pulsators and whose oscillation frequencies were observed by the \kepler\ mission. In this paper we have modeled the asteroseismic and spectroscopic data of the stars, and have determined their convection-zone helium abundances using the signatures left by the He{\sc ii} ionization zone on the mode frequencies. As expected we find that the components of the binary have the same age within uncertainties ($3.087 \pm 0.580$ Gyr and $3.569 \pm 0.912$ Gyr); they also have the same initial helium abundance (Y$_{\mathrm{init}}$=0.253 $\pm$ 0.006 and 0.254 $\pm$ 0.008). Their current metallicity ([Fe/H]=$-0.275 \pm 0.04$ and $-0.285 \pm 0.04$) is also the same within errors. Fits to the signature of the He{\sc ii} acoustic glitch yields current helium abundances of $Y_{\rm A} = 0.224 \pm 0.006$ and $Y_{\rm B} = 0.233 \pm 0.008$ for the two components. Analyzing the complete ensemble of models generated for this investigation we find that both the amplitude and acoustic depth of the glitch signature arising from the second helium ionization zone and the base of the convection zone (CZ) are functions of mass. We show that the acoustic depths of these glitches are positively correlated with each other. The analysis can help us to detect the internal structure and constrain the chemical compositions.

We present Atacama Large Millimeter Array (ALMA) observations at an angular resolution of 0.1-0.2" of the disk surrounding the young Herbig Ae star MWC 758. The data consist of images of the dust continuum emission recorded at 0.88 millimeter, as well as images of the 13CO and C18O J = 3-2 emission lines. The dust continuum emission is characterized by a large cavity of roughly 40 au in radius which might contain a mildly inner warped disk. The outer disk features two bright emission clumps at radii of about 47 and 82 au that present azimuthal extensions and form a double-ring structure. The comparison with radiative transfer models indicates that these two maxima of emission correspond to local increases in the dust surface density of about a factor 2.5 and 6.5 for the south and north clumps, respectively. The optically thick 13CO peak emission, which traces the temperature, and the dust continuum emission, which probes the disk midplane, additionally reveal two spirals previously detected in near-IR at the disk surface. The spirals seen in the dust continuum emission present, however, a slight shift of a few au towards larger radii and one of the spirals crosses the south dust clump. Finally, we present different scenarios in order to explain the complex structure of the disk.

We present the first ALMA dust polarization observations towards the high-mass star-forming regions W51 e2, e8, and W51 North in Band 6 (230 GHz) with a resolution around 0.26" ($\sim5$mpc). Polarized emission in all three sources is clearly detected and resolved. Measured relative polarization levels are between 0.1\% and 10\%. While the absolute polarization shows complicated structures, the relative polarization displays the typical anti-correlation with Stokes $I$, though with a large scatter. Inferred magnetic (B) field morphologies are organized and connected. Detailed substructures are resolved, revealing new features such as cometary-shaped B-field morphologies in satellite cores, symmetrically converging B-field zones, and possibly streamlined morphologies. The local B-field dispersion shows some anti-correlation with the relative polarization. Moreover, lowest polarization percentages together with largest dispersions coincide with B-field convergence zones. We put forward $\sin\omega$, where $\omega$ is the measurable angle between a local B-field orientation and local gravity, as a measure of how effectively the B-field can oppose gravity. Maps of $\sin\omega$ for all three sources show organized structures that suggest a locally varying role of the B-field, with some regions where gravity can largely act unaffectedly, possibly in a network of narrow magnetic channels, and other regions where the B-field can work maximally against gravity.

How important is the magnetic (B-) field when compared to gravity and turbulence in the star-formation process? Does its importance depend on scale and location? We summarize submm dust polarization observations towards the large filamentary infrared dark cloud G34 and towards a dense core in the high-mass star-forming region W51. We detect B-field orientations that are either perpendicular or parallel to the G34 filament axis. These B-field orientations further correlate with local velocity gradients. Towards three cores in G34 we find a varying importance between B-field, gravity, and turbulence that seems to dictate varying types of fragmentation. At highest resolution towards the gravity-dominated collapsing core W51 e2 we resolve new B-field features, such as converging B-field lines and possibly magnetic channels.

We present the results of 2D, moving-mesh, viscous hydrodynamical simulations of an accretion disc around a merging supermassive black hole binary (SMBHB). The simulation is pseudo-Newtonian, with the BHs modeled as point masses with a Paczynski-Wiita potential, and includes viscous heating, shock heating, and radiative cooling. We follow the gravitational inspiral of an equal-mass binary with a component mass M_bh=10^6 M_sun from an initial separation of 60r_g (where r_g=GM_bh/c^2 is the gravitational radius) to the merger. We find that a central, low-density cavity forms around the binary, as in previous work, but that the BHs capture gas from the circumbinary disc and accrete efficiently via their own minidiscs, well after their inspiral outpaces the viscous evolution of the disc. The system remains luminous, displaying strong periodicity at twice the binary orbital frequency throughout the entire inspiral process, all the way to the merger. In the soft X-ray band, the thermal emission is dominated by the inner edge of the circumbinary disc with especially clear periodicity in the early inspiral. By comparison, harder X-ray emission is dominated by the minidiscs, and the light curve is initially more noisy but develops a clear periodicity in the late inspiral stage. This variability pattern should help identify the EM counterparts of SMBHBs detected by the space-based gravitational-wave detector LISA.

Observations by the Atacama Large Millimeter/sub-millimeter Array of the dust continuum and $^{13}$CO(3-2) millimeter emissions of the triple stellar system GG Tau A are analysed, giving evidence for a rotating gas disc and a concentric and coplanar dust ring. The present work complements an earlier analysis (Tang et al. 2016) by exploring detailed properties of the gas disc. A 95% confidence level upper limit of 0.24 arcsec (34 au) is placed on the disc scale height at a distance of 1 arcsec (140 au) from the central stars. Evidence for Keplerian rotation of the gas disc is presented, the rotation velocity reaching ~3.1 kms$^{-1}$ at 1 arcsec from the central stars, and a 99% confidence level upper limit of 9% is placed on a possible in-fall velocity relative contribution. Variations of the intensity across the disc area are studied in detail and confirm the presence of a hot spot in the south-eastern quadrant. However several other significant intensity variations, in particular a depression in the northern direction, are also revealed. Variations of the intensity are found to be positively correlated to variations of the line width. Possible contributions to the measured line width are reviewed, suggesting an increase of the disc temperature and opacity with decreasing distance from the stars.

We investigate the relation between gas and star formation in sub-galactic regions, ~360 pc to ~1.5 kpc in size, within the nearby starburst dwarf NGC4449, in order to separate the underlying relation from the effects of sampling at varying spatial scales. Dust and gas mass surface densities are derived by combining new observations at 1.1 mm, obtained with the AzTEC instrument on the Large Millimeter Telescope, with archival infrared images in the range 8-500 micron from the Spitzer Space Telescope and the Herschel Space Observatory. We extend the dynamic range of our mm (and dust) maps at the faint end, using a correlation between the far-infrared/millimeter colors F(70)/F(1100) [and F(160)/F(1100)] and the mid-infrared color F(8)/F(24) that we establish for the first time for this and other galaxies. Supplementing our data with maps of the extinction-corrected star formation rate (SFR) surface density, we measure both the SFR-molecular gas and the SFR-total gas relations in NGC4449. We find that the SFR-molecular gas relation is described by a power law with exponent that decreases from ~1.5 to ~1.2 for increasing region size, while the exponent of the SFR-total gas relation remains constant with value ~1.5 independent of region size. We attribute the molecular law behavior to the increasingly better sampling of the molecular cloud mass function at larger region sizes; conversely, the total gas law behavior likely results from the balance between the atomic and molecular gas phases achieved in regions of active star formation. Our results indicate a non-linear relation between SFR and gas surface density in NGC4449, similar to what is observed for galaxy samples.

In this study, we present some general features of gamma-ray spectra from dark matter. We find that the spectrum with sharp features could appear in a wide class of dark matter models and mimic the gamma line signals. If all other physical degrees of freedom are heavy or effectively decoupled, the resulting gamma ray from dark matter decay or annihilation would generally have polynomial-type spectra or power-law with positive index. We illustrate our findings in a model-independent framework with generic kinematic analysis. Similar results can also apply for cosmic ray or neutrino cases.

We investigate the particle production from thermal gravitational annihilation in the very early universe, which is an important contribution for particles that might not be in thermal equilibrium or/and only have gravitational interaction, such as dark matter (DM). For particles with spin 0, 1/2 and 1 we calculate the relevant cross sections through gravitational annihilation and give the analytic formulas with full mass-dependent terms. We find that DM with mass between TeV and $10^{16}$GeV could have the relic abundance that fits the observation, with small dependence on its spin. We also discuss the effects of gravitational annihilation from inflatons. Interestingly, contributions from inflatons could be dominant and have the same power dependence on Hubble parameter of inflation as that from vacuum fluctuation. Also, fermion production from inflatons, in comparison to boson, is suppressed by its mass due to helicity selection.

In the light of recent possible tensions in the Hubble constant $H_0$ and the structure growth rate $\sigma_8$ between the Planck and other measurements, we investigate a hidden-charged dark matter (DM) model where DM interacts with hidden chiral fermions, which are charged under the hidden SU(N) and U(1) gauge interactions. The symmetries in this model assure these fermions to be massless. The DM in this model, which is a Dirac fermion and singlet under the hidden SU(N), is also assumed to be charged under the U(1) gauge symmetry, through which it can interact with the chiral fermions. Below the confinement scale of SU(N), the hidden quark condensate spontaneously breaks the U(1) gauge symmetry such that there remains a discrete symmetry, which accounts for the stability of DM. This condensate also breaks a flavor symmetry in this model and Nambu-Goldstone bosons associated with this flavor symmetry appear below the confinement scale. The hidden U(1) gauge boson and hidden quarks/Nambu-Goldstone bosons are components of dark radiation (DR) above/below the confinement scale. These light fields increase the effective number of neutrinos by $\delta N_{\rm eff}\simeq 0.59$ above the confinement scale for $N=2$, resolving the tension in the measurements of the Hubble constant by Planck and Hubble Space Telescope if the confinement scale is $\lesssim 1$ eV. DM and DR continuously scatter with each other via the hidden U(1) gauge interaction, which suppresses the matter power spectrum and results in a smaller structure growth rate. The DM sector couples to the Standard Model sector through the exchange of a real singlet scalar mixing with the Higgs boson, which makes it possible to probe our model in DM direct detection experiments. Variants of this model are also discussed, which may offer alternative ways to investigate this scenario.

Most protostars have luminosities that are fainter than expected from steady accretion over the protostellar lifetime. The solution to this problem may lie in episodic mass accretion -- prolonged periods of very low accretion punctuated by short bursts of rapid accretion. However, the timescale and amplitude for variability at the protostellar phase is almost entirely unconstrained. In "A JCMT/SCUBA-2 Transient Survey of Protostars in Nearby Star Forming Regions", we are monitoring monthly with SCUBA-2 the sub-mm emission in eight fields within nearby (<500 pc) star forming regions to measure the accretion variability of protostars. The total survey area of ~1.6 sq.deg. includes ~105 peaks with peaks brighter than 0.5 Jy/beam (43 associated with embedded protostars or disks) and 237 peaks of 0.125-0.5 Jy/beam (50 with embedded protostars or disks). Each field has enough bright peaks for flux calibration relative to other peaks in the same field, which improves upon the nominal flux calibration uncertainties of sub-mm observations to reach a precision of ~2-3% rms, and also provides quantified confidence in any measured variability. The timescales and amplitudes of any sub-mm variation will then be converted into variations in accretion rate and subsequently used to infer the physical causes of the variability. This survey is the first dedicated survey for sub-mm variability and complements other transient surveys at optical and near-IR wavelengths, which are not sensitive to accretion variability of deeply embedded protostars.

The interaction between a supermassive black hole (SMBH) and the surrounding material is of primary importance in modern astrophysics. The detection of the molecular 2-pc circumnuclear disk (CND) immediately around the Milky Way SMBH, SgrA*, provides an unique opportunity to study SMBH accretion at sub-parsec scales. Our new wide-field CS(2-1) map toward the Galactic center (GC) reveals multiple dense molecular streamers originated from the ambient clouds 20-pc further out, and connecting to the central 2 parsecs of the CND. These dense gas streamers appear to carry gas directly toward the nuclear region and might be captured by the central potential. Our phase-plot analysis indicates that these streamers show a signature of rotation and inward radial motion with progressively higher velocities as the gas approaches the CND and finally ends up co-rotating with the CND. Our results might suggest a possible mechanism of gas feeding the CND from 20 pc around 2 pc in the GC. In this paper, we discuss the morphology and the kinematics of these streamers. As the nearest observable Galactic nucleus, this feeding process may have implications for understanding the processes in extragalactic nuclei.

The protoplanetary disk around HL Tau is so far the youngest candidate of planet formation, and it is still embedded in a protostellar envelope with a size of thousands of au. In this work, we study the gas kinematics in the envelope and its possible influence on the embedded disk. We present our new ALMA cycle 3 observational results of HL Tau in the 13CO (2-1) and C18O (2-1) emission at resolutions of 0.8" (110 au), and we compare the observed velocity pattern with models of different kinds of gas motions. Both the 13CO and C18O emission lines show a central compact component with a size of 2" (280 au), which traces the protoplanetary disk. The disk is clearly resolved and shows a Keplerian motion, from which the protostellar mass of HL Tau is estimated to be 1.8+/-0.3 M$_\odot$, assuming the inclination angle of the disk to be 47 deg from the plane of the sky. The 13CO emission shows two arc structures with sizes of 1000-2000 au and masses of 3E-3 M$_\odot$ connected to the central disk. One is blueshifted and stretches from the northeast to the northwest, and the other is redshifted and stretches from the southwest to the southeast. We find that simple kinematical models of infalling and (counter-)rotating flattened envelopes cannot fully explain the observed velocity patterns in the arc structures. The gas kinematics of the arc structures can be better explained with three-dimensional infalling or outflowing motions. Nevertheless, the observed velocity in the northwestern part of the blueshifted arc structure is ~60-70% higher than the expected free-fall velocity. We discuss two possible origins of the arc structures: (1) infalling flows externally compressed by an expanding shell driven by XZ Tau and (2) outflowing gas clumps caused by gravitational instabilities in the protoplanetary disk around HL Tau.

In supersymmetric axion models, if the gravitino or axino is the lightest SUSY particle (LSP), the other is often the next-to-LSP (NLSP). We investigate the cosmology of such a scenario and point out that the lifetime of the NLSP naturally becomes comparable to the present age of the universe in a viable parameter region. This is a well-motivated example of the so-called decaying dark matter model, which is recently considered as an extension of the $\Lambda$CDM model to relax some cosmological tensions.

We present Submillimeter Array (SMA) observations at 345 GHz towards the intermediate/high-mass cluster-forming region NGC 6334 V. From the dust emission we spatially resolve three dense condensations, the brightest one presenting the typical chemistry of a hot core. The magnetic field (derived from the dust polarized emission) shows a bimodal converging pattern towards the hot core. The molecular emission traces two filamentary structures at two different velocities, separated by 2 km/s, converging to the hot core and following the magnetic field distribution. We compare the velocity field and the magnetic field derived from the SMA observations with MHD simulations of star-forming regions dominated by gravity. This comparison allows us to show how the gas falls in from the larger-scale extended dense core (~0.1 pc) of NGC 6334 V towards the higher-density hot core region (~0.02 pc) through two distinctive converging flows dragging the magnetic field, whose strength seems to have been overcome by gravity.

We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC 1 region. We see that the magnetic field lies approximately perpendicular to the famous 'integral filament' in the densest regions of that filament. Furthermore, we see an 'hour-glass' magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density north-eastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments.

We report the results of ALMA observations of a protoplanetary disk surrounding the Herbig Ae star AB Aurigae. We obtained high-resolution (0.1"; 14 au) images in $^{12}$CO (J=2-1) emission and in dust continuum at the wavelength of 1.3 mm. The continuum emission is detected at the center and at the ring with a radius of $\sim$ 120 au. The CO emission is dominated by two prominent spirals within the dust ring. These spirals are trailing and appear to be about 4 times brighter than their surrounding medium. Their kinematics is consistent with Keplerian rotation at an inclination of 23 degree. The apparent two-arm-spiral pattern is best explained by tidal disturbances created by an unseen companion located at 60--80 au, with dust confined in the pressure bumps created outside this companion orbit. An additional companion at r of 30 au, coinciding with the peak CO brightness and a large pitch angle of the spiral, would help to explain the overall emptiness of the cavity. Alternative mechanisms to excite the spirals are discussed. The origin of the large pitch angle detected here remain puzzling.

Gaseous circumbinary accretion discs provide a promising mechanism to facilitate the mergers of supermassive black holes (SMBHs) in galactic nuclei. We measure the torques exerted on accreting SMBH binaries, using 2D, isothermal, moving-mesh, viscous hydrodynamical simulations of circumbinary accretion discs. Our computational domain includes the entire inner region of the circumbinary disk with the individual black holes (BHs) included as point masses on the grid and a sink prescription to model accretion onto each BH. The BHs each acquire their own well-resolved accretion discs ("minidiscs"). We explore a range of mass removal rates for the sink prescription removing gas from the central regions of the minidiscs. We find that the torque exerted on the binary is primarily gravitational, and dominated by the gas orbiting close behind and ahead of the individual BHs. The torques from the distorted circumbinary disc farther out and from the direct accretion of angular momentum are subdominant. The torques are sensitive to the sink prescription: slower sinks result in more gas accumulating near the BHs and more negative torques, driving the binary to merger more rapidly. For faster sinks, the torques are less negative and eventually turn positive (for unphysically fast sinks). When the minidiscs are modeled as standard alpha discs, our results are insensitive to the choice of sink radius. Scaling the simulations to a binary orbital period tbin = 1yr and background disc accretion rate Mdot = 0.3MEdd in Eddington units, the binary inspirals on a timescale of 3X10^6 years, irrespective of the SMBH masses. For binaries with total mass <10^7Msun, this is shorter than the inspiral time due to gravitational wave (GW) emission alone, implying that gas discs will have a significant impact on the SMBH binary population and can affect the GW signal for Pulsar Timing Arrays.

We propose a novel particle physics model in which vector dark matter (VDM) and dark radiation (DR) originate from the same non-Abelian dark sector. We show an illustrating example where dark $SU(3)$ is spontaneously broken into $SU(2)$ subgroup by the nonzero vacuum expectation value of a complex scalar in fundamental representation of $SU(3)$. The massless gauge bosons associated with the residual unbroken $SU(2)$ constitute DR and help to relieve the tension in Hubble constant measurements between $\textit{Planck}$ and Hubble Space Telescope. In the meantime, massive dark gauge bosons associated with the broken generators are VDM candidates. Intrinsically, this non-Abelian VDM can interact with non-Abelian DR in the cosmic background, which results in a suppressed matter power spectrum and leads to a smaller $\sigma_8$ for structure formation.

By performing non-masked polarization imaging with Subaru/HiCIAO, polarized scattered light from the inner region of the disk around the GG Tau A system was successfully detected in the $H$ band with a spatial resolution of approximately 0.07$\arcsec$, revealing the complicated inner disk structures around this young binary. This paper reports the observation of an arc-like structure to the north of GG Tau Ab and part of a circumstellar structure that is noticeable around GG Tau Aa extending to a distance of approximately 28 AU from the primary star. The speckle noise around GG Tau Ab constrains its disk radius to <13 AU. Based on the size of the circumbinary ring and the circumstellar disk around GG Tau Aa, the semi-major axis of the binary's orbit is likely to be 62 AU. A comparison of the present observations with previous ALMA and near-infrared (NIR) H$_2$ emission observations suggests that the north arc could be part of a large streamer flowing from the circumbinary ring to sustain the circumstellar disks. According to the previous studies, the circumstellar disk around GG Tau Aa has enough mass and can sustain itself for a duration sufficient for planet formation; thus, our study indicates that planets can form within close (separation $\lesssim$ 100 AU) young binary systems.

Motivated by the tensions in the Hubble constant $H_0$ and the structure growth $\sigma_8$ between $Planck$ results and other low redshift measurements, we discuss some cosmological effects of a dark sector model in which dark matter (DM) interacts with fermionic dark radiation (DR) through a light gauge boson (dark photon). Such kind of models are very generic in particle physics with a dark sector with dark gauge symmetries. The effective number of neutrinos is increased by $\delta N_{eff} \sim 0.5$ due to light dark photon and fermionic DR, thereby resolving the conflicts in $H_0$. The elastic scattering between DM and DR induces suppression for DM's density perturbation, but without acoustic oscillations. For weakly-interacting DM around $100$GeV, the new gauge coupling should be $\sim 10^{-4}$ to have sizable effect on matter power spectrum in order to relax the tension in $\sigma_8$.

Reconstructing the evolution history of the dark energy equation of state parameter $w(z)$ directly from observational data is highly valuable in cosmology, since it contains substantial clues in understanding the nature of the accelerated expansion of the Universe. Many works have focused on reconstructing $w(z)$ using Type Ia supernova data, however, only a few studies pay attention to Hubble parameter data. In the present work, we explore the merit of Hubble parameter data and make an attempt to reconstruct $w(z)$ from them through the principle component analysis approach. We find that current Hubble parameter data perform well in reconstructing $w(z)$; though, when compared to supernova data, the data are scant and their quality is worse. Both $\Lambda$CDM and evolving $w(z)$ models can be constrained within $10\%$ at redshifts $z \lesssim 1.5$ and even $5\%$ at redshifts 0.1 $\lesssim$ z $\lesssim$ 1 by using simulated $H(z)$ data of observational quality.

In this study, we investigate a scenario that dark matter (DM) has only gravitational interaction. In the framework of effective field theory of gravity, we find that DM is still stable at tree level even if there is no symmetry to protect its longevity, but could decay into standard model particles due to gravitational loop corrections. The radiative corrections can lead to both higher- and lower-dimensional effective operators. We also first explore how DM can be produced in the early universe. Through gravitational interaction at high temperature, DM is then found to have mass around TeV $\lesssim m_X \lesssim 10^{11}$GeV to get the right relic abundance. When DM decays, it mostly decays into gravitons, which could be tested by current and future CMB experiments. We also estimate the resulting fluxes for cosmic rays, gamma-ray and neutrino.

Cosmic rays are the highest energy particles found in nature. Measurements of the mass composition of cosmic rays between 10^{17} eV and 10^{18} eV are essential to understand whether this energy range is dominated by Galactic or extragalactic sources. It has also been proposed that the astrophysical neutrino signal comes from accelerators capable of producing cosmic rays of these energies. Cosmic rays initiate cascades of secondary particles (air showers) in the atmosphere and their masses are inferred from measurements of the atmospheric depth of the shower maximum, Xmax, or the composition of shower particles reaching the ground. Current measurements suffer from either low precision, or a low duty cycle and a high energy threshold. Radio detection of cosmic rays is a rapidly developing technique, suitable for determination of Xmax with a duty cycle of in principle nearly 100%. The radiation is generated by the separation of relativistic charged particles in the geomagnetic field and a negative charge excess in the shower front. Here we report radio measurements of Xmax with a mean precision of 16 g/cm^2 between 10^{17}-10^{17.5} eV. Because of the high resolution in $Xmax we can determine the mass spectrum and find a mixed composition, containing a light mass fraction of ~80%. Unless the extragalactic component becomes significant already below 10^{17.5} eV, our measurements indicate an additional Galactic component dominating at this energy range.

We investigate possible cosmological effects of interacting scalar radiation and dark matter. After its decoupling, scalar radiation can stream freely as neutrinos or self-interact strongly as perfect fluid, highly depending on the magnitude of its self-couplings. We obtain the general and novel structure for self-scattering rate and compare it with the expansion rate of our Universe. If its trilinear/cubic coupling is non-zero, scalar radiation can be eventually treated as perfect fluid. Possible effects on CMB are also discussed. When this scalar also mediates interaction among dark matter particles, the linear matter power spectrum for large scale structure can be modified differently from other models. We propose to use Debye shielding to avoid the singularity appearing in the scattering between scalar radiation and dark matter.

We aim to unveil the observational imprint of physical mechanisms that govern planetary formation in the young, multiple system GG Tau A. We present ALMA observations of $^{12}$CO and $^{13}$CO 3-2 and 0.9 mm continuum emission with 0.35" resolution. The $^{12}$CO 3-2 emission, found within the cavity of the circumternary dust ring (at radius $< 180$ AU) where no $^{13}$CO emission is detected, confirms the presence of CO gas near the circumstellar disk of GG Tau Aa. The outer disk and the recently detected hot spot lying at the outer edge of the dust ring are mapped both in $^{12}$CO and $^{13}$CO. The gas emission in the outer disk can be radially decomposed as a series of slightly overlapping Gaussian rings, suggesting the presence of unresolved gaps or dips. The dip closest to the disk center lies at a radius very close to the hot spot location at $\sim250-260$~AU. The CO excitation conditions indicate that the outer disk remains in the shadow of the ring. The hot spot probably results from local heating processes. The two latter points reinforce the hypothesis that the hot spot is created by an embedded proto-planet shepherding the outer disk.

We are performing a series of observations with ground-based telescopes toward Planck Galactic cold clumps (PGCCs) in the $\lambda$ Orionis complex in order to systematically investigate the effects of stellar feedback. In the particular case of PGCC G192.32-11.88, we discovered an extremely young Class 0 protostellar object (G192N) and a proto-brown dwarf candidate (G192S). G192N and G192S are located in a gravitationally bound bright-rimmed clump. The velocity and temperature gradients seen in line emission of CO isotopologues indicate that PGCC G192.32-11.88 is externally heated and compressed. G192N probably has the lowest bolometric luminosity ($\sim0.8$ L$_{\sun}$) and accretion rate (6.3$\times10^{-7}$ M$_{\sun}$~yr$^{-1}$) when compared with other young Class 0 sources (e.g. PACS Bright Red sources (PBRs)) in the Orion complex. It has slightly larger internal luminosity ($0.21\pm0.01$ L$_{\sun}$) and outflow velocity ($\sim$14 km~s$^{-1}$) than the predictions of first hydrostatic cores (FHSCs). G192N might be among the youngest Class 0 sources, which are slightly more evolved than a FHSC. Considering its low internal luminosity ($0.08\pm0.01$ L$_{\odot}$) and accretion rate (2.8$\times10^{-8}$ M$_{\sun}$~yr$^{-1}$), G192S is an ideal proto-brown dwarf candidate. The star formation efficiency ($\sim$0.3\%-0.4\%) and core formation efficiency ($\sim$1\%) in PGCC G192.32-11.88 are significantly smaller than in other giant molecular clouds or filaments, indicating that the star formation therein is greatly suppressed due to stellar feedback.

We highlight distinct and systematic observational features of magnetic field morphologies in polarized submm dust continuum. We illustrate this with specific examples and show statistical trends from a sample of 50 star-forming regions.

Ground-based observations at terahertz (THz) frequencies are a newly explorable area of astronomy for the next ten years. We discuss science cases for a first-generation 10-m class THz telescope, focusing on the Greenland Telescope as an example of such a facility. We propose science cases and provide quantitative estimates for each case. The largest advantage of ground-based THz telescopes is their higher angular resolution (~ 4 arcsec for a 10-m dish), as compared to space or airborne THz telescopes. Thus, high-resolution mapping is an important scientific argument. In particular, we can isolate zones of interest for Galactic and extragalactic star-forming regions. The THz windows are suitable for observations of high-excitation CO lines and [N II] 205 um lines, which are scientifically relevant tracers of star formation and stellar feedback. Those lines are the brightest lines in the THz windows, so that they are suitable for the initiation of ground-based THz observations. THz polarization of star-forming regions can also be explored since it traces the dust population contributing to the THz spectral peak. For survey-type observations, we focus on ``sub-THz'' extragalactic surveys, whose uniqueness is to detect galaxies at redshifts z ~ 1--2, where the dust emission per comoving volume is the largest in the history of the Universe. Finally we explore possibilities of flexible time scheduling, which enables us to monitor active galactic nuclei, and to target gamma-ray burst afterglows. For these objects, THz and submillimeter wavelength ranges have not yet been explored.

LOFAR offers the unique capability of observing pulsars across the 10-240 MHz frequency range with a fractional bandwidth of roughly 50%. This spectral range is well-suited for studying the frequency evolution of pulse profile morphology caused by both intrinsic and extrinsic effects: such as changing emission altitude in the pulsar magnetosphere or scatter broadening by the interstellar medium, respectively. The magnitude of most of these effects increases rapidly towards low frequencies. LOFAR can thus address a number of open questions about the nature of radio pulsar emission and its propagation through the interstellar medium. We present the average pulse profiles of 100 pulsars observed in the two LOFAR frequency bands: High Band (120-167 MHz, 100 profiles) and Low Band (15-62 MHz, 26 profiles). We compare them with Westerbork Synthesis Radio Telescope (WSRT) and Lovell Telescope observations at higher frequencies (350 and1400 MHz) in order to study the profile evolution. The profiles are aligned in absolute phase by folding with a new set of timing solutions from the Lovell Telescope, which we present along with precise dispersion measures obtained with LOFAR. We find that the profile evolution with decreasing radio frequency does not follow a specific trend but, depending on the geometry of the pulsar, new components can enter into, or be hidden from, view. Nonetheless, in general our observations confirm the widening of pulsar profiles at low frequencies, as expected from radius-to-frequency mapping or birefringence theories. We offer this catalog of low-frequency pulsar profiles in a user friendly way via the EPN Database of Pulsar Profiles (http://www.epta.eu.org/epndb/).

Most molecular clouds are filamentary or elongated. Among those forming low-mass stars, their long axes tend to be either parallel or perpendicular to the large-scale (10-100 pc) magnetic field (B-field) in the surrounding inter cloud medium. This arises because, along the dynamically dominant B-fields, the competition between self-gravity and turbulent pressure will shape the cloud to be elongated either perpendicular or parallel to the fields. Recent study also suggested that, on the scales of 0.1-0.01 pc, fields are dynamically important within cloud cores forming massive stars. But whether the core field morphologies are inherited from the inter cloud medium or governed by cloud turbulence is under vigorous debate, so is the role played by B-fields in cloud fragmentation at 10 - 0.1 pc scales. Here we report B-field maps covering 100-0.01 pc scales inferred from polarimetric observations of a massive-star forming region, NGC 6334. First, the main filament also lies perpendicular to the ambient field. NGC 6334 hosts young star-forming sites where fields are not severely affected by stellar feedback, and their directions do not change significantly over the entire scale range. This means that the fields are dynamically important. At various scales, we find that the hourglass-shaped field lines are pinched where the gas column density peaks and the field strength is proportional to the 0.4-power of the density. We conclude that B-fields play a crucial role in the fragmentation of NGC 6334.

The recently observed IceCube PeV events could be due to heavy dark matter (DM) decay. In this paper, we propose a simple DM model with extra $U(1)_X$ gauge symmetry and bridge it with standard model particles through heavy right-handed neutrino. The Dirac fermion DM $\chi$ with mass ~5 PeV can dominantly decay into a dark Higgs ($\phi$), the SM Higgs ($h$) and a neutrino ($\nu$). If the lifetime of $\chi$ is ~O($10^{28}$) sec, the resulting neutrino flux can fit data consistently. The neutrino flux from $\chi \rightarrow \phi h \nu$ in our model is softer than the one predicted from $\chi \rightarrow \nu h$, for example. We also discuss a possible mechanism to produce DM with the right relic abundance.

It is generally expected that adding light sterile species would increase the effective number of neutrinos, $N_{eff}$. In this paper we discuss a scenario that $N_{eff}$ can actually decrease due to the neutrino oscillation effect if sterile neutrinos have self-interactions. We specifically focus on the eV mass range, as suggested by the neutrino anomalies. With large self-interactions, sterile neutrinos are not fully thermalized in the early Universe because of the suppressed effective mixing angle or matter effect. As the Universe cools down, flavor equilibrium between active and sterile species can be reached after big bang nucleosynthesis (BBN) epoch, but leading to a decrease of $N_{eff}$. In such a scenario, we also show that the conflict with cosmological mass bounds on the additional sterile neutrinos can be relaxed further when more light species are introduced.

We present the Multifrequency Snapshot Sky Survey (MSSS), the first northern-sky LOFAR imaging survey. In this introductory paper, we first describe in detail the motivation and design of the survey. Compared to previous radio surveys, MSSS is exceptional due to its intrinsic multifrequency nature providing information about the spectral properties of the detected sources over more than two octaves (from 30 to 160 MHz). The broadband frequency coverage, together with the fast survey speed generated by LOFAR's multibeaming capabilities, make MSSS the first survey of the sort anticipated to be carried out with the forthcoming Square Kilometre Array (SKA). Two of the sixteen frequency bands included in the survey were chosen to exactly overlap the frequency coverage of large-area Very Large Array (VLA) and Giant Metrewave Radio Telescope (GMRT) surveys at 74 MHz and 151 MHz respectively. The survey performance is illustrated within the "MSSS Verification Field" (MVF), a region of 100 square degrees centered at J2000 (RA,Dec)=(15h,69deg). The MSSS results from the MVF are compared with previous radio survey catalogs. We assess the flux and astrometric uncertainties in the catalog, as well as the completeness and reliability considering our source finding strategy. We determine the 90% completeness levels within the MVF to be 100 mJy at 135 MHz with 108" resolution, and 550 mJy at 50 MHz with 166" resolution. Images and catalogs for the full survey, expected to contain 150,000-200,000 sources, will be released to a public web server. We outline the plans for the ongoing production of the final survey products, and the ultimate public release of images and source catalogs.

The bright TeV source HESS J1640-465 is positionally coincident with the young SNR G338.3-0.0, and the nearby HESS J1641-463 with TeV gamma-ray emission seems to be closely associated with it. Based on the nonlinear diffusion shock acceleration (NLDSA) model, we explore the emission from these two TeV sources, the particle diffusion is assumed to be different inside and outside the absorbing boundary of the particles accelerated in the SNR shock. The results indicate that (1) the GeV to TeV emission from the region of the HESS J1640-465 is produced as a result of the particle acceleration inside the SNR G338.3-0.0; and (2) the runaway cosmic-ray particles outside the SNR are interacting with nearby dense molecular cloud (MC) at the region of the HESS J1641-463, corresponding $\pi^0$ decay gamma-ray in proton-proton collision contribute to the TeV emission from the HESS J1641-463. Also we investigate the possible X-ray emission in molecular cloud from synchrotron procedure by secondary $e^\pm$ produced through escaped protons interaction with the MC.

PSR B0823+26, a 0.53-s radio pulsar, displays a host of emission phenomena over timescales of seconds to (at least) hours, including nulling, subpulse drifting, and mode-changing. Studying pulsars like PSR B0823+26 provides further insight into the relationship between these various emission phenomena and what they might teach us about pulsar magnetospheres. Here we report on the LOFAR discovery that PSR B0823+26 has a weak and sporadically emitting 'quiet' (Q) emission mode that is over 100 times weaker (on average) and has a nulling fraction forty-times greater than that of the more regularly-emitting 'bright' (B) mode. Previously, the pulsar has been undetected in the Q-mode, and was assumed to be nulling continuously. PSR B0823+26 shows a further decrease in average flux just before the transition into the B-mode, and perhaps truly turns off completely at these times. Furthermore, simultaneous observations taken with the LOFAR, Westerbork, Lovell, and Effelsberg telescopes between 110 MHz and 2.7 GHz demonstrate that the transition between the Q-mode and B-mode occurs within one single rotation of the neutron star, and that it is concurrent across the range of frequencies observed.

We present Higgs-Portal dark matter (DM) models to explain the reported Galactic Center GeV gamma-ray excess. Naive effective theories are inconsistent with direct detection constraint for the relevant parameter range. Simple extended models with dark gauge symmetries can easily accommodate the gamma-ray excess through the Higgs-Portal coupling while satisfying various constraints.

Using observations obtained with the LOw Fequency ARray (LOFAR), the Westerbork Synthesis Radio Telescope (WSRT) and archival Very Large Array (VLA) data, we have traced the radio emission to large scales in the complex source 4C 35.06 located in the core of the galaxy cluster Abell 407. At higher spatial resolution (~4"), the source was known to have two inner radio lobes spanning 31 kpc and a diffuse, low-brightness extension running parallel to them, offset by about 11 kpc (in projection). At 62 MHz, we detect the radio emission of this structure extending out to 210 kpc. At 1.4 GHz and intermediate spatial resolution (~30"), the structure appears to have a helical morphology. We have derived the characteristics of the radio spectral index across the source. We show that the source morphology is most likely the result of at least two episodes of AGN activity separated by a dormant period of around 35 Myr. The AGN is hosted by one of the galaxies located in the cluster core of Abell 407. We propose that it is intermittently active as it moves in the dense environment in the cluster core. Using LOFAR, we can trace the relic plasma from that episode of activity out to greater distances from the core than ever before. Using the the WSRT, we detect HI in absorption against the center of the radio source. The absorption profile is relatively broad (FWHM of 288 km/s), similar to what is found in other clusters. Understanding the duty cycle of the radio emission as well as the triggering mechanism for starting (or restarting) the radio-loud activity can provide important constraints to quantify the impact of AGN feedback on galaxy evolution. The study of these mechanisms at low frequencies using morphological and spectral information promises to bring new important insights in this field.

Context. The LOw Frequency ARray (LOFAR) radio telescope is a giant digital phased array interferometer with multiple antennas distributed in Europe. It provides discrete sets of Fourier components of the sky brightness. Recovering the original brightness distribution with aperture synthesis forms an inverse problem that can be solved by various deconvolution and minimization methods Aims. Recent papers have established a clear link between the discrete nature of radio interferometry measurement and the "compressed sensing" (CS) theory, which supports sparse reconstruction methods to form an image from the measured visibilities. Empowered by proximal theory, CS offers a sound framework for efficient global minimization and sparse data representation using fast algorithms. Combined with instrumental direction-dependent effects (DDE) in the scope of a real instrument, we developed and validated a new method based on this framework Methods. We implemented a sparse reconstruction method in the standard LOFAR imaging tool and compared the photometric and resolution performance of this new imager with that of CLEAN-based methods (CLEAN and MS-CLEAN) with simulated and real LOFAR data Results. We show that i) sparse reconstruction performs as well as CLEAN in recovering the flux of point sources; ii) performs much better on extended objects (the root mean square error is reduced by a factor of up to 10); and iii) provides a solution with an effective angular resolution 2-3 times better than the CLEAN images. Conclusions. Sparse recovery gives a correct photometry on high dynamic and wide-field images and improved realistic structures of extended sources (of simulated and real LOFAR datasets). This sparse reconstruction method is compatible with modern interferometric imagers that handle DDE corrections (A- and W-projections) required for current and future instruments such as LOFAR and SKA

Hidden sector dark matter (DM) models with local dark gauge symmetry make a natural playground for the $\gamma$-ray excess from the galactic center (GC). In this paper, we first discuss in detail the GC $\gamma$-ray excess in scalar dark matter model with local $Z_{3}$ symmetry which was recently proposed by the present authors. Within this model, scalar DM with mass $30-70$GeV is allowed due to newly-opened (semi-)annihilation channels of DM pair into dark Higgs $\phi$ and/or dark photon $Z^{'}$ pair, and can fit the $\gamma$-ray spectrum from the GC. Then we argue that the GC gamma ray excess can be easily accommodated in hidden sector dark matter models where DM is stabilized by local gauge symmetries, due to the presence of dark Higgs (and also dark photon for Abelian dark gauge symmetry). As by-products, the Higgs portal interaction between the dark Higgs and the SM Higgs boson can improve the EW vacuum stability up to Planck scale, and enable the Higgs inflation scenario to have potentially large tensor-to-scalar ratio, $r\sim O(0.01 - 0.1)$.

Positron excess observed by PAMELA, Fermi and AMS02 may be due to dark matter (DM) pair annihilation or decay dominantly into muons. In this paper, we consider a scenario with thermal fermionic DM ($\chi$) with mass $\sim O(1-2)$ TeV decaying into a dark Higgs ($\phi$) and an active neutrino ($\nu_a$) instead of the SM Higgs boson and $\nu_a$. We first present a renormalizable model for this scenario with local dark $U(1)_X$ gauge symmetry, in which the DM $\chi$ can be thermalized by Higgs portal and the gauge kinetic mixing. Assuming the dark Higgs ($\phi$) mass is in the range $2 m_\mu < m_\phi < 2 m_{\pi^0}$, the positron excess can be fit in a natural manner without conflict with constraints from antiproton and gamma ray fluxes or direct detection experiments. Also having such a light dark Higgs, the self-interaction of DM can be enhanced to some extent, and three puzzles in the CDM paradigm can be somewhat relaxed.

We perform imaging and analyses of SMA 1.3 mm continuum, C18O (2-1) and 12CO (2-1) line data of 17 Class 0 and 0/I protostars to study their gas kinematics on a 1,000-AU scale. Continuum and C18O (2-1) emission are detected toward all the sample sources and show central primary components with sizes of ~600-1,500 AU associated with protostars. The velocity gradients in C18O (2-1) have wide ranges of orientations from parallel to perpendicular to the outflows, with magnitudes from ~1 to ~530 km/s/pc. We construct a simple kinematic model to reproduce the observed velocity gradients, estimate the infalling and rotational velocities, and infer the disk radii and the protostellar masses. The inferred disk radii range from <5 AU to >500 AU with estimated protostellar masses from <0.1 Msun to >1 Msun. Our results hint that both large and small disks are possibly present around Class 0 protostars, which could be a sign of disk growth at the Class 0 stage. In addition, the directions of the overall velocity gradients in 7 out of the 17 sources are close to perpendicular to their outflow axes, which is a signature of significant rotational motions. From our model fitting, the specific angular momenta in these sources are estimated to be >2E-4 km/s*pc, suggesting that magnetic braking is unlikely efficient on a 1,000-AU scale in these Class 0 and 0/I sources. In a sub-sample with observed magnetic field orientations, we find no source with large specific angular momenta together with closely aligned magnetic field and outflow axes. This possibly hints that the magnetic field, if originally aligned with the rotational axis, can play a role in removing angular momentum from infalling material at the Class 0 stage. We discuss our results in comparison with theoretical models of collapsing dense cores with and without magnetic fields in the context of disk formation.

Measuring radio emission from air showers offers a novel way to determine properties of the primary cosmic rays such as their mass and energy. Theory predicts that relativistic time compression effects lead to a ring of amplified emission which starts to dominate the emission pattern for frequencies above ~100 MHz. In this article we present the first detailed measurements of this structure. Ring structures in the radio emission of air showers are measured with the LOFAR radio telescope in the frequency range of 110 - 190 MHz. These data are well described by CoREAS simulations. They clearly confirm the importance of including the index of refraction of air as a function of height. Furthermore, the presence of the Cherenkov ring offers the possibility for a geometrical measurement of the depth of shower maximum, which in turn depends on the mass of the primary particle.

Submillimeter dust polarization measurements of a sample of 50 star-forming regions, observed with the SMA and the CSO covering pc-scale clouds to mpc-scale cores, are analyzed in order to quantify the magnetic field importance. The magnetic field misalignment $\delta$ -- the local angle between magnetic field and dust emission gradient -- is found to be a prime observable, revealing distinct distributions for sources where the magnetic field is preferentially aligned with or perpendicular to the source minor axis. Source-averaged misalignment angles $\langle|\delta|\rangle$ fall into systematically different ranges, reflecting the different source-magnetic field configurations. Possible bimodal $\langle|\delta|\rangle$-distributions are found for the separate SMA and CSO samples. Combining both samples broadens the distribution with a wide maximum peak at small $\langle|\delta|\rangle$-values. Assuming the 50 sources to be representative, the prevailing source-magnetic field configuration is one that statistically prefers small magnetic field misalignments $|\delta|$. When interpreting $|\delta|$ together with an MHD force equation, as developed in the framework of the polarization-intensity gradient method, a sample-based scaling fits the magnetic field tension-to-gravity force ratio $\langle\Sigma_B\rangle$ versus $\langle|\delta|\rangle$ with $\langle\Sigma_B\rangle = 0.116 \cdot \exp(0.047\cdot \langle|\delta|\rangle)\pm 0.20$ (mean error), providing a way to estimate the relative importance of the magnetic field, only based on measurable field misalignments $|\delta|$. The force ratio $\Sigma_B$ discriminates systems that are collapsible on average ($\langle \Sigma_B\rangle <1$) from other molecular clouds where the magnetic field still provides enough resistance against gravitational collapse ($\langle \Sigma_B\rangle >1$) (abridged).

Aims. An efficient means of locating calibrator sources for International LOFAR is developed and used to determine the average density of usable calibrator sources on the sky for subarcsecond observations at 140 MHz. Methods. We used the multi-beaming capability of LOFAR to conduct a fast and computationally inexpensive survey with the full International LOFAR array. Sources were pre-selected on the basis of 325 MHz arcminute-scale flux density using existing catalogues. By observing 30 different sources in each of the 12 sets of pointings per hour, we were able to inspect 630 sources in two hours to determine if they possess a sufficiently bright compact component to be usable as LOFAR delay calibrators. Results. Over 40% of the observed sources are detected on multiple baselines between international stations and 86 are classified as satisfactory calibrators. We show that a flat low-frequency spectrum (from 74 to 325 MHz) is the best predictor of compactness at 140 MHz. We extrapolate from our sample to show that the density of calibrators on the sky that are sufficiently bright to calibrate dispersive and non-dispersive delays for the International LOFAR using existing methods is 1.0 per square degree. Conclusions. The observed density of satisfactory delay calibrator sources means that observations with International LOFAR should be possible at virtually any point in the sky, provided that a fast and efficient search using the methodology described here is conducted prior to the observation to identify the best calibrator.

We report the first detection of an hourglass magnetic field aligned with a well-defined outflow-rotation system in a high-mass star-forming region. The observations were performed with Submillimeter Array toward G240.31+0.07, which harbors a massive, flattened, and fragmenting molecular cloud core and a wide-angle bipolar outflow. The polarized dust emission at 0.88 mm reveals a clear hourglass-shaped magnetic field aligned within 20 degree of the outflow axis. Maps of high-density tracing spectral lines, e.g., H13CO+ (4-3), show that the core is rotating about its minor axis, which is also aligned with the magnetic field axis. Therefore, both the magnetic field and kinematic properties observed in this region are surprisingly consistent with the theoretical predictions of the classic paradigm of isolated low-mass star formation. The strength of the magnetic field in the plane of sky is estimated to be about 1.1 mG, resulting in a mass-to-magnetic flux ratio of 1.4 times the critical value and a turbulent to ordered magnetic energy ratio of 0.4. We also find that the specific angular momentum almost linearly decreases from r~0.6 pc to 0.03 pc scales, which is most likely attributed to magnetic braking.

We present a new method of analysis for determining the surface geometry of five protoplanetary disks observed with near-infrared imaging polarimetry using Subaru-HiCIAO. Using as inputs the observed distribution of polarized intensity (PI), disk inclination, assumed properties for dust scattering, and other reasonable approximations, we calculate a differential equation to derive the surface geometry. This equation is numerically integrated along the distance from the star at a given position angle. We show that, using these approximations, the local maxima in the PI distribution of spiral arms (SAO 206462, MWC 758) and rings (2MASS J16042165-2130284, PDS 70) are associated with local concave-up structures on the disk surface. We also show that the observed presence of an inner gap in scattered light still allows the possibility of a disk surface that is parallel to the light path from the star, or a disk that is shadowed by structures in the inner radii. Our analysis for rings does not show the presence of a vertical inner wall as often assumed in studies of disks with an inner gap. Finally, we summarize the implications of spiral and ring structures as potential signatures of ongoing planet formation.

We present LOFAR Low Band observations of the Bootes and 3C295 fields. Our images made at 34, 46, and 62 MHz reach noise levels of 12, 8, and 5 mJy beam$^{-1}$, making them the deepest images ever obtained in this frequency range. In total, we detect between 300 and 400 sources in each of these images, covering an area of 17 to 52 deg$^{2}$. From the observations we derive Euclidean-normalized differential source counts. The 62 MHz source counts agree with previous GMRT 153 MHz and VLA 74 MHz differential source counts, scaling with a spectral index of $-0.7$. We find that a spectral index scaling of $-0.5$ is required to match up the LOFAR 34 MHz source counts. This result is also in agreement with source counts from the 38 MHz 8C survey, indicating that the average spectral index of radio sources flattens towards lower frequencies. We also find evidence for spectral flattening using the individual flux measurements of sources between 34 and 1400 MHz and by calculating the spectral index averaged over the source population. To select ultra-steep spectrum ($\alpha < -1.1$) radio sources, that could be associated with massive high redshift radio galaxies, we compute spectral indices between 62 MHz, 153 MHz and 1.4 GHz for sources in the Bo\"otes field. We cross-correlate these radio sources with optical and infrared catalogues and fit the spectral energy distribution to obtain photometric redshifts. We find that most of these ultra-steep spectrum sources are located in the $ 0.7 \lesssim z \lesssim 2.5$ range.

We propose a ultraviolet complete theory for cold dark matter(CDM) and sterile neutrino that can accommodate both cosmological data and neutrino oscillation experiments at $1\sigma$ level. A new $U(1)_X$ dark gauge symmetry is introduced, and is assumed to be broken at $\sim \mathcal{O}$(MeV) scale. Such a light mediator for DM's self-scattering and scattering-off sterile neutrinos can resolve three controversies for cold DM on small cosmological scales, cusp vs. core, too-big-to-fail and missing satellites problems. We can also accommodate $\sim$ eV scale sterile neutrinos as the hot dark matter(HDM) and can fit some neutrino anomalies from neutrino oscillation experiments within $1\sigma$. Finally the right amount of HDM can make a sizable contribution to dark radiation, and also helps to reconcile the tension between the data on the tensor-to-scalar ratio reported by Planck and BICEP2 Collaborations.

We have conducted two pilot surveys for radio pulsars and fast transients with the Low-Frequency Array (LOFAR) around 140 MHz and here report on the first low-frequency fast-radio burst limit and the discovery of two new pulsars. The first survey, the LOFAR Pilot Pulsar Survey (LPPS), observed a large fraction of the northern sky, ~1.4 x 10^4 sq. deg, with 1-hr dwell times. Each observation covered ~75 sq. deg using 7 independent fields formed by incoherently summing the high-band antenna fields. The second pilot survey, the LOFAR Tied-Array Survey (LOTAS), spanned ~600 sq. deg, with roughly a 5-fold increase in sensitivity compared with LPPS. Using a coherent sum of the 6 LOFAR "Superterp" stations, we formed 19 tied-array beams, together covering 4 sq. deg per pointing. From LPPS we derive a limit on the occurrence, at 142 MHz, of dispersed radio bursts of < 150 /day/sky, for bursts brighter than S > 107 Jy for the narrowest searched burst duration of 0.66 ms. In LPPS, we re-detected 65 previously known pulsars. LOTAS discovered two pulsars, the first with LOFAR or any digital aperture array. LOTAS also re-detected 27 previously known pulsars. These pilot studies show that LOFAR can efficiently carry out all-sky surveys for pulsars and fast transients, and they set the stage for further surveying efforts using LOFAR and the planned low-frequency component of the Square Kilometer Array.

Recent exo-planetary surveys reveal that planets can orbit and survive around binary stars. This suggests that some fraction of young binary systems which possess massive circumbinary disks (CB) may be in the midst of planet formation. However, there are very few CB disks detected. We revisit one of the known CB disks, the UY Aurigae system, and probe 13CO 2-1, C18O 2-1, SO 5(6)-4(5) and 12CO 3-2 line emission and the thermal dust continuum. Our new results confirm the existence of the CB disk. In addition, the circumstellar (CS) disks are clearly resolved in dust continuum at 1.4 mm. The spectral indices between the wavelengths of 0.85 mm and 6 cm are found to be surprisingly low, being 1.6 for both CS disks. The deprojected separation of the binary is 1.26" based on our 1.4 mm continuum data. This is 0.07" (10 AU) larger than in earlier studies. Combining the fact of the variation of UY Aur B in $R$ band, we propose that the CS disk of an undetected companion UY Aur Bb obscures UY Aur Ba. A very complex kinematical pattern inside the CB disk is observed due to a mixing of Keplerian rotation of the CB disk, the infall and outflow gas. The streaming gas accreting from the CB ring toward the CS disks and possible outflows are also identified and resolved. The SO emission is found to be at the bases of the streaming shocks. Our results suggest that the UY Aur system is undergoing an active accretion phase from the CB disk to the CS disks. The UY Aur B might also be a binary system, making the UY Aur a triple system.

We present radio observations of the Moon between $35$ and $80$ MHz to demonstrate a novel technique of interferometrically measuring large-scale diffuse emission extending far beyond the primary beam (global signal) for the first time. In particular, we show that (i) the Moon appears as a negative-flux source at frequencies $35<\nu<80$ MHz since it is `colder' than the diffuse Galactic background it occults, (ii) using the (negative) flux of the lunar disc, we can reconstruct the spectrum of the diffuse Galactic emission with the lunar thermal emission as a reference, and (iii) that reflected RFI (radio-frequency interference) is concentrated at the center of the lunar disc due to specular nature of reflection, and can be independently measured. Our RFI measurements show that (i) Moon-based Cosmic Dawn experiments must design for an Earth-isolation of better than $80$ dB to achieve an RFI temperature $<1$ mK, (ii) Moon-reflected RFI contributes to a dipole temperature less than $20$ mK for Earth-based Cosmic Dawn experiments, (iii) man-made satellite-reflected RFI temperature exceeds $20$ mK if the aggregate cross section of visible satellites exceeds $80$ m$^2$ at $800$ km height, or $5$ m$^2$ at $400$ km height. Currently, our diffuse background spectrum is limited by sidelobe confusion on short baselines (10-15% level). Further refinement of our technique may yield constraints on the redshifted global $21$-cm signal from Cosmic Dawn ($40>z>12$) and the Epoch of Reionization ($12>z>5$).

The Sun is an active source of radio emission which is often associated with energetic phenomena such as solar flares and coronal mass ejections (CMEs). At low radio frequencies (<100 MHz), the Sun has not been imaged extensively because of the instrumental limitations of previous radio telescopes. Here, the combined high spatial, spectral and temporal resolution of the Low Frequency Array (LOFAR) was used to study solar Type III radio bursts at 30-90 MHz and their association with CMEs. The Sun was imaged with 126 simultaneous tied-array beams within 5 solar radii of the solar centre. This method offers benefits over standard interferometric imaging since each beam produces high temporal (83 ms) and spectral resolution (12.5 kHz) dynamic spectra at an array of spatial locations centred on the Sun. LOFAR's standard interferometric output is currently limited to one image per second. Over a period of 30 minutes, multiple Type III radio bursts were observed, a number of which were found to be located at high altitudes (4 solar radii from the solar center at 30 MHz) and to have non-radial trajectories. These bursts occurred at altitudes in excess of values predicted by 1D radial electron density models. The non-radial high altitude Type III bursts were found to be associated with the expanding flank of a CME. The CME may have compressed neighbouring streamer plasma producing larger electron densities at high altitudes, while the non-radial burst trajectories can be explained by the deflection of radial magnetic fields as the CME expanded in the low corona.

Massive stars ($M > 8$ \msun) typically form in parsec-scale molecular clumps that collapse and fragment, leading to the birth of a cluster of stellar objects. We investigate the role of magnetic fields in this process through dust polarization at 870 $\mu$m obtained with the Submillimeter Array (SMA). The SMA observations reveal polarization at scales of $\lsim$ 0.1 pc. The polarization pattern in these objects ranges from ordered hour-glass configurations to more chaotic distributions. By comparing the SMA data with the single dish data at parsec scales, we found that magnetic fields at dense core scales are either aligned within $40^\circ$ of or perpendicular to the parsec-scale magnetic fields. This finding indicates that magnetic fields play an important role during the collapse and fragmentation of massive molecular clumps and the formation of dense cores. We further compare magnetic fields in dense cores with the major axis of molecular outflows. Despite a limited number of outflows, we found that the outflow axis appears to be randomly oriented with respect to the magnetic field in the core. This result suggests that at the scale of accretion disks ($\lsim 10^3$ AU), angular momentum and dynamic interactions possibly due to close binary or multiple systems dominate over magnetic fields. With this unprecedentedly large sample massive clumps, we argue on a statistical basis that magnetic fields play an important role during the formation of dense cores at spatial scale of 0.01 - 0.1 pc in the context of massive star and cluster star formation.

We report on a study of the 2300-2600\AA FeII/FeII* multiplets in the rest-UV spectra of star-forming galaxies at 1.0<z<2.6 as probes of galactic-scale outflows. We extracted a mass-limited sample of 97 galaxies at z~1.0-2.6 from ultra-deep spectra obtained during the GMASS spetroscopic survey in the GOODS South field with the VLT and FORS2. We obtain robust measures of the rest equivalent width of the FeII absorption lines down to a limit of W_r>1.5 \AA and of the FeII* emission lines to W_r>0.5 \AA. Whenever we can measure the systemic redshift of the galaxies from the [OII] emission line, we find that both the FeII and MgII absorption lines are blueshifted, indicative that both species trace gaseous outflows. We also find, however, that the FeII gas has generally lower outflow velocity relative to that of MgII. We investigate the variation of FeII line profiles as a function of the radiative transfer properties of the lines, and find that transitions with higher oscillator strengths are more blueshifted in terms of both line centroids and line wings. We discuss the possibility that FeII lines are suppressed by stellar absorptions. The lower velocities of the FeII lines relative to the MgII doublet, as well as the absence of spatially extended FeII* emission in 2D stacked spectra, suggest that most clouds responsible for the FeII absorption lie close (3~4 kpc) to the disks of galaxies. We show that the FeII/FeII* multiplets offer unique probes of the kinematic structure of galactic outflows.

This study aims to characterise the polarized foreground emission in the ELAIS-N1 field and to address its possible implications for the extraction of the cosmological 21-cm signal from the Low-Frequency Array - Epoch of Reionization (LOFAR-EoR) data. We use the high band antennas of LOFAR to image this region and RM-synthesis to unravel structures of polarized emission at high Galactic latitudes. The brightness temperature of the detected Galactic emission is on average 4 K in polarized intensity and covers the range from -10 to +13rad m^-2 in Faraday depth. The total polarized intensity and polarization angle show a wide range of morphological features. We have also used the Westerbork Synthesis Radio Telescope (WSRT) at 350 MHz to image the same region. The LOFAR and WSRT images show a similar complex morphology, at comparable brightness levels, but their spatial correlation is very low. The fractional polarization at 150 MHz, expressed as a percentage of the total intensity, amounts to 1.5%. There is no indication of diffuse emission in total intensity in the interferometric data, in line with results at higher frequencies. The wide frequency range, good angular resolution and good sensitivity make LOFAR an exquisite instrument for studying Galactic polarized emission at a resolution of 1-2 rad m^-2 in Faraday depth. The different polarised patterns observed at 150 MHz and 350 MHz are consistent with different source distributions along the line of sight wring in a variety of Faraday thin regions of emission. The presence of polarised foregrounds is a serious complication for Epoch of Reionization experiments. To avoid the leakage of polarized emission into total intensity, which can depend on frequency, we need to calibrate the instrumental polarization across the field of view to a small fraction of 1%.

Extensive air showers, induced by high energy cosmic rays impinging on the Earth's atmosphere, produce radio emission that is measured with the LOFAR radio telescope. As the emission comes from a finite distance of a few kilometers, the incident wavefront is non-planar. A spherical, conical or hyperbolic shape of the wavefront has been proposed, but measurements of individual air showers have been inconclusive so far. For a selected high-quality sample of 161 measured extensive air showers, we have reconstructed the wavefront by measuring pulse arrival times to sub-nanosecond precision in 200 to 350 individual antennas. For each measured air shower, we have fitted a conical, spherical, and hyperboloid shape to the arrival times. The fit quality and a likelihood analysis show that a hyperboloid is the best parametrization. Using a non-planar wavefront shape gives an improved angular resolution, when reconstructing the shower arrival direction. Furthermore, a dependence of the wavefront shape on the shower geometry can be seen. This suggests that it will be possible to use a wavefront shape analysis to get an additional handle on the atmospheric depth of the shower maximum, which is sensitive to the mass of the primary particle.

We construct a self-interacting scalar dark matter (DM) model with local discrete $Z_{3}$ symmetry that stabilizes a weak scale scalar dark matter $X$. The model assumes a hidden sector with a local $U(1)_X$ dark gauge symmetry, which is broken spontaneously into $Z_3$ subgroup by nonzero VEV of dark Higgs field $\phi_X$ ($ \langle \phi_X \rangle \neq 0$). Compared with global $Z_3$ DM models, the local $Z_3$ model has two new extra fields: a dark gauge field $Z^{'}$ and a dark Higgs field $\phi$ (a remnant of the $U(1)_X$ breaking). After imposing various constraints including the upper bounds on the spin-independent direct detection cross section and thermal relic density, we find that the scalar DM with mass less than $125$ GeV is allowed in the local $Z_3$ model, in contrary to the global $Z_3$ model. This is due to new channels in the DM pair annihilations open into $Z^{'}$ and $\phi$ in the local $Z_3$ model. Most parts of the newly open DM mass region can be probed by XENON1T and other similar future experiments. Also if $\phi$ is light enough (a few MeV $\lesssim m_\phi \lesssim$ O(100) MeV), it can generate a right size of DM self-interaction and explain the astrophysical small scale structure anomalies. This would lead to exotic decays of Higgs boson into a pair of dark Higgs bosons, which could be tested at LHC and ILC.

We show that the $\mathinner{\mathrm{GeV}}$ scale $\gamma$-ray excess from the direction of the Galactic Center can be naturally explained by the pair annihilation of Abelian vector dark matter (VDM) into a pair of dark Higgs bosons $VV \rightarrow \phi \phi$, followed by the subsequent decay of $\phi$ into $\phi \rightarrow b\bar{b} , \tau \bar{\tau} $. All the processes are described by a renormalizable VDM model with the Higgs portal, which is naturally flavor-dependent. Some parameter space of this scenario can be tested at the near future direct dark matter search experiments such as LUX and XENON1T.

We present the first detection of carbon radio recombination line absorption along the line of sight to Cygnus A. The observations were carried out with the LOw Frequency ARray in the 33 to 57 MHz range. These low frequency radio observations provide us with a new line of sight to study the diffuse, neutral gas in our Galaxy. To our knowledge this is the first time that foreground Milky Way recombination line absorption has been observed against a bright extragalactic background source. By stacking 48 carbon $\alpha$ lines in the observed frequency range we detect carbon absorption with a signal-to-noise ratio of about 5. The average carbon absorption has a peak optical depth of 2$\times$10$^{-4}$, a line width of 10 km s$^{-1}$ and a velocity of +4 km s$^{-1}$ with respect to the local standard of rest. The associated gas is found to have an electron temperature $T_{e}\sim$ 110 K and density $n_{e}\sim$ 0.06 cm$^{-3}$. These properties imply that the observed carbon $\alpha$ absorption likely arises in the cold neutral medium of the Orion arm of the Milky Way. Hydrogen and helium lines were not detected to a 3$\sigma$ peak optical depth limit of 1.5$\times$10$^{-4}$ for a 4 km s$^{-1}$ channel width. Radio recombination lines associated with Cygnus A itself were also searched for, but are not detected. We set a 3$\sigma$ upper limit of 1.5$\times$10$^{-4}$ for the peak optical depth of these lines for a 4 km s$^{-1}$ channel width.

The low frequency array (LOFAR), is the first radio telescope designed with the capability to measure radio emission from cosmic-ray induced air showers in parallel with interferometric observations. In the first $\sim 2\,\mathrm{years}$ of observing, 405 cosmic-ray events in the energy range of $10^{16} - 10^{18}\,\mathrm{eV}$ have been detected in the band from $30 - 80\,\mathrm{MHz}$. Each of these air showers is registered with up to $\sim1000$ independent antennas resulting in measurements of the radio emission with unprecedented detail. This article describes the dataset, as well as the analysis pipeline, and serves as a reference for future papers based on these data. All steps necessary to achieve a full reconstruction of the electric field at every antenna position are explained, including removal of radio frequency interference, correcting for the antenna response and identification of the pulsed signal.

We report on detection of an ordered magnetic field (B field) threading a massive star-forming clump in the molecular cloud G35.2-0.74, using Submillimeter Array observations of polarized dust emission. Thanks to the sensitive and high-angular-resolution observations, we are able to resolve the morphology of the B field in the plane of sky and detect a great turn of 90 degree in the B field direction: Over the northern part of the clump, where a velocity gradient is evident, the B field is largely aligned with the long axis of the clump, whereas in the southern part, where the velocity field appears relatively uniform, the B field is slightly pinched with its mean direction perpendicular to the clump elongation. We suggest that the clump forms as its parent cloud collapses more along the large scale B field. In this process, the northern part carries over most of the angular momentum, forming a fast rotating system, and pulls the B field into a toroidal configuration. In contrast, the southern part is not significantly rotating and the B field remains in a poloidal configuration. A statistical analysis of the observed polarization dispersion yields a B field strength of ~ 1 mG, a turbulent-to-magnetic energy ratio of order unity, and a mass-to-magnetic flux ratio of ~ 2--3 in units of the critical value. Detailed calculations support our hypothesis that the B field in the northern part is being rotationally distorted. Our observations, in conjunction with early single-dish data, suggest that the B field may play a critical role in the formation of the dense clump, whereas rotation and turbulence could also be important in further dynamical evolution of the clump. The observations also provide evidence for a wide-angle outflow driven from a strongly rotating region whose B field is largely toroidal.

We present the general purpose spectral energy distribution (SED) fitting tool SED Analysis Through Markov Chains (SATMC). Utilizing Monte Carlo Markov Chain (MCMC) algorithms, SATMC fits an observed SED to SED templates or models of the user's choice to infer intrinsic parameters, generate confidence levels and produce the posterior parameter distribution. Here we describe the key features of SATMC from the underlying MCMC engine to specific features for handling SED fitting. We detail several test cases of SATMC, comparing results obtained to traditional least-squares methods, which highlight its accuracy, robustness and wide range of possible applications. We also present a sample of submillimetre galaxies that have been fitted using the SED synthesis routine GRASIL as input. In general, these SMGs are shown to occupy a large volume of parameter space, particularly in regards to their star formation rates which range from ~30-3000 M_sun yr^-1 and stellar masses which range from ~10^10-10^12 M_sun. Taking advantage of the Bayesian formalism inherent to SATMC, we also show how the fitting results may change under different parametrizations (i.e., different initial mass functions) and through additional or improved photometry, the latter being crucial to the study of high-redshift galaxies.

Dust polarization observations from the Submillimeter Array (SMA) and the Caltech Submillimeter Observatory (CSO) are analyzed with the goal of providing a general tool to interpret the role of the magnetic field in molecular clouds. Magnetic field and dust emission gradient orientations are observed to show distinct patterns and features. The angle $\delta$ between these two orientations can be interpreted as a magnetic field alignment deviation, assuming the emission gradient orientation to coincide with the density gradient orientation in the magnetohydrodynamics (MHD) force equation. In SMA high-resolution (collapsing) cores, additional symmetry properties in $\delta$ can reveal accretion and outflow zones. All these observational findings suggest the angle $\delta$ to be a relevant quantity that can assess the role of the magnetic field. Indeed, when comparing this angle with the (projection-free) magnetic field significance \Sigma_B (Koch et al. 2012a), it is demonstrated that |\delta| yields an approximation to the change in \Sigma_B. Thus, changes in the magnetic field alignment deviation $\delta$ trace changes in the role of the magnetic field. The angle $\delta$ is observationally straightforward to determine, providing a tool to distinguish between zones of minor or significant magnetic field impact. This is exemplified by the CSO M+0.25+0.01, Mon R2, CO+0.02-0.02, M-0.02-0.07 sources and by the SMA high-resolution data from W51 e2, W51 North, Orion BN/KL and g5.89. Additional CSO sources are analyzed, providing further support of this result. Finally, based on the different features found in our sample of 31 sources in total, covering sizes from large-scale complexes to collapsing cores, a schematic evolutionary scenario is proposed (abridged).

The characteristic outer scale of turbulence and the ratio of the random to ordered components of the magnetic field are key parameters to characterise magnetic turbulence in the interstellar gas, which affects the propagation of cosmic rays within the Galaxy. We provide new constraints to those two parameters. We use the LOw Frequency ARray (LOFAR) to image the diffuse continuum emission in the Fan region at (l,b) (137.0,+7.0) at 80"x70" resolution in the range [146,174] MHz. We detect multi-scale fluctuations in the Galactic synchrotron emission and compute their power spectrum. Applying theoretical estimates and derivations from the literature for the first time, we derive the outer scale of turbulence and the ratio of random to ordered magnetic field from the characteristics of these fluctuations . We obtain the deepest image of the Fan region to date and find diffuse continuum emission within the primary beam. The power spectrum of the foreground synchrotron fluctuations displays a power law behaviour for scales between 100 and 8 arcmin with a slope of (-1.84+/-0.19). We find an upper limit of about 20 pc for the outer scale of the magnetic interstellar turbulence toward the Fan region. We also find a variation of the ratio of random to ordered field as a function of Galactic coordinates, supporting different turbulent regimes. We use power spectra fluctuations from LOFAR as well as earlier GMRT and WSRT observations to constrain the outer scale of turbulence of the Galactic synchrotron foreground, finding a range of plausible values of 10-20 pc. Then, we use this information to deduce lower limits of the ratio of ordered to random magnetic field strength. These are found to be 0.3, 0.3, and 0.5 for the LOFAR, WSRT and GMRT fields considered respectively. Both these constraints are in agreement with previous estimates.

Faint undetected sources of radio-frequency interference (RFI) might become visible in long radio observations when they are consistently present over time. Thereby, they might obstruct the detection of the weak astronomical signals of interest. This issue is especially important for Epoch of Reionisation (EoR) projects that try to detect the faint redshifted HI signals from the time of the earliest structures in the Universe. We explore the RFI situation at 30-163 MHz by studying brightness histograms of visibility data observed with LOFAR, similar to radio-source-count analyses that are used in cosmology. An empirical RFI distribution model is derived that allows the simulation of RFI in radio observations. The brightness histograms show an RFI distribution that follows a power-law distribution with an estimated exponent around -1.5. With several assumptions, this can be explained with a uniform distribution of terrestrial radio sources whose radiation follows existing propagation models. Extrapolation of the power law implies that the current LOFAR EoR observations should be severely RFI limited if the strength of RFI sources remains strong after time integration. This is in contrast with actual observations, which almost reach the thermal noise and are thought not to be limited by RFI. Therefore, we conclude that it is unlikely that there are undetected RFI sources that will become visible in long observations. Consequently, there is no indication that RFI will prevent an EoR detection with LOFAR.

We present high-angular-resolution observations of the massive star forming core DR21(OH) at 880 mum using the Submillimeter Array (SMA). The dense core exhibits an overall velocity gradient in a Keplerian-like pattern, which breaks at the center of the core where SMA 6 and SMA 7 are located. The dust polarization shows a complex magnetic field, compatible with a toroidal configuration. This is in contrast with the large, parsec-scale filament that surrounds the core, where there is a smooth magnetic field. The total magnetic field strengths in the filament and in the core are 0.9 and 2.1 mG, respectively. We found evidence of magnetic field diffusion at the core scales, far beyond the expected value for ambipolar diffusion. It is possible that the diffusion arises from fast magnetic reconnection in the presence of turbulence. The dynamics of the DR 21(OH) core appear to be controlled energetically in equal parts by the magnetic field, magneto-hydrodynamic (MHD) turbulence and the angular momentum. The effect of the angular momentum (this is a fast rotating core) is probably causing the observed toroidal field configuration. Yet, gravitation overwhelms all the forces, making this a clear supercritical core with a mass-to-flux ratio of ~6 times the critical value. However, simulations show that this is not enough for the high level of fragmentation observed at 1000 AU scales. Thus, rotation and outflow feedback is probably the main cause of the observed fragmentation.

A statistical survey of 26 major electron events during the period 2002 February through the end of solar cycle 23 is presented. We have obtained electron solar onset times and the peak flux spectra for each event by fitting to a powerlaw spectrum truncated by an exponential high-energy tail. We also derived the coronal magnetic configurations of the related solar active regions (ARs) from the potential-field source-surface model. It is found that (1) 10 of the 11 well-connected open field-line events are prompt events whose solar onset times coincide with the maxima of flare emission and 13 of the 14 closed field-line events are delayed events. (2) A not-wellconnected open field-line event and one of the closed field-line events are prompt events, they are both associated with large-scale coronal disturbances or dimming. (3)An averaged harder spectrum is found in open field-line events compared with the closed ones. Specifically, the averaged spectral index is of 1.6 +/- 0.3 in open field-line events and of 2.0 +/- 0.4 in closed ones. The spectra of three closed field-line events show infinite rollover energies E0. These correlations clearly establish a significant link between the coronal magnetic field-line topology and the escape of charged particles from the flaring ARs into interplanetary space during the major solar energetic particle events.

An X3.4 solar flare and a fast halo coronal mass ejection (CME) occurred on 2006 December 13, accompanied by a high flux of energetic particles recorded both in near-Earth space and at ground level. Our purpose is to provide evidence of flare acceleration in a major solar energetic particle (SEP) event. We first present observations from ACE/EPAM, GOES, and the Apatity neutron monitor. It is found that the initial particle release time coincides with the flare emission and that the spectrum becomes softer and the anisotropy becomes weaker during particle injection, indicating that the acceleration source changes from a confined coronal site to a widespread interplanetary CME-driven shock. We then describe a comprehensive study of the associated flare active region. By use of imaging data from HINODE/SOT and SOHO/MDI magnetogram, we infer the flare magnetic reconnection rate in the form of the magnetic flux change rate. This correlates in time with the microwave emission, indicating a physical link between the flare magnetic reconnection and the acceleration of nonthermal particles. Combining radio spectrograph data from Huairou/NOAC, Culgoora/IPS, Learmonth/RSTN, and WAVES/WIND leads to a continuous and longlasting radio burst extending from a few GHz down to several kHz. Based on the photospheric vector magnetogram from Huairou/NOAC and the nonlinear force free field (NFFF) reconstruction method, we derive the 3D magnetic field configuration shortly after the eruption. Furthermore, we also compute coronal field lines extending to a few solar radii using a potential-field source-surface (PFSS) model. Both the so-called type III-l burst and the magnetic field configuration suggest that open-field lines extend from the flare active region into interplanetary space, allowing the accelerated and charged particles escape.

LOFAR, the LOw-Frequency ARray, is a new-generation radio interferometer constructed in the north of the Netherlands and across europe. Utilizing a novel phased-array design, LOFAR covers the largely unexplored low-frequency range from 10-240 MHz and provides a number of unique observing capabilities. Spreading out from a core located near the village of Exloo in the northeast of the Netherlands, a total of 40 LOFAR stations are nearing completion. A further five stations have been deployed throughout Germany, and one station has been built in each of France, Sweden, and the UK. Digital beam-forming techniques make the LOFAR system agile and allow for rapid repointing of the telescope as well as the potential for multiple simultaneous observations. With its dense core array and long interferometric baselines, LOFAR achieves unparalleled sensitivity and angular resolution in the low-frequency radio regime. The LOFAR facilities are jointly operated by the International LOFAR Telescope (ILT) foundation, as an observatory open to the global astronomical community. LOFAR is one of the first radio observatories to feature automated processing pipelines to deliver fully calibrated science products to its user community. LOFAR's new capabilities, techniques and modus operandi make it an important pathfinder for the Square Kilometre Array (SKA). We give an overview of the LOFAR instrument, its major hardware and software components, and the core science objectives that have driven its design. In addition, we present a selection of new results from the commissioning phase of this new radio observatory.

Faraday rotation measurements using the current and next generation of low-frequency radio telescopes will provide a powerful probe of astronomical magnetic fields. However, achieving the full potential of these measurements requires accurate removal of the time-variable ionospheric Faraday rotation contribution. We present ionFR, a code that calculates the amount of ionospheric Faraday rotation for a specific epoch, geographic location, and line-of-sight. ionFR uses a number of publicly available, GPS-derived total electron content maps and the most recent release of the International Geomagnetic Reference Field. We describe applications of this code for the calibration of radio polarimetric observations, and demonstrate the high accuracy of its modeled ionospheric Faraday rotations using LOFAR pulsar observations. These show that we can accurately determine some of the highest-precision pulsar rotation measures ever achieved. Precision rotation measures can be used to monitor rotation measure variations - either intrinsic or due to the changing line-of-sight through the interstellar medium. This calibration is particularly important for nearby sources, where the ionosphere can contribute a significant fraction of the observed rotation measure. We also discuss planned improvements to ionFR, as well as the importance of ionospheric Faraday rotation calibration for the emerging generation of low-frequency radio telescopes, such as the SKA and its pathfinders.

Cassiopeia A was observed using the Low-Band Antennas of the LOw Frequency ARray (LOFAR) with high spectral resolution. This allowed a search for radio recombination lines (RRLs) along the line-of-sight to this source. Five carbon-alpha RRLs were detected in absorption between 40 and 50 MHz with a signal-to-noise ratio of > 5 from two independent LOFAR datasets. The derived line velocities (v_LSR ~ -50 km/s) and integrated optical depths (~ 13 s^-1) of the RRLs in our spectra, extracted over the whole supernova remnant, are consistent within each LOFAR dataset and with those previously reported. For the first time, we are able to extract spectra against the brightest hotspot of the remnant at frequencies below 330 MHz. These spectra show significantly higher (15-80 %) integrated optical depths, indicating that there is small-scale angular structure on the order of ~1 pc in the absorbing gas distribution over the face of the remnant. We also place an upper limit of 3 x 10^-4 on the peak optical depths of hydrogen and helium RRLs. These results demonstrate that LOFAR has the desired spectral stability and sensitivity to study faint recombination lines in the decameter band.

Some radio pulsars show clear drifting subpulses, in which subpulses are seen to drift in pulse longitude in a systematic pattern. Here we examine how the drifting subpulses of PSR B0809+74 evolve with time and observing frequency. We show that the subpulse period (P3) is constant on timescales of days, months and years, and between 14-5100 MHz. Despite this, the shapes of the driftbands change radically with frequency. Previous studies have concluded that, while the subpulses appear to move through the pulse window approximately linearly at low frequencies (< 500 MHz), a discrete step of 180 degrees in subpulse phase is observed at higher frequencies (> 820 MHz) near to the peak of the average pulse profile. We use LOFAR, GMRT, GBT, WSRT and Effelsberg 100-m data to explore the frequency-dependence of this phase step. We show that the size of the subpulse phase step increases gradually, and is observable even at low frequencies. We attribute the subpulse phase step to the presence of two separate driftbands, whose relative arrival times vary with frequency - one driftband arriving 30 pulses earlier at 20 MHz than it does at 1380 MHz, whilst the other arrives simultaneously at all frequencies. The drifting pattern which is observed here cannot be explained by either the rotating carousel model or the surface oscillation model, and could provide new insight into the physical processes happening within the pulsar magnetosphere.

The aim of the LOFAR Epoch of Reionization (EoR) project is to detect the spectral fluctuations of the redshifted HI 21cm signal. This signal is weaker by several orders of magnitude than the astrophysical foreground signals and hence, in order to achieve this, very long integrations, accurate calibration for stations and ionosphere and reliable foreground removal are essential. One of the prospective observing windows for the LOFAR EoR project will be centered at the North Celestial Pole (NCP). We present results from observations of the NCP window using the LOFAR highband antenna (HBA) array in the frequency range 115 MHz to 163 MHz. The data were obtained in April 2011 during the commissioning phase of LOFAR. We used baselines up to about 30 km. With about 3 nights, of 6 hours each, effective integration we have achieved a noise level of about 100 microJy/PSF in the NCP window. Close to the NCP, the noise level increases to about 180 microJy/PSF, mainly due to additional contamination from unsubtracted nearby sources. We estimate that in our best night, we have reached a noise level only a factor of 1.4 above the thermal limit set by the noise from our Galaxy and the receivers. Our continuum images are several times deeper than have been achieved previously using the WSRT and GMRT arrays. We derive an analytical explanation for the excess noise that we believe to be mainly due to sources at large angular separation from the NCP.

We present the first high-angular resolution (up to 0.7", ~5000 AU) polarization and thermal dust continuum images toward the massive star-forming region W51 North. The observations were carried out with the Submillimeter Array (SMA) in both the subcompact (SMA-SubC) and extended (SMA-Ext) configurations at a wavelength of 870 micron. W51 North is resolved into four cores (SMA1 to SMA4) in the 870 micron continuum image. The associated dust polarization exhibits more complex structures than seen at lower angular resolution. We analyze the inferred morphologies of the plane-of-sky magnetic field (B_bot) in the SMA1 to SMA4 cores and in the envelope using the SMA-Ext and SMA-SubC data. These results are compared with the B_bot archive images obtained from the CSO and JCMT. A correlation between dust intensity gradient position angles (phi_{nabla I}) and magnetic field position angles (phi_B) is found in the CSO, JCMT and both SMA data sets. This correlation is further analyzed quantitatively. A systematically tighter correlation between phi_{nabla I} and phi_B is found in the cores, whereas the correlation decreases in outside-core regions. Magnetic field-to-gravity force ratio (Sigma_B) maps are derived using the newly developed polarization - intensity gradient method by Koch, Tang & Ho 2012. We find that the force ratios tend to be small (Sigma_B <= 0.5) in the cores in all 4 data sets. In regions outside of the cores, the ratios increase or the field is even dominating gravity (Sigma_B > 1). This possibly provides a physical explanation of the tightening correlation between phi_{nabla I} and phi_B in the cores: the more the B field lines are dragged and aligned by gravity, the tighter the correlation is. Finally, we propose a schematic scenario for the magnetic field in W51 North to interpret the four polarization observations at different physical scales.

M87 is a giant elliptical galaxy located in the centre of the Virgo cluster, which harbours a supermassive black hole of mass 6.4x10^9 M_sun, whose activity is responsible for the extended (80 kpc) radio lobes that surround the galaxy. The energy generated by matter falling onto the central black hole is ejected and transferred to the intra-cluster medium via a relativistic jet and morphologically complex systems of buoyant bubbles, which rise towards the edges of the extended halo. Here we present the first observations made with the new Low-Frequency Array (LOFAR) of M87 at frequencies down to 20 MHz. Images of M87 were produced at low radio frequencies never explored before at these high spatial resolution and dynamic range. To disentangle different synchrotron models and place constraints on source magnetic field, age and energetics, we also performed a detailed spectral analysis of M87 extended radio-halo using these observations together with archival data. We do not find any sign of new extended emissions; on the contrary the source appears well confined by the high pressure of the intra-cluster medium. A continuous injection of relativistic electrons is the model that best fits our data, and provides a scenario in which the lobes are still supplied by fresh relativistic particles from the active galactic nuclei. We suggest that the discrepancy between the low-frequency radio-spectral slope in the core and in the halo implies a strong adiabatic expansion of the plasma as soon as it leaves the core area. The extended halo has an equipartition magnetic field strength of ~10 uG, which increases to ~13 uG in the zones where the particle flows are more active. The continuous injection model for synchrotron ageing provides an age for the halo of ~40 Myr, which in turn provides a jet kinetic power of 6-10x10^44 erg/s.

Aims: This paper discusses the spectral occupancy for performing radio astronomy with the Low-Frequency Array (LOFAR), with a focus on imaging observations. Methods: We have analysed the radio-frequency interference (RFI) situation in two 24-h surveys with Dutch LOFAR stations, covering 30-78 MHz with low-band antennas and 115-163 MHz with high-band antennas. This is a subset of the full frequency range of LOFAR. The surveys have been observed with a 0.76 kHz / 1 s resolution. Results: We measured the RFI occupancy in the low and high frequency sets to be 1.8% and 3.2% respectively. These values are found to be representative values for the LOFAR radio environment. Between day and night, there is no significant difference in the radio environment. We find that lowering the current observational time and frequency resolutions of LOFAR results in a slight loss of flagging accuracy. At LOFAR's nominal resolution of 0.76 kHz and 1 s, the false-positives rate is about 0.5%. This rate increases approximately linearly when decreasing the data frequency resolution. Conclusions: Currently, by using an automated RFI detection strategy, the LOFAR radio environment poses no perceivable problems for sensitive observing. It remains to be seen if this is still true for very deep observations that integrate over tens of nights, but the situation looks promising. Reasons for the low impact of RFI are the high spectral and time resolution of LOFAR; accurate detection methods; strong filters and high receiver linearity; and the proximity of the antennas to the ground. We discuss some strategies that can be used once low-level RFI starts to become apparent. It is important that the frequency range of LOFAR remains free of broadband interference, such as DAB stations and windmills.

The circumstellar disk of AB Aurigae has garnered strong attention owing to the apparent existence of spirals at a relatively young stage and also the asymmetric disk traced in thermal dust emission. However, the physical conditions of the spirals are still not well understood. The origin of the asymmetric thermal emission is unclear. We observed the disk at 230 GHz (1.3 mm) in both the continuum and the spectral line ^12CO J=2-1 with IRAM 30-m, the Plateau de Bure interferometer, and the Submillimeter Array to sample all spatial scales from 0.37" to about 50". To combine the data obtained from these telescopes, several methods and calibration issues were checked and discussed. The 1.3 mm continuum (dust) emission is resolved into inner disk and outer ring. Molecular gas at high velocities traced by the CO line is detected next to the stellar location. The inclination angle of the disk is found to decrease toward the center. On a larger scale, based on the intensity weighted dispersion and the integrated intensity map of ^12CO J=2-1, four spirals are identified, where two of them are also detected in the near infrared. The total gas mass of the 4 spirals (M_spiral) is 10^-7 < M_spiral < 10^-5 M_sun, which is 3 orders of magnitude smaller than the mass of the gas ring. Surprisingly, the CO gas inside the spiral is apparently counter-rotating with respect to the CO disk, and it only exhibits small radial motion. The wide gap, the warped disk, and the asymmetric dust ring suggest that there is an undetected companion with a mass of 0.03 M_sun at a radius of 45 AU. Although an hypothetical fly-by cannot be ruled out, the most likely explanation of the AB Aurigae system may be inhomogeneous accretion well above or below the main disk plane from the remnant envelope, which can explain both the rotation and large-scale motions detected with the 30-m image.

Dispersion in the interstellar medium is a well known phenomenon that follows a simple relationship, which has been used to predict the time delay of dispersed radio pulses since the late 1960s. We performed wide-band simultaneous observations of four pulsars with LOFAR (at 40-190 MHz), the 76-m Lovell Telescope (at 1400 MHz) and the Effelsberg 100-m Telescope (at 8000 MHz) to test the accuracy of the dispersion law over a broad frequency range. In this paper we present the results of these observations which show that the dispersion law is accurate to better than 1 part in 100000 across our observing band. We use this fact to constrain some of the properties of the ISM along the line-of-sight and use the lack of any aberration or retardation effects to determine upper limits on emission heights in the pulsar magnetosphere. We also discuss the effect of pulse profile evolution on our observations, and the implications that it could have for precision pulsar timing projects such as the detection of gravitational waves with pulsar timing arrays.

Abell 2256 is one of the best known examples of a galaxy cluster hosting large-scale diffuse radio emission that is unrelated to individual galaxies. It contains both a giant radio halo and a relic, as well as a number of head-tail sources and smaller diffuse steep-spectrum radio sources. The origin of radio halos and relics is still being debated, but over the last years it has become clear that the presence of these radio sources is closely related to galaxy cluster merger events. Here we present the results from the first LOFAR Low band antenna (LBA) observations of Abell 2256 between 18 and 67 MHz. To our knowledge, the image presented in this paper at 63 MHz is the deepest ever obtained at frequencies below 100 MHz in general. Both the radio halo and the giant relic are detected in the image at 63 MHz, and the diffuse radio emission remains visible at frequencies as low as 20 MHz. The observations confirm the presence of a previously claimed ultra-steep spectrum source to the west of the cluster center with a spectral index of -2.3 \pm 0.4 between 63 and 153 MHz. The steep spectrum suggests that this source is an old part of a head-tail radio source in the cluster. For the radio relic we find an integrated spectral index of -0.81 \pm 0.03, after removing the flux contribution from the other sources. This is relatively flat which could indicate that the efficiency of particle acceleration at the shock substantially changed in the last \sim 0.1 Gyr due to an increase of the shock Mach number. In an alternative scenario, particles are re-accelerated by some mechanism in the downstream region of the shock, resulting in the relatively flat integrated radio spectrum. In the radio halo region we find indications of low-frequency spectral steepening which may suggest that relativistic particles are accelerated in a rather inhomogeneous turbulent region.

Reconnection of the self-generated magnetic fields in laser-plasma interaction was first investigated experimentally by Nilson {\it et al.} [Phys. Rev. Lett. 97, 255001 (2006)] by shining two laser pulses a distance apart on a solid target layer. An elongated current sheet (CS) was observed in the plasma between the two laser spots. In order to more closely model magnetotail reconnection, here two side-by-side thin target layers, instead of a single one, are used. It is found that at one end of the elongated CS a fan-like electron outflow region including three well-collimated electron jets appears. The ($>1$ MeV) tail of the jet energy distribution exhibits a power-law scaling. The enhanced electron acceleration is attributed to the intense inductive electric field in the narrow electron dominated reconnection region, as well as additional acceleration as they are trapped inside the rapidly moving plasmoid formed in and ejected from the CS. The ejection also induces a secondary CS.

Recent high-angular-resolution (up to 0.7") dust polarization observations toward star forming regions are summarized. With the Sub-Millimeter Array, the emission from the dense structures is traced and resolved. The detected magnetic field morphologies vary from hourglass-like structures to isolated patches depending on the evolutionary stage of the source. These observed features have also served as a testbed to develop new analysis methods, with a particular focus on quantifying the role of the magnetic field in the star formation process.

Dust polarization observational results are analyzed for the high-mass star formation region W51 from the largest parent cloud ($\sim$ 2~pc, JCMT) to the large-scale envelope ($\sim$ 0.5~pc, BIMA) down to the collapsing core e2 ($\sim$ 60~mpc, SMA). Magnetic field and dust emission gradient orientations reveal a correlation which becomes increasingly more tight with higher resolution. The previously developed polarization - intensity gradient method (Koch et al. 2012) is applied in order to quantify the magnetic field significance. This technique provides a way to estimate the local magnetic field force compared to gravity without the need of any mass or field strength measurements, solely making use of measured angles which reflect the geometrical imprint of the various forces. All three data sets clearly show regions with distinct features in the field-to-gravity force ratio. Azimuthally averaged radial profiles of this force ratio reveal a transition from a field dominance at larger distances to a gravity dominance closer to the emission peaks. Normalizing these profiles to a characteristic core scale points toward self-similarity. Furthermore, the polarization intensity-gradient method is linked to the mass-to-flux ratio, providing a new approach to estimate the latter one without mass and field strength inputs. A transition from a magnetically supercritical to a subcritical state as a function of distance from the emission peak is found for the e2 core. Finally, based on the measured radius-dependent field-to-gravity force ratio we derive a modified star formation efficiency with a diluted gravity force. Compared to a standard (free-fall) efficiency, the observed field is capable of reducing the efficiency down to 10\% or less.

Dust polarization orientations in molecular clouds often tend to be close to tangential to the Stokes $I$ dust continuum emission contours. The magnetic field and the emission gradient orientations, therefore, show some correlation. A method is proposed, which -- in the framework of ideal magneto-hydrodynamics (MHD) -- connects the measured angle between magnetic field and emission gradient orientations to the total field strength. The approach is based on the assumption that a change in emission intensity (gradient) is a measure for the resulting direction of motion in the MHD force equation. In particular, this new method leads to maps of position-dependent magnetic field strength estimates. When evaluating the field curvature and the gravity direction locally on a map, the method can be generalized to arbitrary cloud shapes. The technique is applied to high-resolution ($\sim0\farcs7$) Submillimeter Array polarization data of the collapsing core W51 e2. A tentative $\sim 7.7$~mG field strength is found when averaging over the entire core. The analysis further reveals some structures and an azimuthally averaged radial profile $\sim r^{-1/2}$ for the field strength. Maximum values close to the center are around $19$~mG. The currently available observations lack higher resolution data to probe the innermost part of the core where the largest field strength is expected from the method. Application regime and limitations of the method are discussed. As a further important outcome of this technique, the local significance of the magnetic field force compared to the other forces can be quantified in a model-independent way, from measured angles only. Finally, the method can potentially also be expanded and applied to other objects (besides molecular clouds) with measurements that reveal the field morphology, as e.g. Faraday rotation measurements in galaxies.

We have imaged the extremely high-velocity outflowing gas in CO (2-1) and (3-2) associated with the shell-like ultracompact HII region G5.89-0.39 at a resolution of ~3" (corresponding to ~4000 AU) with the Submillimeter Array. The integrated high-velocity (>45 km/s) CO emission reveals at least three blueshifted lobes and two redshifted lobes. These lobes belong to two outflows, one oriented N-S, the other NW-SE. The NW-SE outflow is likely identical to the previously detected Br_gamma outflow. Furthermore, these outflow lobes all clearly show a Hubble-like kinematic structure. For the first time, we estimate the temperature of the outflowing gas as a function of velocity with the large velocity gradient calculations. Our results reveal a clear increasing trend of temperature with gas velocity. The observational features of the extremely high-velocity gas associated with G5.89-0.39 qualitatively favor the jet-driven bow shock model.

When an ultra-high energy neutrino or cosmic ray strikes the Lunar surface a radio-frequency pulse is emitted. We plan to use the LOFAR radio telescope to detect these pulses. In this work we propose an efficient trigger implementation for LOFAR optimized for the observation of short radio pulses.

We report the discovery of large amounts of cold (T ~ 10^4 K), chemically young gas in an overdensity of galaxies at redshift z ~ 1.6 in the Great Observatories Origins Deep Survey southern field (GOODS-S). The gas is identified thanks to the ultra-strong Mg II absorption features it imprints in the rest-frame UV spectra of galaxies in the background of the overdensity. There is no evidence that the optically-thick gas is part of any massive galaxy (i.e. M_star > 4x10^9 M_sun), but rather is associated with the overdensity; less massive and fainter galaxies (25.5 < z_850 < 27.5 mag) have too large an impact parameter to be causing ultra-strong absorption systems, based on our knowledge of such systems. The lack of corresponding Fe II absorption features, not detected even in co-added spectra, suggests that the gas is chemically more pristine than the ISM and outflows of star-forming galaxies at similar redshift, including those in the overdensity itself, and comparable to the most metal-poor stars in the Milky Way halo. A crude estimate of the projected covering factor of the high-column density gas (N_H >~ 10^20 cm-2) based on the observed fraction of galaxies with ultra-strong absorbers is C_F ~ 0.04. A broad, continuum absorption profile extending to the red of the interstellar Mg II absorption line by <~ 2000 km/s is possibly detected in two independent co-added spectra of galaxies of the overdensity, consistent with a large-scale infall motion of the gas onto the overdensity and its galaxies. Overall, these findings provides the first tentative evidence of accretion of cold, chemically young gas onto galaxies at high redshift, possibly feeding their star formation activity. The fact that the galaxies are members of a large structure, as opposed to field galaxies, might play a significant role in our ability to detect the accreting gas.

We report here our search for molecular outflows from young very low mass stars and brown dwarfs in Taurus and rho Ophiuchi. Using the Submillimeter Array (SMA) and the Combined Array for Research in Millimeter-wave Astronomy (CARMA), we have observed 4 targets at 1.3 mm wavelength (230 GHz) to search for CO J=2-->1 outflows. A young very low mass star MHO 5 (in Taurus) with an estimated mass of 90 M_J, which is just above the hydrogen-burning limit, shows two gas lobes that are likely outflows. While the CO map of MHO 5 does not show a clear structure of outflow, possibly due to environment gas, its position-velocity diagram indicates two distinct blue- and redshifted components. We therefore conclude that they are components of a bipolar molecular outflow from MHO 5. We estimate an outflow mass of 7.0 x 10^-5 M_Sun and a mass-loss rate of 9.0 x 10^-10 M_Sun. These values are over two orders of magnitude smaller than the typical ones for T Tauri stars and somewhat weaker than those we have observed in the young brown dwarf ISO-Oph 102 of 60 M_J in rho Ophiuchi. This makes MHO~5 the first young very low mass star showing a bipolar molecular outflow in Taurus. The detection boosts the scenario that very low mass objects form like low-mass stars but in a version scaled down by a factor of over 100.

Low frequency radio waves, while challenging to observe, are a rich source of information about pulsars. The LOw Frequency ARray (LOFAR) is a new radio interferometer operating in the lowest 4 octaves of the ionospheric "radio window": 10-240MHz, that will greatly facilitate observing pulsars at low radio frequencies. Through the huge collecting area, long baselines, and flexible digital hardware, it is expected that LOFAR will revolutionize radio astronomy at the lowest frequencies visible from Earth. LOFAR is a next-generation radio telescope and a pathfinder to the Square Kilometre Array (SKA), in that it incorporates advanced multi-beaming techniques between thousands of individual elements. We discuss the motivation for low-frequency pulsar observations in general and the potential of LOFAR in addressing these science goals. We present LOFAR as it is designed to perform high-time-resolution observations of pulsars and other fast transients, and outline the various relevant observing modes and data reduction pipelines that are already or will soon be implemented to facilitate these observations. A number of results obtained from commissioning observations are presented to demonstrate the exciting potential of the telescope. This paper outlines the case for low frequency pulsar observations and is also intended to serve as a reference for upcoming pulsar/fast transient science papers with LOFAR.

We present the Spitzer spectroscopic mapping observations toward the central kpc region of the nearby elliptical LINER galaxy, NGC 4278, by using the Infrared Spectrograph (IRS). These observations reveal rich mid-IR emission features of extended ionized gas, warm molecular hydrogen and dust. Different phases of gas and dust are closely related and belong to a same elongated feature. We further study properties of multi-phase dusty gas to uncover the underlying mechanism of ionization and excitation. The band ratio and intensity of PAH features in the central region might reflect modified size distribution resulted from selective destruction. H$_{2}$ S(0)-S(7) pure rotational lines of molecular hydrogen show excessive intensity and moderately high excitation temperature comparing with photon dissociation region (PDR). A strong and extended [SiII] emission line is detected, which could be a sign of reduced depletion of silicon in interstellar dust. We also discover an extended high ionization region associated with enhanced H$_{2}$ S(1) emission. We conclude that a shock-heating component is required to account for observed emission characteristics, which could be triggered by cloud-cloud interactions during accretion of cold gas from the large HI disk.

Context. Asteroseismology is an effcient tool not only for testing stellar structure and evolutionary theory but also constraining the parameters of stars for which solar-like oscillations are detected, presently. As an important southern asteroseismic target, Tau Ceti, is a metal-poor star. The main features of the oscillations and some frequencies of ? Ceti have been identified. Many scientists propose to comprehensively observe this star as part of the Stellar Observations Network Group. Aims. Our goal is to obtain the optimal model and reliable fundamental parameters for the metal-poor star Tau Ceti by combining all non-asteroseismic observations with these seismological data. Methods. Using the Yale stellar evolution code (YREC), a grid of stellar model candidates that fall within all the error boxes in the HR diagram have been constructed, and both the model frequencies and large- and small- frequency separations are calculated using the Guenther's stellar pulsation code. The \chi2c minimization is performed to identify the optimal modelling parameters that reproduce the observations within their errors. The frequency corrections of near-surface effects to the calculated frequencies using the empirical law, as proposed by Kjeldsen and coworkers, are applied to the models. Results. We derive optimal models, corresponding to masses of about 0.775 - 0.785 M? and ages of about 8 - 10 Gyr. Furthermore, we find that the quantities derived from the non-asteroseismic observations (effective temperature and luminosity) acquired spectroscopically are more accurate than those inferred from interferometry for ? Ceti, because our optimal models are in the error boxes B and C, which are derived from spectroscopy results.

We present sensitive high angular resolution (0.57$''$-0.78$''$) SO, SO$_2$, CO, C$_2$H$_5$OH, HC$_3$N, and HCOCH$_2$OH line observations at millimeter and submillimeter wavelengths of the young O-type protostar W51 North made with the Submillimeter Array (SMA). We report the presence of a large (of about 8000 AU) and hot molecular circumstellar disk around this object, which connects the inner dusty disk with the molecular ring or toroid reported recently, and confirms the existence of a single bipolar outflow emanating from this object. The molecular emission from the large disk is observed in layers with the transitions characterized by high excitation temperatures in their lower energy states (up to 1512 K) being concentrated closer to the central massive protostar. The molecular emission from those transitions with low or moderate excitation temperatures are found in the outermost parts of the disk and exhibits an inner cavity with an angular size of around 0.7$''$. We modeled all lines with a Local Thermodynamic Equilibrium (LTE) synthetic spectra. A detail study of the kinematics of the molecular gas together with a LTE model of a circumstellar disk shows that the innermost parts of the disk are also Keplerian plus a contracting velocity. The emission of the HCOCH$_2$OH reveals the possible presence of a warm ``companion'' located to the northeast of the disk, however its nature is unclear. The emission of the SO and SO$_2$ is observed in the circumstellar disk as well as in the outflow. We suggest that the massive protostar W51 North appears to be in a phase before the presence of a Hypercompact or an Ultracompact HII (HC/UCHII) region, and propose a possible sequence on the formation of the massive stars.

Polarization data from high mass star formation regions (W51 e2/e8, Orion BN/KL) are used to derive statistical properties of the plane of sky projected magnetic field. Structure function and auto-correlation function are calculated for observations with various resolutions from the BIMA and SMA interferometers, covering a range in physical scales from $\sim 70$~mpc to $\sim 2.1$~mpc. Results for the magnetic field turbulent dispersion, its turbulent to mean field strength ratio and the large-scale polarization angle correlation length are presented as a function of the physical scale at the star formation sites. Power law scaling relations emerge for some of these physical quantities. The turbulent to mean field strength ratio is found to be close to constant over the sampled observing range, with a hint of a decrease toward smaller scales, indicating that the role of magnetic field and turbulence is evolving with physical scale. A statistical method is proposed to separate large and small scale correlations from an initial ensemble of polarization segments. This also leads to a definition of a turbulent polarization angle correlation length.

We present observational results of the thermal dust continuum emission and its linear polarization in one of the nearest massive star-forming sites Orion BN/KL in Orion Molecular Cloud-1. The observations were carried out with the Submillimeter Array. With an angular resolution of 1" (~2 mpc; 480 AU), we have detected and resolved the densest cores near the BN/KL region. At a wavelength of ~870 micron, the polarized dust emission can be used to trace the structure of the magnetic field in this star-forming core. The dust continuum appears to arise from a V-shaped region, with a cavity nearly coincident with the center of the explosive outflows observed on larger scales. The position angles (P.A.s) of the observed polarization vary significantly by a total of about 90 degree but smoothly, i.e., curl-like, across the dust ridges. Such a polarization pattern can be explained with dust grains being magnetically aligned instead of mechanically with outflows, since the latter mechanism would cause the P.A.s to be parallel to the direction of the outflow, i.e., radial-like. The magnetic field projected in the plane of sky is therefore derived by rotating the P.A.s of the polarization by 90 degree. We find an azimuthally symmetric structure in the overall magnetic field morphology, with the field directions pointing toward 2.5" west to the center of the explosive outflows. We also find a preferred symmetry plane at a P.A. of 36 degree, which is perpendicular to the mean magnetic field direction (120 degree) of the 0.5 pc dust ridge. Two possible interpretations of the origin of the observed magnetic field structure are discussed.

The brightest continuum source in the Orion Molecular Cloud-3 region (OMC-3), MMS 6, was observed with the Very Large Array (VLA), the Nobeyama Millimeter Array (NMA), and the Submillimeter Array (SMA). Our data were supplemented by near- to mid-infrared archival data taken by Spitzer Space Telescope. The compact continuum source, MMS 6-main, was detected with an H_2 mass of 3.0 Msun with a size of 510 AU. Despite its compact and well condensed appearance, neither clear CO outflow, radio jet, nor infrared sources (at a wave-length shorter than 8 um) were detected at MMS 6-main even with the present high-spatial resolution and high-sensitivity observations. The derived H_2 column density, 2.6x10^25 cm^-2, corresponds to a visual extinction of A_v~15000 mag., and the derived number density is at least two orders of magnitude higher than for the other OMC-2/3 continuum sources. The volume density profile of the source was estimated to have a power-law index of 2 or steeper down to a radius of ~450 AU. The time scale to form a protostar at the center or the time scale elapsed after its formation is estimated to be 830 to 7600 yr. This is much shorter than the typical lifetime of the Class 0/I protostars, which is ~10^(4-5) yr, suggesting that MMS 6-main is probably in either the earliest stage of the proto-stellar core or in the latest stage of the pre-stellar phase.

We present the Spitzer Space Telescope Infrared Array Camera (IRAC) observations for a sample of local elliptical galaxies to study later stages of AGN activities. A sample of 36 elliptical galaxies is selected from the Palomar spectroscopic survey. We detect nuclear non-stellar infrared emission in 9 of them. There is unambiguous evidence of circumnuclear dust in these 9 galaxies in their optical images. We also find a remarkable correlation between the infrared excess emission and the nuclear radio/X-ray emission, suggesting that infrared excess emission is tightly related with nuclear activity. Possible origin of infrared excess emission from hot dust heated by the central AGN is supported by spectral indices of IR excess emission.

We report our observational results of 870 $\mu$m continuum emission and its linear polarization in the massive star formation site W51 e2/e8. Inferred from the linear polarization maps, the magnetic field in the plane of sky (B$_{\bot}$) is traced with an angular resolution of 0$\farcs$7 with the Submillimeter Array (SMA). Whereas previous BIMA observations with an angular resolution of 3$\arcsec$ (0.1 pc) showed a uniform B field, our revealed B$_{\bot}$ morphology is hourglass-like in the collapsing core near the Ultracompact H II region e2 and also possibly in e8. The decrease in polarization near the continuum peak seen at lower angular resolution is apparently due to the more complex structures at smaller scales. In e2, the pinched direction of the hourglass-like B field morphology is parallel to the plane of the ionized accretion flow traced by H53$\alpha$, suggesting that the massive stars are formed via processes similar to the low mass stars, i.e. accretion through a disk, except that the mass involved is much larger. Furthermore, our finding that the resolved collapsing cores in e2 and e8 lie within one subcritical 0.5 pc envelope supports the scenario of \textit{magnetic fragmentation} via ambipolar diffusion. We therefore suggest that magnetic fields control the dynamical evolution of the envelope and cores in W51 e2 and e8.

We report high angular resolution (3") Submillimeter Array (SMA) observations of the molecular cloud associated with the Ultra-Compact HII region G5.89-0.39. Imaged dust continuum emission at 870 micron reveals significant linear polarization. The position angles (PAs) of the polarization vary enormously but smoothly in a region of 2x10^4 AU. Based on the distribution of the PAs and the associated structures, the polarized emission can be separated roughly into two components. The component "x" is associated with a well defined dust ridge at 870 micron, and is likely tracing a compressed B field. The component "o" is located at the periphery of the dust ridge and is probably from the original B field associated with a pre-existing extended structure. The global B field morphology in G5.89, as inferred from the PAs, is clearly disturbed by the expansion of the HII region and the molecular outflows. Using the Chandrasekhar-Fermi method, we estimate from the smoothness of the field structures that the B field strength in the plane of sky can be no more than 2-3 mG. We then compare the energy densities in the radiation, the B field, and the mechanical motions as deduced from the C^17O 3-2 line emission. We conclude that the B field structures are already overwhelmed and dominated by the radiation, outflows, and turbulence from the newly formed massive stars.

Studying the earliest stages in the birth of stars is crucial for understanding how they form. Brown dwarfs with masses between that of stars and planets are not massive enough to maintain stable hydrogen-burning fusion reactions during most of their lifetime. Their origins are subject to much debate in recent literature because their masses are far below the typical mass where core collapse is expected to occur. We present the first confirmed evidence that brown dwarfs undergo a phase of molecular outflow that is typical of young stars. Using the Submillimeter Array, we have obtained a map of a bipolar molecular outflow from a young brown dwarf. We estimate an outflow mass of 1.6 x 10^-4 M_Sun and a mass-loss rate of 1.4 x 10^-9 M_Sun. These values are over two orders of magnitude smaller than the typical ones for T Tauri stars. From our millimiter continuum data and our own analysis of Spitzer infrared photometry, we estimate that the brown dwarf has a disk with a mass of 8 x 10^-3 M_Sun and an outer disk radius of 80 AU. Our results demonstrate that the bipolar molecular outflow operates down to planetary masses, occurring in brown dwarfs as a scaled-down version of the universal process seen in young stars.

Nonthermal X-rays and very high-energy (VHE) $\gamma$-rays have been detected from the supernova remnant (SNR) RX J1713.7-3946, and especially the recent observations with the \textit{Suzaku} satellite clearly reveal a spectral cutoff in the X-ray spectrum, which directly relates to the cutoff of the energy spectrum of the parent electrons. However, whether the origin of the VHE $\gamma$-rays from the SNR is hadronic or leptonic is still in debate. We studied the multi-band nonthermal emission from RX J1713.7-3946 based on a semi-analytical approach to the nonlinear shock acceleration process by including the contribution of the accelerated electrons to the nonthermal radiation. The results show that the multi-band observations on RX J1713.7-3946 can be well explained in the model with appropriate parameters and the TeV $\gamma$-rays have hadronic origin, i.e., they are produced via proton-proton (p-p) interactions as the relativistic protons accelerated at the shock collide with the ambient matter.

Taking into account the helium and metal diffusion, we explore the possible evolutionary status and perform seismic analysis of MOST target: the star $\epsilon$ Eridani. We adopt the different input parameters to construct the models by fitting the available observational constraints: e.g., $T_{eff}$, $L$, $R$, $[Fe/H]$. From computation, we obtain the average large spacings of $\epsilon$ Eridani about $194\pm 1 \mu$Hz. The age of the diffused models has been found to be about 1 Gyr, which is younger than one determined previously by models without diffusion. We found that the effect of pure helium diffusion on the internal structure of the young low-mass star is slight, but the metal diffusion influence is obvious. The metal diffusion leads the models to have much higher temperature in the radiation interior, correspondingly the higher sound speed in the interior of the model, thereby the larger frequency and spacings.

The basic intent of this paper is to model 70 Ophiuchi A using the latest asteroseismic observations as complementary constraints and to determine the fundamental parameters of the star. Additionally, we propose a new quantity to lift the degeneracy between the initial chemical composition and stellar age. Using the Yale stellar evolution code (YREC7), we construct a series of stellar evolutionary tracks for the mass range $M$ = 0.85 -- 0.93 $M_{\odot}$ with different composition $Y_{i}$ (0.26 -- 0.30) and $Z_{i}$ (0.017 -- 0.023). Along these tracks, we select a grid of stellar model candidates that fall within the error box in the HR diagram to calculate the theoretical frequencies, the large- and small- frequency separations using the Guenther's stellar pulsation code. Following the asymptotic formula of stellar $p$-modes, we define a quantity $r_{01}$ which is correlated with stellar age. Also, we test it by theoretical adiabatic frequencies of many models. Many detailed models of 70 Ophiuchi A have been listed in Table 3. By combining all non-asteroseismic observations available for 70 Ophiuchi A with these seismological data, we think that Model 60, Model 125 and Model 126, listed in Table 3, are the optimum models presently. Meanwhile, we predict that the radius of this star is about 0.860 -- 0.865 $R_{\odot}$ and the age is about 6.8 -- 7.0 Gyr with mass 0.89 -- 0.90 $M_{\odot}$. Additionally, we prove that the new quantity $r_{01}$ can be a useful indicator of stellar age.

Asteroseismology, as a tool to use the indirect information contained in stellar oscillations to probe the stellar interiors, is an active field of research presently. Stellar age, as a fundamental property of star apart from its mass, is most difficult to estimate. In addition, the estimating of stellar age can provide the chance to study the time evolution of astronomical phenomena. In our poster, we summarize our previous work and further present a method to determine age of low-mass main-sequence star.

We test whether there is a relation between the observed tidal interactions and Seyfert activity by imaging in HI twenty inactive galaxies at the same spatial resolution and detection threshold as the Seyfert sample. This control sample of inactive galaxies were closely matched in Hubble type, range in size and inclination, and have roughly comparable galaxy optical luminosity to the Seyfert galaxies. We find that only ~15% of the galaxies in our control sample are disturbed in HI, whereas the remaining ~85% show no disturbances whatsoever in HI. Even at a spatial resolution of ~10 kpc, none of the latter galaxies show appreciable HI disturbances reminiscent of tidal features. In a companion paper (Kuo et al. 2008), we report results from the first systematic imaging survey of Seyfert galaxies in atomic hydrogen (HI) gas. We find that only ~28% of the eighteen Seyfert galaxies in that sample are visibly disturbed in optical starlight. By contrast, ~94% of the same Seyfert galaxies are disturbed spatially and usually also kinematically in HI gas on galactic scales of >~20 kpc. In at least ~67% and up to perhaps ~94% of cases, the observed disturbances can be traced to tidal interactions with neighboring galaxies detected also in HI. The dramatic contrast between the observed prevalence of HI disturbances in the Seyfert and control samples implicates tidal interactions in initiating events that lead to luminous Seyfert activity in a large fraction of local disk galaxies.

No mechanisms have hitherto been conclusively demonstrated to be responsible for initiating optically-luminous nuclear (Seyfert) activity in local disk galaxies. Only a small minority of such galaxies are visibly disturbed in optical starlight, with the observed disturbances being at best marginally stronger than those found in matched samples of inactive galaxies. Here, we report the first systematic study of an optically-selected sample of twenty-three active galaxies in atomic hydrogen (HI) gas, which is the most sensitive and enduring tracer known of tidal interactions. Eighteen of these galaxies are (generally) classified as Seyferts, with over half (and perhaps all) having [OIII] luminosities within two orders of magnitude of Quasi-Stellar Objects. Only ~28% of these Seyfert galaxies are visibly disturbed in optical DSS2 images. By contrast, ~94% of the same galaxies are disturbed in HI, in nearly all cases not just spatially but also kinematically on galactic (>~20 kpc) scales. In at least ~67% and perhaps up to ~94% of cases, the observed HI disturbances can be traced to tidal interactions with neighboring galaxies detected also in HI. The majority of these neighboring galaxies have projected separations of <~ 100 kpc and differ in radial velocities by <~100 km/s from their respective Seyfert galaxies, and many have optical luminosities ranging from the Small to Large Magellanic Clouds. In a companion paper, we show that only ~15% of a matched control sample of inactive galaxies display comparable HI disturbances. Our results suggest that: i) most Seyfert galaxies (with high nuclear luminosities) have experienced tidal interactions in the recent past; ii) in most cases, these tidal interactions are responsible for initiating events that lead to their nuclear activity.

An X17.2 solar flare occurred on 2003 October 28, accompanied by multi-wavelength emissions and a high flux of relativistic particles observed at 1AU. We present the analytic results of the TRACE, SOHO, RHESSI, ACE, GOES, hard X-ray (INTEGRAL satellite), radio (Onderejov radio telescope), and neutron monitor data. It is found that the inferred magnetic reconnection electric field correlates well with the hard X-ray, gamma-ray, and neutron emission at the Sun. Thus the flare's magnetic reconnection probably makes a crucial contribution to the prompt relativistic particles, which could be detected at 1 AU. Since the neutrons were emitted a few minutes before the injection of protons and electrons, we propose a magnetic-field evolution configuration to explain this delay. We do not exclude the effect of CME-driven shock, which probably plays an important role in the delayed gradual phase of solar energetic particles.

In large gradual solar energetic particle (SEP) events, especially the ground-level enhancement (GLE) events, where and how energetic particles are accelerated is still a problem. By using imaging data from TRACE, Yohkoh/HXT, SOHO/MDI and SOHO/EIT, along with the data from the GOES, Apatity NM, and SOHO/LASCO CME catalog, the evolution of the X5.7 two-ribbon flare and the associated SEP event on 14 July 2000 are studied. It is found that the magnetic reconnection in this event consists of two parts, and the induced electric field Erec is temporally correlated with the evolution of hard X-ray and gamma-ray emission. In particular, the first hard X-ray and gamma-ray emission peak occurred at 10:22 UT, corresponding to the magnetic reconnection in the western part of the flare ribbons and the maximum Erec of 9.5 V/cm; the second emission peak at 10:27 UT, corresponding to the eastern part and the maximum Erec of 13.0 V/cm. We also analyze the SEP injection profiles as functions of time and CME-height, and find two-component injection which may result from different acceleration mechanisms. A reasonable conclusion is that reconnection electric field makes a crucial contribution to the acceleration of relativistic particles and to the impulsive component of the large gradual SEP event, while CME-driven shocks play a dominant role in the gradual component.

The purpose of this work is to search the Yangbajing neutron monitor data obtained between 1998 October and 2000 June for solar neutrons associated with solar flares. Using the onset times of 166 BATSE-detected flares with the GOES peak flux (1 -- 8 \AA) higher than $1.0 \times 10^{-5}$ $\mathrm{Wm^{-2}}$, we prepare for each flare a light curve of the Yangbajing neutron monitor, spanning $\pm$ 1.5 hours from the BATSE onset time. Based on the light curves, a systematic search for solar neutrons in energies above 100 MeV from the 166 flares was performed. No statistically significant signals due to solar neutrons were found in the present work. Therefore, we put upper limits on the $>$ 100 MeV solar-neutron flux for 18 events consisting of 2 X and 16 M class flares. The calculation assumed a power-law shaped neutron energy spectrum and three types of neutron emission profiles at the Sun. Compared with the other positive neutron detections associated with X-class flares, typical 95% confidence level upper limits for the two X-class flares are found to be comparable to the lowest and second lowest neutron fluxes at the top of the atmosphere.In addition, the upper limits for M-class flares scatter in the range of $10^{-2}$ to 1 neutrons $\mathrm{cm^{-2}s^{-1}}$. This provides the first upper limits on the solar-neutron flux from M-class solar flares, using space observatories as well as ground-based neutron monitors.

The consortium for Galactic studies with the Arecibo L-band Feed Array (ALFA) is conducting a neutral hydrogen (HI) survey of the whole Arecibo sky (declination range from -1 to 38 degrees), with high angular (3.5') and velocity resolution (0.2 km/s). The precursor observations with ALFA of a region in the Galactic anti-center reveal numerous isolated, small (a few pc in size), and cold (T_k<400 K) HI clouds at low negative velocities, distinctly separated from the HI disk emission (`low-velocity clouds', LVCs). These clouds are most likely located in the transition region between the Galactic disk and halo (at scale heights of 60--900 pc), yet they have properties of typical cold neutral clouds. LVCs are colder and, most likely, smaller and less massive than Lockman's clouds in the disk/halo interface region of the inner Galaxy. Our observations demonstrate that the cloudy structure of the interface region is most likely a general phenomenon, not restricted to the inner Galaxy. LVCs have sizes and radial velocities in agreement with the expectations for clouds formed in low-temperature fountain flows, although we measure a factor of ten higher HI column densities. Alternatively, LVCs could represent the final stages of the infalling intergalactic material in the on-going construction of the Galaxy. In the same dataset at higher negative velocities, we have discovered a `companion' HI cloud located 50' southwest of CHVC186+19-114. CHVC186+19-114 is a typical compact high velocity cloud (HVC) with a well-defined core/envelope structure. The companion cloud has a diameter of only 7'x9', and is one of the smallest HVCs known, most likely stripped from the main cloud through the interactions with the halo medium.

The 352 MHz Westerbork In the Southern Hemisphere (WISH) survey is the southern extension of the WENSS, covering 1.60 sr between -9 < DEC < -26 to a limiting flux density of ~18 mJy (5sigma). Due to the very low elevation of the observations, the survey has a much lower resolution in declination than in right ascension (54" x 54"cosec(DEC)). A correlation with the 1.4 GHz NVSS shows that the positional accuracy is less constrained in declination than in right ascension, but there is no significant systematic error. We present a source list containing 73570 sources. We correlate this WISH catalogue with the NVSS to construct a sample of faint Ultra Steep Spectrum (USS) sources, which is accessible for follow-up studies with large optical telescopes in the southern hemisphere. This sample is aimed at increasing the number of known high redshift radio galaxies to allow detailed follow-up studies of these massive galaxies and their environments in the early Universe.

An FRII radio source with an angular extension of 18.4 arcminutes has been discovered in the Westerbork Northern Sky Survey that is being carried out at 325 MHz. This source is identified with Mrk 1498 at $z=0.056$ indicating that its projected extension is 1.6 \h50\ Mpc. The optical spectrum of Mrk 1498 has the typical emission lines of a narrow line radio galaxy; H$\alpha$ is the only permitted line that has a faint broad component. The discovery shows that the good sensitivity of WENSS at low frequencies to low surface brightness features at scales of a few minutes, is uniquely suited to provide well defined and relatively large samples of giant radio galaxies. The future definition of such samples will be important for our understanding of giant radio galaxies.