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.

M64 is a nearby spiral galaxy with a Type III anti-truncation component. To trace the origin of the Type III component, we present Hubble Space Telescope/Advanced Camera for Surveys $F606W/F814W$ photometry of resolved stars in a field located in the outer disk ($2.5' \lesssim r \lesssim 6.5'$) of M64. At $r\approx 5.5'$ (7 kpc) to the east, we discover a new metal-poor globular cluster ($R_{\rm eff}=5.73\pm0.02$ pc and $M_V=-9.54\pm0.09$ mag), M64-GC1. This is the first globular cluster found in M64. The color-magnitude diagram (CMD) of the resolved stars in M64-GC1 is well matched by 12 Gyr isochrones with [Fe/H] $=-1.5\pm0.2$, showing that this cluster belongs to a halo. The CMD of the resolved stars in the entire ACS field shows two distinguishable red giant branches (RGBs): a curved metal-rich RGB and a vertical metal-poor RGB. The metal-rich RGB represents an old metal-rich ([Fe/H] $\approx -0.4$) disk population. In contrast, the CMD of the metal-poor RGB stars is very similar to the CMD of M64-GC1, showing that the metal-poor RGB represents a halo population. The radial number density profile of the metal-rich RGB stars is described by an exponential disk law, while the profile of the metal-poor RGB stars is described by a de Vaucouleurs's law. From these, we conclude that the origin of the Type III component in M64 is a halo which has a much lower metallicity than a disk or bulge population.

We present results of continuum and spectral line observations with ALMA and 22 GHz water (H$_2$O) maser observations using KaVA and VERA toward a high-mass star-forming region, G25.82-0.17. Multiple 1.3 mm continuum sources are revealed, indicating the presence of young stellar objects (YSOs) at different evolutionary stages, namely an ultra-compact HII region, G25.82-E, a high-mass young stellar object (HM-YSO), G25.82-W1, and starless cores, G25.82-W2 and G25.82-W3. Two SiO outflows, at N-S and SE-NW orientations, are identified. The CH$_3$OH 8$_{-1}$-7$_{0}$ E line, known to be a class I CH$_3$OH maser at 229 GHz is also detected showing a mixture of thermal and maser emission. Moreover, the H$_2$O masers are distributed in a region ~0.25" shifted from G25.82-W1. The CH$_3$OH 22$_{4}$-21$_{5}$ E line shows a compact ring-like structure at the position of G25.82-W1 with a velocity gradient, indicating a rotating disk or envelope. Assuming Keplerian rotation, the dynamical mass of G25.82-W1 is estimated to be $>$25 M$_{\odot}$ and the total mass of 20 M$_\odot$-84 M$_\odot$ is derived from the 1.3 mm continuum emission. The driving source of the N-S SiO outflow is G25.82-W1 while that of the SE-NW SiO outflow is uncertain. Detection of multiple high-mass starless$/$protostellar cores and candidates without low-mass cores implies that HM-YSOs could form in individual high-mass cores as predicted by the turbulent core accretion model. If this is the case, the high-mass star formation process in G25.82 would be consistent with a scaled-up version of low-mass star formation.

We studied a flare onset process in terms of stability of a three-dimensional (3D) magnetic field in active region 12371 producing an eruptive M6.5 flare in 2015 June 22. In order to reveal the 3D magnetic structure, we first extrapolated the 3D coronal magnetic fields based on time series of the photospheric vector magnetic fields under a nonlinear force-free field (NLFFF) approximation. The NLFFFs nicely reproduced the observed sigmoidal structure which is widely considered to be preeruptive magnetic configuration. In particular, we found that the sigmoid is composed of two branches of sheared arcade loops. On the basis of the NLFFFs, we investigated the sheared arcade loops to explore the onset process of the eruptive flare using three representative magnetohydrodynamic instabilities: the kink, torus, and double arc instabilities (DAI). The DAI, recently proposed by Ishiguro & Kusano, is a double arc loop that can be more easily destabilized than a torus loop. Consequently, the NLFFFs are found to be quite stable against the kink and torus instabilities. However, the sheared arcade loops formed prior to the flare possibly become unstable against the DAI. As a possible scenario for the onset process of the M6.5 flare, we suggest a three-step process: (1) double arc loops are formed by the sheared arcade loops through the tether-cutting reconnection during an early phase of the flare, (2) the DAI contributes to the expansion of destabilized double arc loops, and (3) finally, the torus instability makes the full eruption.

We report the discovery of ultra-diffuse galaxies (UDGs) in Abell 370 ($z=0.375$). We find 46 UDGs in Abell 370 from the images of the Hubble Frontier Fields (HFF). Most UDGs are low-luminosity red sequence galaxies, while a few of them are blue UDGs. We estimate the abundance of UDGs in Abell 370, $N(\rm UDG)=644\pm104$. Combining these results with those of Abell S1063 ($z=0.348$) and Abell 2744 ($z=0.308$) \citep{Lee17}, we derive a mean radial number density profile of UDGs in the three clusters. The number density profiles of UDGs and bright galaxies show a discrepancy in the central region of the clusters: the profile of UDGs shows a flattening as clustercentric distance decreases, while that of bright galaxies shows a continuous increase. This implies that UDGs are prone to disruption in the central region of the clusters. The relation between the abundance of UDGs and virial masses of their host systems is described by a power-law with an index nearly one: $N({\rm UDG})\propto M_{200}^{0.99\pm0.05}$ for $M_{200}>10^{13}~M_{\odot}$. We estimate approximately dynamical masses of UDGs using the fundamental manifold method, and find that most UDGs have dwarf-like masses $(M_{200}<10^{11}$ $M_{\odot})$. This implies that most UDGs have a dwarf-like origin and a small number of them could be failed $L^{*}$ galaxies. These results suggest that multiple origins may contribute to the formation and evolution of UDGs in massive galaxy clusters.

The BICEP/Keck (BK) experiment aims to detect the imprint of primordial gravitational waves in the Cosmic Microwave Background polarization, which would be direct evidence of the inflation theory. While the tensor-to-scalar ratio has been constrained to be r_0.05 < 0.06 at 95% c.l., further improvements on this upper limit are hindered by polarized Galactic foreground emissions and removal of gravitational lensing polarization. The 30/40 GHz receiver of the BICEP Array (BA) will deploy at the end of 2019 and will constrain the synchrotron foreground with unprecedented accuracy within the BK sky patch. We will show the design of the 30/40 GHz detectors and test results summarizing its performance. The low optical and atmospheric loading at these frequencies requires our TES detectors to have low saturation power in order to be photon-noise dominated. To realize the low thermal conductivity required from a 250 mK base temperature, we developed new bolometer leg designs. We will present the relevant measured detector parameters: G, Tc, Rn, Psat , and spectral bands, and noise spectra. We achieved a per bolometer NEP including all noise components of 2.07E-17 W/sqrt(Hz), including an anticipated photon noise level 1.54E-17 W/sqrt(Hz).

The BICEP/Keck experiment (BK) is a series of small-aperture refracting telescopes observing degree-scale Cosmic Microwave Background (CMB) polarization from the South Pole in search of a primordial $B$-mode signature. This $B$-mode signal arises from primordial gravitational waves interacting with the CMB, and has amplitude parametrized by the tensor-to-scalar ratio $r$. Since 2016, BICEP3 and the Keck Array have been observing with 4800 total antenna-coupled transition-edge sensor detectors, with frequency bands spanning 95, 150, 220, and 270 GHz. Here we present the optical performance of these receivers from 2016 to 2019, including far-field beams measured in situ with an improved chopped thermal source and instrument spectral response measured with a field-deployable Fourier Transform Spectrometer. As a pair differencing experiment, an important systematic that must be controlled is the differential beam response between the co-located, orthogonally polarized detectors. We generate per-detector far-field beam maps and the corresponding differential beam mismatch that is used to estimate the temperature-to-polarization leakage in our CMB maps and to give feedback on detector and optics fabrication. The differential beam parameters presented here were estimated using improved low-level beam map analysis techniques, including efficient removal of non-Gaussian noise as well as improved spatial masking. These techniques help minimize systematic uncertainty in the beam analysis, with the goal of constraining the bias on $r$ induced by temperature-to-polarization leakage to be subdominant to the statistical uncertainty. This is essential as we progress to higher detector counts in the next generation of CMB experiments.

Families of cosmic inflation models predict a primordial gravitational-wave background that imprints B-mode polarization pattern in the Cosmic Microwave Background (CMB). High sensitivity instruments with wide frequency coverage and well-controlled systematic errors are needed to constrain the faint B-mode amplitude. We have developed antenna-coupled Transition Edge Sensor (TES) arrays for high-sensitivity polarized CMB observations over a wide range of millimeter-wave bands. BICEP Array, the latest phase of the BICEP/Keck experiment series, is a multi-receiver experiment designed to search for inflationary B-mode polarization to a precision $\sigma$(r) between 0.002 and 0.004 after 3 full years of observations, depending on foreground complexity and the degree of lensing removal. We describe the electromagnetic design and measured performance of BICEP Array low-frequency 40-GHz detector, their packaging in focal plane modules, and optical characterization including efficiency and beam matching between polarization pairs. We summarize the design and simulated optical performance, including an approach to improve the optical efficiency due to mismatch losses. We report the measured beam maps for a new broad-band corrugation design to minimize beam differential ellipticity between polarization pairs caused by interactions with the module housing frame, which helps minimize polarized beam mismatch that converts CMB temperature to polarization ($T \rightarrow P$) anisotropy in CMB maps.

Branches of cosmic inflationary models, such as slow-roll inflation, predict a background of primordial gravitational waves that imprints a unique odd-parity B-mode pattern in the Cosmic Microwave Background (CMB) at amplitudes that are within experimental reach. The BICEP/Keck (BK) experiment targets this primordial signature, the amplitude of which is parameterized by the tensor-to-scalar ratio r, by observing the polarized microwave sky through the exceptionally clean and stable atmosphere at the South Pole. B-mode measurements require an instrument with exquisite sensitivity, tight control of systematics, and wide frequency coverage to disentangle the primordial signal from the Galactic foregrounds. BICEP Array represents the most recent stage of the BK program, and comprises four BICEP3-class receivers observing at 30/40, 95, 150 and 220/270 GHz. The 30/40 GHz receiver will be deployed at the South Pole during the 2019/2020 austral summer. After 3 full years of observations with 30,000+ detectors, BICEP Array will measure primordial gravitational waves to a precision $\sigma(r)$ between 0.002 and 0.004, depending on foreground complexity and the degree of lensing removal. In this paper we give an overview of the instrument, highlighting the design features in terms of cryogenics, magnetic shielding, detectors and readout architecture as well as reporting on the integration and tests that are ongoing with the first receiver at 30/40 GHz.

We describe an on-sky demonstration of a microwave-multiplexing readout system in one of the receivers of the Keck Array, a polarimetry experiment observing the cosmic microwave background at the South Pole. During the austral summer of 2018-2019, we replaced the time-division multiplexing readout system with microwave-multiplexing components including superconducting microwave resonators coupled to radio-frequency superconducting quantum interference devices at the sub-Kelvin focal plane, coaxial-cable plumbing and amplification between room temperature and the cold stages, and a SLAC Microresonator Radio Frequency system for the warm electronics. In the range 5-6 GHz, a single coaxial cable reads out 528 channels. The readout system is coupled to transition-edge sensors, which are in turn coupled to 150-GHz slot-dipole phased-array antennas. Observations began in April 2019, and we report here on an initial characterization of the system performance.

Ultra-light axions (ULAs) with mass less than 10^-20 eV have interesting behaviors that may contribute to either dark energy or dark matter at different epochs of the Universe. Its properties can be explored by cosmological observations, such as expansion history of the Universe, cosmic large-scale structure, cosmic microwave background, etc. In this work, we study the ULAs with a mass around 10^-33 eV, which means the ULA field still rolls slowly at present with the equation of state w=-1 as dark energy. In order to investigate the mass and other properties of this kind of ULA field, we adopt the measurements of Type Ia supernova (SN Ia), baryon acoustic oscillation (BAO), and Hubble parameter H(z). The Markov Chain Monte Carlo (MCMC) technique is employed to perform the constraints on the parameters. Finally, by exploring four cases of the model, we find that the mass of this ULA field is about 3x10^-33 eV if assuming the initial axion field phi_i=M_pl. We also investigate a general case by assuming phi_i< M_pl and find that the fitting results of phi_i/M_pl are consistent with or close to 1 for the datasets we use.

Precision measurements of cosmic microwave background (CMB) polarization require extreme control of instrumental systematics. In a companion paper we have presented cosmological constraints from observations with the BICEP2 and Keck Array experiments up to and including the 2015 observing season (BK15), resulting in the deepest CMB polarization maps to date and a statistical sensitivity to the tensor-to-scalar ratio of $\sigma(r) = 0.020$. In this work we characterize the beams and constrain potential systematic contamination from main beam shape mismatch at the three BK15 frequencies (95, 150, and 220 GHz). Far-field maps of 7,360 distinct beam patterns taken from 2010-2015 are used to measure differential beam parameters and predict the contribution of temperature-to-polarization leakage to the BK15 B-mode maps. In the multifrequency, multicomponent likelihood analysis that uses BK15, Planck, and WMAP maps to separate sky components, we find that adding this predicted leakage to simulations induces a bias of $\Delta r = 0.0027 \pm 0.0019$. Future results using higher-quality beam maps and improved techniques to detect such leakage in CMB data will substantially reduce this uncertainty, enabling the levels of systematics control needed for BICEP Array and other experiments that plan to definitively probe large-field inflation.

We present the strongest constraints to date on anisotropies of cosmic microwave background (CMB) polarization rotation derived from 150 GHz data taken by the BICEP2/Keck Array CMB experiments up to and including the 2014 observing season (BK14). The definition of the polarization angle in BK14 maps has gone through self-calibration in which the overall angle is adjusted to minimize the observed TB and EB power spectra. After this procedure, the QU maps lose sensitivity to a uniform polarization rotation but are still sensitive to anisotropies of polarization rotation. This analysis places constraints on the anisotropies of polarization rotation, which could be generated by CMB photons interacting with axionlike pseudoscalar fields or Faraday rotation induced by primordial magnetic fields. The sensitivity of BK14 maps ($\sim 3\mu$K-arcmin) makes it possible to reconstruct anisotropies of the polarization rotation angle and measure their angular power spectrum much more precisely than previous attempts. Our data are found to be consistent with no polarization rotation anisotropies, improving the upper bound on the amplitude of the rotation angle spectrum by roughly an order of magnitude compared to the previous best constraints. Our results lead to an order of magnitude better constraint on the coupling constant of the Chern-Simons electromagnetic term $g_{a\gamma}\leq 7.2\times 10^{-2}/H_I$ (95% confidence) than the constraint derived from the B-mode spectrum, where $H_I$ is the inflationary Hubble scale. This constraint leads to a limit on the decay constant of $10^{-6}\lesssim f_a/M_{\rm pl}$ at mass range of $10^{-33}< m_a< 10^{-28}$ eV for $r=0.01$, assuming $g_{a\gamma}\sim\alpha/(2\pi f_a)$ with $\alpha$ denoting the fine structure constant. The upper bound on the amplitude of the primordial magnetic fields is 30nG (95% confidence) from the polarization rotation anisotropies.

Recently, spiral wave patterns (SWPs) have been detected in 3 minute oscillations of sunspot umbrae, but the nature of this phenomenon has remained elusive. We present a theoretical model that interprets the observed SWPs as the superposition of two different azimuthal modes of slow magnetoacoustic waves driven below the surface in an untwisted and non-rotating magnetic cylinder. We apply this model to SWPs of the line-of-sight (LOS) velocity in a pore observed by the Fast Imaging Solar Spectrograph installed at the 1.6 m Goode Solar Telescope. One- and two-armed SWPs were identified in instantaneous amplitudes of LOS Doppler velocity maps of 3 minute oscillations. The associated oscillation periods are about 160 s, and the durations are about 5 minutes. In our theoretical model, the observed spiral structures are explained by the superposition of non-zero azimuthal modes driven 1600 km below the photosphere in the pore. The one-armed SWP is produced by the slow-body sausage (m = 0) and kink (m = 1) modes, and the two-armed SWP is formed by the slow-body sausage (m = 0) and fluting (m = 2) modes of the magnetic flux tube forming the pore.

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 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 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.

Intensive Monitoring Survey of Nearby Galaxies (IMSNG) is a high cadence observation program monitoring nearby galaxies with high probabilities of hosting supernovae (SNe). IMSNG aims to constrain the SN explosion mechanism by inferring sizes of SN progenitor systems through the detection of the shock-heated emission that lasts less than a few days after the SN explosion. To catch the signal, IMSNG utilizes a network of 0.5-m to 1-m class telescopes around the world and monitors the images of 60 nearby galaxies at distances D < 50 Mpc to a cadence as short as a few hours. The target galaxies are bright in near-ultraviolet (NUV) with M_NUV < -18.4 AB mag and have high probabilities of hosting SNe (0.06 SN/yr per galaxy). With this strategy, we expect to detect the early light curves of 3.4 SNe per year to a depth of R ~ 19.5 mag, enabling us to detect the shock-heated emission from a progenitor star with a radius as small as 0.1 R_sun. The accumulated data will be also useful for studying faint features around the target galaxies and other science projects. So far, 18 SNe have occurred in our target fields (16 in IMSNG galaxies) over 5 years, confirming our SN rate estimate of 0.06 SN/yr per galaxy.

_NGC 4589, a bright E2 merger-remnant galaxy, hosts the peculiar fast and faint calcium-rich Type Ib supernova (SN) SN 2005cz. The progenitor of Ca-rich SNe Ib has been controversial: it could be a) a young massive star with 6-12 M$\odot$ in a binary system, or b) an old low-mass star in a binary system that was kicked out from the galaxy center. Moreover, previous distance estimates for this galaxy have shown a large spread, ranging from 20 Mpc to 60 Mpc. Thus, using archival $Hubble$ $Space$ $Telescope$/ACS $F435W$, $F555W$, and $F814W$ images, we search for star clusters in NGC 4589 in order to help resolve these issues. We find a small population of young star clusters with $25<V\leq27$ ($-7.1<M_V\leq-5.1$) mag and age $< 1$ Gyr in the central region at $R<0.5'$ ($<3.8$ kpc), thus supporting the massive-star progenitor scenario for SN 2005cz. In addition to young star clusters, we also find a large population of old globular clusters. In contrast to previous results in the literature, we find that the color distribution of the globular clusters is clearly bimodal. The turnover (Vega) magnitude in the $V$-band luminosity functions of the blue (metal-poor) globular clusters is determined to be $V_0{(\rm max)}=24.40\pm0.10$ mag. We derive the total number of globular clusters, $N_{\rm GC} =640\pm50$, and the specific frequency, $S_N =1.7\pm0.2$. Adopting a calibration for the metal-poor globular clusters, $M_V({\rm max})=-7.66\pm0.14$ mag, we derive a distance to this galaxy: $(m-M)_0=32.06\pm0.10({\rm ran})\pm0.15({\rm sys})$ ($d=25.8\pm2.2$ Mpc).

The BICEP and Keck Array experiments are a suite of small-aperture refracting telescopes observing the microwave sky from the South Pole. They target the degree-scale B-mode polarization signal imprinted in the Cosmic Microwave Background (CMB) by primordial gravitational waves. Such a measurement would shed light on the physics of the very early universe. While BICEP2 observed for the first time a B-mode signal at 150 GHz, higher frequencies from the Planck satellite showed that it could be entirely due to the polarized emission from Galactic dust, though uncertainty remained high. Keck Array has been observing the same region of the sky for several years, with an increased detector count, producing the deepest polarized CMB maps to date. New detectors at 95 GHz were installed in 2014, and at 220 GHz in 2015. These observations enable a better constraint of galactic foreground emissions, as presented here. In 2015, BICEP2 was replaced by BICEP3, a 10 times higher throughput telescope observing at 95 GHz, while Keck Array is now focusing on higher frequencies. In the near future, BICEP Array will replace Keck Array, and will allow unprecedented sensitivity to the gravitational wave signal. High resolution observations from the South Pole Telescope (SPT) will also be used to remove the lensing contribution to B-modes.

We present results from an analysis of all data taken by the BICEP2/Keck CMB polarization experiments up to and including the 2015 observing season. This includes the first Keck Array observations at 220 GHz and additional observations at 95 & 150 GHz. The $Q/U$ maps reach depths of 5.2, 2.9 and 26 $\mu$K$_{cmb}$ arcmin at 95, 150 and 220 GHz respectively over an effective area of $\approx 400$ square degrees. The 220 GHz maps achieve a signal-to-noise on polarized dust emission approximately equal to that of Planck at 353 GHz. We take auto- and cross-spectra between these maps and publicly available WMAP and Planck maps at frequencies from 23 to 353 GHz. We evaluate the joint likelihood of the spectra versus a multicomponent model of lensed-$\Lambda$CDM+$r$+dust+synchrotron+noise. The foreground model has seven parameters, and we impose priors on some of these using external information from Planck and WMAP derived from larger regions of sky. The model is shown to be an adequate description of the data at the current noise levels. The likelihood analysis yields the constraint $r_{0.05}<0.07$ at 95% confidence, which tightens to $r_{0.05}<0.06$ in conjunction with Planck temperature measurements and other data. The lensing signal is detected at $8.8 \sigma$ significance. Running maximum likelihood search on simulations we obtain unbiased results and find that $\sigma(r)=0.020$. These are the strongest constraints to date on primordial gravitational waves.

Bicep Array is a cosmic microwave background (CMB) polarization experiment that will begin observing at the South Pole in early 2019. This experiment replaces the five Bicep2 style receivers that compose the Keck Array with four larger Bicep3 style receivers observing at six frequencies from 30 to 270GHz. The 95GHz and 150GHz receivers will continue to push the already deep Bicep/Keck CMB maps while the 30/40GHz and 220/270GHz receivers will constrain the synchrotron and galactic dust foregrounds respectively. Here we report on the design and performance of the Bicep Array instruments focusing on the mount and cryostat systems.

BICEP Array is a degree-scale Cosmic Microwave Background (CMB) experiment that will search for primordial B-mode polarization while constraining Galactic foregrounds. BICEP Array will be comprised of four receivers to cover a broad frequency range with channels at 30/40, 95, 150 and 220/270 GHz. The first low-frequency receiver will map synchrotron emission at 30 and 40 GHz and will deploy to the South Pole at the end of 2019. In this paper, we give an overview of the BICEP Array science and instrument, with a focus on the detector module. We designed corrugations in the metal frame of the module to suppress unwanted interactions with the antenna-coupled detectors that would otherwise deform the beams of edge pixels. This design reduces the residual beam systematics and temperature-to-polarization leakage due to beam steering and shape mismatch between polarized beam pairs. We report on the simulated performance of single- and wide-band corrugations designed to minimize these effects. Our optimized design alleviates beam differential ellipticity caused by the metal frame to about 7% over 57% bandwidth (25 to 45 GHz), which is close to the level due the bare antenna itself without a metal frame. Initial laboratory measurements are also presented.

Targeting faint polarization patterns arising from Primordial Gravitational Waves in the Cosmic Microwave Background requires excellent observational sensitivity. Optical elements in small aperture experiments such as Bicep3 and Keck Array are designed to optimize throughput and minimize losses from transmission, reflection and scattering at millimeter wavelengths. As aperture size increases, cryostat vacuum windows must withstand larger forces from atmospheric pressure and the solution has often led to a thicker window at the expense of larger transmission loss. We have identified a new candidate material for the fabrication of vacuum windows: with a tensile strength two orders of magnitude larger than previously used materials, woven high-modulus polyethylene could allow for dramatically thinner windows, and therefore significantly reduced losses and higher sensitivity. In these proceedings we investigate the suitability of high-modulus polyethylene windows for ground-based CMB experiments, such as current and future receivers in the Bicep/Keck Array program. This includes characterizing their optical transmission as well as their mechanical behavior under atmospheric pressure. We find that such ultra-thin materials are promising candidates to improve the performance of large-aperture instruments at millimeter wavelengths, and outline a plan for further tests ahead of a possible upcoming field deployment of such a science-grade window.

BICEP Array is the newest multi-frequency instrument in the BICEP/Keck Array program. It is comprised of four 550 mm aperture refractive telescopes observing the polarization of the cosmic microwave background (CMB) at 30/40, 95, 150 and 220/270 GHz with over 30,000 detectors. We present an overview of the receiver, detailing the optics, thermal, mechanical, and magnetic shielding design. BICEP Array follows BICEP3's modular focal plane concept, and upgrades to 6" wafer to reduce fabrication with higher detector count per module. The first receiver at 30/40 GHz is expected to start observing at the South Pole during the 2019-20 season. By the end of the planned BICEP Array program, we project $\sigma(r) \sim 0.003$, assuming current modeling of polarized Galactic foreground and depending on the level of delensing that can be achieved with higher resolution maps from the South Pole Telescope.

BICEP3 is a 520mm aperture on-axis refracting telescope observing the polarization of the cosmic microwave background (CMB) at 95GHz in search of the B-mode signal originating from inflationary gravitational waves. BICEP3's focal plane is populated with modularized tiles of antenna-coupled transition edge sensor (TES) bolometers. BICEP3 was deployed to the South Pole during 2014-15 austral summer and has been operational since. During the 2016-17 austral summer, we implemented changes to optical elements that lead to better noise performance. We discuss this upgrade and show the performance of BICEP3 at its full mapping speed from the 2017 and 2018 observing seasons. BICEP3 achieves an order-of-magnitude improvement in mapping speed compared to a Keck 95GHz receiver. We demonstrate $6.6\mu K\sqrt{s}$ noise performance of the BICEP3 receiver.

For icosahedral inflation, we compute the tensor modes' two-point function in the presence of higher derivative corrections, and show that in general this features anisotropies that are aligned with the underlying icosahedral structure. The effect is small within the regime of validity of the effective theory, but it is the leading contribution to a mixed correlator between the two different helicities. We also estimate the magnitude of a similar effect for a mixed scalar-tensor correlator, whose detailed computation we leave for future work. Along the way, we clarify a number of aspects of the spin decomposition of generic icosahedral-invariant physical quantities.

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.

NGC 4993 hosts a binary neutron star merger emitting gravitational waves and electromagnetic waves, GW170817/GRB 170817A. The distance to this galaxy is not well established. We select the globular cluster candidates from the Hubble Space Telescope/ACS F606W images of NGC 4993 in the archive, using the structural parameters of the detected sources. The radial number density distribution of these candidates shows a significant central concentration around the galaxy center at the galactocentric distance $r<50''$, showing that they are mostly the members of NGC 4993. Also the luminosity function of these candidates is fit well by a Gaussian function. Therefore the selected candidates at $r<50''$ are mostly considered to be globular clusters in NGC 4993. We derive an extinction-corrected turnover Vega magnitude in the luminosity function of the globular clusters at $20''<r<50''$, F606W (max)$_0= 25.36\pm0.08$ ($V_0 =25.52\pm0.11$)} mag. Adopting the calibration of the turnover magnitudes of the globular clusters, $M_V({\rm max})=-7.58\pm0.11$, we derive a distance to NGC 4993, $d=41.65\pm3.00$ Mpc ($(m-M)_0=33.10\pm0.16$). The systematic error of this method can be as large as $\pm0.3$ mag. This value is consistent with the previous distance estimates based on the fundamental plane relation and the gravitational wave method in the literature. The distance in this study can be used to constrain the values of the parameters including the inclination angle of the binary system in the models of gravitational wave analysis.

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.

We study reheating in $\alpha$-attractor models of inflation in which the inflaton couples to other scalars or fermions. We show that the parameter space contains viable regions in which the inflaton couplings to radiation can be determined from the properties of CMB temperature fluctuations, in particular the spectral index. This may be the only way to measure these fundamental microphysical parameters, which shaped the universe by setting the initial temperature of the hot big bang and contain important information about the embedding of a given model of inflation into a more fundamental theory of physics. The method can be applied to other models of single field inflation.

We present analysis of near-infrared, high-resolution spectroscopy towards the Flat-spectrum YSO (Young Stellar Object) ESO H$\alpha$ 279a (1.5 solar mass) in the Serpens star forming region, at the distance of 429 pc. Using the Immersion GRating INfrared Spectrometer (IGRINS, R=45,000), we detect emission lines originating from the accretion channel flow, jet, and inner disk. Specifically, we identify hydrogen Brackett series recombination, [Fe II], [Fe III], [Fe IV], Ca I, Na I, H2, H2O and CO overtone emission lines. By modeling five bands of CO overtone emission lines, and the symmetric double-peaked line profile for Na I emission lines, we find that ESO H$\alpha$ 279a has an actively accreting Keplerian disk. From our Keplerian disk model, we find that Na I emission lines originate between 0.04 AU and 1.00 AU, while CO overtone emission lines are from the outer part of disk, in the range between 0.22 AU and 3.00 AU. It reveals that the neutral atomic Na gas is a good tracer of the innermost region of the actively accreting disk. We derive a mass accretion rate of 2-10x10^{-7} M_solar/yr from the measured Br_gamma emission luminosity of 1.78(+-0.31)x10^{31} erg/s.

We present the detection of a large population of ultra diffuse galaxies (UDGs) in two massive galaxy clusters, Abell S1063 at $z=0.348$ and Abell 2744 at $z=0.308$, based on F814W and F105W images in the Hubble Frontier Fields Program. We find 47 and 40 UDGs in Abell S1063 and Abell 2744, respectively. Color-magnitude diagrams of the UDGs show that they are mostly located at the faint end of the red sequence. From the comparison with simple stellar population models, we estimate their stellar mass to range from $10^8$ to $10^9 M_\odot$. Radial number density profiles of the UDGs show a turnover or a flattening in the central region at $r<100$ kpc. We estimate the total masses of the UDGs using the galaxy scaling relations. A majority of the UDGs have total masses, $M_{200} = 10^{10}$ to $10^{11}~M_\odot$, and only a few of them have total masses, $M_{200} = 10^{11}$ to $10^{12}~M_\odot$. The total number of UDGs within the virial radius is estimated to be N(UDG)$=770\pm114$ for Abell S1063, and N(UDG)$=814\pm122$ for Abell 2744. Combining these results with data in the literature, we fit the relation between the total numbers of UDGs and the masses of their host systems for $M_{200}>10^{13} M_\odot$ with a power law, N(UDG) $= M_{200}^{1.05\pm0.09}$. These results suggest that a majority of the UDGs have a dwarf galaxy origin, while only a small number of the UDGs are massive $L_*$ galaxies that failed to form a normal population of stars.

We present the results of very long baseline interferometry (VLBI) observations of gamma-ray bright blazar S5 0716+714 using the Korean VLBI Network (KVN) at the 22, 43, 86, and 129 GHz bands, as part of the Interferometric Monitoring of Gamma-ray Bright AGNs (iMOGABA) KVN key science program. Observations were conducted in 29 sessions from January 16, 2013 to March 1, 2016, with the source being detected and imaged at all available frequencies. In all epochs, the source was compact on the milliarcsecond (mas) scale, yielding a compact VLBI core dominating the synchrotron emission on these scales. Based on the multi-wavelength data between 15 GHz (Owens Valley Radio Observatory) and 230 GHz (Submillimeter Array), we found that the source shows multiple prominent enhancements of the flux density at the centimeter (cm) and millimeter (mm) wavelengths, with mm enhancements leading cm enhancements by -16$\pm$8 days. The turnover frequency was found to vary between 21 to 69GHz during our observations. By assuming a synchrotron self-absorption model for the relativistic jet emission in S5 0716+714, we found the magnetic field strength in the mas emission region to be $\le$5 mG during the observing period, yielding a weighted mean of 1.0$\pm$0.6 mG for higher turnover frequencies (e.g., >45 GHz).

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.

Class I methanol masers are found to be good tracers of the interaction between outflows from massive young stellar objects with their surrounding media. Although polarization observations of Class II methanol masers have been able to provide information about magnetic fields close to the central (proto)stars, polarization observations of Class I methanol masers are rare, especially at 44 and 95 GHz. We present the results of linear polarization observations of 39 Class I methanol maser sources at 44 and 95 GHz. These two lines are observed simultaneously with one of the 21 m Korean VLBI Network telescopes in single-dish mode. Approximately 60% of the observed sources have fractional polarizations of a few percent in at least one transition. This is the first reported detection of linear polarization of the 44 GHz methanol maser. The two maser transitions show similar polarization properties, indicating that they trace similar magnetic environments, although the fraction of the linear polarization is slightly higher at 95 GHz. We discuss the association between the directions of polarization angles and outflows. We also discuss some targets having different polarization properties at both lines, including DR21(OH) and G82.58+0.20, which show the 90{\deg} polarization angle flip at 44 GHz.

Neutral atomic hydrogen (HI) gas in interstellar space is largely organized into filaments, loops, and shells, the most prominent of which are "supershells". These gigantic structures requiring $\gtrsim 3 \times 10^{52}$ erg to form are generally thought to be produced by either the explosion of multiple supernovae (SNe) in OB associations or alternatively by the impact of high-velocity clouds (HVCs) falling to the Galactic disk. Here we report the detection of a kiloparsec (kpc)-size supershell in the outskirts of the Milky Way with the compact HVC 040+01$-$282 (hereafter CHVC040) at its geometrical center using the "Inner-Galaxy Arecibo L-band Feed Array" HI 21-cm survey data. The morphological and physical properties of both objects suggest that CHVC040, which is either a fragment of a nearby disrupted galaxy or a cloud originated from an intergalactic accreting flow, collided with the disk $\sim 5$ Myrs ago to form the supershell. Our result shows that some compact HVCs can survive their trip through the Galactic halo and inject energy and momentum into the Milky Way disk.

BICEP3, the latest telescope in the BICEP/Keck program, started science observations in March 2016. It is a 550mm aperture refractive telescope observing the polarization of the cosmic microwave background at 95 GHz. We show the focal plane design and detector performance, including spectral response, optical efficiency and preliminary sensitivity of the upgraded BICEP3. We demonstrate 9.72$\mu$K$\sqrt{\textrm{s}}$ noise performance of the BICEP3 receiver.

BICEP3 is a small-aperture refracting cosmic microwave background (CMB) telescope designed to make sensitive polarization maps in pursuit of a potential B-mode signal from inflationary gravitational waves. It is the latest in the BICEP/Keck Array series of CMB experiments at the South Pole, which has provided the most stringent constraints on inflation to date. For the 2016 observing season, BICEP3 was outfitted with a full suite of 2400 optically coupled detectors operating at 95 GHz. In these proceedings we report on the far field beam performance using calibration data taken during the 2015-2016 summer deployment season in situ with a thermal chopped source. We generate high-fidelity per-detector beam maps, show the array-averaged beam profile, and characterize the differential beam response between co-located, orthogonally polarized detectors which contributes to the leading instrumental systematic in pair differencing experiments. We find that the levels of differential pointing, beamwidth, and ellipticity are similar to or lower than those measured for BICEP2 and Keck Array. The magnitude and distribution of BICEP3's differential beam mismatch - and the level to which temperature-to-polarization leakage may be marginalized over or subtracted in analysis - will inform the design of next-generation CMB experiments with many thousands of detectors.

BICEP3 is a 520 mm aperture, compact two-lens refractor designed to observe the polarization of the cosmic microwave background (CMB) at 95 GHz. Its focal plane consists of modularized tiles of antenna-coupled transition edge sensors (TESs), similar to those used in BICEP2 and the Keck Array. The increased per-receiver optical throughput compared to BICEP2/Keck Array, due to both its faster f/1.7 optics and the larger aperture, more than doubles the combined mapping speed of the BICEP/Keck program. The BICEP3 receiver was recently upgraded to a full complement of 20 tiles of detectors (2560 TESs) and is now beginning its second year of observation (and first science season) at the South Pole. We report on its current performance and observing plans. Given its high per-receiver throughput while maintaining the advantages of a compact design, BICEP3-class receivers are ideally suited as building blocks for a 3rd-generation CMB experiment, consisting of multiple receivers spanning 35 GHz to 270 GHz with total detector count in the tens of thousands. We present plans for such an array, the new "BICEP Array" that will replace the Keck Array at the South Pole, including design optimization, frequency coverage, and deployment/observing strategies.

A linear polarization field on the sphere can be uniquely decomposed into an E-mode and a B-mode component. These two components are analytically defined in terms of spin-2 spherical harmonics. Maps that contain filtered modes on a partial sky can also be decomposed into E-mode and B-mode components. However, the lack of full sky information prevents orthogonally separating these components using spherical harmonics. In this paper, we present a technique for decomposing an incomplete map into E and B-mode components using E and B eigenmodes of the pixel covariance in the observed map. This method is found to orthogonally define E and B in the presence of both partial sky coverage and spatial filtering. This method has been applied to the BICEP2 and the Keck Array maps and results in reducing E to B leakage from LCDM E-modes to a level corresponding to a tensor-to-scalar ratio of $r<1\times10^{-4}$.

We calculate the production rate of singlet fermions from the decay of neutral or charged scalar fields in a hot plasma. We find that there are considerable thermal corrections when the temperature of the plasma exceeds the mass of the decaying scalar. We give analytic expressions for the temperature-corrected production rates in the regime where the decay products are relativistic. We also study the regime of non-relativistic decay products numerically. Our results can be used to determine the abundance and momentum distribution of Dark Matter particles produced in scalar decays. The inclusion of thermal corrections helps to improve predictions for the free streaming of the Dark Matter particles, which is crucial to test the compatibility of a given model with cosmic structure formation. With some modifications, our results may be generalised to the production of other Dark Matter candidates in scalar decays.

We present measurements of polarization lensing using the 150 GHz maps which include all data taken by the BICEP2 & Keck Array CMB polarization experiments up to and including the 2014 observing season (BK14). Despite their modest angular resolution ($\sim 0.5^\circ$), the excellent sensitivity ($\sim 3\mu$K-arcmin) of these maps makes it possible to directly reconstruct the lensing potential using only information at larger angular scales ($\ell\leq 700$). From the auto-spectrum of the reconstructed potential we measure an amplitude of the spectrum to be $A^{\phi\phi}_{\rm L}=1.15\pm 0.36$ (Planck $\Lambda$CDM prediction corresponds to $A^{\phi\phi}_{\rm L}=1$), and reject the no-lensing hypothesis at 5.8$\sigma$, which is the highest significance achieved to date using an EB lensing estimator. Taking the cross-spectrum of the reconstructed potential with the Planck 2015 lensing map yields $A^{\phi\phi}_{\rm L}=1.13\pm 0.20$. These direct measurements of $A^{\phi\phi}_{\rm L}$ are consistent with the $\Lambda$CDM cosmology, and with that derived from the previously reported BK14 B-mode auto-spectrum ($A^{\rm BB}_{\rm L}=1.20\pm 0.17$). We perform a series of null tests and consistency checks to show that these results are robust against systematics and are insensitive to analysis choices. These results unambiguously demonstrate that the B-modes previously reported by BICEP / Keck at intermediate angular scales ($150\lesssim\ell\lesssim 350$) are dominated by gravitational lensing. The good agreement between the lensing amplitudes obtained from the lensing reconstruction and B-mode spectrum starts to place constraints on any alternative cosmological sources of B-modes at these angular scales.

BICEP3 is a $550~mm$ aperture telescope with cold, on-axis, refractive optics designed to observe at the $95~GHz$ band from the South Pole. It is the newest member of the BICEP/Keck family of inflationary probes specifically designed to measure the polarization of the cosmic microwave background (CMB) at degree-angular scales. BICEP3 is designed to house 1280 dual-polarization pixels, which, when fully-populated, totals to $\sim$9$\times$ the number of pixels in a single Keck $95~GHz$ receiver, thus further advancing the BICEP/Keck program's $95~GHz$ mapping speed. BICEP3 was deployed during the austral summer of 2014-2015 with 9 detector tiles, to be increased to its full capacity of 20 in the second season. After instrument characterization measurements were taken, CMB observation commenced in April 2015. Together with multi-frequency observation data from Planck, BICEP2, and the Keck Array, BICEP3 is projected to set upper limits on the tensor-to-scalar ratio to $r$ $\lesssim 0.03$ at $95\%$ C.L..

We generalize the model of solid inflation to an anisotropic cosmic solid. Barring fine tunings, the observed isotropy of the cosmological background and of the scalar two-point function isolate the icosahedral group as the only possible symmetry group of such a solid. In such a case, higher-point correlation functions---starting with the three-point one---are naturally maximally anisotropic, which makes the standard detection strategies highly inefficient and calls for a dedicated analysis of CMB data. The tensor two-point function can also be highly anisotropic, but only in the presence of sizable higher-derivative couplings.

Scalar fields appear in many theories beyond the Standard Model of particle physics. In the early universe, they are exposed to extreme conditions, including high temperature and rapid cosmic expansion. Understanding their behavior in this environment is crucial to understand the implications for cosmology. We calculate the finite temperature effective action for the field expectation value in two particularly important cases, for damped oscillations near the ground state and for scalar fields with a flat potential. We find that the behavior in both cases can in good approximation be described by a complex valued effective potential that yields Markovian equations of motion. Near the potential minimum, we recover the solution to the well-known Langevin equation. For large field values we find a very different behavior, and our result for the damping coefficient differs from the expressions frequently used in the literature. We illustrate our results in a simple scalar model, for which we give analytic approximations for the effective potential and damping coefficient. We also provide various expressions for loop integrals at finite temperature that are useful for future calculations in other models.

BICEP2 and the Keck Array are polarization-sensitive microwave telescopes that observe the cosmic microwave background (CMB) from the South Pole at degree angular scales in search of a signature of inflation imprinted as B-mode polarization in the CMB. BICEP2 was deployed in late 2009, observed for three years until the end of 2012 at 150 GHz with 512 antenna-coupled transition edge sensor bolometers, and has reported a detection of B-mode polarization on degree angular scales. The Keck Array was first deployed in late 2010 and will observe through 2016 with five receivers at several frequencies (95, 150, and 220 GHz). BICEP2 and the Keck Array share a common optical design and employ the field-proven BICEP1 strategy of using small-aperture, cold, on-axis refractive optics, providing excellent control of systematics while maintaining a large field of view. This design allows for full characterization of far-field optical performance using microwave sources on the ground. Here we describe the optical design of both instruments and report a full characterization of the optical performance and beams of BICEP2 and the Keck Array at 150 GHz.

A fraction of the first generation of stars in the early Universe may be very massive ($\gtrsim 300~\mathrm{M_\odot}$) as they form in metal-free environments. Formation of black holes from these stars can be accompanied by supermassive collapsars to produce long gamma-ray bursts of a unique type having a very high total energy ($\sim 10^{54}~\mathrm{erg}$) as recently suggested by several authors. We present new stellar evolution models of very massive Population III stars including the effect of rotation to provide theoretical constraints on super-collapsar progenitors. We find that the angular momentum condition for super-collapsar can be fulfilled if magnetic torques are ignored, in which case Eddington-Sweet circulations play the dominant role for the transport of angular momentum. We further find that the initial mass range for super-collapsar progenitors would be limited to $300~\mathrm{M_\odot} \lesssim M \lesssim 700~\mathrm{M_\odot}$. However, all of our very massive star models of this mass range end their lives as red supergiants rather than blue supergiants, in good agreement with most of the previous studies. The predicted final fate of these stars is either a jet-powered type IIP supernova or an ultra-long, relatively faint gamma-ray transient, depending on the initial amount of angular momentum.

We investigate a new scenario of dark matter production in a bouncing universe, in which dark matter was produced completely out of equilibrium in the contracting as well as expanding phase. We explore possibilities of using dark matter as a probe of the bouncing universe, focusing on the relationship between a critical temperature of the bouncing universe and the present relic abundance of dark matter.

Force-freeness of a solar magnetic field is a key to reconstructing invisible coronal magnetic structure of an emerging flux region on the Sun where active phenomena such as flares and coronal mass ejections frequently occur. We have performed magnetohydrodynamic (MHD) simulations which are adjusted to investigate force-freeness of an emerging magnetic field by using the virial theorem. Our focus is on how the force-free range of an emerging flux region develops and how it depends on the twist of a pre-emerged magnetic field. As an emerging flux region evolves, the upper limit of the force-free range continuously increases while the lower limit is asymptotically reduced to the order of a photospheric pressure scale height above the solar surface. As the twist becomes small the lower limit increases and then seems to be saturated. We also discuss the applicability of the virial theorem to an evolving magnetic structure on the Sun.

We calculate the relaxation rate of a scalar field in a plasma of other scalars and fermions with gauge interactions using thermal quantum field theory. It yields the rate of cosmic reheating and thereby determines the temperature of the "hot big bang" in inflationary cosmology. The total rate originates from various processes, including decays and inverse decays as well as Landau damping by scatterings. It involves quantum statistical effects and off-shell transport. Its temperature dependence can be highly nontrivial, making it impossible to express the reheating temperature in terms of the model parameters in a simple way. We pay special attention to the temperature dependence of the phase space due to the modified dispersion relations in the plasma. We find that it can have a drastic effect on the efficiency of perturbative reheating, which depends on the way particles in the primordial plasma interact. For some interactions thermal masses can effectively close the phase space for the dominant dissipative processes and thereby impose an upper bound on the reheating temperature. In other cases they open up new channels of dissipation, hence increase the reheating temperature. At high temperatures we find that the universe can even be heated through couplings to fermions, which are often assumed to be negligible due to Pauli-blocking. These effects may also be relevant for baryogenesis, dark matter production, the fate of moduli and in scenarios of warm inflation.

Oxygen-rich Asymptotic Giant Branch (AGB) stars can be intense emitters of SiO ($v$=1 and 2, J=1$\rightarrow$0) and H$_2$O maser lines at 43 and 22 GHz, respectively. VLBI observations of the maser emission provide a unique tool to probe the innermost layers of the circumstellar envelopes in AGB stars. Nevertheless, the difficulties in achieving astrometrically aligned \water\ and $v$=1 and $v$=2 SiO maser maps have traditionally limited the physical constraints that can be placed on the SiO maser pumping mechanism. We present phase referenced simultaneous spectral-line VLBI images for the SiO $v$=1 and $v$=2, J=1$\rightarrow$0, and H$_2$O maser emission around the AGB star R\,LMi, obtained from the Korean VLBI Network (KVN). The simultaneous multi-channel receivers of the KVN offer great possibilities for astrometry in the frequency domain. With this facility we have produced images with bona-fide absolute astrometric registration between high frequency maser transitions of different species to provide the positions of the \water\ maser emission, and the centre of the SiO maser emission, and hence reducing the uncertainty in the proper motion for R\,LMi by an order of magnitude over that from Hipparcos. This is the first successful demonstration of source frequency phase referencing for mm-VLBI spectral-line observations and also where the ratio between the frequencies is not an integer.

BICEP3 is a 550 mm-aperture refracting telescope for polarimetry of radiation in the cosmic microwave background at 95 GHz. It adopts the methodology of BICEP1, BICEP2 and the Keck Array experiments - it possesses sufficient resolution to search for signatures of the inflation-induced cosmic gravitational-wave background while utilizing a compact design for ease of construction and to facilitate the characterization and mitigation of systematics. However, BICEP3 represents a significant breakthrough in per-receiver sensitivity, with a focal plane area 5$\times$ larger than a BICEP2/Keck Array receiver and faster optics ($f/1.6$ vs. $f/2.4$). Large-aperture infrared-reflective metal-mesh filters and infrared-absorptive cold alumina filters and lenses were developed and implemented for its optics. The camera consists of 1280 dual-polarization pixels; each is a pair of orthogonal antenna arrays coupled to transition-edge sensor bolometers and read out by multiplexed SQUIDs. Upon deployment at the South Pole during the 2014-15 season, BICEP3 will have survey speed comparable to Keck Array 150 GHz (2013), and will significantly enhance spectral separation of primordial B-mode power from that of possible galactic dust contamination in the BICEP2 observation patch.

The Korean VLBI Network (KVN) is a new mm-VLBI dedicated array with capability for simultaneous observations at multiple frequencies, up to 129 GHz. The innovative multi-channel receivers present significant benefits for astrometric measurements in the frequency domain. The aim of this work is to verify the astrometric performance of the KVN using a comparative study with the VLBA, a well established instrument. For that purpose, we carried out nearly contemporaneous observations with the KVN and the VLBA, at 14/7 mm, in April 2013. The KVN observations consisted of simultaneous dual frequency observations, while the VLBA used fast frequency switching observations. We used the Source Frequency Phase Referencing technique for the observational and analysis strategy. We find that having simultaneous observations results in a superior performance for compensation of all atmospheric terms in the observables, in addition to offering other significant benefits for astrometric analysis. We have compared the KVN astrometry measurements to those from the VLBA. We find that the structure blending effects introduce dominant systematic astrometric shifts and these need to be taken into account. We have tested multiple analytical routes to characterize the impact of the low resolution effects for extended sources in the astrometric measurements. The results from the analysis of KVN and full VLBA datasets agree within 2-$\sigma$ of the thermal error estimate. We interpret the discrepancy as arising from the different resolutions. We find that the KVN provides astrometric results with excellent agreement, within 1-$\sigma$, when compared to a VLBA configuration which has a similar resolution. Therefore this comparative study verifies the astrometric performance of KVN using SFPR at 14/7 mm, and validates the KVN as an astrometric instrument.

The Korean very-long-baseline interferometry (VLBI) network (KVN) and VLBI Exploration of Radio Astrometry (VERA) Array (KaVA) is the first international VLBI array dedicated to high-frequency (23 and 43 GHz bands) observations in East Asia. Here, we report the first imaging observations of three bright active galactic nuclei (AGNs) known for their complex morphologies: 4C 39.25, 3C 273, and M 87. This is one of the initial result of KaVA early science. Our KaVA images reveal extended outflows with complex substructure such as knots and limb brightening, in agreement with previous Very Long Baseline Array (VLBA) observations. Angular resolutions are better than 1.4 and 0.8 milliarcsecond at 23 GHz and 43 GHz, respectively. KaVA achieves a high dynamic range of ~1000, more than three times the value achieved by VERA. We conclude that KaVA is a powerful array with a great potential for the study of AGN outflows, at least comparable to the best existing radio interferometric arrays.

The Korean very-long-baseline interferometry (VLBI) network (KVN) and VLBI Exploration of Radio Astrometry (VERA) Array (KaVA) is the first international VLBI array dedicated to high-frequency (23 and 43 GHz bands) observations in East Asia. Here, we report the first imaging observations of three bright active galactic nuclei (AGNs) known for their complex morphologies: 4C 39.25, 3C 273, and M 87. This is one of the initial result of KaVA early science. Our KaVA images reveal extended outflows with complex substructure such as knots and limb brightening, in agreement with previous Very Long Baseline Array (VLBA) observations. Angular resolutions are better than 1.4 and 0.8 milliarcsecond at 23 GHz and 43 GHz, respectively. KaVA achieves a high dynamic range of ~1000, more than three times the value achieved by VERA. We conclude that KaVA is a powerful array with a great potential for the study of AGN outflows, at least comparable to the best existing radio interferometric arrays.

We have carried out the first very long baseline interferometry (VLBI) imaging of 44 GHz class I methanol maser (7_{0}-6_{1}A^{+}) associated with a millimeter core MM2 in a massive star-forming region IRAS 18151-1208 with KaVA (KVN and VERA Array), which is a newly combined array of KVN (Korean VLBI Network) and VERA (VLBI Exploration of Radio Astrometry). We have succeeded in imaging compact maser features with a synthesized beam size of 2.7 milliarcseconds x 1.5 milliarcseconds (mas). These features are detected at a limited number of baselines within the length of shorter than approximately 650 km corresponding to 100 Mlambda in the uv-coverage. The central velocity and the velocity width of the 44 GHz methanol maser are consistent with those of the quiescent gas rather than the outflow traced by the SiO thermal line. The minimum component size among the maser features is ~ 5 mas x 2 mas, which corresponds to the linear size of ~ 15 AU x 6 AU assuming a distance of 3 kpc. The brightness temperatures of these features range from ~ 3.5 x 10^{8} to 1.0 x 10^{10} K, which are higher than estimated lower limit from a previous Very Large Array observation with the highest spatial resolution of ~ 50 mas. The 44 GHz class I methanol maser in IRAS 18151-1208 is found to be associated with the MM2 core, which is thought to be less evolved than another millimeter core MM1 associated with the 6.7 GHz class II methanol maser.

We have carried out the first very long baseline interferometry (VLBI) imaging of 44 GHz class I methanol maser (7_{0}-6_{1}A^{+}) associated with a millimeter core MM2 in a massive star-forming region IRAS 18151-1208 with KaVA (KVN and VERA Array), which is a newly combined array of KVN (Korean VLBI Network) and VERA (VLBI Exploration of Radio Astrometry). We have succeeded in imaging compact maser features with a synthesized beam size of 2.7 milliarcseconds x 1.5 milliarcseconds (mas). These features are detected at a limited number of baselines within the length of shorter than approximately 650 km corresponding to 100 Mlambda in the uv-coverage. The central velocity and the velocity width of the 44 GHz methanol maser are consistent with those of the quiescent gas rather than the outflow traced by the SiO thermal line. The minimum component size among the maser features is ~ 5 mas x 2 mas, which corresponds to the linear size of ~ 15 AU x 6 AU assuming a distance of 3 kpc. The brightness temperatures of these features range from ~ 3.5 x 10^{8} to 1.0 x 10^{10} K, which are higher than estimated lower limit from a previous Very Large Array observation with the highest spatial resolution of ~ 50 mas. The 44 GHz class I methanol maser in IRAS 18151-1208 is found to be associated with the MM2 core, which is thought to be less evolved than another millimeter core MM1 associated with the 6.7 GHz class II methanol maser.

We search for fast-expanding Hi shells associated with Galactic supernova remnants (SNRs) in the longitude range l \approx 32\arcdeg to 77\arcdeg using 21-cm line data from the Inner-Galaxy Arecibo L-band Feed Array (I-GALFA) Hi survey. Among the 39 known Galactic SNRs in this region, we find such Hi shells in four SNRs: W44, G54.4-0.3, W51C, and CTB 80. All four were previously identified in low-resolution surveys, and three of those (excluding G54.4-0.3) were previously studied with the Arecibo telescope. A remarkable new result, however, is the detection of Hi emission at both very high positive and negative velocities in W44 from the receding and approaching parts of the Hi expanding shell, respectively. This is the first detection of both sides of an expanding shell associated with an SNR in Hi 21-cm emission. The high-resolution I-GALFA survey data also reveal a prominent expanding Hi shell with high circular symmetry associated with G54.4-0.3. We explore the physical characteristics of four SNRs and discuss what differentiates them from other SNRs in the survey area. We conclude that these four SNRs are likely the remnants of core-collapse supernovae interacting with a relatively dense (>~1 cm^-3) ambient medium, and we discuss the visibility of SNRs in the Hi 21-cm line.

This paper presents the catalog of correlated flux densities in three ranges of baseline projection lengths of 637 sources from a 43 GHz (Q-band) survey observed with the Korean VLBI Network. Of them, 623 sources have not been observed before at Q-band with VLBI. The goal of this work in the early science phase of the new VLBI array is twofold: to evaluate the performance of the new instrument that operates in a frequency range of 22-129 GHz and to build a list of objects that can be used as targets and as calibrators. We have observed the list of 799 target sources with declinations down to -40 degrees. Among them, 724 were observed before with VLBI at 22 GHz and had correlated flux densities greater than 200 mJy. The overall detection rate is 78%. The detection limit, defined as the minimum flux density for a source to be detected with 90% probability in a single observation, was in a range of 115-180 mJy depending on declination. However, some sources as weak as 70 mJy have been detected. Of 623 detected sources, 33 objects are detected for the first time in VLBI mode. We determined their coordinates with the median formal uncertainty 20 mas. The results of this work set the basis for future efforts to build the complete flux-limited sample of extragalactic sources at frequencies 22 GHz and higher at 3/4 of the celestial sphere.

A survey of the Milky Way disk and the Magellanic System at the wavelengths of the 21-cm atomic hydrogen (HI) line and three 18-cm lines of the OH molecule will be carried out with the Australian Square Kilometre Array Pathfinder telescope. The survey will study the distribution of HI emission and absorption with unprecedented angular and velocity resolution, as well as molecular line thermal emission, absorption, and maser lines. The area to be covered includes the Galactic plane (|b|< 10deg) at all declinations south of delta = +40deg, spanning longitudes 167deg through 360deg to 79deg at b=0deg, plus the entire area of the Magellanic Stream and Clouds, a total of 13,020 square degrees. The brightness temperature sensitivity will be very good, typically sigma_T ~ 1 K at resolution 30arcsec and 1 km/s. The survey has a wide spectrum of scientific goals, from studies of galaxy evolution to star formation, with particular contributions to understanding stellar wind kinematics, the thermal phases of the interstellar medium, the interaction between gas in the disk and halo, and the dynamical and thermal states of gas at various positions along the Magellanic Stream.

We present the results of HI 21-cm line observations to explore the nature of the high-velocity (HV) HI gas at l ~ 173{\circ} . In low-resolution Hi surveys this HV gas appears as faint, wing-like, HI emission that extends to velocities beyond those allowed by Galactic rotation. We designate this feature FVW (Forbidden Velocity Wing) 172.8+1.5. Our high-resolution (3.'4) Arecibo HI observations show that FVW 172.8+1.5 is composed of knots, filaments, and ring-like structures distributed over an area a few degrees in extent. These HV HI emission features are confined within the limits of the HII complex G173+1.5, which is composed of five Sharpless HII regions distributed along a radio continuum loop of size 4{\circ}.4 {\times} 3{\circ}.4, or ~ 138 pc {\times} 107 pc, at a distance of 1.8 kpc. G173+1.5 is one of the largest star-forming regions in the outer Galaxy. We demonstrate that the HV HI gas is well correlated with the radio continuum loop and that the two seem to trace an expanding shell. The expansion velocity of the shell is large (55 km s-1) suggesting that it represents a supernova-remnant (SNR). We derive physical parameters for the shell and show these to be consistent with the object being a SNR. We also detect hot X-ray emitting gas inside the HII complex by analyzing the ROSAT all-sky X-ray background survey data. This also supports the SNR interpretation. We conclude that the HV HI gas and the X-rays are most likely the products of a supernova explosion(s) within the HII complex, possibly in a cluster that triggered the formation of these HII regions.

We report simultaneous multi-frequency observing performance at 22 and 43 GHz of the 21-m shaped-Cassegrain radio telescopes of the Korean VLBI Network (KVN). KVN is the first millimeter-dedicated VLBI network in Korea having a maximum baseline length of 480 km. It currently operates at 22 and 43 GHz and planed to operate in four frequency bands, 22, 43, 86, and 129 GHz. The unique quasioptics of KVN enable simultaneous multi-frequency observations based on efficient beam filtering and accuarate antenna-beam alignment at 22 and 43 GHz. We found that the offset of the beams is within <5 arcseconds over all pointing directions of antenna. The dual polarization, cooled HEMT receivers at 22 and 43 GHz result in receiver noise temperatures less than 40 K at 21.25-23.25 GHz and 80 K at 42.11-44.11 GHz. The pointing accuracies have been measured to be 3 arcseconds in azimuth and elevation for all antennas. The measured aperture efficiencies are 65%(K)/67%(Q), 62%(K)/59%(Q), and 66%(K)/60%(Q) for the three KVN antennas, KVNYS, KVNUS, and KVNTN, respectively. The main-beam efficiencies are measured to be 50%(K)/52%(Q), 48%(K)/50%(Q), and 50%(K)/47%(Q) for KVNYS, KVNUS, and KVNTN, respectively. The estimated Moon efficiencies are 77%(K)/90%(Q), 74%(K)/79%(Q), and 80%(K)/86%(Q) for KVNYS, KVNUS, KVNTN, respectively. The elevation dependence of the aperture efficiencies is quite flat for elevations > 20 degrees.

The I-GALFA survey is mapping all the HI in the inner Galactic disk visible to the Arecibo 305m telescope within 10 degrees of the Galactic plane (longitudes of 32 to 77 deg at b=0 deg). The survey, which will obtain 1.3 million independent spectra, became possible with the installation of the 7-beam Arecibo L-Band Feed Array (ALFA) receiver in 2004. ALFA's 3.4 arcmin resolution and tremendous sensitivity offer a great opportunity to observe the fine details of HI in the Galaxy. The I-GALFA survey began in 2008 May and will be completed in 2009 September. The data will be made publicly available when the calibrated and gridded cubes are completed. Further information on the I-GALFA project may be found at www.naic.edu/~galfa.

We describe string frame bouncing universe scenarios involving the creation and annihilation of a non-BPS D9-brane in type IIA superstring theory. We find several classes of solutions, in which the bounce is driven by the tachyon dynamics of the non-BPS brane. The metric and the dilaton are consistently described in terms of the lowest order effective action. The bounce solutions interpolate between contracting and expanding pre-big bang (or post-big bang) solutions. The singular behavior of our bounce solutions is the same as that of the pre-big bang (or post-big bang) solution. Upon adding a simple dilaton potential the asymptotic curvature singularity is removed but the dilaton still grows without bound. Such a potential may result from $\alpha'$ corrections in the open string sector.

In the present work the primordial Big-Bang Nucleosynthesis (BBN) and weakly interacting massive particle (WIMP) dark matter are discussed in a certain class of modified gravitational theories, namely $f(R) \sim R^n$ gravity. The new gravitational model is characterized by a single parameter $n$. First we determine the conditions under which the theoretical predictions for the $^{4}$He abundance are in agreement with the observations. More precisely, during BBN the physics is known and all the parameters are known. The only free parameter to be constrained is the power $n$ related to the new gravitational model. After that, for cold dark matter we use the value of $n$ determined from the BBN considerations and determine how the mass of the dark matter particle is related to the annihilation cross section in order for the cold dark matter constraint to be satisfied.

The new ekpyrotic scenario attempts to solve the singularity problem by involving violation of the null energy condition in a model which combines the ekpyrotic/cyclic scenario with the ghost condensate theory and the curvaton mechanism of production of adiabatic perturbations of metric. The Lagrangian of this theory, as well as of the ghost condensate model, contains a term with higher derivatives, which was added to the theory to stabilize its vacuum state. We found that this term may affect the dynamics of the cosmological evolution. Moreover, after a proper quantization, this term results in the existence of a new ghost field with negative energy, which leads to a catastrophic vacuum instability. We explain why one cannot treat this dangerous term as a correction valid only at small energies and momenta below some UV cut-off, and demonstrate the problems arising when one attempts to construct a UV completion of this theory.

Cosmological scenarios with k-essence are invoked in order to explain the observed late-time acceleration of the universe. These scenarios avoid the need for fine-tuned initial conditions (the "coincidence problem") because of the attractor-like dynamics of the k-essence field \phi. It was recently shown that all k-essence scenarios with Lagrangians p=L(X)/\phi^2, necessarily involve an epoch where perturbations of \phi propagate faster than light (the "no-go theorem"). We carry out a comprehensive study of attractor-like cosmological solutions ("trackers") involving a k-essence scalar field \phi and another matter component. The result of this study is a complete classification of k-essence Lagrangians that admit asymptotically stable tracking solutions, among all Lagrangians of the form p=K(\phi)L(X) . Using this classification, we select the class of models that describe the late-time acceleration and avoid the coincidence problem through the tracking mechanism. An analogous "no-go theorem" still holds for this class of models, indicating the existence of a superluminal epoch. In the context of k-essence cosmology, the superluminal epoch does not lead to causality violations. We discuss the implications of superluminal signal propagation for possible causality violations in Lorentz-invariant field theories.

We present the results of a search for faint HI 21-cm emission line wings at velocities forbidden by Galactic rotation in the Galactic plane using the Leiden/Dwingeloo HI Survey data and the HI Southern Galactic Plane Survey data. These ``forbidden-velocity wings (FVWs)'' appear as protruding excessive emission in comparison with their surroundings in limited (< 2 deg) spatial regions over velocity extent more than ~20 km/s in large-scale (l-v) diagrams. Their high-velocities imply that there should be some dynamical phenomena associated. We have identified 87 FVWs. We present their catalog, and discuss their distribution and statistical properties. We found that 85% of FVWs are not coincident with known supernova remnants (SNRs), galaxies, or high-velocity clouds. Their natures are currently unknown. We suspect that many of them are fast-moving HI shells and filaments associated with the oldest SNRs that are essentially invisible except via their HI line emission. We discuss other possible origins.

Although some 20--30,000 supernova remnants (SNRs) are expected to exist in the Milky Way, only about 230 are presently known. This implies that most SNRs are ``missing''. Recently, we proposed that small ($\simlt 1^\circ$), faint, high-velocity features seen in large-scale 21-cm line surveys of atomic hydrogen ({\sc Hi}) in the Galactic plane could be examples of such {\it missing} old SNRs. Here we report on high-resolution \schi observations of one such candidate, FVW 190.2+1.1, which is revealed to be a rapidly expanding ($\sim 80$ \kms) shell. The parameters of this shell seem only consistent with FVW 190.2+1.1 being the remnant of a SN explosion that occurred in the outermost fringes of the Galaxy some $\sim 3\times 10^5$ yr ago. This shell is not seen in any other wave band suggesting that it represents the oldest type of SNR, that which is essentially invisible except via its \schi line emission. FVW 190.2+1.1 is one of a hundred "forbidden-velocity wings" (FVWs) recently identified in the Galactic plane, and our discovery suggests that many of these are likely to be among the oldest SNRs. We discuss the possible link between FVWs and fast-moving atomic clouds in the Galaxy.

We construct an anomaly free supersymmetric U(1)' model with a secluded U(1)'-breaking sector. We study the one-loop effective potential at finite temperature, and show that there exists a strong enough first order electroweak phase transition for electroweak baryogenesis (EWBG) because of the large trilinear term $A_h h S H_d H_u$ in the tree-level Higgs potential. Unlike in the MSSM, the lightest stop can be very heavy. We consider the non-local EWBG mechanism in the thin wall regime, and find that within uncertainties the observed baryon number can be generated from the $\tau$ lepton contribution, with the secluded sector playing an essential role. The chargino and neutralino contributions and the implications for the Z' mass and electric dipole moments are briefly discussed.

We have searched for HI 21 cm line emission from shocked atomic gas associated with southern supernova remnants (SNRs) using data from the Southern Galactic Plane Survey. Among the 97 sources studied, we have detected 10 SNRs with high-velocity HI emission confined to the SNR. The large velocity and the spatial confinement suggest that the emission is likely from the gas accelerated by the SN blast wave. We also detected 22 SNRs which show HI emission significantly brighter than the surrounding regions over a wide ($>10$\kms) velocity interval. The association with these SNRs is less certain. We present the parameters and maps of the excess emission in these SNRs. We discuss in some detail the ten individual SNRs with associated high-velocity HI emission.

We estimate the number of old, radiative supernova remnants (SNRs) detectable in HI 21-cm emission line in the Galaxy. According to our result, the contamination due to the Galactic background HI emission limits the number of visible SNRs to 270, or 9% of the total HI SNRs. The telescope sensitivity further limits the number. We compare the result with observations to find that the observed number (< 25) of HI SNRs is much less than the expected. A plausible explanation is that previous observational studies, which were made toward the SNRs identified mostly in radio continuum, missed most of the HI SNRs because they are too faint to be visible in radio continuum. We propose that the faint, extended HI 21-cm emission line wings protruding from the Galactic background HI emission in large-scale (l, v) diagrams could be possible candidates for HI SNRs, although our preliminary result shows that their number is considerably less than the expected in the inner Galaxy. We conclude that a possible explanation for the small number of \schi SNRs in the inner Galaxy is that the interstellar space there is largely filled with a very tenuous gas as in the three-phase ISM model.

We briefly summarize the survey of HI 21 cm emission lines to search for shocked atomic gas associated with Galactic supernova remnants (SNRs) in the southern sky. For G347.3-0.5, we discuss the distance to the SNR and the implications of the HI results.