We present sensitive 850$\mu$m imaging of the COSMOS field using 640hr of new and archival observations taken with SCUBA-2 at the East Asian Observatory's James Clerk Maxwell Telescope. The SCUBA-2 COSMOS survey (S2COSMOS) achieves a median noise level of $\sigma_{850\mu{\mathrm{m}}}$=1.2mJy/beam over an area of 1.6 sq. degree (MAIN; HST/ACS footprint), and $\sigma_{850\mu{\mathrm{m}}}$=1.7mJy/beam over an additional 1 sq. degree of supplementary (SUPP) coverage. We present a catalogue of 1020 and 127 sources detected at a significance level of >4$\sigma$ and >4.3$\sigma$ in the MAIN and SUPP regions, respectively, corresponding to a uniform 2% false-detection rate. We construct the single-dish 850$\mu$m number counts at $S_{850}$>2mJy and show that these S2COSMOS counts are in agreement with previous single-dish surveys, demonstrating that degree-scale fields are sufficient to overcome the effects of cosmic variance in the $S_{850}$=2-10mJy population. To investigate the properties of the galaxies identified by S2COSMOS sources we measure the surface density of near-infrared-selected galaxies around their positions and identify an average excess of 2.0$\pm$0.2 galaxies within a 13$''$ radius (~100kpc at $z$~2). The bulk of these galaxies represent near-infrared-selected SMGs and/or spatially-correlated sources and lie at a median photometric redshift of $z$=2.0$\pm$0.1. Finally, we perform a stacking analysis at sub-millimeter and far-infrared wavelengths of stellar-mass-selected galaxies ($M_{\star}$=10$^{10}$-10$^{12}{\rm M_{\odot}}$) from $z$=0-4, obtaining high-significance detections at 850um in all subsets (SNR=4-30), and investigate the relation between far-infrared luminosity, stellar mass, and the peak wavelength of the dust SED. The publication of this survey adds a new deep, uniform sub-millimeter layer to the wavelength coverage of this well-studied COSMOS field.

We report a search for astronomical neutrinos in the energy region from several GeV to TeV in the direction of the blazar TXS0506+056 using the Super-Kamiokande detector following the detection of a 100 TeV neutrino from the same location by the IceCube collaboration. Using Super-Kamiokande neutrino data across several data samples observed from April 1996 to February 2018 we have searched for both a total excess above known backgrounds across the entire period as well as localized excesses on smaller time scales in that interval. No significant excess nor significant variation in the observed event rate are found in the blazar direction. Upper limits are placed on the electron and muon neutrino fluxes at 90\% confidence level as $6.03 \times 10^{-7}$ and $4.52 \times 10^{-7}$ to $9.26 \times 10^{-10}$ [${\rm erg}/{\rm cm}^2/{\rm s}$], respectively.

We identify multi-wavelength counterparts to 1,147 submillimeter sources from the S2COSMOS SCUBA-2 survey of the COSMOS field by employing a recently developed radio$+$machine-learning method trained on a large sample of ALMA-identified submillimeter galaxies (SMGs), including 260 SMGs identified in the AS2COSMOS pilot survey. In total, we identify 1,222 optical/near-infrared(NIR)/radio counterparts to the 897 S2COSMOS submillimeter sources with S$_{850}$>1.6mJy, yielding an overall identification rate of ($78\pm9$)%. We find that ($22\pm5$)% of S2COSMOS sources have multiple identified counterparts. We estimate that roughly 27% of these multiple counterparts within the same SCUBA-2 error circles very likely arise from physically associated galaxies rather than line-of-sight projections by chance. The photometric redshift of our radio$+$machine-learning identified SMGs ranges from z=0.2 to 5.7 and peaks at $z=2.3\pm0.1$. The AGN fraction of our sample is ($19\pm4$)%, which is consistent with that of ALMA SMGs in the literature. Comparing with radio/NIR-detected field galaxy population in the COSMOS field, our radio+machine-learning identified counterparts of SMGs have the highest star-formation rates and stellar masses. These characteristics suggest that our identified counterparts of S2COSMOS sources are a representative sample of SMGs at z<3. We employ our machine-learning technique to the whole COSMOS field and identified 6,877 potential SMGs, most of which are expected to have submillimeter emission fainter than the confusion limit of our S2COSMOS surveys (S$_{850}$<1.5mJy). We study the clustering properties of SMGs based on this statistically large sample, finding that they reside in high-mass dark matter halos ($(1.2\pm0.3)\times10^{13}\,h^{-1}\,\rm M_{\odot}$), which suggests that SMGs may be the progenitors of massive ellipticals we see in the local Universe.

Supernova detection is a major objective of the Super-Kamiokande (SK) experiment. In the next stage of SK (SK-Gd), gadolinium (Gd) sulfate will be added to the detector, which will improve the ability of the detector to identify neutrons. A core-collapse supernova will be preceded by an increasing flux of neutrinos and anti-neutrinos, from thermal and weak nuclear processes in the star, over a timescale of hours; some of which may be detected at SK-Gd. This could provide an early warning of an imminent core-collapse supernova, hours earlier than the detection of the neutrinos from core collapse. Electron anti-neutrino detection will rely on inverse beta decay events below the usual analysis energy threshold of SK, so Gd loading is vital to reduce backgrounds while maximising detection efficiency. Assuming normal neutrino mass ordering, more than 200 events could be detected in the final 12 hours before core collapse for a 15-25 solar mass star at around 200 pc, which is representative of the nearest red supergiant to Earth, $\mathrm{\alpha}$Ori (Betelgeuse). At a statistical false alarm rate of 1 per century, detection could be up to 10 hours before core collapse, and a pre-supernova star could be detected by SK-Gd up to 600 pc away. A pre-supernova alert could be provided to the astrophysics community following gadolinium loading.

We present CO J=4-3 line and 3 mm dust continuum observations of a 100 kpc-scale filamentary Ly{\alpha} nebula (SSA22 LAB18) at z=3.1 using the Atacama Large Millimeter/submillimeter Array (ALMA). We detected the CO J=4-3 line at a systemic z(CO)=3.093 {\pm} 0.001 at 11 {\sigma} from one of the ALMA continuum sources associated with the Ly{\alpha} filament. We estimated the CO J=4-3 luminosity of L'CO(4-3)=(2.3\pm0.2)x10^9 K km s^{-1} pc^2 for this CO source, which is one order of magnitude smaller than those of typical z>1 dusty star-forming galaxies (DSFGs) of similar far-infrared luminosity L(IR)~10^{12} Lsun. We derived a molecular gas mass of Mgas=(4.4^{+0.9}_{-0.6})x10^9 Msun and a star-formation rate of SFR=270\pm160 Msun yr^{-1}. We also estimated a gas depletion time of {\tau}(dep)=17\pm10 Myr, being shorter than those of typical DSFGs. It is suggested that this source is in a transition phase from DSFG to a gas-poor, early-type galaxy. From ALMA to Herschel multi-band dust continuum observations, we measured a dust emissivity index {\beta}=2.3\pm0.2, which is similar to those of local gas-poor, early-type galaxies. Such a high {\beta} can be reproduced by specific chemical compositions for interstellar dust at the submillimeter wavelengths from recent laboratory experiments. ALMA CO and multi-band dust continuum observations can constrain the evolutionary stage of high-redshift galaxies through {\tau}(dep) and {\beta}, and thus we can investigate dust chemical compositions even in the early Universe.

To search for dust-obscured star-formation activity in the early Universe, it is essential to obtain a deep and wide submillimeter/millimeter map. The advent of the Atacama Large Millimeter/submillimeter Array (ALMA) has enabled us to obtain such maps at sufficiently high spatial resolution to be free from source confusion. We present a new 1.1mm map obtained by ALMA in the SSA22 field. SSA22 contains a remarkable proto-cluster at z=3.09 and is therefore an ideal region to investigate the role of large-scale cosmic web on dust-obscured star formation. The typical 1sigma depth of our map is 73 uJy/beam at a 0.5 resolution; combined with earlier, archived observations, we map an area of 20 arcmin^2 (71 comoving Mpc^2 at z=3.09). Within the combined survey area we have detected 35 sources at a signal-to-noise ratio (S/N) >5, with flux densities, S1.1mm=0.43--5.6 mJy, equivalent to star-formation rates of >=100--1000 Msun/yr at z=3.09, for a Chabrier initial mass function; of these, 17 are new detections. The cumulative number counts show a factor 3--5x excess compared to blank fields. The excess suggests enhanced dust-enshrouded star-formation activity in the proto-cluster on a 10 comoving Mpc scale, indicating accelerated galaxy evolution in this overdense region.

We report the results of a neutrino search in Super-Kamiokande for coincident signals with the first detected gravitational wave produced by a binary neutron star merger, GW170817, which was followed by a short gamma-ray burst, GRB170817A, and a kilonova/macronova. We searched for coincident neutrino events in the range from 3.5 MeV to $\sim$100 PeV, in a time window $\pm$500 seconds around the gravitational wave detection time, as well as during a 14-day period after the detection. No significant neutrino signal was observed for either time window. We calculated 90% confidence level upper limits on the neutrino fluence for GW170817. From the upward-going-muon events in the energy region above 1.6 GeV, the neutrino fluence limit is $16.0^{+0.7}_{-0.6}$ ($21.3^{+1.1}_{-0.8}$) cm$^{-2}$ for muon neutrinos (muon antineutrinos), with an error range of $\pm5^{\circ}$ around the zenith angle of NGC4993, and the energy spectrum is under the assumption of an index of $-2$. The fluence limit for neutrino energies less than 100 MeV, for which the emission mechanism would be different than for higher-energy neutrinos, is also calculated. It is $6.6 \times 10^7$ cm$^{-2}$ for anti-electron neutrinos under the assumption of a Fermi-Dirac spectrum with average energy of 20 MeV.

We report the Australia Telescope Compact Array 16 cm observation of CIZA J1358.9-4750. Recent X-ray studies imply that this galaxy cluster is composed of merging, binary clusters. Using the EW367 configuration, we found no significant diffuse radio emission in and around the cluster. An upper limit of the total radio power at 1.4 GHz is $\sim 1.1 \times 10^{22}$ Watt/Hz in $30$ square arcminutes which is a typical size of radio relics. It is known that an empirical relation holds between the total radio power and X-ray luminosity of the host cluster. The upper limit is about one order of magnitude lower than the power expected from the relation. Very young ($\sim 70$ Myr) shocks with low Mach numbers ($\sim 1.3$), which are often seen at an early stage of merger simulations, are suggested by the previous X-ray observation. The shocks may generate cosmic-ray electrons with a steep energy spectrum, which is consistent with non-detection of bright ($>10^{23}$ Watt/Hz) relic in this 16 cm band observation. Based on the assumption of energy equipartition, the upper limit gives the magnetic-field strength below $0.68f(D_{\rm los}/{\rm 1~Mpc})^{-1}(\gamma_{\rm min}/200)^{-1}$ $\mu$G, where $f$ is the cosmic-ray total energy density over the cosmic-ray electron energy density, $D_{\rm los}$ is the depth of the shock wave along the sightline and $\gamma_{\rm min}$ is the lower cutoff Lorentz factor of the cosmic-ray electron energy spectrum.

We present the results of a search in the Super-Kamiokande (SK) detector for excesses of neutrinos with energies above a few GeV that are in the direction of the track events reported in IceCube. Data from all SK phases (SK-I through SK-IV) were used, spanning a period from April 1996 to April 2016 and corresponding to an exposure of 225 kilotonne-years . We considered the 14 IceCube track events from a data set with 1347 livetime days taken from 2010 to 2014. We use Poisson counting to determine if there is an excess of neutrinos detected in SK in a 10 degree search cone (5 degrees for the highest energy data set) around the reconstructed direction of the IceCube event. No significant excess was found in any of the search directions we examined. We also looked for coincidences with a recently reported IceCube multiplet event. No events were detected within a $\pm$ 500 s time window around the first detected event, and no significant excess was seen from that direction over the lifetime of SK.

The present paper investigates the temperature structure of the X-ray emitting plasma in the core of the Perseus cluster using the 1.8--20.0 keV data obtained with the Soft X-ray Spectrometer (SXS) onboard the Hitomi Observatory. A series of four observations were carried out, with a total effective exposure time of 338 ks and covering a central region $\sim7'$ in diameter. The SXS was operated with an energy resolution of $\sim$5 eV (full width at half maximum) at 5.9 keV. Not only fine structures of K-shell lines in He-like ions but also transitions from higher principal quantum numbers are clearly resolved from Si through Fe. This enables us to perform temperature diagnostics using the line ratios of Si, S, Ar, Ca, and Fe, and to provide the first direct measurement of the excitation temperature and ionization temperature in the Perseus cluster. The observed spectrum is roughly reproduced by a single temperature thermal plasma model in collisional ionization equilibrium, but detailed line ratio diagnostics reveal slight deviations from this approximation. In particular, the data exhibit an apparent trend of increasing ionization temperature with increasing atomic mass, as well as small differences between the ionization and excitation temperatures for Fe, the only element for which both temperatures can be measured. The best-fit two-temperature models suggest a combination of 3 and 5 keV gas, which is consistent with the idea that the observed small deviations from a single temperature approximation are due to the effects of projection of the known radial temperature gradient in the cluster core along the line of sight. Comparison with the Chandra/ACIS and the XMM-Newton/RGS results on the other hand suggests that additional lower-temperature components are present in the ICM but not detectable by Hitomi SXS given its 1.8--20 keV energy band.

Pulsar magnetospheres have strong magnetic fields and large amounts of plasma. The structures of these magnetospheres are studied using force-free electrodynamics. To understand pulsar magnetospheres, discussions must include their outer region. However, force-free electrodynamics is limited in it does not handle dissipation. Therefore, a resistive pulsar magnetic field model is needed. To break the ideal magnetohydrodynamic (MHD) condition $E \cdot B = 0$, Ohm's law is used. In this work, I introduce resistivity depending upon the distance from the star and obtain a self-consistent steady state by time integration. Poloidal current circuits form in the magnetosphere while the toroidal magnetic-field region expands beyond the light cylinder and the Poynting flux radiation appears. High electric resistivity causes a large space scale poloidal current circuit and the magnetosphere radiates a larger Poynting flux than the linear increase outside of the light cylinder radius. The formed poloidal current circuit has width, which grows with the electric conductivity. This result contributes to a more concrete dissipative pulsar magnetosphere model.

We study the heating mechanisms and Ly{\alpha} escape fractions of 35 Ly{\alpha} blobs (LABs) at z = 3.1 in the SSA22 field. Dust continuum sources have been identified in 11 of the 35 LABs, all with star formation rates (SFRs) above 100 Msun/yr. Likely radio counterparts are detected in 9 out of 29 investigated LABs. The detection of submm dust emission is more linked to the physical size of the Ly{\alpha} emission than to the Ly{\alpha} luminosities of the LABs. A radio excess in the submm/radio detected LABs is common, hinting at the presence of active galactic nuclei. Most radio sources without X-ray counterparts are located at the centers of the LABs. However, all X-ray counterparts avoid the central regions. This may be explained by absorption due to exceptionally large column densities along the line-of-sight or by LAB morphologies, which are highly orientation dependent. The median Ly{\alpha} escape fraction is about 3\% among the submm-detected LABs, which is lower than a lower limit of 11\% for the submm-undetected LABs. We suspect that the large difference is due to the high dust attenuation supported by the large SFRs, the dense large-scale environment as well as large uncertainties in the extinction corrections required to apply when interpreting optical data.

In this publication we continue the work started in Quintero Noda et al. (2017) examining this time a numerical simulation of a magnetic flux tube concentration. Our goal is to study if the physical phenomena that take place in it, in particular, the magnetic pumping, leaves a specific imprint on the examined spectral lines. We find that the profiles from the interior of the flux tube are periodically dopplershifted following an oscillation pattern that is also reflected in the amplitude of the circular polarization signals. In addition, we analyse the properties of the Stokes profiles at the edges of the flux tube discovering the presence of linear polarization signals for the Ca II lines, although they are weak with an amplitude around 0.5% of the continuum intensity. Finally, we compute the response functions to perturbations in the longitudinal field and we estimate the field strength using the weak field approximation. Our results indicate that the height of formation of the spectral lines changes during the magnetic pumping process which makes the interpretation of the inferred magnetic field strength and its evolution more difficult. These results complement those from previous works demonstrating the capabilities and limitations of the 850 nm spectrum for chromospheric Zeeman polarimetry in a very dynamic and complex atmosphere.

Solar "magnetic tornadoes" are produced by rotating magnetic field structures that extend from the upper convection zone and the photosphere to the corona of the Sun. Recent studies show that such rotating features are an integral part of atmospheric dynamics and occur on a large range of spatial scales. A systematic statistical study of magnetic tornadoes is a necessary next step towards understanding their formation and their role for the mass and energy transport in the solar atmosphere. For this purpose, we have developed a new automatic detection method for chromospheric swirls, i.e. the observable signature of solar tornadoes or, more generally, chromospheric vortex flows and rotating motions. Unlike the previous studies that relied on visual inspections, our new method combines a line integral convolution (LIC) imaging technique and a scalar quantity which represents a vortex flow on a two-dimensional plane. We have tested two detection algorithms, based on the enhanced vorticity and vorticity strength quantities, by applying them to 3D numerical simulations of the solar atmosphere with CO5BOLD. We conclude that the vorticity strength method is superior compared to the enhanced vorticity method in all aspects. Applying the method to a numerical simulation of the solar atmosphere revealed very abundant small-scale, short-lived chromospheric vortex flows that had not been found by visual inspection before.

We report the result from observations conducted with the Atacama Large Millimeter/submillimeter Array (ALMA) to detect [CII] 158 um fine structure line emission from galaxies embedded in one of the most spectacular Lyman-alpha blobs (LABs) at z=3.1, SSA22-LAB1. Of three dusty star-forming galaxies previously discovered by ALMA 860 um dust continuum survey toward SSA22-LAB1, we detected the [CII] line from one, LAB1-ALMA3 at z=3.0993+/-0.0004. No line emission was detected, associated with the other ALMA continuum sources or from three rest-frame UV/optical selected z_spec~3.1 galaxies within the field of view. For LAB1-ALMA3, we find relatively bright [CII] emission compared to the infrared luminosity (L_[CII]/L_[CII]) and an extremely high [CII] 158 um and [NII] 205 um emission line ratio (L_[CII]/L_[NII]>55). The relatively strong [CII] emission may be caused by abundant photodissociation regions and sub-solar metallicity, or by shock heating. The origin of the unusually strong [CII] emission could be causally related to the location within the giant LAB, although the relationship between extended Lyman-alpha emission and ISM conditions of associated galaxies is yet to be understand.

We present results from a deep 2'x3' (comoving scale of 3.7 Mpc x 5.5 Mpc at z=3) survey at 1.1 mm taken with the Atacama Large Millimeter/submillimeter Array (ALMA) in the SSA22 field. We observe the core region of a z = 3.09 protocluster, achieving a typical rms sensitivity of 60 micro-Jy/beam at a spatial resolution of 0".7. We detect 18 robust ALMA sources at a signal-to-noise ratio (SNR) > 5. Comparison between the ALMA map and a 1.1 mm map taken with the AzTEC camera on the Atacama Submillimeter Telescope Experiment (ASTE) indicates that three submillimeter sources discovered by the AzTEC/ASTE survey are resolved into eight individual submillimeter galaxies (SMGs) by ALMA. At least ten of our 18 ALMA SMGs have spectroscopic redshifts of z = 3.09, placing them in the protocluster. This shows that a number of dusty starburst galaxies are forming simultaneously in the core of the protocluster. The nine brightest ALMA SMGs with SNR > 10 have a median intrinsic angular size of 0".32+0".13-0".06 (2.4+1.0-0.4 physical kpc at z = 3.09), which is consistent with previous size measurements of SMGs in other fields. As expected the source counts show a possible excess compared to the counts in the general fields at S_1.1mm >= 1.0 mJy due to the protocluster. Our contiguous mm mapping highlights the importance of large-scale structures on the formation of dusty starburst galaxies.

Using radiation magnetohydrodynamic simulations of the solar atmospheric layers from the upper convection zone to the lower corona, we investigate the self-consistent excitation of slow magneto-acoustic body waves (slow modes) in a magnetic flux concentration. We find that the convective downdrafts in the close surroundings of a two-dimensional flux slab "pump" the plasma inside it in the downward direction. This action produces a downflow inside the flux slab, which encompasses ever higher layers, causing an upwardly propagating rarefaction wave. The slow mode, excited by the adiabatic compression of the downflow near the optical surface, travels along the magnetic field in the upward direction at the tube speed. It develops into a shock wave at chromospheric heights, where it dissipates, lifts the transition region, and produces an offspring in the form of a compressive wave that propagates further into the corona. In the wake of downflows and propagating shock waves, the atmosphere inside the flux slab in the chromosphere and higher tends to oscillate with a period of $\nu\approx 4$~mHz. We conclude that this process of "magnetic pumping" is a most plausible mechanism for the direct generation of longitudinal chromospheric and coronal compressive waves within magnetic flux concentrations, and it may provide an important heat source in the chromosphere. It may also be responsible for certain types of dynamic fibrils.

We present a Herschel/SPIRE survey of three protoclusters at z=2-3 (2QZCluster, HS1700, SSA22). Based on the SPIRE colours (S350/S250 and S500/S350) of 250 $\mu$m sources, we selected high redshift dusty star-forming galaxies potentially associated with the protoclusters. In the 2QZCluster field, we found a 4-sigma overdensity of six SPIRE sources around 4.5' (~2.2 Mpc) from a density peak of H$\alpha$ emitters at z=2.2. In the HS1700 field, we found a 5-sigma overdensity of eight SPIRE sources around 2.1' (~1.0 Mpc) from a density peak of LBGs at z=2.3. We did not find any significant overdensities in SSA22 field, but we found three 500 $\mu$m sources are concentrated 3' (~1.4 Mpc) east to the LAEs overdensity. If all the SPIRE sources in these three overdensities are associated with protoclusters, the inferred star-formation rate densities are 10$^3$-10$^4$ times higher than the average value at the same redshifts. This suggests that dusty star-formation activity could be very strongly enhanced in z~2-3 protoclusters. Further observations are needed to confirm the redshifts of the SPIRE sources and to investigate what processes enhance the dusty star-formation activity in z~2-3 protoclusters.

We observed three $\gamma$-ray bursts related to thunderclouds in winter using the prototype of anti-neutrino detector PANDA made of 360-kg plastic scintillator deployed at Ohi Power Station at the coastal area of the Japan Sea. The maximum rate of the events which deposited the energy higher than $3\,$MeV was $(5.5 \pm 0.1) \times 10^2 {\rm /s}$. Monte Carlo simulation showed that electrons with approximately monochromatic energy falling downwards from altitudes of order $100\,$m roughly produced the observed total energy spectra of the bursts. It is supposed that secondary cosmic-ray electrons, which act as seed, were accelerated in electric field of thunderclouds and multiplied by relativistic runaway electron avalanche. We actually found that the $\gamma$-rays of the bursts entered into the detector from the direction close to the zenith. The direction stayed constant during the burst within the detector resolution. In addition, taking advantage of the delayed coincidence detection of the detector, we found neutron events in one of the bursts at the maximum rate of $\sim 14\pm5\,{\rm /s}$.

Abell 548W, one of the galaxy clusters located in the Abell 548 region, has about an order of magnitude lower X-ray luminosity compared to ordinal clusters in view of the well known intracluster medium (ICM) temperature vs X-ray luminosity (kT-L_X) relation. The cluster hosts a pair of diffuse radio sources to the north west and north, both about 10' apart from the cluster center. They are candidate radio relics, frequently associated with merging clusters. A Suzaku deep observation with exposure of 84.4 ks was performed to search signatures for merging in this cluster. The XIS detectors successfully detected the ICM emission out to 16' from the cluster center. The temperature is ~3.6 keV around its center, and ~2 keV at the outermost regions. The hot region (~6 keV) aside the relic candidates shifted to the cluster center reported by XMM-Newton was not seen in the Suzaku data, although its temperature of 3.6 keV itself is higher than the average temperature of 2.5 keV around the radio sources. In addition, a signature of a cool (kT ~0.9 keV) component was found around the north west source. A marginal temperature jump at its outer-edge was also found, consistent with the canonical idea of shock acceleration origin of the radio relics. The cluster has among the highest central entropy of ~400 keV cm^2 and is one of the so-called low surface brightness clusters. Taking into account the fact that its shape itself is relatively circular and smooth and also its temperature structure is nearly flat, possible scenarios for merging is discussed.

The Atacama Large Millimeter/submillimeter Array (ALMA) is a new powerful tool for observing the Sun at high spatial, temporal, and spectral resolution. These capabilities can address a broad range of fundamental scientific questions in solar physics. The radiation observed by ALMA originates mostly from the chromosphere - a complex and dynamic region between the photosphere and corona, which plays a crucial role in the transport of energy and matter and, ultimately, the heating of the outer layers of the solar atmosphere. Based on first solar test observations, strategies for regular solar campaigns are currently being developed. State-of-the-art numerical simulations of the solar atmosphere and modeling of instrumental effects can help constrain and optimize future observing modes for ALMA. Here we present a short technical description of ALMA and an overview of past efforts and future possibilities for solar observations at submillimeter and millimeter wavelengths. In addition, selected numerical simulations and observations at other wavelengths demonstrate ALMA's scientific potential for studying the Sun for a large range of science cases.

We report the results of $1^{\prime}.5 \times3^{\prime}$ mapping at 1.1~mm with the Atacama Large Millimeter/submillimeter Array (ALMA) toward the central region of the $z=3.09$ SSA22 protocluster. By combining our source catalog with archival spectroscopic redshifts, we find that eight submillimeter galaxies (SMGs) with flux densities, $S_{\rm 1.1~mm}=0.7-6.4$~mJy ($L_{\rm IR}\sim10^{12.1}-10^{13.1}L_\odot$) are at $z=3.08-3.10$. Not only are these SMGs members of the protocluster but they in fact reside within the node at the junction of the 50 Mpc-scale filamentary three-dimensional structure traced by Lyman-$\alpha$ emitters (LAEs) in this field. The eight SMGs account for a star formation rate density (SFRD) $\sim$10 $M_\odot$ yr$^{-1}$ Mpc$^{-3}$ in the node, which is two orders of magnitudes higher than the global SFRD at this redshift. We find that four of the eight SMGs host a X-ray luminous active galactic nuclei (AGN). Our results suggest that the vigorous star formation activity and the growth of super massive black holes (SMBHs) occurred simultaneously in the densest regions at $z\sim3$, which may correspond to the most active historical phase of the massive galaxy population found in the core of the clusters in the present universe. Two SMGs are associated with Lyman-$\alpha$ blobs (LABs), implying that the two populations coexist in high density environments for a few cases.

CIZA J1358.9-4750 is a nearby (z = 0.074) pair of clusters of galaxies located close to the Galactic plane. It consists of two X-Ray extended humps at north-west and south-east separated by 14 arcmin (~ 1.2 Mpc), and an X-Ray bright bridge-like structure in between. With Suzaku, the south-east hump was shown to have a temperature of 5.6 keV and the north-west one 4.6 keV. Neither humps exhibit significant central cool component. The bridge region has a temperature higher than 9 keV at the maximum, and this hot region is distributed almost orthogonal to the bridge axis in agreement with the shock heating seen in numerical simulations at an early phase of a head-on major merger. This resemblance is supported by good positional coincidence between the X-Ray peaks and cD galaxies associated with each cluster. In a short exposure XMM-Newton image, a significant intensity jump was found at a position where the Suzaku-measured temperature exhibits a steep gradient. These properties indicate the presence of a shock discontinuity. The Mach number is estimated to be 1.32 from the temperature difference across the identified shock front, which gives the colliding velocity of approximately 1800 km/s. From optical redshifts of the member galaxies, the two clusters are indicated to be merging nearly on the sky plane. Thus, CIZA J1358.9-4750 is considered as a valuable nearby example of early-phase merger with a clear shock feature.

The Solar Simulations for the Atacama Large Millimeter Observatory Network (SSALMON) was initiated in 2014 in connection with two ALMA development studies. The Atacama Large Millimeter/submillimeter Array (ALMA) is a powerful new tool, which can also observe the Sun at high spatial, temporal, and spectral resolution. The international SSALMONetwork aims at coordinating the further development of solar observing modes for ALMA and at promoting scientific opportunities for solar physics with particular focus on numerical simulations, which can provide important constraints for the observing modes and can aid the interpretation of future observations. The radiation detected by ALMA originates mostly in the solar chromosphere - a complex and dynamic layer between the photosphere and corona, which plays an important role in the transport of energy and matter and the heating of the outer layers of the solar atmosphere. Potential targets include active regions, prominences, quiet Sun regions, flares. Here, we give a brief overview over the network and potential science cases for future solar observations with ALMA.

We present a model of the magnetosphere around an oscillating neutron star. The electromagnetic fields are numerically solved by modeling electric charge and current induced by the stellar torsional mode, with particular emphasis on outgoing radiation passing through the magnetosphere. The current is modeled using Ohm's law, whereby an increase in conductivity results in an increase in the induced current. As a result, the fields are drastically modified, and energy flux is thereby enhanced. This behavior is however localized in the vicinity of the surface since the induced current disappears outwardly in our model, in which the exterior is assumed to gradually approach a vacuum.

We investigate mass losses via stellar winds from sun-like main sequence stars with a wide range of activity levels. We perform forward-type magnetohydrodynamical numerical experiments for Alfven wave-driven stellar winds with a wide range of the input Poynting flux from the photosphere. Increasing the magnetic field strength and the turbulent velocity at the stellar photosphere from the current solar level, the mass loss rate rapidly increases at first owing to the suppression of the reflection of the Alfven waves. The surface materials are lifted up by the magnetic pressure associated with the Alfven waves, and the cool dense chromosphere is intermittently extended to 10 -- 20 % of the stellar radius. The dense atmospheres enhance the radiative losses and eventually most of the input Poynting energy from the stellar surface escapes by the radiation. As a result, there is no more sufficient energy remained for the kinetic energy of the wind; the stellar wind saturates in very active stars, as observed in Wood et al. The saturation level is positively correlated with B_{r,0}f_0, where B_{r,0} and f_0 are the magnetic field strength and the filling factor of open flux tubes at the photosphere. If B_{r,0}f_0 is relatively large >~ 5 G, the mass loss rate could be as high as 1000 times. If such a strong mass loss lasts for ~ 1 billion years, the stellar mass itself is affected, which could be a solution to the faint young sun paradox. We derive a Reimers-type scaling relation that estimates the mass loss rate from the energetics consideration of our simulations. Finally, we derive the evolution of the mass loss rates, \dot{M} t^{-1.23}, of our simulations, combining with an observed time evolution of X-ray flux from sun-like stars, which is shallower than \dot{M} t^{-2.33+/-0.55} in Wood et al.(2005).

We propose the experiments on the collision of laser light and high intensity electromagnetic pulses generated by relativistic flying mirrors, with electron bunches produced by a conventional accelerator and with laser wake field accelerated electrons for studying extreme field limits in the nonlinear interaction of electromagnetic waves. The regimes of dominant radiation reaction, which completely changes the electromagnetic wave-matter interaction, will be revealed in the laser plasma experiments. This will result in a new powerful source of ultra short high brightness gamma-ray pulses. A possibility of the demonstration of the electron-positron pair creation in vacuum in a multi-photon processes can be realized. This will allow modeling under terrestrial laboratory conditions neutron star magnetospheres, cosmological gamma ray bursts and the Leptonic Era of the Universe.

We report major performance test results of the Atacama Compact Array (ACA) 7-m and 12-m antennas of ALMA (Atacama Large Millimeter/submillimeter Array). The four major performances of the ACA antennas are all-sky pointing (to be not more than 2.0 arcsec), offset pointing (to be < 0.6 arcsec) surface accuracy (< 25(20) micrometer for 12(7)m-antenna), stability of path-length (15 micrometer over 3 min), and high servo capability (6 degrees/s for Azimuth and 3 degrees/s for Elevation). The high performance of the ACA antenna has been extensively evaluated at the Site Erection Facility area at an altitude of about 2900 meters. Test results of pointing performance, surface performance, and fast motion capability are demonstrated.

Relativistic effects in the oscillatory damping of magnetic disturbances near two-dimensional X-points are investigated. By taking into account displacement current, we study new features of extremely magnetized systems, in which the Alfv\'en velocity is almost the speed of light. The frequencies of the least-damped mode are calculated using linearized relativistic MHD equations for wide ranges of the Lundquist number S and the magnetization parameter $\sigma$. These timescales approach constant values in the large resistive limit: the oscillation time becomes a few times the light crossing time, irrespective of $\sigma$, and the decay time is proportional to $\sigma$ and therefore is longer for a highly magnetized system.

From radiation magnetohydrodynamic simulations of the solar atmosphere we find a new mechanism for the excitation of longitudinal slow modes within magnetic flux concentrations. We find that the convective downdrafts in the immediate surroundings of magnetic elements are responsible for the excitation of slow modes. The coupling between the external downdraft and the plasma motion internal to the flux concentration is mediated by the inertial forces of the downdraft that act on the magnetic flux concentration. These forces, in conjunction with the downward movement, pump the internal atmosphere in the downward direction, which entails a fast downdraft in the photospheric and chromospheric layers of the magnetic element. Subsequent to the transient pumping phase, the atmosphere rebounds, causing a slow mode traveling along the magnetic flux concentration in the upward direction. It develops into a shock wave in chromospheric heights, possibly capable of producing some kind of dynamic fibril. We propose an observational detection of this process.

We examine the recoil velocity induced by the superposition of the magnetic dipole and quadrupole radiation from a pulsar/magnetar born with rapid rotation. The resultant velocity depends on not the magnitude, but rather the ratio of the two moments and their geometrical configuration. The model does not necessarily lead to high spatial velocity for a magnetar with a strong magnetic field, which is consistent with the recent observational upper bound. The maximum velocity predicted with this model is slightly smaller than that of observed fast-moving pulsars.

We determine the spin of a supermassive black hole in the context of discseismology by comparing newly detected quasi-periodic oscillations (QPOs) of radio emission in the Galactic centre, Sagittarius A* (Sgr A*), as well as infrared and X-ray emissions with those of the Galactic black holes. We find that the spin parameters of black holes in Sgr A* and in Galactic X-ray sources have a unique value of $\approx 0.44$ which is smaller than the generally accepted value for supermassive black holes, suggesting evidence for the angular momentum extraction of black holes during the growth of supermassive black holes. Our results demonstrate that the spin parameter approaches the equilibrium value where spin-up via accretion is balanced by spin-down via the Blandford-Znajek mechanism regardless of its initial spin. We anticipate that measuring the spin of black holes by using QPOs will open a new window for exploring the evolution of black holes in the Universe.

We report the detection of radio QPOs with structure changes using the Very Long Baseline Array (VLBA) at 43 GHz. We found conspicuous patterned changes of the structure with P = 16.8 +- 1.4, 22.2 +- 1.4, 31.2 +- 1.5, 56.4 +- 6 min, very roughly in a 3:4:6:10 ratio. The first two periods show a rotating one-arm structure, while the P = 31.4 min shows a rotating 3-arm structure, as if viewed edge-on. At the central 50 microasec the P = 56.4 min period shows a double amplitude variation of those in its surroundings. Spatial distributions of the oscillation periods suggest that the disk of SgrA* is roughly edge-on, rotating around an axis with PA = -10 degree. Presumably, the observed VLBI images of SgrA* at 43 GHz retain several features of the black hole accretion disk of SgrA* in spite of being obscured and broadened by scattering of surrounding plasma.

By solving radiative transfer equations, we examine three-dimensional radiative properties of a magnetohydrodynamic accretion flow model confronting with the observed spectrum of Sgr A*, in the vicinity of supermassive black hole at the Galactic centre. As a result, we find that the core of radio emission is larger than the size of the event horizon shadow and its peak location is shifted from the gravitational centre. We also find that the self-absorbed synchrotron emissions by the superposition of thermal electrons within a few tens of the Schwartzschild radius can account for low-frequency spectra below the critical frequency $\nu_{c}\approx 10^{12}$ Hz. Above the critical frequency, the synchrotron self-Compton emission by thermal electrons can account for variable emissions in recent near-infrared observations. In contrast to the previous study by Ohsuga et al. (2005), we found that the X-ray spectra by Bremsstrahlung emission of thermal electrons for the different mass accretion rates can be consistent with both the flaring state and the quiescent state of Sgr A* observed by {\it Chandra}.

Black-hole accretion systems are known to possess several distinct modes (or spectral states), such as low/hard state, high/soft state, and so on. Since the dynamics of the corresponding flows is distinct, theoretical models were separately discussed for each state. We here propose a unified model based on our new, global, two-dimensional radiation-magnetohydrodynamic simulations. By controlling a density normalization we could for the first time reproduce three distinct modes of accretion flow and outflow with one numerical code. When the density is large (model A), a geometrically thick, very luminous disk forms, in which photon trapping takes place. When the density is moderate (model B), the accreting gas can effectively cool by emitting radiation, thus generating a thin disk, i.e., the soft-state disk. When the density is too low for radiative cooling to be important (model C), a disk becomes hot, thick, and faint; i.e., the hard-state disk. The magnetic energy is amplified within the disk up to about twice, 30%, and 20% of the gas energy in models A, B, and C, respectively. Notably, the disk outflows with helical magnetic fields, which are driven either by radiation pressure force or magnetic pressure force, are ubiquitous in any accretion modes. Finally, our simulations are consistent with the phenomenological alpha-viscosity prescription, that is, the disk viscosity is proportional to the pressure.

We calculate the radiation spectrum and its time variability of the black hole accretion disk-corona system based on the three-dimensional magnetohydrodynamic simulation. In explaining the spectral properties of active galactic nuclei (AGNs), it is often assumed that they consist of a geometrically thin, optically thick disk and hot, optically thin corona surrounding the thin disk. As for a model of the corona, we adopt the simulation data of three-dimensional, non-radiative MHD accretion flows calculated by Kato and coworkers, while for a thin disk we assume a standard type disk. We perform Monte Carlo radiative transfer simulations in the corona, taking into account the Compton scattering of soft photons from the thin disk by hot thermal electrons and coronal irradiation heating of the thin disk, which emits blackbody radiation. By adjusting the density parameter of the MHD coronal flow, we can produce the emergent spectra which are consistent with those of typical Seyfert galaxies. Moreover, we find rapid time variability in X-ray emission spectra, originating from the density fluctuation produced by the magnetorotational instability in the MHD corona. The features of reflection component including iron fluorescent line emission are also briefly discussed.

In spite of a large number of global three-dimensional (3D) magneto-hydrodynamical (MHD) simulations of accretion flows and jets being made recently, their astrophysical relevance for realistic situations is not well known. In order to examine to what extent the simulated MHD flows can account for the observed spectral energy distribution (SED) of Sagittarius A* (Sgr A*), for the first time we calculate the emergent spectra from 3D MHD flows in a wide range of wavelengths (from radio to X-ray) by solving the 3D radiative transfer equations. We use the simulation data by Kato, Mineshige, and Shibata (2004) and perform Monte Carlo radiative transfer simulations, in which synchrotron emission/absorption, free-free emission/absorption, and Compton/inverse Compton scattering are taken into account. We assume two temperature plasmas and calculate electron temperatures by solving the electron energy equation. Only thermal electrons are considered. It is found that the 3D MHD flow generally over-produces X-rays by means of bremsstrahlung radiation from the regions at large radii. A flatter density profile, r^{-a} with a<1, than that of the advection-dominated accretion flow (ADAF), r^{-3/2}, is the main reason for this. If we restrict the size of the emission region to be as small as 10r_s, where r_s is the Schwarzschild radius, the MHD model can reproduce the basic features of the observed SED of Sgr A* during its flaring state. Yet, the spectrum in the quiescent state remains to be understood. We also calculate the time-dependent spectral changes, finding that the fluxes fluctuate in a wide range of the frequency and the flux at each wavelength does not always vary coherently.

We study the radial and vertical oscillations in the three-dimensional magnetohydrodynamic (MHD) accretion flows around black holes. General relativistic effects are taken into account by using pseudo-Newtonian potential. We find that the structure of MHD flows is changed at $3.8 r_{\rm s}\leq r\leq 6.3 r_{\rm s}$ and the two pairs of quasi-periodic oscillations (twin QPOs) are excited in that region. The time evolution of the power spectrum density (PSD) indicates that these twin QPOs are most likely to be produced by the resonance between the Keplerian frequency $\nu_{\rm K}$ and the epicyclic frequency $\nu_{\kappa}$. The PSD shows that the lower peak frequency $\nu_{\rm l}$ corresponds to $\nu_{\rm K}$, while the upper peak frequency $\nu_{\rm u}$ corresponds to $\nu_{\rm K} + \nu_{\kappa}$. The ratio of two peak frequencies is close to 2:3. The results provide the first direct evidence for the excitation of the resonant disk oscillation in the MHD accretion flows.

We present the results of 2.5-dimensional resistive magnetohydrodynamic (MHD) simulations of the magnetic interaction between a weakly magnetized neutron star and its accretion disk. General relativistic effects are simulated by using the pseudo-Newtonian potential. We find that well-collimated jets traveling along the rotation axis of the disk are formed by the following mechanism: (1) The magnetic loops connecting the neutron star and the disk are twisted due to the differential rotation between the neutron star and the disk. (2) Twist injection from the disk initiates expansion of the loop. (3) The expanding magnetic loops create a magnetic tower in which accelerated disk material travel as collimated bipolar jets. The propagation speed of the working surface of the jet is the order of 10% of the speed of light ($\sim 0.1c$). (4) Magnetic reconnection taking place inside the expanding magnetic loops injects hot bubbles intermittently into the magnetic tower. The ejection speed of the bubble is the order of the local Alfv\'{e}n speed of the launching point and $\sim 0.2c$ in our simulations. (5) The hot bubbles moving inside the tower catch up with the working surface of the jet. High energy electrons created by the magnetic reconnection are a plausible source of radio emission. Our model can explain the formation process of a narrow jet from a weakly magnetized ($|{\boldmath$B_{*}$}|\le 10^{9}$ gauss) neutron star and the correlation between radio flares of the core and of the lobe observed in Sco X-1.

We present three-dimensional (3-D) magnetohydrodynamical (MHD) simulations of radiatively inefficient accretion flow around black holes. General relativistic effects are simulated by using the pseudo-Newtonian potential. We start calculations with a rotating torus threaded by localized poloidal magnetic fields with plasma beta, a ratio of the gas pressure to the magnetic pressure, $\beta =10$ and 100. When the bulk of torus material reaches the innermost region close to a central black hole, a magnetically driven jet emerges. This magnetic jet is derived by vertically inflating toroidal fields (`magnetic tower') and has a two-component structure: low-$\beta$ ($\lsim 1$) plasmas threaded with poloidal (vertical) fields are surrounded by that with toroidal fields. The collimation width of the jet depends on external pressure, pressure of ambient medium; the weaker the external pressure is, the wider and the longer-lasting becomes the jet. Unless the external pressure is negligible, the bipolar jet phase ceases after several dynamical timescales at the original torus position and a subsequent quasi-steady state starts. The black hole is surrounded by quasi-spherical zone with highly inhomogeneous structure in which toroidal fields are dominant except near the rotation axis. Mass accretion takes place mainly along the equatorial plane. Comparisons with other MHD simulation results and observational implications are discussed.

One hundred and forty six long-period red variable stars in the Large Magellanic Cloud (LMC) from the three year MOA project database were analysed. A careful periodic analysis was performed on these stars and a catalogue of their magnitudes, colours, periods and amplitudes is presented. We convert our blue and red magnitudes to $K$ band values using 19 oxygen-rich stars. A group of red short-period stars separated from the Mira sequence has been found on a (log P, K) diagram. They are located at the short period side of the Mira sequence consistent with the work of Wood and Sebo (1996). There are two interpretations for such stars; a difference in pulsation mode or a difference in chemical composition. We investigated the properties of these stars together with their colour, amplitude and periodicity. We conclude that they have small amplitudes and less regular variability. They are likely to be higher mode pulsators. A large scatter has been also found on the long period side of the (log P, K) diagram. This is possibly a systematic spread given that the blue band of our photometric system covers both standard B and V bands and affects carbon-rich stars.

We carried out two dimensional high-resolution magnetohydrodynamic (MHD) simulations of an accretion disk around a weakly magnetized neutron star. General relativistic effects are taken into account by using pseudo-Newtonian potential of Pacz\'{y}nski, B., P. J. Witta (1980). When magnetic loops connect the neutron star and the accretion disk, the twist injection from the disk or from the rotating neutron star triggers expansion of the loops. Since the expanding magnetic loops prevent inflow toward the magnetic poles of the neutron star, disk matter accumulates on the boundary between the magnetosphere and the disk. Magnetic reconnection taking place in the loops creates a channel along which the disk matter can accrete and unloads the magnetosphere. This process produces quasi-periodic variation of the accretion flow in the innermost region of the disk. We found two kinds of oscillations. One is the magnetospheric oscillation regulated by magnetic reconnection. The other is the radial disk oscillation. The typical frequency of the oscillations is 100 Hz to 2 kHz. Furthermore, we predict that QPO sources inevitably accompany X-ray flares by magnetic reconnection and bipolar outflows of hot X-ray emitting plasma similar to the optical jets in protostars.