Showing 1–50 of 235 results for author: Nakagawa, T

A measure of departure from symmetry via the Fisher-Rao distance for contingency tables

Authors: Wataru Urasaki, Go Kawamitsu, Tomoyuki Nakagawa, Kouji Tahata

Abstract: A measure of asymmetry is a quantification method that allows for the comparison of categorical evaluations before and after treatment effects or among different target populations, irrespective of sample size. We focus on square contingency tables that summarize survey results between two time points or cohorts, represented by the same categorical variables. We propose a measure to evaluate the d… ▽ More A measure of asymmetry is a quantification method that allows for the comparison of categorical evaluations before and after treatment effects or among different target populations, irrespective of sample size. We focus on square contingency tables that summarize survey results between two time points or cohorts, represented by the same categorical variables. We propose a measure to evaluate the degree of departure from a symmetry model using cosine similarity. This proposal is based on the Fisher-Rao distance, allowing asymmetry to be interpreted as a geodesic distance between two distributions. Various measures of asymmetry have been proposed, but visualizing the relationship of these quantification methods on a two-dimensional plane demonstrates that the proposed measure provides the geometrically simplest and most natural quantification. Moreover, the visualized figure indicates that the proposed method for measuring departures from symmetry is less affected by very few cells with extreme asymmetry. A simulation study shows that for square contingency tables with an underlying asymmetry model, our method can directly extract and quantify only the asymmetric structure of the model, and can more sensitively detect departures from symmetry than divergence-type measures. △ Less

Submitted 21 October, 2024; originally announced October 2024.

MSC Class: 62H17; 62H20

Extended Japanese Commonsense Morality Dataset with Masked Token and Label Enhancement

Authors: Takumi Ohashi, Tsubasa Nakagawa, Hitoshi Iyatomi

Abstract: Rapid advancements in artificial intelligence (AI) have made it crucial to integrate moral reasoning into AI systems. However, existing models and datasets often overlook regional and cultural differences. To address this shortcoming, we have expanded the JCommonsenseMorality (JCM) dataset, the only publicly available dataset focused on Japanese morality. The Extended JCM (eJCM) has grown from the… ▽ More Rapid advancements in artificial intelligence (AI) have made it crucial to integrate moral reasoning into AI systems. However, existing models and datasets often overlook regional and cultural differences. To address this shortcoming, we have expanded the JCommonsenseMorality (JCM) dataset, the only publicly available dataset focused on Japanese morality. The Extended JCM (eJCM) has grown from the original 13,975 sentences to 31,184 sentences using our proposed sentence expansion method called Masked Token and Label Enhancement (MTLE). MTLE selectively masks important parts of sentences related to moral judgment and replaces them with alternative expressions generated by a large language model (LLM), while re-assigning appropriate labels. The model trained using our eJCM achieved an F1 score of 0.857, higher than the scores for the original JCM (0.837), ChatGPT one-shot classification (0.841), and data augmented using AugGPT, a state-of-the-art augmentation method (0.850). Specifically, in complex moral reasoning tasks unique to Japanese culture, the model trained with eJCM showed a significant improvement in performance (increasing from 0.681 to 0.756) and achieved a performance close to that of GPT-4 Turbo (0.787). These results demonstrate the validity of the eJCM dataset and the importance of developing models and datasets that consider the cultural context. △ Less

Submitted 12 October, 2024; originally announced October 2024.

Journal ref: Proceedings of the 33rd ACM International Conference on Information and Knowledge Management (CIKM '24), October 21--25, 2024, Boise, ID, USA

Systematic Study of the Inner Structure of Molecular Tori in Nearby U/LIRGs using Velocity Decomposition of CO Rovibrational Absorption Lines

Authors: Shusuke Onishi, Takao Nakagawa, Shunsuke Baba, Kosei Matsumoto, Naoki Isobe, Mai Shirahata, Hiroshi Terada, Tomonori Usuda, Shinki Oyabu

Abstract: Determining the inner structure of the molecular torus around an active galactic nucleus is essential for understanding its formation mechanism. However, spatially resolving the torus is difficult because of its small size. To probe the clump conditions in the torus, we therefore perform the systematic velocity-decomposition analyses of the gaseous CO rovibrational absorption lines (… ▽ More Determining the inner structure of the molecular torus around an active galactic nucleus is essential for understanding its formation mechanism. However, spatially resolving the torus is difficult because of its small size. To probe the clump conditions in the torus, we therefore perform the systematic velocity-decomposition analyses of the gaseous CO rovibrational absorption lines ($v=0\to 1,ΔJ=\pm 1$) at $λ\sim 4.67 \mathrm{μ{m}}$ observed toward four (ultra)luminous infrared galaxies using the high-resolution ($R\sim 5000\text{--}10000$) spectroscopy from the Subaru Telescope. We find that each transition has two to five distinct velocity components with different line-of-sight (LOS) velocities ($V_\mathrm{LOS}\sim -240\text{--}+100\mathrm{km\,s^{-1}}$) and dispersions ($σ_V\sim 15\text{--}190\mathrm{km\,s^{-1}}$); i.e., the components (a), (b), ..., beginning with the broadest one in each target, indicating that the tori have clumpy structures. By assuming a hydrostatic disk ($σ_V\propto R_\mathrm{rot}^{-0.5}$), we find that the tori have dynamic inner structures, with the innermost component (a) outflowing with velocity $|V_\mathrm{LOS}|\sim 160\text{--}240\mathrm{km\,s^{-1}}$, and the outer components (b) and (c) outflowing more slowly or infalling with $|V_\mathrm{LOS}|\lesssim 100\mathrm{km\,s^{-1}}$. In addition, we find that the innermost component (a) can be attributed to collisionally excited hot ($\gtrsim 530$K) and dense ($n_\mathrm{H_2}\gtrsim 10^6\mathrm{cm^{-3}}$) clumps, based on the level populations. Conversely, the outer component (b) can be attributed to cold ($\sim 30\text{--}140$K) clumps radiatively excited by a far-infrared-to-submillimeter background with a brightness temperature higher than $\sim 20\text{--}400$K. These observational results demonstrate the clumpy and dynamic structure of tori in the presence of background radiation. △ Less

Submitted 9 October, 2024; originally announced October 2024.

Comments: 43 pages, 15 figures, accepted for publication in ApJ

Visualization for departures from symmetry with the power-divergence-type measure in two-way contingency tables

Authors: Wataru Urasaki, Tomoyuki Nakagawa, Jun Tsuchida, Kouji Tahata

Abstract: When the row and column variables consist of the same category in a two-way contingency table, it is specifically called a square contingency table. Since it is clear that the square contingency tables have an association structure, a primary objective is to examine symmetric relationships and transitions between variables. While various models and measures have been proposed to analyze these stru… ▽ More When the row and column variables consist of the same category in a two-way contingency table, it is specifically called a square contingency table. Since it is clear that the square contingency tables have an association structure, a primary objective is to examine symmetric relationships and transitions between variables. While various models and measures have been proposed to analyze these structures understanding changes between two variables in behavior at two-time points or cohorts, it is also necessary to require a detailed investigation of individual categories and their interrelationships, such as shifts in brand preferences. This paper proposes a novel approach to correspondence analysis (CA) for evaluating departures from symmetry in square contingency tables with nominal categories, using a power-divergence-type measure. The approach ensures that well-known divergences can also be visualized and, regardless of the divergence used, the CA plot consists of two principal axes with equal contribution rates. Additionally, the scaling is independent of sample size, making it well-suited for comparing departures from symmetry across multiple contingency tables. Confidence regions are also constructed to enhance the accuracy of the CA plot. △ Less

Submitted 30 September, 2024; originally announced October 2024.

Narrowing band gap chemically and physically: Conductive dense hydrocarbon

Authors: Takeshi Nakagawa, Caoshun Zhang, Kejun Bu, Philip Dalladay-Simpson, Martina Vrankić, Sarah Bolton, Dominique Laniel, Dong Wang, Akun Liang, Hirofumi Ishii, Nozomu Hiraoka, Gaston Garbarino, Angelika D. Rosa, Qingyang Hu, Xujie Lü, Ho-kwang Mao, Yang Ding

Abstract: Band gap energy of an organic molecule can be reduced by intermolecular interaction enhancement, and thus, certain polycyclic aromatic hydrocarbons (PAHs), which are insulators with wide band gaps, are expected to undergo insulator-metal transitions by simple compression. Such a pressure-induced electronic transition can be exploited to transform non-metallic organic materials into states featurin… ▽ More Band gap energy of an organic molecule can be reduced by intermolecular interaction enhancement, and thus, certain polycyclic aromatic hydrocarbons (PAHs), which are insulators with wide band gaps, are expected to undergo insulator-metal transitions by simple compression. Such a pressure-induced electronic transition can be exploited to transform non-metallic organic materials into states featuring intriguing electronic characteristics such as high-temperature superconductivity. Numerous attempts have been made to metalize various small PAHs, but so far only pressure-induced amorphization well below the megabar region was observed. The wide band gap energy of the small PAHs and low chemical stability under simple compression are the bottlenecks. We have investigated the band gap energy evolution and the crystal structural compression of the large PAH molecules, where the band gap energy is significantly reduced by increasing the number of π-electrons and improved chemical stability with fully benzenoid molecular structure. Herein, we present a pressure-induced transition in dicoronylene, C48H20, an insulator at ambient conditions that transforms into a semi-metallic state above 23.0 GPa with a three-order-of-magnitude reduction in resistivity. In-situ UV-visible absorption, transport property measurement, Raman spectroscopy, X-ray diffraction and density functional theory calculations were performed to provide tentative explanations to the alterations in its electronic structure at high pressure. The discovery of an electronic transition at pressures well below the megabar is a promising step towards realization of a single component purely hydrocarbon molecular metal in the near future. △ Less

Submitted 18 September, 2024; originally announced September 2024.

Comments: 22 pages, 5 figures

Diffusion Driven Transient Hydrogenation in Metal Superhydrides at Extreme Conditions

Authors: Yishan Zhou, Yunhua Fu, Meng Yang, Israel Osmond, Rajesh Jana, Takeshi Nakagawa, Owen Moulding, Jonathan Buhot, Sven Friedemann, Dominique Laniel, Thomas Meier

Abstract: In recent years, metal hydride research has become one of the driving forces of the high-pressure community, as it is believed to hold the key to superconductivity close to ambient temperature. While numerous novel metal hydride compounds have been reported and extensively investigated for their superconducting properties, little attention has been focused on the atomic and electronic states of hy… ▽ More In recent years, metal hydride research has become one of the driving forces of the high-pressure community, as it is believed to hold the key to superconductivity close to ambient temperature. While numerous novel metal hydride compounds have been reported and extensively investigated for their superconducting properties, little attention has been focused on the atomic and electronic states of hydrogen, the main ingredient in these novel compounds. Here, we present combined $^{1}H$- and $^{139}La$-NMR data on lanthanum superhydrides, $LaH_{x}$, ($x = 10.2 - 11.1$), synthesized after laser heating at pressures above 160 GPa. Strikingly, we found hydrogen to be in a highly diffusive state at room temperature, with diffusion coefficients in the order of $10^{-6}~cm^2s^{-1}$. We found that this diffusive state of hydrogen results in a dynamic de-hydrogenation of the sample over the course of several weeks, approaching a composition similar to its precursor materials. Quantitative measurements demonstrate that the synthesized superhydrides continuously decompose over time. Transport measurements underline this conclusion as superconducting critical temperatures were found to decrease significantly over time as well. This observation sheds new light on formerly unanswered questions on the long-term stability of metal superhydrides. △ Less

Submitted 23 August, 2024; originally announced August 2024.

Hexagonal to Monoclinic Phase Transition in Dense Hydrogen Phase III Detected by High-Pressure NMR

Authors: Meng Yang, Yishan Zhou, Rajesh Jana, Takeshi Nakagawa, Yunhua Fu, Thomas Meier

Abstract: Conclusive crystal structure determination of the high pressure phases of hydrogen remains elusive due to lack of core electrons and vanishing wave vectors, rendering standard high-pressure experimental methods moot. Ab-initio DFT calculations have shown that structural polymorphism might be solely resolvable using high-resolution nuclear magnetic resonance (NMR) spectroscopy at mega-bar pressures… ▽ More Conclusive crystal structure determination of the high pressure phases of hydrogen remains elusive due to lack of core electrons and vanishing wave vectors, rendering standard high-pressure experimental methods moot. Ab-initio DFT calculations have shown that structural polymorphism might be solely resolvable using high-resolution nuclear magnetic resonance (NMR) spectroscopy at mega-bar pressures, however technical challenges have precluded such experiments thus far. Here, we present in-situ high-pressure high-resolution NMR experiments in hydrogen phase III between 181 GPa and 208 GPa at room temperature. Our spectra suggest that at lower pressures phase III adopts a hexagonal P6122 crystal structure, transitioning into a monoclinic C2/c phase at about 197 GPa. The high resolution spectra are in excellent agreement with earlier structural and spectral predictions and underline the possibility of a subtle P6122 to C2/c phase transition in hydrogen phase III. These experiments show the importance of a combination of ab-initio calculations and low-Z sensitive spectral probes in high-pressure science in elucidating the structural complexity of the most abundant element in our universe. △ Less

Submitted 27 July, 2024; originally announced July 2024.

Robust VAEs via Generating Process of Noise Augmented Data

Authors: Hiroo Irobe, Wataru Aoki, Kimihiro Yamazaki, Yuhui Zhang, Takumi Nakagawa, Hiroki Waida, Yuichiro Wada, Takafumi Kanamori

Abstract: Advancing defensive mechanisms against adversarial attacks in generative models is a critical research topic in machine learning. Our study focuses on a specific type of generative models - Variational Auto-Encoders (VAEs). Contrary to common beliefs and existing literature which suggest that noise injection towards training data can make models more robust, our preliminary experiments revealed th… ▽ More Advancing defensive mechanisms against adversarial attacks in generative models is a critical research topic in machine learning. Our study focuses on a specific type of generative models - Variational Auto-Encoders (VAEs). Contrary to common beliefs and existing literature which suggest that noise injection towards training data can make models more robust, our preliminary experiments revealed that naive usage of noise augmentation technique did not substantially improve VAE robustness. In fact, it even degraded the quality of learned representations, making VAEs more susceptible to adversarial perturbations. This paper introduces a novel framework that enhances robustness by regularizing the latent space divergence between original and noise-augmented data. Through incorporating a paired probabilistic prior into the standard variational lower bound, our method significantly boosts defense against adversarial attacks. Our empirical evaluations demonstrate that this approach, termed Robust Augmented Variational Auto-ENcoder (RAVEN), yields superior performance in resisting adversarial inputs on widely-recognized benchmark datasets. △ Less

Submitted 26 July, 2024; originally announced July 2024.

General Bayesian inference for causal effects using covariate balancing procedure

Authors: Shunichiro Orihara, Tomotaka Momozaki, Tomoyuki Nakagawa

Abstract: In observational studies, the propensity score plays a central role in estimating causal effects of interest. The inverse probability weighting (IPW) estimator is commonly used for this purpose. However, if the propensity score model is misspecified, the IPW estimator may produce biased estimates of causal effects. Previous studies have proposed some robust propensity score estimation procedures.… ▽ More In observational studies, the propensity score plays a central role in estimating causal effects of interest. The inverse probability weighting (IPW) estimator is commonly used for this purpose. However, if the propensity score model is misspecified, the IPW estimator may produce biased estimates of causal effects. Previous studies have proposed some robust propensity score estimation procedures. However, these methods require considering parameters that dominate the uncertainty of sampling and treatment allocation. This study proposes a novel Bayesian estimating procedure that necessitates probabilistically deciding the parameter, rather than deterministically. Since the IPW estimator and propensity score estimator can be derived as solutions to certain loss functions, the general Bayesian paradigm, which does not require the considering the full likelihood, can be applied. Therefore, our proposed method only requires the same level of assumptions as ordinary causal inference contexts. The proposed Bayesian method demonstrates equal or superior results compared to some previous methods in simulation experimentss, and is also applied to real data, namely the Whitehall dataset. △ Less

Submitted 30 September, 2024; v1 submitted 14 April, 2024; originally announced April 2024.

Comments: covariate balancing, general Bayes, inverse probability weighting, M-estimator, propensity score

Balmer Decrement Anomalies in Galaxies at z ~ 6 Found by JWST Observations: Density-Bounded Nebulae or Excited H I Clouds?

Authors: Hiroto Yanagisawa, Masami Ouchi, Kimihiko Nakajima, Hidenobu Yajima, Hiroya Umeda, Shunsuke Baba, Takao Nakagawa, Minami Nakane, Akinori Matsumoto, Yoshiaki Ono, Yuichi Harikane, Yuki Isobe, Yi Xu, Yechi Zhang

Abstract: We investigate the physical origins of the Balmer decrement anomalies in GS-NDG-9422 (Cameron et al. 2023) and RXCJ2248-ID (Topping et al. 2024) galaxies at $z\sim 6$ whose $\mathrm{H}α/\mathrm{H}β$ values are significantly smaller than $2.7$, the latter of which also shows anomalous $\mathrm{H}γ/\mathrm{H}β$ and $\mathrm{H}δ/\mathrm{H}β$ values beyond the errors. Because the anomalous Balmer decr… ▽ More We investigate the physical origins of the Balmer decrement anomalies in GS-NDG-9422 (Cameron et al. 2023) and RXCJ2248-ID (Topping et al. 2024) galaxies at $z\sim 6$ whose $\mathrm{H}α/\mathrm{H}β$ values are significantly smaller than $2.7$, the latter of which also shows anomalous $\mathrm{H}γ/\mathrm{H}β$ and $\mathrm{H}δ/\mathrm{H}β$ values beyond the errors. Because the anomalous Balmer decrements are not reproduced under the Case B recombination, we explore the nebulae with the optical depths smaller and larger than the Case B recombination by physical modeling. We find two cases quantitatively explaining the anomalies; 1) density-bounded nebulae that are opaque only up to around Ly$γ$-Ly8 transitions and 2) ionization-bounded nebulae partly/fully surrounded by optically-thick excited H{\sc i} clouds. The case of 1) produces more H$β$ photons via Ly$γ$ absorption in the nebulae, requiring fine tuning in optical depth values, while this case helps ionizing photon escape for cosmic reionization. The case of 2) needs the optically-thick excited H{\sc i} clouds with $N_2\simeq 10^{12}-10^{13}$ $\mathrm{cm^{-2}}$, where $N_2$ is the column density of the hydrogen atom with the principal quantum number of $n=2$. Interestingly, the high $N_2$ values qualitatively agree with the recent claims for GS-NDG-9422 with the strong nebular continuum requiring a number of $2s$-state electrons and for RXCJ2248-ID with the dense ionized regions likely coexisting with the optically-thick clouds. While the physical origin of the optically-thick excited H{\sc i} clouds is unclear, these results may suggest gas clouds with excessive collisional excitation caused by an amount of accretion and supernovae in the high-$z$ galaxies. △ Less

Submitted 25 August, 2024; v1 submitted 29 March, 2024; originally announced March 2024.

Comments: 13 pages, 7 figures, accepted for publication in ApJ

Effects of galaxy environment on merger fraction

Authors: W. J. Pearson, D. J. D. Santos, T. Goto, T. -C. Huang, S. J. Kim, H. Matsuhara, A. Pollo, S. C. -C. Ho, H. S. Hwang, K. Małek, T. Nakagawa, M. Romano, S. Serjeant, L. Suelves, H. Shim, G. J. White

Abstract: Aims. In this work, we intend to examine how environment influences the merger fraction, from the low density field environment to higher density groups and clusters. We also aim to study how the properties of a group or cluster, as well as the position of a galaxy in the group or cluster, influences the merger fraction. Methods. We identified galaxy groups and clusters in the North Ecliptic Pol… ▽ More Aims. In this work, we intend to examine how environment influences the merger fraction, from the low density field environment to higher density groups and clusters. We also aim to study how the properties of a group or cluster, as well as the position of a galaxy in the group or cluster, influences the merger fraction. Methods. We identified galaxy groups and clusters in the North Ecliptic Pole using a friends-of-friends algorithm and the local density. Once identified, we determined the central galaxies, group radii, velocity dispersions, and group masses of these groups and clusters. Merging systems were identified with a neural network as well as visually. With these, we examined how the merger fraction changes as the local density changes for all galaxies as well as how the merger fraction changes as the properties of the groups or clusters change. Results. We find that the merger fraction increases as local density increases and decreases as the velocity dispersion increases, as is often found in literature. A decrease in merger fraction as the group mass increases is also found. We also find groups with larger radii have higher merger fractions. The number of galaxies in a group does not influence the merger fraction. Conclusions. The decrease in merger fraction as group mass increases is a result of the link between group mass and velocity dispersion. Hence, this decrease of merger fraction with increasing mass is a result of the decrease of merger fraction with velocity dispersion. The increasing relation between group radii and merger fraction may be a result of larger groups having smaller velocity dispersion at a larger distance from the centre or larger groups hosting smaller, infalling groups with more mergers. However, we do not find evidence of smaller groups having higher merger fractions. △ Less

Submitted 18 March, 2024; originally announced March 2024.

Comments: 15 pages, 10 figures, 8 tables, 2 appendices, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 686, A94 (2024)

Spectroscopic Constraints on the Mid-Infrared Attenuation Curve: I -- Attenuation Model using PAH Emission

Authors: Thomas S. -Y. Lai, J. D. T. Smith, Els Peeters, Henrik W. W. Spoon, Shunsuke Baba, Masatoshi Imanishi, Takao Nakagawa

Abstract: We introduce a novel model to spectroscopically constrain the mid-infrared (MIR) extinction/attenuation curve from 3--17 um, using Polycyclic Aromatic Hydrocarbon (PAH) emission drawn from an AKARI-Spitzer extragalactic cross-archival dataset. Currently proposed MIR extinction curves vary significantly in their slopes toward the near-infrared, and the variation of the strengths and shapes of the 9… ▽ More We introduce a novel model to spectroscopically constrain the mid-infrared (MIR) extinction/attenuation curve from 3--17 um, using Polycyclic Aromatic Hydrocarbon (PAH) emission drawn from an AKARI-Spitzer extragalactic cross-archival dataset. Currently proposed MIR extinction curves vary significantly in their slopes toward the near-infrared, and the variation of the strengths and shapes of the 9.7 um and 18 um silicate absorption features make MIR spectral modeling and interpretation challenging, particularly for heavily obscured galaxies. By adopting the basic premise that PAH bands have relatively consistent intrinsic ratios within dusty starbursting galaxies, we can, for the first time, empirically determine the overall shape of the MIR attenuation curve by measuring the differential attenuation at specific PAH wavelengths. Our attenuation model shows PAH emission in most (U)LIRGs is unambiguously subjected to attenuation, and we find strong evidence that PAH bands undergo differential attenuation as obscuration increases. Compared to pre-existing results, the MIR attenuation curve derived from the model favors relatively gray continuum absorption from 3-8 $μ$m and silicate features with intermediate strength at 9.7 um but with stronger than typical 18 um opacity. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: 15 pages, 5 figures. Paper re-submitted to ApJ after addressing the first round of comments by the referee

Direct excitation of Kelvin waves on quantized vortices

Authors: Yosuke Minowa, Yuki Yasui, Tomo Nakagawa, Sosuke Inui, Makoto Tsubota, Masaaki Ashida

Abstract: Helices and spirals, prevalent across various systems, play a crucial role in characterizing symmetry, describing dynamics, and imparting unique functionalities, attributed to their inherent simplicity and chiral nature. A helical excitation on a quantized vortex, an example of a one-dimensional topological defect, emerges as a Nambu-Goldstone mode following spontaneous symmetry breaking, known as… ▽ More Helices and spirals, prevalent across various systems, play a crucial role in characterizing symmetry, describing dynamics, and imparting unique functionalities, attributed to their inherent simplicity and chiral nature. A helical excitation on a quantized vortex, an example of a one-dimensional topological defect, emerges as a Nambu-Goldstone mode following spontaneous symmetry breaking, known as a Kelvin wave. Kelvin waves play a vital role in energy dissipation within inviscid quantum fluids. However, deliberately exciting Kelvin waves has proven to be challenging. Here, we introduce a controlled method for exciting Kelvin waves on a quantized vortex in superfluid helium-4. We used a charged nanoparticle, oscillated by a time-varying electric field, to stimulate Kelvin waves on the vortex. A major breakthrough in our research is the confirmation of the helical nature of Kelvin waves through three-dimensional image reconstruction, providing visual evidence of their complex dynamics. Additionally, we determined the dispersion relation and the phase velocity of the Kelvin wave and identified the vorticity direction, enhancing our understanding of quantum fluid behavior. This work elucidates the dynamics of Kelvin waves and pioneers a novel approach for manipulating and observing quantized vortices in three dimensions, thereby opening new avenues for exploring quantum fluidic systems. △ Less

Submitted 26 February, 2024; originally announced February 2024.

Comments: 18 pages, 4 figures

Unveiling a Novel Metal-to-Metal Transition in LuH2: Critically Challenging Superconductivity Claims in Lutetium Hydrides

Authors: Dong Wang, Ningning Wang, Caoshun Zhang, Chunsheng Xia, Weicheng Guo, Xia Yin, Kejun Bu, Takeshi Nakagawa, Jianbo Zhang, Federico Gorelli, Philip Dalladay-Simpson, Thomas Meier, Xujie Lü, Liling Sun, Jinguang Cheng, Qiaoshi Zeng, Yang Ding, Ho-kwang Mao

Abstract: Following the recent report by Dasenbrock-Gammon et al. (2023) of near-ambient superconductivity in nitrogen-doped lutetium trihydride (LuH3-δNε), significant debate has emerged surrounding the composition and interpretation of the observed sharp resistance drop. Here, we meticulously revisit these claims through comprehensive characterization and investigations. We definitively identify the repor… ▽ More Following the recent report by Dasenbrock-Gammon et al. (2023) of near-ambient superconductivity in nitrogen-doped lutetium trihydride (LuH3-δNε), significant debate has emerged surrounding the composition and interpretation of the observed sharp resistance drop. Here, we meticulously revisit these claims through comprehensive characterization and investigations. We definitively identify the reported material as lutetium dihydride (LuH2), resolving the ambiguity surrounding its composition. Under similar conditions (270-295 K and 1-2 GPa), we replicate the reported sharp decrease in electrical resistance with a 30% success rate, aligning with Dasenbrock-Gammon et al.'s observations. However, our extensive investigations reveal this phenomenon to be a novel, pressure-induced metal-to-metal transition intrinsic to LuH2, distinct from superconductivity. Intriguingly, nitrogen doping exerts minimal impact on this transition. Our work not only elucidates the fundamental properties of LuH2 and LuH3 but also critically challenges the notion of superconductivity in these lutetium hydride systems. These findings pave the way for future research on lutetium hydride systems while emphasizing the crucial importance of rigorous verification in claims of ambient temperature superconductivity. △ Less

Submitted 28 January, 2024; v1 submitted 25 January, 2024; originally announced January 2024.

Journal ref: Matter Radiat. Extremes 9, 037401 (2024)

Semiparametric Copula Estimation for Spatially Correlated Multivariate Mixed Outcomes: Analyzing Visual Sightings of Fin Whales from Line Transect Survey

Authors: Tomotaka Momozaki, Tomoyuki Nakagawa, Shonosuke Sugasawa, Hiroko Kato Solvang

Abstract: Multivariate data having both continuous and discrete variables is known as mixed outcomes and has widely appeared in a variety of fields such as ecology, epidemiology, and climatology. In order to understand the probability structure of multivariate data, the estimation of the dependence structure among mixed outcomes is very important. However, when location information is equipped with multivar… ▽ More Multivariate data having both continuous and discrete variables is known as mixed outcomes and has widely appeared in a variety of fields such as ecology, epidemiology, and climatology. In order to understand the probability structure of multivariate data, the estimation of the dependence structure among mixed outcomes is very important. However, when location information is equipped with multivariate data, the spatial correlation should be adequately taken into account; otherwise, the estimation of the dependence structure would be severely biased. To solve this issue, we propose a semiparametric Bayesian inference for the dependence structure among mixed outcomes while eliminating spatial correlation. To this end, we consider a hierarchical spatial model based on the rank likelihood and a latent multivariate Gaussian process. We develop an efficient algorithm for computing the posterior using the Markov Chain Monte Carlo. We also provide a scalable implementation of the model using the nearest-neighbor Gaussian process under large spatial datasets. We conduct a simulation study to validate our proposed procedure and demonstrate that the procedure successfully accounts for spatial correlation and correctly infers the dependence structure among outcomes. Furthermore, the procedure is applied to a real example collected during an international synoptic krill survey in the Scotia Sea of the Antarctic Peninsula, which includes sighting data of fin whales (Balaenoptera physalus), and the relevant oceanographic data. △ Less

Submitted 19 December, 2023; originally announced December 2023.

Comments: 23 pages, 5 figures

MSC Class: 62F15 (Primary); 62H11 (Secondary)

Nonparametric Bayesian Adjustment of Unmeasured Confounders in Cox Proportional Hazards Models

Authors: Shunichiro Orihara, Shonosuke Sugasawa, Tomohiro Ohigashi, Tomoyuki Nakagawa, Masataka Taguri

Abstract: In observational studies, unmeasured confounders present a crucial challenge in accurately estimating desired causal effects. To calculate the hazard ratio (HR) in Cox proportional hazard models for time-to-event outcomes, two-stage residual inclusion and limited information maximum likelihood are typically employed. However, these methods are known to entail difficulty in terms of potential bias… ▽ More In observational studies, unmeasured confounders present a crucial challenge in accurately estimating desired causal effects. To calculate the hazard ratio (HR) in Cox proportional hazard models for time-to-event outcomes, two-stage residual inclusion and limited information maximum likelihood are typically employed. However, these methods are known to entail difficulty in terms of potential bias of HR estimates and parameter identification. This study introduces a novel nonparametric Bayesian method designed to estimate an unbiased HR, addressing concerns that previous research methods have had. Our proposed method consists of two phases: 1) detecting clusters based on the likelihood of the exposure and outcome variables, and 2) estimating the hazard ratio within each cluster. Although it is implicitly assumed that unmeasured confounders affect outcomes through cluster effects, our algorithm is well-suited for such data structures. The proposed Bayesian estimator has good performance compared with some competitors. △ Less

Submitted 1 October, 2024; v1 submitted 4 December, 2023; originally announced December 2023.

Comments: general Bayes, invalid instrumental variable, Mendelian randomization, UK Biobank, weak instrumental variable

Structure responsible for the superconducting state in La3Ni2O7 at high pressure and low temperature conditions

Authors: Luhong Wang, Yan Li, Shengyi Xie, Fuyang Liu, Hualei Sun, Caoxin Huang, Yang Gao, Takeshi Nakagawa, Boyang Fu, Bo Dong, Zhenhui Cao, Runze Yu, Saori I. Kawaguchi, Hirokazu Kadobayashi, Meng Wang, Changqing Jin, Ho-kwang Mao, Haozhe Liu

Abstract: Very recently, a new superconductor with Tc = 80 K was reported in nickelate (La3Ni2O7) at around 15 - 40 GPa conditions (Nature, 621, 493, 2023) [1], which is the second type of unconventional superconductor, beside the cuprates, with Tc above liquid nitrogen temperature. However, the phase diagram plotted in this report was mostly based on the transport measurement at low temperature and high pr… ▽ More Very recently, a new superconductor with Tc = 80 K was reported in nickelate (La3Ni2O7) at around 15 - 40 GPa conditions (Nature, 621, 493, 2023) [1], which is the second type of unconventional superconductor, beside the cuprates, with Tc above liquid nitrogen temperature. However, the phase diagram plotted in this report was mostly based on the transport measurement at low temperature and high pressure conditions, and the assumed corresponding X-ray diffraction (XRD) results was carried out at room temperature. This encouraged us to carry out in situ high pressure and low temperature synchrotron XRD experiments to determine which phase is responsible for the high Tc state. In addition to the phase transition from orthorhombic Amam structure to orthorhombic Fmmm structure, a tetragonal phase with space group of I4/mmm was discovered when the sample was compressed to 19 GPa at 40 K where the superconductivity takes palce in La3Ni2O7. The calculations based on this tetragonal structure reveal that the electronic states approached to the Fermi energy were mainly dominated by the eg orbitals (3dz2 and 3dx2-y2) of Ni atoms, which are located in the oxygen octahedral crystal field. The correlation between Tc and this structural evolution, especially Ni-O octahedra regularity and the in-plane Ni-O-Ni bonding angles, are analyzed. This work sheds new lights to identify what is the most likely phase responsible for superconductivity in the double layered nickelate. △ Less

Submitted 21 November, 2023; v1 submitted 15 November, 2023; originally announced November 2023.

Self-consistent dust and non-LTE line radiative transfer with SKIRT

Authors: Kosei Matsumoto, Peter Camps, Maarten Baes, Frederik De Ceuster, Keiichi Wada, Takao Nakagawa, Kentaro Nagamine

Abstract: We introduce Monte Carlo-based non-LTE line radiative transfer calculations in the 3D dust radiative transfer code SKIRT, which was originally set up as a dust radiative transfer code. By doing so, we develop a generic and powerful 3D radiative transfer code that can self-consistently generate spectra with molecular and atomic lines against the underlying continuum. We test the accuracy of the non… ▽ More We introduce Monte Carlo-based non-LTE line radiative transfer calculations in the 3D dust radiative transfer code SKIRT, which was originally set up as a dust radiative transfer code. By doing so, we develop a generic and powerful 3D radiative transfer code that can self-consistently generate spectra with molecular and atomic lines against the underlying continuum. We test the accuracy of the non-LTE line radiative transfer module in the extended SKIRT code using standard benchmarks. We find excellent agreement between the SKIRT results, the published benchmark results, and results obtained using the ray-tracing non-LTE line radiative transfer code MAGRITTE, which validates our implementation. We apply the extended SKIRT code on a 3D hydrodynamic simulation of a dusty AGN torus model and generate multi-wavelength images with CO rotational-line spectra against the underlying dust continuum. We find that the low-J CO emission traces the geometrically thick molecular torus, whereas the higher-J CO lines originate from the gas with high kinetic temperature located in the innermost regions of the torus. Comparing the calculations with and without dust radiative transfer, we find that higher-J CO lines are slightly attenuated by the surrounding cold dust when seen edge-on. This shows that atomic and molecular lines can experience attenuation, an effect that is particularly important for transitions at mid- and near-infrared wavelengths. Therefore, our self-consistent dust and non-LTE line radiative transfer calculations can help interpret the observational data from Herschel, ALMA, and JWST. △ Less

Submitted 5 September, 2023; originally announced September 2023.

Comments: 12 pages, 9 figures, Accepted for publication in A&A

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

Dynamics of pinned quantized vortices in superfluid $^4$He in a microelectromechanical oscillator

Authors: Tomo Nakagawa, Makoto Tsubota, Keegan Gunther, Yoonseok Lee

Abstract: We numerically studied the vortex dynamics at zero temperature in superfluid $^4$He confined between two parallel rough solid boundaries, one of which oscillates in a shear mode. This study was motivated by the experimental work by Barquist $et$ $al.$ which employed a microelectromechanical systems (MEMS) oscillator operating in superfluid $^4$He at a near-zero temperature. Their experiments sugge… ▽ More We numerically studied the vortex dynamics at zero temperature in superfluid $^4$He confined between two parallel rough solid boundaries, one of which oscillates in a shear mode. This study was motivated by the experimental work by Barquist $et$ $al.$ which employed a microelectromechanical systems (MEMS) oscillator operating in superfluid $^4$He at a near-zero temperature. Their experiments suggest that the motion of the MEMS oscillator is damped by quantized vortices. In our study, we postulated that this damping effect was closely associated with vortex pinning phenomena and developed pinning models. Our primary objective is to understand the vortex dynamics in the presence of pinning and to provide insight into the experimental observations regarding the damping mechanism. We confirmed that Kelvin waves were excited in the pinned vortices when the oscillation frequency of the solid boundary matched with the mode frequency of the Kelvin wave. Additionally, we examined the formation and evolution of vortex tangles between the boundaries. The vortex tangle was suppressed in the presence of pinning, while the absence of pinning allowed to form well developed vortex tangle resulting in turbulence. Finally, by evaluating the tension of pinned vortices we extracted the damping force acting on the solid boundaries. △ Less

Submitted 25 October, 2023; v1 submitted 23 August, 2023; originally announced August 2023.

Comments: 11 pages + supplement, 17 figures

Journal ref: Phys. Rev. B 108, 144110 (2023)

Geometric Mean Type of Proportional Reduction in Variation Measure for Two-Way Contingency Tables

Authors: Wataru Urasaki, Yuki Wada, Tomoyuki Nakagawa, Kouji Tahata, Sadao Tomizawa

Abstract: In a two-way contingency table analysis with explanatory and response variables, the analyst is interested in the independence of the two variables. However, if the test of independence does not show independence or clearly shows a relationship, the analyst is interested in the degree of their association. Various measures have been proposed to calculate the degree of their association, one of whi… ▽ More In a two-way contingency table analysis with explanatory and response variables, the analyst is interested in the independence of the two variables. However, if the test of independence does not show independence or clearly shows a relationship, the analyst is interested in the degree of their association. Various measures have been proposed to calculate the degree of their association, one of which is the proportional reduction in variation (PRV) measure which describes the PRV from the marginal distribution to the conditional distribution of the response. The conventional PRV measures can assess the association of the entire contingency table, but they can not accurately assess the association for each explanatory variable. In this paper, we propose a geometric mean type of PRV (geoPRV) measure that aims to sensitively capture the association of each explanatory variable to the response variable by using a geometric mean, and it enables analysis without underestimation when there is partial bias in cells of the contingency table. Furthermore, the geoPRV measure is constructed by using any functions that satisfy specific conditions, which has application advantages and makes it possible to express conventional PRV measures as geometric mean types in special cases. △ Less

Submitted 5 July, 2023; originally announced July 2023.

Comments: 16 pages

ALMA ACA detection of submillimeter emission associated with the west hot spot of the radio galaxy Pictor A

Authors: Naoki Isobe, Hiroshi Nagai, Motoki Kino, Shunsuke Baba, Takao Nakagawa, Yuji Sunada, Makoto Tashiro

Abstract: In order to investigate the far-infrared excess detected from the west hot spot of the radio galaxy Pictor A with the Herschel observatory, a submillimeter photometry is performed with the Atacama Compact Array (ACA) of the Atacama Large Millimeter/submillimeter Array at Band 8 with the reference frequency of 405 GHz. A submillimeter source is discovered at the radio peak of the hot spot. Because… ▽ More In order to investigate the far-infrared excess detected from the west hot spot of the radio galaxy Pictor A with the Herschel observatory, a submillimeter photometry is performed with the Atacama Compact Array (ACA) of the Atacama Large Millimeter/submillimeter Array at Band 8 with the reference frequency of 405 GHz. A submillimeter source is discovered at the radio peak of the hot spot. Because the 405 GHz flux density of the source, $80.7\pm3.1$ mJy, agrees with the extrapolation of the synchrotron radio spectrum, the far-infrared excess is suggested to exhibit no major contribution at the ACA band. In contrast, by subtracting the power-law spectrum tightly constrained by the radio and ACA data, the significance of the excess in the Herschel band is well confirmed. No diffuse submillimeter emission is detected within the ACA field of view, and thus, the excess is ascribed to the west hot spot itself. In comparison to the previous estimate based on the Herschel data, the relative contribution of the far-infrared excess is reduced by a factor of $\sim 1.5$. The spectrum of the excess below the far-infrared band is determined to be harder than that of the diffusive shock acceleration. This strengthens the previous interpretation that the excess originates via the magnetic turbulence in the substructures within the hot spot. The ACA data are utilized to evaluate the magnetic field strength of the excess and of diffuse radio structure associated to the hot spot. △ Less

Submitted 26 June, 2023; originally announced June 2023.

Comments: 10 pages,4 figures, 3 tables, accepted for ApJ

Existence of density functions for SDEs driven by pure-jump processes

Authors: Takuya Nakagawa, Ryoichi Suzuki

Abstract: We verify the existence of density functions of the running maximum of a stochastic differential equation (SDE) driven by a Brownian motion and a non-truncated stable process. This is proved by the existence of density functions of the running maximum of Wiener-Poisson functionals resulting from Bismut's approach to Malliavin calculus for jump processes. We verify the existence of density functions of the running maximum of a stochastic differential equation (SDE) driven by a Brownian motion and a non-truncated stable process. This is proved by the existence of density functions of the running maximum of Wiener-Poisson functionals resulting from Bismut's approach to Malliavin calculus for jump processes. △ Less

Submitted 5 June, 2023; originally announced June 2023.

Comments: 21 pages

MSC Class: 60H10 60G52 60H07

Warm Molecular Gas in the Central Parsecs of the Buried Nucleus of NGC 4418 Traced with the Fundamental CO Ro-vibrational Absorptions

Authors: Youichi Ohyama, Shusuke Onishi, Takao Nakagawa, Kosei Matsumoto, Naoki Isobe, Mai Shirahata, Shunsuke Baba, Kazushi Sakamoto

Abstract: We investigated the inner buried nucleus of a nearby luminous infrared galaxy NGC 4418 using high-resolution spectroscopy of fundamental carbon monoxide (CO) ro-vibrational absorptions around $4.67 μ$m for the first time. This method allowed us to examine the physical and kinematical properties in the hot inner region of this nucleus. We detected a series of both very deep (partly saturated)… ▽ More We investigated the inner buried nucleus of a nearby luminous infrared galaxy NGC 4418 using high-resolution spectroscopy of fundamental carbon monoxide (CO) ro-vibrational absorptions around $4.67 μ$m for the first time. This method allowed us to examine the physical and kinematical properties in the hot inner region of this nucleus. We detected a series of both very deep (partly saturated) $^{12}$CO and moderately deep (optically thin) $^{13}$CO absorption lines and inferred a large column density ($N_\mathrm{H2}=(5\pm3)\times10^{23}$ cm$^{-2}$ in front of the $5 μ$m photosphere) of warm ($T_\mathrm{ex}\simeq170$ K) molecular gas by assuming an isothermal plane-parallel slab illuminated by a compact background MIR-emitting source. We modeled that the warm CO absorber almost covers the central heating source and that it is an inner layer around the $5 μ$m photosphere (at $r=$several pc) of a compact shroud of gas and dust ($d\sim100$ pc). The width of the absorption lines ($110$ km s$^{-1}$) and their small deviation from the systemic velocity ($<10$ km s$^{-1}$) are consistent with a warm and turbulent layer with little bulk motion in the radial direction. △ Less

Submitted 31 May, 2023; v1 submitted 17 May, 2023; originally announced May 2023.

Comments: 18 pages, 3 figures, 2 tables, accepted for publication in ApJ

Robustness of Bayesian ordinal response model against outliers via divergence approach

Authors: Tomotaka Momozaki, Tomoyuki Nakagawa

Abstract: Ordinal response model is a popular and commonly used regression for ordered categorical data in a wide range of fields such as medicine and social sciences. However, it is empirically known that the existence of ``outliers'', combinations of the ordered categorical response and covariates that are heterogeneous compared to other pairs, makes the inference with the ordinal response model unreliabl… ▽ More Ordinal response model is a popular and commonly used regression for ordered categorical data in a wide range of fields such as medicine and social sciences. However, it is empirically known that the existence of ``outliers'', combinations of the ordered categorical response and covariates that are heterogeneous compared to other pairs, makes the inference with the ordinal response model unreliable. In this article, we prove that the posterior distribution in the ordinal response model does not satisfy the posterior robustness with any link functions, i.e., the posterior cannot ignore the influence of large outliers. Furthermore, to achieve robust Bayesian inference in the ordinal response model, this article defines general posteriors in the ordinal response model with two robust divergences (the density-power and $γ$-divergences) based on the framework of the general posterior inference. We also provide an algorithm for generating posterior samples from the proposed posteriors. The robustness of the proposed methods against outliers is clarified from the posterior robustness and the index of robustness based on the Fisher-Rao metric. Through numerical experiments on artificial data and two real datasets, we show that the proposed methods perform better than the ordinary bayesian methods with and without outliers in the data for various link functions. △ Less

Submitted 12 May, 2023; originally announced May 2023.

Comments: 29 pages

Denoising Cosine Similarity: A Theory-Driven Approach for Efficient Representation Learning

Authors: Takumi Nakagawa, Yutaro Sanada, Hiroki Waida, Yuhui Zhang, Yuichiro Wada, Kōsaku Takanashi, Tomonori Yamada, Takafumi Kanamori

Abstract: Representation learning has been increasing its impact on the research and practice of machine learning, since it enables to learn representations that can apply to various downstream tasks efficiently. However, recent works pay little attention to the fact that real-world datasets used during the stage of representation learning are commonly contaminated by noise, which can degrade the quality of… ▽ More Representation learning has been increasing its impact on the research and practice of machine learning, since it enables to learn representations that can apply to various downstream tasks efficiently. However, recent works pay little attention to the fact that real-world datasets used during the stage of representation learning are commonly contaminated by noise, which can degrade the quality of learned representations. This paper tackles the problem to learn robust representations against noise in a raw dataset. To this end, inspired by recent works on denoising and the success of the cosine-similarity-based objective functions in representation learning, we propose the denoising Cosine-Similarity (dCS) loss. The dCS loss is a modified cosine-similarity loss and incorporates a denoising property, which is supported by both our theoretical and empirical findings. To make the dCS loss implementable, we also construct the estimators of the dCS loss with statistical guarantees. Finally, we empirically show the efficiency of the dCS loss over the baseline objective functions in vision and speech domains. △ Less

Submitted 19 April, 2023; originally announced April 2023.

Towards Understanding the Mechanism of Contrastive Learning via Similarity Structure: A Theoretical Analysis

Authors: Hiroki Waida, Yuichiro Wada, Léo Andéol, Takumi Nakagawa, Yuhui Zhang, Takafumi Kanamori

Abstract: Contrastive learning is an efficient approach to self-supervised representation learning. Although recent studies have made progress in the theoretical understanding of contrastive learning, the investigation of how to characterize the clusters of the learned representations is still limited. In this paper, we aim to elucidate the characterization from theoretical perspectives. To this end, we con… ▽ More Contrastive learning is an efficient approach to self-supervised representation learning. Although recent studies have made progress in the theoretical understanding of contrastive learning, the investigation of how to characterize the clusters of the learned representations is still limited. In this paper, we aim to elucidate the characterization from theoretical perspectives. To this end, we consider a kernel-based contrastive learning framework termed Kernel Contrastive Learning (KCL), where kernel functions play an important role when applying our theoretical results to other frameworks. We introduce a formulation of the similarity structure of learned representations by utilizing a statistical dependency viewpoint. We investigate the theoretical properties of the kernel-based contrastive loss via this formulation. We first prove that the formulation characterizes the structure of representations learned with the kernel-based contrastive learning framework. We show a new upper bound of the classification error of a downstream task, which explains that our theory is consistent with the empirical success of contrastive learning. We also establish a generalization error bound of KCL. Finally, we show a guarantee for the generalization ability of KCL to the downstream classification task via a surrogate bound. △ Less

Submitted 1 April, 2023; originally announced April 2023.

Effects of Thermal Modification on the Flexure Properties, Fracture Energy, and Hardness of Western Hemlock

Authors: Troy Nakagawa, Erik Poulin, Talbot Rueppel, Zhisong Chen, Juliet Swinea, Mark O'Brien, Guy Houser, Geoffrey Wood, Malloree Weinheimer, Pouria Bahmani, Peter Stynoski, Marco Salviato

Abstract: This study investigates the effect of thermal modification on the flexural properties, transverse fracture energy, and hardness of western hemlock, a material which is finding increasing applications in construction. Flexure tests on specimens featuring longitudinal and transverse grains showed that thermal modification at 167C slightly improves the flexural modulus and strength and leads to less… ▽ More This study investigates the effect of thermal modification on the flexural properties, transverse fracture energy, and hardness of western hemlock, a material which is finding increasing applications in construction. Flexure tests on specimens featuring longitudinal and transverse grains showed that thermal modification at 167C slightly improves the flexural modulus and strength and leads to less statistical variability compared to unmodified samples. On the other hand, the fracture and Janka hardness tests revealed a more pronounced brittleness of the thermally modified samples. In fact, the total mode I fracture energy of modified Single Edge Notch Bending (SENB) samples was about 47% lower for radial-longitudinal systems and 60% lower for tangential-longitudinal systems. Similarly, the average Janka hardness in the tangential, radial, and transverse planes was 8.5%, 3.9%, and 9.4% lower in the modified specimens, respectively. The results presented in this work show that thermal modification can have a significant effect on the fracturing behavior of western hemlock and its energy dissipation capabilities. For design, this must be taken into serious consideration as these properties significantly influence the damage tolerance of this wood in the presence of stress concentrations such as e.g., those induced in bolted joints and cut outs. Fracture energy and hardness are also strongly correlated to ballistic performance. △ Less

Submitted 31 March, 2023; originally announced April 2023.

Comments: 27 pages, 25 figures

Direct Imaging Explorations for Companions around Mid-Late M Stars from the Subaru/IRD Strategic Program

Authors: Taichi Uyama, Charles Beichman, Masayuki Kuzuhara, Markus Janson, Takayuki Kotani, Dimitri Mawet, Bun'ei Sato, Motohide Tamura, Hiroyuki Tako Ishikawa, Bryson Cale, Thayne Currie, Hiroki Harakawa, Thomas Henning, Teruyuki Hirano, Klaus Hodapp, Yasunori Hori, Masato Ishizuka, Shane Jacobson, Yui Kasagi, Eiichiro Kokubo, Mihoko Konishi, Tomoyuki Kudo, Takashi Kurokawa, Nobuhiko Kusakabe, Jungmi Kwon , et al. (10 additional authors not shown)

Abstract: The Subaru telescope is currently performing a strategic program (SSP) using the high-precision near-infrared (NIR) spectrometer IRD to search for exoplanets around nearby mid/late-M~dwarfs via radial velocity (RV) monitoring. As part of the observing strategy for the exoplanet survey, signatures of massive companions such as RV trends are used to reduce the priority of those stars. However, this… ▽ More The Subaru telescope is currently performing a strategic program (SSP) using the high-precision near-infrared (NIR) spectrometer IRD to search for exoplanets around nearby mid/late-M~dwarfs via radial velocity (RV) monitoring. As part of the observing strategy for the exoplanet survey, signatures of massive companions such as RV trends are used to reduce the priority of those stars. However, this RV information remains useful for studying the stellar multiplicity of nearby M~dwarfs. To search for companions around such ``deprioritized" M~dwarfs, we observed 14 IRD-SSP targets using Keck/NIRC2 observations with pyramid wavefront sensing at NIR wavelengths, leading to high sensitivity to substellar-mass companions within a few arcseconds. We detected two new companions (LSPM~J1002+1459~B and LSPM~J2204+1505~B) and two new candidates that are likely companions (LSPM~J0825+6902~B and LSPM~J1645+0444~B) as well as one known companion. Including two known companions resolved by the IRD fiber injection module camera, we detected seven (four new) companions at projected separations between $\sim2-20$~au in total. A comparison of the colors with the spectral library suggests that LSPM~J2204+1505~B and LSPM~J0825+6902~B are located at the boundary between late-M and early-L spectral types. Our deep high-contrast imaging for targets where no bright companions were resolved did not reveal any additional companion candidates. The NIRC2 detection limits could constrain potential substellar-mass companions ($\sim10-75\ M_{\rm Jup}$) at 10~au or further. The failure with Keck/NIRC2 around the IRD-SSP stars having significant RV trends makes these objects promising targets for further RV monitoring or deeper imaging with JWST to search for smaller-mass companions below the NIRC2 detection limits. △ Less

Submitted 26 February, 2023; originally announced February 2023.

Comments: 16 pages, 8 figures, accepted for publication in AJ

LEALLA: Learning Lightweight Language-agnostic Sentence Embeddings with Knowledge Distillation

Authors: Zhuoyuan Mao, Tetsuji Nakagawa

Abstract: Large-scale language-agnostic sentence embedding models such as LaBSE (Feng et al., 2022) obtain state-of-the-art performance for parallel sentence alignment. However, these large-scale models can suffer from inference speed and computation overhead. This study systematically explores learning language-agnostic sentence embeddings with lightweight models. We demonstrate that a thin-deep encoder ca… ▽ More Large-scale language-agnostic sentence embedding models such as LaBSE (Feng et al., 2022) obtain state-of-the-art performance for parallel sentence alignment. However, these large-scale models can suffer from inference speed and computation overhead. This study systematically explores learning language-agnostic sentence embeddings with lightweight models. We demonstrate that a thin-deep encoder can construct robust low-dimensional sentence embeddings for 109 languages. With our proposed distillation methods, we achieve further improvements by incorporating knowledge from a teacher model. Empirical results on Tatoeba, United Nations, and BUCC show the effectiveness of our lightweight models. We release our lightweight language-agnostic sentence embedding models LEALLA on TensorFlow Hub. △ Less

Submitted 26 December, 2023; v1 submitted 16 February, 2023; originally announced February 2023.

Comments: EACL 2023 main conference; LEALLA models: https://www.kaggle.com/models/google/lealla (modified url in v2 of this paper)

Effects of Average Number of Platelets Through the Thickness and Platelet Width on the Mechanical Properties of Discontinuous Fiber Composites

Authors: Seunghyun Ko, Troy Nakagawa, Zhisong Chen, William B. Avery, Ebonni J. Adams, Matthew R. Soja, Michael H. Larson, Chul Y. Park, Jinkyu Yang, Marco Salviato

Abstract: In this study, we experimentally and numerically investigate the evolution of the tensile material properties of Discontinuous Fiber Composites (DFCs) with an increasing average number of platelets through the thickness for two different platelet widths. The results show that both the number of platelets and the platelet width have significant effects on the tensile modulus and strength. We find t… ▽ More In this study, we experimentally and numerically investigate the evolution of the tensile material properties of Discontinuous Fiber Composites (DFCs) with an increasing average number of platelets through the thickness for two different platelet widths. The results show that both the number of platelets and the platelet width have significant effects on the tensile modulus and strength. We find that not only the average mechanical properties but also their coefficients of variation change according to the different DFC mesostructures. To understand the relationship between material morphology at the mesoscale and corresponding material properties, we developed a random platelet mesostructure generation algorithm combined with explicit finite element models. Leveraging the computational tools, we find that moduli and strength increase with increasing average number of platelets through the thickness. The increasing trend continues until reaching an asymptotic limit at about 45 layers through the thickness for the narrow platelets and 27 layers for the square platelets. In the study, we address the importance of having accurate simulations of the mesostructure to match not only the average modulus and strength but also their associated coefficients of variation. We show that it is possible to accurately predict the tensile material properties of DFCs, including their B-basis design values. This is a quintessential condition for the adoption of DFCs in structural applications. △ Less

Submitted 6 February, 2023; originally announced February 2023.

Report number: E-0223

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

Comment on Comparison of Holocene temperature reconstructions based on GISP2 multiple-gas-isotope measurements by Doring and Leuenberger (2022)

Authors: Takuro Kobashi, Tosiyuki Nakaegawa

Abstract: Doring and Lueneberger (Doring and Leuenberger, 2022) conducted an important work on Greenland temperature reconstruction from occluded air in an ice core for the Holocene by comparing various reconstructed temperatures. From their analyses, they conclude that their d15N temperature reconstruction (d15N method) is the best method among the compared methods, and Kobashi et al. temperature reconstru… ▽ More Doring and Lueneberger (Doring and Leuenberger, 2022) conducted an important work on Greenland temperature reconstruction from occluded air in an ice core for the Holocene by comparing various reconstructed temperatures. From their analyses, they conclude that their d15N temperature reconstruction (d15N method) is the best method among the compared methods, and Kobashi et al. temperature reconstruction with d15Nexcess (Kobashi et al., 2017) should not be used for climate interpretation because d40Ar and d15Nexcess data are not accurate enough. However, they failed to recognize important issues to reach their conclusions. Therefore, we comment on those points to make it clear that their conclusions are erroneous. △ Less

Submitted 16 November, 2022; originally announced November 2022.

Journal ref: Quaternary Science Reviews, Volume 298,107707, 2022

Efficient search for superspecial hyperelliptic curves of genus four with automorphism group containing $\mathbb{Z}_6$

Authors: Momonari Kudo, Tasuku Nakagawa, Tsuyoshi Takagi

Abstract: In arithmetic and algebraic geometry, superspecial (s.sp.\ for short) curves are one of the most important objects to be studied, with applications to cryptography and coding theory. If $g \geq 4$, it is not even known whether there exists such a curve of genus $g$ in general characteristic $p > 0$, and in the case of $g=4$, several computational approaches to search for those curves have been pro… ▽ More In arithmetic and algebraic geometry, superspecial (s.sp.\ for short) curves are one of the most important objects to be studied, with applications to cryptography and coding theory. If $g \geq 4$, it is not even known whether there exists such a curve of genus $g$ in general characteristic $p > 0$, and in the case of $g=4$, several computational approaches to search for those curves have been proposed. In the genus-$4$ hyperelliptic case, Kudo-Harashita proposed a generic algorithm to enumerate all s.sp.\ curves, and recently Ohashi-Kudo-Harashita presented an algorithm specific to the case where automorphism group contains the Klein 4-group. In this paper, we propose an algorithm with complexity $\tilde{O}(p^4)$ in theory but $\tilde{O}(p^3)$ in practice to enumerate s.sp.\ hyperelliptic curves of genus 4 with automorphism group containing the cyclic group of order $6$. By executing the algorithm over Magma, we enumerate those curves for $p$ up to $1000$. We also succeeded in finding a s.sp.\ hyperelliptic curve of genus $4$ in every $p$ with $p \equiv 2 \pmod{3}$. As a theoretical result, we classify hyperelliptic curves of genus $4$ in terms of automorphism groups in the appendix. △ Less

Submitted 26 October, 2022; originally announced October 2022.

Robustness against outliers in ordinal response model via divergence approach

Authors: Tomotaka Momozaki, Tomoyuki Nakagawa

Abstract: This study deals with the problem of outliers in ordinal response model, which is a regression on ordered categorical data as the response variable. ``Outlier" means that the combination of ordered categorical data and its covariates is heterogeneous compared to other pairs. Although the ordinal response model is important for data analysis in various fields such as medicine and social sciences, i… ▽ More This study deals with the problem of outliers in ordinal response model, which is a regression on ordered categorical data as the response variable. ``Outlier" means that the combination of ordered categorical data and its covariates is heterogeneous compared to other pairs. Although the ordinal response model is important for data analysis in various fields such as medicine and social sciences, it is known that the maximum likelihood method with probit, logit, log-log and complementary log-log link functions, which are often used, is strongly affected by outliers, and statistical analysts are forced to limit their analysis when there may be outliers in the data. To solve this problem, this paper provides inference methods with two robust divergences (the density-power and $γ$-divergences). We also derive influence functions for the proposed methods and show conditions on the link function for them to be bounded and to redescendence. Since the commonly used link functions satisfy these conditions, the analyst can perform robust and flexible analysis with our methods. In addition, and this is a result that further highlights our contributions, we show that the influence function in the maximum likelihood method does not have redescendence for any link function in the ordinal response model. Through numerical experiments using artificial and two real data, we show that the proposed methods perform better than the maximum likelihood method with and without outliers in the data for various link functions. △ Less

Submitted 27 December, 2022; v1 submitted 24 September, 2022; originally announced September 2022.

Comments: 30 pages

Expressions Causing Differences in Emotion Recognition in Social Networking Service Documents

Authors: Tsubasa Nakagawa, Shunsuke Kitada, Hitoshi Iyatomi

Abstract: It is often difficult to correctly infer a writer's emotion from text exchanged online, and differences in recognition between writers and readers can be problematic. In this paper, we propose a new framework for detecting sentences that create differences in emotion recognition between the writer and the reader and for detecting the kinds of expressions that cause such differences. The proposed f… ▽ More It is often difficult to correctly infer a writer's emotion from text exchanged online, and differences in recognition between writers and readers can be problematic. In this paper, we propose a new framework for detecting sentences that create differences in emotion recognition between the writer and the reader and for detecting the kinds of expressions that cause such differences. The proposed framework consists of a bidirectional encoder representations from transformers (BERT)-based detector that detects sentences causing differences in emotion recognition and an analysis that acquires expressions that characteristically appear in such sentences. The detector, based on a Japanese SNS-document dataset with emotion labels annotated by both the writer and three readers of the social networking service (SNS) documents, detected "hidden-anger sentences" with AUC = 0.772; these sentences gave rise to differences in the recognition of anger. Because SNS documents contain many sentences whose meaning is extremely difficult to interpret, by analyzing the sentences detected by this detector, we obtained several expressions that appear characteristically in hidden-anger sentences. The detected sentences and expressions do not convey anger explicitly, and it is difficult to infer the writer's anger, but if the implicit anger is pointed out, it becomes possible to guess why the writer is angry. Put into practical use, this framework would likely have the ability to mitigate problems based on misunderstandings. △ Less

Submitted 3 September, 2022; v1 submitted 30 August, 2022; originally announced August 2022.

Comments: 5 pages, 3 figures. Accepted at the 31st ACM International Conference on Information and Knowledge Management (CIKM '22) as a short paper

Journal ref: Proceedings of the 31st ACM International Conference on Information and Knowledge Management (CIKM'22), October 17--21, 2022, Atlanta, GA, USA

Simulations of the Spectral Resolving Power of a Compact Space-Borne Immersion-Echelle Spectrometer Using Mid-Infrared Wave Tracing

Authors: Satoshi Itoh, Daisuke Ishihara, Takehiko Wada, Takao Nakagawa, Shinki Oyabu, Hidehiro Kaneda, Yasuhiro Hirahara, the SMI consortium

Abstract: We performed wave-optics-based numerical simulations at mid-infrared wavelengths to investigate how the presence or absence of entrance slits and optical aberrations affect the spectral resolving power $R$ of a compact, high-spectral-resolving-power spectrometer containing an immersion-echelle grating. We tested three cases of telescope aberration (aberration-free, astigmatism and spherical aberra… ▽ More We performed wave-optics-based numerical simulations at mid-infrared wavelengths to investigate how the presence or absence of entrance slits and optical aberrations affect the spectral resolving power $R$ of a compact, high-spectral-resolving-power spectrometer containing an immersion-echelle grating. We tested three cases of telescope aberration (aberration-free, astigmatism and spherical aberration), assuming the aberration budget of the Space Infrared Telescope for Cosmology and Astrophysics (SPICA), which has a 20-$\mathrm{μm}$-wavelength diffraction limit. In cases with a slit, we found that the value of $R$ at around 10--20 $\mathrm{μm}$ is approximately independent of the assumed aberrations, which is significantly different from the prediction of geometrical optics. Our results also indicate that diffraction from the slit improves $R$ by enlarging the effective illuminated area on the grating window and that this improvement decreases at short wavelengths. For the slit-less cases, we found that the impact of aberrations on $R$ can be roughly estimated using the Strehl ratio. △ Less

Submitted 6 June, 2022; v1 submitted 6 June, 2022; originally announced June 2022.

Comments: Accepted for publication in the Journal of Astronomical Telescopes, Instruments, and Systems (JATIS)

Probing dynamics and thermal properties inside molecular tori with CO rovibrational absorption lines

Authors: Kosei Matsumoto, Takao Nakagawa, Keiichi Wada, Shunsuke Baba, Shusuke Onishi, Taisei Uzuo, Naoki Isobe, Yuki Kudoh

Abstract: A recent hydrodynamic model, "radiation-driven fountain model" (Wada et al. 2016), presented a dynamical picture that active galactic nuclei (AGNs) tori sustain their geometrical thickness by gas circulation around AGNs, and previous papers confirmed that this picture is consistent with multi-wavelength observations of nearby Seyfert galaxies. Recent near-infrared observations implied that CO rovi… ▽ More A recent hydrodynamic model, "radiation-driven fountain model" (Wada et al. 2016), presented a dynamical picture that active galactic nuclei (AGNs) tori sustain their geometrical thickness by gas circulation around AGNs, and previous papers confirmed that this picture is consistent with multi-wavelength observations of nearby Seyfert galaxies. Recent near-infrared observations implied that CO rovibrational absorption lines ($ΔJ=\pm1$, $v=0-1$, $λ\sim 4.7$ $\mathrm{μm}$) could probe physical properties of the inside tori. However, the origin of the CO absorption lines has been under debate. In this paper, we investigate the origin of the absorption lines and conditions for detecting them by performing line radiative transfer calculations based on the radiation-driven fountain model. We find that CO rovibrational absorption lines are detected at inclination angles $θ_\mathrm{obs} = 50-80$ $^{\circ}$. At the inclination angle $θ_\mathrm{obs} = 77$ $^{\circ}$, we observe multi-velocity components: inflow ($v_\mathrm{LOS}=30$ $\mathrm{kms^{-1}}$), systemic ($v_\mathrm{LOS}=0 \, \mathrm{kms^{-1}}$), and outflows ($v_\mathrm{LOS}=-75,\, -95,$ and $-105$ $\mathrm{kms^{-1}}$). The inflow and outflow components ($v_\mathrm{LOS}= 30$ and $-95$ $\mathrm{kms^{-1}}$) are collisionally excited at the excitation temperature of $186$ and $380$ K up to $J=12$ and $4$, respectively. The inflow and outflow components originate from the accreting gas on the equatorial plane at $1.5$ pc from the AGN center and the outflowing gas driven by AGN radiation pressure at $1.0$ pc, respectively. These results suggest that CO rovibrational absorption lines can provide us with the velocities and kinetic temperatures of the inflow and outflow in the inner a-few-pc regions of AGN tori, and the observations can probe the gas circulation inside the tori. △ Less

Submitted 1 June, 2022; originally announced June 2022.

Comments: 19 pages, 17 figures, 3 tables, accepted by The Astrophysical Journal

A Super-Earth Orbiting Near the Inner Edge of the Habitable Zone around the M4.5-dwarf Ross 508

Authors: Hiroki Harakawa, Takuya Takarada, Yui Kasagi, Teruyuki Hirano, Takayuki Kotani, Masayuki Kuzuhara, Masashi Omiya, Hajime Kawahara, Akihiko Fukui, Yasunori Hori, Hiroyuki Tako Ishikawa, Masahiro Ogihara, John Livingston, Timothy D. Brandt, Thayne Currie, Wako Aoki, Charles A. Beichman, Thomas Henning, Klaus Hodapp, Masato Ishizuka, Hideyuki Izumiura, Shane Jacobson, Markus Janson, Eiji Kambe, Takanori Kodama , et al. (24 additional authors not shown)

Abstract: We report the near-infrared radial-velocity (RV) discovery of a super-Earth planet on a 10.77-day orbit around the M4.5 dwarf Ross 508 ($J_\mathrm{mag}=9.1$). Using precision RVs from the Subaru Telescope IRD (InfraRed Doppler) instrument, we derive a semi-amplitude of $3.92^{+0.60}_{-0.58}$ ${\rm m\,s}^{-1}$, corresponding to a planet with a minimum mass… ▽ More We report the near-infrared radial-velocity (RV) discovery of a super-Earth planet on a 10.77-day orbit around the M4.5 dwarf Ross 508 ($J_\mathrm{mag}=9.1$). Using precision RVs from the Subaru Telescope IRD (InfraRed Doppler) instrument, we derive a semi-amplitude of $3.92^{+0.60}_{-0.58}$ ${\rm m\,s}^{-1}$, corresponding to a planet with a minimum mass $m \sin i = 4.00^{+0.53}_{-0.55}\ M_{\oplus}$. We find no evidence of significant signals at the detected period in spectroscopic stellar activity indicators or MEarth photometry. The planet, Ross 508 b, has a semimajor-axis of $0.05366^{+0.00056}_{-0.00049}$ au. This gives an orbit-averaged insolation of $\approx$1.4 times the Earth's value, placing Ross 508 b near the inner edge of its star's habitable zone. We have explored the possibility that the planet has a high eccentricity and its host is accompanied by an additional unconfirmed companion on a wide orbit. Our discovery demonstrates that the near-infrared RV search can play a crucial role to find a low-mass planet around cool M dwarfs like Ross 508. △ Less

Submitted 24 May, 2022; originally announced May 2022.

Comments: Accepted for publication in PASJ (May 23, 2022)

Generalized Cramér's coefficient via $f$-divergence for contingency tables

Authors: Wataru Urasaki, Tomoyuki Nakagawa, Tomotaka Momozaki, Sadao Tomizawa

Abstract: Various measures in two-way contingency table analysis have been proposed to express the strength of association between row and column variables in contingency tables. Tomizawa et al. (2004) proposed more general measures, including Cramér's coefficient, using the power-divergence. In this paper, we propose measures using the $f$-divergence that has a wider class than the power-divergence. Unlike… ▽ More Various measures in two-way contingency table analysis have been proposed to express the strength of association between row and column variables in contingency tables. Tomizawa et al. (2004) proposed more general measures, including Cramér's coefficient, using the power-divergence. In this paper, we propose measures using the $f$-divergence that has a wider class than the power-divergence. Unlike statistical hypothesis tests, these measures provide quantification of the association structure in contingency tables. The contribution of our study is proving that a measure applying a function that satisfies the condition of the $f$-divergence has desirable properties for measuring the strength of association in contingency tables. With this contribution, we can easily construct a new measure using a divergence that has essential properties for the analyst. For example, we conducted numerical experiments with a measure applying the $θ$-divergence. Furthermore, we can give further interpretation of the association between the row and column variables in the contingency table, which could not be obtained with the conventional one. We also show a relationship between our proposed measures and the correlation coefficient in the bivariate normal distribution of latent variables in the contingency tables. △ Less

Submitted 6 July, 2023; v1 submitted 25 April, 2022; originally announced April 2022.

Comments: 36 pages

NeuRecover: Regression-Controlled Repair of Deep Neural Networks with Training History

Authors: Shogo Tokui, Susumu Tokumoto, Akihito Yoshii, Fuyuki Ishikawa, Takao Nakagawa, Kazuki Munakata, Shinji Kikuchi

Abstract: Systematic techniques to improve quality of deep neural networks (DNNs) are critical given the increasing demand for practical applications including safety-critical ones. The key challenge comes from the little controllability in updating DNNs. Retraining to fix some behavior often has a destructive impact on other behavior, causing regressions, i.e., the updated DNN fails with inputs correctly h… ▽ More Systematic techniques to improve quality of deep neural networks (DNNs) are critical given the increasing demand for practical applications including safety-critical ones. The key challenge comes from the little controllability in updating DNNs. Retraining to fix some behavior often has a destructive impact on other behavior, causing regressions, i.e., the updated DNN fails with inputs correctly handled by the original one. This problem is crucial when engineers are required to investigate failures in intensive assurance activities for safety or trust. Search-based repair techniques for DNNs have potentials to tackle this challenge by enabling localized updates only on "responsible parameters" inside the DNN. However, the potentials have not been explored to realize sufficient controllability to suppress regressions in DNN repair tasks. In this paper, we propose a novel DNN repair method that makes use of the training history for judging which DNN parameters should be changed or not to suppress regressions. We implemented the method into a tool called NeuRecover and evaluated it with three datasets. Our method outperformed the existing method by achieving often less than a quarter, even a tenth in some cases, number of regressions. Our method is especially effective when the repair requirements are tight to fix specific failure types. In such cases, our method showed stably low rates (<2%) of regressions, which were in many cases a tenth of regressions caused by retraining. △ Less

Submitted 4 March, 2022; v1 submitted 28 February, 2022; originally announced March 2022.

Extremely Buried Nucleus of IRAS 17208$-$0014 Observed at Sub-Millimeter and Near-Infrared Wavelengths

Authors: Shunsuke Baba, Masatoshi Imanishi, Takuma Izumi, Taiki Kawamuro, Dieu D. Nguyen, Takao Nakagawa, Naoki Isobe, Shusuke Onishi, Kosei Matsumoto

Abstract: The ultraluminous infrared galaxy IRAS 17208$-$0014 is a late-stage merger that hosts a buried active galactic nucleus (AGN). To investigate its nuclear structure, we performed high spatial resolution ($\sim0.\!\!^{\prime\prime}04\sim32\,\mathrm{pc}$) Atacama Large Millimeter/submillimeter Array (ALMA) observations in Band 9 ($\sim$450\,\micron\ or $\sim$660\,GHz), along with near-infrared AKARI s… ▽ More The ultraluminous infrared galaxy IRAS 17208$-$0014 is a late-stage merger that hosts a buried active galactic nucleus (AGN). To investigate its nuclear structure, we performed high spatial resolution ($\sim0.\!\!^{\prime\prime}04\sim32\,\mathrm{pc}$) Atacama Large Millimeter/submillimeter Array (ALMA) observations in Band 9 ($\sim$450\,\micron\ or $\sim$660\,GHz), along with near-infrared AKARI spectroscopy in 2.5--5.0\,\micron. The Band 9 dust continuum peaks at the AGN location, and toward this position CO($J$=6--5) and CS($J$=14--13) are detected in absorption. Comparison with non-local thermal equilibrium calculations indicates that, within the central beam ($r\sim20\,\mathrm{pc}$), there exists a concentrated component that is dense ($10^7\,\mathrm{cm}^{-2}$) and warm ($>$200\,K) and has a large column density ($N_\mathrm{H_2}>10^{23}\,\mathrm{cm}^{-2}$). The AKARI spectrum shows deep and broad CO ro-vibrational absorption at 4.67\,\micron. Its band profile is well reproduced with a similarly dense and large column but hotter ($\sim$1000\,K) gas. The region observed through absorption in the near-infrared is highly likely in the nuclear direction, as in the sub-millimeter, but with a narrower beam including a region closer to the nucleus. The central component is considered to possess a hot structure where vibrationally excited HCN emission originates. The most plausible heating source for the gas is X-rays from the AGN. The AKARI spectrum does not show other AGN signs in 2.5--4\,\micron, but this absence may be usual for AGNs buried in a hot mid-infrared core. Besides, based on our ALMA observations, we relate various nuclear structures of IRAS 17208$-$0014 that have been proposed in the literature. △ Less

Submitted 22 February, 2022; originally announced February 2022.

Comments: 28 pages, 16 figures, accepted for publication in ApJ

Effects of Out Time on the Mechanical and Fracture Properties of Chopped Fiber Composites Made From Repurposed Aerospace Prepreg Scrap and Waste

Authors: Troy Nakagawa, Seunghyun Ko, Cory Slaughter, Talal Abdullah, Guy Houser, Marco Salviato

Abstract: In this study, the effects of prepreg out time on the mechanical and fracture properties of Discontinuous Fiber Composites (DFCs) are investigated. Carbon fiber prepregs are aged at 0x, 1x, 2x, and 3x the out life in an environmental chamber at constant temperature and humidity. Degree of cure is measured via Differential Scanning Calorimetry (DSC) while tension, compression, and shear tests are p… ▽ More In this study, the effects of prepreg out time on the mechanical and fracture properties of Discontinuous Fiber Composites (DFCs) are investigated. Carbon fiber prepregs are aged at 0x, 1x, 2x, and 3x the out life in an environmental chamber at constant temperature and humidity. Degree of cure is measured via Differential Scanning Calorimetry (DSC) while tension, compression, and shear tests are performed to investigate the effects that aging has on these mechanical properties. For the first time, Mode I intra-laminar fracture and its size effect are also investigated by means of fracture tests on geometrically-scaled Single Edge Notch Tension (SENT) specimens. From the tension, compression, and shear experiments it is seen that the out time has no effect on the elastic moduli. However, the strength increases with increasing age of the specimens for all the loading conditions. The percent increase compared to the non-aged material ranges from 15% to 33%. This is likely caused by plasticization of the matrix with age, allowing for higher energy absorption. More complex trends are reported for the SENT specimens for all the sizes. It is found that the fracture energy and characteristic length initially decrease with age, and then finally increase for the longest out time. This trend is owed to two factors with countering effects on the fracture behavior: 1) the increase of the average number of platelets with increasing aging due to increase in resin viscosity, and 2) the plasticization of the matrix with aging. The results from this study suggest that Discontinuous Fiber Composites (DFCs) made from reused materials can have equal, if not better, performance than non-aged DFCs. The experimental data presented in this work can be used as a baseline to design DFC composite components made from repurposed prepreg scrap and waste. △ Less

Submitted 21 December, 2021; originally announced December 2021.

Infrared Absorption and its Sources of CdZnTe at Cryogenic Temperature

Authors: Hiroshi Maeshima, Kosei Matsumoto, Yasuhiro Hirahara, Takao Nakagawa, Ryoichi Koga, Yusuke Hanamura, Takehiko Wada, Koichi Nagase, Shinki Oyabu, Toyoaki Suzuki, Takuma Kokusho, Hidehiro Kaneda, Daichi Ishikawa

Abstract: To reveal the infrared absorption causes in the wavelength region between electronic and lattice absorptions, we measured the temperature dependence of the absorption coefficient of $p$-type low-resistivity ($\sim 10^2~{\rm Ωcm}$) CdZnTe crystals. We measured the absorption coefficients of CdZnTe crystals in four-wavelength bands ($λ=6.45$, 10.6, 11.6, 15.1$~μ$m) over the temperature range of… ▽ More To reveal the infrared absorption causes in the wavelength region between electronic and lattice absorptions, we measured the temperature dependence of the absorption coefficient of $p$-type low-resistivity ($\sim 10^2~{\rm Ωcm}$) CdZnTe crystals. We measured the absorption coefficients of CdZnTe crystals in four-wavelength bands ($λ=6.45$, 10.6, 11.6, 15.1$~μ$m) over the temperature range of $T=8.6-300$ K with an originally developed system. The CdZnTe absorption coefficient was measured to be $α=0.3-0.5$ ${\rm cm^{-1}}$ at $T=300$ K and $α=0.4-0.9$ ${\rm cm^{-1}}$ at $T=8.6$ K in the investigated wavelength range. With an absorption model based on transitions of free holes and holes trapped at an acceptor level, we conclude that the absorption due to free holes at $T=150-300$ K and that due to trapped-holes at $T<50$ K are dominant absorption causes in CdZnTe. We also discuss a method to predict the CdZnTe absorption coefficient at cryogenic temperature based on the room-temperature resistivity. △ Less

Submitted 18 December, 2021; originally announced December 2021.

Comments: Accepted for publication in Journal of Electronic Material, 30 pages, 6 figures, 5 tables

Elemental abundances of nearby M dwarfs based on high-resolution near-infrared spectra obtained by the Subaru/IRD survey: Proof of concept

Authors: Hiroyuki Tako Ishikawa, Wako Aoki, Teruyuki Hirano, Takayuki Kotani, Masayuki Kuzuhara, Masashi Omiya, Yasunori Hori, Eiichiro Kokubo, Tomoyuki Kudo, Takashi Kurokawa, Nobuhiko Kusakabe, Norio Narita, Jun Nishikawa, Masahiro Ogihara, Akitoshi Ueda, Thayne Currie, Thomas Henning, Yui Kasagi, Jared R. Kolecki, Jungmi Kwon, Masahiro N. Machida, Michael W. McElwain, Takao Nakagawa, Sebastien Vievard, Ji Wang , et al. (2 additional authors not shown)

Abstract: Detailed chemical analyses of M dwarfs are scarce but necessary to constrain the formation environment and internal structure of planets being found around them. We present elemental abundances of 13 M dwarfs (2900 < Teff < 3500 K) observed in the Subaru/IRD planet search project. They are mid-to-late M dwarfs whose abundance of individual elements has not been well studied. We use the high-resolu… ▽ More Detailed chemical analyses of M dwarfs are scarce but necessary to constrain the formation environment and internal structure of planets being found around them. We present elemental abundances of 13 M dwarfs (2900 < Teff < 3500 K) observed in the Subaru/IRD planet search project. They are mid-to-late M dwarfs whose abundance of individual elements has not been well studied. We use the high-resolution (~70,000) near-infrared (970-1750 nm) spectra to measure the abundances of Na, Mg, Si, K, Ca, Ti, V, Cr, Mn, Fe, and Sr by the line-by-line analysis based on model atmospheres, with typical errors ranging from 0.2 dex for [Fe/H] to 0.3-0.4 dex for other [X/H]. We measure radial velocities from the spectra and combine them with Gaia astrometry to calculate the Galactocentric space velocities UVW. The resulting [Fe/H] values agree with previous estimates based on medium-resolution K-band spectroscopy, showing a wide distribution of metallicity (-0.6 < [Fe/H] < +0.4). The abundance ratios of individual elements [X/Fe] are generally aligned with the solar values in all targets. While the [X/Fe] distributions are comparable to those of nearby FGK stars, most of which belong to the thin disk population, the most metal-poor object, GJ 699, could be a thick disk star. The UVW velocities also support this. The results raise the prospect that near-infrared spectra of M dwarfs obtained in the planet search projects can be used to grasp the trend of elemental abundances and Galactic stellar population of nearby M dwarfs. △ Less

Submitted 30 November, 2021; originally announced December 2021.

Comments: 31 pages, 16 figures, accepted for publication in AJ

Electron g-factor determined for quantum dot circuit fabricated from (110)-oriented GaAs quantum well

Authors: T. Nakagawa, S. Lamoureux, T. Fujita, J. Ritzmann, A. Ludwig, A. D. Wieck, A. Oiwa, M. Korkusinski, A. Sachrajda, D. G. Austing, L. Gaudreau

Abstract: The choice of substrate orientation for semiconductor quantum dot circuits offers opportunities for tailoring spintronic properties such as g-factors for specific functionality. In this letter, we demonstrate the operation of a few-electron double quantum dot circuit fabricated from a (110)-oriented GaAs quantum well. We estimate the in-plane electron g-factor from the profile of the enhanced inte… ▽ More The choice of substrate orientation for semiconductor quantum dot circuits offers opportunities for tailoring spintronic properties such as g-factors for specific functionality. In this letter, we demonstrate the operation of a few-electron double quantum dot circuit fabricated from a (110)-oriented GaAs quantum well. We estimate the in-plane electron g-factor from the profile of the enhanced inter-dot tunneling (leakage) current near zero magnetic field. Spin-blockade due to Pauli exclusion can block inter-dot tunneling. However, this blockade becomes inactive due to hyperfine interaction mediated spin flip-flop processes between electron spin states and the nuclear spin of the host material. The g-factor of absolute value ~0.1 found for a magnetic field parallel to the direction [11(bar)0], is approximately a factor of four lower than that for comparable circuits fabricated from material grown on widely-employed standard (001) GaAs substrates, and is in line with reported values determined by purely optical means for quantum well structures grown on (110) GaAs substrates. △ Less

Submitted 14 November, 2021; originally announced November 2021.

Comments: 17 pages, 4 figures

Dynamic Spatio-temporal Zero-inflated Poisson Models for Predicting Capelin Distribution in the Barents Sea

Authors: Shonosuke Sugasawa, Tomoyuki Nakagawa, Hiroko Kato Solvang, Sam Subbey, Salah Alrabeei

Abstract: We consider modeling and prediction of Capelin distribution in the Barents sea based on zero-inflated count observation data that vary continuously over a specified survey region. The model is a mixture of two components; a one-point distribution at the origin and a Poisson distribution with spatio-temporal intensity, where both intensity and mixing proportions are modeled by some auxiliary variab… ▽ More We consider modeling and prediction of Capelin distribution in the Barents sea based on zero-inflated count observation data that vary continuously over a specified survey region. The model is a mixture of two components; a one-point distribution at the origin and a Poisson distribution with spatio-temporal intensity, where both intensity and mixing proportions are modeled by some auxiliary variables and unobserved spatio-temporal effects. The spatio-temporal effects are modeled by a dynamic linear model combined with the predictive Gaussian process. We develop an efficient posterior computational algorithm for the model using a data augmentation strategy. The performance of the proposed model is demonstrated through simulation studies, and an application to the number of Capelin caught in the Barents sea from 2014 to 2019. △ Less

Submitted 19 October, 2022; v1 submitted 1 November, 2021; originally announced November 2021.

Comments: 25 pages

Dust dissipation timescales in the intermediate and outer regions of protoplanetary disks

Authors: Hiroshi Maeshima, Takao Nakagawa, Takuya Kojima, Satoshi Takita, Jungmi Kwon

Abstract: Dust and gas in protoplanetary disks dissipate as central stars evolve. In order to estimate the dust dissipation timescales in the protoplanetary disks, we stacked the WISE 12, 22, and the AKARI 90 $μ$m survey images of known T Tauri stars and derived the average fluxes, well below the survey flux limit in the 90 $μ$m band. We classified 4,783 T Tauri stars into three age groups, which are young… ▽ More Dust and gas in protoplanetary disks dissipate as central stars evolve. In order to estimate the dust dissipation timescales in the protoplanetary disks, we stacked the WISE 12, 22, and the AKARI 90 $μ$m survey images of known T Tauri stars and derived the average fluxes, well below the survey flux limit in the 90 $μ$m band. We classified 4,783 T Tauri stars into three age groups, which are young ($<$2 Myr), mid-age (2-6 Myr), and old ($>$6 Myr) groups, and stacked the WISE 12 and 22 and the AKARI 90 $μ$m images in each group. The photometry of the stacked image shows the flux decay timescales of 1.4$\pm$0.2, 1.38$\pm$0.05, and 1.4$^{+0.6}_{-0.5}$ Myr in the 12, 22, and 90 $μ$m bands, respectively. In optically thin disks with one-solar luminosity central stars, the 12 and 22 $μ$m fluxes are attributed to the emission from the intermediate ($\sim$ 1 au) region and the 90 $μ$m flux corresponds to that from the outer ($\sim$ 10 au) region in the disk. We hence conclude the dust dissipation timescale is $τ_{\rm med,dust}\sim$1.4 Myr in the intermediate disks and is $τ_{\rm outer, dust}=$1.4$^{+0.6}_{-0.5}$ Myr in the outer disks. The dust-dissipation time difference between the outer and intermediate disks is $Δτ_{\rm dust} = τ_{\rm outer,dust} -τ_{\rm med,dust}=0.0^{+0.6}_{-0.5}$ Myr, indicating that the dust in the intermediate and outer disks dissipates on almost the same timescale. △ Less

Submitted 4 October, 2021; originally announced October 2021.

Comments: 29 pages, 11 figures, 5 tables, accepted for publication in PASJ

Improving the accuracy of estimating indexes in contingency tables using Bayesian estimators

Authors: Tomotaka Momozaki, Koji Cho, Tomoyuki Nakagawa, Sadao Tomizawa

Abstract: In contingency table analysis, one is interested in testing whether a model of interest (e.g., the independent or symmetry model) holds using goodness-of-fit tests. When the null hypothesis where the model is true is rejected, the interest turns to the degree to which the probability structure of the contingency table deviates from the model. Many indexes have been studied to measure the degree of… ▽ More In contingency table analysis, one is interested in testing whether a model of interest (e.g., the independent or symmetry model) holds using goodness-of-fit tests. When the null hypothesis where the model is true is rejected, the interest turns to the degree to which the probability structure of the contingency table deviates from the model. Many indexes have been studied to measure the degree of the departure, such as the Yule coefficient and Cramér coefficient for the independence model, and Tomizawa's symmetry index for the symmetry model. The inference of these indexes is performed using sample proportions, which are estimates of cell probabilities, but it is well-known that the bias and mean square error (MSE) values become large without a sufficient number of samples. To address the problem, this study proposes a new estimator for indexes using Bayesian estimators of cell probabilities. Assuming the Dirichlet distribution for the prior of cell probabilities, we asymptotically evaluate the value of MSE when plugging the posterior means of cell probabilities into the index, and propose an estimator of the index using the Dirichlet hyperparameter that minimizes the value. Numerical experiments show that when the number of samples per cell is small, the proposed method has smaller values of bias and MSE than other methods of correcting estimation accuracy. We also show that the values of bias and MSE are smaller than those obtained by using the uniform and Jeffreys priors. △ Less

Submitted 25 October, 2023; v1 submitted 20 September, 2021; originally announced September 2021.

Comments: 19 pages, 6 figures

Anomalous Hydrogen Recombination-Line Ratios in Ultraluminous Infrared Galaxies

Authors: Kenichi Yano, Shunsuke Baba, Takao Nakagawa, Matthew Malkan, Naoki Isobe, Mai Shirahata, Ryosuke Doi, Vanshree Bhalotia

Abstract: We conducted systematic observations of the H I Br$α$ (4.05 $μ$m) and Br$β$ (2.63 $μ$m) lines in 52 nearby ($z<0.3$) ultraluminous infrared galaxies (ULIRGs) with AKARI. Among 33 ULIRGs wherein the lines are detected, three galaxies show anomalous Br$β$/Br$α$ line ratios ($\sim1.0$), which are significantly higher than those for case B (0.565). Our observations also show that ULIRGs have a tendenc… ▽ More We conducted systematic observations of the H I Br$α$ (4.05 $μ$m) and Br$β$ (2.63 $μ$m) lines in 52 nearby ($z<0.3$) ultraluminous infrared galaxies (ULIRGs) with AKARI. Among 33 ULIRGs wherein the lines are detected, three galaxies show anomalous Br$β$/Br$α$ line ratios ($\sim1.0$), which are significantly higher than those for case B (0.565). Our observations also show that ULIRGs have a tendency to exhibit higher Br$β$/Br$α$ line ratios than those observed in Galactic H II regions. The high Br$β$/Br$α$ line ratios cannot be explained by a combination of dust extinction and case B since dust extinction reduces the ratio. We explore possible causes for the high Br$β$/Br$α$ line ratios and show that the observed ratios can be explained by a combination of an optically thick Br$α$ line and an optically thin Br$β$ line. We simulated the H II regions in ULIRGs with the Cloudy code, and our results show that the high Br$β$/Br$α$ line ratios can be explained by high-density conditions, wherein the Br$α$ line becomes optically thick. To achieve a column density large enough to make the Br$α$ line optically thick within a single H II region, the gas density must be as high as $n\sim10^8$ $\mathrm{cm}^{-3}$. We therefore propose an ensemble of H II regions, in each of which the Br$α$ line is optically thick, to explain the high Br$β$/Br$α$ line ratio. △ Less

Submitted 14 September, 2021; originally announced September 2021.

Comments: 28 pages, 15 figures, accepted for publication in ApJ

Study of the Inner Structure of the Molecular Torus in IRAS 08572+3915 NW with Velocity Decomposition of CO Rovibrational Absorption Lines

Authors: Shusuke Onishi, Takao Nakagawa, Shunsuke Baba, Kosei Matsumoto, Naoki Isobe, Mai Shirahata, Hiroshi Terada, Tomonori Usuda, Shinki Oyabu

Abstract: Understanding the inner structure of the clumpy molecular torus surrounding the active galactic nucleus is essential in revealing the forming mechanism. However, spatially resolving the torus is difficult because of its size of a few parsecs. Thus, to probe the clump conditions in the torus, we performed the velocity decomposition of the CO rovibrational absorption lines (… ▽ More Understanding the inner structure of the clumpy molecular torus surrounding the active galactic nucleus is essential in revealing the forming mechanism. However, spatially resolving the torus is difficult because of its size of a few parsecs. Thus, to probe the clump conditions in the torus, we performed the velocity decomposition of the CO rovibrational absorption lines ($Δ{v}=0\to 1,\ Δ{J}=\pm 1$) at $λ\sim 4.67\,\mathrm{μ{m}}$ observed toward an ultraluminous infrared galaxy IRAS 08572+3915 NW with the high-resolution spectroscopy ($R\sim 10{,}000$) of Subaru Telescope. Consequently, we found that each transition had two outflowing components, i.e., (a) and (b), both at approximately $\sim -160\,\mathrm{km\,s^{-1}}$, but with broad and narrow widths, and an inflowing component, i.e., (c), at approximately $\sim +100\,\mathrm{km\,s^{-1}}$, which were attributed to the torus. The ratios of the velocity dispersions of each component lead to those of the rotating radii around the black hole of $R_\mathrm{rot,a}:R_\mathrm{rot,b}:R_\mathrm{rot,c}\approx 1:5:17$, indicating the torus where clumps are outflowing in the inner regions and inflowing in the outer regions if a hydrostatic disk with $σ_V\propto R_\mathrm{rot}^{-0.5}$ is assumed. Based on the kinetic temperature of components (a) and (b) of $\sim 720\,\mathrm{K}$ and $\sim 25\,\mathrm{K}$ estimated from the level population, the temperature gradient is $T_\mathrm{kin}\propto R_\mathrm{rot}^{-2.1}$. Magnetohydrodynamic models with large density fluctuations of two orders of magnitude or more are necessary to reproduce this gradient. △ Less

Submitted 11 October, 2021; v1 submitted 10 August, 2021; originally announced August 2021.

Comments: Accepted for publication in ApJ, 27 pages, 13 figures, 6 tables, Minor changes after proof (v3)