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

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

We present the jet kinematics of the flat spectrum radio quasar (FSRQ) 4C +21.35 using time-resolved KaVA very long baseline interferometry array radio maps obtained from September 2014 to July 2016. During two out of three observing campaigns, observations were performed bi-weekly at 22 and 43 GHz quasi-simultaneously. At 22 GHz, we identified three jet components near the core with apparent speeds up to (14.4+/-2.1)c. The timing of the ejection of a new component detected in 2016 is consistent with a gamma-ray flare in November 2014. At 43 GHz, we found four inner jet (<3 mas) components with speeds from (3.5+/-1.4)c to (6.8+/-1.5)c. Jet component speeds tend to be higher with increasing distances from the core. We compared our data with archival Very Long Baseline Array (VLBA) data from the Boston University (BU) 43 GHz and the Monitoring Of Jets in Active galactic nuclei with VLBA Experiments (MOJAVE) 15.4 GHz monitoring programs. Whereas MOJAVE data and our data are in good agreement, jet speeds obtained from the BU Program data in the same time period are about twice as high as the ones we obtain from the KaVA data. The discrepancy at 43 GHz indicates that radio arrays with different angular resolution identify and trace different jet features even when the data are obtained at the same frequency and at the same time. The flux densities of jet components decay exponentially, in agreement with a synchrotron cooling time scale of about 1 year. Using known electron Lorentz factor values (about 9,000), we estimate the magnetic field strength to be around 1-3 micro-Tesla. When adopting a jet viewing angle of 5 degrees, the intrinsic jet speed is of order 0.99c.

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

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

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

We performed a molecular line survey of 82 high-mass protostellar objects in search for inflow signatures associated with high-mass star formation. Using the H$^{13}$CO$^+$ (1$-$0) line as an optically thin tracer, we detected a statistically significant excess of blue asymmetric line profiles in the HCO$^+$ (1$-$0) transition, but nonsignificant excesses in the HCO$^+$ (3$-$2) and H$_2$CO (2$_{12}$$-$1$_{11}$) transitions. The negative blue excess for the HCN (3$-$2) transition suggests that the line profiles are affected by dynamics other than inflow motion. The HCO$^+$ (1$-$0) transition thus seems to be the suitable tracer of inflow motions in high-mass star-forming regions, as previously suggested. We found 27 inflow candidates that have at least one blue asymmetric profile and no red asymmetric profile, and derived the inflow velocities to be 0.23$-$2.00 km s$^{-1}$ for 20 of them using a simple two-layer radiative transfer model. Our sample is divided into two groups in different evolutionary stages. The blue excess of the group in relatively earlier evolutionary stages was estimated to be slightly higher than that of the other in the HCO$^+$ (1$-$0) transition.

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

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

Analysis of all-sky Planck submillimetre observations and the IRAS 100um data has led to the detection of a population of Galactic cold clumps. The clumps can be used to study star formation and dust properties in a wide range of Galactic environments. Our aim is to measure dust spectral energy distribution (SED) variations as a function of the spatial scale and the wavelength. We examine the SEDs at large scales using IRAS, Planck, and Herschel data. At smaller scales, we compare with JCMT/SCUBA-2 850um maps with Herschel data that are filtered using the SCUBA-2 pipeline. Clumps are extracted using the Fellwalker method and their spectra are modelled as modified blackbody functions. According to IRAS and Planck data, most fields have dust colour temperatures T_C ~ 14-18K and opacity spectral index values of beta=1.5-1.9. The clumps/cores identified in SCUBA-2 maps have T~ 13K and similar beta values. There are some indications of the dust emission spectrum becoming flatter at wavelengths longer than 500um. In fits involving Planck data, the significance is limited by the uncertainty of the corrections for CO line contamination. The fits to the SPIRE data give a median beta value slightly above 1.8. In the joint SPIRE and SCUBA-2 850um fits the value decreases to beta ~1.6. Most of the observed T-beta anticorrelation can be explained by noise. The typical submillimetre opacity spectral index beta of cold clumps is found to be ~1.7. This is above the values of diffuse clouds but lower than in some previous studies of dense clumps. There is only tentative evidence of T-beta anticorrelation and beta decreasing at millimetre wavelengths.

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.

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.

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 report the initial results of our high-cadence monitoring program on the radio jet in the active galaxy M87, obtained by the KVN and VERA Array (KaVA) at 22 GHz. This is a pilot study that preceded a larger KaVA-M87 monitoring program, which is currently ongoing. The pilot monitoring was mostly performed every two to three weeks from December 2013 to June 2014, at a recording rate of 1 Gbps, obtaining the data for a total of 10 epochs. We successfully obtained a sequence of good quality radio maps that revealed the rich structure of this jet from <~1 mas to 20 mas, corresponding to physical scales (projected) of ~0.1-2 pc (or ~140-2800 Schwarzschild radii). We detected superluminal motions at these scales, together with a trend of gradual acceleration. The first evidence for such fast motions and acceleration near the jet base were obtained from recent VLBA studies at 43 GHz, and the fact that very similar kinematics are seen at a different frequency and time with a different instrument suggests these properties are fundamental characteristics of this jet. This pilot program demonstrates that KaVA is a powerful VLBI array for studying the detailed structural evolution of the M87 jet and also other relativistic jets.

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.

A remarkable prediction of the Standard Model is that, in the absence of corrections lifting the energy density, the Higgs potential becomes negative at large field values. If the Higgs field samples this part of the potential during inflation, the negative energy density may locally destabilize the spacetime. We use numerical simulations of the Einstein equations to study the evolution of inflation-induced Higgs fluctuations as they grow towards the true (negative-energy) minimum. These simulations show that forming a single patch of true vacuum in our past light cone during inflation is incompatible with the existence of our Universe; the boundary of the true vacuum region grows outward in a causally disconnected manner from the crunching interior, which forms a black hole. We also find that these black hole horizons may be arbitrarily elongated---even forming black strings---in violation of the hoop conjecture. By extending the numerical solution of the Fokker-Planck equation to the exponentially suppressed tails of the field distribution at large field values, we derive a rigorous correlation between a future measurement of the tensor-to-scalar ratio and the scale at which the Higgs potential must receive stabilizing corrections in order for the Universe to have survived inflation until today.

We observed 146 Galactic clumps in HCN (4-3) and CS (7-6) with the Atacama Submillimeter Telescope Experiment (ASTE) 10-m telescope. A tight linear relationship between star formation rate and gas mass traced by dust continuum emission was found for both Galactic clumps and the high redshift (z>1) star forming galaxies (SFGs), indicating a constant gas depletion time of ~100 Myr for molecular gas in both Galactic clumps and high z SFGs. However, low z galaxies do not follow this relation and seem to have a longer global gas depletion time. The correlations between total infrared luminosities (L_TIR) and molecular line luminosities (L'_mol) of HCN (4-3) and CS (7-6) are tight and sublinear extending down to clumps with LTIR~10^{3} L_sun. These correlations become linear when extended to external galaxies. A bimodal behavior in the LTIR--L'mol correlations was found for clumps with different dust temperature, luminosity-to-mass ratio, and sigmaline/sigmavir. Such bimodal behavior may be due to evolutionary effects. The slopes of LTIR--L'mol correlations become more shallow as clumps evolve. We compared our results with lower J transition lines in wu et al. (2010). The correlations between clump masses and line luminosities are close to linear for low effective excitation density tracers but become sublinear for high effective excitation density tracers for clumps with LTIR larger than LTIR~10^4.5 Lsun. High effective excitation density tracers cannot linearly trace the total clump masses, leading to a sublinear correlations for both Mclump-L'mol and LTIR-L'mol relations.

The Chandrasekhar-Fermi method is a powerful technique for estimating the strength of the mean magnetic field projected on the plane of the sky. In this paper, we present a technique for improving the Chandrasekhar-Fermi method, in which we take into account the averaging effect arising from independent eddies along the line of sight . In the conventional Chandrasekhar-Fermi method, the strength of fluctuating magnetic field divided by $\sqrt{4 \pi \bar{\rho}}$, where $\bar{\rho}$ is average density, is assumed to be comparable to the line-of-sight velocity dispersion. This however is not true when the driving scale of turbulence $L_f$, i.e. the outer scale of turbulence, is smaller than the size of the system along the line of sight $L_{los}$. In fact, the conventional Chandrasekhar-Fermi method over-estimates the strength of the mean plane-of-the-sky magnetic field by a factor of $\sim \sqrt{ L_{los}/L_f}$. We show that the standard deviation of centroid velocities divided by the average line-of-sight velocity dispersion is a good measure of $\sqrt{ L_{los}/L_f}$, which enables us to propose a modified Chandrasekhar-Fermi method.

The highest brightness temperature ever observed are from "nanoshots" from the Crab pulsar which we argue could be the signature of bursts of vacuum $e^{\pm}$ pair production. If so this would be the first time the astronomical Schwinger effect has been observed. These "Schwinger sparks" would be an intermittent but extremely powerful, $\sim 10^3 L_{\astrosun}$, 10 PeV $e^{\pm}$ accelerator in the heart of the Crab. These nanosecond duration sparks are generated in a volume less than $1 m^3$ and the existence of such sparks has implications for the small scale structure of the magnetic field of young pulsars such as the Crab. This mechanism may also play a role in producing other enigmatic bright short radio transients such as fast radio bursts.

Blue CDM-photon isocurvature perturbations are attractive in terms of observability and may be typical from the perspective of generic mass relations in supergravity. We present and apply three theorems useful for blue isocurvature perturbations arising from linear spectator scalar fields. In the process, we give a more precise formula for the blue spectrum associated with the axion model of 0904.3800, which can in a parametric corner give a factor of O(10) correction. We explain how a conserved current associated with Peccei-Quinn symmetry plays a crucial role and explicitly plot several example spectra including the breaks in the spectra. We also resolve a little puzzle arising from a naive multiplication of isocurvature expression that sheds light on the gravitational imprint of the adiabatic perturbations on the fields responsible for blue isocurvature fluctuations.

Turbulence is ubiquitous in astrophysical fluids such as the interstellar medium (ISM) and the intracluster medium (ICM). In turbulence studies, it is customary to assume that fluid is driven on a single scale. However, in astrophysical fluids, there can be many different driving mechanisms that act on different scales. If there are multiple energy-injection scales, the process of energy cascade and turbulence dynamo will be different compared with the case of single energy-injection scale. In this work, we perform three-dimensional incompressible/compressible magnetohydrodynamic (MHD) turbulence simulations. We drive turbulence in Fourier space in two wavenumber ranges, 2\$\leq k \leq \sqrt{12}\$ (large-scale) and 15 \$\lesssim k \lesssim \$ 26 (small-scale). We inject different amount of energy in each range by changing the amplitude of forcing in the range. We present the time evolution of the kinetic and magnetic energy densities and discuss the turbulence dynamo in the presence of energy injections at two scales. We show how kinetic, magnetic and density spectra are affected by the two-scale energy injections and we discuss the observational implications. In the case \$\epsilon_L < \epsilon_S\$, where \$\epsilon_L\$ and \$\epsilon_S\$ are energy-injection rates at the large and small scales, respectively, our results show that even a tiny amount of large-scale energy injection can significantly change the properties of turbulence. On the other hand, when \$\epsilon_L \gtrsim \epsilon_S\$, the small-scale driving does not influence the turbulence statistics much unless \$\epsilon_L \sim \epsilon_S\$.

Isocurvature perturbations in the inflationary literature typically involve quantum fluctuations of bosonic field degrees of freedom. In this work, we consider isocurvature perturbations from fermionic quantum fluctuations during inflation. When a stable massive fermion is coupled to a non-conformal sector different from the scalar metric perturbations, observably large amplitude scale invariant isocurvature perturbations can be generated. In addition to the computation of the isocurvature two-point function, an estimate of the local non-Gaussianities is also given and found to be promising for observations in a corner of the parameter space. The results provide a new class of cosmological probes for theories with stable massive fermions. On the technical side, we explicitly renormalize the composite operator in curved spacetime and show that gravitational Ward identities play an important role in suppressing certain contributions to the fermionic isocurvature perturbations.

Sources of isocurvature perturbations and large non-Gaussianities include field degrees of freedom whose vacuum expectation values are smaller than the expansion rate of inflation. The inhomogeneities in the energy density of such fields are quadratic in the fields to leading order in the inhomogeneity expansion. Although it is often assumed that such isocurvature perturbations and inflaton-driven curvature perturbations are uncorre- lated, this is not obvious from a direct computational point of view due to the form of the minimal gravitational interactions. We thus compute the irreducible gravitational contributions to the quadratic isocurvature-curvature cross-correlation. We find a small but non-decaying cross-correlation, which in principle serves as a consistency prediction of this large class of isocurvature perturbations. We apply our cross-correlation result to two dark matter isocurvature perturbation scenarios: QCD axions and WIMPZILLAs. On the technical side, we utilize a gravita- tional Ward identity in a novel manner to demonstrate the gauge invariance of the computation. Furthermore, the detailed computation is interpreted in terms of a soft-{\zeta} theorem and a gravitational Ward identity. Finally, we also identify explicitly all the counterterms that are necessary for renormalizing the isocurvature perturbation composite operator in inflationary cosmological backgrounds.

Axions and WIMPZILLAs are well-motivated dark matter candidates with interesting cos- mological implications, such as isocurvature perturbations and non-Gaussianities. How- ever, these predicted implications in the literature are based on the assumption that the cross-correlation between curvature and CDM isocurvature is negligible. This contribution discusses the cross-correlation in the axion and the WIMPZILLA scenarios. Particularly, it is shown that the gravitational Ward identity associated with diffeomorphism invariance plays an important role in this cross-correlation calculation confirming the assumption.

Gravitational particle production naturally occurs during the transition from the inflationary phase to the non-inflationary phase. If the particles are stable and very weakly interacting, they are natural nonthermal dark matter candidates. We show that such nonthermal dark matter particles can produce local non-Gaussianities large enough to be observed by ongoing and near future experiments without being in conflict with the existing isocurvature bounds. Of particular interest is the fact that these particles can be observable through local non-Gaussianities even when they form a very small fraction of the total dark matter content.

Turbulence dynamo deals with amplification of a seed magnetic field in a turbulent medium and has been studied mostly for uniform or spatially homogeneous seed magnetic fields. However, some astrophysical processes (e.g. jets from active galaxies, galactic winds, or ram-pressure stripping in galaxy clusters) can provide localized seed magnetic fields. In this paper, we numerically study amplification of localized seed magnetic fields in a turbulent medium. Throughout the paper, we assume that driving scale of turbulence is comparable to the size of the system. Our findings are as follows. First, turbulence can amplify a localized seed magnetic field very efficiently. The growth rate of magnetic energy density is as high as that for a uniform seed magnetic field. This result implies that a magnetic field ejected from an astrophysical object can be a viable source of magnetic field in a cluster. Second, the localized seed magnetic field disperses and fills the whole system very fast. If turbulence in a system (e.g. a galaxy cluster or a filament) is driven at large scales, we expect that it takes a few large-eddy turnover times for magnetic field to fill the whole system. Third, growth and turbulence diffusion of a localized seed magnetic field are also fast in high magnetic Prandtl number turbulence. Fourth, even in decaying turbulence, a localized seed magnetic field can ultimately fill the whole system. Although the dispersal rate of magnetic field is not fast in purely decaying turbulence, it can be enhanced by an additional forcing.

We revisit the gravitational production of massive Dirac fermions in inflationary cosmology with a focus on clarifying the analytic computation of the particle number density in both the large and the small mass regimes. For the case in which the masses of the gravitationally produced fermions are small compared to the Hubble expansion rate at the end of inflation, we obtain a universal result for the number density that is nearly independent of the details of the inflationary model. The result is identical to the case of conformally coupled scalars up to an overall multiplicative factor of order unity for reasons other than just counting the fermionic degrees of freedom.