Search

Fast radio bursts (FRBs) are millisecond-duration radio waves from the Universe. Even though more than 50 physical models have been proposed, the origin and physical mechanism of FRB emissions are still unknown. The classification of FRBs is one of the primary approaches to understanding their mechanisms, but previous studies classified conventionally using only a few observational parameters, such as fluence and duration, which might be incomplete. To overcome this problem, we use an unsupervised machine-learning model, the Uniform Manifold Approximation and Projection (UMAP) to handle seven parameters simultaneously, including amplitude, linear temporal drift, time duration, central frequency, bandwidth, scaled energy, and fluence. We test the method for homogeneous 977 sub-bursts of FRB 20121102A detected by the Arecibo telescope. Our machine-learning analysis identified five distinct clusters, suggesting the possible existence of multiple different physical mechanisms responsible for the observed FRBs from the FRB 20121102A source. The geometry of the emission region and the propagation effect of FRB signals could also make such distinct clusters. This research will be a benchmark for future FRB classifications when dedicated radio telescopes such as the Square Kilometer Array (SKA) or Bustling Universe Radio Survey Telescope in Taiwan (BURSTT) discover more FRBs than before., Comment: Accepted for publication for PASA

Polycyclic aromatic hydrocarbon (PAH) dust emission has been proposed as an effective extinction-independent star formation rate (SFR) indicator in the mid-infrared (MIR), but this may depend on conditions in the interstellar medium. The coverage of the AKARI/Infrared Camera (IRC) allows us to study the effects of metallicity, starburst intensity, and active galactic nuclei on PAH emission in galaxies with $f_{\nu}(L18W)\lesssim 19$ AB mag. Observations include follow-up, rest-frame optical spectra of 443 galaxies within the AKARI North Ecliptic Pole survey that have IRC detections from 7-24 $\mu$m. We use optical emission line diagnostics to infer SFR based on H$\alpha$ and [O II]$\lambda\lambda 3726,3729$ emission line luminosities. The PAH 6.2 $\mu$m and PAH 7.7 $\mu$m luminosities ($L(PAH\ 6.2\ \mu m)$ and $L(PAH\ 7.7\ \mu m)$, respectively) derived using multi-wavelength model fits are consistent with those derived from slitless spectroscopy within 0.2 dex. $L(PAH\ 6.2\ \mu m)$ and $L(PAH\ 7.7\ \mu m)$ correlate linearly with the 24 $\mu$m-dust corrected H$\alpha$ luminosity only for normal, star-forming ``main-sequence" galaxies. Assuming multi-linear correlations, we quantify the additional dependencies on metallicity and starburst intensity, which we use to correct our PAH SFR calibrations at $0

The nature of the interaction between active galactic nuclei (AGNs) and their host galaxies remains an unsolved question. Therefore, conducting an AGN census is valuable to AGN research. Nevertheless, a significant fraction of AGNs are obscured by their environment, which blocks UV and optical emissions due to the dusty torus surrounding the central supermassive black hole (SMBH). To overcome this challenge, mid-infrared (IR) surveys have emerged as a valuable tool for identifying obscured AGNs, as the obscured light is re-emitted in this range. With its high sensitivity, the James Webb Space Telescope (JWST) uncovered more fainter objects than previous telescopes. By applying the SED fitting, this work investigates AGN candidates in JWST Cosmic Evolution Early Release Science (CEERS) fields. We identified 42 candidates, 30 of them are classified as composites ($0.2\leq f_{\rm AGN, IR}< 0.5$), and 12 of them are AGNs ($f_{\rm AGN, IR}\geq 0.5$). We report the AGN luminosity contributions and AGN number fractions as a function of redshift and total infrared luminosity, showing that previously reported increasing relations are not apparent in our sample due to the sample size. We also extend the previous results on ultra-luminous infrared galaxies (ULIRGs, $L_{\rm TIR}\geq 10^{12} L_{\odot}$) to less luminous AGNs, highlighting the power of JWST., Comment: 15 pages, 20 figures, 4 tables. Accepted for publication in MNRAS. The 3 min summary: https://www.youtube.com/watch?v=mWUebbgUOh8

Mid-infrared (MIR) light from galaxies is sensitive to dust-obscured star-formation activities because it traces the characteristic emission of dust heated by young, massive stars. By constructing the MIR luminosity functions (LFs), we are able to quantify the overall dusty star formation history and the evolution of galaxies over cosmic time. In this work, we report the first rest-frame MIR LFs at 7.7, 10, 12.8, 15, 18, and 21 $\mu$m as well as the total IR LF from the James Webb Space Telescope (JWST) Cosmic Evolution Early Release Science (CEERS) survey. We identify 506 galaxies at $z=0-5.1$ in the CEERS survey that also have optical photometry from the Hubble Space Telescope. With the unprecedented sensitivity of the JWST, we probe the faintest end of the LFs at $z=0-1$ down to $L^* \sim 10^7 L_\odot$, $\sim 2$ orders of magnitude fainter than those from the previous generation of IR space telescopes. Our findings connect well with and continue the faint end of the MIR LFs from the deepest observations in past works. As a proxy of star formation history, we present the MIR-based luminosity density up to $z\simeq4.0$, marking the first probe of the early Universe by JWST MIRI., Comment: 22 pages, 22 figures, 7 tables. Accepted for publication in MNRAS. A summary video can be found at https://youtu.be/TRb6bjmGfOU

It has been an unanswered question how many dusty galaxies have been undetected from the state-of-the-art observational surveys. JWST enables us to detect faint IR galaxies that have prominent polycyclic aromatic hydrocarbon (PAH) features in the mid-IR wavelengths. PAH is a valuable tracer of star formation and dust properties in the mid-infrared wavelength. The JWST Cosmic Evolution Early Release Science (CEERS) fields provide us with wavelength coverage from 7.7 to 21 $\mu$m using six photometric bands of the mid-infrared instrument (MIRI). We have identified galaxies dominated by mid-IR emission from PAHs, termed PAH galaxies. From our multi-band photometry catalogue, we selected ten PAH galaxies displaying high flux ratios of $\log(S_{15}/S_{10}) > 0.8$. The SED fitting analysis indicates that these galaxies are star-forming galaxies with total IR luminosities of $10^{10}$ $\sim$ $10^{11.5}$ $L_{\odot}$ at z $\sim 1$. The morphology of PAH galaxies does not show any clear signatures of major merging or interaction within the MIRI resolution. The majority of them are on the star-formation main sequence at $z \sim 1$. Our result demonstrates that JWST can detect PAH emissions from normal star-forming galaxies at $z \sim 1$, in addition to ultra-luminous infrared galaxies (ULIRGs) or luminous infrared galaxies (LIRGs)., Comment: 12 pages, 20 figures, 4 tables. Accepted by MNRAS. A summary video is at https://www.youtube.com/watch?v=UtPaVTFM4f8&ab_channel=NTHUCosmology

With the advent of the James Webb Space Telescope (JWST), extra-galactic source count studies were conducted down to sub-microJy in the mid-infrared (MIR), which is several tens of times fainter than what the previous-generation infrared (IR) telescopes achieved in the MIR. In this work, we aim to interpret the JWST source counts and constrain cosmic star-formation history (CSFH) and black hole accretion history (BHAH). We employ the backward evolution of local luminosity functions (LLFs) of galaxies to reproduce the observed source counts from sub-microJy to a few tens of mJy in the MIR bands of the JWST. The shapes of the LLFs at the MIR bands are determined using the model templates of the spectral energy distributions (SEDs) for five representative galaxy types (star-forming galaxies, starbursts, composite, AGN type 2 and 1). By simultaneously fitting our model to all the source counts in the six MIR bands, along with the previous results, we determine the best-fit evolutions of MIR LFs for each of the five galaxy types, and subsequently estimate the CSFH and BHAH. Thanks to the JWST, our estimates are based on several tens of times fainter MIR sources, the existence of which was merely an extrapolation in previous studies., Comment: 15 pages, 12 figures, published in MNRAS, https://doi.org/10.1093/mnras/stad3499. A summary video is https://youtu.be/Md6wragrYyM

Fast Radio Bursts (FRBs) are mysterious bursts in the millisecond timescale at radio wavelengths. Currently, there is little understanding about the classification of repeating FRBs, based on difference in physics, which is of great importance in understanding their origin. Recent works from the literature focus on using specific parameters to classify FRBs to draw inferences on the possible physical mechanisms or properties of these FRB subtypes. In this study, we use publicly available 1652 repeating FRBs from FRB121102 detected with the Five-hundred-meter Aperture Spherical Telescope (FAST), and studied them with an unsupervised machine learning model. By fine-tuning the hyperparameters of the model, we found that there is an indication for four clusters from the bursts of FRB121102 instead of the two clusters ("Classical" and "Atypical") suggested in the literature. Wherein, the "Atypical" cluster can be further classified into three sub-clusters with distinct characteristics. Our findings show that the clustering result we obtained is more comprehensive not only because our study produced results which are consistent with those in the literature but also because our work uses more physical parameters to create these clusters. Overall, our methods and analyses produced a more holistic approach in clustering the repeating FRBs of FRB121102., Comment: 24 pages, 14 figures, accepted for publication in MNRAS. For summary video, please see https://www.youtube.com/watch?v=wYx6t2G__84&list=PLOpYDs2PkYlYIiKDjDz6r6aKXcXdJZXYb&index=13&ab_channel=NCHUAstronomy

We present a distant T$-$type brown dwarf candidate at $\approx2.55$ kpc discovered in the Cosmic Evolution Early Release Science (CEERS) fields by James Webb Space Telescope (JWST) NIRCam. In addition to the superb sensitivity, we utilised 7 filters from JWST in near-IR and thus is advantageous in finding faint, previously unseen brown dwarfs. From the model spectra in new JWST/NIRCam filter wavelengths, the selection criteria of F115W-F277W$<$-0.8 and F277W-F444W$>$1.1 were chosen to target the spectrum features of brown dwarfs having temperatures from 500K to 1300K. Searching through the data from Early Release Observations (ERO) and Early Release Science (ERS), we find 1 promising candidate in the CEERS field. The result of SED fitting suggested an early T spectral type with a low effective temperature of T$_\text{eff}\approx$1300K, the surface gravity of $\log{g}\approx5.25\text{cm s}^{-2}$, and an eddy diffusion parameter of logK$_{zz}\approx7\text{cm}^2 \text{s}^{-1}$, which indicates an age of $\approx$1.8Gyr and a mass of $\approx0.05$M$_{\odot}$. In contrast to typically found T$-$dwarf within several hundred parsecs, the estimated distance of the source is $\approx2.55$kpc, showing the JWST's power to extend the search to a much larger distance., Comment: 5 pages, 6 figures and 1 table; accepted for publication in MNRAS; A summary video is available at https://youtu.be/PQW79tuS0mI

Understanding dark matter is one of the most urgent questions in modern physics. A very interesting candidate is primordial black holes (PBHs; Carr2016). For the mass ranges of $< 10^{-16} M_{\odot}$ and $> 100 M_{\odot}$, PBHs have been ruled out. However, they are still poorly constrained in the mass ranges of $10^{-16} - 100 M_{\odot}$ (Belotsky et al. 2019). Fast radio bursts (FRBs) are millisecond flashes of radio light of unknown origin mostly from outside the Milky Way. Due to their short timescales, gravitationally lensed FRBs, which are yet to be detected, have been proposed as a useful probe for constraining the presence of PBHs in the mass window of $< 100M_{\odot}$ (Mu\~noz et al. 2016). Up to now, the most successful project in finding FRBs has been CHIME. Due to its large field of view (FoV), CHIME is detecting at least 600 FRBs since 2018. However, none of them is confirmed to be gravitationally lensed (Leung et al. 2022). Taiwan plans to build a new telescope, BURSTT dedicated to detecting FRBs. Its survey area will be 25 times greater than CHIME. BURSTT can localize all of these FRBs through very-long-baseline interferometry (VLBI). We estimate the probability to find gravitationally lensed FRBs, based on the scaled redshift distribution from the latest CHIME catalog and the lensing probability function from Mu\~noz et al. (2016). BURSTT-2048 can detect ~ 24 lensed FRBs out of ~ 1,700 FRBs per annum. With BURSTT's ability to detect nanosecond FRBs, we can constrain PBHs to form a part of dark matter down to $10^{-4}M_{\odot}$., Comment: Accepted for publication in ApJ. A summary video is available at this https://youtu.be/yivrtvuMDHE

Fast radio bursts (FRBs) are astronomical transients with millisecond timescales. Although most of the FRBs are not observed to repeat, a few of them are detected to repeat more than hundreds of times. There exist a large variety of physical properties among these bursts, suggesting heterogeneous mechanisms of FRBs. In this paper, we conduct a categorisation on the extremely frequently repeating FRB 20201124A with the assistance of machine learning, as such techniques have the potential to use subtle differences and correlations that humans are unaware of to better classify bursts. The research is carried out by applying the unsupervised Uniform Manifold Approximation and Projection (UMAP) model on the FRB 20201124A data provided by Five-hundred-meter Aperture Spherical radio Telescope (FAST). The algorithm eventually categorises the bursts into three clusters. In addition to the two categories in previous work based on waiting time, a new way for categorisation has been found. The three clusters are either high energy, high frequency, or low frequency, reflecting the distribution of FRB energy and frequency. Importantly, a similar machine learning result is found in another frequently repeating FRB20121102A, implying a common mechanism among this kind of FRB. This work is one of the first steps towards the systematical categorisation of the extremely frequently repeating FRBs., Comment: 8 pages, 9 figures, accepted for publication in MNRAS. For summary video, please see https://www.youtube.com/watch?v=1PPCMKcrIf4&list=PLOpYDs2PkYlYIiKDjDz6r6aKXcXdJZXYb&index=14&ab_channel=NCHUAstronomy

We report on 167 infrared (IR) galaxies selected by AKARI and IRAS and detected in the Atacama Cosmology Telescope (ACT) Data Release 5 (DR5) sky maps at the 98, 150 and 220 GHz frequency bands. Of these detections, 134 (80%) of the millimeter counterparts are first-time identifications with ACT. We expand the previous ACT extragalactic source catalogs, by including new 98 GHz detections and measurements from ACT DR5. We also report flux density measurements at the 98, 150, and 220 GHz frequency bands. We compute $\alpha_{98-150}$, $\alpha_{98-220}$, and $\alpha_{150-220}$ millimeter-wave spectral indices and far-IR to millimeter-wave spectral indices between 90 micron and 98, 150, and 220 GHz. We specify the galaxy type, based on $\alpha_{150-220}$. We combine publicly available multiwavelength data-including ultraviolet, optical, near-IR, mid-IR, far-IR, and the millimeter measurements obtained in this work-and perform spectral energy distribution (SED) fitting with CIGALE. With the radio emission decomposition advantage of CIGALE V2022.0, we identify the origins of the millimeter emissions for 69 galaxies in our sample. Our analysis also shows that millimeter data alone indicates the need for a radio synchrotron component in the SEDs that are produced by active galactic nuclei (AGNs) and/or star formation. We present SEDs and measured physical properties of these galaxies, such as the dust luminosity, AGN luminosity, the total IR luminosity, and the ratio of the IR and radio luminosity. We quantify the relationships between the total IR luminosity and the millimeter-band luminosities, which can be used in the absence of SED analysis., Comment: 29 pages, 6 tables, 11 figures, 1 appendix with 1 further figure. Accepted for publication in ApJS

Fast radio bursts (FRBs) are intense bursts of radio emission with durations of milliseconds. Although researchers have found them happening frequently all over the sky, they are still in the dark to understand what causes the phenomena because the existing radio observatories have encountered certain challenges during the discovery of FRB progenitors. The construction of Bustling Universe Radio Survey Telescope in Taiwan (BURSTT) is being proposed to solve these challenges. We simulate mock Galactic FRB-like events by applying a range of spatial distributions, pulse widths and luminosity functions. The effect of turbulent Interstellar Medium (ISM) on the detectability of FRB-like events within the Milky Way plane is considered to estimate the dispersion measure and pulse scattering of mock events. We evaluate the fraction of FRB-like events in the Milky Way that are detectable by BURSTT and compare the result with those by Survey for Transient Astronomical Radio Emission 2 (STARE2) and Galactic Radio Explorer (GReX). We find that BURSTT could increase the detection rate by more than two orders of magnitude compared with STARE2 and GReX, depending on the slope of luminosity function of the events. We also investigate the influence of the specifications of BURSTT on its detection improvement. This leads to the fact that greatly higher sensitivity and improved coverage of the Milky Way plane have significant effects on the detection improvement of BURSTT. We find that the upgrade version of BURSTT, BURSTT-2048 could increase the detection rate of faint Galactic FRB-like events by a factor of 3., Comment: 9 pages, 9 figures, 2 tables, accepted for publication in the Monthly Notices of the Royal Astronomical Society

The Hubble constant ($H_{0}$) is a measurement to describe the expansion rate of the Universe in the current era. However, there is a $4.4\sigma$ discrepancy between the measurements from the early Universe and the late Universe. In this research, we propose a model-free and distance-free method to constrain $H_{0}$. Combining Friedman-Lema\^itre-Robertson-Walker cosmology with geometrical relation of the proper motion of extragalactic jets, the lower limit ($H_{\rm 0,min}$) of $H_{0}$ can be determined using only three cosmology-free observables: the redshifts of the host galaxies, as well as the approaching and receding angular velocities of radio jets. Using these, we propose to use the Kolmogorov-Smirnov test (K-S test) between cumulative distribution functions of $H_{\rm 0,min}$ to differentiate cosmology. We simulate 100, 200, and 500 extragalactic jets with 3 levels of accuracy of the proper motion ($\mu_{a}$ and $\mu_{r}$), at $10\%$, $5\%$, and $1\%$, corresponding to the accuracies of the current and future radio interferometers. We perform K-S tests between the simulated samples as theoretical distributions with different $H_{0}$ and power-law index of velocity distribution of jets and mock observational data. Our result suggests increasing sample sizes leads to tighter constraints on both power-law index and the Hubble constant at moderate accuracy (i.e., $10\%$ and $5\%$) while at $1\%$ accuracy, increasing sample sizes leads to tighter constraints on power-law index more. Improving accuracy results in better constraints in the Hubble constant compared with the power-law index in all cases but it alleviates the degeneracy., Comment: 13 pages, 14 figures, accepted for publication in MNRAS

Addressing how strong UV radiation affects galaxy formation is central to understanding their evolution. The quenching of star formation via strong UV radiation (from starbursts or AGN) has been proposed in various scenes to solve certain astrophysical problems. Around luminous sources, some evidence of decreased star formation has been found but is limited to a handful of individual cases. No direct, conclusive evidence on the actual role of strong UV radiation in quenching star formation has been found. Here we present statistical evidence of decreased number density of faint (AB magnitude $\geq$ 24.75 mag) Ly\alpha emitters (LAEs) around bright (AB magnitude < 24.75 mag) LAEs even when the radius goes up to 10 pMpc for $z$ $\simeq$ 5.7 LAEs. A similar trend is found for z $\simeq$ 6.6 LAEs but only within 1 pMpc radius from the bright LAEs. We use a large sample of 1077 (962) LAEs at $z$ $\simeq$ 5.7 ($z$ $\simeq$ 6.6) selected in total areas of 14 (21) deg$^2$ with Subaru/Hyper Suprime-Cam narrow-band data, and thus, the result is of statistical significance for the first time at these high redshift ranges. A simple analytical calculation indicates that the radiation from the central LAE is not enough to suppress LAEs with AB mag $\geq$ 24.75 mag around them, suggesting additional physical mechanisms we are unaware of are at work. Our results clearly show that the environment is at work for the galaxy formation at $z$ $\sim$ 6 in the Universe., Comment: Accepted for publication at MNRAS

Source counts -- the number density of sources as a function of flux density -- represent one of the fundamental metrics in observational cosmology due to their straightforward and simple nature. It is an important tool that provides information on galaxy formation and evolution. Source counting is a direct measurement. Compared to advanced analyzes that require more observational input such as luminosity/mass functions, it is less affected by any cosmological parameter assumptions or any errors propagated from luminosities. In this study, we present source counts at the six mid-infrared bands, i.e., 7.7, 10, 12.8, 15, 18, and 21 $\mu$m from the MIR instrument of the James Webb Space Telescope (JWST). Contrasted with the infrared source counts achieved by prior generations of infrared space telescopes, our source counts delve up to $\sim$100 times deeper, showcasing the exceptional sensitivity of the JWST, and aligning with the model predictions based on preceding observations. In a follow-up study, we utilize our source counts to establish a new IR galaxy population evolutionary model that provides a physical interpretation., Comment: 10 pages, 24 figures, accepted for publication in MNRAS, https://youtu.be/xlX25O8T89g

We present mid-infrared galaxy number counts based on the Early Release Observations obtained by the James Webb Space Telescope (JWST) at 7.7-, 10- and 15-$\mu$m (F770W, F1000W and F1500W, respectively) bands of the Mid-Infrared Instrument (MIRI). Due to the superior sensitivity of JWST, the 80 percent completeness limits reach 0.32, 0.79 and 2.0 $\mu$Jy in F770W, F1000W and F1500W filters, respectively, i.e., $\sim$100 times deeper than previous space infrared telescopes such as Spitzer or AKARI. The number counts reach much deeper than the broad bump around $0.05\sim0.5$ mJy due to polycyclic aromatic hydrocarbon (PAH) emissions. An extrapolation towards fainter flux from the evolutionary models in the literature agrees amazingly well with the new data, where the extrapolated faint-end of infrared luminosity functions combined with the cosmic star-formation history to higher redshifts can reproduce the deeper number counts by JWST. Our understanding of the faint infrared sources has been confirmed by the observed data due to the superb sensitivity of JWST., Comment: 6 pages, 8 figures. Accepted for publication in MNRAS

A correlation between the intrinsic energy and the burst duration of non-repeating fast radio bursts (FRBs) has been reported. If it exists, the correlation can be used to estimate intrinsic energy from the duration, and thus can provide us with a new distance measure for cosmology. However, the correlation suffered from small number statistics (68 FRBs) and was not free from contamination by latent repeating populations, which might not have such a correlation. How to separate/exclude the repeating bursts from the mixture of all different types of FRBs is essential to see this property. Using a much larger sample from the new FRB catalogue (containing 536 FRBs) recently released by the CHIME/FRB project, combined with a new classification method developed based on unsupervised machine learning, we carried out further scrutiny of the relation. We found that there is a weak correlation between the intrinsic energy and duration for non-repeating FRBs at z < 0.3 with Kendall's tau correlation coefficient of 0.239 and significance of 0.001 (statistically significant), whose slope looks similar to that of gamma-ray bursts. This correlation becomes weaker and insignificant at higher redshifts (z > 0.3), possibly due to the lack of the faint FRBs at high-z and/or the redshift evolution of the correlation. The scattering time in the CHIME/FRB catalogue shows an intriguing trend: it varies along the line obtained from linear fit on the energy versus duration plane between these two parameters. A possible cosmological application of the relation must wait for faint FRBs at high-z., Comment: 9 pages, 7 figures, accepted for publication in MNRAS

Fast Radio Bursts (FRBs) are bright millisecond-duration radio transients that appear about 1,000 times per day, all-sky, for a fluence threshold 5 Jy ms at 600 MHz. The FRB radio-emission physics and the compact objects involved in these events are subjects of intense active debate. To better constrain source models, the Bustling Universe Radio Survey Telescope in Taiwan (BURSTT) is optimized to discover and localize a large sample of rare, high-fluence, nearby FRBs. This is the population most amenable to multi-messenger, multi-wavelength follow-up, allowing deeper understanding of source mechanisms. BURSTT will provide horizon-to-horizon sky coverage with a half power field-of-view (FoV) of $\sim$10$^{4}$ deg$^{2}$, a 400 MHz effective bandwidth between 300-800 MHz, and sub-arcsecond localization, made possible using outrigger stations hundreds to thousands of km from the main array. Initially, BURSTT will employ 256 antennas. After tests of various antenna designs and optimization of system performance we plan to expand to 2048 antennas. We estimate that BURSTT-256 will detect and localize $\sim$100 bright ($\geq$100 Jy ms) FRBs per year. Another advantage of BURSTT's large FoV and continuous operation will be greatly enhanced monitoring of FRBs for repetition. The current lack of sensitive all-sky observations likely means that many repeating FRBs are currently cataloged as single-event FRBs., Comment: Published version

A1689-zD1 is one of the most distant galaxies, discovered with the aid of gravitational lensing, providing us with an important opportunity to study galaxy formation in the very early Universe. In this study, we report the detection of [C II]158$\mu$m and [O III]88$\mu$m emission lines of A1689-zD1 in the ALMA Bands 6 and 8. We measure the redshift of this galaxy as $z_{\rm{sys}}=7.133\pm0.005$ based on the [C II] and [O III] emission lines, consistent with that adopted by Bakx et al. (2021). The observed $L_{[\rm{O\,III]}}/L_{[\rm{C\,II]}}$ ratio is $2.09\pm0.09$, higher than most of the local galaxies, but consistent with other $z\sim7$ galaxies. The moderate-spatial resolution of ALMA data provided us with a precious opportunity to investigate spatial variation of $L_{[\rm{O\,III]}}/L_{[\rm{C\,II]}}$. In contrast to the average value of 2.09, we find a much higher $L_{[\rm{O\,III]}}/L_{[\rm{C\,II]}}$ of $\sim 7$ at the center of the galaxy. This spatial variation of $L_{[\rm{O\,III]}}/L_{[\rm{C\,II]}}$ was seldom reported for other high-z galaxies. It is also interesting that the peak of the ratio does not overlap with optical peaks. Possible physical reasons include a central AGN, shock heating from merging, and starburst. Our moderate-spatial resolution data also reveals that in addition to the observed two clumps shown in previous HST images, there is a redshifted segment to the west of the northern optical clump. Such a structure is consistent with previous claims that A1689-zD1 is a merging galaxy, but with the northern redshifted part being some ejected materials, or that the northern redshifted materials being from a third more highly obscured region of the galaxy., Comment: 11 pages, 5 figures, accepted for publication in ApJ. A video summary is available at https://youtu.be/VkR9FlclAp4

Fast radio bursts (FRBs) are mysterious millisecond pulses in radio, most of which originate from distant galaxies. Revealing the origin of FRBs is becoming central in astronomy. The redshift evolution of the FRB energy function, i.e., the number density of FRB sources as a function of energy, provides important implications for the FRB progenitors. Here we show the energy functions of FRBs selected from the recently released Canadian Hydrogen Intensity Mapping Experiment (CHIME) catalogue using the $V_{\rm max}$ method. The $V_{\rm max}$ method allows us to measure the redshift evolution of the energy functions as it is without any prior assumption on the evolution. We use a homogeneous sample of 164 non-repeating FRB sources, which are about one order of magnitude larger than previously investigated samples. The energy functions of non-repeating FRBs show Schechter function-like shapes at $z\lesssim1$. The energy functions and volumetric rates of non-repeating FRBs decrease towards higher redshifts similar to the cosmic stellar-mass density evolution: there is no significant difference between the non-repeating FRB rate and cosmic stellar-mass density evolution with a 1\% significance threshold, whereas the cosmic star-formation rate scenario is rejected with a more than 99\% confidence level. Our results indicate that the event rate of non-repeating FRBs is likely controlled by old populations rather than young populations which are traced by the cosmic star-formation rate density. This suggests old populations such as old neutron stars and black holes as more likely progenitors of non-repeating FRBs., Comment: Accepted for publication in MNRAS. A summary video is available at https://youtu.be/qbUctrFHLAs

Einstein's General Relativity (GR) is the basis of modern astronomy and astrophysics. Testing the validity of basic assumptions of GR is important. In this work, we test a possible violation of the Weak Equivalence Principle (WEP), i.e., there might be a time-lag between photons of different frequencies caused by the effect of gravitational fields if the speeds of photons are slightly different at different frequencies. We use Fast Radio Bursts (FRBs) , which are astronomical transients with millisecond timescales detected in the radio frequency range. Being at cosmological distances, accumulated time delay of FRBs can be caused by the plasma in between an FRB source and an observer, and by gravitational fields in the path of the signal. We segregate the delay due to dispersion and gravitational field using the post-Newtonian formalism (PPN) parameter $\Delta \gamma$, which defines the space-curvature due to gravity by a unit test mass. We did not detect any time-delay from FRBs but obtained tight constraints on the upper limit of $\Delta \gamma$. For FRB20181117C with $z = 1.83 \pm 0.28$ and $\nu_{obs}$ = $676.5\,{\rm MHz}$, the best possible constraint is obtained at log($\Delta \gamma$) = $-21.58 ^{+0.10}_{-0.12}$ and log($\Delta \gamma$/$r_{\rm E}$) = $-21.75 ^{+0.10}_{-0.14}$, respectively, where $r_{\rm E}$ is the energy ratio of two photons of the same FRB signal. This constraint is about one order of magnitude better than the previous constraint obtained with FRBs, and five orders tighter than any constraint obtained using other cosmological sources., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society

Fast radio bursts (FRBs) are astronomical transients with millisecond timescales occurring at cosmological distances. The observed time lag between different energies of each FRB is well described by the inverse-square law of the observed frequency, i.e., dispersion measure. Therefore, FRBs provide one of the ideal laboratories to test Einstein's weak equivalence principle (WEP): the hypothetical time lag between photons with different energies under a gravitational potential. If WEP is violated, such evidence should be exposed within the observational uncertainties of dispersion measures, unless the WEP violation also depends on the inverse-square of the observed frequency. In this work, we constrain the difference of gamma parameters ($\Delta\gamma$) between photons with different energies using the observational uncertainties of FRB dispersion measures, where $\Delta\gamma=0$ for Einstein's general relativity. Adopting the averaged 'Shapiro time delay' for cosmological sources, FRB 121002 at $z=1.6\pm0.3$ and FRB 180817.J1533+42 at $z=1.0\pm0.2$ place the most stringent constraints of $\log\Delta\gamma<-20.8\pm0.1$ and $\log(\Delta\gamma/r_{E}) < -20.9\pm0.2$, respectively, where $r_{E}$ is the energy ratio between the photons. The former is about three orders of magnitude lower than those of other astrophysical sources in previous works under the same formalization of the Shapiro time delay while the latter is comparable to the tightest constraint so far., Comment: Accepted for publication in Physical Review D. A summary video is available at https://youtu.be/Rt7tXN02wpQ

The origins of fast radio bursts (FRBs), astronomical transients with millisecond timescales, remain unknown. One of the difficulties stems from the possibility that observed FRBs could be heterogeneous in origin; as some of them have been observed to repeat, and others have not. Due to limited observing periods and telescope sensitivities, some bursts may be misclassified as non-repeaters. Therefore, it is important to clearly distinguish FRBs into repeaters and non-repeaters, to better understand their origins. In this work, we classify repeaters and non-repeaters using unsupervised machine learning, without relying on expensive monitoring observations. We present a repeating FRB recognition method based on the Uniform Manifold Approximation and Projection (UMAP). The main goals of this work are to: (i) show that the unsupervised UMAP can classify repeating FRB population without any prior knowledge about their repetition, (ii) evaluate the assumption that non-repeating FRBs are contaminated by repeating FRBs, and (iii) recognise the FRB repeater candidates without monitoring observations and release a corresponding catalogue. We apply our method to the Canadian Hydrogen Intensity Mapping Experiment Fast Radio Burst (CHIME/FRB) database. We found that the unsupervised UMAP classification provides a repeating FRB completeness of 95 per cent and identifies 188 FRB repeater source candidates from 474 non-repeater sources. This work paves the way to a new classification of repeaters and non-repeaters based on a single epoch observation of FRBs., Comment: 10 pages, 9 figures, accepted for publication in MNRAS. For summary video, please see https://youtu.be/fWfvfFPDhcQ

We study the galaxy merger fraction and its dependence on star formation mode in the5.4 square degrees of the North Ecliptic Pole-Wide field. We select 6352 galaxies withAKARI 9{\mu}m detections, and identify mergers among them using the Gini coefficientand M20derived from the Subaru/HSC optical images. We obtain the total infraredluminosity and star formation rate of galaxies using the spectral energy distributiontemplates based on one band, AKARI 9{\mu}m. We classify galaxies into three differentstar formation modes (i.e. starbursts, main sequence, and quiescent galaxies) andcalculate the merger fractions for each. We find that the merger fractions of galaxiesincrease with redshift atz<0.6. The merger fractions of starbursts are higher thanthose of main sequence and quiescent galaxies in all redshift bins. We also examinethe merger fractions of far-infrared detected galaxies which have at least one detectionfromHerschel/SPIRE. We find thatHerscheldetected galaxies have higher mergerfraction compared to non-Herscheldetected galaxies, and bothHerscheldetected andnon-Herscheldetected galaxies show clearly different merger fractions depending onthe star formation modes., Comment: 13 pages, 11 figures, Accepted for publication by MNRAS

How does the environment affect active galactic nucleus (AGN) activity? We investigated this question in an extinction-free way, by selecting 1120 infrared galaxies in the $AKARI$ North Ecliptic Pole Wide field at redshift $z$ $\leq$ 1.2. A unique feature of the $AKARI$ satellite is its continuous 9-band infrared (IR) filter coverage, providing us with an unprecedentedly large sample of IR spectral energy distributions (SEDs) of galaxies. By taking advantage of this, for the first time, we explored the AGN activity derived from SED modelling as a function of redshift, luminosity, and environment. We quantified AGN activity in two ways: AGN contribution fraction (ratio of AGN luminosity to the total IR luminosity), and AGN number fraction (ratio of number of AGNs to the total galaxy sample). We found that galaxy environment (normalised local density) does not greatly affect either definitions of AGN activity of our IRG/LIRG samples (log ${\rm L}_{\rm TIR}$ $\leq$ 12). However, we found a different behavior for ULIRGs (log ${\rm L}_{\rm TIR}$ $>$ 12). At our highest redshift bin (0.7 $\lesssim$ z $\lesssim$ 1.2), AGN activity increases with denser environments, but at the intermediate redshift bin (0.3 $\lesssim$ z $\lesssim$ 0.7), the opposite is observed. These results may hint at a different physical mechanism for ULIRGs. The trends are not statistically significant (p $\geq$ 0.060 at the intermediate redshift bin, and p $\geq$ 0.139 at the highest redshift bin). Possible different behavior of ULIRGs is a key direction to explore further with future space missions (e.g., $JWST$, $Euclid$, $SPHEREx$)., Comment: Accepted for publication in MNRAS. A summary video is available at https://www.youtube.com/watch?v=7Y_a0kJkLI4&ab_channel=NthuCosmology

The search for extraterrestrial intelligence (SETI) has been conducted for nearly 60 years. A Dyson Sphere, a spherical structure that surrounds a star and transports its radiative energy outward as an energy source for an advanced civilisation, is one of the main targets of SETI. In this study, we discuss whether building a Dyson Sphere around a black hole is effective. We consider six energy sources: (i) the cosmic microwave background, (ii) the Hawking radiation, (iii) an accretion disk, (iv) Bondi accretion, (v) a corona, and (vi) relativistic jets. To develop future civilisations (for example, a Type II civilisation), $4\times10^{26}\,{\rm W}$($1\,{\rm L_{\odot}}$) is expected to be needed. Among (iii) to (vi), the largest luminosity can be collected from an accretion disk, reaching $10^{5}\,{\rm L_{\odot}}$, enough to maintain a Type II civilisation. Moreover, if a Dyson Sphere collects not only the electromagnetic radiation but also other types of energy (e.g., kinetic energy) from the jets, the total collected energy would be approximately 5 times larger. Considering the emission from a Dyson Sphere, our results show that the Dyson Sphere around a stellar-mass black hole in the Milky Way ($10\,\rm kpc$ away from us) is detectable in the ultraviolet$(\rm 10-400\,{\rm nm)}$, optical$(\rm 400-760\,{\rm nm)}$, near-infrared($\rm 760\,{\rm nm}-5\,{\rm \mu m}$), and mid-infrared($\rm 5-40\,{\rm \mu m}$) wavelengths via the waste heat radiation using current telescopes such as Galaxy Evolution Explorer Ultraviolet Sky Surveys. Performing model fitting to observed spectral energy distributions and measuring the variability of radial velocity may help us to identify these possible artificial structures., Comment: This paper has been accepted for publication in MNRAS

Context. The North Ecliptic Pole (NEP) field provides a unique set of panchromatic data, well suited for active galactic nuclei (AGN) studies. Selection of AGN candidates is often based on mid-infrared (MIR) measurements. Such method, despite its effectiveness, strongly reduces a catalog volume due to the MIR detection condition. Modern machine learning techniques can solve this problem by finding similar selection criteria using only optical and near-infrared (NIR) data. Aims. Aims of this work were to create a reliable AGN candidates catalog from the NEP field using a combination of optical SUBARU/HSC and NIR AKARI/IRC data and, consequently, to develop an efficient alternative for the MIR-based AKARI/IRC selection technique. Methods. A set of supervised machine learning algorithms was tested in order to perform an efficient AGN selection. Best of the models were formed into a majority voting scheme, which used the most popular classification result to produce the final AGN catalog. Additional analysis of catalog properties was performed in form of the spectral energy distribution (SED) fitting via the CIGALE software. Results. The obtained catalog of 465 AGN candidates (out of 33 119 objects) is characterized by 73% purity and 64% completeness. This new classification shows consistency with the MIR-based selection. Moreover, 76% of the obtained catalog can be found only with the new method due to the lack of MIR detection for most of the new AGN candidates. Training data, codes and final catalog are available via the github repository. Final AGN candidates catalog will be also available via the CDS service after publication., Comment: Accepted for publication in Astronomy and Astrophysics

Revealing the cosmic reionisation history is at the frontier of extragalactic astronomy. The power spectrum of the cosmic microwave background (CMB) polarisation can be used to constrain the reionisation history. Here we propose a CMB-independent method using fast radio bursts (FRBs) to directly measure the ionisation fraction of the intergalactic medium (IGM) as a function of redshift. FRBs are new astronomical transients with millisecond timescales. Their dispersion measure (DM$_{\rm IGM}$) is an indicator of the amount of ionised material in the IGM. Since the differential of DM$_{\rm IGM}$ against redshift is proportional to the ionisation fraction, our method allows us to directly measure the reionisation history without any assumption on its functional shape. As a proof of concept, we constructed mock non-repeating FRB sources to be detected with the Square Kilometre Array, assuming three different reionisation histories with the same optical depth of Thomson scattering. We considered three cases of redshift measurements: (A) spectroscopic redshift for all mock data, (B) spectroscopic redshift for 10% of mock data, and (C) redshift estimated from an empirical relation of FRBs between their time-integrated luminosity and rest-frame intrinsic duration. In all cases, the reionisation histories are consistently reconstructed from the mock FRB data using our method. Our results demonstrate the capability of future FRBs in constraining the reionisation history., Comment: Accepted for publication in MNRAS. A summary video is available at https://youtu.be/uis6h_cBnpE

To understand the cosmic accretion history of supermassive black holes, separating the radiation from active galactic nuclei (AGNs) and star-forming galaxies (SFGs) is critical. However, a reliable solution on photometrically recognising AGNs still remains unsolved. In this work, we present a novel AGN recognition method based on Deep Neural Network (Neural Net; NN). The main goals of this work are (i) to test if the AGN recognition problem in the North Ecliptic Pole Wide (NEPW) field could be solved by NN; (ii) to shows that NN exhibits an improvement in the performance compared with the traditional, standard spectral energy distribution (SED) fitting method in our testing samples; and (iii) to publicly release a reliable AGN/SFG catalogue to the astronomical community using the best available NEPW data, and propose a better method that helps future researchers plan an advanced NEPW database. Finally, according to our experimental result, the NN recognition accuracy is around 80.29% - 85.15%, with AGN completeness around 85.42% - 88.53% and SFG completeness around 81.17% - 85.09%., Comment: 12 pages, 12 figures

The $AKARI$ space infrared telescope has performed near- to mid-infrared (MIR) observations on the North Ecliptic Pole Wide (NEPW) field (5.4 deg$^2$) for about one year. $AKARI$ took advantage of its continuous nine photometric bands, compared with NASA's $Spitzer$ and WISE space telescopes, which had only four filters with a wide gap in the MIR. The $AKARI$ NEPW field lacked deep and homogeneous optical data, limiting the use of nearly half of the IR sources for extra-galactic studies owing to the absence of photometric redshifts (photo-zs). To remedy this, we have recently obtained deep optical imaging over the NEPW field with 5 bands ($g$, $r$, $i$, $z$, and $Y$) of the Hyper Suprime-Camera (HSC) on the Subaru 8m telescope. We optically identify AKARI-IR sources along with supplementary $Spitzer$ and WISE data as well as pre-existing optical data. In this work, we derive new photo-zs using a $\chi^2$ template-fitting method code ($Le$ $Phare$) and reliable photometry from 26 selected filters including HSC, $AKARI$, CFHT, Maidanak, KPNO, $Spitzer$ and WISE data. We take 2026 spectroscopic redshifts (spec-z) from all available spectroscopic surveys over the NEPW to calibrate and assess the accuracy of the photo-zs. At z < 1.5, we achieve a weighted photo-z dispersion of $\sigma_{\Delta{z/(1+z)}}$ = 0.053 with $\eta$ = 11.3% catastrophic errors., Comment: 20 pages, 13 figures, accepted for publication in MNRAS. For summary video, please see http://youtu.be/hjNJRCoBIgg

The north ecliptic pole (NEP) field is a natural deep field location for many satellite observations. It has been targeted manytimes since it was surveyed by the AKARI space telescope with its unique wavelength coverage from the near- to mid-infrared(mid-IR). Many follow-up observations have been carried out and made this field one of the most frequently observed areas witha variety of facilities, accumulating abundant panchromatic data from X-ray to radio wavelength range. Recently, a deep opticalsurvey with the Hyper Suprime-Cam (HSC) at the Subaru telescope covered the NEP-Wide (NEPW) field, which enabled us toidentify faint sources in the near- and mid-IR bands, and to improve the photometric redshift (photo-z) estimation. In this work,we present newly identified AKARI sources by the HSC survey, along with multi-band photometry for 91,861 AKARI sourcesobserved over the NEPW field. We release a new band-merged catalogue combining various photometric data from GALEXUV to the submillimetre (sub-mm) bands (e.g., Herschel/SPIRE, JCMT/SCUBA-2). About 20,000 AKARI sources are newlymatched to the HSC data, most of which seem to be faint galaxies in the near- to mid-infrared AKARI bands. This cataloguemotivates a variety of current research, and will be increasingly useful as recently launched (eROSITA/ART-XC) and futurespace missions (such as JWST, Euclid, and SPHEREx) plan to take deep observations in the NEP field., Comment: 17 pages, 12 Figures, accepted in MNRAS, (in press)

In order to understand the interaction between the central black hole and the whole galaxy or their co-evolution history along with cosmic time, a complete census of active galactic nuclei (AGN) is crucial. However, AGNs are often missed in optical, UV and soft X-ray observations since they could be obscured by gas and dust. A mid-infrared (mid-IR) survey supported by multiwavelength data is one of the best ways to find obscured AGN activities because it suffers less from extinction. Previous large IR photometric surveys, e.g., $WISE$ and $Spitzer$, have gaps between the mid-IR filters. Therefore, star forming galaxy (SFG)-AGN diagnostics in the mid-IR were limited. The $AKARI$ satellite has a unique continuous 9-band filter coverage in the near to mid-IR wavelengths. In this work, we take advantage of the state-of-the-art spectral energy distribution (SED) modelling software, CIGALE, to find AGNs in mid-IR. We found 126 AGNs in the NEP-Wide field with this method. We also investigate the energy released from the AGN as a fraction of the total IR luminosity of a galaxy. We found that the AGN contribution is larger at higher redshifts for a given IR luminosity. With the upcoming deep IR surveys, e.g., $JWST$, we expect to find more AGNs with our method., Comment: 14 pages, 18 figures. Accepted for publication in MNRAS. For associated video, please see https://www.youtube.com/watch?v=B01jL0Bol9Q&feature=emb_logo

Fast radio bursts (FRBs) are millisecond transients of unknown origin(s) occurring at cosmological distances. Here we, for the first time, show time-integrated-luminosity functions and volumetric occurrence rates of non-repeating and repeating FRBs against redshift. The time-integrated-luminosity functions of non-repeating FRBs do not show any significant redshift evolution. The volumetric occurrence rates are almost constant during the past $\sim$10 Gyr. The nearly-constant rate is consistent with a flat trend of cosmic stellar-mass density traced by old stellar populations. Our findings indicate that the occurrence rate of non-repeating FRBs follows the stellar-mass evolution of long-living objects with $\sim$Gyr time scales, favouring e.g. white dwarfs, neutron stars, and black holes, as likely progenitors of non-repeating FRBs. In contrast, the occurrence rates of repeating FRBs may increase towards higher redshifts in a similar way to the cosmic star formation-rate density or black hole accretion-rate density if the slope of their luminosity function does not evolve with redshift. Short-living objects with $\lesssim$ Myr time scales associated with young stellar populations (or their remnants, e.g., supernova remnants, young pulsars, and magnetars) or active galactic nuclei might be favoured as progenitor candidates of repeating FRBs., Comment: Accepted for publication in MNRAS. A video summary is available at https://youtu.be/8CpiJepkBeo The catalogue used in this work is available at http://www.phys.nthu.edu.tw/~tetsuya/

We report a promising candidate for a distant supercluster at z $\sim 1.1$ in the Dark Energy Survey Science Verification data. We examine smoothed semi-3D galaxy density maps in various photo-z slices. Among several overdense regions, in this work we report the most significant one as having a $3\sigma$ overdensity at a redshift of $\sim1.1$, over a $\sim160$ Mpc scale, much larger than the regular cluster scale (several Mpc). The shape of the supercluster is not circular in the sky projection. Therefore, we regard the point of maximum overdensity as the center for quantitative measurements. Two different estimates suggest the mass of the supercluster to be $1.37\substack{+1.31 \\ -0.79} \times 10^{17} M_{\odot}$, more than one order more massive than regular galaxy clusters. Except for protosuperclusters found with emission-line galaxies, this could be the most distant supercluster to date defined by regular galaxies. A spectroscopic confirmation would make this a very interesting object for cosmology. We discuss the possible implications of such a massive structure for $\Lambda$CDM cosmology.

Fast radio bursts (FRBs) are mysterious extragalactic radio signals. Revealing their origin is one of the central foci in modern astronomy. Previous studies suggest that occurrence rates of non-repeating and repeating FRBs could be controlled by the cosmic stellar-mass density (CSMD) and star formation-rate density (CSFRD), respectively. The Square Kilometre Array (SKA) is one of the best future instruments to address this subject due to its high sensitivity and high-angular resolution. Here, we predict the number of FRBs to be detected with the SKA. In contrast to previous predictions, we estimate the detections of non-repeating and repeating FRBs separately, based on latest observational constraints on their physical properties including the spectral indices, FRB luminosity functions, and their redshift evolutions. We consider two cases of redshift evolution of FRB luminosity functions following either the CSMD or CSFRD. At $z\gtrsim2$, $z\gtrsim6$ and $z\gtrsim10$, non-repeating FRBs will be detected with the SKA at a rate of $\sim10^{4}$, $\sim10^{2}$, and $\sim10$ (sky$^{-1}$ day$^{-1}$), respectively, if their luminosity function follows the CSMD evolution. At $z\gtrsim1$, $z\gtrsim2$, and $z\gtrsim4$, sources of repeating FRBs will be detected at a rate of $\sim10^{3}$, $\sim10^{2}$, and $\lesssim10$ (sky$^{-1}$ day$^{-1}$), respectively, assuming that the redshift evolution of their luminosity function is scaled with the CSFRD. These numbers could change by about one order of magnitude depending on the assumptions on the CSMD and CSFRD. In all cases, abundant FRBs will be detected by the SKA, which will further constrain the luminosity functions and number density evolutions., Comment: Accepted for publication in MNRAS. A summary is available at https://www.youtube.com/watch?v=c-30UghzCO8

Gamma-Ray Bursts (GRBs) can be a promising tracer of cosmic star-formation rate history (CSFRH). In order to reveal the CSFRH using GRBs, it is important to understand whether they are biased tracers or not. For this purpose, it is crucial to understand properties of GRB host galaxies, in comparison to field galaxies. In this work, we report ALMA far-infrared (FIR) observations of six $z\sim2$ IR-bright GRB host galaxies, which are selected for the brightness in IR. Among them, four host galaxies are detected for the first time in the rest-frame FIR. In addition to the ALMA data, we collected multi-wavelength data from previous studies for the six GRB host galaxies. Spectral energy distribution (SED) fitting analyses were performed with \texttt{CIGALE} to investigate physical properties of the host galaxies, and to test whether active galactic nucleus (AGN) and radio components are required or not. Our results indicate that the best-fit templates of five GRB host galaxies do not require an AGN component, suggesting the absence of AGNs. One GRB host galaxy, 080207, shows a very small AGN contribution. While derived stellar masses of the three host galaxies are mostly consistent with those in previous studies, interestingly the value of star-formation rates (SFRs) of all six GRB hosts are inconsistent with previous studies. Our results indicate the importance of rest-frame FIR observations to correctly estimate SFRs by covering thermal emission from cold dust heated by star formation., Comment: This paper has been accepted for publication in MNRAS

We present the physical properties of AKARI sources without optical counterparts in optical images from the Hyper Suprime-Cam (HSC) on the Subaru telescope. Using the AKARI infrared (IR) source catalog and HSC optical catalog, we select 583 objects that do not have HSC counterparts in the AKARI North Ecliptic Pole (NEP) wide survey field ($\sim 5$ deg$^{2}$). Because the HSC limiting magnitude is deep ($g_{\rm AB}$ $\sim 28.6$), these are good candidates for extremely red star-forming galaxies (SFGs) and/or active galactic nuclei (AGNs), possibly at high redshifts. We compile multi-wavelength data out to 500 $\mu$m and use it for Spectral Energy Distribution (SED) fitting with CIGALE to investigate the physical properties of AKARI galaxies without optical counterparts. We also compare their physical quantities with AKARI mid-IR selected galaxies with HSC counterparts. The estimated redshifts of AKARI objects without HSC counterparts range up to $z\sim 4$, significantly higher than that of AKARI objects with HSC counterparts. We find that: (i) 3.6 $-$ 4.5 $\mu$m color, (ii) AGN luminosity, (iii) stellar mass, (iv) star formation rate, and (v) $V$-band dust attenuation in the interstellar medium of AKARI objects without HSC counterparts are systematically larger than those of AKARI objects with counterparts. These results suggest that our sample includes luminous, heavily dust-obscured SFGs/AGNs at $z\sim 1-4$ that are missed by previous optical surveys, providing very interesting targets for the coming James Webb Space Telescope era., Comment: 17 pages, 13 figures, and 2 tables, accepted for publication in ApJ