The "ATOMS," standing for {\it ALMA Three-millimeter Observations of Massive Star-forming regions}, survey has observed 146 active star forming regions with ALMA Band 3, aiming to systematically investigate the spatial distribution of various dense gas tracers in a large sample of Galactic massive clumps, to study the roles of stellar feedback in star formation, and to characterize filamentary structures inside massive clumps. In this work, the observations, data analysis, and example science of the "ATOMS" survey are presented, using a case study for the G9.62+0.19 complex. Toward this source, some transitions, commonly assumed to trace dense gas, including CS $J = 2-1$, HCO$^+$ $J = 1-0$ and HCN $J = 1-0$, are found to show extended gas emission in low density regions within the clump; less than 25\% of their emission is from dense cores. SO, CH$_3$OH, H$^{13}$CN and HC$_3$N show similar morphologies in their spatial distributions and reveal well the dense cores. Widespread narrow SiO emission is present (over $\sim$1 pc), which may be caused by slow shocks from large--scale colliding flows or H{\sc ii} regions. Stellar feedback from an expanding H{\sc ii} region has greatly reshaped the natal clump, significantly changed the spatial distribution of gas, and may also account for the sequential high-mass star formation in the G9.62+0.19 complex. The ATOMS survey data can be jointly analyzed with other survey data, e.g., "MALT90", "Orion B", "EMPIRE", "ALMA\_IMF", and "ALMAGAL", to deepen our understandings of "dense gas" star formation scaling relations and massive proto-cluster formation.

We report studies of the relationships between the total bolometric luminosity ($L_{\rm bol}$ or $L_{\rm TIR}$) and the molecular line luminosities of $J=1-0$ transitions of H$^{13}$CN, H$^{13}$CO$^+$, HCN, and HCO$^+$ with data obtained from ACA observations in the "ATOMS" survey of 146 active Galactic star forming regions. The correlations between $L_{\rm bol}$ and molecular line luminosities $L'_{\rm mol}$ of the four transitions all appear to be approximately linear. Line emission of isotopologues shows as large scatters in $L_{\rm bol}$-$L'_{\rm mol}$ relations as their main line emission. The log($L_{\rm bol}$/$L'_{\rm mol}$) for different molecular line tracers have similar distributions. The $L_{\rm bol}$-to-$L'_{\rm mol}$ ratios do not change with galactocentric distances ($R_{\rm GC}$) and clump masses ($M_{\rm clump}$). The molecular line luminosity ratios (HCN-to-HCO$^+$, H$^{13}$CN-to-H$^{13}$CO$^+$, HCN-to-H$^{13}$CN and HCO$^+$-to-H$^{13}$CO$^+$) all appear constant against $L_{\rm bol}$, dust temperature ($T_{\rm d}$), $M_{\rm clump}$ and $R_{\rm GC}$. Our studies suggest that both the main lines and isotopologue lines are good tracers of the total masses of dense gas in Galactic molecular clumps. The large optical depths of main lines do not affect the interpretation of the slopes in star formation relations. We find that the mean star formation efficiency (SFE) of massive Galactic clumps in the "ATOMS" survey is reasonably consistent with other measures of the SFE for dense gas, even those using very different tracers or examining very different spatial scales.

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.

In the episodic accretion scenario, a large fraction of the protostellar mass accretes during repeated and large bursts of accretion. Since outbursts on protostars are typically identified at specific wavelengths, interpreting these outbursts requires converting this change in flux to a change in total luminosity. The Class I young stellar object EC 53 in the Serpens Main cloud has undergone repeated increases in brightness at 850 $\mu$m that are likely caused by bursts of accretion. In this study, we perform two- and three-dimensional continuum radiative transfer modeling to quantify the internal luminosity rise in EC 53 that corresponds to the factor of $\sim$1.5 enhancement in flux at 850 $\mu$m. We model the spectral energy distribution and radial intensity profile in both the quiescent and outburst phases. The internal luminosity in the outburst phase is $\sim 3.3$ times brighter than the luminosity in the quiescent phase. The radial intensity profile analysis demonstrates that the detected sub-mm flux variation of EC 53 comes from the heated envelope by the accretion burst. We also find that the role of external heating of the EC 53 envelope by the interstellar radiation field is insignificant.

In the third catalog of active galactic nuclei detected by the $Fermi$ Large Area Telescope Clean (3LAC) sample, there are 402 blazars candidates of uncertain type (BCU). The proposed analysis will help to evaluate the potential optical classification flat spectrum radio quasars (FSRQs) versus BL Lacertae (BL Lacs) objects of BCUs, which can help to understand which is the most elusive class of blazar hidden in the Fermi sample. By studying the 3LAC sample, we found some critical values of $\gamma$-ray photon spectral index ($\Gamma_{\rm ph}$), variability index (VI) and radio flux (${\rm F_R}$) of the sources separate known FSRQs and BL Lac objects. We further utilize those values to defined an empirical "high-confidence" candidate zone that can be used to classify the BCUs. Within such a zone ($\Gamma_{\rm ph}<2.187$, log${\rm F_R}<2.258$ and ${ \rm logVI <1.702}$), we found that 120 BCUs can be classified BL Lac candidates with a higher degree of confidence (with a misjudged rate $<1\%$). Our results suggest that an empirical "high confidence" diagnosis is possible to distinguish the BL Lacs from the Fermi observations based on only on the direct observational data of $\Gamma_{\rm ph}$, VI and ${\rm F_R}$.

We scrutinise the widely studied minimal scotogenic model of dark matter and radiative neutrino mass from the requirement of a strong first order electroweak phase transition (EWPT) and observable gravitational waves at future planned space based experiments. The scalar DM scenario is similar to inert scalar doublet extension of standard model where a strong first order EWPT favours the low mass regime of DM, a narrow fine tuned region currently allowed from collider and direct detection bounds. In the fermion DM scenario, we get newer region of parameter space which favours strong first order EWPT as the restriction on mass ordering within inert scalar doublet gets relaxed. While such leptophilic fermion DM remains safe from stringent direct detection bounds, newly allowed low mass regime of charged scalar can leave tantalising signatures at colliders. While we get such new region of parameter space satisfying DM relic, strong first order EWPT with detectable gravitational waves, light neutrino mass and other relevant constraints, we also improve upon previous analysis in similar model by incorporating appropriate resummation effects in effective finite temperature potential.

In a new release of the fourth \emph{Fermi} Large Area Telescope source catalog (4FGL) in September (URL: https://fermi.gsfc.nasa.gov/ssc/data/access/lat/8yr_catalog/), 5065 $\gamma$-ray sources are reported, including 3207 active galactic nuclei (AGNs), 239 pulsars, 190 other sources and 1429 unassociated sources. We employ two different supervised machine learning classifiers, combined with the direct observation parameters given by the 4FGL fits (URL: https://fermi.gsfc.nasa.gov/ssc/data/access/lat/8yr_catalog/gll_psc_v20.fit), to search for sources potentially classified as AGNs and pulsars in the 1429 unassociated sources. In order to reduce the error caused by the large difference in the sizes of samples, we divide the classification process into two separate steps in order to identify the AGNs and the pulsars. First, we select the identified AGNs from all of the samples, and then select the identified pulsars from the remaining. Using the 4FGL sources associated or identified as AGNs, pulsars, and other sources with the features selected through the K-S test, we trained, optimized, and tested our classifier models. Then, the models are applied to classify the 1429 unassociated sources. According to the direct calculation results of the two classifiers, we show the sensitivity, specificity, accuracy in each step, and the class of unassociated sources given by each classifier. The accuracy obtained in the first step is approximately $95\%$; in the second step, the obtained overall accuracy is approximately $80\%$. Combining the results of the two classifiers, we predict that there are 674 AGN-type candidates, 86 pulsar-type candidates, 177 other types of $\gamma$-ray candidates, and 492 of uncertain type.

We propose a new mechanism for generating matter-antimatter asymmetry via the interference of tree-level diagrams only. We first derive a general result that a nonzero $CP$-asymmetry can be generated via at least two sets of interfering tree-level diagrams involving either $2\rightarrow2$ or $1\rightarrow {\cal N}$ (with ${\cal N}\geq3$) processes. We illustrate this point in a simple TeV-scale extension of the Standard Model with an inert Higgs doublet and right-handed neutrinos, along with an electroweak triplet scalar field. The imaginary part needed for the required $CP$-asymmetry comes from the trilinear coupling of the inert doublet with the triplet scalar. Small Majorana neutrino masses are generated by both the scotogenic and type-II seesaw mechanisms. The real part of the neutral component of the inert doublet serves as a cold dark matter candidate. The evolutions of the dark matter relic density and the baryon asymmetry are intimately related in this scenario.

We show that the weakly interacting massive particle (WIMP) scenario of proton-philic spin-dependent inelastic dark matter can still provide a viable explanation of the observed DAMA effect in compliance with the constraints from other experiments. We also show that, although the COSINE-100 collaboration has recently tested the DAMA effect using the same target material, for the time being the comparison between DAMA and COSINE-100 still depends on the particle-physics model.

The recently published fourth Fermi Large Area Telescope source catalog (4FGL) reports 5065 gamma-ray sources in terms of direct observational gamma-ray properties. Among the sources, the largest population is the Active Galactic Nuclei (AGN), which consists of 3137 blazars, 42 radio galaxies, and 28 other AGNs. The blazar sample comprises 694 flat-spectrum radio quasars (FSRQs), 1131 BL Lac-type objects (BL Lacs), and 1312 blazar candidates of an unknown type (BCUs). The classification of blazars is difficult using optical spectroscopy given the limited knowledge with respect to their intrinsic properties, and the limited availability of astronomical observations. To overcome these challenges, machine learning algorithms are being investigated as alternative approaches. Using the 4FGL catalog, a sample of 3137 Fermi blazars with 23 parameters is systematically selected. Three established supervised machine learning algorithms (random forests (RFs), support vector machines (SVMs), artificial neural networks (ANNs)) are employed to general predictive models to classify the BCUs. We analyze the results for all of the different combinations of parameters. Interestingly, a previously reported trend the use of more parameters leading to higher accuracy is not found. Considering the least number of parameters used, combinations of eight, 12 or 10 parameters in the SVM, ANN, or RF generated models achieve the highest accuracy (Accuracy $\simeq$ 91.8\%, or $\simeq$ 92.9\%). Using the combined classification results from the optimal combinations of parameters, 724 BL Lac type candidates and 332 FSRQ type candidates are predicted; however, 256 remain without a clear prediction.

We analyze the DAMA/Libra-phase2 modulation result using a basis of 16 effective operators describing the WIMP interaction with photons, gluons and quarks up to mass dimension seven. For each operator we fix the effective theory at the scale of 2 GeV and parametrize WIMP-quark interactions in terms of two independent couplings for up-type and down-type quarks. We discuss the connection with the non-relativistic limit of the effective theory in terms of operators invariant by Galilean transformations, and the impact of the ensuing momentum-dependent Wilson coefficients due to long-range interactions or light meson poles. Most relativistic operators yield a good fit of the DAMA modulation effect, although for parameters excluded by the constraints from XENON1T and PICO60.

We discuss the possibility of realising a two-component dark matter (DM) scenario where the two DM candidates differ from each other by virtue of their production mechanism in the early universe. One of the DM candidates is thermally generated in a way similar to the weakly interacting massive particle (WIMP) paradigm where the DM abundance is governed by its freeze-out while the other candidate is produced only from non-thermal contributions similar to freeze-in mechanism. We discuss this in a minimal extension of the standard model where light neutrino masses arise radiatively in a way similar to the scotogenic models with DM particles going inside the loop. The lepton asymmetry is generated at the same time from WIMP DM annihilations as well as partially from the mother particle for non-thermal DM. This can be achieved while satisfying the relevant experimental bounds, and keeping the scale of leptogenesis or the thermal DM mass as low as 3 TeV, well within present experimental reach. In contrast to the TeV scale thermal DM mass, the non-thermal DM can be as low as a few keV, giving rise to the possibility of a sub-dominant warm dark matter (WDM) component that can have interesting consequences on structure formation. The model also has tantalizing prospects of being detected at ongoing direct detection experiments as well as the ones looking for charged lepton flavour violating process like $\mu \rightarrow e \gamma$.

We started a multi-scale analysis of G202.3+2.5, an intertwined filamentary region of Monoceros OB1. In Paper I, we examined the distributions of dense cores and protostars and found enhanced star formation (SF) activity in the junction region of the filaments. In this second paper, we aim to unveil the connections between the core and filament evolutions, and between the filament dynamics and the global evolution of the cloud. We characterise the gas dynamics and energy balance using Herschel and WISE observations and molecular tracers observed with the IRAM 30m and TRAO 14m telescopes. The velocity field of the cloud is examined and velocity-coherent structures are put in perspective with the cloud environment. Two main velocity components (VCs) are revealed, well separated in the north and merged around the location of intense N2H+ emission where Paper I found the peak of SF activity. The relative position of the two VCs along the sightline, and the velocity gradient in N2H+ emission imply that the VCs have been undergoing collision for ~10^5 yrs. The dense gas where N2H+ is detected is interpreted as the compressed region between the two filaments, which corresponds to a high mass inflow rate of ~1e-3 Msun/yr and possibly leads to an increase in its SF efficiency. We identify a protostar in the junction region that possibly powers two crossed intermittent outflows. We show that the HII region around the nearby cluster NCG 2264 is still expanding and its role in the collision is examined. However, we cannot rule out the idea that the collision arises mostly from the global collapse of the cloud. The (sub-)filament-scale observables examined in this paper reveal a collision between G202.3+2.5 sub-structures and its probable role in feeding the cores in the junction region. One must now characterise the cloud morphology, its fragmentation, and magnetic field, all at high resolution.

Current theories and models attempt to explain star formation globally, from core scales to giant molecular cloud scales. A multi-scale observational characterisation of an entire molecular complex is necessary to constrain them. We investigate star formation in G202.3+2.5, a ~10x3 pc sub-region of the Monoceros OB1 cloud with a complex morphology harbouring interconnected filamentary structures. We aim to connect the evolution of cores and filaments in G202.3+2.5 with the global evolution of the cloud and to identify the engines of the cloud dynamics. In this first paper, the star formation activity is evaluated by surveying the distributions of dense cores and protostars, and their evolutionary state, as characterised using both infrared observations from the Herschel and WISE telescopes and molecular line observations with the IRAM 30-m telescope. We find ongoing star formation in the whole cloud, with a local peak in star formation activity around the centre of G202.3+2.5 where a chain of massive cores (10-50 Msun) forms a massive ridge (>150 Msun). All evolutionary stages from starless cores to Class II protostars are found in G202.3+2.5, including a possibly starless, large column density (8x10^{22} cm-2), and massive (52 Msun) core. All the core-scale observables examined in this paper point to an enhanced star formation activity centred on the junction between the three main branches of the ramified structure of G202.3+2.5. This suggests that the increased star-formation activity results from the convergence of these branches. To further investigate the origin of this enhancement, it is now necessary to extend the analysis to larger scales, in order to examine the relationship between cores, filaments and their environment. We address these points through the analysis of the dynamics of G202.3+2.5 in a joint paper.

Assuming for Weakly Interacting Massive Particles (WIMPs) a Maxwellian velocity distribution in the Galaxy we provide an assessment of the sensitivity of existing Dark Matter (DM) direct detection (DD) experiments to operators up to dimension 7 of the relativistic effective field theory describing dark matter interactions with quarks and gluons . In particular we focus on a systematic approach, including an extensive set of experiments and large number of couplings, both exceeding for completeness similar analyses in the literature. The relativistic effective theory requires to fix one coupling for each quark flavor, so in principle for each different combination the bounds should be recalculated starting from direct detection experimental data. To address this problem we propose an approximate model-independent procedure that allows to directly calculate the bounds for any combination of couplings in terms of model-independent limits on the Wilson coefficients of the non-relativistic theory expressed in terms of the WIMP mass and of the neutron-to-proton coupling ratio $c^n/c^p$. We test the result of the approximate procedure against that of a full calculation, and discuss its possible pitfalls and limitations. We also provide a simple interpolating interface in Python that allows to apply our method quantitatively.

This white paper gives a brief summary of Galactic continuum surveys with the next generation 850 micron camera at the James Clerk Maxwell Telescope (JCMT) in the next decade. This new camera will have mapping speeds at least 10-20 times faster than the present SCUBA-2 camera, and will enable deep (<10 mJy/beam) and extensive continuum surveys such as a wider (-5 deg.<l<240 deg., |b|<2-5 deg.) Galactic Plane survey, a larger Gould Belt survey, and a follow-up survey of all Planck compact objects visible from the northern hemisphere. These surveys will provide a complete census of molecular clouds, filaments, and dense cores across the Galaxy, vital for studying star formation in various environments.

Assuming a standard Maxwellian for the WIMP velocity distribution, we obtain the bounds from null WIMP search results of 59.5 days of COSINE-100 data on the DAMA/LIBRA-phase2 modulation effect within the context of the non-relativistic effective theory of WIMP-nucleus scattering. Here, we systematically assume that one of the effective operators allowed by Galilean invariance dominates in the effective Hamiltonian of a spin-1/2 dark matter (DM) particle. We find that, although DAMA/LIBRA and COSINE-100 use the same sodium-iodide target, the comparison of the two results still depends on the particle-physics model. This is mainly due to two reasons: i) the WIMP signal spectral shape; ii) the expected modulation fractions, when the upper bound on the time-averaged rate in COSINE-100 is converted into a constraint on the annual modulation component in DAMA/LIBRA. We find that the latter effect is the dominant one. For several effective operators the expected modulation fractions are larger than in the standard spin-independent or spin-dependent interaction cases. As a consequence, compatibility between the modulation effect observed in DAMA/LIBRA and the null result from COSINE-100 is still possible for several non-relativistic operators. At low WIMP masses such relatively high values of the modulation fractions arise because COSINE-100 is mainly sensitive to WIMP-sodium scattering events, due to the higher threshold compared to DAMA/LIBRA. A next COSINE analysis is expected to have a full sensitivity for the 5$\sigma$ region of DAMA/LIBRA.

In the early stages of star formation, a protostar is deeply embedded in an optically thick envelope such that it is not directly observable. Variations in the protostellar accretion rate, however, will cause luminosity changes that are reprocessed by the surrounding envelope and are observable at submillimeter wavelengths. We searched for submillimeter flux variability toward 12 Planck Galactic Cold Clumps detected by the James Clerk Maxwell Telescope (JCMT)-SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE) survey. These observations were conducted at 850 um using the JCMT/SCUBA-2. Each field was observed three times over about 14 months between 2016 April and 2017 June. We applied a relative flux calibration and achieved a calibration uncertainty of ~ 3.6% on average. We identified 136 clumps across 12 fields and detected four sources with flux variations of ~ 30%. For three of these sources, the variations appear to be primarily due to large-scale contamination, leaving one plausible candidate. The flux change of the candidate may be associated with low- or intermediate-mass star formation assuming a distance of 1.5 kpc, although we cannot completely rule out the possibility that it is a random deviation. Further studies with dedicated monitoring would provide a better understanding of the detailed relationship between submillimeter flux and accretion rate variabilities while enhancing the search for variability in star-forming clumps farther away than the Gould Belt.

Young stars exhibit variability due to changes in the gas accretion rate onto them, an effect that should be quite significant in the early stages of their formation. As protostars are embedded within their natal cloud, this variability may only be inferred through long wavelength observations. We perform radiative transfer simulations of young stellar objects (YSOs) formed in hydrodynamical simulations, varying the structure and luminosity properties in order to estimate the long-wavelength, sub-mm and mm, variations of their flux. We find that the flux increase due to an outburst event depends on the protostellar structure and is more prominent at sub-mm wavelengths than at mm wavelengths; e.g. a factor of 40 increase in the luminosity of the young protostar leads to a flux increase of a factor of 10 at 250 micron but only a factor of 2.5 at 1.3 mm. We find that the interstellar radiation field dilutes the flux increase but that this effect may be avoided if resolution permits the monitoring of the inner regions of a YSO, where the heating is primarily due to protostellar radiation. We also confirm that the bolometric temperature and luminosity of outbursting protostars may result in an incorrect classification of their evolutionary stage.

Accretion onto protostars may occur in sharp bursts. Accretion bursts during the embedded phase of young protostars are probably most intense, but can only be inferred indirectly through long-wavelength observations. We perform radiative transfer calculations for young stellar objects (YSOs) formed in hydrodynamic simulations to predict the long wavelength, sub-mm and mm, flux responses to episodic accretion events, taking into account heating from the young protostar and from the interstellar radiation field. We find that the flux increase due to episodic accretion events is more prominent at sub-mm wavelengths than at mm wavelengths; e.g. a factor of ~570 increase in the luminosity of the young protostar leads to a flux increase of a factor of 47 at 250 micron but only a factor of 10 at 1.3 mm. Heating from the interstellar radiation field may reduce further the flux increase observed at longer wavelengths. We find that during FU Ori-type outbursts the bolometric temperature and luminosity may incorrectly classify a source as a more evolved YSO, due to a larger fraction of the radiation of the object being emitted at shorter wavelengths

The flat spectrum radio quasar 4C 38.41 showed a significant increase of its radio flux density during the period 2012 March - 2015 August which correlates with gamma-ray flaring activity. Multi-frequency simultaneous VLBI observations were conducted as part of the interferometric monitoring of gamma-ray bright active galactic nuclei (iMOGABA) program and supplemented with additional monitoring observations at various bands across the electromagnetic spectrum. The epochs of the maxima for the two largest gamma-ray flares coincide with the ejection of two respective new VLBI components and the evolution of the physical properties seem to be in agreement with the shock-in-jet model. Derived synchrotron self absorption magnetic fields, of the order of 0.1 mG, do not seem to dramatically change during the flares, and are much smaller, by a factor 10,000, than the estimated equipartition magnetic fields, indicating that the source of the flare may be associated with a particle dominated emitting region.

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 report the first high-resolution (sub-arcminute) large-scale mapping $^{12}$CO and $^{13}$CO observations of the molecular clouds associated with the giant outer Galaxy HII region CTB~102 (KR 1). These observations were made using a newly commissioned receiver system on the 13.7-m radio telescope at the Taeduk Radio Astronomy Observatory. Our observations show that the molecular clouds have a spatial extent of $60 \times 35$ pc and a total mass of $10^{4.8} - 10^{5.0}$ M$_\odot$. Infrared data from WISE and 2MASS were used to identify and classify the YSO population associated with ongoing star formation activity within the molecular clouds. We directly detect 18 class I/class II YSOs and six transition disk objects. Moving away from the HII region, there is an age/class gradient consistent with sequential star formation. The infrared and molecular-line data were combined to estimate the star formation efficiency (SFE) of the entire cloud as well as the SFE for various sub-regions of the cloud. We find that the overall SFE is between $\sim5 - 10$%, consistent with previous observations of giant molecular clouds. One of the sub-regions, region 1a, is a clear outlier, with a SFE of 17 $-$ 35% on a 5 pc spatial scale. This high SFE is more typical for much smaller (sub-pc scale) star-forming cores, and we think region 1a is likely an embedded massive protocluster.

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.

A long standing problem in the study of Active Galactic Nuclei (AGNs) is that the observed VLBI core is in fact a blending of the actual AGN core (classically defined by the $\tau=1$ surface) and the upstream regions of the jet or optically thin emitting region flows. This blending may cause some biases towards the observables of the core, such as its flux density, size or brightness temperature, which may lead to misleading interpretation of the derived quantities and physics. We study the effects of such blending under the view of the Korean VLBI Network (KVN) for a sample of AGNs at 43 GHz by comparing their observed properties with observations with the Very Large Baseline Array (VLBA). Our results suggest that the observed core sizes are a factor $\sim11$ larger than these of VLBA, which is similar to the factor expected by considering the different resolutions of the two facilities. We suggest the use of this factor to consider blending effects in KVN measurements. Other parameters, such as flux density or brightness temperature, seem to possess a more complicated dependence.

Assuming for Weakly Interacting Massive Particles (WIMPs) a Maxwellian velocity distribution in the Galaxy we explore in a systematic way the relative sensitivity of an extensive set of existing and projected Dark Matter (DM) direct detection experiments to each of the 14 couplings that parameterize the most general non-relativistic (NR) effective Hamiltonian allowed by Galilean invariance for a contact interaction driving the elastic scattering off nuclei of WIMPs of spin 1/2. We perform our analysis in terms of two free parameters: the WIMP mass $m_{\chi}$ and the ratio between the WIMP-neutron and the WIMP-proton couplings $c^n/c^p$. We include the modified signal spectral shape due to non-standard interactions when it is needed in the determination of the bound, such as in the case of background subtraction or of the application of the optimal-interval method. For each coupling, in the $m_{\chi}$-$c^n/c^p$ plane we provide contour plots of the most stringent 90 % C.L. bound on the WIMP-nucleon cross section and show the experiment providing it. We also introduce NRDD_constraints, a simple interpolating code written in Python that allows to obtain the numerical value of the bound as a function of the WIMP mass $m_{\chi}$ and of the coupling ratio $c^n/c^p$ for each NR coupling. We find that 9 experiments out of the 14 present Dark Matter searches considered in our analysis provide the most stringent bound on some of the effective couplings for a given choice of $(m_{\chi},c^n/c^p)$: this is evidence of the complementarity of different target nuclei and/or different combinations of count-rates and energy thresholds when the search of DM is extended to a wide range of possible interactions.

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.

There are still some important unanswered questions about the detailed particle acceleration and escape occurring during the quiescent epoches. As a result, the particle distribution that is adopted in the blazar quiescent spectral model have numerous unconstrained shapes. To help remedy this problem, we introduce a analytical particle transport model to reproduce quiescent broadband spectral energy distribution of blazar. In this model, the exact electron distribution is solved from a generalized transport equation that contains the terms describing first-order and secondary-order \emph{Fermi} acceleration, escape of particle due to both the advection and spatial diffusion, energy losses due to synchrotron emission and inverse-Compton scattering of an assumed soft photon field. We suggest that the advection is a significant escape mechanism in blazar jet. We find that in our model the advection process tends to harden the particle distribution, which enhances the high energy components of resulting synchrotron and synchrotron self-Comptom spectrum from jet. Our model is able to roughly reproduce the observed spectra of extreme BL Lac object 1ES 0414+009 with reasonable assumptions about the physical parameters.

In the third catalog of active galactic nuclei detected by the Fermi-LAT (3LAC) Clean Sample, there are 402 blazars candidates of uncertain type (BCU). Due to the limitations of astronomical observation or intrinsic properties, it is difficult to classify blazars using optical spectroscopy. The potential classification of BCUs using machine learning algorithms is essential. Based on the 3LAC Clean Sample, we collect 1420 Fermi blazars with 8 parameters of {\gamma}-ray photon spectral index, radio flux, flux density, curve significance, the integral photon flux in 100 to 300 MeV, 0.3 to 1 GeV, 10 to 100 GeV and variability index. Here, we apply 4 different supervised machine learning (SML) algorithms (\emph{Decision trees, Random forests, support vector machines and Mclust Gaussian finite mixture models}) to evaluate the classification of BCUs based on the direct observational properties. All the 4 methods can perform exceedingly well with a more accuracy and can effective forecast the classification of Fermi BCUs. The evaluating results show the results of these methods (SML) are valid and robust, where, about 1/4 sources are FSRQs and 3/4 are BL Lacs in 400 BCUs, which are consistent with some other recent results. Although a number of factors influence the accuracy of SML, the results are stable at a fixed ratio 1:3 between FSRQs and BL Lacs, which suggests that the SML can provides an effective method to evaluate the potential classification of BCUs. Among the 4 methods, Mclust Gaussian Mixture Modelling has the highest accuracy for our training sample (4/5, seed=123).

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 show that the Weakly Interacting Massive Particle scenario of proton-philic spin-dependent inelastic Dark Matter (pSIDM) can still provide a viable explanation of the observed DAMA modulation amplitude in compliance with the constraints from other experiments after the release of the DAMA/LIBRA-phase2 data and including the recent bound from COSINE-100, that uses the same $NaI$ target of DAMA. The pSIDM scenario provided a viable explanation of DAMA/LIBRA--phase1 both for a Maxwellian WIMP velocity distribution and in a halo-independent approach. At variance with DAMA/LIBRA-phase1, for which the modulation amplitudes showed an isolated maximum at low energy, the DAMA/LIBRA-phase2 spectrum is compatible to a monotonically decreasing one. Moreover, due to its lower threshold, it is sensitive to WIMP-iodine interactions at low WIMP masses. Due to the combination of these two effects pSIDM can now explain the yearly modulation observed by DAMA/LIBRA only when the WIMP velocity distribution departs from a standard Maxwellian. In this case the WIMP mass $m_{\chi}$ and mass splitting $\delta$ fall in the approximate ranges 7 GeV $\lesssim m_{\chi}\lesssim$ 17 GeV and 18 keV$\lesssim\delta\lesssim$29 keV. The recent COSINE-100 bound is naturally evaded in the pSDIM scenario due to its large expected modulation fractions.

We re-examine the case of anapole dark matter as an explanation for the DAMA annual modulation in light of the DAMA/LIBRA-phase2 results and improved upper limits from other DM searches. If the WIMP velocity distribution is assumed to be a Maxwellian, anapole dark matter is unable to provide an explanation of the DAMA modulation compatible with the other searches. Nevertheless, anapole dark matter provides a better fit to the DAMA-phase2 modulation data than an isoscalar spin-independent interaction, due to its magnetic coupling with sodium targets. A halo-independent analysis shows that explaining the DAMA modulation above 2 keVee in terms of anapole dark matter is basically impossible in face of the other null results, while the DAMA/LIBRA-phase2 modulation measurements below 2 keVee are marginally allowed. We conclude that in light of current measurements, anapole dark matter does not seem to be a viable explanation for the totality of the DAMA modulation.

How stellar feedback from high-mass stars (e.g., H{\sc ii} regions) influences the surrounding interstellar medium and regulates new star formation is still unclear. To address this question, we observed the G9.62+0.19 complex in 850 $\mu$m continuum with the JCMT/POL-2 polarimeter. An ordered magnetic field has been discovered in its youngest clump, the G9.62 clump. The magnetic field strength is determined to be $\sim$1 mG. Magnetic field plays a larger role than turbulence in supporting the clump. However, the G9.62 clump is still unstable against gravitational collapse even if thermal, turbulent, and magnetic field support are taken into account all together. The magnetic field segments in the outskirts of the G9.62 clump seem to point toward the clump center, resembling a dragged-in morphology, indicating that the clump is likely undergoing magnetically-regulated global collapse. However, The magnetic field in its central region is aligned with the shells of the photodissociation regions (PDRs) and is approximately parallel to the ionization (or shock) front, indicating that the magnetic field therein is likely compressed by the expanding H{\sc ii} regions that formed in the same complex.

Recent Fermi-Large Area Telescope (LAT) light curves indicate an active $\gamma$-ray state spanning about five months from 2016 June to 2016 October in the BL Lac object 1749+096 (OT 081). During this period, we find two notable $\gamma$-ray events: an exceptionally strong outburst followed by a significant enhancement (local peak). In this study, we analyze multi-waveband light curves (radio, optical, X-ray, and $\gamma$-ray) plus very-long baseline interferometry (VLBI) data to investigate the nature of the $\gamma$-ray events. The $\gamma$-ray outburst coincides with flux maxima at longer wavelengths. We find a spectral hardening of the $\gamma$-ray photon index during the $\gamma$-ray outburst. The photon index shows a transition from a softer-when-brighter to a harder-when-brighter trend at around 1.8 $\times$ $10^{-7}$ ph cm$^{-2}$ s$^{-1}$. We see indication that both the $\gamma$-ray outburst and the subsequent enhancement precede the propagation of a polarized knot in a region near the VLBI core. The highest polarized intensity, 230\,mJy, and an electric vector position angle rotation, by $\sim$32$^{\circ}$, are detected about 12 days after the $\gamma$-ray outburst. We conclude that both $\gamma$-ray events are caused by the propagation of a disturbance in the mm-wave core.

The DAMA/LIBRA collaboration has recently released updated results from their search for the annual modulation signal expected from Dark Matter (DM) scattering in their NaI detectors. We have fitted the updated DAMA result for the modulation amplitudes in terms of a Weakly Interacting Massive Particle (WIMP) signal, parameterizing the interaction with nuclei in terms of the most general effective Lagrangian for a WIMP particle spin up to 1/2, systematically assuming dominance of one of the 14 possible interaction terms, and assuming for the WIMP velocity distribution a standard Maxwellian. We find that most of the couplings of the non-relativistic effective Hamiltonian can provide a better fit compared to the standard Spin Independent interaction case, and with a reduced fine-tuning of the three parameters (WIMP mass, WIMP-nucleon effective cross-section and ratio between the WIMP-neutron and the WIMP-proton couplings). Moreover, effective models for which the cross section depends explicitly on the WIMP incoming velocity can provide a better fit of the DAMA data at large values of $m_{\chi}$ compared to the standard velocity-independent cross-section due to a different phase of the modulation amplitudes. All the best fit solutions are in tension with exclusion plots of both XENON1T and PICO60.

We propose a WIMP baryogensis achieved by the annihilation of non-thermally produced WIMPs from decay of heavy particles, which can result in low reheating temerature. Dark matter (DM) can be produced non-thermally during a reheating period created by the decay of long-lived heavy particle, and subsequently re-annihilate to lighter particles even after the thermal freeze-out. The re-annihilation of DM provides the observed baryon asymmetry as well as the correct relic density of DM. We investigate how wahout effects can affect the generation of the baryon asymmetry and study a model suppressing them. In this scenario, we find that DM can be heavy enough and its annihilation cross section can also be larger than that adopted in the usual thermal WIMP baryogenesis.

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.

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.

Based on the 850 $\mu$m dust continuum data from SCUBA-2 at James Clerk Maxwell Telescope (JCMT), we compare overall properties of Planck Galactic Cold Clumps (PGCCs) in the $\lambda$ Orionis cloud to those of PGCCs in the Orion A and B clouds. The Orion A and B clouds are well known active star-forming regions, while the $\lambda$ Orionis cloud has a different environment as a consequence of the interaction with a prominent OB association and a giant Hii region. PGCCs in the $\lambda$ Orionis cloud have higher dust temperatures ($Td=16.13\pm0.15$ K) and lower values of dust emissivity spectral index ($ \beta=1.65\pm0.02$) than PGCCs in the Orion A (Td=13.79$\pm 0.21$K, $\beta=2.07\pm0.03$) and Orion B ($Td=13.82\pm0.19$K, $\beta=1.96\pm0.02$) clouds. We find 119 sub-structures within the 40 detected PGCCs and identify them as cores. Of total 119 cores, 15 cores are discovered in the $\lambda$ Orionis cloud, while 74 and 30 cores are found in the Orion A and B clouds, respectively. The cores in the $\lambda$ Orionis cloud show much lower mean values of size R=0.08 pc, column density N(H2)=$(9.5\pm1.2) \times 10^{22}$ cm$^{-2}$, number density n(H2)=$(2.9 \pm 0.4)\times10^{5}$ cm$^{-3}$, and mass $M_{core}$=$1.0\pm0.3$ M$_{\odot}$ compared to the cores in the Orion A (R=0.11pc, $N(H2)=(2.3\pm0.3) \times 10^{23}$ cm$^{-2}$, n(H2)=$(3.8\pm0.5) \times 10^{5}$cm$^{-3}$, and $M_{core}$=$2.4 \pm 0.3$ M$_{\odot}$) and Orion B (R=0.16pc, N(H2)=$(3.8 \pm 0.4) \times 10^{23}$cm$^{-2}$, n(H2)=$(15.6\pm1.8)\times10^{5}$ cm$^{-3}$, and $M_{core}$= $2.7\pm0.3$ M$_{\odot}$) clouds. These core properties in the $\lambda$ Orionis cloud can be attributed to the photodissociation and external heating by the nearby Hii region, which may prevent the PGCCs from forming gravitationally bound structures and eventually disperse them. These results support the idea of negative stellar feedback on core formation.

We study the linear polarization of the radio cores of eight blazars simultaneously at 22, 43, and 86 GHz with observations obtained by the Korean VLBI Network (KVN) in three epochs between late 2016 and early 2017 in the frame of the Plasma-physics of Active Galactic Nuclei (PAGaN) project. We investigate the Faraday rotation measure (RM) of the cores; the RM is expected to increase with observing frequency if core positions depend on frequency due to synchrotron self-absorption. We find a systematic increase of RMs at higher observing frequencies in our targets. The RM--$\nu$ relations follow power-laws with indices distributed around 2, indicating conically expanding outflows serving as Faraday rotating media. Comparing our KVN data with contemporaneous optical polarization data from the Steward Observatory for a few sources, we find indication that the increase of RM with frequency saturates at frequencies of a few hundreds GHz. This suggests that blazar cores are physical structures rather than simple $\tau=1$ surfaces. A single region, e.g. a recollimation shock, might dominate the jet emission downstream of the jet launching region. We detect a sign change in the observed RMs of CTA 102 on a time scale of $\approx$1 month, which might be related to new superluminal components emerging from its core undergoing acceleration/deceleration and/or bending. We see indication for quasars having higher core RMs than BL Lac objects, which could be due to denser inflows/outflows in quasars.

The flat spectrum radio quasar 1633+382 (4C~38.41) showed a significant increase of its radio flux density during the period 2012 March - 2015 August which correlates with gamma-ray flaring activity. Multi-frequency simultaneous VLBI observations were conducted as part of the interferometric monitoring of gamma-ray bright active galactic nuclei (iMOGABA) program and supplemented with additional radio monitoring observations with the OVRO 40m telescope, the Boston University VLBI program, and the Submillimeter Array. The epochs of the maxima for the two largest gamma-ray flares coincide with the ejection of two respective new VLBI components. Analysis of the spectral energy distribution indicates a higher turnover frequency after the flaring events. The evolution of the flare in the turnover frequency - turnover flux density plane probes the adiabatic losses in agreement with the shock-in-jet model. Derived synchrotron self absorption magnetic fields, of the order of 0.1mG, do not seem to dramatically change during the flares, and are much smaller, by a factor $10^4$, than the estimated equipartition magnetic fields, indicating that the source of the flare may be associated with a particle dominated emitting region.

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.

The curvature of the $\gamma$-ray spectrum in blazars may reflect the intrinsic distribution of the emitting electron distribution, which will further give some information on the possible acceleration and cooling processes in the emitting region. The $\gamma$-ray spectra of Fermi blazars are normally fitted either by a single power-law (PL) or a log-normal (call Logarithmic Parabola, LP) form. The possible reason for this differnece is not unclear. We statistically explore this issue based on the different observational properties of 1419 Fermi blazars in the 3LAC Clean sample. We find that the $\gamma$-ray flux (100 MeV-100 GeV) and variability index follow bimodal distributions for PL and LP blazars, where $\gamma$-ray flux and variability index show {a positive correlation}. However, the distributions of the $\gamma$-ray luminosity and redshift follow a unimodal distribution. Our results suggest that the bimodal distribution of $\gamma$-ray flux for LP and PL blazars may be not intrinsic and all blazars may have an intrinsic curved $\gamma$-ray spectrum and the PL spectrum is just caused by the fitting effect due to the less photons.

3C 84 (NGC 1275) is a well-studied mis-aligned Active Galactic Nucleus (AGN), which has been active in Gamma rays since at least 2008. We have monitored the source at four wavelengths (14 mm, 7 mm, 3 mm and 2 mm) using the Korean VLBI network (KVN) since 2013 as part of the interferometric monitoring of $\gamma$-ray bright AGN (iMOGABA) program. 3C 84 exhibits bright radio emission both near the central supermassive black hole (SMBH) feature known as C1 and from a moving feature located to the south known as C3. Other facilities have also detected these short-term variations above a slowly rising trend at shorter wavelengths, such as in Gamma ray and 1 mm total intensity light-curves. We find that the variations in the $\gamma$ rays and 1 mm total intensity light-curves are correlated, with the $\gamma$ rays leading and lagging the radio emission. Analysis of the 2 mm KVN data shows that both the Gamma rays and 1 mm total intensity short-term variations are better correlated with the SMBH region than C3, likely placing the short-term variations in C1. We interpret the emission as being due to the random alignment of spatially separated emission regions. We place the slowly rising trend in C3, consistent with previous results. Additionally, we report that since mid-2015, a large mm-wave radio flare has been occurring in C3, with a large Gamma ray flare coincident with the onset of this flare at all radio wavelengths.

We study the centimeter- to millimeter-wavelength synchrotron spectrum of the core of the radio galaxy M87 at $\lesssim0.8\,{\rm mas}~\sim110R_{s}$ spatial scales using four years of fully simultaneous, multi-frequency VLBI data obtained by the Korean VLBI Network (KVN). We find a core spectral index $\alpha$ of $\gtrsim-0.37$ ($S\propto \nu^{+\alpha}$) between 22GHz and 129GHz. By combining resolution-matched flux measurements from the Very Long Baseline Array (VLBA) at 15GHz and taking the Event Horizon Telescope (EHT) 230GHz core flux measurements in epochs 2009 and 2012 as lower limits, we find evidence of a nearly flat core spectrum across 15GHz and 129GHz, which could naturally connect the 230GHz VLBI core flux. The extremely flat spectrum is a strong indication that the jet base does not consist of a simple homogeneous plasma, but of inhomogeneous multi-energy components, with at least one component with the turn-over frequency $\gtrsim100$GHz. The spectral shape can be qualitatively explained if both the strongly (compact, optically thick at $>$100GHz) and the relatively weakly magnetized (more extended, optically thin at $<$100GHz) plasma components are colocated in the footprint of the relativistic jet.

We analyze results from the first eighteen months of monthly sub-mm monitoring of eight star-forming regions in the JCMT Transient Survey. In our search for stochastic variability in 1643 bright peaks, only the previously identified source, EC53, shows behavior well above the expected measurement uncertainty. Another four sources, two disks and two protostars, show moderately-enhanced standard deviations in brightness, as expected for stochastic variables. For the two protostars, this apparent variability is the result of single epochs that are much brighter than the mean. In our search for secular brightness variations that are linear in time, we measure the fractional brightness change per year for 150 bright peaks, fifty of which are protostellar. The ensemble distribution of slopes is well fit by a normal distribution with sigma ~ 0.023. Most sources are not rapidly brightening or fading in the sub-mm. Comparison against time-randomized realizations shows that the width of the distribution is dominated by the uncertainty in the individual brightness measurements of the sources. A toy model for secular variability reveals that an underlying Gaussian distribution of linear fractional brightness change sigma = 0.005 would be unobservable in the present sample, whereas an underlying distribution with sigma = 0.02 is ruled out. Five protostellar sources, 10% of the protostellar sample, are found to have robust secular measures deviating from a constant flux. The sensitivity to secular brightness variations will improve significantly with a larger time sample, with a factor of two improvement expected by the conclusion of our 36-month survey.

Atmospheric propagation effects at millimeter wavelengths can significantly alter the phases of radio signals and reduce the coherence time, putting tight constraints on high frequency Very Long Baseline Interferometry (VLBI) observations. In previous works, it has been shown that non-dispersive (e.g. tropospheric) effects can be calibrated with the frequency phase transfer (FPT) technique. The coherence time can thus be significantly extended. Ionospheric effects, which can still be significant, remain however uncalibrated after FPT as well as the instrumental effects. In this work, we implement a further phase transfer between two FPT residuals (i.e. so-called FPT-square) to calibrate the ionospheric effects based on their frequency dependence. We show that after FPT-square, the coherence time at 3 mm can be further extended beyond 8~hours, and the residual phase errors can be sufficiently canceled by applying the calibration of another source, which can have a large angular separation from the target (>20 deg) and significant temporal gaps. Calibrations for all-sky distributed sources with a few calibrators are also possible after FPT-square. One of the strengths and uniqueness of this calibration strategy is the suitability for high-frequency all-sky survey observations including very weak sources. We discuss the introduction of a pulse calibration system in the future to calibrate the remaining instrumental effects and allowing the possibility of imaging the source structure at high frequencies with FPT-square, where all phases are fully calibrated without involving any additional sources.

We present multi-frequency simultaneous VLBI radio observations of the flat spectrum radio quasar 1633+382 (4C~38.41) as part of the interferometric monitoring of gamma-ray bright active galactic nuclei (iMOGABA) program combined with additional observations in radio, optical, X-rays and $\gamma-$rays carried out between the period 2012 March - 2015 August. The monitoring of this source reveals a significant long-lived increase in its activity since approximately two years in the radio bands, which correlates with a similar increase on all other bands from sub-millimeter to $\gamma-$rays. A significant correlation is also found between radio fluxes and simultaneous spectral indices during this period. The study of the discrete correlation function (DCF) indicates time lags smaller than the $\sim40$ days uncertainties among both radio bands and also high-energy bands, and a time lag of $\sim$70 days, with $\gamma-$rays leading radio. We interpret that the high-energy and radio fluxes are arising from different emitting regions, located at $1\pm12$ and $40\pm13$ pc from the central engine respectively.

We fit the spectral energy distributions (SEDs) of 46 GeV - TeV BL Lac objects in the frame of leptonic one-zone synchrotron self-Compton (SSC) model and investigate the physical properties of these objects. We use the Markov Chain Monte Carlo (MCMC) method to obtain the basic parameters, such as magnetic field (B), the break energy of the relativistic electron distribution ($\gamma'_{\rm{b}}$) and the electron energy spectral index. Based on the modeling results, we support the following scenarios on GeV-TeV BL Lac objects: (1) Some sources have large Doppler factors, implying other radiation mechanism should be considered. (2) Comparing with FSRQs, GeV-TeV BL Lac objects have weaker magnetic field and larger Doppler factor, which cause the ineffective cooling and shift the SEDs to higher bands. Their jet powers are around $4.0\times 10^{45}~\rm{ erg\cdot s}^{-1}$, comparing with radiation power, $5.0\times 10^{42}~\rm{ erg\cdot s}^{-1}$, indicating that only a small fraction of jet power is transformed into the emission power. (3) For some BL Lacs with large Doppler factors, their jet components could have two substructures, e.g., the fast core and the slow sheath. For most GeV-TeV BL Lacs, Kelvin-Helmholtz instabilities are suppressed by their higher magnetic fields, leading few micro-variability or intro-day variability in the optical bands. (4) Combined with a sample of FSRQs, an anti-correlation between the peak luminosity $L_{\rm {pk}}$ and the peak frequency $\nu_{\rm {pk}}$ is obtained, favoring the blazar sequence scenario. In addition, an anti-correlation between the jet power $P_{\rm {jet}}$ and the break Lorentz factor $\gamma_{\rm {b}}$ also supports the blazar sequence.

We present a study of the inexplicit connection between radio jet activity and gamma-ray emission of BL Lacertae (BL Lac; 2200+420). We analyze the long-term millimeter activity of BL Lac via interferometric observations with the Korean VLBI Network (KVN) obtained at 22, 43, 86, and 129 GHz simultaneously over three years (from January 2013 to March 2016); during this time, two gamma-ray outbursts (in November 2013 and March 2015) can be seen in gamma-ray light curves obtained from Fermi observations. The KVN radio core is optically thick at least up to 86 GHz; there is indication that it might be optically thin at higher frequencies. To first order, the radio light curves decay exponentially over the time span covered by our observations, with decay timescales of 411+/-85 days, 352+/-79 days, 310+/-57 days, and 283+/-55 days at 22, 43, 86, and 129 GHz, respectively. Assuming synchrotron cooling, a cooling time of around one year is consistent with magnetic field strengths B~2microT and electron Lorentz factors gamma~10,000. Taking into account that our formal measurement errors include intrinsic variability and thus over-estimate the statistical uncertainties, we find that the decay timescale tau scales with frequency nu like tau~nu^{-0.2}. This relation is much shallower than the one expected from opacity effects (core shift), but in agreement with the (sub-)mm radio core being a standing recollimation shock. We do not find convincing radio flux counterparts to the gamma$ray outbursts. The spectral evolution is consistent with the `generalized shock model' of Valtaoja et al. (1992). A temporary increase in the core opacity and the emergence of a knot around the time of the second gamma-ray event indicate that this gamma-ray outburst might be an `orphan' flare powered by the `ring of fire' mechanism.

The V380 Ori NE bipolar outflow was imaged in the SiO and CO J = 1 - 0 lines, and dense cores in L1641 were observed in the 2.0-0.89 mm continuum. The highly collimated SiO jet shows point-symmetric oscillation patterns in both position and velocity, which suggests that the jet axis is precessing and the driving source may belong to a non-coplanar binary system. By considering the position and velocity variabilities together, accurate jet parameters were derived. The protostellar system is viewed nearly edge-on, and the jet has a flow speed of 35 km/s and a precession period of 1600 years. The CO outflow length gives a dynamical timescale of 6300 years, and the protostar must be extremely young. The inferred binary separation of 6-70 au implies that this protobinary system may have been formed through the disk instability process. The continuum spectra of L1641 dense cores indicate that the emission comes from dust, and the fits with modified blackbody functions give emissivity power indices of beta = 0.3-2.2. The emissivity index shows a positive correlation with the molecular line width, but no strong correlation with bolometric luminosity or temperature. V380 Ori NE has a particularly low value of beta = 0.3, which tentatively suggests the presence of millimeter-sized dust grains. Because the dust growth takes millions of years, much longer than the protostellar age, this core may have produced large grains in the starless core stage. HH 34 MMS and HH 147 MMS also have low emissivity indices.

Investigating variability at the earliest stages of low-mass star formation is fundamental in understanding how a protostar assembles mass. While many simulations of protostellar disks predict non-steady accretion onto protostars, deeper investigation requires robust observational constraints on the frequency and amplitude of variability events characterised across the observable SED. In this study, we develop methods to robustly analyse repeated observations of an area of the sky for submillimetre variability in order to determine constraints on the magnitude and frequency of deeply embedded protostars. We compare \mbox{850 $\mu$m} JCMT Transient Survey data with archival JCMT Gould Belt Survey data to investigate variability over 2-4 year timescales. Out of 175 bright, independent emission sources identified in the overlapping fields, we find 7 variable candidates, 5 of which we classify as \textit{Strong} and the remaining 2 as \textit{Extended} to indicate the latter are associated with larger-scale structure. For the \textit{Strong} variable candidates, we find an average fractional peak brightness change per year of |4.0|\% yr$^{-1}$ with a standard deviation of $2.7\%\mathrm{\:yr}^{-1}$. In total, 7\% of the protostars associated with \mbox{850 $\mu$m} emission in our sample show signs of variability. Four of the five \textit{Strong} sources are associated with a known protostar. The remaining source is a good follow-up target for an object that is anticipated to contain an enshrouded, deeply embedded protostar. In addition, we estimate the \mbox{850 $\mu$m} periodicity of the submillimetre variable source, EC 53, to be \mbox{567 $\pm$ 32 days} based on the archival Gould Belt Survey data.

During the protostellar phase of stellar evolution, accretion onto the star is expected to be variable, but this suspected variability has been difficult to detect because protostars are deeply embedded. In this paper, we describe a sub-mm luminosity burst of the Class I protostar EC 53 in Serpens Main, the first variable found during our dedicated JCMT/SCUBA-2 monitoring program of eight nearby star-forming regions. EC 53 remained quiescent for the first 6 months of our survey, from February to August 2016. The sub-mm emission began to brighten in September 2016, reached a peak brightness of $1.5$ times the faint state, and has been decaying slowly since February 2017. The change in sub-mm brightness is interpreted as dust heating in the envelope, generated by a luminosity increase of the protostar of a factor of $\ge 4$. The 850~$\mu$m lightcurve resembles the historical $K$-band lightcurve, which varies by a factor of $\sim 6$ with a 543 period and is interpreted as accretion variability excited by interactions between the accretion disk and a close binary system. The predictable detections of accretion variability observed at both near-infrared and sub-mm wavelengths make the system a unique test-bed, enabling us to capture the moment of the accretion burst and to study the consequences of the outburst on the protostellar disk and envelope.

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.

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.

We employ a single-zone leptonic jet model, with synchrotron, synchrotron self-Compton (SSC) and external Compton (EC) process, to reproduce the quasi-simultaneous multi-wavelength spectral energy distributions in active and quiescent states of the narrow-line gamma-ray-loud radio source GB 1310+487. In the case of the EC process, the external seed photons from both broad line region (BLR) and dust torus are considered by assuming that the gamma-ray emission region is located at the outside boundary of the BLR and inside the dust torus. Comparing the energy density of external photon fields $U_{\rm BLR}$ obtained by model fitting with that constrained from the BLR observations. We find that the location of the gamma-ray emitting region of GB 1310+487 can be tightly constrained at the outer edge of the BLR (the dissipation distance of the $\gamma$-ray emission region from central black hole $r_{\rm diss} \sim {\rm~a~few~times~of}~R_{\rm BLR}$). The ratio of magnetic energy and emitting-electron energy in the radiation blob ($\epsilon_{B}=L_{B}/L_{\rm e}$) is gradually increased from Flare 1, Flare 2 to Post-Flare, where the magnetic energy increase and matter energy decrease. These results suggest that the conversion of the magnetic field and the matter (radiation electrons) energy and the location of the $\gamma$-ray emission region (or ambient photon field) may play an important role in different radiation states of GB 1310+487.

We present results of single-epoch very long baseline interferometry (VLBI) observations of gamma-ray bright active galactic nuclei (AGNs) using the Korean VLBI Network (KVN) at 22, 43, 86, and 129~GHz bands, which are part of a KVN key science program, Interferometric Monitoring of Gamma-ray Bright AGNs (iMOGABA). We selected a total of 34 radio-loud AGNs of which 30 sources are gamma-ray bright AGNs with flux densities of $>6\times10^{-10}$~ph~cm$^{-2}$~s$^{-1}$. Single-epoch multi-frequency VLBI observations of the target sources were conducted during a 24-hr session on 2013 November 19 and 20. All observed sources were detected and imaged at all frequency bands with or without a frequency phase transfer technique which enabled the imaging of 12 faint sources at 129~GHz, except for one source. Many of the target sources are resolved on milliarcsecond scales, yielding a core-jet structure with the VLBI core dominating the synchrotron emission on the milliarcsecond scale. CLEAN flux densities of the target sources are 0.43-28~Jy, 0.32-21~Jy, 0.18-11~Jy, and 0.35-8.0~Jy in the 22, 43, 86, and 129~GHz bands, respectively.

We present the results of a broadband radio-to-GeV observing campaign organized to get a better understanding of the radiation processes responsible for the $\gamma$-ray flares observed in 3C 279. The total intensity and polarization observations of the source were carried out between December 28, 2013 and January 03, 2014 using the Fermi-LAT, Swift-XRT, Swift-UVOT, and KVN telescopes. A prominent flare observed in the optical/near-UV passbands was found to be correlated with a concurrent $\gamma$-ray flare at a confidence level $>$95$\%$, which suggests a co-spatial origin of the two. Moreover, the flaring activity in the two regimes was accompanied by no significant spectral variations. A peak in the X-ray light curve coincides with the peaks of the fractional polarization curves at 43 and 86 GHz radio bands. No prominent variation was noticed for the total intensity and the electric vector position angle (EVPA) observations at radio bands during this period. We noticed a possible hint of steepening of the radio spectrum with an increase in percentage polarization, which suggests that the radio polarization variations could be simply due to a spectral change. In a simple scenario, the correlated optical/$\gamma$-ray flares could be caused by the same population of emitting particles. The coincidence of the increase in radio polarization with the X-ray flux supports the picture that X-rays are produced via inverse-Compton scattering of radio photons. The observed fractional variability for the $\gamma$-ray flare $\sim$0.23 does not exceed that in the optical regime, which is inconsistent with what we usually observe for 3C 279; it could be due to different dependencies of the magnetic field and the external radiation field energy density profiles along the jet.

We present an ordinary case in the momentum diffusion equation for the electron spectrum evolution and investigate the energy spectra and time dependent properties of flare in Mrk 421 in the frame of a time-dependent one-zone synchrotron self-Compton model. In this model, electrons are accelerated to extra-relativistic energy through the strong magnetic turbulence and evolve with time, and non-thermal photons are produced by both synchrotron and inverse Comtpon scattering off synchrotron photons. Moreover, non-thermal photons during a pre-flare are produced by the relativistic electrons in the steady state and those during a flare are produced by the electrons whose injection rate is changed during some time interval. We apply the model to the energy spectra and time-dependent properties of flares in Mrk 421 by reproducing the pre-burst spectrum of the source and varying the injection rate in the Bohm diffusion (q=1) and the hard sphere approximation (q=2) case, respectively. Our results show that Bohm diffusion case leads to hard photon spectra, and hard sphere case approximation seems to reproduce the energy spectra and the time dependent properties of flare still better.

We investigate MeV-GeV $\gamma$-ray outbursts in 3C 273 in the frame of a time-dependent one-zone synchrotron self-Compton (SSC) model. In this model, electrons are accelerated to extra-relativistic energy through the stochastic particle acceleration and evolve with the time, nonthermal photons are produced by both synchrotron and inverse Compton scattering of synchrotron photons. Moreover, nonthermal photons during a quiescent are produced by the relativistic electrons in the steady state and those during a outburst are produced by the electrons whose injection rate is changed at some time interval. We apply the model to two exceptionally luminous $\gamma$-ray outbursts observed by the Fermi-LAT from 3C 273 in September, 2009 and obtain the multi-wavelength spectra during the quiescent and during the outburst states, respectively. Our results show that the time-dependent properties of outbursts can be reproduced by adopting the appropriate injection rate function of the electron population.

Context. Very high-energy (VHE) $\gamma$-ray measurements of distant TeV blazars can be nicely explained by TeV spectra induced by ultra high-energy cosmic rays. Aims. We develop a model for a plausible origin of hard spectra in distant TeV blazars. Methods. In the model, the TeV emission in distant TeV blazars is dominated by two mixed components. The first is the internal component with the photon energy around 1 TeV produced by inverse Compton scattering of the relativistic electrons on the synchrotron photons (SSC) with a correction for extragalactic background light absorbtion and the other is the external component with the photon energy more than 1 TeV produced by the cascade emission from high-energy protons propagating through intergalactic space. Results. Assuming suitable model parameters, we apply the model to observed spectra of distant TeV blazars of 1ES 0229+200. Our results show that 1) the observed spectrum properties of 1ES 0229+200, especially the TeV $\gamma$-ray tail of the observed spectra, could be reproduced in our model and 2) an expected TeV $\gamma$-ray spectrum with photon energy $>$1 TeV of 1ES 0229+200 should be comparable with the 50-hour sensitivity goal of the Cherenkov Telescope Array (CTA) and the differential sensitivity curve for the one-year observation with the Large High Altitude Air Shower Observatory (LHAASO). Conclusions. We argue that strong evidence for the Bethe-Heitler cascades along the line of sight as a plausible origin of hard spectra in distant TeV blazars could be obtained from VHE observations with CTA, LHAASO, HAWC, and HiSCORE.

PKS 1424+240 is a distant very high energy gamma-ray BL Lac object with redshift $z=0.601$. It was found that pure synchrotron self-Compton (SSC) process normally need extreme input parameters (e.g., very low magnetic field intensity and extraordinarily large Doppler factor) to explain its multi-wavelength spectral energy distributions (SEDs). To avoid the extreme model parameters, different models have been proposed (e.g., two-zone SSC model or lepto-hadronic model). In this work, we employ the traditional one-zone leptonic model after including a weak external Compton component to re-explore the simultaneous multi-wavelength SEDs of PKS 1424+240 in both high (2009) and low (2013) states. We find that the input parameters of magnetic field and Doppler factor are roughly consistent with those of other BL Lacs if a weak external photon field from either broad line region (BLR) or the dust torus. However, the required energy density of seed photons from BLR or torus is about 3 orders of magnitude less than that constrained in luminous quasars (e.g., flat-spectrum radio quasars, FSRQs). This result suggests that the BLR/torus in BL Lacs is much weaker than that of luminous FSRQs (but not fully disappear), and the inverse-Compton of external photons from BLR/torus may still play a role even in high synchrotron peaked blazars.

We present results of long-term multi-wavelength polarization observations of the powerful blazar 3C~279 after its $\gamma$-ray flare on 2013~December 20. We followed up this flare with single-dish polarization observations using two 21-m telescopes of the Korean VLBI Network. Observations carried out weekly from 2013~December~25 to 2015~January~11, at 22~GHz, 43~GHz, 86~GHz simultaneously, as part of the Monitoring Of GAmma-ray Bright AGN (MOGABA) program. We measured 3C~279 total flux densities of 22--34~Jy at 22~GHz, 15--28~Jy (43~GHz), and 10--21~Jy (86~GHz), showing mild variability of $\leq 50\,\%$ over the period of our observations. The spectral index between 22~GHz and 86~GHz ranged from $-0.13$ to $-0.36$. Linear polarization angles were 27$^{\circ}$--38$^{\circ}$, 30$^{\circ}$--42$^{\circ}$, and 33$^{\circ}$--50$^{\circ}$ at 22~GHz, 43~GHz, and 86~GHz, respectively. The degree of linear polarization was in the range of 6--12\,\%, and slightly decreased with time at all frequencies. We investigated Faraday rotation and depolarization of the polarized emission at 22--86~GHz, and found Faraday rotation measures (RM) of $-300$ to $-1200$~rad~m$^{-2}$ between 22~GHz and 43~GHz, and $-800$ to $-5100$~rad~m$^{-2}$ between 43~GHz and 86~GHz. The RM values follow a power law with a mean power law index $a$ of $2.2$, implying that the polarized emission at these frequencies travels through a Faraday screen in or near the jet. We conclude that the regions emitting polarized radio emission may be different from the region responsible for the 2013 December $\gamma$-ray flare and are maintained by the dominant magnetic field perpendicular to the direction of the radio jet at milliarcsecond scales.

We report first results from KVN and VERA Array (KaVA) VLBI observations obtained in the frame of our Plasma-physics of Active Galactic Nuclei (PAGaN) project. We observed eight selected AGN at 22 and 43 GHz in single polarization (LCP) between March 2014 and April 2015. Each source was observed for 6 to 8 hours per observing run to maximize the $uv$ coverage. We obtained a total of 15 deep high-resolution images permitting the identification of individual circular Gaussian jet components and three spectral index maps of BL Lac, 3C 111 and 3C 345 from simultaneous dual-frequency observations. The spectral index maps show trends in agreement with general expectations -- flat core and steep jets -- while the actual value of the spectral index for jets shows indications for a dependence on AGN type. We analyzed the kinematics of jet components of BL Lac and 3C 111, detecting superluminal proper motions with maximum apparent speeds of about $5c$. This constrains the lower limits of the intrinsic component velocities to $\sim0.98c$ and the upper limits of the angle between jet and line of sight to $\sim$20$\deg$. In agreement with global jet expansion, jet components show systematically larger diameters $d$ at larger core distances $r$, following the global relation $d\approx0.2r$, albeit within substantial scatter.

The Interferometric Monitoring of Gamma-ray Bright Active galactic nuclei (iMOGABA) program provides not only simultaneous multifrequency observations of bright gamma-ray detected active galactic nuclei (AGN), but also covers the highest Very Large Baseline Interferometry (VLBI) frequencies ever being systematically monitored, up to 129 GHz. However, observation and imaging of weak sources at the highest observed frequencies is very challenging. In the second paper in this series, we evaluate the viability of the frequency phase transfer technique to iMOGABA in order to obtain larger coherence time at the higher frequencies of this program (86 and 129 GHz) and image additional sources that were not detected using standard techniques. We find that this method is applicable to the iMOGABA program even under non-optimal weather conditions.

We report the first detection of linearly polarized emission at an observing wavelength of 350 $\mu$m from the radio-loud Active Galactic Nucleus 3C 279. We conducted polarization observations for 3C 279 using the SHARP polarimeter in the Caltech Submillimeter Observatory on 2014 March 13 and 14. For the first time, we detected the linear polarization with the degree of polarization of 13.3\%$\pm$3.4\% ($3.9\sigma$) and the Electric Vector Position Angle (EVPA) of 34.7$^\circ\pm5.6^\circ$. We also observed 3C~279 simultaneously at 13, 7, and 3.5 mm in dual polarization with the Korean very long baseline interferometry (VLBI) Network on 2014 March 6 (single dish) and imaged in milliarcsecond (mas) scales at 13, 7, 3.5, and 2.3 mm on March 22 (VLBI). We found that the degree of linear polarization increases from 10\% to 13\% at 13 mm to 350 $\mu$m and the EVPAs at all observing frequencies are parallel within $<10^\circ$ to the direction of the jet at mas scale, implying that the integrated magnetic fields are perpendicular to the jet in the innermost regions. We also found that the Faraday rotation measures RM are in a range of $-6.5\times10^2 \sim -2.7\times10^3$ rad m$^{-2}$ between 13 and 3.5 mm, and are scaled as a function of wavelength: $|{\rm RM}|\propto\lambda^{-2.2}$. These results indicate that the mm and sub-mm polarization emission are generated in the compact jet within 1~mas scale and affected by a Faraday screen in or in the close proximity of the jet.

In the light of the new IceCube data on the (yet unidentified) astrophysical neutrino flux in the PeV and sub-PeV range, we present an update on the status of decaying dark matter interpretation of the events. In particular, we develop further the angular distribution analysis and discuss the perspectives for diagnostics. By performing various statistical tests (maximum likelihood, Kolmogorov-Smirnov and Anderson-Darling tests) we conclude that currently the data show a mild preference (below the two sigma level) for the angular distribution expected from dark matter decay vs. the isotropic distribution foreseen for a conventional astrophysical flux of extragalactic origin. Also, we briefly develop some general considerations on heavy dark matter model building and on the compatibility of the expected energy spectrum of decay products with the IceCube data, as well as with existing bounds from gamma-rays. Alternatively, assuming that the IceCube data originate from conventional astrophysical sources, we derive bounds on both decaying and annihilating dark matter for various final states. The lower limits on heavy dark matter lifetime improve by up to an order of magnitude with respect to existing constraints, definitively making these events---even if astrophysical in origin---an important tool for astroparticle physics studies.

We fit the (quasi-)simultaneous multi-waveband spectral energy distributions (SEDs) for a sample of low-synchrotron-peaked (LSP) blazars with a one-zone leptonic model. The seed photons that predominantly come from broad line region (BLR) and infrared (IR) molecular torus are considered respectively in external Compton process. We find that the modeling with IR seed photons is systematically better than that with BLR photons based on a $\chi^2$ test, which suggest that $\gamma$-ray emitting region most possibly stay outside the BLR. The minimum electron Lorentz factor, $\gamma_{\rm min}$, is constrained from the modeling for these LSP blazars with good soft X-ray data (ranges from 5 to 160 with a median value of 55), which plays a key role in jet power estimation. Assuming one-to-one ratio of proton and electron, we find that the jet power for LSP blazars is systematically higher than that of FR II radio galaxies at given 151 MHz radio luminosity, $L_{\rm 151MHz}$, even though FR IIs are regarded as same as LSP blazars in unification scheme except the jet viewing angle. The possible reason is that there are some $e^{\pm}$ pairs in the jet of these blazars. If this is the case, we find the number density of $e^{\pm}$ pairs should be several times higher than that of $e^{-}-p$ pairs by assuming the jet power is the same for LSP blazars and FR IIs at given $L_{\rm 151MHz}$.

We started two observing programs with the Korean VLBI Network (KVN) monitoring changes in the flux density and polarization of relativistic jets in gamma-ray bright AGNs simultaneously at 22, 43, 86, 129 GHz. One is a single-dish weekly-observing program in dual polarization with KVN 21-m diameter radio telescopes beginning in 2011 May. The other is a VLBI monthly-observing program with the three-element VLBI network at an angular resolution range of 1.0--9.2 mas beginning in 2012 December. The monitoring observations aim to study correlation of variability in gamma-ray with that in radio flux density and polarization of relativistic jets when they flare up. These observations enable us to study the origin of the gamma-ray flares of AGNs.

We present an analysis of HCO+ (J=3-2) and H13CO+ (J=3-2) observations of the massive (M ~ 20 Msun) submm/IR source IRAS 01202+6133 located on the periphery of the HII region KR 120 (Sh 2-187). The HCO+ line profile has a classic blue-asymmetric shape with the optically thin H13CO+ line peaking at the position expected if the HCO+ line arises from a combination of self-absorption and infall motion. We have modified existing analytic radiative transfer models to allow for the fitting of submm/mm line profiles that have both self-absorption features and optically thin wings and applied these models to our HCO+ spectrum of IRAS 01202+6133. We conclude that it is a young Class I YSO with a substantial envelope undergoing slow infall and having some outflow motions. The young age of the HII region rules out a "collect and collapse" scenario. While we cannot eliminate the possibility that IRAS 01202+6133 formed spontaneously at its current location, considering its early evolutionary state and its proximity to the HII region we think that the formation of IRAS 01202+6133 was triggered by the expansion of KR 120 (Sh 2-187).

We consider the Higgs inflation in the extension of the Standard Model with two Higgs doublets coupled to gravity non-minimally. In the presence of an approximate global U(1) symmetry in the Higgs sector, both radial and angular modes of neutral Higgs bosons drive inflation where large non-Gaussianity is possible from appropriate initial conditions on the angular mode. We also discuss the case with single-field inflation for which the U(1) symmetry is broken to a Z_2 subgroup. We show that inflationary constraints, perturbativity and stability conditions restrict the parameter space of the Higgs quartic couplings at low energy in both multi- and single-field cases. Focusing on the inert doublet models where Z_2 symmetry remains unbroken at low energy, we show that the extra neutral Higgs boson can be a dark matter candidate consistent with the inflationary constraints. The doublet dark matter is always heavy in multi-field inflation while it can be light due to the suppression of the co-annihilation in single-field inflation. The implication of the extra quartic couplings on the vacuum stability bound is also discussed in the light of the recent LHC limits on the Higgs mass.

We present a comprehensive abundance analysis of 20 elements for 16 new low-metallicity stars from the Chemical Abundances of Stars in the Halo (CASH) project. The abundances have been derived from both Hobby-Eberly Telescope High Resolution Spectrograph snapshot spectra (R~15,000) and corresponding high-resolution (R~35,000) Magellan MIKE spectra. The stars span a metallicity range from [Fe/H] from -2.9 to -3.9, including four new stars with [Fe/H]<-3.7. We find four stars to be carbon-enhanced metal-poor (CEMP) stars, confirming the trend of increasing [C/Fe] abundance ratios with decreasing metallicity. Two of these objects can be classified as CEMP-no stars, adding to the growing number of these objects at [Fe/H]<-3. We also find four neutron-capture enhanced stars in the sample, one of which has [Eu/Fe] of 0.8 with clear r-process signatures. These pilot sample stars are the most metal-poor ([Fe/H]<-3.0) of the brightest stars included in CASH and are used to calibrate a newly-developed, automated stellar parameter and abundance determination pipeline. This code will be used for the entire ~500 star CASH snapshot sample. We find that the pipeline results are statistically identical for snapshot spectra when compared to a traditional, manual analysis from a high-resolution spectrum.

Noncommutative geometry may be a starting point to a quantum gravity. We study the influence of the spacetime noncommutative parameter on the strong field gravitational lensing in the noncommutative Schwarzschild black-hole spacetime and obtain the angular position and magnification of the relativistic images. Supposing that the gravitational field of the supermassive central object of the galaxy described by this metric, we estimate the numerical values of the coefficients and observables for strong gravitational lensing. Comparing to the Reissner-Norstr\"{om} black hole, we find that the influences of the spacetime noncommutative parameter is similar to those of the charge, just these influences are much smaller. This may offer a way to distinguish a noncommutative black hole from a Reissner-Norstr\"{om} black hole, and may probe the spacetime noncommutative constant $\vartheta$ \cite{foot} by the astronomical instruments in the future.

We present the high resolution spectra of the D type symbiotic stars V1016 Cygni and HM Sagittae obtained with the Bohyunsan Optical Echelle Spectrograph (BOES), and investigate the double-peaked asymmetric profiles of the Raman scattered O VI 6825. By adopting a wind accretion disk model, we assume that the O VI emission region is described by a Keplerian thin disk. The Raman scattering occurs in a neutral region near the giant, taking in the form of a slow stellar wind, part of which is ionized by the strong UV radiation from the hot white dwarf. Using a Monte Carlo technique, we compute the line profiles that are modulated by the slow spherical stellar wind from the giant component with the ionization front approximated by a hyperboloid. In order to account for the asymmetry and the existence of a central dip in the profiles, we add an O VI resonance scattering region between the hot white dwarf and the giant star which hinders the incidence of slightly blue O VI photons upon the H I region. Overall good fits to the observed data are obtained from our model, which lends support to the accretion disk emission model in these objects. The best fitting parameters for V1016 Cyg are $v_o=30{\rm km s^{-1}}$, $v_\infty=11{\rm km s^{-1}}$, and $v_{c}=10{\rm km s^{-1}}$, where $v_o$, $v_\infty$ and $v_{c}$ are the velocity of the outer disk rim, the terminal velocity of the giant wind, and the velocity component of the resonance scattering O VI region along the binary axis, respectively. Similar fitting parameters $v_o=27{\rm km s^{-1}}$, $v_\infty=10{\rm km s^{-1}}$ and $v_{c}=9{\rm km s^{-1}}$ are obtained for HM Sge. We also investigate the effect of a hot spot in a disk that is well known in accretion disks in cataclysmic variables.

To maximize the number of planet detections, current microlensing follow-up observations are focusing on high-magnification events which have a higher chance of being perturbed by central caustics. In this paper, we investigate the properties of central caustics and the perturbations induced by them. We derive analytic expressions of the location, size, and shape of the central caustic as a function of the star-planet separation, $s$, and the planet/star mass ratio, $q$, under the planetary perturbative approximation and compare the results with those based on numerical computations. While it has been known that the size of the planetary caustic is \propto \sqrt{q}, we find from this work that the dependence of the size of the central caustic on $q$ is linear, i.e., \propto q, implying that the central caustic shrinks much more rapidly with the decrease of $q$ compared to the planetary caustic. The central-caustic size depends also on the star-planet separation. If the size of the caustic is defined as the separation between the two cusps on the star-planet axis (horizontal width), we find that the dependence of the central-caustic size on the separation is \propto (s+1/s). While the size of the central caustic depends both on $s$ and q, its shape defined as the vertical/horizontal width ratio, R_c, is solely dependent on the planetary separation and we derive an analytic relation between R_c and s. Due to the smaller size of the central caustic combined with much more rapid decrease of its size with the decrease of q, the effect of finite source size on the perturbation induced by the central caustic is much more severe than the effect on the perturbation induced by the planetary caustic. Abridged.

Various methods have been proposed to search for extrasolar planets. Compared to the other methods, microlensing has unique applicabilities to the detections of Earth-mass and free-floating planets. However, the microlensing method is seriously flawed by the fact that the masses of the detected planets cannot be uniquely determined. Recently, Gould, Gaudi, & Han introduced an observational setup that enables one to resolve the mass degeneracy of the Earth-mass planets. The setup requires a modest adjustment to the orbit of an already proposed Microlensing planet-finder satellite combined with ground-based observations. In this paper, we show that a similar observational setup can also be used for the mass determinations of free-floating planets with masses ranging from ~0.1 M_J to several Jupiter masses. If the proposed observational setup is realized, the future lensing surveys will play important roles in the studies of Earth-mass and free-floating planets, which are the populations of planets that have not been previously probed.