Showing 1–50 of 354 results for author: Kim, M

Eddington Ratios of Dust-obscured Quasars at $z \lesssim 1$: Evidence Supporting Dust-obscured Quasars as Young Quasars

Authors: Dohyeong Kim, Yongjung Kim, Myungshin Im, Eilat Glikman, Minjin Kim, Tanya Urrutia, Gu Lim

Abstract: Dust-obscured quasars have been suspected as the intermediate stage galaxies between merger-driven star-forming galaxies and unobscured quasars. This merger-driven galaxy evolution scenario suggests that dust-obscured quasars exhibit higher Eddington ratios ($λ_{\rm Edd}$) than those of unobscured quasars. However, their high dust obscuration poses challenges to accurately measuring their… ▽ More Dust-obscured quasars have been suspected as the intermediate stage galaxies between merger-driven star-forming galaxies and unobscured quasars. This merger-driven galaxy evolution scenario suggests that dust-obscured quasars exhibit higher Eddington ratios ($λ_{\rm Edd}$) than those of unobscured quasars. However, their high dust obscuration poses challenges to accurately measuring their $λ_{\rm Edd}$ using commonly employed bolometric luminosity ($L_{\rm bol}$) and black hole (BH) mass ($M_{\rm BH}$) estimators based on the ultraviolet (UV) or optical luminosity. Recently, Kim et al. (2023) established new estimators for $L_{\rm bol}$ and $M_{\rm BH}$ based on mid-infrared (MIR) continuum luminosity ($L_{\rm MIR}$), which are less affected by dust obscuration. These estimators enable the study of a large number of dust-obscured quasars across a wide redshift range. In this study, we measure the $λ_{\rm Edd}$ values of 30 dust-obscured quasars at $z \lesssim 1$, the largest sample size to date, using the $L_{\rm MIR}$-based $L_{\rm bol}$ and $M_{\rm BH}$ estimators. Our findings reveal that dust-obscured quasars exhibit significantly higher $λ_{\rm Edd}$ values compared to unobscured quasars. Moreover, we confirm that the enhanced $λ_{\rm Edd}$ values of dust-obscured quasars maintain consistency across the redshift span of 0 to 1. Our results strongly support the picture that dust-obscured quasars are in the earlier stage than unobscured quasars in the merger-driven galaxy evolutionary track. △ Less

Submitted 7 August, 2024; v1 submitted 6 August, 2024; originally announced August 2024.

Comments: 15 pages, 13 figures, Accepted for publication in A&A

The characterisation of water ice in debris discs: implications for JWST scattered light observations

Authors: Minjae Kim, Grant M. Kennedy, Veronica Roccatagliata

Abstract: Water ice plays a crucial role throughout the different stages of planetary evolution and is abundant in the Universe. However, its presence and nature in debris discs of exoplanetary systems are not yet strongly established observationally. In this study, we quantify and discuss the impact of ice parameters such as volume fraction ${\mathcal{F}}_{\rm ice}$, blow-out grain size, size distribution,… ▽ More Water ice plays a crucial role throughout the different stages of planetary evolution and is abundant in the Universe. However, its presence and nature in debris discs of exoplanetary systems are not yet strongly established observationally. In this study, we quantify and discuss the impact of ice parameters such as volume fraction ${\mathcal{F}}_{\rm ice}$, blow-out grain size, size distribution, and its phase on the observational appearance of debris discs, considering the diverse nature of these systems around stellar spectral types ranging from A to M. Our findings reveal that the prominent ice features at approximately 2.7 and 3.3\,$μ$m depend on both the water ice fraction ${\mathcal{F}}_{\rm ice}$ and the scattering angle, with backscattering geometries yielding the most prominent signatures. When the phase function is considered and data are not background limited, strong forward and backward scattering (near edge-on discs) are expected to yield the strongest detections in images/spectra for A or F-type stars, while scattering angle matters less for later type stars. The Fresnel peak at 3.1\,$μ$m serves as a viable discriminant for the transitional phase (crystalline/amorphous), while simultaneously constraining the water ice temperature. For JWST imaging, we find that the F356W and F444W filter combination is most effective for constraining the grain size distribution, while the F356W and F277W filter combination provides better constraints on the ice fraction ${\mathcal{F}}_{\rm ice}$ in debris discs. However, degeneracy between the grain size distribution and ice fraction when using photometric flux ratios means that obtaining robust constraints will likely require more than two filters, or spectroscopic data. △ Less

Submitted 6 August, 2024; originally announced August 2024.

Comments: 22 pages, 20 figures

Red Type-1 Quasars after Cosmic Noon and Impact on $L_{\rm UV}$-related Quasar Statistics

Authors: Yongjung Kim, Dohyeong Kim, Myungshin Im, Minjin Kim

Abstract: Over the past decades, nearly a million quasars have been explored to shed light on the evolution of supermassive black holes and galaxies. The ultraviolet-to-optical spectra of type-1 quasars particularly offer insights into their black hole activities. Recent findings, however, raise questions about the prevalence of red type-1 quasars of which colors might be due to dust-obscuration and their p… ▽ More Over the past decades, nearly a million quasars have been explored to shed light on the evolution of supermassive black holes and galaxies. The ultraviolet-to-optical spectra of type-1 quasars particularly offer insights into their black hole activities. Recent findings, however, raise questions about the prevalence of red type-1 quasars of which colors might be due to dust-obscuration and their potential influence on luminosity-related properties of quasars. We examine the fraction of red type-1 quasars within the redshift range of $0.68\leq z < 2.20$, applying spectral energy distribution (SED) fitting using optical-to-MIR photometric data of Sloan Digital Sky Survey Data Release 14 quasars. Approximately 10\,\% of the type-1 quasars exhibit red colors suggestive of dust obscuration. There is an association between the brightness of the MIR luminosity and a higher fraction of red type-1 quasars, albeit with negligible redshift evolution. By employing $E(B-V)$ values from the SED fitting, we obtained dereddened luminosity of the red type-1 quasars and reassess the quasar luminosity function (QLF) and black hole mass ($M_{\rm BH}$) estimates. Result shows a modest increase in the number density of bright quasars, linking to more flatten bright-end slope of QLFs, while $M_{\rm BH}$ adjustments are minimal. Current SDSS selections in optical could miss a significant population of heavily dust-obscured quasars. As future MIR surveys like SPHEREx expand, they may reveal enough obscured quasars to prompt a more profound revision of fundamental quasar properties. △ Less

Submitted 8 August, 2024; v1 submitted 6 August, 2024; originally announced August 2024.

Comments: 19 pages, 12 figures, 1 table. Accepted for publication in ApJ

Swift-BAT GUANO follow-up of gravitational-wave triggers in the third LIGO-Virgo-KAGRA observing run

Authors: Gayathri Raman, Samuele Ronchini, James Delaunay, Aaron Tohuvavohu, Jamie A. Kennea, Tyler Parsotan, Elena Ambrosi, Maria Grazia Bernardini, Sergio Campana, Giancarlo Cusumano, Antonino D'Ai, Paolo D'Avanzo, Valerio D'Elia, Massimiliano De Pasquale, Simone Dichiara, Phil Evans, Dieter Hartmann, Paul Kuin, Andrea Melandri, Paul O'Brien, Julian P. Osborne, Kim Page, David M. Palmer, Boris Sbarufatti, Gianpiero Tagliaferri , et al. (1797 additional authors not shown)

Abstract: We present results from a search for X-ray/gamma-ray counterparts of gravitational-wave (GW) candidates from the third observing run (O3) of the LIGO-Virgo-KAGRA (LVK) network using the Swift Burst Alert Telescope (Swift-BAT). The search includes 636 GW candidates received in low latency, 86 of which have been confirmed by the offline analysis and included in the third cumulative Gravitational-Wav… ▽ More We present results from a search for X-ray/gamma-ray counterparts of gravitational-wave (GW) candidates from the third observing run (O3) of the LIGO-Virgo-KAGRA (LVK) network using the Swift Burst Alert Telescope (Swift-BAT). The search includes 636 GW candidates received in low latency, 86 of which have been confirmed by the offline analysis and included in the third cumulative Gravitational-Wave Transient Catalogs (GWTC-3). Targeted searches were carried out on the entire GW sample using the maximum--likelihood NITRATES pipeline on the BAT data made available via the GUANO infrastructure. We do not detect any significant electromagnetic emission that is temporally and spatially coincident with any of the GW candidates. We report flux upper limits in the 15-350 keV band as a function of sky position for all the catalog candidates. For GW candidates where the Swift-BAT false alarm rate is less than 10$^{-3}$ Hz, we compute the GW--BAT joint false alarm rate. Finally, the derived Swift-BAT upper limits are used to infer constraints on the putative electromagnetic emission associated with binary black hole mergers. △ Less

Submitted 13 July, 2024; originally announced July 2024.

Comments: 50 pages, 10 figures, 4 tables

Development of MMC-based lithium molybdate cryogenic calorimeters for AMoRE-II

Authors: A. Agrawal, V. V. Alenkov, P. Aryal, H. Bae, J. Beyer, B. Bhandari, R. S. Boiko, K. Boonin, O. Buzanov, C. R. Byeon, N. Chanthima, M. K. Cheoun, J. S. Choe, S. Choi, S. Choudhury, J. S. Chung, F. A. Danevich, M. Djamal, D. Drung, C. Enss, A. Fleischmann, A. M. Gangapshev, L. Gastaldo, Y. M. Gavrilyuk, A. M. Gezhaev , et al. (84 additional authors not shown)

Abstract: The AMoRE collaboration searches for neutrinoless double beta decay of $^{100}$Mo using molybdate scintillating crystals via low temperature thermal calorimetric detection. The early phases of the experiment, AMoRE-pilot and AMoRE-I, have demonstrated competitive discovery potential. Presently, the AMoRE-II experiment, featuring a large detector array with about 90 kg of $^{100}$Mo isotope, is und… ▽ More The AMoRE collaboration searches for neutrinoless double beta decay of $^{100}$Mo using molybdate scintillating crystals via low temperature thermal calorimetric detection. The early phases of the experiment, AMoRE-pilot and AMoRE-I, have demonstrated competitive discovery potential. Presently, the AMoRE-II experiment, featuring a large detector array with about 90 kg of $^{100}$Mo isotope, is under construction.This paper discusses the baseline design and characterization of the lithium molybdate cryogenic calorimeters to be used in the AMoRE-II detector modules. The results from prototype setups that incorporate new housing structures and two different crystal masses (316 g and 517 - 521 g), operated at 10 mK temperature, show energy resolutions (FWHM) of 7.55 - 8.82 keV at the 2.615 MeV $^{208}$Tl $γ$ line, and effective light detection of 0.79 - 0.96 keV/MeV. The simultaneous heat and light detection enables clear separation of alpha particles with a discrimination power of 12.37 - 19.50 at the energy region around $^6$Li(n, $α$)$^3$H with Q-value = 4.785 MeV. Promising detector performances were demonstrated at temperatures as high as 30 mK, which relaxes the temperature constraints for operating the large AMoRE-II array. △ Less

Submitted 16 July, 2024; originally announced July 2024.

Submillimeter and Mid-Infrared Variability of Young Stellar Objects in the M17SWex Intermediate-Mass Star-Forming Region

Authors: Geumsook Park, Doug Johnstone, Carlos Contreras Pena, Jeong-Eun Lee, Sheng-Yuan Liu, Gregory Herczeg, Steve Mairs, Zhiwei Chen, Jennifer Hatchell, Kee-Tae Kim, Mi-Ryang Kim, Keping Qiu, Yao-Te Wang, Xu Zhang, The JCMT Transient Team

Abstract: We present a comprehensive analysis of young stellar object (YSO) variability within the M17 Southwest Extension (M17 SWex), using 3.5 years of monitoring data from the JCMT Transient Survey at sub-millimeter (sub-mm) and 9 years from the NEOWISE mission at mid-infrared (mid-IR). Our study encompasses observations of 147 bright sub-mm peaks identified within our deep JCMT co-added map as well as 1… ▽ More We present a comprehensive analysis of young stellar object (YSO) variability within the M17 Southwest Extension (M17 SWex), using 3.5 years of monitoring data from the JCMT Transient Survey at sub-millimeter (sub-mm) and 9 years from the NEOWISE mission at mid-infrared (mid-IR). Our study encompasses observations of 147 bright sub-mm peaks identified within our deep JCMT co-added map as well as 156 YSOs in NEOWISE W1 and 179 in W2 that were previously identified in Spitzer surveys. We find three robust sub-mm variables: two are candidate YSOs and one is a likely extragalactic source. At mid-IR wavelengths, our analysis reveals secular and stochastic variability in 47 YSOs, with the highest fraction of secular variability occurring at the earliest evolutionary stage. This is similar to what has previously been observed for low-mass YSO variability within the Gould Belt. However, we observe less overall variability in M17SWex at both the sub-mm and mid-IR. We suspect that this lower fraction is due to the greater distance to M17 SWex. Our findings showcase the utility of multi-wavelength observations to better capture the complex variability phenomena inherent to star formation processes and demonstrate the importance of years-long monitoring of a diverse selection of star-forming environments. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Comments: Accepted for Publication in The Astronomical Journal

The Focal-plane Actualized Shifted Technique Realized for a Shack Hartmann Wavefront Sensor (fastrSHWFS)

Authors: Benjamin L. Gerard, Aaron Lemmer, Bautista R. Fernandez, Xiaoxing Xia, Cesar Laguna, Mike Kim, Stephen Mark Ammons, Brian Bauman, Lisa Poyneer

Abstract: Astronomical adaptive optics (AO) is a critical approach to enable ground-based diffraction-limited imaging and high contrast science, with the potential to enable habitable exoplanet imaging on future extremely large telescopes. However, AO systems must improve significantly to enable habitable exoplanet imaging. Time lag between the end of an exposure and end of deformable mirror commands being… ▽ More Astronomical adaptive optics (AO) is a critical approach to enable ground-based diffraction-limited imaging and high contrast science, with the potential to enable habitable exoplanet imaging on future extremely large telescopes. However, AO systems must improve significantly to enable habitable exoplanet imaging. Time lag between the end of an exposure and end of deformable mirror commands being applied in an AO loop is now the dominant error term in many extreme AO systems (e.g., Poyneer et al. 2016), and within that lag component detector read time is becoming non-negligible (e.g., Cetre et al. 2018). This term will decrease as faster detector readout capabilities are developed by vendors. In complement, we have developed a modified Shack Hartmann Wavefront Sensor (SHWFS) to address this problem called the Focal-plane Actualized Shifted Technique Realized for a SHWFS (fastrSHWFS). The novelty of this design is to replace the usual lenslet array with a bespoke pupil-plane phase mask that redistributes the spot pattern on the detector into a rectangular array with a custom aspect ratio (in an extreme case, if the detector size can accommodate it, the array can be a single line). We present the fastrSHWFS concept and preliminary laboratory tests. For some detectors and AO systems, the fastrSHWFS technique can decrease the read time per frame compared to a regular SHWFS by up to 30x, supporting the goal of reduced AO lag needed to eventually enable habitable exoplanet imaging. △ Less

Submitted 14 June, 2024; originally announced June 2024.

Comments: submitted to SPIE Astronomical Telescopes and Instrumentation, paper number 13097-283

The Size-luminosity Relation of the AGN Torus Determined from the Comparison between Optical and Mid-infrared Variability

Authors: Minjin Kim, Suyeon Son, Luis C. Ho

Abstract: We investigate the optical variability of low-redshift ($0.15< z\leq0.4$) active galactic nuclei using the multi-epoch data from the Zwicky Transient Facility. We find that a damped random walk model well describes the ensemble structure function in the $g$ band. Consistent with previous studies, more luminous active galactic nuclei tend to have a steeper structure function at a timescale less tha… ▽ More We investigate the optical variability of low-redshift ($0.15< z\leq0.4$) active galactic nuclei using the multi-epoch data from the Zwicky Transient Facility. We find that a damped random walk model well describes the ensemble structure function in the $g$ band. Consistent with previous studies, more luminous active galactic nuclei tend to have a steeper structure function at a timescale less than the break timescale and smaller variability amplitude. By comparing the structure functions in the optical with the mid-infrared obtained from the Wide-field Infrared Survey Explorer, we derive the size of the dusty torus using a toy model for the geometry of the torus. The size of the torus positively correlates with the luminosity of the active nucleus, following a relation that agrees well with previous studies based on reverberation mapping. This result demonstrates that the structure function method can be used as a powerful and highly efficient tool to examine the size of the torus. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 7 pages, 5 figures. Accepted for publication in A&A

The quantum-refractive evolution of polarization states in pulsar emission

Authors: Dong-Hoon Kim, Chul Min Kim, Sang Pyo Kim

Abstract: Highly magnetized neutron stars have quantum refraction effects on pulsar emission due to the non-linearity of the quantum electrodynamics (QED) action. In this paper, we investigate the evolution of the polarization states under the quantum refraction effects combined with the frequency dependence of pulsar emission; we solve a system of evolution equations of the Stokes vector, where the birefri… ▽ More Highly magnetized neutron stars have quantum refraction effects on pulsar emission due to the non-linearity of the quantum electrodynamics (QED) action. In this paper, we investigate the evolution of the polarization states under the quantum refraction effects combined with the frequency dependence of pulsar emission; we solve a system of evolution equations of the Stokes vector, where the birefringent vector, in which such effects are encoded, acts on the Stokes vector. At a fixed frequency of emission, depending on the magnitude of the birefringent vector, dominated mostly by the magnetic field strength, the evolution of the Stokes vector largely exhibits three different patterns: (i) monotonic, or (ii) half-oscillatory, or (iii) highly oscillatory behaviours. These features are understood and confirmed by means of approximate analytical solutions to the evolution equations. △ Less

Submitted 9 June, 2024; originally announced June 2024.

Comments: 14 pages, 9 figures

Modeling the refractive index profile n(z) of polar ice for ultra-high energy neutrino experiments

Authors: S. Ali, P. Allison, S. Archambault, J. J. Beatty, D. Z. Besson, A. Bishop, P. Chen, Y. C. Chen, B. A. Clark, W. Clay, A. Connolly, K. Couberly, L. Cremonesi, A. Cummings, P. Dasgupta, R. Debolt, S. de Kockere, K. D. de Vries, C. Deaconu, M. A. DuVernois, J. Flaherty, E. Friedman, R. Gaior, P. Giri, J. Hanson , et al. (45 additional authors not shown)

Abstract: We develop an in-situ index of refraction profile using the transit time of radio signals broadcast from an englacial transmitter to 2-5 km distant radio-frequency receivers, deployed at depths up to 200 m. Maxwell's equations generally admit two ray propagation solutions from a given transmitter, corresponding to a direct path (D) and a refracted path (R); the measured D vs. R (dt(D,R)) timing di… ▽ More We develop an in-situ index of refraction profile using the transit time of radio signals broadcast from an englacial transmitter to 2-5 km distant radio-frequency receivers, deployed at depths up to 200 m. Maxwell's equations generally admit two ray propagation solutions from a given transmitter, corresponding to a direct path (D) and a refracted path (R); the measured D vs. R (dt(D,R)) timing differences provide constraints on the index of refraction profile near South Pole, where the Askaryan Radio Array (ARA) neutrino observatory is located. We constrain the refractive index profile by simulating D and R ray paths via ray tracing and comparing those to measured dt(D,R) signals. Using previous ice density data as a proxy for n(z), we demonstrate that our data strongly favors a glaciologically-motivated three-phase densification model rather than a single exponential scale height model. Simulations show that the single exponential model overestimates ARA neutrino sensitivity compared to the three-phase model. △ Less

Submitted 11 June, 2024; v1 submitted 2 June, 2024; originally announced June 2024.

MHONGOOSE discovery of a gas-rich low-surface brightness galaxy in the Dorado Group

Authors: F. M. Maccagni, W. J. G. de Blok, P. E. Mancera Piña, R. Ragusa, E. Iodice, M. Spavone, S. McGaugh, K. A. Oman, T. A. Oosterloo, B. S. Koribalski, M. Kim, E. A. K. Adams, P. Amram, A. Bosma, F. Bigiel, E. Brinks, L. Chemin, F. Combes, B. Gibson, J. Healy, B. W. Holwerda, G. I. G. Józsa, P. Kamphuis, D. Kleiner, S. Kurapati , et al. (6 additional authors not shown)

Abstract: We present the discovery of a low-mass gas-rich low-surface brightness galaxy in the Dorado Group, at a distance of 17.7 Mpc. Combining deep MeerKAT 21-cm observations from the MeerKAT HI Observations of Nearby Galactic Objects: Observing Southern Emitters (MHONGOOSE) survey with deep photometric images from the VST Early-type Galaxy Survey (VEGAS) we find a stellar and neutral atomic hydrogen (HI… ▽ More We present the discovery of a low-mass gas-rich low-surface brightness galaxy in the Dorado Group, at a distance of 17.7 Mpc. Combining deep MeerKAT 21-cm observations from the MeerKAT HI Observations of Nearby Galactic Objects: Observing Southern Emitters (MHONGOOSE) survey with deep photometric images from the VST Early-type Galaxy Survey (VEGAS) we find a stellar and neutral atomic hydrogen (HI) gas mass of $M_\star = 2.23\times10^6$ M$_\odot$ and $M_{\rm HI}=1.68\times10^6$ M$_\odot$, respectively. This low-surface brightness galaxy is the lowest mass HI detection found in a group beyond the Local Universe ($D\gtrsim 10$ Mpc). The dwarf galaxy has the typical overall properties of gas-rich low surface brightness galaxies in the Local group, but with some striking differences. Namely, the MHONGOOSE observations reveal a very low column density ($\sim 10^{18-19}$ cm$^{-2}$) HI disk with asymmetrical morphology possibly supported by rotation and higher velocity dispersion in the centre. There, deep optical photometry and UV-observations suggest a recent enhancement of the star formation. Found at galactocentric distances where in the Local Group dwarf galaxies are depleted of cold gas (at $390$ projected-kpc distance from the group centre), this galaxy is likely on its first orbit within the Dorado group. We discuss the possible environmental effects that may have caused the formation of the HI disk and the enhancement of star formation, highlighting the short-lived phase (a few hundreds of Myr) of the gaseous disk, before either SF or hydrodynamical forces will deplete the gas of the galaxy. △ Less

Submitted 27 May, 2024; originally announced May 2024.

Comments: 18 pages, 14 figures, 5 tables; accepted for publication in A&A

Observation of Gravitational Waves from the Coalescence of a $2.5\text{-}4.5~M_\odot$ Compact Object and a Neutron Star

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, I. Abouelfettouh, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, D. Agarwal, M. Agathos, M. Aghaei Abchouyeh, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, S. Akçay, T. Akutsu, S. Albanesi, R. A. Alfaidi, A. Al-Jodah , et al. (1771 additional authors not shown)

Abstract: We report the observation of a coalescing compact binary with component masses $2.5\text{-}4.5~M_\odot$ and $1.2\text{-}2.0~M_\odot$ (all measurements quoted at the 90% credible level). The gravitational-wave signal GW230529_181500 was observed during the fourth observing run of the LIGO-Virgo-KAGRA detector network on 2023 May 29 by the LIGO Livingston Observatory. The primary component of the so… ▽ More We report the observation of a coalescing compact binary with component masses $2.5\text{-}4.5~M_\odot$ and $1.2\text{-}2.0~M_\odot$ (all measurements quoted at the 90% credible level). The gravitational-wave signal GW230529_181500 was observed during the fourth observing run of the LIGO-Virgo-KAGRA detector network on 2023 May 29 by the LIGO Livingston Observatory. The primary component of the source has a mass less than $5~M_\odot$ at 99% credibility. We cannot definitively determine from gravitational-wave data alone whether either component of the source is a neutron star or a black hole. However, given existing estimates of the maximum neutron star mass, we find the most probable interpretation of the source to be the coalescence of a neutron star with a black hole that has a mass between the most massive neutron stars and the least massive black holes observed in the Galaxy. We provisionally estimate a merger rate density of $55^{+127}_{-47}~\text{Gpc}^{-3}\,\text{yr}^{-1}$ for compact binary coalescences with properties similar to the source of GW230529_181500; assuming that the source is a neutron star-black hole merger, GW230529_181500-like sources constitute about 60% of the total merger rate inferred for neutron star-black hole coalescences. The discovery of this system implies an increase in the expected rate of neutron star-black hole mergers with electromagnetic counterparts and provides further evidence for compact objects existing within the purported lower mass gap. △ Less

Submitted 26 July, 2024; v1 submitted 5 April, 2024; originally announced April 2024.

Comments: 45 pages (10 pages author list, 13 pages main text, 1 page acknowledgements, 13 pages appendices, 8 pages bibliography), 17 figures, 16 tables. Update to match version published in The Astrophysical Journal Letters. Data products available from https://zenodo.org/records/10845779

Report number: LIGO-P2300352

Journal ref: ApJL 970, L34 (2024)

Ultralight vector dark matter search using data from the KAGRA O3GK run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi , et al. (1778 additional authors not shown)

Abstract: Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we prese… ▽ More Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we present the result of a search for $U(1)_{B-L}$ gauge boson DM using the KAGRA data from auxiliary length channels during the first joint observation run together with GEO600. By applying our search pipeline, which takes into account the stochastic nature of ultralight DM, upper bounds on the coupling strength between the $U(1)_{B-L}$ gauge boson and ordinary matter are obtained for a range of DM masses. While our constraints are less stringent than those derived from previous experiments, this study demonstrates the applicability of our method to the lower-mass vector DM search, which is made difficult in this measurement by the short observation time compared to the auto-correlation time scale of DM. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 20 pages, 5 figures

Report number: LIGO-P2300250

Filamentary Network and Magnetic Field Structures Revealed with BISTRO in the High-Mass Star-Forming Region NGC2264 : Global Properties and Local Magnetogravitational Configurations

Authors: Jia-Wei Wang, Patrick M. Koch, Seamus D. Clarke, Gary Fuller, Nicolas Peretto, Ya-Wen Tang, Hsi-Wei Yen, Shih-Ping Lai, Nagayoshi Ohashi, Doris Arzoumanian, Doug Johnstone, Ray Furuya, Shu-ichiro Inutsuka, Chang Won Lee, Derek Ward-Thompson, Valentin J. M. Le Gouellec, Hong-Li Liu, Lapo Fanciullo, Jihye Hwang, Kate Pattle, Frédérick Poidevin, Mehrnoosh Tahani, Takashi Onaka, Mark G. Rawlings, Eun Jung Chung , et al. (132 additional authors not shown)

Abstract: We report 850 $μ$m continuum polarization observations toward the filamentary high-mass star-forming region NGC 2264, taken as part of the B-fields In STar forming Regions Observations (BISTRO) large program on the James Clerk Maxwell Telescope (JCMT). These data reveal a well-structured non-uniform magnetic field in the NGC 2264C and 2264D regions with a prevailing orientation around 30 deg from… ▽ More We report 850 $μ$m continuum polarization observations toward the filamentary high-mass star-forming region NGC 2264, taken as part of the B-fields In STar forming Regions Observations (BISTRO) large program on the James Clerk Maxwell Telescope (JCMT). These data reveal a well-structured non-uniform magnetic field in the NGC 2264C and 2264D regions with a prevailing orientation around 30 deg from north to east. Field strengths estimates and a virial analysis for the major clumps indicate that NGC 2264C is globally dominated by gravity while in 2264D magnetic, gravitational, and kinetic energies are roughly balanced. We present an analysis scheme that utilizes the locally resolved magnetic field structures, together with the locally measured gravitational vector field and the extracted filamentary network. From this, we infer statistical trends showing that this network consists of two main groups of filaments oriented approximately perpendicular to one another. Additionally, gravity shows one dominating converging direction that is roughly perpendicular to one of the filament orientations, which is suggestive of mass accretion along this direction. Beyond these statistical trends, we identify two types of filaments. The type-I filament is perpendicular to the magnetic field with local gravity transitioning from parallel to perpendicular to the magnetic field from the outside to the filament ridge. The type-II filament is parallel to the magnetic field and local gravity. We interpret these two types of filaments as originating from the competition between radial collapsing, driven by filament self-gravity, and the longitudinal collapsing, driven by the region's global gravity. △ Less

Submitted 23 January, 2024; originally announced January 2024.

Comments: Accepted for publication in the Astrophysical Journal. 43 pages, 32 figures, and 4 tables (including Appendix)

Quantum refraction effects in pulsar radio emission

Authors: Dong-Hoon Kim, Chul Min Kim, Sang Pyo Kim

Abstract: Highly magnetized neutron stars exhibit the vacuum non-linear electrodynamics effects, which can be well described using the one-loop effective action for quantum electrodynamics. In this context, we study the propagation and polarization of pulsar radio emission, based on the post-Maxwellian Lagrangian from the Heisenberg-Euler-Schwinger action. Given the refractive index obtained from this Lagra… ▽ More Highly magnetized neutron stars exhibit the vacuum non-linear electrodynamics effects, which can be well described using the one-loop effective action for quantum electrodynamics. In this context, we study the propagation and polarization of pulsar radio emission, based on the post-Maxwellian Lagrangian from the Heisenberg-Euler-Schwinger action. Given the refractive index obtained from this Lagrangian, we determine the leading-order corrections to both the propagation and polarization vectors due to quantum refraction via perturbation analysis. In addition, the effects on the orthogonality between the propagation and polarization vectors and the Faraday rotation angle, all due to quantum refraction are investigated. Furthermore, from the dual refractive index and the associated polarization modes, we discuss quantum birefringence, with the optical phenomenology analogous to its classical counterpart. △ Less

Submitted 11 June, 2024; v1 submitted 8 January, 2024; originally announced January 2024.

Comments: 14 pages, 5 figures

Journal ref: Published in Monthly Notices of the Royal Astronomical Society, 531, 2148 on 20 May 2024

The JCMT Transient Survey: Six-Year Summary of 450/850\,$μ$m Protostellar Variability and Calibration Pipeline Version 2.0

Authors: Steve Mairs, Seonjae Lee, Doug Johnstone, Colton Broughton, Jeong-Eun Lee, Gregory J. Herczeg, Graham S. Bell, Zhiwei Chen, Carlos Contreras-Peña, Logan Francis, Jennifer Hatchell, Mi-Ryang Kim, Sheng-Yuan Liu, Geumsook Park, Keping Qiu, Yao-Te Wang, Xu Zhang, The JCMT Transient Team

Abstract: The JCMT Transient Survey has been monitoring eight Gould Belt low-mass star-forming regions since December 2015 and six somewhat more distant intermediate-mass star-forming regions since February 2020 with SCUBA-2 on the JCMT at \ShortS and \LongS and with an approximately monthly cadence. We introduce our Pipeline v2 relative calibration procedures for image alignment and flux calibration across… ▽ More The JCMT Transient Survey has been monitoring eight Gould Belt low-mass star-forming regions since December 2015 and six somewhat more distant intermediate-mass star-forming regions since February 2020 with SCUBA-2 on the JCMT at \ShortS and \LongS and with an approximately monthly cadence. We introduce our Pipeline v2 relative calibration procedures for image alignment and flux calibration across epochs, improving on our previous Pipeline v1 by decreasing measurement uncertainties and providing additional robustness. These new techniques work at both \LongS and \ShortNS, where v1 only allowed investigation of the \LongS data. Pipeline v2 achieves better than $0.5^{\prime\prime}$ relative image alignment, less than a tenth of the submillimeter beam widths. The v2 relative flux calibration is found to be 1\% at \LongS and $<5$\% at \ShortNS. The improvement in the calibration is demonstrated by comparing the two pipelines over the first four years of the survey and recovering additional robust variables with v2. Using the full six years of the Gould Belt survey the number of robust variables increases by 50\,\%, and at \ShortS we identify four robust variables, all of which are also robust at \LongNS. The multi-wavelength light curves for these sources are investigated and found to be consistent with the variability being due to dust heating within the envelope in response to accretion luminosity changes from the central source. △ Less

Submitted 7 January, 2024; originally announced January 2024.

Comments: Accepted for Publication in the The Astrophysical Journal. DOI link to data will become public after the proof stage is complete

Impact of discontinuous grain size distributions on the spectral energy distribution of debris disks

Authors: Minjae Kim, Sebastian Wolf

Abstract: The collisional evolution of debris disks is expected to result in a characteristic wavy pattern of the grain size distributions, i.e., an under/overabundance of particles of specific sizes. This perturbed grain size distribution potentially leaves characteristic patterns in the spectral energy distribution (SED) of the disk system. We aim to quantify and understand the specific influence of disco… ▽ More The collisional evolution of debris disks is expected to result in a characteristic wavy pattern of the grain size distributions, i.e., an under/overabundance of particles of specific sizes. This perturbed grain size distribution potentially leaves characteristic patterns in the spectral energy distribution (SED) of the disk system. We aim to quantify and understand the specific influence of discontinuous particle size distributions on the appearance of debris disks. For this purpose, we consider dust emission models based on two different grain size distributions, i.e., once with a single and once with a broken power law. We compare the spectral index $α$ ($F_ν\,\propto ν^{\rm{\,α}}$) in the case of a continuous grain size distribution with that of a discontinuous grain size distribution. We perform this comparison for central stars with different spectral types and two different disk structures (e.g., slim and broad debris dust rings). Within the considered parameter space, we find a characteristic difference between the spectral slopes of the SED in the different scenarios. More specifically, the overabundance of small grains leads to a steeper slope in the far-infrared/submillimeter regime, while the spectral index in the mm regime is hardly affected. On the other hand, the underabundance of medium-sized grains results in a slight steepening of the far-infrared slope of SED, while its primary effect is on the mm slope of SED, causing it to become shallower. We also find that the impact of an overabundance of small dust particles is more pronounced than that of an underabundance of medium-sized dust particles. Furthermore, we find that the difference between the spectral indices for the two different grain size distributions is largest for debris disks around brighter central stars and broader disks. △ Less

Submitted 23 December, 2023; originally announced December 2023.

Comments: 14 pages, 10 figures

The JWST Early Release Science Program for Direct Observations of Exoplanetary Systems V: Do Self-Consistent Atmospheric Models Represent JWST Spectra? A Showcase With VHS 1256 b

Authors: Simon Petrus, Niall Whiteford, Polychronis Patapis, Beth A. Biller, Andrew Skemer, Sasha Hinkley, Genaro Suárez, Anna Lueber, Paulina Palma-Bifani, Jordan M. Stone, Johanna M. Vos, Caroline V. Morley, Pascal Tremblin, Benjamin Charnay, Christiane Helling, Brittany E. Miles, Aarynn L. Carter, Jason J. Wang, Markus Janson, Eileen C. Gonzales, Ben Sutlieff, Kielan K. W. Hoch, Mickaël Bonnefoy, Gaël Chauvin, Olivier Absil , et al. (97 additional authors not shown)

Abstract: The unprecedented medium-resolution (R~1500-3500) near- and mid-infrared (1-18um) spectrum provided by JWST for the young (140+/-20Myr) low-mass (12-20MJup) L-T transition (L7) companion VHS1256b gives access to a catalogue of molecular absorptions. In this study, we present a comprehensive analysis of this dataset utilizing a forward modelling approach, applying our Bayesian framework, ForMoSA. W… ▽ More The unprecedented medium-resolution (R~1500-3500) near- and mid-infrared (1-18um) spectrum provided by JWST for the young (140+/-20Myr) low-mass (12-20MJup) L-T transition (L7) companion VHS1256b gives access to a catalogue of molecular absorptions. In this study, we present a comprehensive analysis of this dataset utilizing a forward modelling approach, applying our Bayesian framework, ForMoSA. We explore five distinct atmospheric models to assess their performance in estimating key atmospheric parameters: Teff, log(g), [M/H], C/O, gamma, fsed, and R. Our findings reveal that each parameter's estimate is significantly influenced by factors such as the wavelength range considered and the model chosen for the fit. This is attributed to systematic errors in the models and their challenges in accurately replicating the complex atmospheric structure of VHS1256b, notably the complexity of its clouds and dust distribution. To propagate the impact of these systematic uncertainties on our atmospheric property estimates, we introduce innovative fitting methodologies based on independent fits performed on different spectral windows. We finally derived a Teff consistent with the spectral type of the target, considering its young age, which is confirmed by our estimate of log(g). Despite the exceptional data quality, attaining robust estimates for chemical abundances [M/H] and C/O, often employed as indicators of formation history, remains challenging. Nevertheless, the pioneering case of JWST's data for VHS1256b has paved the way for future acquisitions of substellar spectra that will be systematically analyzed to directly compare the properties of these objects and correct the systematics in the models. △ Less

Submitted 31 January, 2024; v1 submitted 6 December, 2023; originally announced December 2023.

Comments: 32 pages, 16 figures, 6 tables, 2 appendices

The Seoul National University AGN Monitoring Project III: H$β$ lag measurements of 32 luminous AGNs and the high-luminosity end of the size--luminosity relation

Authors: Jong-Hak Woo, Shu Wang, Suvendu Rakshit, Hojin Cho, Donghoon Son, Vardha N. Bennert, Elena Gallo, Edmund Hodges-Kluck, Tommaso Treu, Aaron J. Barth, Wanjin Cho, Adi Foord, Jaehyuk Geum, Hengxiao Guo, Yashashree Jadhav, Yiseul Jeon, Kyle M. Kabasares, Won-Suk Kang, Changseok Kim, Minjin Kim, Tae-Woo Kim, Huynh Anh N. Le, Matthew A. Malkan, Amit Kumar Mandal, Daeseong Park , et al. (6 additional authors not shown)

Abstract: We present the main results from a long-term reverberation mapping campaign carried out for the Seoul National University Active Galactic Nuclei (AGN) Monitoring Project. High-quality data were obtained during 2015-2021 for 32 luminous AGNs (i.e., continuum luminosity in the range of $10^{44-46}$ erg s$^{-1}$) at a regular cadence, of 20-30 days for spectroscopy and 3-5 days for photometry. We obt… ▽ More We present the main results from a long-term reverberation mapping campaign carried out for the Seoul National University Active Galactic Nuclei (AGN) Monitoring Project. High-quality data were obtained during 2015-2021 for 32 luminous AGNs (i.e., continuum luminosity in the range of $10^{44-46}$ erg s$^{-1}$) at a regular cadence, of 20-30 days for spectroscopy and 3-5 days for photometry. We obtain time lag measurements between the variability in the H$β$ emission and the continuum for 32 AGNs; twenty-five of those have the best lag measurements based on our quality assessment, examining correlation strength, and the posterior lag distribution. Our study significantly increases the current sample of reverberation-mapped AGNs, particularly at the moderate to high luminosity end. Combining our results with literature measurements, we derive a H$β$ broad line region size--luminosity relation with a shallower slope than reported in the literature. For a given luminosity, most of our measured lags are shorter than the expectation, implying that single-epoch black hole mass estimators based on previous calibrations could suffer large systematic uncertainties. △ Less

Submitted 26 November, 2023; originally announced November 2023.

Comments: Accepted by ApJ; 39 pages, 22 figures

Photometry of the Didymos system across the DART impact apparition

Authors: Nicholas Moskovitz, Cristina Thomas, Petr Pravec, Tim Lister, Tom Polakis, David Osip, Theodore Kareta, Agata Rożek, Steven R. Chesley, Shantanu P. Naidu, Peter Scheirich, William Ryan, Eileen Ryan, Brian Skiff, Colin Snodgrass, Matthew M. Knight, Andrew S. Rivkin, Nancy L. Chabot, Vova Ayvazian, Irina Belskaya, Zouhair Benkhaldoun, Daniel N. Berteşteanu, Mariangela Bonavita, Terrence H. Bressi, Melissa J. Brucker , et al. (56 additional authors not shown)

Abstract: On 26 September 2022, the Double Asteroid Redirection Test (DART) spacecraft impacted Dimorphos, the satellite of binary near-Earth asteroid (65803) Didymos. This demonstrated the efficacy of a kinetic impactor for planetary defense by changing the orbital period of Dimorphos by 33 minutes (Thomas et al. 2023). Measuring the period change relied heavily on a coordinated campaign of lightcurve phot… ▽ More On 26 September 2022, the Double Asteroid Redirection Test (DART) spacecraft impacted Dimorphos, the satellite of binary near-Earth asteroid (65803) Didymos. This demonstrated the efficacy of a kinetic impactor for planetary defense by changing the orbital period of Dimorphos by 33 minutes (Thomas et al. 2023). Measuring the period change relied heavily on a coordinated campaign of lightcurve photometry designed to detect mutual events (occultations and eclipses) as a direct probe of the satellite's orbital period. A total of 28 telescopes contributed 224 individual lightcurves during the impact apparition from July 2022 to February 2023. We focus here on decomposable lightcurves, i.e. those from which mutual events could be extracted. We describe our process of lightcurve decomposition and use that to release the full data set for future analysis. We leverage these data to place constraints on the post-impact evolution of ejecta. The measured depths of mutual events relative to models showed that the ejecta became optically thin within the first ~1 day after impact, and then faded with a decay time of about 25 days. The bulk magnitude of the system showed that ejecta no longer contributed measurable brightness enhancement after about 20 days post-impact. This bulk photometric behavior was not well represented by an HG photometric model. An HG1G2 model did fit the data well across a wide range of phase angles. Lastly, we note the presence of an ejecta tail through at least March 2023. Its persistence implied ongoing escape of ejecta from the system many months after DART impact. △ Less

Submitted 3 November, 2023; originally announced November 2023.

Comments: 52 pages, 5 tables, 9 figures, accepted to PSJ

Scaling slowly rotating asteroids by stellar occultations

Authors: A. Marciniak, J. Ďurech, A. Choukroun, J. Hanuš, W. Ogłoza, R. Szakáts, L. Molnár, A. Pál, F. Monteiro, E. Frappa, W. Beisker, H. Pavlov, J. Moore, R. Adomavičienė, R. Aikawa, S. Andersson, P. Antonini, Y. Argentin, A. Asai, P. Assoignon, J. Barton, P. Baruffetti, K. L. Bath, R. Behrend, L. Benedyktowicz , et al. (154 additional authors not shown)

Abstract: As evidenced by recent survey results, majority of asteroids are slow rotators (P>12 h), but lack spin and shape models due to selection bias. This bias is skewing our overall understanding of the spins, shapes, and sizes of asteroids, as well as of their other properties. Also, diameter determinations for large (>60km) and medium-sized asteroids (between 30 and 60 km) often vary by over 30% for m… ▽ More As evidenced by recent survey results, majority of asteroids are slow rotators (P>12 h), but lack spin and shape models due to selection bias. This bias is skewing our overall understanding of the spins, shapes, and sizes of asteroids, as well as of their other properties. Also, diameter determinations for large (>60km) and medium-sized asteroids (between 30 and 60 km) often vary by over 30% for multiple reasons. Our long-term project is focused on a few tens of slow rotators with periods of up to 60 hours. We aim to obtain their full light curves and reconstruct their spins and shapes. We also precisely scale the models, typically with an accuracy of a few percent. We used wide sets of dense light curves for spin and shape reconstructions via light-curve inversion. Precisely scaling them with thermal data was not possible here because of poor infrared data: large bodies are too bright for WISE mission. Therefore, we recently launched a campaign among stellar occultation observers, to scale these models and to verify the shape solutions, often allowing us to break the mirror pole ambiguity. The presented scheme resulted in shape models for 16 slow rotators, most of them for the first time. Fitting them to stellar occultations resolved previous inconsistencies in size determinations. For around half of the targets, this fitting also allowed us to identify a clearly preferred pole solution, thus removing the ambiguity inherent to light-curve inversion. We also address the influence of the uncertainty of the shape models on the derived diameters. Overall, our project has already provided reliable models for around 50 slow rotators. Such well-determined and scaled asteroid shapes will, e.g. constitute a solid basis for density determinations when coupled with mass information. Spin and shape models continue to fill the gaps caused by various biases. △ Less

Submitted 13 October, 2023; originally announced October 2023.

Comments: Accepted to Astronomy & Astrophysics. 12 pages + appendices

Journal ref: A&A 679, A60 (2023)

The Structure Function of Mid-infrared Variability in Low-redshift Active Galactic Nuclei

Authors: Suyeon Son, Minjin Kim, Luis C. Ho

Abstract: Using the multi-epoch mid-infrared (MIR) photometry from the Wide-field Infrared Survey Explorer spanning a baseline of $\sim10$ yr, we extensively investigate the MIR variability of nearby active galactic nuclei (AGNs) at $0.15 < z < 0.4$. We find that the ensemble structure function in the W1 band ($3.4\ μ$m) can be modeled with a broken power law. Type 1 AGNs tend to exhibit larger variability… ▽ More Using the multi-epoch mid-infrared (MIR) photometry from the Wide-field Infrared Survey Explorer spanning a baseline of $\sim10$ yr, we extensively investigate the MIR variability of nearby active galactic nuclei (AGNs) at $0.15 < z < 0.4$. We find that the ensemble structure function in the W1 band ($3.4\ μ$m) can be modeled with a broken power law. Type 1 AGNs tend to exhibit larger variability amplitudes than type 2 AGNs, possibly due to the extinction by the torus. The variability amplitude is inversely correlated with the AGN luminosity, consistent with a similar relation known in the optical. Meanwhile, the slope of the power law increases with AGN luminosity. This trend can be attributed to the fact that the inner radius of the torus is proportional to the AGN luminosity, as expected from the size$-$luminosity relation of the torus. Interestingly, low-luminosity type 2 AGNs, unlike low-luminosity type 1 AGNs, tend to exhibit smaller variability amplitude than do high-luminosity AGNs. We argue that either low-luminosity type 2 AGNs have distinctive central structures due to their low luminosity or their MIR brightness is contaminated by emission from the cold dust in the host galaxy. Our findings suggest that the AGN unification scheme may need to be revised. We find that the variability amplitude of dust-deficient AGNs is systematically larger than that of normal AGNs, supporting the notion that the hot and warm dust in dust-deficient AGNs may be destroyed and reformed according to the strength of the ultraviolet radiation from the accretion disk. △ Less

Submitted 10 October, 2023; v1 submitted 8 October, 2023; originally announced October 2023.

Comments: Accepted for publication in ApJ, 17 pages, 13 figures

Protoplanetary and debris disks in the $η$ Chamaeleontis Association: A sub-millimeter survey obtained with APEX/LABOCA

Authors: V. Roccatagliata, A. Sicilia-Aguilar, M. Kim, J. Campbell-White, M. Fang, S. J. Murphy, S. Wolf, W. A. Lawson, Th. Henning, J. Bouwman

Abstract: Nearby associations are ideal regions to study coeval samples of protoplanetary and debris disks down to late M-type stars. Those aged 5-10,Myrs, where most of the disk should have already dissipated forming planets, are of particular interest. We present the first complete study of both protoplanetary and debris disks in a young region, using the $η$ Chamaeleontis association as a test bench to s… ▽ More Nearby associations are ideal regions to study coeval samples of protoplanetary and debris disks down to late M-type stars. Those aged 5-10,Myrs, where most of the disk should have already dissipated forming planets, are of particular interest. We present the first complete study of both protoplanetary and debris disks in a young region, using the $η$ Chamaeleontis association as a test bench to study the cold disk content. We obtained sub-millimeter data for the entire core population down to late M-type stars, plus a few halo members. We performed a continuum sub-millimeter survey with APEX/LABOCA of all the core populations of $η$ Cha association. Disk properties have been derived by modeling protoplanetary and debris disks using RADMC 2D and DMS, respectively. We find that protoplanetary disks in $η$ Cha typically have holes with radii of the order of 0.01 to 0.03 AU, while ring-like emission from the debris disks is located between 20 au and 650 au from the central star. The parallaxes and Gaia eDR3 photometry, in combination with the PARSEC and COLIBRI isochrones, enable us to confirm an age of $η$ Cha between 7 and 9 Myrs. In general, the disk mass seems insufficient to support accretion over a long time, even for the lowest mass accretors, a clear difference compared with other regions and also a sign that the mass budget is further underestimated. We do not find a correlation between the stellar masses, accretion rates, and disk masses, although this could be due to sample issues. We confirm that the presence of inner holes is not enough to stop accretion unless accompanied by dramatic changes to the total disk mass content. Comparing $η$ Cha with other regions at different ages, we find that the physical processes responsible for debris disks (e.g., dust growth, dust trapping) efficiently act in less than 5 Myrs. △ Less

Submitted 4 October, 2023; originally announced October 2023.

Comments: 20 pages, 9 figures, 7 tables

Estimators of Bolometric Luminosity and Black Hole Mass with Mid-infrared Continuum Luminosities for Dust-obscured Quasars: Prevalence of Dust-obscured SDSS Quasars

Authors: Dohyeong Kim, Myungshin Im, Minjin Kim, Yongjung Kim, Suhyun Shin, Hyunjin Shim, Hyunmi Song

Abstract: We present bolometric luminosity ($L_{\rm bol}$) and black hole (BH) mass ($M_{\rm BH}$) estimators based on mid-infrared (MIR) continuum luminosity (hereafter, $L_{\rm MIR}$) that are measured from infrared (IR) photometric data. The $L_{\rm MIR}$-based estimators are relatively immune from dust extinction effects, hence they can be used for dust-obscured quasars. To derive the $L_{\rm bol}$ and… ▽ More We present bolometric luminosity ($L_{\rm bol}$) and black hole (BH) mass ($M_{\rm BH}$) estimators based on mid-infrared (MIR) continuum luminosity (hereafter, $L_{\rm MIR}$) that are measured from infrared (IR) photometric data. The $L_{\rm MIR}$-based estimators are relatively immune from dust extinction effects, hence they can be used for dust-obscured quasars. To derive the $L_{\rm bol}$ and $M_{\rm BH}$ estimators, we use unobscured quasars selected from the Sloan Digital Sky Survey (SDSS) quasar catalog, which have wide ranges of $L_{\rm bol}$ ($10^{44.62}$--$10^{46.16}$\,$\rm erg\,s^{-1}$) and $M_{\rm BH}$ ($10^{7.14}$--$10^{9.69}$\,$M_{\odot}$). We find empirical relations between (i) continuum luminosity at 5100\,$\rm{Å}$ (hereafter, L5100) and $L_{\rm MIR}$; (ii) $L_{\rm bol}$ and $L_{\rm MIR}$. Using these relations, we derive the $L_{\rm MIR}$-based $L_{\rm bol}$ and $M_{\rm BH}$ estimators. We find that our estimators allow the determination of $L_{\rm bol}$ and $M_{\rm BH}$ at an accuracy of $\sim$0.2\,dex against the fiducial estimates based on the optical properties of the unobscured quasars. We apply the $L_{\rm MIR}$-based estimators to SDSS quasars at $z \lesssim 0.5$ including obscured ones. The ratios of $L_{\rm bol}$ from the $L_{\rm MIR}$-based estimators to those from the optical luminosity-based estimators become larger with the amount of the dust extinction, and a non-negligible fraction ($\sim$15\,\%) of the SDSS quasars exhibits ratios greater than 1.5. This result suggests that dust extinction can significantly affect physical parameter derivations even for SDSS quasars, and that dust extinction needs to be carefully taken into account when deriving quasar properties. △ Less

Submitted 2 October, 2023; originally announced October 2023.

Comments: 17 pages, 10 figures, 3 tables, Published in ApJ

Journal ref: 2023ApJ...954..156K

Schwinger Pair Production and Vacuum Birefringence around High Magnetized Neutron Stars

Authors: Chul Min Kim, Sang Pyo Kim

Abstract: Highly magnetized neutron stars have magnetic fields of order of the critical field and can lead to measurable QED effects. We consider the Goldreich-Julian pulsar model with supercritical magnetic fields, induced subcritical electric fields, and a period of milliseconds. We then study the strong field physics, such as Schwinger pair production and the vacuum birefringence including the wrench eff… ▽ More Highly magnetized neutron stars have magnetic fields of order of the critical field and can lead to measurable QED effects. We consider the Goldreich-Julian pulsar model with supercritical magnetic fields, induced subcritical electric fields, and a period of milliseconds. We then study the strong field physics, such as Schwinger pair production and the vacuum birefringence including the wrench effect, whose X-ray polarimetry will be observed in future space missions. △ Less

Submitted 30 August, 2023; originally announced August 2023.

Comments: 13 pages, 6 figures; Talk at the Fifth Zeldovich meeting, an international conference in honor of Ya. B. Zeldovich held in Yerevan, Armenia on June 12--16, 2023

A Joint Fermi-GBM and Swift-BAT Analysis of Gravitational-Wave Candidates from the Third Gravitational-wave Observing Run

Authors: C. Fletcher, J. Wood, R. Hamburg, P. Veres, C. M. Hui, E. Bissaldi, M. S. Briggs, E. Burns, W. H. Cleveland, M. M. Giles, A. Goldstein, B. A. Hristov, D. Kocevski, S. Lesage, B. Mailyan, C. Malacaria, S. Poolakkil, A. von Kienlin, C. A. Wilson-Hodge, The Fermi Gamma-ray Burst Monitor Team, M. Crnogorčević, J. DeLaunay, A. Tohuvavohu, R. Caputo, S. B. Cenko , et al. (1674 additional authors not shown)

Abstract: We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM on-board triggers and sub-threshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses,… ▽ More We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM on-board triggers and sub-threshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses, the Targeted Search and the Untargeted Search, we investigate whether there are any coincident GRBs associated with the GWs. We also search the Swift-BAT rate data around the GW times to determine whether a GRB counterpart is present. No counterparts are found. Using both the Fermi-GBM Targeted Search and the Swift-BAT search, we calculate flux upper limits and present joint upper limits on the gamma-ray luminosity of each GW. Given these limits, we constrain theoretical models for the emission of gamma-rays from binary black hole mergers. △ Less

Submitted 25 August, 2023; originally announced August 2023.

Calibration and Physics with ARA Station 1: A Unique Askaryan Radio Array Detector

Authors: M. F. H Seikh, D. Z. Besson, S. Ali, P. Allison, S. Archambault, J. J. Beatty, A. Bishop, P. Chen, Y. C. Chen, B. A. Clark, W. Clay, A. Connolly, K. Couberly, L. Cremonesi, A. Cummings, P. Dasgupta, R. Debolt, S. De Kockere, K. D. de Vries, C. Deaconu, M. A. DuVernois, J. Flaherty, E. Friedman, R. Gaior, P. Giri , et al. (48 additional authors not shown)

Abstract: The Askaryan Radio Array Station 1 (A1), the first among five autonomous stations deployed for the ARA experiment at the South Pole, is a unique ultra-high energy neutrino (UHEN) detector based on the Askaryan effect that uses Antarctic ice as the detector medium. Its 16 radio antennas (distributed across 4 strings, each with 2 Vertically Polarized (VPol), 2 Horizontally Polarized (HPol) receivers… ▽ More The Askaryan Radio Array Station 1 (A1), the first among five autonomous stations deployed for the ARA experiment at the South Pole, is a unique ultra-high energy neutrino (UHEN) detector based on the Askaryan effect that uses Antarctic ice as the detector medium. Its 16 radio antennas (distributed across 4 strings, each with 2 Vertically Polarized (VPol), 2 Horizontally Polarized (HPol) receivers), and 2 strings of transmitting antennas (calibration pulsers, CPs), each with 1 VPol and 1 HPol channel, are deployed at depths less than 100 m within the shallow firn zone of the 2.8 km thick South Pole (SP) ice. We apply different methods to calibrate its Ice Ray Sampler second generation (IRS2) chip for timing offset and ADC-to-Voltage conversion factors using a known continuous wave input signal to the digitizer, and achieve a precision of sub-nanoseconds. We achieve better calibration for odd, compared to even samples, and also find that the HPols under-perform relative to the VPol channels. Our timing calibrated data is subsequently used to calibrate the ADC-to-Voltage conversion as well as precise antenna locations, as a precursor to vertex reconstruction. The calibrated data will then be analyzed for UHEN signals in the final step of data compression. The ability of A1 to scan the firn region of SP ice sheet will contribute greatly towards a 5-station analysis and will inform the design of the planned IceCube Gen-2 radio array. △ Less

Submitted 14 August, 2023; originally announced August 2023.

Comments: 10 pages

Journal ref: PoS ICRC2023 (2023) 1163

Cometary dust collected by MIDAS on board Rosetta II. Particle shape descriptors and pristineness evaluation

Authors: M. Kim, T. Mannel, J. Lasue, A. Longobardo, M. S. Bentley, R. Moissl, the MIDAS team

Abstract: The MIDAS (Micro-Imaging Dust Analysis System) atomic force microscope on board the Rosetta comet orbiter investigated and measured the 3D topography of a few hundred nm to tens of $μ$m sized dust particles of 67P/Churyumov-Gerasimenko with resolutions down to a few nanometers, giving insights into the physical processes of our early Solar System. We analyze the shapes of the cometary dust particl… ▽ More The MIDAS (Micro-Imaging Dust Analysis System) atomic force microscope on board the Rosetta comet orbiter investigated and measured the 3D topography of a few hundred nm to tens of $μ$m sized dust particles of 67P/Churyumov-Gerasimenko with resolutions down to a few nanometers, giving insights into the physical processes of our early Solar System. We analyze the shapes of the cometary dust particles collected by MIDAS on the basis of a recently updated particle catalog with the aim to determine which structural properties remained pristine. We develop a set of shape descriptors and metrics such as aspect ratio, elongation, circularity, convexity, and particle surface/volume distribution, which can be used to describe the distribution of particle shapes. Furthermore, we compare the structure of the MIDAS dust particles and the clusters in which the particles were deposited to those found in previous laboratory experiments and by Rosetta/COSIMA. Finally, we combine our findings to calculate a pristineness score for MIDAS particles and determine the most pristine particles and their properties. We find that the morphological properties of all cometary dust particles at the micrometer scale are surprisingly homogeneous despite originating from diverse cometary environments (e.g., different collection targets that are associated with cometary activities/source regions and collection velocities/periods). We next find that the types of clusters found by MIDAS show good agreement with those defined by previous laboratory experiments, however, there are some differences to those found by Rosetta/COSIMA. Based on our result, we rate 19 out of 1082 MIDAS particles at least moderately pristine, i.e., they are not substantially flattened by impact, not fragmented, and/or not part of a fragmentation cluster. △ Less

Submitted 17 August, 2023; v1 submitted 10 August, 2023; originally announced August 2023.

Comments: 40 pages, 31 figures, 1 online table

Journal ref: A&A 678, A179 (2023)

On the origin of optical and IR emission lines in star forming galaxies

Authors: Mariela Martínez-Paredes, Gustavo Bruzual, Christophe Morisset, Minsun Kim, Marcio Meléndez, Luc Binette

Abstract: Combining the {\sc Cloudy} photoionization code with updated stellar population synthesis results, we simultaneously model the MIR $\neiii/\neii$ vs. $\oiv/\neiii$, the MIR-FIR $\neiii/\neii$ vs. $\oiv/\oiii$ and the classical BPT diagnostic diagrams. We focus on the properties of optically classified \hii\,galaxies that lie in the normal star forming zone in the MIR diagnostic diagram. We find th… ▽ More Combining the {\sc Cloudy} photoionization code with updated stellar population synthesis results, we simultaneously model the MIR $\neiii/\neii$ vs. $\oiv/\neiii$, the MIR-FIR $\neiii/\neii$ vs. $\oiv/\oiii$ and the classical BPT diagnostic diagrams. We focus on the properties of optically classified \hii\,galaxies that lie in the normal star forming zone in the MIR diagnostic diagram. We find that a small fraction of our models lie in this zone, but most of them correspond to the lowest explored metallicity, \zstar\,=\,0.0002, at age $\sim1$ Gyr. This value of \zstar\,is, by far, lower than the values derived for these galaxies from optical emission lines, suggesting that the far-UV emission produced by post-AGB stars (a.k.a. HOLMES, hot low-mass evolved stars) is NOT the source of ionization. Instead, shock models can easily reproduce this part of the MIR diagram. We suggest that it is likely that some of these galaxies have been misclassified and that in them, shocks, produced by a weak AGN-outflow, could be an important source of ionizaton. Using a subset of our models, we derive a new demarcation line for the maximal contribution of retired galaxies in the BPT diagram. This demarcation line allows for a larger contamination from the neighbouring AGN-dominated region. Considering the importance of disentangling the different ionising mechanisms in weak or deeply obscured systems, new observational efforts to classify galaxies both in the optical and IR are required to better constrain this kind of models and understand their evolutionary paths. △ Less

Submitted 10 August, 2023; originally announced August 2023.

Comments: Accepted for publication in MNRAS

Search for Eccentric Black Hole Coalescences during the Third Observing Run of LIGO and Virgo

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi , et al. (1750 additional authors not shown)

Abstract: Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effect… ▽ More Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effects of eccentricity. Here, we present observational results for a waveform-independent search sensitive to eccentric black hole coalescences, covering the third observing run (O3) of the LIGO and Virgo detectors. We identified no new high-significance candidates beyond those that were already identified with searches focusing on quasi-circular binaries. We determine the sensitivity of our search to high-mass (total mass $M>70$ $M_\odot$) binaries covering eccentricities up to 0.3 at 15 Hz orbital frequency, and use this to compare model predictions to search results. Assuming all detections are indeed quasi-circular, for our fiducial population model, we place an upper limit for the merger rate density of high-mass binaries with eccentricities $0 < e \leq 0.3$ at $0.33$ Gpc$^{-3}$ yr$^{-1}$ at 90\% confidence level. △ Less

Submitted 7 August, 2023; originally announced August 2023.

Comments: 24 pages, 5 figures

Report number: LIGO-P2300080

Photometric Selection of Unobscured QSOs in the Ecliptic Poles: KMTNet in the South Field and Pan-STARRS in the North Field

Authors: Woowon Byun, Minjin Kim, Yun-Kyeong Sheen, Dongseob Lee, Luis C. Ho, Jongwan Ko, Kwang-Il Seon, Hyunjin Shim, Dohyeong Kim, Yongjung Kim, Joon Hyeop Lee, Hyunjin Jeong, Jong-Hak Woo, Woong-Seob Jeong, Byeong-Gon Park, Sang Chul Kim, Yongseok Lee, Sang-Mok Cha, Hyunmi Song, Donghoon Son, Yujin Yang

Abstract: We search for quasi-stellar objects (QSOs) in a wide area of the south ecliptic pole (SEP) field, which has been and will continue to be intensively explored through various space missions. For this purpose, we obtain deep broadband optical images of the SEP field covering an area of $\sim$$14.5\times14.5$ deg$^2$ with the Korea Microlensing Telescope Network. The 5$σ$ detection limits for point s… ▽ More We search for quasi-stellar objects (QSOs) in a wide area of the south ecliptic pole (SEP) field, which has been and will continue to be intensively explored through various space missions. For this purpose, we obtain deep broadband optical images of the SEP field covering an area of $\sim$$14.5\times14.5$ deg$^2$ with the Korea Microlensing Telescope Network. The 5$σ$ detection limits for point sources in the $BVRI$ bands are estimated to be $\sim$22.59, 22.60, 22.98, and 21.85 mag, respectively. Utilizing data from Wide-field Infrared Survey Explorer, unobscured QSO candidates are selected among the optically point-like sources using the mid-infrared (MIR) and optical-MIR colors. To further refine our selection and eliminate any contamination not adequately removed by the color-based selection, we perform the spectral energy distribution fitting with archival photometric data ranging from optical to MIR. As a result, we identify a total of 2,383 unobscured QSO candidates in the SEP field. We also apply a similar method to the north ecliptic pole field using the Pan-STARRS data and obtain a similar result of identifying 2,427 candidates. The differential number count per area of our QSO candidates is in good agreement with those measured from spectroscopically confirmed ones in other fields. Finally, we compare the results with the literature and discuss how this work will be implicated in future studies, especially with the upcoming space missions. △ Less

Submitted 28 July, 2023; originally announced July 2023.

Comments: 14 pages, 9 figures, accepted for publication in ApJS

The Seoul National University AGN Monitoring Project IV: H$α$ reverberation mapping of 6 AGNs and the H$α$ Size-Luminosity Relation

Authors: Hojin Cho, Jong-Hak Woo, Shu Wang, Donghoon Son, Jaejin Shin, Suvendu Rakshit, Aaron J. Barth, Vardha N. Bennert, Elena Gallo, Edmund Hodges-Kluck, Tommaso Treu, Hyun-Jin Bae, Wanjin Cho, Adi Foord, Jaehyuk Geum, Yashashree Jadhav, Yiseul Jeon, Kyle M. Kabasares, Daeun Kang, Wonseok Kang, Changseok Kim, Donghwa Kim, Minjin Kim, Taewoo Kim, Huynh Anh N. Le , et al. (7 additional authors not shown)

Abstract: The broad line region (BLR) size-luminosity relation has paramount importance for estimating the mass of black holes in active galactic nuclei (AGNs). Traditionally, the size of the H$β$ BLR is often estimated from the optical continuum luminosity at 5100\angstrom{} , while the size of the H$α$ BLR and its correlation with the luminosity is much less constrained. As a part of the Seoul National Un… ▽ More The broad line region (BLR) size-luminosity relation has paramount importance for estimating the mass of black holes in active galactic nuclei (AGNs). Traditionally, the size of the H$β$ BLR is often estimated from the optical continuum luminosity at 5100\angstrom{} , while the size of the H$α$ BLR and its correlation with the luminosity is much less constrained. As a part of the Seoul National University AGN Monitoring Project (SAMP) which provides six-year photometric and spectroscopic monitoring data, we present our measurements of the H$α$ lags of 6 high-luminosity AGNs. Combined with the measurements for 42 AGNs from the literature, we derive the size-luminosity relations of H$α$ BLR against broad H$α$ and 5100\angstrom{} continuum luminosities. We find the slope of the relations to be $0.61\pm0.04$ and $0.59\pm0.04$, respectively, which are consistent with the \hb{} size-luminosity relation. Moreover, we find a linear relation between the 5100\angstrom{} continuum luminosity and the broad H$α$ luminosity across 7 orders of magnitude. Using these results, we propose a new virial mass estimator based on the H$α$ broad emission line, finding that the previous mass estimates based on the scaling relations in the literature are overestimated by up to 0.7 dex at masses lower than $10^7$~M$_{\odot}$. △ Less

Submitted 29 June, 2023; originally announced June 2023.

Comments: Accepted for publication in ApJ (Jun. 25th, 2023). 21 pages, 12 figures

Deficit of Hot Dust in Low-redshift Active Galactic Nuclei

Authors: Suyeon Son, Minjin Kim, Luis C. Ho

Abstract: We assemble a broad-band spectral energy distribution (SED) ranging from optical to mid-infrared of nearby active galactic nuclei at $z < 0.4$. SED fitting analysis is performed using semi-empirical templates derived from Palomar-Green quasars to classify the sample into normal, warm-dust-deficient (WDD), and hot-dust-deficient (HDD) AGNs. Kolmogorov-Smirnov tests reveal that HDD AGNs exhibit, on… ▽ More We assemble a broad-band spectral energy distribution (SED) ranging from optical to mid-infrared of nearby active galactic nuclei at $z < 0.4$. SED fitting analysis is performed using semi-empirical templates derived from Palomar-Green quasars to classify the sample into normal, warm-dust-deficient (WDD), and hot-dust-deficient (HDD) AGNs. Kolmogorov-Smirnov tests reveal that HDD AGNs exhibit, on average higher AGN luminosity than normal and WDD AGNs. HDD fraction, on the other hand, is only weakly correlated with black hole mass and inversely correlated with Eddington ratio. By fixing the other parameters, we conclude that the HDD fraction is primarily connected with the AGN luminosity. It implies that there is a causal connection between the covering factor of the hot dust component and AGN luminosity, possibly due to the sublimation of the innermost dust or the thickening of the intervening gas in the broad-line region. Analysis of the outflow properties traced by the wing of [O III]$\lambda5007$ suggests that outflows may be related to the formation and maintenance of the hot dust component. Finally, we demonstrate through comparison with previous studies that the classification of HDD AGNs requires careful subtraction of the host galaxy light. △ Less

Submitted 22 June, 2023; originally announced June 2023.

Comments: Accepted for publication in ApJ

2021 occultations and transits of Linus orbiting (22) Kalliope: I. Polygonal and `cliptracing' algorithm

Authors: M. Brož, J. Ďurech, M. Ferrais, H. -J. Lee, M. -J. Kim, D. -G. Roh, H. -S. Yim, E. Jehin, A. Burdanov, J. de Wit, P. Fatka, J. Hanuš, B. Carry

Abstract: The satellite Linus orbiting the main-belt asteroid (22) Kalliope exhibited occultation and transit events in late 2021. A photometric campaign was organized and observations were taken by the TRAPPIST-South, SPECULOOS-Artemis, OWL-Net, and BOAO telescopes, with the goal to constrain models of this system. Our dynamical model is complex, with multipoles (up to the order $\ell = 2$), internal tides… ▽ More The satellite Linus orbiting the main-belt asteroid (22) Kalliope exhibited occultation and transit events in late 2021. A photometric campaign was organized and observations were taken by the TRAPPIST-South, SPECULOOS-Artemis, OWL-Net, and BOAO telescopes, with the goal to constrain models of this system. Our dynamical model is complex, with multipoles (up to the order $\ell = 2$), internal tides, and external tides. The model was constrained by astrometry (spanning 2001--2021), occultations, adaptive-optics imaging, calibrated photometry, as well as relative photometry. Our photometric model was substantially improved. A new precise (${<}\,0.1\,{\rm mmag}$) light curve algorithm was implemented, based on polygon intersections, which are computed exactly -- by including partial eclipses and partial visibility of polygons. Moreover, we implemented a `cliptracing' algorithm, based again on polygon intersections, in which partial contributions to individual pixels are computed exactly. Both synthetic light curves and synthetic images are then very smooth. Based on our combined solution, we confirmed the size of Linus, $(28\pm 1)\,{\rm km}$. However, this solution exhibits some tension between the light curves and the PISCO speckle-interferometry dataset. In most solutions, Linus is darker than Kalliope, with the albedos $A_{\rm w} = 0.40$ vs. $0.44$. This is confirmed on deconvolved images. A~detailed revision of astrometric data allowed us to revise also the $J_2 \equiv -C_{20}$ value of Kalliope. Most importantly, a~homogeneous body is excluded. For a differentiated body, two solutions exist: low-oblateness ($C_{20} \simeq -0.12$), with a~spherical iron core, and alternatively, high-oblateness ($C_{20} \simeq -0.22$) with an elongated iron core. These correspond to the low- and high-energy collisions, respectively, studied by means of SPH simulations in our previous work. △ Less

Submitted 7 June, 2023; originally announced June 2023.

Comments: Astronomy and Astrophysics, accepted

Journal ref: A&A 676, A60 (2023)

The JCMT BISTRO Survey: Studying the Complex Magnetic Field of L43

Authors: Janik Karoly, Derek Ward-Thompson, Kate Pattle, David Berry, Anthony Whitworth, Jason Kirk, Pierre Bastien, Tao-Chung Ching, Simon Coude, Jihye Hwang, Woojin Kwon, Archana Soam, Jia-Wei Wang, Tetsuo Hasegawa, Shih-Ping Lai, Keping Qiu, Doris Arzoumanian, Tyler L. Bourke, Do-Young Byun, Huei-Ru Vivien Chen, Wen Ping Chen, Mike Chen, Zhiwei Chen, Jungyeon Cho, Minho Choi , et al. (133 additional authors not shown)

Abstract: We present observations of polarized dust emission at 850 $μ$m from the L43 molecular cloud which sits in the Ophiuchus cloud complex. The data were taken using SCUBA-2/POL-2 on the James Clerk Maxwell Telescope as a part of the BISTRO large program. L43 is a dense ($N_{\rm H_2}\sim 10^{22}$-10$^{23}$ cm$^{-2}$) complex molecular cloud with a submillimetre-bright starless core and two protostellar… ▽ More We present observations of polarized dust emission at 850 $μ$m from the L43 molecular cloud which sits in the Ophiuchus cloud complex. The data were taken using SCUBA-2/POL-2 on the James Clerk Maxwell Telescope as a part of the BISTRO large program. L43 is a dense ($N_{\rm H_2}\sim 10^{22}$-10$^{23}$ cm$^{-2}$) complex molecular cloud with a submillimetre-bright starless core and two protostellar sources. There appears to be an evolutionary gradient along the isolated filament that L43 is embedded within, with the most evolved source closest to the Sco OB2 association. One of the protostars drives a CO outflow that has created a cavity to the southeast. We see a magnetic field that appears to be aligned with the cavity walls of the outflow, suggesting interaction with the outflow. We also find a magnetic field strength of up to $\sim$160$\pm$30 $μ$G in the main starless core and up to $\sim$90$\pm$40 $μ$G in the more diffuse, extended region. These field strengths give magnetically super- and sub-critical values respectively and both are found to be roughly trans-Alfvénic. We also present a new method of data reduction for these denser but fainter objects like starless cores. △ Less

Submitted 22 May, 2023; v1 submitted 18 May, 2023; originally announced May 2023.

Comments: Accepted for publication in ApJ. 23 pages, 9 figures (7 main text, 2 appendix)

Coronal Heating as Determined by the Solar Flare Frequency Distribution Obtained by Aggregating Case Studies

Authors: James Paul Mason, Alexandra Werth, Colin G. West, Allison A. Youngblood, Donald L. Woodraska, Courtney Peck, Kevin Lacjak, Florian G. Frick, Moutamen Gabir, Reema A. Alsinan, Thomas Jacobsen, Mohammad Alrubaie, Kayla M. Chizmar, Benjamin P. Lau, Lizbeth Montoya Dominguez, David Price, Dylan R. Butler, Connor J. Biron, Nikita Feoktistov, Kai Dewey, N. E. Loomis, Michal Bodzianowski, Connor Kuybus, Henry Dietrick, Aubrey M. Wolfe , et al. (977 additional authors not shown)

Abstract: Flare frequency distributions represent a key approach to addressing one of the largest problems in solar and stellar physics: determining the mechanism that counter-intuitively heats coronae to temperatures that are orders of magnitude hotter than the corresponding photospheres. It is widely accepted that the magnetic field is responsible for the heating, but there are two competing mechanisms th… ▽ More Flare frequency distributions represent a key approach to addressing one of the largest problems in solar and stellar physics: determining the mechanism that counter-intuitively heats coronae to temperatures that are orders of magnitude hotter than the corresponding photospheres. It is widely accepted that the magnetic field is responsible for the heating, but there are two competing mechanisms that could explain it: nanoflares or Alfvén waves. To date, neither can be directly observed. Nanoflares are, by definition, extremely small, but their aggregate energy release could represent a substantial heating mechanism, presuming they are sufficiently abundant. One way to test this presumption is via the flare frequency distribution, which describes how often flares of various energies occur. If the slope of the power law fitting the flare frequency distribution is above a critical threshold, $α=2$ as established in prior literature, then there should be a sufficient abundance of nanoflares to explain coronal heating. We performed $>$600 case studies of solar flares, made possible by an unprecedented number of data analysts via three semesters of an undergraduate physics laboratory course. This allowed us to include two crucial, but nontrivial, analysis methods: pre-flare baseline subtraction and computation of the flare energy, which requires determining flare start and stop times. We aggregated the results of these analyses into a statistical study to determine that $α= 1.63 \pm 0.03$. This is below the critical threshold, suggesting that Alfvén waves are an important driver of coronal heating. △ Less

Submitted 9 May, 2023; originally announced May 2023.

Comments: 1,002 authors, 14 pages, 4 figures, 3 tables, published by The Astrophysical Journal on 2023-05-09, volume 948, page 71

Physical properties of the slow-rotating near-Earth asteroid (2059) Baboquivari from one apparition

Authors: Orhan Erece, Irek Khamitov, Murat Kaplan, Yucel Kilic, Hee-Jae Lee, Myung-Jin Kim, Ilfan F. Bikmaev, Rustem I. Gumerov, Eldar N. Irtuganov

Abstract: In this study, we carried out photometric, spectroscopic, and for the first time, polarimetric observations of the Amor-type near-Earth asteroid (2059) Baboquivari. Our findings represent the first reliable determination of Baboquivari's physical properties. We used data from a 1m-class telescope (T100) along with ALCDEF data for photometric analyses and a 1.5-m-class telescope (RTT150) for polari… ▽ More In this study, we carried out photometric, spectroscopic, and for the first time, polarimetric observations of the Amor-type near-Earth asteroid (2059) Baboquivari. Our findings represent the first reliable determination of Baboquivari's physical properties. We used data from a 1m-class telescope (T100) along with ALCDEF data for photometric analyses and a 1.5-m-class telescope (RTT150) for polarimetric, spectroscopic, and additional photometric observations. We obtained the synodic rotation period of Baboquivari as 129.93 +/- 2.31 hours and the standard phase function parameters H and G as 16.05 +/- 0.05, 0.22 +/- 0.02, respectively. Our colour index (V-R) measurement of 0.45 +/- 0.02 is consistent with spectroscopic observations, indicating an S (or sub-S) spectral type. Using the polarimetric and spectroscopic data, we found that the geometric albedo is 0.15 +/- 0.03, and the spectral type is Sq. Based on the estimated albedo and absolute magnitude, Baboquivari has an effective diameter of 2.12 +/- 0.21 km. Due to the scattered data in the light curve, its slow rotation and location among the NEAs suggest that Baboquivari may be a non-principal axis (NPA) rotator. △ Less

Submitted 9 May, 2023; originally announced May 2023.

Comments: 17 pages, 9 figures, accepted to publish in Planetary and Space Science (2023)

First BISTRO observations of the dark cloud Taurus L1495A-B10: the role of the magnetic field in the earliest stages of low-mass star formation

Authors: Derek Ward-Thompson, Janik Karoly, Kate Pattle, Anthony Whitworth, Jason Kirk, David Berry, Pierre Bastien, Tao-Chung Ching, Simon Coude, Jihye Hwang, Woojin Kwon, Archana Soam, Jia-Wei Wang, Tetsuo Hasegawa, Shih-Ping Lai, Keping Qiu, Doris Arzoumanian, Tyler L. Bourke, Do-Young Byun, Huei-Ru Vivien Chen, Wen Ping Chen, Mike Chen, Zhiwei Chen, Jungyeon Cho, Minho Choi , et al. (133 additional authors not shown)

Abstract: We present BISTRO Survey 850 μm dust emission polarisation observations of the L1495A-B10 region of the Taurus molecular cloud, taken at the JCMT. We observe a roughly triangular network of dense filaments. We detect 9 of the dense starless cores embedded within these filaments in polarisation, finding that the plane-of-sky orientation of the core-scale magnetic field lies roughly perpendicular to… ▽ More We present BISTRO Survey 850 μm dust emission polarisation observations of the L1495A-B10 region of the Taurus molecular cloud, taken at the JCMT. We observe a roughly triangular network of dense filaments. We detect 9 of the dense starless cores embedded within these filaments in polarisation, finding that the plane-of-sky orientation of the core-scale magnetic field lies roughly perpendicular to the filaments in almost all cases. We also find that the large-scale magnetic field orientation measured by Planck is not correlated with any of the core or filament structures, except in the case of the lowest-density core. We propose a scenario for early prestellar evolution that is both an extension to, and consistent with, previous models, introducing an additional evolutionary transitional stage between field-dominated and matter-dominated evolution, observed here for the first time. In this scenario, the cloud collapses first to a sheet-like structure. Uniquely, we appear to be seeing this sheet almost face-on. The sheet fragments into filaments, which in turn form cores. However, the material must reach a certain critical density before the evolution changes from being field-dominated to being matter-dominated. We measure the sheet surface density and the magnetic field strength at that transition for the first time and show consistency with an analytical prediction that had previously gone untested for over 50 years (Mestel 1965). △ Less

Submitted 23 February, 2023; originally announced February 2023.

Comments: 14 pages, 5 figures. ApJ accepted

Cometary dust collected by MIDAS on board Rosetta. I. Dust particle catalog and statistics

Authors: M. Kim, T. Mannel, P. D. Boakes, M. S. Bentley, A. Longobardo, H. Jeszenszky, R. Moissl, the MIDAS team

Abstract: We aim to catalog all dust particles collected and analyzed by MIDAS, together with their main statistical properties such as size, height, basic shape descriptors, and collection time. Furthermore, we aim to present the scientific results that can be extracted from the catalog (e.g., the size distribution and statistical characteristics of cometary dust particles). The existing MIDAS particle cat… ▽ More We aim to catalog all dust particles collected and analyzed by MIDAS, together with their main statistical properties such as size, height, basic shape descriptors, and collection time. Furthermore, we aim to present the scientific results that can be extracted from the catalog (e.g., the size distribution and statistical characteristics of cometary dust particles). The existing MIDAS particle catalog has been greatly improved by a careful re-analysis of the AFM images, leading to the addition of more dust particles and a detailed description of the particle properties. The catalog documents all images of identified dust particles and includes a variety of derived information tabulated one record per particle. Furthermore, the best image of each particle was chosen for subsequent studies. Finally, we created dust coverage maps and clustering maps of the MIDAS collection targets and traced any possible fragmentation of collected particles with a detailed algorithm. The revised MIDAS catalog includes 3523 MIDAS particles in total, where 1857 particles are expected to be usable for further analysis (418 scans of particles before perihelion + 1439 scans of particles after perihelion, both after the removal of duplicates), ranging from about 40 nm to about 8 $μ$m in size. The mean value of the equivalent radius derived from the 2D projection of the particles is 0.91 ${\pm}$ 0.79 $μ$m. A slightly improved equivalent radius based on the particle's volume coincides in the range of uncertainties with a value of 0.56 ${\pm}$ 0.45 $μ$m. We note that those sizes and all following MIDAS particle size distributions are expected to be influenced by the fragmentation of MIDAS particles upon impact on the collection targets. Furthermore, fitting the slope of the MIDAS particle size distribution with a power law of a r ${^b}$ yields an index b of ${\sim}$ -1.67 to -1.88. △ Less

Submitted 17 March, 2023; v1 submitted 21 February, 2023; originally announced February 2023.

Comments: 28 pages, 5 figure, 2 online tables

Report number: 673, A129

Journal ref: A&A 673, A129 (2023)

Flybys in debris disk systems with Gaia eDR3

Authors: Leonardo Bertini, Veronica Roccatagliata, Minjae Kim

Abstract: We aim to observationally and statistically constrain the influence of flybys in the formation and evolution of debris disks. We compiled a sample of 254 debris disks with ages between 2 Myr and 8 Gyr that are either part of an association or isolated, drawing the binary and planetary companions of the systems mainly from the literature. Using the Gaia eDR3 astrometric data and radial velocities o… ▽ More We aim to observationally and statistically constrain the influence of flybys in the formation and evolution of debris disks. We compiled a sample of 254 debris disks with ages between 2 Myr and 8 Gyr that are either part of an association or isolated, drawing the binary and planetary companions of the systems mainly from the literature. Using the Gaia eDR3 astrometric data and radial velocities of our sample, as well as all the sources in a specific region of the sky, we reconstructed the relative linear motions in the last 5 Myr and made predictions for the next 2 Myr. Relating the Hill radius of each debris disk system and the closest distances reached by the two sources, we defined the flyby events in terms of position and time. We find that in the period between the last 5 Myrs and the next 2 Myrs, 90% of the analyzed systems have experienced at least a close flyby, while 7% of them have experienced flybys at distances greater than 0.5R Hill. In particular, 75% of them have experienced at least one past close encounter and 36% multiple past close encounters. From the sub-sample of resolved debris disk (41 out of 94), 80% of the analyzed systems experience at least an encounter within 0.8 pc. From the subsample of 10 debris disks with planets, half of these systems do show misalignments between disk and planet, stirring, or asymmetries. Systems with a misalignment between the planetary orbit and the disk do indeed experience at least one flyby event. In particular, when the planet orbits have a difference with the disk inclination higher than about 20 degree, as in the case of HD 38529, we find that multiple close encounters have taken place in the last 5 Myr, as theoretically predicted. The high incidence of encounters, particularly close encounters, experienced by the systems in the last 5 Myr suggests the fundamental impact of flybys on the evolution of debris disks. △ Less

Submitted 13 February, 2023; originally announced February 2023.

Comments: 9 pages, 2 figures

Journal ref: A&A 671, L2 (2023)

JCMT BISTRO Observations: Magnetic Field Morphology of Bubbles Associated with NGC 6334

Authors: Mehrnoosh Tahani, Pierre Bastien, Ray S. Furuya, Kate Pattle, Doug Johnstone, Doris Arzoumanian, Yasuo Doi, Tetsuo Hasegawa, Shu-ichiro Inutsuka, Simon Coudé, Laura Fissel, Michael Chun-Yuan Chen, Frédérick Poidevin, Sarah Sadavoy, Rachel Friesen, Patrick M. Koch, James Di Francesco, Gerald H. Moriarty-Schieven, Zhiwei Chen, Eun Jung Chung, Chakali Eswaraiah, Lapo Fanciullo, Tim Gledhill, Valentin J. M. Le Gouellec, Thiem Hoang , et al. (120 additional authors not shown)

Abstract: We study the HII regions associated with the NGC 6334 molecular cloud observed in the sub-millimeter and taken as part of the B-fields In STar-forming Region Observations (BISTRO) Survey. In particular, we investigate the polarization patterns and magnetic field morphologies associated with these HII regions. Through polarization pattern and pressure calculation analyses, several of these bubbles… ▽ More We study the HII regions associated with the NGC 6334 molecular cloud observed in the sub-millimeter and taken as part of the B-fields In STar-forming Region Observations (BISTRO) Survey. In particular, we investigate the polarization patterns and magnetic field morphologies associated with these HII regions. Through polarization pattern and pressure calculation analyses, several of these bubbles indicate that the gas and magnetic field lines have been pushed away from the bubble, toward an almost tangential (to the bubble) magnetic field morphology. In the densest part of NGC 6334, where the magnetic field morphology is similar to an hourglass, the polarization observations do not exhibit observable impact from HII regions. We detect two nested radial polarization patterns in a bubble to the south of NGC 6334 that correspond to the previously observed bipolar structure in this bubble. Finally, using the results of this study, we present steps (incorporating computer vision; circular Hough Transform) that can be used in future studies to identify bubbles that have physically impacted magnetic field lines. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: Accepted for publication in Astrophysical Journal (ApJ)

Metallicity-PAH Relation of MIR-selected Star-forming Galaxies in AKARI North Ecliptic Pole-wide Survey

Authors: Hyunjin Shim, Ho Seong Hwang, Woong-Seob Jeong, Yoshiki Toba, Minjin Kim, Dohyeong Kim, Hyunmi Song, Tetsuya Hashimoto, Takao Nakagawa, Ambra Nanni, William J. Pearson, Toshinobu Takagi

Abstract: We investigate the variation in the mid-infrared spectral energy distributions of 373 low-redshift ($z<0.4$) star-forming galaxies, which reflects a variety of polycyclic aromatic hydrocarbon (PAH) emission features. The relative strength of PAH emission is parameterized as $q_\mathrm{PAH}$, which is defined as the mass fraction of PAH particles in the total dust mass. With the aid of continuous m… ▽ More We investigate the variation in the mid-infrared spectral energy distributions of 373 low-redshift ($z<0.4$) star-forming galaxies, which reflects a variety of polycyclic aromatic hydrocarbon (PAH) emission features. The relative strength of PAH emission is parameterized as $q_\mathrm{PAH}$, which is defined as the mass fraction of PAH particles in the total dust mass. With the aid of continuous mid-infrared photometric data points covering 7-24$μ$m and far-infrared flux densities, $q_\mathrm{PAH}$ values are derived through spectral energy distribution fitting. The correlation between $q_\mathrm{PAH}$ and other physical properties of galaxies, i.e., gas-phase metallicity ($12+\mathrm{log(O/H)}$), stellar mass, and specific star-formation rate (sSFR) are explored. As in previous studies, $q_\mathrm{PAH}$ values of galaxies with high metallicity are found to be higher than those with low metallicity. The strength of PAH emission is also positively correlated with the stellar mass and negatively correlated with the sSFR. The correlation between $q_\mathrm{PAH}$ and each parameter still exists even after the other two parameters are fixed. In addition to the PAH strength, the application of metallicity-dependent gas-to-dust mass ratio appears to work well to estimate gas mass that matches the observed relationship between molecular gas and physical parameters. The result obtained will be used to calibrate the observed PAH luminosity-total infrared luminosity relation, based on the variation of MIR-FIR SED, which is used in the estimation of hidden star formation. △ Less

Submitted 8 December, 2022; originally announced December 2022.

Comments: 19 pages, 9 figures, AJ, in press

The JCMT BISTRO-2 Survey: Magnetic Fields of the Massive DR21 Filament

Authors: Tao-Chung Ching, Keping Qiu, Di Li, Zhiyuan Ren, Shih-Ping Lai, David Berry, Kate Pattle, Ray Furuya, Derek Ward-Thompson, Doug Johnstone, Patrick M. Koch, Chang Won Lee, Thiem Hoang, Tetsuo Hasegawa, Woojin Kwon, Pierre Bastien, Chakali Eswaraiah, Jia-Wei Wang, Kyoung Hee Kim, Jihye Hwang, Archana Soam, A-Ran Lyo, Junhao Liu, Valentin J. M. Le Gouellec, Doris Arzoumanian , et al. (132 additional authors not shown)

Abstract: We present 850 $μ$m dust polarization observations of the massive DR21 filament from the B-fields In STar-forming Region Observations (BISTRO) survey, using the POL-2 polarimeter and the SCUBA-2 camera on the James Clerk Maxwell Telescope. We detect ordered magnetic fields perpendicular to the parsec-scale ridge of the DR21 main filament. In the sub-filaments, the magnetic fields are mainly parall… ▽ More We present 850 $μ$m dust polarization observations of the massive DR21 filament from the B-fields In STar-forming Region Observations (BISTRO) survey, using the POL-2 polarimeter and the SCUBA-2 camera on the James Clerk Maxwell Telescope. We detect ordered magnetic fields perpendicular to the parsec-scale ridge of the DR21 main filament. In the sub-filaments, the magnetic fields are mainly parallel to the filamentary structures and smoothly connect to the magnetic fields of the main filament. We compare the POL-2 and Planck dust polarization observations to study the magnetic field structures of the DR21 filament on 0.1--10 pc scales. The magnetic fields revealed in the Planck data are well aligned with those of the POL-2 data, indicating a smooth variation of magnetic fields from large to small scales. The plane-of-sky magnetic field strengths derived from angular dispersion functions of dust polarization are 0.6--1.0 mG in the DR21 filament and $\sim$ 0.1 mG in the surrounding ambient gas. The mass-to-flux ratios are found to be magnetically supercritical in the filament and slightly subcritical to nearly critical in the ambient gas. The alignment between column density structures and magnetic fields changes from random alignment in the low-density ambient gas probed by Planck to mostly perpendicular in the high-density main filament probed by JCMT. The magnetic field structures of the DR21 filament are in agreement with MHD simulations of a strongly magnetized medium, suggesting that magnetic fields play an important role in shaping the DR21 main filament and sub-filaments. △ Less

Submitted 4 December, 2022; originally announced December 2022.

Comments: 26 pages, 13 figures, ApJ accepted

WALLABY Pilot Survey: HI gas kinematics of galaxy pairs in cluster environment

Authors: Shin-Jeong Kim, Se-Heon Oh, Jing Wang, Lister Staveley-Smith, Bärbel S. Koribalski, Minsu Kim, Hye-Jin Park, Shinna Kim, Kristine Spekkens, Tobias Westmeier, O. Ivy Wong, Gerhardt R. Meurer, Peter Kamphuis., Barbara Catinella, Kristen B. W. McQuinn, Frank Bigiel, Benne W. Holwerda, Jonghwan Rhee, Karen Lee-Waddell, Nathan Deg, Lourdes Verdes-Montenegro, Bi-Qing For, Juan P. Madrid, Helga Dénes, Ahmed Elagali

Abstract: We examine the HI gas kinematics of galaxy pairs in two clusters and a group using Australian Square Kilometre Array Pathfinder (ASKAP) WALLABY pilot survey observations. We compare the HI properties of galaxy pair candidates in the Hydra I and Norma clusters, and the NGC 4636 group, with those of non-paired control galaxies selected in the same fields. We perform HI profile decomposition of the s… ▽ More We examine the HI gas kinematics of galaxy pairs in two clusters and a group using Australian Square Kilometre Array Pathfinder (ASKAP) WALLABY pilot survey observations. We compare the HI properties of galaxy pair candidates in the Hydra I and Norma clusters, and the NGC 4636 group, with those of non-paired control galaxies selected in the same fields. We perform HI profile decomposition of the sample galaxies using a tool, {\sc baygaud} which allows us to de-blend a line-of-sight velocity profile with an optimal number of Gaussian components. We construct HI super-profiles of the sample galaxies via stacking of their line profiles after aligning the central velocities. We fit a double Gaussian model to the super-profiles and classify them as kinematically narrow and broad components with respect to their velocity dispersions. Additionally, we investigate the gravitational instability of HI gas disks of the sample galaxies using Toomre Q parameters and HI morphological disturbances. We investigate the effect of the cluster environment on the HI properties of galaxy pairs by dividing the cluster environment into three subcluster regions (i.e., outskirts, infalling and central regions). We find that the denser cluster environment (i.e., infalling and central regions) is likely to impact the HI gas properties of galaxies in a way of decreasing the amplitude of the kinematically narrow HI gas ($M_{\rm{narrow}}^{\rm{HI}}$/$M_{\rm{total}}^{\rm{HI}}$), and increasing the Toomre Q values of the infalling and central galaxies. This tendency is likely to be more enhanced for galaxy pairs in the cluster environment. △ Less

Submitted 28 November, 2022; originally announced November 2022.

Comments: Accepted for publication in MNRAS, 25 pages, 20 figures. Supplementary material included (20 pages)

The JCMT BISTRO Survey: A Spiral Magnetic Field in a Hub-filament Structure, Monoceros R2

Authors: Jihye Hwang, Jongsoo Kim, Kate Pattle, Chang Won Lee, Patrick M. Koch, Doug Johnstone, Kohji Tomisaka, Anthony Whitworth, Ray S. Furuya, Ji-hyun Kang, A-Ran Lyo, Eun Jung Chung, Doris Arzoumanian, Geumsook Park, Woojin Kwon, Shinyoung Kim, Motohide Tamura, Jungmi Kwon, Archana Soam, Ilseung Han, Thiem Hoang, Kyoung Hee Kim, Takashi Onaka, Eswaraiah Chakali, Derek Ward-Thompson , et al. (135 additional authors not shown)

Abstract: We present and analyze observations of polarized dust emission at 850 $μ$m towards the central 1 pc $\times$ 1 pc hub-filament structure of Monoceros R2 (Mon R2). The data are obtained with SCUBA-2/POL-2 on the James Clerk Maxwell Telescope (JCMT) as part of the BISTRO (B-fields in Star-forming Region Observations) survey. The orientations of the magnetic field follow the spiral structure of Mon R… ▽ More We present and analyze observations of polarized dust emission at 850 $μ$m towards the central 1 pc $\times$ 1 pc hub-filament structure of Monoceros R2 (Mon R2). The data are obtained with SCUBA-2/POL-2 on the James Clerk Maxwell Telescope (JCMT) as part of the BISTRO (B-fields in Star-forming Region Observations) survey. The orientations of the magnetic field follow the spiral structure of Mon R2, which are well-described by an axisymmetric magnetic field model. We estimate the turbulent component of the magnetic field using the angle difference between our observations and the best-fit model of the underlying large-scale mean magnetic field. This estimate is used to calculate the magnetic field strength using the Davis-Chandrasekhar-Fermi method, for which we also obtain the distribution of volume density and velocity dispersion using a column density map derived from $Herschel$ data and the C$^{18}$O ($J$ = 3-2) data taken with HARP on the JCMT, respectively. We make maps of magnetic field strengths and mass-to-flux ratios, finding that magnetic field strengths vary from 0.02 to 3.64 mG with a mean value of 1.0 $\pm$ 0.06 mG, and the mean critical mass-to-flux ratio is 0.47 $\pm$ 0.02. Additionally, the mean Alfvén Mach number is 0.35 $\pm$ 0.01. This suggests that in Mon R2, magnetic fields provide resistance against large-scale gravitational collapse, and magnetic pressure exceeds turbulent pressure. We also investigate the properties of each filament in Mon R2. Most of the filaments are aligned along the magnetic field direction and are magnetically sub-critical. △ Less

Submitted 13 December, 2022; v1 submitted 12 October, 2022; originally announced October 2022.

Comments: This paper is accepted to the ApJ

KMTNet Nearby Galaxy Survey: Overview and a Survey Description

Authors: Woowon Byun, Yun-Kyeong Sheen, Kwang-Il Seon, Luis C. Ho, Joon Hyeop Lee, Hyunjin Jeong, Sang Chul Kim, Byeong-Gon Park, Yongseok Lee, Sang-Mok Cha, Minjin Kim

Abstract: Recently, there has been an increasing demand for deep imaging surveys to investigate the history of the mass assembly of galaxies in detail by examining the remnants of mergers and accretions, both of which have very low surface brightness (LSB). In addition, the nature of star formation in LSB regions, such as galaxy outer disks, is also an intriguing topic in terms of understanding the physical… ▽ More Recently, there has been an increasing demand for deep imaging surveys to investigate the history of the mass assembly of galaxies in detail by examining the remnants of mergers and accretions, both of which have very low surface brightness (LSB). In addition, the nature of star formation in LSB regions, such as galaxy outer disks, is also an intriguing topic in terms of understanding the physical mechanisms of disk evolution. To address these issues, this study conducts a survey project, called the Korea Microlensing Telescope Network (KMTNet) Nearby Galaxy Survey to construct a deep imaging data set of nearby galaxies in the southern hemisphere using KMTNet. It provides deep and wide-field images with a field-of-view of $\sim$12 deg$^2$ for 13 nearby galaxies drawn from the Carnegie-Irvine Galaxy Survey catalog, in optical broadbands ($BRI$) and an H$α$ narrowband. Through a dedicated data reduction, the surface brightness limit in 10$^{\prime\prime}\times10^{\prime\prime}$ boxes was found to reach as deep as $μ_{1σ}\sim29$-31 mag arcsec$^{-2}$ in the optical broadbands and $f_{1σ}\sim1$-$2\times 10^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ in the H$α$ narrowband. To conclude the paper, several possible scientific applications for this data set are described. △ Less

Submitted 2 September, 2022; originally announced September 2022.

Comments: 32 pages, 21 figures, accepted for publication in PASP

Global HI Properties of Galaxies via Super-profile Analysis

Authors: Minsu Kim, Se-Heon Oh

Abstract: We present a new method which constructs an HI super-profile of a galaxy which is based on profile decomposition analysis. The decomposed velocity profiles of an HI data cube with an optimal number of Gaussian components are co-added after being aligned in velocity with respect to their centroid velocities. This is compared to the previous approach where no prior profile decomposition is made for… ▽ More We present a new method which constructs an HI super-profile of a galaxy which is based on profile decomposition analysis. The decomposed velocity profiles of an HI data cube with an optimal number of Gaussian components are co-added after being aligned in velocity with respect to their centroid velocities. This is compared to the previous approach where no prior profile decomposition is made for the velocity profiles being stacked. The S/N improved super-profile is useful for deriving the galaxy's global HI properties like velocity dispersion and mass from observations which do not provide sufficient surface brightness sensitivity for the galaxy. As a practical test, we apply our new method to 64 high-resolution HI data cubes of nearby galaxies in the local Universe which are taken from THINGS and LITTLE THINGS. In addition, we also construct two additional HI super-profiles of the sample galaxies using symmetric and all velocity profiles of the cubes whose centroid velocities are determined from Hermite $h_3$ polynomial fitting, respectively. We find that the HI super-profiles constructed using the new method have narrower cores and broader wings in shape than the other two super-profiles. This is mainly due to the effect of either asymmetric velocity profiles' central velocity bias or the removal of asymmetric velocity profiles in the previous methods on the resulting HI super-profiles. We discuss how the shapes ($σ_{\rm{n}}/σ_{\rm{b}}$, $A_{\rm{n}}/A_{\rm{b}}$, and $A_{\rm{n}}/A_{\rm{tot}}$) of the new HI super-profiles which are measured from a double Gaussian fit are correlated with star formation rates of the sample galaxies and are compared with those of the other two super-profiles. △ Less

Submitted 4 November, 2022; v1 submitted 1 September, 2022; originally announced September 2022.

Comments: 24 pages, 13 figures, published in JKAS

Journal ref: JKAS Vol.55 No.5 pp.149-172 (2022)

A Mid-infrared Flare in the Seyfert Galaxy NGC 3786: A Changing-look Event Triggered by an Obscured Tidal Disruption Event?

Authors: Suyeon Son, Minjin Kim, Luis C. Ho, Dohyeong Kim, Taehyun Kim

Abstract: We report an exceptional mid-infrared flare in the Seyfert 1.8 NGC 3786. In the multi-epoch data from the Wide-field Infrared Survey Explorer, the nuclear mid-infrared brightness of NGC 3786 appears to vary substantially up to $0.5-0.8$ mag around mid-2020. However, there is no evidence of significant variation in the corresponding light curve of the optical band from the Zwicky Transient Facility… ▽ More We report an exceptional mid-infrared flare in the Seyfert 1.8 NGC 3786. In the multi-epoch data from the Wide-field Infrared Survey Explorer, the nuclear mid-infrared brightness of NGC 3786 appears to vary substantially up to $0.5-0.8$ mag around mid-2020. However, there is no evidence of significant variation in the corresponding light curve of the optical band from the Zwicky Transient Facility. This implies that the flare may have been heavily obscured by nuclear dust. Through follow-up spectroscopic observations with Gemini-North after the flare, we find that broad emission lines in ${\rm Pa}α$ and ${\rm Pa}β$ newly appear, while the broad ${\rm H}β$ emission is marginally detected in the post-flare spectrum. In addition, their central wavelengths are systematically redshifted up to 900 km s$^{-1}$ with respect to the narrow emission lines. This reveals that the flare is associated with the changing-look phenomenon from type 1.8 to type 1. Based on these findings, we argue that the flare is likely to originate from an obscured tidal disruption event, although extreme variation in the accretion rate may not be ruled out completely. △ Less

Submitted 17 August, 2022; originally announced August 2022.

Comments: Accepted to ApJ

Reflectivity of Venus' dayside disk during the 2020 observation campaign: outcomes and future perspectives

Authors: Yeon Joo Lee, Antonio García Muñoz, Atsushi Yamazaki, Eric Quémerais, Stefano Mottola, Stephan Hellmich, Thomas Granzer, Gilles Bergond, Martin Roth, Eulalia Gallego-Cano, Jean-Yves Chaufray, Rozenn Robidel, Go Murakami, Kei Masunaga, Murat Kaplan, Orhan Erece, Ricardo Hueso, Petr Kabáth, Magdaléna Špoková, Agustín Sánchez-Lavega, Myung-Jin Kim, Valeria Mangano, Kandis-Lea Jessup, Thomas Widemann, Ko-ichiro Sugiyama , et al. (6 additional authors not shown)

Abstract: We performed a unique Venus observation campaign to measure the disk brightness of Venus over a broad range of wavelengths in August and September 2020. The primary goal of the campaign is to investigate the absorption properties of the unknown absorber in the clouds. The secondary goal is to extract a disk mean SO$_2$ gas abundance, whose absorption spectral feature is entangled with that of the… ▽ More We performed a unique Venus observation campaign to measure the disk brightness of Venus over a broad range of wavelengths in August and September 2020. The primary goal of the campaign is to investigate the absorption properties of the unknown absorber in the clouds. The secondary goal is to extract a disk mean SO$_2$ gas abundance, whose absorption spectral feature is entangled with that of the unknown absorber at the ultraviolet (UV) wavelengths. A total of 3 spacecraft and 6 ground-based telescopes participated in this campaign, covering the 52 to 1700~nm wavelength range. After careful evaluation of the observational data, we focused on the data sets acquired by 4 facilities. We accomplished our primary goal by analyzing the reflectivity spectrum of the Venus disk over the 283-800 nm wavelengths. Considerable absorption is present in the 350-450 nm range, for which we retrieved the corresponding optical depth by the unknown absorber. The result shows a consistent wavelength dependence of the relative optical depth with that at low latitudes during the Venus flyby by MESSENGER in 2007 (Pérez-Hoyos et al. 2018), which was expected because the overall disk reflectivity is dominated by low latitudes. Last, we summarize the experience obtained during this first campaign that should enable us to accomplish our second goal in future campaigns. △ Less

Submitted 27 July, 2022; originally announced July 2022.

Comments: Accepted for publication in PSJ, 45 pages, 13 figures

The Infrared Medium-deep Survey. IX. Discovery of Two New $z\sim6$ Quasars and Space Density down to $M_{1450}\sim-23.5$ mag

Authors: Yongjung Kim, Myungshin Im, Yiseul Jeon, Minjin Kim, Linhua Jiang, Suhyun Shin, Changsu Choi, Minhee Hyun, Hyunsung D. Jun, Dohyeong Kim, Duho Kim, Jae-Woo Kim, Ji Hoon Kim, Bumhyun Lee, Seong-Kook Lee, Juan Molina, Soojong Pak, Won-Kee Park, Yoon Chan Taak, Yongmin Yoon

Abstract: We present the result of the Infrared Medium-deep Survey (IMS) $z\sim6$ quasar survey, using the combination of the IMS near-infrared images and the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) optical images. The traditional color-selection method results in 25 quasar candidates over $86$ deg$^{2}$. We introduce the corrected Akaike Information Criterion (AICc) with the high-redshift qua… ▽ More We present the result of the Infrared Medium-deep Survey (IMS) $z\sim6$ quasar survey, using the combination of the IMS near-infrared images and the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) optical images. The traditional color-selection method results in 25 quasar candidates over $86$ deg$^{2}$. We introduce the corrected Akaike Information Criterion (AICc) with the high-redshift quasar and late-type star models to prioritize the candidates efficiently. Among the color-selected candidates, seven plausible candidates finally passed the AICc selection of which three are known quasars at $z\sim6$. The follow-up spectroscopic observations for the remaining four candidates were carried out, and we confirmed that two out of four are $z\sim6$ quasars. With this complete sample, we revisited the quasar space density at $z\sim6$ down to $M_{1450}\sim-23.5$ mag. Our result supports the low quasar space density at the luminosity where the quasar's ultraviolet ionizing emissivity peaks, favoring a minor contribution of quasars to the cosmic reionization. △ Less

Submitted 27 July, 2022; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: 21 pages, 13 figures, 4 tables. Accepted for publication in AJ