Showing 1–50 of 97 results for author: Oh, B

Inferring Cosmological Parameters on SDSS via Domain-Generalized Neural Networks and Lightcone Simulations

Authors: Jun-Young Lee, Ji-hoon Kim, Minyong Jung, Boon Kiat Oh, Yongseok Jo, Songyoun Park, Jaehyun Lee, Yuan-Sen Ting, Ho Seong Hwang

Abstract: We present a proof-of-concept simulation-based inference on $Ω_{\rm m}$ and $σ_{8}$ from the SDSS BOSS LOWZ NGC catalog using neural networks and domain generalization techniques without the need of summary statistics. Using rapid lightcone simulations, ${\rm L{\scriptsize -PICOLA}}$, mock galaxy catalogs are produced that fully incorporate the observational effects. The collection of galaxies is… ▽ More We present a proof-of-concept simulation-based inference on $Ω_{\rm m}$ and $σ_{8}$ from the SDSS BOSS LOWZ NGC catalog using neural networks and domain generalization techniques without the need of summary statistics. Using rapid lightcone simulations, ${\rm L{\scriptsize -PICOLA}}$, mock galaxy catalogs are produced that fully incorporate the observational effects. The collection of galaxies is fed as input to a point cloud-based network, ${\texttt{Minkowski-PointNet}}$. We also add relatively more accurate ${\rm G{\scriptsize ADGET}}$ mocks to obtain robust and generalizable neural networks. By explicitly learning the representations which reduces the discrepancies between the two different datasets via the semantic alignment loss term, we show that the latent space configuration aligns into a single plane in which the two cosmological parameters form clear axes. Consequently, during inference, the SDSS BOSS LOWZ NGC catalog maps onto the plane, demonstrating effective generalization and improving prediction accuracy compared to non-generalized models. Results from the ensemble of 25 independently trained machines find $Ω_{\rm m}=0.339 \pm 0.056$ and $σ_{8}=0.801 \pm 0.061$, inferred only from the distribution of galaxies in the lightcone slices without relying on any indirect summary statistics. A single machine that best adapts to the ${\rm G{\scriptsize ADGET}}$ mocks yields a tighter prediction of $Ω_{\rm m}=0.282 \pm 0.014$ and $σ_{8}=0.786 \pm 0.036$. We emphasize that adaptation across multiple domains can enhance the robustness of the neural networks in observational data. △ Less

Submitted 3 September, 2024; originally announced September 2024.

Comments: 28 pages, 14 figures; Accepted for publication in ApJ

DUNE Phase II: Scientific Opportunities, Detector Concepts, Technological Solutions

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, M. Adamowski, D. Adams, M. Adinolfi, C. Adriano, A. Aduszkiewicz, J. Aguilar, F. Akbar, K. Allison, S. Alonso Monsalve, M. Alrashed, A. Alton, R. Alvarez, T. Alves, H. Amar, P. Amedo, J. Anderson, C. Andreopoulos, M. Andreotti , et al. (1347 additional authors not shown)

Abstract: The international collaboration designing and constructing the Deep Underground Neutrino Experiment (DUNE) at the Long-Baseline Neutrino Facility (LBNF) has developed a two-phase strategy toward the implementation of this leading-edge, large-scale science project. The 2023 report of the US Particle Physics Project Prioritization Panel (P5) reaffirmed this vision and strongly endorsed DUNE Phase I… ▽ More The international collaboration designing and constructing the Deep Underground Neutrino Experiment (DUNE) at the Long-Baseline Neutrino Facility (LBNF) has developed a two-phase strategy toward the implementation of this leading-edge, large-scale science project. The 2023 report of the US Particle Physics Project Prioritization Panel (P5) reaffirmed this vision and strongly endorsed DUNE Phase I and Phase II, as did the European Strategy for Particle Physics. While the construction of the DUNE Phase I is well underway, this White Paper focuses on DUNE Phase II planning. DUNE Phase-II consists of a third and fourth far detector (FD) module, an upgraded near detector complex, and an enhanced 2.1 MW beam. The fourth FD module is conceived as a "Module of Opportunity", aimed at expanding the physics opportunities, in addition to supporting the core DUNE science program, with more advanced technologies. This document highlights the increased science opportunities offered by the DUNE Phase II near and far detectors, including long-baseline neutrino oscillation physics, neutrino astrophysics, and physics beyond the standard model. It describes the DUNE Phase II near and far detector technologies and detector design concepts that are currently under consideration. A summary of key R&D goals and prototyping phases needed to realize the Phase II detector technical designs is also provided. DUNE's Phase II detectors, along with the increased beam power, will complete the full scope of DUNE, enabling a multi-decadal program of groundbreaking science with neutrinos. △ Less

Submitted 22 August, 2024; originally announced August 2024.

Report number: FERMILAB-TM-2833-LBNF

Demonstrating a Bi-directional Asymmetric Frequency Conversion in Nonlinear Phononic Crystals

Authors: Yeongtae Jang, Beomseok Oh, Eunho Kim, Junsuk Rho

Abstract: Beyond the constraints of conservative systems, altering wave propagation frequency emerges as a crucial factor across diverse physical domains. This Letter demonstrates bi-directional asymmetric frequency conversion -- either upward or downward -- depending on the excitation direction in the elastic domain, moving beyond uni-directional approaches. We numerically and experimentally demonstrate it… ▽ More Beyond the constraints of conservative systems, altering wave propagation frequency emerges as a crucial factor across diverse physical domains. This Letter demonstrates bi-directional asymmetric frequency conversion -- either upward or downward -- depending on the excitation direction in the elastic domain, moving beyond uni-directional approaches. We numerically and experimentally demonstrate its practical realization in a model system of cylindrical beam crystals, a type of granular crystal characterized by intrinsic local resonance. This novel wave transport mechanism operates through the interplay of nonlinear contact, spatial asymmetry, and the coupling of local resonance. Thanks to the proposed highly tunable architecture, we demonstrate various ways to manipulate wave transport, including tunable frequency conversion. Given that the local resonance we employ exemplifies avoided crossings (i.e., a strong coupling effect), our work may inspire investigations into diverse physical nonlinear domains that support material/structural resonance. △ Less

Submitted 17 August, 2024; originally announced August 2024.

First Measurement of the Total Inelastic Cross-Section of Positively-Charged Kaons on Argon at Energies Between 5.0 and 7.5 GeV

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, M. Adamowski, D. Adams, M. Adinolfi, C. Adriano, A. Aduszkiewicz, J. Aguilar, F. Akbar, K. Allison, S. Alonso Monsalve, M. Alrashed, A. Alton, R. Alvarez, T. Alves, H. Amar, P. Amedo, J. Anderson, C. Andreopoulos, M. Andreotti , et al. (1341 additional authors not shown)

Abstract: ProtoDUNE Single-Phase (ProtoDUNE-SP) is a 770-ton liquid argon time projection chamber that operated in a hadron test beam at the CERN Neutrino Platform in 2018. We present a measurement of the total inelastic cross section of charged kaons on argon as a function of kaon energy using 6 and 7 GeV/$c$ beam momentum settings. The flux-weighted average of the extracted inelastic cross section at each… ▽ More ProtoDUNE Single-Phase (ProtoDUNE-SP) is a 770-ton liquid argon time projection chamber that operated in a hadron test beam at the CERN Neutrino Platform in 2018. We present a measurement of the total inelastic cross section of charged kaons on argon as a function of kaon energy using 6 and 7 GeV/$c$ beam momentum settings. The flux-weighted average of the extracted inelastic cross section at each beam momentum setting was measured to be 380$\pm$26 mbarns for the 6 GeV/$c$ setting and 379$\pm$35 mbarns for the 7 GeV/$c$ setting. △ Less

Submitted 1 August, 2024; originally announced August 2024.

Report number: CERN-EP-2024-211, FERMILAB-PUB-24-0216-V

Supernova Pointing Capabilities of DUNE

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, M. Adamowski, D. Adams, M. Adinolfi, C. Adriano, A. Aduszkiewicz, J. Aguilar, B. Aimard, F. Akbar, K. Allison, S. Alonso Monsalve, M. Alrashed, A. Alton, R. Alvarez, T. Alves, H. Amar, P. Amedo, J. Anderson, D. A. Andrade , et al. (1340 additional authors not shown)

Abstract: The determination of the direction of a stellar core collapse via its neutrino emission is crucial for the identification of the progenitor for a multimessenger follow-up. A highly effective method of reconstructing supernova directions within the Deep Underground Neutrino Experiment (DUNE) is introduced. The supernova neutrino pointing resolution is studied by simulating and reconstructing electr… ▽ More The determination of the direction of a stellar core collapse via its neutrino emission is crucial for the identification of the progenitor for a multimessenger follow-up. A highly effective method of reconstructing supernova directions within the Deep Underground Neutrino Experiment (DUNE) is introduced. The supernova neutrino pointing resolution is studied by simulating and reconstructing electron-neutrino charged-current absorption on $^{40}$Ar and elastic scattering of neutrinos on electrons. Procedures to reconstruct individual interactions, including a newly developed technique called ``brems flipping'', as well as the burst direction from an ensemble of interactions are described. Performance of the burst direction reconstruction is evaluated for supernovae happening at a distance of 10 kpc for a specific supernova burst flux model. The pointing resolution is found to be 3.4 degrees at 68% coverage for a perfect interaction-channel classification and a fiducial mass of 40 kton, and 6.6 degrees for a 10 kton fiducial mass respectively. Assuming a 4% rate of charged-current interactions being misidentified as elastic scattering, DUNE's burst pointing resolution is found to be 4.3 degrees (8.7 degrees) at 68% coverage. △ Less

Submitted 14 July, 2024; originally announced July 2024.

Comments: 25 pages, 16 figures

Report number: FERMILAB-PUB-24-0319-LBNF

Leading Whitespaces of Language Models' Subword Vocabulary Poses a Confound for Calculating Word Probabilities

Authors: Byung-Doh Oh, William Schuler

Abstract: Word-by-word conditional probabilities from Transformer-based language models are increasingly being used to evaluate their predictions over minimal pairs or to model the incremental processing difficulty of human readers. In this paper, we argue that there is a confound posed by the subword tokenization scheme of such language models, which has gone unaddressed thus far. This is due to the fact t… ▽ More Word-by-word conditional probabilities from Transformer-based language models are increasingly being used to evaluate their predictions over minimal pairs or to model the incremental processing difficulty of human readers. In this paper, we argue that there is a confound posed by the subword tokenization scheme of such language models, which has gone unaddressed thus far. This is due to the fact that tokens in the subword vocabulary of most language models have leading whitespaces and therefore do not naturally define stop probabilities of words. We first prove that this can result in word probabilities that sum to more than one, thereby violating the axiom that $\mathsf{P}(Ω) = 1$. This property results in a misallocation of word-by-word surprisal, where the unacceptability of the current 'end of word' is incorrectly carried over to the next word. Additionally, language models' such implicit prediction of word boundaries is incongruous with psycholinguistic experiments where human subjects directly observe upcoming word boundaries. We present a simple decoding technique to reaccount the probability of the trailing whitespace into that of the current word, which resolves this confound. As a case study, we show that this results in significantly different estimates of garden-path effects in transitive/intransitive sentences, where a comma is strongly expected before the critical word. △ Less

Submitted 16 June, 2024; originally announced June 2024.

HyperCLOVA X Technical Report

Authors: Kang Min Yoo, Jaegeun Han, Sookyo In, Heewon Jeon, Jisu Jeong, Jaewook Kang, Hyunwook Kim, Kyung-Min Kim, Munhyong Kim, Sungju Kim, Donghyun Kwak, Hanock Kwak, Se Jung Kwon, Bado Lee, Dongsoo Lee, Gichang Lee, Jooho Lee, Baeseong Park, Seongjin Shin, Joonsang Yu, Seolki Baek, Sumin Byeon, Eungsup Cho, Dooseok Choe, Jeesung Han , et al. (371 additional authors not shown)

Abstract: We introduce HyperCLOVA X, a family of large language models (LLMs) tailored to the Korean language and culture, along with competitive capabilities in English, math, and coding. HyperCLOVA X was trained on a balanced mix of Korean, English, and code data, followed by instruction-tuning with high-quality human-annotated datasets while abiding by strict safety guidelines reflecting our commitment t… ▽ More We introduce HyperCLOVA X, a family of large language models (LLMs) tailored to the Korean language and culture, along with competitive capabilities in English, math, and coding. HyperCLOVA X was trained on a balanced mix of Korean, English, and code data, followed by instruction-tuning with high-quality human-annotated datasets while abiding by strict safety guidelines reflecting our commitment to responsible AI. The model is evaluated across various benchmarks, including comprehensive reasoning, knowledge, commonsense, factuality, coding, math, chatting, instruction-following, and harmlessness, in both Korean and English. HyperCLOVA X exhibits strong reasoning capabilities in Korean backed by a deep understanding of the language and cultural nuances. Further analysis of the inherent bilingual nature and its extension to multilingualism highlights the model's cross-lingual proficiency and strong generalization ability to untargeted languages, including machine translation between several language pairs and cross-lingual inference tasks. We believe that HyperCLOVA X can provide helpful guidance for regions or countries in developing their sovereign LLMs. △ Less

Submitted 13 April, 2024; v1 submitted 2 April, 2024; originally announced April 2024.

Comments: 44 pages; updated authors list and fixed author names

Performance of a modular ton-scale pixel-readout liquid argon time projection chamber

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, M. Adamowski, D. Adams, M. Adinolfi, C. Adriano, A. Aduszkiewicz, J. Aguilar, B. Aimard, F. Akbar, K. Allison, S. Alonso Monsalve, M. Alrashed, A. Alton, R. Alvarez, T. Alves, H. Amar, P. Amedo, J. Anderson, D. A. Andrade , et al. (1340 additional authors not shown)

Abstract: The Module-0 Demonstrator is a single-phase 600 kg liquid argon time projection chamber operated as a prototype for the DUNE liquid argon near detector. Based on the ArgonCube design concept, Module-0 features a novel 80k-channel pixelated charge readout and advanced high-coverage photon detection system. In this paper, we present an analysis of an eight-day data set consisting of 25 million cosmi… ▽ More The Module-0 Demonstrator is a single-phase 600 kg liquid argon time projection chamber operated as a prototype for the DUNE liquid argon near detector. Based on the ArgonCube design concept, Module-0 features a novel 80k-channel pixelated charge readout and advanced high-coverage photon detection system. In this paper, we present an analysis of an eight-day data set consisting of 25 million cosmic ray events collected in the spring of 2021. We use this sample to demonstrate the imaging performance of the charge and light readout systems as well as the signal correlations between the two. We also report argon purity and detector uniformity measurements, and provide comparisons to detector simulations. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 47 pages, 41 figures

Report number: FERMILAB-PUB-24-0073-LBNF

The AGORA High-resolution Galaxy Simulations Comparison Project IV: Halo and Galaxy Mass Assembly in a Cosmological Zoom-in Simulation at $z\le2$

Authors: Santi Roca-Fàbrega, Ji-hoon Kim, Joel R. Primack, Minyong Jung, Anna Genina, Loic Hausammann, Hyeonyong Kim, Alessandro Lupi, Kentaro Nagamine, Johnny W. Powell, Yves Revaz, Ikkoh Shimizu, Clayton Strawn, Héctor Velázquez, Tom Abel, Daniel Ceverino, Bili Dong, Thomas R. Quinn, Eun-jin Shin, Alvaro Segovia-Otero, Oscar Agertz, Kirk S. S. Barrow, Corentin Cadiou, Avishai Dekel, Cameron Hummels , et al. (3 additional authors not shown)

Abstract: In this fourth paper from the AGORA Collaboration, we study the evolution down to redshift $z=2$ and below of a set of cosmological zoom-in simulations of a Milky Way mass galaxy by eight of the leading hydrodynamic simulation codes. We also compare this CosmoRun suite of simulations with dark matter-only simulations by the same eight codes. We analyze general properties of the halo and galaxy at… ▽ More In this fourth paper from the AGORA Collaboration, we study the evolution down to redshift $z=2$ and below of a set of cosmological zoom-in simulations of a Milky Way mass galaxy by eight of the leading hydrodynamic simulation codes. We also compare this CosmoRun suite of simulations with dark matter-only simulations by the same eight codes. We analyze general properties of the halo and galaxy at $z=4$ and 3, and before the last major merger, focusing on the formation of well-defined rotationally-supported disks, the mass-metallicity relation, the specific star formation rate, the gas metallicity gradients, and the non-axisymmetric structures in the stellar disks. Codes generally converge well to the stellar-to-halo mass ratios predicted by semi-analytic models at $z\sim$2. We see that almost all the hydro codes develop rotationally-supported structures at low redshifts. Most agree within 0.5 dex with the observed MZR at high and intermediate redshifts, and reproduce the gas metallicity gradients obtained from analytical models and low-redshift observations. We confirm that the inter-code differences in the halo assembly history reported in the first paper of the collaboration also exist in CosmoRun, making the code-to-code comparison more difficult. We show that such differences are mainly due to variations in code-dependent parameters that control the time-stepping strategy of the gravity solver. We find that variations in the early stellar feedback can also result in differences in the timing of the low-redshift mergers. All the simulation data down to $z=2$ and the auxiliary data will be made publicly available. △ Less

Submitted 9 February, 2024; originally announced February 2024.

Comments: Submitted to ApJ. Visit the AGORA Collaboration website (www.agorasimulations.org <http://www.agorasimulations.org/>) for more information

The AGORA High-resolution Galaxy Simulations Comparison Project. V: Satellite Galaxy Populations In A Cosmological Zoom-in Simulation of A Milky Way-mass Halo

Authors: Minyong Jung, Santi Roca-Fàbrega, Ji-hoon Kim, Anna Genina, Loic Hausammann, Hyeonyong Kim, Alessandro Lupi, Kentaro Nagamine, Johnny W. Powell, Yves Revaz, Ikkoh Shimizu, Héctor Velázquez, Daniel Ceverino, Joel R. Primack, Thomas R. Quinn, Clayton Strawn, Tom Abel, Avishai Dekel, Bili Dong, Boon Kiat Oh, Romain Teyssier

Abstract: We analyze and compare the satellite halo populations at $z\sim2$ in the high-resolution cosmological zoom-in simulations of a $10^{12}\,{\rm M}_{\odot}$ target halo ($z=0$ mass) carried out on eight widely-used astrophysical simulation codes ({\sc Art-I}, {\sc Enzo}, {\sc Ramses}, {\sc Changa}, {\sc Gadget-3}, {\sc Gear}, {\sc Arepo-t}, and {\sc Gizmo}) for the {\it AGORA} High-resolution Galaxy… ▽ More We analyze and compare the satellite halo populations at $z\sim2$ in the high-resolution cosmological zoom-in simulations of a $10^{12}\,{\rm M}_{\odot}$ target halo ($z=0$ mass) carried out on eight widely-used astrophysical simulation codes ({\sc Art-I}, {\sc Enzo}, {\sc Ramses}, {\sc Changa}, {\sc Gadget-3}, {\sc Gear}, {\sc Arepo-t}, and {\sc Gizmo}) for the {\it AGORA} High-resolution Galaxy Simulations Comparison Project. We use slightly different redshift epochs near $z=2$ for each code (hereafter ``$z\sim2$') at which the eight simulations are in the same stage in the target halo's merger history. After identifying the matched pairs of halos between the {\it CosmoRun} simulations and the DMO simulations, we discover that each {\it CosmoRun} halo tends to be less massive than its DMO counterpart. When we consider only the halos containing stellar particles at $z\sim2$, the number of satellite {\it galaxies} is significantly fewer than that of dark matter halos in all participating {\it AGORA} simulations, and is comparable to the number of present-day satellites near the Milky Way or M31. The so-called ``missing satellite problem' is fully resolved across all participating codes simply by implementing the common baryonic physics adopted in {\it AGORA} and the stellar feedback prescription commonly used in each code, with sufficient numerical resolution ($\lesssim100$ proper pc at $z=2$). We also compare other properties such as the stellar mass$-$halo mass relation and the mass$-$metallicity relation. Our work highlights the value of comparison studies such as {\it AGORA}, where outstanding problems in galaxy formation theory are studied simultaneously on multiple numerical platforms. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: Accepted for publication to the ApJ, 19 pages, 10 figures. Visit the AGORA Collaboration website (www.agorasimulations.org) for more information. For summary video, please see https://www.youtube.com/watch?v=HShRJYvPe1k

Journal ref: ApJ 964 123 (2024)

The AGORA High-resolution Galaxy Simulations Comparison Project. VI. Similarities and Differences in the Circumgalactic Medium

Authors: Clayton Strawn, Santi Roca-Fàbrega, Joel R. Primack, Ji-hoon Kim, Anna Genina, Loic Hausammann, Hyeonyong Kim, Alessandro Lupi, Kentaro Nagamine, Johnny W. Powell, Yves Revaz, Ikkoh Shimizu, Héctor Velázquez, Tom Abel, Daniel Ceverino, Bili Dong, Minyong Jung, Thomas R. Quinn, Eun-jin Shin, Kirk S. S. Barrow, Avishai Dekel, Boon Kiat Oh, Nir Mandelker, Romain Teyssier, Cameron Hummels , et al. (4 additional authors not shown)

Abstract: We analyze the circumgalactic medium (CGM) for eight commonly-used cosmological codes in the AGORA collaboration. The codes are calibrated to use identical initial conditions, cosmology, heating and cooling, and star formation thresholds, but each evolves with its own unique code architecture and stellar feedback implementation. Here, we analyze the results of these simulations in terms of the str… ▽ More We analyze the circumgalactic medium (CGM) for eight commonly-used cosmological codes in the AGORA collaboration. The codes are calibrated to use identical initial conditions, cosmology, heating and cooling, and star formation thresholds, but each evolves with its own unique code architecture and stellar feedback implementation. Here, we analyze the results of these simulations in terms of the structure, composition, and phase dynamics of the CGM. We show properties such as metal distribution, ionization levels, and kinematics are effective tracers of the effects of the different code feedback and implementation methods, and as such they can be highly divergent between simulations. This is merely a fiducial set of models, against which we will in the future compare multiple feedback recipes for each code. Nevertheless, we find that the large parameter space these simulations establish can help disentangle the different variables that affect observable quantities in the CGM, e.g., showing that abundances for ions with higher ionization energy are more strongly determined by the simulation's metallicity, while abundances for ions with lower ionization energy are more strongly determined by the gas density and temperature. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: 32 pages, 17 figures. Published in ApJ

Journal ref: ApJ 962 29 (2024)

Frequency Explains the Inverse Correlation of Large Language Models' Size, Training Data Amount, and Surprisal's Fit to Reading Times

Authors: Byung-Doh Oh, Shisen Yue, William Schuler

Abstract: Recent studies have shown that as Transformer-based language models become larger and are trained on very large amounts of data, the fit of their surprisal estimates to naturalistic human reading times degrades. The current work presents a series of analyses showing that word frequency is a key explanatory factor underlying these two trends. First, residual errors from four language model families… ▽ More Recent studies have shown that as Transformer-based language models become larger and are trained on very large amounts of data, the fit of their surprisal estimates to naturalistic human reading times degrades. The current work presents a series of analyses showing that word frequency is a key explanatory factor underlying these two trends. First, residual errors from four language model families on four corpora show that the inverse correlation between model size and fit to reading times is the strongest on the subset of least frequent words, which is driven by excessively accurate predictions of larger model variants. Additionally, training dynamics reveal that during later training steps, all model variants learn to predict rare words and that larger model variants do so more accurately, which explains the detrimental effect of both training data amount and model size on fit to reading times. Finally, a feature attribution analysis demonstrates that larger model variants are able to accurately predict rare words based on both an effectively longer context window size as well as stronger local associations compared to smaller model variants. Taken together, these results indicate that Transformer-based language models' surprisal estimates diverge from human-like expectations due to the superhumanly complex associations they learn for predicting rare words. △ Less

Submitted 3 February, 2024; originally announced February 2024.

Comments: EACL 2024

Doping Liquid Argon with Xenon in ProtoDUNE Single-Phase: Effects on Scintillation Light

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, M. Adamowski, D. Adams, M. Adinolfi, C. Adriano, A. Aduszkiewicz, J. Aguilar, B. Aimard, F. Akbar, K. Allison, S. Alonso Monsalve, M. Alrashed, A. Alton, R. Alvarez, H. Amar Es-sghir, P. Amedo, J. Anderson, D. A. Andrade, C. Andreopoulos , et al. (1297 additional authors not shown)

Abstract: Doping of liquid argon TPCs (LArTPCs) with a small concentration of xenon is a technique for light-shifting and facilitates the detection of the liquid argon scintillation light. In this paper, we present the results of the first doping test ever performed in a kiloton-scale LArTPC. From February to May 2020, we carried out this special run in the single-phase DUNE Far Detector prototype (ProtoDUN… ▽ More Doping of liquid argon TPCs (LArTPCs) with a small concentration of xenon is a technique for light-shifting and facilitates the detection of the liquid argon scintillation light. In this paper, we present the results of the first doping test ever performed in a kiloton-scale LArTPC. From February to May 2020, we carried out this special run in the single-phase DUNE Far Detector prototype (ProtoDUNE-SP) at CERN, featuring 720 t of total liquid argon mass with 410 t of fiducial mass. A 5.4 ppm nitrogen contamination was present during the xenon doping campaign. The goal of the run was to measure the light and charge response of the detector to the addition of xenon, up to a concentration of 18.8 ppm. The main purpose was to test the possibility for reduction of non-uniformities in light collection, caused by deployment of photon detectors only within the anode planes. Light collection was analysed as a function of the xenon concentration, by using the pre-existing photon detection system (PDS) of ProtoDUNE-SP and an additional smaller set-up installed specifically for this run. In this paper we first summarize our current understanding of the argon-xenon energy transfer process and the impact of the presence of nitrogen in argon with and without xenon dopant. We then describe the key elements of ProtoDUNE-SP and the injection method deployed. Two dedicated photon detectors were able to collect the light produced by xenon and the total light. The ratio of these components was measured to be about 0.65 as 18.8 ppm of xenon were injected. We performed studies of the collection efficiency as a function of the distance between tracks and light detectors, demonstrating enhanced uniformity of response for the anode-mounted PDS. We also show that xenon doping can substantially recover light losses due to contamination of the liquid argon by nitrogen. △ Less

Submitted 2 August, 2024; v1 submitted 2 February, 2024; originally announced February 2024.

Comments: 36 pages, 20 figures. Corrected author list; corrected typos across paper and polished text

Report number: CERN-EP-2024-024; FERMILAB-PUB-23-0819-LBNF

New Beam Dynamics Code for Cyclotron Analysis

Authors: G-H. Kim, H-J. Cho, B-H. Oh, G-R. Hahn, M. Chung, S. Park, S. Shin

Abstract: This paper describes the beam dynamic simulation with transfer matrix method for cyclotron. Starting from a description on the equation of motion in the cyclotron, lattice functions were determined from transfer matrix method and the solutions for the 2nd-order nonlinear Hamiltonian were introduced and used in phase space particle tracking. Based on the description of beam dynamics in the cyclotro… ▽ More This paper describes the beam dynamic simulation with transfer matrix method for cyclotron. Starting from a description on the equation of motion in the cyclotron, lattice functions were determined from transfer matrix method and the solutions for the 2nd-order nonlinear Hamiltonian were introduced and used in phase space particle tracking. Based on the description of beam dynamics in the cyclotron, simulation code was also developed for cyclotron design. △ Less

Submitted 19 January, 2024; originally announced January 2024.

Comments: arXiv admin note: substantial text overlap with arXiv:1807.01397

The DUNE Far Detector Vertical Drift Technology, Technical Design Report

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, M. Adamowski, D. Adams, M. Adinolfi, C. Adriano, A. Aduszkiewicz, J. Aguilar, B. Aimard, F. Akbar, K. Allison, S. Alonso Monsalve, M. Alrashed, A. Alton, R. Alvarez, H. Amar, P. Amedo, J. Anderson, D. A. Andrade, C. Andreopoulos , et al. (1304 additional authors not shown)

Abstract: DUNE is an international experiment dedicated to addressing some of the questions at the forefront of particle physics and astrophysics, including the mystifying preponderance of matter over antimatter in the early universe. The dual-site experiment will employ an intense neutrino beam focused on a near and a far detector as it aims to determine the neutrino mass hierarchy and to make high-precisi… ▽ More DUNE is an international experiment dedicated to addressing some of the questions at the forefront of particle physics and astrophysics, including the mystifying preponderance of matter over antimatter in the early universe. The dual-site experiment will employ an intense neutrino beam focused on a near and a far detector as it aims to determine the neutrino mass hierarchy and to make high-precision measurements of the PMNS matrix parameters, including the CP-violating phase. It will also stand ready to observe supernova neutrino bursts, and seeks to observe nucleon decay as a signature of a grand unified theory underlying the standard model. The DUNE far detector implements liquid argon time-projection chamber (LArTPC) technology, and combines the many tens-of-kiloton fiducial mass necessary for rare event searches with the sub-centimeter spatial resolution required to image those events with high precision. The addition of a photon detection system enhances physics capabilities for all DUNE physics drivers and opens prospects for further physics explorations. Given its size, the far detector will be implemented as a set of modules, with LArTPC designs that differ from one another as newer technologies arise. In the vertical drift LArTPC design, a horizontal cathode bisects the detector, creating two stacked drift volumes in which ionization charges drift towards anodes at either the top or bottom. The anodes are composed of perforated PCB layers with conductive strips, enabling reconstruction in 3D. Light-trap-style photon detection modules are placed both on the cryostat's side walls and on the central cathode where they are optically powered. This Technical Design Report describes in detail the technical implementations of each subsystem of this LArTPC that, together with the other far detector modules and the near detector, will enable DUNE to achieve its physics goals. △ Less

Submitted 5 December, 2023; originally announced December 2023.

Comments: 425 pages; 281 figures Central editing team: A. Heavey, S. Kettell, A. Marchionni, S. Palestini, S. Rajogopalan, R. J. Wilson

Report number: Fermilab Report no: TM-2813-LBNF

Merger Tree-based Galaxy Matching: A Comparative Study Across Different Resolutions

Authors: Minyong Jung, Ji-hoon Kim, Boon Kiat Oh, Sungwook E. Hong, Jaehyun Lee, Juhan Kim

Abstract: We introduce a novel halo/galaxy matching technique between two cosmological simulations with different resolutions, which utilizes the positions and masses of halos along their subhalo merger tree. With this tool, we conduct a study of resolution biases through the {\it galaxy-by-galaxy} inspection of a pair of simulations that have the same simulation configuration but different mass resolutions… ▽ More We introduce a novel halo/galaxy matching technique between two cosmological simulations with different resolutions, which utilizes the positions and masses of halos along their subhalo merger tree. With this tool, we conduct a study of resolution biases through the {\it galaxy-by-galaxy} inspection of a pair of simulations that have the same simulation configuration but different mass resolutions, utilizing a suite of {\sc IllustrisTNG} simulations to assess the impact on galaxy properties. We find that, with the subgrid physics model calibrated for TNG100-1, subhalos in TNG100-1 (high resolution) have $\lesssim0.5$ dex higher stellar masses than their counterparts in the TNG100-2 (low-resolution). It is also discovered that the subhalos with $M_{\mathrm{gas}}\sim10^{8.5}\,{\rm M}_\odot$ in TNG100-1 have $\sim0.5$ dex higher gas mass than those in TNG100-2. The mass profiles of the subhalos reveal that the dark matter masses of subhalos in TNG100-2 converge well with those from TNG100-1, except within 4 kpc of the resolution limit. The differences in stellar mass and hot gas mass are most pronounced in the central region. We exploit machine learning to build a correction mapping for the physical quantities of subhalos from low- to high-resolution simulations (TNG300-1 and TNG100-1), which enables us to find an efficient way to compile a high-resolution galaxy catalog even from a low-resolution simulation. Our tools can easily be applied to other large cosmological simulations, testing and mitigating the resolution biases of their numerical codes and subgrid physics models. △ Less

Submitted 19 March, 2024; v1 submitted 4 December, 2023; originally announced December 2023.

Comments: 24 pages, 14 figures; Accepted for publication in to ApJ

Journal ref: ApJ 965 156 (2024)

Dense stellar clump formation driven by strong quasar winds in the FIRE cosmological hydrodynamic simulations

Authors: Jonathan Mercedes-Feliz, Daniel Anglés-Alcázar, Boon Kiat Oh, Christopher C. Hayward, Rachel K. Cochrane, Alexander J. Richings, Claude-André Faucher-Giguère, Sarah Wellons, Bryan A. Terrazas, Jorge Moreno, Kung Yi Su, Philip F. Hopkins

Abstract: We investigate the formation of dense stellar clumps in a suite of high-resolution cosmological zoom-in simulations of a massive, star forming galaxy at $z \sim 2$ under the presence of strong quasar winds. Our simulations include multi-phase ISM physics from the Feedback In Realistic Environments (FIRE) project and a novel implementation of hyper-refined accretion disk winds. We show that powerfu… ▽ More We investigate the formation of dense stellar clumps in a suite of high-resolution cosmological zoom-in simulations of a massive, star forming galaxy at $z \sim 2$ under the presence of strong quasar winds. Our simulations include multi-phase ISM physics from the Feedback In Realistic Environments (FIRE) project and a novel implementation of hyper-refined accretion disk winds. We show that powerful quasar winds can have a global negative impact on galaxy growth while in the strongest cases triggering the formation of an off-center clump with stellar mass ${\rm M}_{\star}\sim 10^{7}\,{\rm M}_{\odot}$, effective radius ${\rm R}_{\rm 1/2\,\rm Clump}\sim 20\,{\rm pc}$, and surface density $Σ_{\star} \sim 10^{4}\,{\rm M}_{\odot}\,{\rm pc}^{-2}$. The clump progenitor gas cloud is originally not star-forming, but strong ram pressure gradients driven by the quasar winds (orders of magnitude stronger than experienced in the absence of winds) lead to rapid compression and subsequent conversion of gas into stars at densities much higher than the average density of star-forming gas. The AGN-triggered star-forming clump reaches ${\rm SFR} \sim 50\,{\rm M}_{\odot}\,{\rm yr}^{-1}$ and $Σ_{\rm SFR} \sim 10^{4}\,{\rm M}_{\odot}\,{\rm yr}^{-1}\,{\rm kpc}^{-2}$, converting most of the progenitor gas cloud into stars in $\sim$2\,Myr, significantly faster than its initial free-fall time and with stellar feedback unable to stop star formation. In contrast, the same gas cloud in the absence of quasar winds forms stars over a much longer period of time ($\sim$35\,Myr), at lower densities, and losing spatial coherency. The presence of young, ultra-dense, gravitationally bound stellar clumps in recently quenched galaxies could thus indicate local positive feedback acting alongside the strong negative impact of powerful quasar winds, providing a plausible formation scenario for globular clusters. △ Less

Submitted 17 April, 2024; v1 submitted 30 October, 2023; originally announced October 2023.

Comments: 14 pages, 12 figures. Accepted for publication in MNRAS; doi:10.1093/mnras/stae1021

Primitively $2$-universal senary integral quadratic forms

Authors: Byeong-Kweon Oh, Jongheun Yoon

Abstract: For a positive integer $m$, a (positive definite integral) quadratic form is called primitively $m$-universal if it primitively represents all quadratic forms of rank $m$. It was proved in arXiv:2202.13573 that there are exactly $107$ equivalence classes of primitively $1$-universal quaternary quadratic forms. In this article, we prove that the minimal rank of primitively $2$-universal quadratic f… ▽ More For a positive integer $m$, a (positive definite integral) quadratic form is called primitively $m$-universal if it primitively represents all quadratic forms of rank $m$. It was proved in arXiv:2202.13573 that there are exactly $107$ equivalence classes of primitively $1$-universal quaternary quadratic forms. In this article, we prove that the minimal rank of primitively $2$-universal quadratic forms is six, and there are exactly $201$ equivalence classes of primitively $2$-universal senary quadratic forms. △ Less

Submitted 2 September, 2023; originally announced September 2023.

Comments: 35 pages

MSC Class: 11E12; 11E20

Isolations of the sum of two squares from its proper subforms

Authors: Jangwon Ju, Daejun Kim, Kyoungmin Kim, Mingyu Kim, Byeong-Kweon Oh

Abstract: For a (positive definite and integral) quadratic form $f$, a quadratic form is said to be {\it an isolation of $f$ from its proper subforms} if it represents all proper subforms of $f$, but not $f$ itself. It was proved that the minimal rank of isolations of the square quadratic form $x^2$ is three, and there are exactly $15$ ternary diagonal isolations of $x^2$. Recently, it was proved that any q… ▽ More For a (positive definite and integral) quadratic form $f$, a quadratic form is said to be {\it an isolation of $f$ from its proper subforms} if it represents all proper subforms of $f$, but not $f$ itself. It was proved that the minimal rank of isolations of the square quadratic form $x^2$ is three, and there are exactly $15$ ternary diagonal isolations of $x^2$. Recently, it was proved that any quaternary quadratic form cannot be an isolation of the sum of two squares $I_2=x^2+y^2$, and there are quinary isolations of $I_2$. In this article, we prove that there are at most $231$ quinary isolations of $I_2$, which are listed in Table $1$. Moreover, we prove that $14$ quinary quadratic forms with dagger mark in Table $1$ are isolations of $I_2$. △ Less

Submitted 9 August, 2023; originally announced August 2023.

Comments: 14 pages

MSC Class: 11E12; 11E20; 11E25

Token-wise Decomposition of Autoregressive Language Model Hidden States for Analyzing Model Predictions

Authors: Byung-Doh Oh, William Schuler

Abstract: While there is much recent interest in studying why Transformer-based large language models make predictions the way they do, the complex computations performed within each layer have made their behavior somewhat opaque. To mitigate this opacity, this work presents a linear decomposition of final hidden states from autoregressive language models based on each initial input token, which is exact fo… ▽ More While there is much recent interest in studying why Transformer-based large language models make predictions the way they do, the complex computations performed within each layer have made their behavior somewhat opaque. To mitigate this opacity, this work presents a linear decomposition of final hidden states from autoregressive language models based on each initial input token, which is exact for virtually all contemporary Transformer architectures. This decomposition allows the definition of probability distributions that ablate the contribution of specific input tokens, which can be used to analyze their influence on model probabilities over a sequence of upcoming words with only one forward pass from the model. Using the change in next-word probability as a measure of importance, this work first examines which context words make the biggest contribution to language model predictions. Regression experiments suggest that Transformer-based language models rely primarily on collocational associations, followed by linguistic factors such as syntactic dependencies and coreference relationships in making next-word predictions. Additionally, analyses using these measures to predict syntactic dependencies and coreferent mention spans show that collocational association and repetitions of the same token largely explain the language models' predictions on these tasks. △ Less

Submitted 2 June, 2023; v1 submitted 17 May, 2023; originally announced May 2023.

Comments: ACL 2023

Transformer-Based Language Model Surprisal Predicts Human Reading Times Best with About Two Billion Training Tokens

Authors: Byung-Doh Oh, William Schuler

Abstract: Recent psycholinguistic studies have drawn conflicting conclusions about the relationship between the quality of a language model and the ability of its surprisal estimates to predict human reading times, which has been speculated to be due to the large gap in both the amount of training data and model capacity across studies. The current work aims to consolidate these findings by evaluating surpr… ▽ More Recent psycholinguistic studies have drawn conflicting conclusions about the relationship between the quality of a language model and the ability of its surprisal estimates to predict human reading times, which has been speculated to be due to the large gap in both the amount of training data and model capacity across studies. The current work aims to consolidate these findings by evaluating surprisal estimates from Transformer-based language model variants that vary systematically in the amount of training data and model capacity on their ability to predict human reading times. The results show that surprisal estimates from most variants with contemporary model capacities provide the best fit after seeing about two billion training tokens, after which they begin to diverge from humanlike expectations. Additionally, newly-trained smaller model variants reveal a 'tipping point' at convergence, after which the decrease in language model perplexity begins to result in poorer fits to human reading times. These results suggest that the massive amount of training data is mainly responsible for the poorer fit achieved by surprisal from larger pre-trained language models, and that a certain degree of model capacity is necessary for Transformer-based language models to capture humanlike expectations. △ Less

Submitted 22 October, 2023; v1 submitted 22 April, 2023; originally announced April 2023.

Comments: Findings of the Association for Computational Linguistics: EMNLP 2023

Why Does Surprisal From Larger Transformer-Based Language Models Provide a Poorer Fit to Human Reading Times?

Authors: Byung-Doh Oh, William Schuler

Abstract: This work presents a detailed linguistic analysis into why larger Transformer-based pre-trained language models with more parameters and lower perplexity nonetheless yield surprisal estimates that are less predictive of human reading times. First, regression analyses show a strictly monotonic, positive log-linear relationship between perplexity and fit to reading times for the more recently releas… ▽ More This work presents a detailed linguistic analysis into why larger Transformer-based pre-trained language models with more parameters and lower perplexity nonetheless yield surprisal estimates that are less predictive of human reading times. First, regression analyses show a strictly monotonic, positive log-linear relationship between perplexity and fit to reading times for the more recently released five GPT-Neo variants and eight OPT variants on two separate datasets, replicating earlier results limited to just GPT-2 (Oh et al., 2022). Subsequently, analysis of residual errors reveals a systematic deviation of the larger variants, such as underpredicting reading times of named entities and making compensatory overpredictions for reading times of function words such as modals and conjunctions. These results suggest that the propensity of larger Transformer-based models to 'memorize' sequences during training makes their surprisal estimates diverge from humanlike expectations, which warrants caution in using pre-trained language models to study human language processing. △ Less

Submitted 22 December, 2022; originally announced December 2022.

Comments: Transactions of the Association for Computational Linguistics (pre-MIT Press publication version)

Entropy- and Distance-Based Predictors From GPT-2 Attention Patterns Predict Reading Times Over and Above GPT-2 Surprisal

Authors: Byung-Doh Oh, William Schuler

Abstract: Transformer-based large language models are trained to make predictions about the next word by aggregating representations of previous tokens through their self-attention mechanism. In the field of cognitive modeling, such attention patterns have recently been interpreted as embodying the process of cue-based retrieval, in which attention over multiple targets is taken to generate interference and… ▽ More Transformer-based large language models are trained to make predictions about the next word by aggregating representations of previous tokens through their self-attention mechanism. In the field of cognitive modeling, such attention patterns have recently been interpreted as embodying the process of cue-based retrieval, in which attention over multiple targets is taken to generate interference and latency during retrieval. Under this framework, this work first defines an entropy-based predictor that quantifies the diffuseness of self-attention, as well as distance-based predictors that capture the incremental change in attention patterns across timesteps. Moreover, following recent studies that question the informativeness of attention weights, we also experiment with alternative methods for incorporating vector norms into attention weights. Regression experiments using predictors calculated from the GPT-2 language model show that these predictors deliver a substantially better fit to held-out self-paced reading and eye-tracking data over a rigorous baseline including GPT-2 surprisal. Additionally, the distance-based predictors generally demonstrated higher predictive power, with effect sizes of up to 6.59 ms per standard deviation on self-paced reading times (compared to 2.82 ms for surprisal) and 1.05 ms per standard deviation on eye-gaze durations (compared to 3.81 ms for surprisal). △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: EMNLP 2022

Growth of A Massive Black Hole Via Tidal Disruption Accretion

Authors: Seungjae Lee, Ji-hoon Kim, Boon Kiat Oh

Abstract: Stars that are tidally disrupted by the massive black hole (MBH) may contribute significantly to the growth of the MBH, especially in dense nuclear star clusters (NSCs). Yet, this tidal disruption accretion (TDA) of stars onto the MBH has largely been overlooked compared to the gas accretion (GA) channel in most numerical experiments until now. In this work, we implement a black hole growth channe… ▽ More Stars that are tidally disrupted by the massive black hole (MBH) may contribute significantly to the growth of the MBH, especially in dense nuclear star clusters (NSCs). Yet, this tidal disruption accretion (TDA) of stars onto the MBH has largely been overlooked compared to the gas accretion (GA) channel in most numerical experiments until now. In this work, we implement a black hole growth channel via TDA in the high-resolution adaptive mesh refinement code Enzo to investigate its influence on a MBH seed's early evolution. We find that a MBH seed grows rapidly from $10^3\,\mathrm{M}_\odot$ to $\gtrsim 10^6\,\mathrm{M}_\odot$ in 200\,Myrs in some of the tested simulations. Compared to a MBH seed that grows only via GA, TDA can enhance the MBH's growth rate by up to more than an order of magnitude. However, as predicted, TDA mainly helps the early growth of the MBH (from $10^{3-4}\,\mathrm{M}_\odot$ to $\lesssim10^{5}\,\mathrm{M}_\odot$) while the later evolution is generally dominated by GA. We also observe that the star formation near the MBH is suppressed when TDA is most active, sometimes with a visible cavity in gas (of size $\sim$ a few pc) created in the vicinity of the MBH. It is because the MBH may grow expeditiously with both GA and TDA, and the massive MBH could consume its neighboring gas faster than being replenished by gas inflows. Our study demonstrates the need to consider different channels of black hole accretion that may provide clues for the existence of supermassive black holes at high redshifts. △ Less

Submitted 28 November, 2022; v1 submitted 4 November, 2022; originally announced November 2022.

Comments: 17 pages, and 10 figures

The fate of baryons in counterfactual universes

Authors: Boon Kiat Oh, John A. Peacock, Sadegh Khochfar, Britton D. Smith

Abstract: We present results from nine simulations that compare the standard $Λ$ Cold Dark Matter cosmology ($Λ$CDM) with counterfactual universes, for approximately $100\,{\rm Gyr}$ using the Enzo simulation code. We vary the value of $Λ$ and the fluctuation amplitude to explore the effect on the evolution of the halo mass function (HMF), the intergalactic medium (IGM) and the star formation history (SFH).… ▽ More We present results from nine simulations that compare the standard $Λ$ Cold Dark Matter cosmology ($Λ$CDM) with counterfactual universes, for approximately $100\,{\rm Gyr}$ using the Enzo simulation code. We vary the value of $Λ$ and the fluctuation amplitude to explore the effect on the evolution of the halo mass function (HMF), the intergalactic medium (IGM) and the star formation history (SFH). The distinct peak in star formation rate density (SFRD) and its subsequent decline are both affected by the interplay between gravitational attraction and the accelerating effects of $Λ$. The IGM cools down more rapidly in models with a larger $Λ$ and also with a lower $σ_8$, reflecting the reduced SFRD associated with these changes -- although changing $σ_8$ is not degenerate with changing $Λ$, either regarding the thermal history of the IGM or the SFH. However, these induced changes to the IGM or ionizing background have little impact on the calculated SFRD. We provide fits for the evolution of the SFRD in these different universes, which we integrate over time to derive an asymptotic star formation efficiency. Together with Weinberg's uniform prior on $Λ$, the estimated probability of observers experiencing a value of $Λ$ no greater than the observed value is 13%, substantially larger than some alternative estimates. Within the Enzo model framework, then, observer selection within a multiverse is able to account statistically for the small value of the cosmological constant, although $Λ$ in our universe does appear to be at the low end of the predicted range. △ Less

Submitted 19 September, 2022; originally announced September 2022.

Comments: 18 pages, 10 figures, MNRAS: MN-22-1662-MJ.R1

Machine-guided Exploration and Calibration of Astrophysical Simulations

Authors: Boon Kiat Oh, Hongjun An, Eun-jin Shin, Ji-hoon Kim, Sungwook E. Hong

Abstract: We apply a novel method with machine learning to calibrate sub-grid models within numerical simulation codes to achieve convergence with observations and between different codes. It utilizes active learning and neural density estimators. The hyper parameters of the machine are calibrated with a well-defined projectile motion problem. Then, using a set of 22 cosmological zoom simulations, we tune t… ▽ More We apply a novel method with machine learning to calibrate sub-grid models within numerical simulation codes to achieve convergence with observations and between different codes. It utilizes active learning and neural density estimators. The hyper parameters of the machine are calibrated with a well-defined projectile motion problem. Then, using a set of 22 cosmological zoom simulations, we tune the parameters of a popular star formation and feedback model within Enzo to match simulations. The parameters that are adjusted include the star formation efficiency, coupling of thermal energy from stellar feedback, and volume into which the energy is deposited. This number translates to a factor of more than three improvements over manual calibration. Despite using fewer simulations, we obtain a better agreement to the observed baryon makeup of a Milky-Way (MW) sized halo. Switching to a different strategy, we improve the consistency of the recommended parameters from the machine. Given the success of the calibration, we then apply the technique to reconcile metal transport between grid-based and particle-based simulation codes using an isolated galaxy. It is an improvement over manual exploration while hinting at a less known relation between the diffusion coefficient and the metal mass in the halo region. The exploration and calibration of the parameters of the sub-grid models with a machine learning approach is concluded to be versatile and directly applicable to different problems. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Primitively universal quaternary quadratic forms

Authors: Jangwon Ju, Daejun Kim, Kyoungmin Kim, Mingyu Kim, Byeong-Kweon Oh

Abstract: A (positive definite and integral) quadratic form $f$ is said to be $\textit{universal}$ if it represents all positive integers, and is said to be $\textit{primitively universal}$ if it represents all positive integers primitively. We also say $f$ is $\textit{primitively almost universal}$ if it represents almost all positive integers primitively. Conway and Schneeberger proved (see [1]) that ther… ▽ More A (positive definite and integral) quadratic form $f$ is said to be $\textit{universal}$ if it represents all positive integers, and is said to be $\textit{primitively universal}$ if it represents all positive integers primitively. We also say $f$ is $\textit{primitively almost universal}$ if it represents almost all positive integers primitively. Conway and Schneeberger proved (see [1]) that there are exactly $204$ equivalence classes of universal quaternary quadratic forms. Recently, Earnest and Gunawardana proved in [4] that among $204$ equivalence classes of universal quaternary quadratic forms, there are exactly $152$ equivalence classes of primitively almost universal quaternary quadratic forms. In this article, we prove that there are exactly $107$ equivalence classes of primitively universal quaternary quadratic forms. We also determine the set of all positive integers that are not primitively represented by each of the remaining $152-107=45$ equivalence classes of primitively almost universal quaternary quadratic forms. △ Less

Submitted 28 February, 2022; originally announced February 2022.

Comments: 21 pages

MSC Class: 11E12; 11E20

Can we recover an integral quadratic form by representing all its subforms?

Authors: Wai Kiu Chan, Byeong-Kweon Oh

Abstract: Let $\mathfrak o$ be the ring of integers of a totally real number field. If $f$ is a quadratic form over $\mathfrak o$ and $g$ is another quadratic form over $\mathfrak o$ which represents all proper subforms of $f$, does $g$ represent $f$? We show that if $g$ is indefinite, then $g$ indeed represents $f$. However, when $f$ is positive definite and indecomposable, then there exists a $g$ which re… ▽ More Let $\mathfrak o$ be the ring of integers of a totally real number field. If $f$ is a quadratic form over $\mathfrak o$ and $g$ is another quadratic form over $\mathfrak o$ which represents all proper subforms of $f$, does $g$ represent $f$? We show that if $g$ is indefinite, then $g$ indeed represents $f$. However, when $f$ is positive definite and indecomposable, then there exists a $g$ which represents all proper subforms of $f$ but not $f$ itself. Along the way we give a new characterization of positive definite decomposable quadratic forms over $\mathfrak o$ and a number-field generalization of the finiteness theorem of representations of quadratic forms by quadratic forms over $\mathbb Z$ which asserts that given any infinite set $\mathscr S$ of classes of positive definite integral quadratic forms over $\mathfrak o$ of a fixed rank, there exists a finite subset $\mathscr S_0$ of $\mathscr S$ with the property that a positive definite quadratic form over $\mathfrak o$ represents all classes in $\mathscr S$ if and only if it represents all classes in $\mathscr S_0$. △ Less

Submitted 22 September, 2023; v1 submitted 21 January, 2022; originally announced January 2022.

Quantum computational study of chloride ion attack on chloromethane for chemical accuracy and quantum noise effects with UCCSD and k-UpCCGSD ansatzes

Authors: Hocheol Lim, Hyeon-Nae Jeon, June-Koo Rhee, Byungdu Oh, Kyoung Tai No

Abstract: Quantum computing is expected to play an important role in solving the problem of huge computational costs in various applications by utilizing the collective properties of quantum states, including superposition, interference, and entanglement, to perform computations. Quantum mechanical (QM) methods are candidates for various applications and can provide accurate absolute energy calculations in… ▽ More Quantum computing is expected to play an important role in solving the problem of huge computational costs in various applications by utilizing the collective properties of quantum states, including superposition, interference, and entanglement, to perform computations. Quantum mechanical (QM) methods are candidates for various applications and can provide accurate absolute energy calculations in structure-based methods. QM methods are powerful tools for describing reaction pathways and their potential energy surfaces (PESs). In this study, we applied quantum computing to describe the PES of the bimolecular nucleophilic substitution (SN2) reaction between chloromethane and chloride ions. We performed noiseless and noise simulations using quantum algorithms and compared the accuracy and noise effects of the ansatzes. In noiseless simulations, the results from UCCSD and k-UpCCGSD are similar to those of full configurational interaction (FCI) with the same active space, which indicates that quantum algorithms can describe the PES of the SN2 reaction. In noise simulations, UCCSD is more susceptible to quantum noise than k-UpCCGSD. Therefore, k-UpCCGSD can serve as an alternative to UCCSD to reduce quantum noisy effects in the noisy intermediate-scale quantum era, and k-UpCCGSD is sufficient to describe the PES of the SN2 reaction in this work. The results showed the applicability of quantum computing to the SN2 reaction pathway and provided valuable information for structure-based molecular simulations with quantum computing. △ Less

Submitted 31 December, 2021; originally announced December 2021.

Comments: 16 pages, 5 figures

Quantum readout error mitigation via deep learning

Authors: Jihye Kim, Byungdu Oh, Yonuk Chong, Euyheon Hwang, Daniel K. Park

Abstract: Quantum computing devices are inevitably subject to errors. To leverage quantum technologies for computational benefits in practical applications, quantum algorithms and protocols must be implemented reliably under noise and imperfections. Since noise and imperfections limit the size of quantum circuits that can be realized on a quantum device, developing quantum error mitigation techniques that d… ▽ More Quantum computing devices are inevitably subject to errors. To leverage quantum technologies for computational benefits in practical applications, quantum algorithms and protocols must be implemented reliably under noise and imperfections. Since noise and imperfections limit the size of quantum circuits that can be realized on a quantum device, developing quantum error mitigation techniques that do not require extra qubits and gates is of critical importance. In this work, we present a deep learning-based protocol for reducing readout errors on quantum hardware. Our technique is based on training an artificial neural network with the measurement results obtained from experiments with simple quantum circuits consisting of singe-qubit gates only. With the neural network and deep learning, non-linear noise can be corrected, which is not possible with the existing linear inversion methods. The advantage of our method against the existing methods is demonstrated through quantum readout error mitigation experiments performed on IBM five-qubit quantum devices. △ Less

Submitted 7 December, 2021; originally announced December 2021.

Comments: 10 pages, 10 figures

Journal ref: New Journal of Physics 24 073009 (2022)

The rank of new regular quadratic forms

Authors: Mingyu Kim, Byeong-Kweon Oh

Abstract: A (positive definite and integral) quadratic form $f$ is called regular if it represents all integers that are locally represented. It is known that there are only finitely many regular ternary quadratic forms up to isometry. However, there are infinitely many equivalence classes of regular quadratic forms of rank $n$ for any integer $n$ greater than or equal to $4$. A regular quadratic form $f$ i… ▽ More A (positive definite and integral) quadratic form $f$ is called regular if it represents all integers that are locally represented. It is known that there are only finitely many regular ternary quadratic forms up to isometry. However, there are infinitely many equivalence classes of regular quadratic forms of rank $n$ for any integer $n$ greater than or equal to $4$. A regular quadratic form $f$ is called new if there does not exist a proper subform $g$ of $f$ such that the set of integers that are represented by $g$ is equal to the set of integers that are represented by $f$. In this article, we prove that the rank of any new regular quadratic form is bounded by an absolute constant. △ Less

Submitted 19 November, 2021; originally announced November 2021.

MSC Class: 11E12; 11E20

How Metals Are Transported In And Out Of A Galactic Disk: Dependence On The Hydrodynamic Schemes In Numerical Simulations

Authors: Eun-jin Shin, Ji-hoon Kim, Boon Kiat Oh

Abstract: Metallicity is a fundamental probe for understanding the baryon physics in a galaxy. Since metals are intricately associated with radiative cooling, star formation, and feedback, reproducing the observed metal distribution through numerical experiments will provide a prominent way to examine our understandings of galactic baryon physics. In this study, we analyze the dependence of the galactic met… ▽ More Metallicity is a fundamental probe for understanding the baryon physics in a galaxy. Since metals are intricately associated with radiative cooling, star formation, and feedback, reproducing the observed metal distribution through numerical experiments will provide a prominent way to examine our understandings of galactic baryon physics. In this study, we analyze the dependence of the galactic metal distribution on the numerical schemes and quantify the differences in the metal mixing among modern galaxy simulation codes (the mesh-based code Enzo and the particle-based codes Gadget-2 and Gizmo-PSPH). In particular, we examine different stellar feedback strengths and an explicit metal diffusion scheme in particle-based codes, as a way to alleviate the well-known discrepancy in metal transport between mesh-based and particle-based simulations. We demonstrate that a sufficient number of gas particles are needed in the gas halo to properly investigate the metal distribution therein. Including an explicit metal diffusion scheme does not significantly affect the metal distribution in the galactic disk but does change the amount of low-metallicity gas in the hot-diffuse halo. We also find that the spatial distribution of metals depends strongly on how the stellar feedback is modeled. We demonstrate that the previously reported discrepancy in metals between mesh-based and particle-based simulations can be mitigated with our proposed prescription, enabling these simulations to be reliably utilized in the study of metals in galactic halos and the circumgalactic medium. △ Less

Submitted 8 June, 2021; originally announced June 2021.

Comments: 18 pages including 12 figures. Accepted for publication in ApJ. For associated movies and information, see http://ejshin0.github.io/page/category/visualization2.html and http://physics.snu.ac.kr/cosmo/index/research.html

Isolations of cubic lattices from their proper sublattices

Authors: Byeong-Kweon Oh

Abstract: A (positive definite and integral) quadratic form is called {\it an isolation} of a quadratic form $f$ if it represents all subforms of $f$ except for $f$ itself. The minimum rank of isolations of a quadratic form $f$ is denoted, if it exists, by $\text{Iso}(f)$. In this article, we show that $\text{Iso}(I_2)=5$ and $\text{Iso}(I_3)=6$, where $I_n=x_1^2+\dots+x_n^2$ is the sum of $n$ squares for a… ▽ More A (positive definite and integral) quadratic form is called {\it an isolation} of a quadratic form $f$ if it represents all subforms of $f$ except for $f$ itself. The minimum rank of isolations of a quadratic form $f$ is denoted, if it exists, by $\text{Iso}(f)$. In this article, we show that $\text{Iso}(I_2)=5$ and $\text{Iso}(I_3)=6$, where $I_n=x_1^2+\dots+x_n^2$ is the sum of $n$ squares for any positive integer $n$. After proving that there always exists an isolation of $I_n$ for any positive integer $n$, we provide some explicit lower and upper bounds for $\text{Iso}(I_n)$. In particular, we show that $\text{Iso}(I_n) \in Ω(n^{\frac32-ε})$ for any $ε>0$. △ Less

Submitted 9 April, 2021; originally announced April 2021.

MSC Class: 11E12; 11E25

Tight universal quadratic forms

Authors: Mingyu Kim, Byeong-Kweon Oh

Abstract: For a positive integer $n$, let $\mathcal T(n)$ be the set of all integers greater than or equal to $n$. An integral quadratic form $f$ is called tight $\mathcal T(n)$-universal if the set of nonzero integers that are represented by $f$ is exactly $\mathcal T(n)$. The smallest possible rank over all tight $\mathcal T(n)$-universal quadratic forms is defined by $t(n)$. In this article, we find all… ▽ More For a positive integer $n$, let $\mathcal T(n)$ be the set of all integers greater than or equal to $n$. An integral quadratic form $f$ is called tight $\mathcal T(n)$-universal if the set of nonzero integers that are represented by $f$ is exactly $\mathcal T(n)$. The smallest possible rank over all tight $\mathcal T(n)$-universal quadratic forms is defined by $t(n)$. In this article, we find all tight $\mathcal T(n)$-universal diagonal quadratic forms. We also prove that $t(n) \in Ω(\log_2(n)) \cap O(\sqrt{n})$. Explicit lower and upper bounds for $t(n)$ will be provided for some small integer $n$. △ Less

Submitted 6 April, 2021; originally announced April 2021.

Evolving beyond z=0: insights about the future of stars and the intergalactic medium

Authors: Boon Kiat Oh, John A. Peacock, Sadegh Khochfar, Britton D. Smith

Abstract: We present results from seven cosmological simulations that have been extended beyond the present era as far as redshift $z=-0.995$ or $t\approx96\,{\rm Gyr}$, using the Enzo simulation code. We adopt the calibrated star formation and feedback prescriptions from our previous work on reproducing the Milky Way with Enzo with modifications to the simulation code, chemistry and cooling library. We the… ▽ More We present results from seven cosmological simulations that have been extended beyond the present era as far as redshift $z=-0.995$ or $t\approx96\,{\rm Gyr}$, using the Enzo simulation code. We adopt the calibrated star formation and feedback prescriptions from our previous work on reproducing the Milky Way with Enzo with modifications to the simulation code, chemistry and cooling library. We then consider the future behaviour of the halo mass function (HMF), the equation of state (EOS) of the IGM, and the cosmic star formation history (SFH). Consistent with previous work, we find a freeze-out in the HMF at $z\approx-0.6$. The evolution of the EOS of the IGM presents an interesting case study of the cosmological coincidence problem, where there is a sharp decline in the IGM temperature immediately after $z=0$. For the SFH, the simulations produce a peak and a subsequent decline into the future. However, we do find a turnaround in the SFH after $z\approx-0.98$ in some simulations, probably due to the limitations of the criteria used for star formation. By integrating the SFH in time up to $z=-0.92$, the simulation with the best spatial resolution predicts an asymptotic total stellar mass that is very close to that obtained from extrapolating the fit of the observed SFR. Lastly, we investigate the future evolution of the partition of baryons within a Milky Way-sized galaxy, using both a zoom and a box simulation. Despite vastly different resolutions, these simulations predict individual haloes containing an equal fraction of baryons in stars and gas at the time of freeze-out ($t\approx30\,{\rm Gyr}$). △ Less

Submitted 3 March, 2021; originally announced March 2021.

Sharp isoperimetric inequalities for infinite plane graphs with bounded vertex and face degrees

Authors: Byung-Geun Oh

Abstract: We provide sharp bounds for the isoperimetric constants of infinite plane graphs (tessellations) with bounded vertex and face degrees. For example, if $G$ is a plane graph satisfying the inequalities $p_1 \leq \mbox{deg}\ v \leq p_2$ for $v \in V(G)$ and $q_1 \leq \mbox{deg}\ f \leq q_2$ for $f \in F(G)$, where $p_1, p_2, q_1$, and $q_2$ are natural numbers such that $1/p_i + 1/q_i \leq 1/2$,… ▽ More We provide sharp bounds for the isoperimetric constants of infinite plane graphs (tessellations) with bounded vertex and face degrees. For example, if $G$ is a plane graph satisfying the inequalities $p_1 \leq \mbox{deg}\ v \leq p_2$ for $v \in V(G)$ and $q_1 \leq \mbox{deg}\ f \leq q_2$ for $f \in F(G)$, where $p_1, p_2, q_1$, and $q_2$ are natural numbers such that $1/p_i + 1/q_i \leq 1/2$, $i=1,2$, then we show that \[ Φ(p_1, q_1) \leq \inf_S \frac{|\partial S|}{|V(S)|} \leq Φ(p_2, q_2), \] where the infimum is taken over all finite nonempty subgraphs $S \subset G$, $\partial S$ is the set of edges connecting $S$ to $G \setminus S$, and $Φ(p,q)$ is defined by \[ Φ(p, q) = (p-2) \sqrt{1 - \frac{4}{(p-2)(q-2)}}. \] For $p_1=3$ this gives an affirmative answer to a conjecture by Lawrencenko, Plummer, and Zha from 2002, and for general $p_i$ and $q_i$ our result fully resolves a question posed in the book by Lyons and Peres from 2016, where they extended the conjecture of Lawrencenko et al. to the above form. △ Less

Submitted 19 August, 2024; v1 submitted 9 September, 2020; originally announced September 2020.

Comments: 33 pages, 10 figures. Some sections are removed from the previous version

MSC Class: 05C10; 05B45; 05C63; 52C20

Minimal universality criterion sets on the representations of quadratic forms

Authors: Kyoungmin Kim, Jeongwon Lee, Byeong-Kweon Oh

Abstract: For a set $S$ of (positive definite and integral) quadratic forms with bounded rank, a quadratic form $f$ is called $S$-universal if it represents all quadratic forms in $S$. A subset $S_0$ of $S$ is called an $S$-universality criterion set if any $S_0$-universal quadratic form is $S$-universal. We say $S_0$ is minimal if there does not exist a proper subset of $S_0$ that is an $S$-universality cr… ▽ More For a set $S$ of (positive definite and integral) quadratic forms with bounded rank, a quadratic form $f$ is called $S$-universal if it represents all quadratic forms in $S$. A subset $S_0$ of $S$ is called an $S$-universality criterion set if any $S_0$-universal quadratic form is $S$-universal. We say $S_0$ is minimal if there does not exist a proper subset of $S_0$ that is an $S$-universality criterion set. In this article, we study various properties of minimal universality criterion sets. In particular, we show that for `most' binary quadratic forms $f$, minimal $S$-universality criterion sets are unique in the case when $S$ is the set of all subforms of the binary form $f$. △ Less

Submitted 8 September, 2020; originally announced September 2020.

Comments: 19 pages

MSC Class: 11E12; 11E20

Dark Matter Deficient Galaxies Produced Via High-velocity Galaxy Collisions In High-resolution Numerical Simulations

Authors: Eun-jin Shin, Minyong Jung, Goojin Kwon, Ji-hoon Kim, Joohyun Lee, Yongseok Jo, Boon Kiat Oh

Abstract: The recent discovery of diffuse dwarf galaxies that are deficient in dark matter appears to challenge the current paradigm of structure formation in our Universe. We describe the numerical experiments to determine if the so-called dark matter deficient galaxies (DMDGs) could be produced when two gas-rich, dwarf-sized galaxies collide with a high relative velocity of $\sim 300\,{\rm kms^{-1}}$. Usi… ▽ More The recent discovery of diffuse dwarf galaxies that are deficient in dark matter appears to challenge the current paradigm of structure formation in our Universe. We describe the numerical experiments to determine if the so-called dark matter deficient galaxies (DMDGs) could be produced when two gas-rich, dwarf-sized galaxies collide with a high relative velocity of $\sim 300\,{\rm kms^{-1}}$. Using idealized high-resolution simulations with both mesh-based and particle-based gravito-hydrodynamics codes, we find that DMDGs can form as high-velocity galaxy collisions separate dark matter from the warm disk gas which subsequently is compressed by shock and tidal interaction to form stars. Then using a large simulated universe IllustrisTNG, we discover a number of high-velocity galaxy collision events in which DMDGs are expected to form. However, we did not find evidence that these types of collisions actually produced DMDGs in the TNG100-1 run. We argue that the resolution of the numerical experiment is critical to realize the "collision-induced" DMDG formation scenario. Our results demonstrate one of many routes in which galaxies could form with unconventional dark matter fractions. △ Less

Submitted 26 July, 2020; v1 submitted 20 July, 2020; originally announced July 2020.

Comments: 18 pages including 13 figures. Accepted for publication in ApJ. For associated movies and information, see https://ejshin0.github.io/page/category/visualizations.html and https://physics.snu.ac.kr/cosmo/index/research.html#dmdg

The Cassels heights of cyclotomic integers

Authors: James McKee, Byeong-Kweon Oh, Chris Smyth

Abstract: We study the set $\mathscr C$ of mean square values of the moduli of the conjugates of cyclotomic integers $β$. For its $k$th derived set $\mathscr C^{(k)}$, we show that $\mathscr C^{(k)}=(k+1)\mathscr C\,\, (k\ge 0)$, so that also ${\mathscr C}^{(k)}+{\mathscr C}^{(\ell)}={\mathscr C}^{(k+\ell+1)}\,\,(k,\ell\ge 0)$. We also calculate the order type of $\mathscr C$, and show that it is the same a… ▽ More We study the set $\mathscr C$ of mean square values of the moduli of the conjugates of cyclotomic integers $β$. For its $k$th derived set $\mathscr C^{(k)}$, we show that $\mathscr C^{(k)}=(k+1)\mathscr C\,\, (k\ge 0)$, so that also ${\mathscr C}^{(k)}+{\mathscr C}^{(\ell)}={\mathscr C}^{(k+\ell+1)}\,\,(k,\ell\ge 0)$. We also calculate the order type of $\mathscr C$, and show that it is the same as that of the set of PV numbers. Furthermore, we describe precisely the restricted set $\mathscr C_p$ where the $β$ are confined to the ring $\mathbb Z[ω_p]$, where $p$ is an odd prime and $ω_p$ is a primitive $p$th root of unity. In order to do this, we prove that both of the quadratic polynomials $a^2+ab+b^2+c^2+a+b+c$ and $a^2+b^2+c^2+ab+bc+ca+a+b+c$ are universal. △ Less

Submitted 1 July, 2020; originally announced July 2020.

Comments: 13 pages

MSC Class: 11D85; 11R18

Prime-universal diagonal quadratic forms

Authors: Jangwon Ju, Daejun Kim, Kyoungmin Kim, Mingyu Kim, Byeong-Kweon Oh

Abstract: A (positive definite and integral) quadratic form is said to be $\textit{prime-universal}$ if it represents all primes. Recently, Doyle and Williams in [2] classified all prime-universal diagonal ternary quadratic forms, and all prime-universal diagonal quaternary quadratic forms under two conjectures proposed by themselves. In this article, we classify all prime-universal diagonal quadratic forms… ▽ More A (positive definite and integral) quadratic form is said to be $\textit{prime-universal}$ if it represents all primes. Recently, Doyle and Williams in [2] classified all prime-universal diagonal ternary quadratic forms, and all prime-universal diagonal quaternary quadratic forms under two conjectures proposed by themselves. In this article, we classify all prime-universal diagonal quadratic forms regardless of ranks. Furthermore, we prove, so called, $67$-Theorem for a diagonal quadratic form to be prime-universal. △ Less

Submitted 26 June, 2020; originally announced June 2020.

Comments: 14 pages

MSC Class: 11E12; 11E20

On the exceptional sets of integral quadratic forms

Authors: Wai Kiu Chan, Byeong-Kweon Oh

Abstract: A collection $\mathcal S$ of equivalence classes of positive definite integral quadratic forms in $n$ variables is called an $n$-exceptional set if there exists a positive definite integral quadratic form which represents all equivalence classes of positive definite integral quadratic forms in $n$ variables except those in $\mathcal S$. We show that, among other results, for any given positive int… ▽ More A collection $\mathcal S$ of equivalence classes of positive definite integral quadratic forms in $n$ variables is called an $n$-exceptional set if there exists a positive definite integral quadratic form which represents all equivalence classes of positive definite integral quadratic forms in $n$ variables except those in $\mathcal S$. We show that, among other results, for any given positive integers $m$ and $n$, there is always an $n$-exceptional set of size $m$ and there are only finitely many of them. △ Less

Submitted 25 March, 2020; originally announced March 2020.

MSC Class: 11E12; 11E25

Some criteria for circle packing types and combinatorial Gauss-Bonnet Theorem

Authors: Byung-Geun Oh

Abstract: We investigate criteria for circle packing(CP) types of disk triangulation graphs embedded into simply connected domains in $ \mathbb{C}$. In particular, by studying combinatorial curvature and the combinatorial Gauss-Bonnet theorem involving boundary turns, we show that a disk triangulation graph is CP parabolic if \[ \sum_{n=1}^\infty \frac{1}{\sum_{j=0}^{n-1} (k_j +6)} = \infty, \] where… ▽ More We investigate criteria for circle packing(CP) types of disk triangulation graphs embedded into simply connected domains in $ \mathbb{C}$. In particular, by studying combinatorial curvature and the combinatorial Gauss-Bonnet theorem involving boundary turns, we show that a disk triangulation graph is CP parabolic if \[ \sum_{n=1}^\infty \frac{1}{\sum_{j=0}^{n-1} (k_j +6)} = \infty, \] where $k_n$ is the degree excess sequence defined by \[ k_n = \sum_{v \in B_n} (\mbox{deg}\, v - 6) \] for combinatorial balls $B_n$ of radius $n$ and centered at a fixed vertex. It is also shown that the simple random walk on a disk triangulation graph is recurrent if \[ \sum_{n=1}^\infty \frac{1}{\sum_{j=0}^{n-1} (k_j +6)+\sum_{j=0}^{n} (k_j +6)} = \infty. \] These criteria are sharp, and generalize a conjecture by He and Schramm in their paper from 1995, which was later proved by Repp in 2001. We also give several criteria for CP hyperbolicity, one of which generalizes a theorem of He and Schramm, and present a necessary and sufficient condition for CP types of layered circle packings generalizing and confirming a criterion given by Siders in 1998. △ Less

Submitted 8 August, 2021; v1 submitted 16 March, 2020; originally announced March 2020.

Comments: 45 pages, 19 figures; to appear in TAMS

MSC Class: 52C15; 05B40; 05C10

Calibration of a star formation and feedback model for cosmological simulations with Enzo

Authors: Boon Kiat Oh, Britton D. Smith, John A. Peacock, Sadegh Khochfar

Abstract: We present results from seventy-one zoom simulations of a Milky Way-sized (MW) halo, exploring the parameter space for a widely-used star formation and feedback model in the {\tt Enzo} simulation code. We propose a novel way to match observations, using functional fits to the observed baryon makeup over a wide range of halo masses. The model MW galaxy is calibrated using three parameters: the star… ▽ More We present results from seventy-one zoom simulations of a Milky Way-sized (MW) halo, exploring the parameter space for a widely-used star formation and feedback model in the {\tt Enzo} simulation code. We propose a novel way to match observations, using functional fits to the observed baryon makeup over a wide range of halo masses. The model MW galaxy is calibrated using three parameters: the star formation efficiency $\left(f_*\right)$, the efficiency of thermal energy from stellar feedback $\left(ε\right)$ and the region into which feedback is injected $\left(r\ {\rm and}\ s\right)$. We find that changing the amount of feedback energy affects the baryon content most significantly. We then identify two sets of feedback parameter values that are both able to reproduce the baryonic properties for haloes between $10^{10}\,\mathrm{M_\odot}$ and $10^{12}\,\mathrm{M_\odot}$. We can potentially improve the agreement by incorporating more parameters or physics. If we choose to focus on one property at a time, we can obtain a more realistic halo baryon makeup. We show that the employed feedback prescription is insensitive to dark matter mass resolution between $10^5\,{\rm M_\odot}$ and $10^7\,{\rm M_\odot}$. Contrasting both star formation criteria and the corresponding combination of optimal feedback parameters, we also highlight that feedback is self-consistent: to match the same baryonic properties, with a relatively higher gas to stars conversion efficiency, the feedback strength required is lower, and vice versa. Lastly, we demonstrate that chaotic variance in the code can cause deviations of approximately 10\% and 25\% in the stellar and baryon mass in simulations evolved from identical initial conditions. △ Less

Submitted 7 February, 2020; originally announced February 2020.

Representations of finite number of quadratic forms with same rank

Authors: Daejun Kim, Byeong-Kweon Oh

Abstract: Let $m, n$ be positive integers with $m\le n$. Let $κ(m,n)$ be the largest integer $k$ such that for any (positive definite and integral) quadratic forms $f_1,\ldots,f_k$ of rank $m$, there exists a quadratic form of rank $n$ that represents $f_i$ for any $i$ with $1\le i \le k$. In this article, we determine the number $κ(m,n)$ for any integer $m$ with $1\le m\le 8$, except for the cases when… ▽ More Let $m, n$ be positive integers with $m\le n$. Let $κ(m,n)$ be the largest integer $k$ such that for any (positive definite and integral) quadratic forms $f_1,\ldots,f_k$ of rank $m$, there exists a quadratic form of rank $n$ that represents $f_i$ for any $i$ with $1\le i \le k$. In this article, we determine the number $κ(m,n)$ for any integer $m$ with $1\le m\le 8$, except for the cases when $(m,n)=(3,5)$ and $(4,6)$. In the exceptional cases, it will be proved that $1\le κ(3,5), \ κ(4,6)\le 2$. We also discuss some related topics. △ Less

Submitted 30 December, 2019; originally announced December 2019.

Comments: 13 pages

A sum of three nonunit squares of integers

Authors: Daejun Kim, Jeongwon Lee, Byeong-Kweon Oh

Abstract: We say a positive integer is a sum of three nonunit squares if it is a sum of three squares of integers other than one. In this article, we find all integers which are sums of three nonunit squares assuming that the Generalized Riemann Hypothesis(GRH) holds. As applications, we find all integers, under the GRH only when $k=3$, which are sums of $k$ nonzero triangular numbers, sums of $k$ nonzero g… ▽ More We say a positive integer is a sum of three nonunit squares if it is a sum of three squares of integers other than one. In this article, we find all integers which are sums of three nonunit squares assuming that the Generalized Riemann Hypothesis(GRH) holds. As applications, we find all integers, under the GRH only when $k=3$, which are sums of $k$ nonzero triangular numbers, sums of $k$ nonzero generalized pentagonal numbers, and sums of $k$ nonzero generalized octagonal numbers, respectively for any integer $k\ge 3$. △ Less

Submitted 11 September, 2019; originally announced September 2019.

Comments: 19 pages

MSC Class: 11E12; 11E20; 11E25

Observation of a partially rotating superfluid of exciton-polariton

Authors: Daegwang Choi, Min Park, Byoung Yong Oh, Min-Sik Kwon, Suk In Park, Sooseok Kang, Jin Dong Song, Yong-Hoon Cho, Hyoungsoon Choi

Abstract: Rotation of a container holding a viscous fluid forms a vortex which grows with increasing angular velocity. A superfluid, however, is intrinsically different from these normal fluids because its rotation is quantized. Even if a container of superfluid is rotating, the fluid itself remains still until a critical velocity is reached. Beyond the critical velocity, all the particles conspire to sudde… ▽ More Rotation of a container holding a viscous fluid forms a vortex which grows with increasing angular velocity. A superfluid, however, is intrinsically different from these normal fluids because its rotation is quantized. Even if a container of superfluid is rotating, the fluid itself remains still until a critical velocity is reached. Beyond the critical velocity, all the particles conspire to suddenly pick up an angular momentum of $\hbar$ each and forms a quantized vortex. As a result, a superfluid is known to increase its rotation by a total angular momentum of $N\hbar$. In this letter, we show that exciton-polariton superfluid can split into an irrotational part and a rotational part. The relative ratio between the two states can be controlled by either pump beam's power or spot size. Consequently, angular momentum of exciton-polariton superfluid can be tuned from zero to $N\hbar$ continuously. This striking observation sets the stage for studying non-equilibrium properties of a superfluid with exciton-polaritons. △ Less

Submitted 20 October, 2020; v1 submitted 5 September, 2019; originally announced September 2019.

Comments: A follow-up study on the two-states in the exciton-polariton superfluid reported in this article revealed that the two states occur sequentially and not concurrently. The apparent coexistence of the two states, thus, is not due to partial rotation of a superfluid. This new finding alters the conclusion of this article significantly so as to necessitate a withdrawal instead of an update

Quadratic forms with a strong regularity property on the representations of squares

Authors: Kyoungmin Kim, Byeong-Kweon Oh

Abstract: A (positive definite and non-classic integral) quadratic form is called strongly $s$-regular if it satisfies a strong regularity property on the number of representations of squares of integers. In this article, we prove that for any integer $k \ge 2$, there are only finitely many isometry classes of strongly $s$-regular quadratic forms with rank $k$ if the minimum of the nonzero squares that are… ▽ More A (positive definite and non-classic integral) quadratic form is called strongly $s$-regular if it satisfies a strong regularity property on the number of representations of squares of integers. In this article, we prove that for any integer $k \ge 2$, there are only finitely many isometry classes of strongly $s$-regular quadratic forms with rank $k$ if the minimum of the nonzero squares that are represented by them is fixed. △ Less

Submitted 4 September, 2019; originally announced September 2019.

Comments: 14 pages

MSC Class: 11E12; 11E20

DAPAS : Denoising Autoencoder to Prevent Adversarial attack in Semantic Segmentation

Authors: Seungju Cho, Tae Joon Jun, Byungsoo Oh, Daeyoung Kim

Abstract: Nowadays, Deep learning techniques show dramatic performance on computer vision area, and they even outperform human. But it is also vulnerable to some small perturbation called an adversarial attack. This is a problem combined with the safety of artificial intelligence, which has recently been studied a lot. These attacks have shown that they can fool models of image classification, semantic segm… ▽ More Nowadays, Deep learning techniques show dramatic performance on computer vision area, and they even outperform human. But it is also vulnerable to some small perturbation called an adversarial attack. This is a problem combined with the safety of artificial intelligence, which has recently been studied a lot. These attacks have shown that they can fool models of image classification, semantic segmentation, and object detection. We point out this attack can be protected by denoise autoencoder, which is used for denoising the perturbation and restoring the original images. We experiment with various noise distributions and verify the effect of denoise autoencoder against adversarial attack in semantic segmentation. △ Less

Submitted 7 April, 2020; v1 submitted 14 August, 2019; originally announced August 2019.

Comments: Accepted to be published in: 2020 International Joint Conference on Neural Networks (IJCNN), Glasgow, July 19--24, 2020

Regular ternary triangular forms

Authors: Mingyu Kim, Byeong-Kweon Oh

Abstract: An integer of the form $T_x=\frac{x(x+1)}2$ for some positive integer $x$ is called a triangular number. A ternary triangular form $aT_{x}+bT_{y}+cT_{z}$ for positive integers $a,b$ and $c$ is called regular if it represents every positive integer that is locally represented. In this article, we prove that there are exactly 49 primitive regular ternary triangular forms. An integer of the form $T_x=\frac{x(x+1)}2$ for some positive integer $x$ is called a triangular number. A ternary triangular form $aT_{x}+bT_{y}+cT_{z}$ for positive integers $a,b$ and $c$ is called regular if it represents every positive integer that is locally represented. In this article, we prove that there are exactly 49 primitive regular ternary triangular forms. △ Less

Submitted 8 March, 2019; originally announced March 2019.

Comments: 28 pages

MSC Class: 11E12; 11E20

Gradient-Free Learning Based on the Kernel and the Range Space

Authors: Kar-Ann Toh, Zhiping Lin, Zhengguo Li, Beomseok Oh, Lei Sun

Abstract: In this article, we show that solving the system of linear equations by manipulating the kernel and the range space is equivalent to solving the problem of least squares error approximation. This establishes the ground for a gradient-free learning search when the system can be expressed in the form of a linear matrix equation. When the nonlinear activation function is invertible, the learning prob… ▽ More In this article, we show that solving the system of linear equations by manipulating the kernel and the range space is equivalent to solving the problem of least squares error approximation. This establishes the ground for a gradient-free learning search when the system can be expressed in the form of a linear matrix equation. When the nonlinear activation function is invertible, the learning problem of a fully-connected multilayer feedforward neural network can be easily adapted for this novel learning framework. By a series of kernel and range space manipulations, it turns out that such a network learning boils down to solving a set of cross-coupling equations. By having the weights randomly initialized, the equations can be decoupled and the network solution shows relatively good learning capability for real world data sets of small to moderate dimensions. Based on the structural information of the matrix equation, the network representation is found to be dependent on the number of data samples and the output dimension. △ Less

Submitted 26 October, 2018; originally announced October 2018.

Comments: The idea of kernel and range projection was first introduced in the IEEE/ACIS ICIS conference which was held in Singapore in June 2018. This article presents a full development of the method supported by extensive numerical results