Showing 1–50 of 404 results for author: Yoshida, N

Fearless Asynchronous Communications with Timed Multiparty Session Protocols

Authors: Ping Hou, Nicolas Lagaillardie, Nobuko Yoshida

Abstract: Session types using affinity and exception handling mechanisms have been developed to ensure the communication safety of protocols implemented in concurrent and distributed programming languages. Nevertheless, current affine session types are inadequate for specifying real-world asynchronous protocols, as they are usually imposed by time constraints which enable timeout exceptions to prevent indef… ▽ More Session types using affinity and exception handling mechanisms have been developed to ensure the communication safety of protocols implemented in concurrent and distributed programming languages. Nevertheless, current affine session types are inadequate for specifying real-world asynchronous protocols, as they are usually imposed by time constraints which enable timeout exceptions to prevent indefinite blocking while awaiting valid messages. This paper proposes the first formal integration of affinity, time constraints, timeouts, and time-failure handling based on multiparty session types for supporting reliability in asynchronous distributed systems. With this theory, we statically guarantee that asynchronous timed communication is deadlock-free, communication safe, while being fearless -- never hindered by timeout errors or abrupt terminations. To implement our theory, we introduce a Rust toolchain designed to facilitate the implementation of safe affine timed protocols. Our toolchain leverages generic types and the time library to handle timed communications, integrated with optional types for affinity. We evaluate our approach by extending diverse examples from the literature to incorporate time and timeouts, demonstrating that our solution incurs negligible overhead compared with an untimed implementation. We also showcase the correctness by construction of our approach by implementing various real-world use cases, including a remote data protocol from the Internet of Remote Things domain, as well as protocols from real-time systems like Android motion sensors and smartwatches. △ Less

Submitted 12 July, 2024; v1 submitted 27 June, 2024; originally announced June 2024.

Comments: ECOOP 2024

JWST, ALMA, and Keck Spectroscopic Constraints on the UV Luminosity Functions at z~7-14: Clumpiness and Compactness of the Brightest Galaxies in the Early Universe

Authors: Yuichi Harikane, Akio K. Inoue, Richard S. Ellis, Masami Ouchi, Yurina Nakazato, Naoki Yoshida, Yoshiaki Ono, Fengwu Sun, Riku A. Sato, Seiji Fujimoto, Nobunari Kashikawa, Derek J. McLeod, Pablo G. Perez-Gonzalez, Marcin Sawicki, Yuma Sugahara, Yi Xu, Satoshi Yamanaka, Adam C. Carnall, Fergus Cullen, James S. Dunlop, Eiichi Egami, Norman Grogin, Yuki Isobe, Anton M. Koekemoer, Nicolas Laporte , et al. (10 additional authors not shown)

Abstract: We present the number densities and physical properties of the bright galaxies spectroscopically confirmed at $z\sim7-14$. Our sample is composed of 53 galaxies at $z_\mathrm{spec}\sim7-14$, including recently-confirmed galaxies at $z_\mathrm{spec}=12.34-14.32$ with JWST, as well as new confirmations at $z_\mathrm{spec}=6.583-7.643$ with $-24< M_\mathrm{UV}< -21$ mag using ALMA and Keck. Our JWST/… ▽ More We present the number densities and physical properties of the bright galaxies spectroscopically confirmed at $z\sim7-14$. Our sample is composed of 53 galaxies at $z_\mathrm{spec}\sim7-14$, including recently-confirmed galaxies at $z_\mathrm{spec}=12.34-14.32$ with JWST, as well as new confirmations at $z_\mathrm{spec}=6.583-7.643$ with $-24< M_\mathrm{UV}< -21$ mag using ALMA and Keck. Our JWST/NIRSpec observations have also revealed that very bright galaxy candidates at $z\sim10-13$ identified from ground-based telescope images before JWST are passive galaxies at $z\sim3-4$, emphasizing the necessity of strict screening and spectroscopy in the selection of the brightest galaxies at $z>10$. The UV luminosity functions derived from these spectroscopic results are consistent with a double power-law function, showing tensions with theoretical models at the bright end. To understand the origin of the overabundance of bright galaxies, we investigate their morphologies using JWST/NIRCam high-resolution images obtained in various surveys including PRIMER and COSMOS-Web. We find that $\sim70\%$ of the bright galaxies at $z\sim7$ exhibit clumpy morphologies with multiple sub-components, suggesting merger-induced starburst activity, which is consistent with SED fitting results showing bursty star formation histories. At $z\gtrsim10$, bright galaxies are classified into two types of galaxies; extended ones with weak high-ionization emission lines, and compact ones with strong high-ionization lines including NIV]$λ$1486, indicating that at least two different processes (e.g., merger-induced starburst and compact star formation/AGN) are shaping the physical properties of the brightest galaxies at $z\gtrsim10$ and are responsible for their overabundance. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 32 pages, 22 figures, submitted to ApJ

Effect of Random Learning Rate: Theoretical Analysis of SGD Dynamics in Non-Convex Optimization via Stationary Distribution

Authors: Naoki Yoshida, Shogo Nakakita, Masaaki Imaizumi

Abstract: We consider a variant of the stochastic gradient descent (SGD) with a random learning rate and reveal its convergence properties. SGD is a widely used stochastic optimization algorithm in machine learning, especially deep learning. Numerous studies reveal the convergence properties of SGD and its simplified variants. Among these, the analysis of convergence using a stationary distribution of updat… ▽ More We consider a variant of the stochastic gradient descent (SGD) with a random learning rate and reveal its convergence properties. SGD is a widely used stochastic optimization algorithm in machine learning, especially deep learning. Numerous studies reveal the convergence properties of SGD and its simplified variants. Among these, the analysis of convergence using a stationary distribution of updated parameters provides generalizable results. However, to obtain a stationary distribution, the update direction of the parameters must not degenerate, which limits the applicable variants of SGD. In this study, we consider a novel SGD variant, Poisson SGD, which has degenerated parameter update directions and instead utilizes a random learning rate. Consequently, we demonstrate that a distribution of a parameter updated by Poisson SGD converges to a stationary distribution under weak assumptions on a loss function. Based on this, we further show that Poisson SGD finds global minima in non-convex optimization problems and also evaluate the generalization error using this method. As a proof technique, we approximate the distribution by Poisson SGD with that of the bouncy particle sampler (BPS) and derive its stationary distribution, using the theoretical advance of the piece-wise deterministic Markov process (PDMP). △ Less

Submitted 23 June, 2024; originally announced June 2024.

Comments: 28 pages

The Supersonic Project: Early Star Formation with the Streaming Velocity

Authors: William Lake, Claire E. Williams, Smadar Naoz, Federico Marinacci, Blakesley Burkhart, Mark Vogelsberger, Naoki Yoshida, Gen Chiaki, Avi Chen, Yeou S. Chiou

Abstract: At high redshifts ($z\gtrsim12$), the relative velocity between baryons and dark matter (the so-called streaming velocity) significantly affects star formation in low-mass objects. Streaming substantially reduces the abundance of low-mass gas objects while simultaneously allowing for the formation of supersonically-induced gas objects (SIGOs) and their associated star clusters outside of dark matt… ▽ More At high redshifts ($z\gtrsim12$), the relative velocity between baryons and dark matter (the so-called streaming velocity) significantly affects star formation in low-mass objects. Streaming substantially reduces the abundance of low-mass gas objects while simultaneously allowing for the formation of supersonically-induced gas objects (SIGOs) and their associated star clusters outside of dark matter halos. Here, we present a study of the population-level effects of streaming on star formation within both halos and SIGOs in a set of simulations with and without streaming. Notably, we find that streaming actually enhances star formation within individual halos of all masses at redshifts between $z=12$ and $z=20$. This is demonstrated both as an increased star formation rate per object as well as an enhancement of the Kennicutt-Schmidt relation for objects with streaming. We find that our simulations are consistent with some observations at high redshift, but on a population level, they continue to under-predict star formation relative to the majority of observations. However, simulations of overdense regions (both with and without streaming) agree with observations, suggesting a strategy for extracting information about the overdensity and streaming velocity in a given survey volume in future observations. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Comments: 16 pages, 8 figures

Separation and Encodability in Mixed Choice Multiparty Sessions (Technical Report)

Authors: Kirstin Peters, Nobuko Yoshida

Abstract: Multiparty session types (MP) are a type discipline for enforcing the structured, deadlock-free communication of concurrent and message-passing programs. Traditional MP have a limited form of choice in which alternative communication possibilities are offered by a single participant and selected by another. Mixed choice multiparty session types (MCMP) extend the choice construct to include both se… ▽ More Multiparty session types (MP) are a type discipline for enforcing the structured, deadlock-free communication of concurrent and message-passing programs. Traditional MP have a limited form of choice in which alternative communication possibilities are offered by a single participant and selected by another. Mixed choice multiparty session types (MCMP) extend the choice construct to include both selections and offers in the same choice. This paper first proposes a general typing system for a mixed choice synchronous multiparty session calculus, and prove type soundness, communication safety, and deadlock-freedom. Next we compare expressiveness of nine subcalcli of MCMP-calculus by examining their encodability (there exists a good encoding from one to another) and separation (there exists no good encoding from one calculus to another). We prove 8 new encodablity results and 20 new separation results. In summary, MCMP is strictly more expressive than classical multiparty sessions (MP) and mixed choice in mixed sessions. This contrasts earlier results where mixed sessions do not add any expressiveness to non-mixed fundamental sessions, shedding a light on expressiveness of multiparty mixed choice. △ Less

Submitted 7 June, 2024; v1 submitted 13 May, 2024; originally announced May 2024.

Comments: Technical report of the paper Separation and Encodability in Mixed Choice Multiparty Sessions by Kirstin Peters and Nobuko Yoshida at LICS'24

Neural network prediction of model parameters for strong lensing samples from Hyper Suprime-Cam Survey

Authors: Priyanka Gawade, Anupreeta More, Surhud More, Akisato Kimura, Alessandro Sonnenfeld, Masamune Oguri, Naoki Yoshida

Abstract: Galaxies that cause the strong gravitational lensing of background galaxies provide us crucial information about the distribution of matter around them. Traditional modelling methods that analyse such strong lenses are both time and resource consuming, require sophisticated lensing codes and modelling expertise. To study the large lens population expected from imaging surveys such as LSST, we need… ▽ More Galaxies that cause the strong gravitational lensing of background galaxies provide us crucial information about the distribution of matter around them. Traditional modelling methods that analyse such strong lenses are both time and resource consuming, require sophisticated lensing codes and modelling expertise. To study the large lens population expected from imaging surveys such as LSST, we need fast and automated analysis methods. In this work, we build and train a simple convolutional neural network with an aim to rapidly predict model parameters of gravitational lenses. We focus on the most important lens mass model parameters, namely, the Einstein radius, the axis ratio and the position angle of the major axis of the mass distribution. The network is trained on a variety of simulated data with an increasing degree of realism and shows satisfactory performance on simulated test data. The trained network is then applied to the real sample of galaxy-scale candidate lenses from the Subaru HSC, a precursor survey to LSST. Unlike the simulated lenses, we do not have the ground truth for the real lenses. Therefore, we have compared our predictions with those from YattaLens, a lens modelling pipeline. Additionally, we also compare the parameter predictions for 10 HSC lenses that were also studied by other conventional modelling methods. These comparisons show a fair quantitative agreement on the Einstein radius, although the axis ratio and the position angle from the network as well as the individual modelling methods, seem to have systematic uncertainties beyond the quoted errors. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: 10 pages, 9 figures

Dark Matter Mounds: towards a realistic description of dark matter overdensities around black holes

Authors: Gianfranco Bertone, A. Renske A. C. Wierda, Daniele Gaggero, Bradley J. Kavanagh, Marta Volonteri, Naoki Yoshida

Abstract: Dark matter overdensities around black holes can be searched for by looking at the characteristic imprint they leave on the gravitational waveform of binary black hole mergers. Current theoretical predictions of the density profile of dark matter overdensities are based on highly idealised formation scenarios, in which black holes are assumed to grow adiabatically from an infinitesimal seed to the… ▽ More Dark matter overdensities around black holes can be searched for by looking at the characteristic imprint they leave on the gravitational waveform of binary black hole mergers. Current theoretical predictions of the density profile of dark matter overdensities are based on highly idealised formation scenarios, in which black holes are assumed to grow adiabatically from an infinitesimal seed to their final mass, compressing dark matter cusps at the center of galactic halos into very dense `spikes'. These scenarios were suitable for dark matter indirect detection studies, since annihilating dark matter cannot reach very high densities, but they fail to capture the dark matter distribution in the innermost regions where the gravitational wave signal is produced. We present here a more realistic formation scenario where black holes form from the collapse of supermassive stars, and follow the evolution of the dark matter density as the supermassive star grows and collapses to a black hole. We show that in this case dark matter forms shallower `mounds', instead of `spikes', on scales comparable with the size of the supermassive stars originating them. We discuss the implications for the detectability of these systems. △ Less

Submitted 12 April, 2024; originally announced April 2024.

Comments: 8 pages, 4 figures

Three Subtyping Algorithms for Binary Session Types and their Complexity Analyses

Authors: Thien Udomsrirungruang, Nobuko Yoshida

Abstract: Session types are a type discipline for describing and specifying communication behaviours of concurrent processes. Session subtyping, firstly introduced by Gay and Hole, is widely used for enlarging typability of session programs. This paper gives the complexity analysis of three algorithms for subtyping of synchronous binary session types. First, we analyse the complexity of the algorithm from t… ▽ More Session types are a type discipline for describing and specifying communication behaviours of concurrent processes. Session subtyping, firstly introduced by Gay and Hole, is widely used for enlarging typability of session programs. This paper gives the complexity analysis of three algorithms for subtyping of synchronous binary session types. First, we analyse the complexity of the algorithm from the original paper, which is based on an inductive tree search. We then introduce its optimised version, which improves the complexity, but is still exponential against the size of the two types. Finally, we propose a new quadratic algorithm based on a graph search using the concept of XYZW-simulation, recently introduced by Silva et al. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Comments: In Proceedings PLACES 2024, arXiv:2404.03712. A full version of this paper which contains complete definitions and examples is available at arXiv:2402.06988

Journal ref: EPTCS 401, 2024, pp. 49-60

Redshift-dependent galaxy formation efficiency at z=5-13 in the FirstLight simulations

Authors: Daniel Ceverino, Yurina Nakazato, Naoki Yoshida, Ralf Klessen, Simon Glover

Abstract: Current models of the formation of first galaxies predict low masses and faint objects at extremely high redshifts, z=9-15. However, the first observations of this epoch indicate a higher-than-expected number of bright (sometimes massive) galaxies. Numerical simulations can help to elucidate the mild evolution of the bright end of the UV luminosity function and they can provide the link between th… ▽ More Current models of the formation of first galaxies predict low masses and faint objects at extremely high redshifts, z=9-15. However, the first observations of this epoch indicate a higher-than-expected number of bright (sometimes massive) galaxies. Numerical simulations can help to elucidate the mild evolution of the bright end of the UV luminosity function and they can provide the link between the evolution of bright galaxies and variations of the galaxy formation efficiency across different redshifts. We use the FirstLight database of 377 zoom-in cosmological simulations of a mass-complete sample of galaxies. Mock luminosities are estimated by a dust model constrained by current observations of an evolution of the beta-MUV relation at high-z. FirstLight contains a high number of bright galaxies, MUV<-20, consistent with current data at z=6-13. The evolution of the UV cosmic density is driven by the evolution of the galaxy efficiency and the relation between MUV and halo mass. The efficiency of galaxy formation increases significantly with redshift at a fixed halo mass because galactic halos at extremely high redshifts convert gas into stars at a higher rate than at lower redshifts. The high gas densities in galaxies at z>9 enable these high efficiencies. Our simulations predict higher number densities of massive galaxies, Ms=10^9 Msun, than other models with constant efficiency. Cosmological simulations of galaxy formation with self-consistent models of star formation and feedback can reproduce the different regimes of galaxy formation across cosmic history. △ Less

Submitted 3 April, 2024; originally announced April 2024.

Comments: 10 pages, 10 figures, submitted to A&A

RIOJA. Complex Dusty Starbursts in a Major Merger B14-65666 at z=7.15

Authors: Yuma Sugahara, Javier Álvarez-Márquez, Takuya Hashimoto, Luis Colina, Akio K. Inoue, Luca Costantin, Yoshinobu Fudamoto, Ken Mawatari, Yi W. Ren, Santiago Arribas, Tom J. L. C. Bakx, Carmen Blanco-Prieto, Daniel Ceverino, Alejandro Crespo Gómez, Masato Hagimoto, Takeshi Hashigaya, Rui Marques-Chaves, Hiroshi Matsuo, Yurina Nakazato, Miguel Pereira-Santaella, Yoichi Tamura, Mitsutaka Usui, Naoki Yoshida

Abstract: We present JWST NIRCam imaging of B14-65666 ("Big Three Dragons"), a bright Lyman-break galaxy system ($M_\text{UV}=-22.5$ mag) at $z=7.15$. The high angular resolution of NIRCam reveals the complex morphology of two galaxy components: galaxy E has a compact core (E-core), surrounded by diffuse, extended, rest-frame optical emission, which is likely to be tidal tails; and galaxy W has a clumpy and… ▽ More We present JWST NIRCam imaging of B14-65666 ("Big Three Dragons"), a bright Lyman-break galaxy system ($M_\text{UV}=-22.5$ mag) at $z=7.15$. The high angular resolution of NIRCam reveals the complex morphology of two galaxy components: galaxy E has a compact core (E-core), surrounded by diffuse, extended, rest-frame optical emission, which is likely to be tidal tails; and galaxy W has a clumpy and elongated morphology with a blue UV slope ($β_\text{UV}=-2.2\pm0.1$). The flux excess, F356W$-$F444W, peaks at the E-core ($1.05^{+0.08}_{-0.09}$ mag), tracing the presence of strong [OIII] 4960,5008 Å emission. ALMA archival data show that the bluer galaxy W is brighter in dust continua than the redder galaxy E, while the tails are bright in [OIII] 88 $\mathrm{μm}$. The UV/optical and sub-mm SED fitting confirms that B14-65666 is a major merger in a starburst phase as derived from the stellar mass ratio (3:1 to 2:1) and the star-formation rate, $\simeq1$ dex higher than the star-formation main sequence at the same redshift. The galaxy E is a dusty ($A_\text{V}=1.2\pm0.1$ mag) starburst with a possible high dust temperature ($\ge63$-$68$ K). The galaxy W would have a low dust temperature ($\le27$-$33$ K) or patchy stellar-and-dust geometry, as suggested from the infrared excess (IRX) and $β_\text{UV}$ diagram. The high optical-to-FIR [OIII] line ratio of the E-core shows its lower gas-phase metallicity ($\simeq0.2$ Z$_{\odot}$) than the galaxy W. These results agree with a scenario where major mergers disturb morphology and induce nuclear dusty starbursts triggered by less-enriched inflows. B14-65666 shows a picture of complex stellar buildup processes during major mergers in the epoch of reionization. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Comments: 18 pages, 6 figures, 4 tables. Submitted to ApJ

Log-rank test with coarsened exact matching

Authors: Tomoya Baba, Nakahiro Yoshida

Abstract: It is of special importance in the clinical trial to compare survival times between the treatment group and the control group. Propensity score methods with a logistic regression model are often used to reduce the effects of confounders. However, the modeling of complex structures between the covariates, the treatment assignment and the survival time is difficult. In this paper, we consider coarse… ▽ More It is of special importance in the clinical trial to compare survival times between the treatment group and the control group. Propensity score methods with a logistic regression model are often used to reduce the effects of confounders. However, the modeling of complex structures between the covariates, the treatment assignment and the survival time is difficult. In this paper, we consider coarsened exact matching (CEM), which does not need any parametric models, and we propose the weighted log-rank statistic based on CEM. We derive asymptotic properties of the weighted log-rank statistic, such as the weak convergence to a Gaussian process in Skorokhod space, in particular the asymptotic normality, under the null hypothesis and the consistency of the log-rank test. Simulation studies show that the log-rank statistic based on CEM is more robust than the log-rank statistic based on the propensity score. △ Less

Submitted 26 March, 2024; v1 submitted 24 March, 2024; originally announced March 2024.

The transition operator of a random walk perturbated by sparse potentials

Authors: Takuya Mine, Nobuo Yoshida

Abstract: We consider an operator $P_V=(1+V)P$ on $\ell^2(Z^d)$, where $P$ is the transition operator of a symmetric irreducible random walk, and $V$ is a ``sparse'' potential. We first characterize the essential spectra of this operator. Secondly, we prove that all the eigenfunctions which correspond to discrete spectra decay exponentially fast. Thirdly, we give a sufficient condition for this operator to… ▽ More We consider an operator $P_V=(1+V)P$ on $\ell^2(Z^d)$, where $P$ is the transition operator of a symmetric irreducible random walk, and $V$ is a ``sparse'' potential. We first characterize the essential spectra of this operator. Secondly, we prove that all the eigenfunctions which correspond to discrete spectra decay exponentially fast. Thirdly, we give a sufficient condition for this operator to have an absolute spectral gap at the right edge of the spectra. Finally, as an application of the absolute spectral gap and the exponential decay of the eigenfunctions, we prove a limit theorem for the random walk under the Gibbs measure associated to the potential $V$. △ Less

Submitted 12 March, 2024; originally announced March 2024.

How to Save My Gas Fees: Understanding and Detecting Real-world Gas Issues in Solidity Programs

Authors: Mengting He, Shihao Xia, Boqin Qin, Nobuko Yoshida, Tingting Yu, Linhai Song, Yiying Zhang

Abstract: The execution of smart contracts on Ethereum, a public blockchain system, incurs a fee called gas fee for its computation and data-store consumption. When programmers develop smart contracts (e.g., in the Solidity programming language), they could unknowingly write code snippets that unnecessarily cause more gas fees. These issues, or what we call gas wastes, could lead to significant monetary was… ▽ More The execution of smart contracts on Ethereum, a public blockchain system, incurs a fee called gas fee for its computation and data-store consumption. When programmers develop smart contracts (e.g., in the Solidity programming language), they could unknowingly write code snippets that unnecessarily cause more gas fees. These issues, or what we call gas wastes, could lead to significant monetary waste for users. Yet, there have been no systematic examination of them or effective tools for detecting them. This paper takes the initiative in helping Ethereum users reduce their gas fees in two important steps: we conduct the first empirical study on gas wastes in popular smart contracts written in Solidity by understanding their root causes and fixing strategies; we then develop a static tool, PeCatch, to effectively detect gas wastes with simple fixes in Solidity programs based on our study findings. Overall, we make seven insights and four suggestions from our gas-waste study, which could foster future tool development, language improvement, and programmer awareness, and develop eight gas-waste checkers, which pinpoint 383 previously unknown gas wastes from famous Solidity libraries. △ Less

Submitted 5 March, 2024; originally announced March 2024.

RISMiCal: A software package to perform fast RISM/3D-RISM calculations

Authors: Yutaka Maruyama, Nori Yoshida

Abstract: Solvent plays an essential role in a variety of chemical, physical, and biological processes that occur in the solution phase. The reference interaction site model (RISM) and its three-dimensional extension (3D-RISM) serve as powerful computational tools for modeling solvation effects in chemical reactions, biological functions, and structure formations. We present the RISM integrated calculator (… ▽ More Solvent plays an essential role in a variety of chemical, physical, and biological processes that occur in the solution phase. The reference interaction site model (RISM) and its three-dimensional extension (3D-RISM) serve as powerful computational tools for modeling solvation effects in chemical reactions, biological functions, and structure formations. We present the RISM integrated calculator (RISMiCal) program package, which is based on RISM and 3D-RISM theories with fast GPU code. RISMiCal has been developed as an integrated RISM/3D-RISM program that has interfaces with external programs such as Gaussian16, GAMESS, and Tinker. Fast 3D-RISM programs for single- and multi-GPU codes written in CUDA would enhance the availability of these hybrid methods because they require the performance of many computationally expensive 3D-RISM calculations. We expect that our package can be widely applied for chemical and biological processes in solvent. The RISMiCal package is available at https://rismical-dev.github.io. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Asymptotic expansion of the drift estimator for the fractional Ornstein-Uhlenbeck process

Authors: Ciprian A. Tudor, Nakahiro Yoshida

Abstract: We present an asymptotic expansion formula of an estimator for the drift coefficient of the fractional Ornstein-Uhlenbeck process. As the machinery, we apply the general expansion scheme for Wiener functionals recently developed by the authors [26]. The central limit theorem in the principal part of the expansion has the classical scaling T^{1/2}. However, the asymptotic expansion formula is a com… ▽ More We present an asymptotic expansion formula of an estimator for the drift coefficient of the fractional Ornstein-Uhlenbeck process. As the machinery, we apply the general expansion scheme for Wiener functionals recently developed by the authors [26]. The central limit theorem in the principal part of the expansion has the classical scaling T^{1/2}. However, the asymptotic expansion formula is a complex in that the order of the correction term becomes the classical T^{-1/2} for H in (1/2,5/8), but T^{4H-3} for H in [5/8, 3/4). △ Less

Submitted 3 April, 2024; v1 submitted 1 March, 2024; originally announced March 2024.

MSC Class: 62M09; 60F05; 62H12

Less is More Revisited

Authors: Nobuko Yoshida, Ping Hou

Abstract: Multiparty session types (MPST) provide a type discipline where a programmer or architect specifies a whole view of communications as a global protocol, and each distributed program is locally type-checked against its end-point projection. After 10 years from the birth of MPST, Scalas and Yoshida discovered that the proofs of type safety in the literature which use the end-point projection with me… ▽ More Multiparty session types (MPST) provide a type discipline where a programmer or architect specifies a whole view of communications as a global protocol, and each distributed program is locally type-checked against its end-point projection. After 10 years from the birth of MPST, Scalas and Yoshida discovered that the proofs of type safety in the literature which use the end-point projection with mergeability are flawed. After this paper, researchers wrongly believed that the end-point projection (with mergeability) was unsound. We correct this misunderstanding, proposing a new general proof technique for type soundness of multiparty session $π$-calculus, which uses an association relation between a global type and its end-point projection. △ Less

Submitted 2 May, 2024; v1 submitted 26 February, 2024; originally announced February 2024.

Quasi-likelihood analysis for adaptive estimation of a degenerate diffusion process

Authors: Arnaud Gloter, Nakahiro Yoshida

Abstract: The adaptive quasi-likelihood analysis is developed for a degenerate diffusion process. Asymptotic normality and moment convergence are proved for the quasi-maximum likelihood estimators and quasi-Bayesian estimators, in the adaptive scheme. The adaptive quasi-likelihood analysis is developed for a degenerate diffusion process. Asymptotic normality and moment convergence are proved for the quasi-maximum likelihood estimators and quasi-Bayesian estimators, in the adaptive scheme. △ Less

Submitted 7 June, 2024; v1 submitted 23 February, 2024; originally announced February 2024.

A merger-driven scenario for clumpy galaxy formation in the epoch of reionization: Physical properties of clumps in the FirstLight simulation

Authors: Yurina Nakazato, Daniel Ceverino, Naoki Yoshida

Abstract: Recent JWST observations with superb angular resolution have revealed the existence of clumpy galaxies at high redshift through the detection of rest-frame optical emission lines. We use the FirstLight simulation to study the properties of (sub-)galactic clumps that are bright in [OIII] 5007$\mathrm{\mathring{A}}$ line with flux greater than $\sim 10^{-18} \, {\rm erg\, s^{-1}\, cm^{-2}}$, to be d… ▽ More Recent JWST observations with superb angular resolution have revealed the existence of clumpy galaxies at high redshift through the detection of rest-frame optical emission lines. We use the FirstLight simulation to study the properties of (sub-)galactic clumps that are bright in [OIII] 5007$\mathrm{\mathring{A}}$ line with flux greater than $\sim 10^{-18} \, {\rm erg\, s^{-1}\, cm^{-2}}$, to be detected by JWST. For 62 simulated galaxies that have stellar masses of $(0.5-6) \times 10^{10} \, M_\odot$ at $z=5$, we find clumps in 1828 snapshots in the redshift range $z = 9.5-5.5$. The clumps are identified by the surface density of star formation rate. About one-tenth of the snapshots show the existence of clumpy systems with two or more components. Most of the clumps are formed by mergers and can be characterized by their ages; central clumps dominated by stellar populations older than 50 Myr, and off-centered clumps dominated by younger stellar populations with specific star formation rates of $\sim 50 \, {\rm Gyr^{-1}}$. The latter type of young clumps is formed from gas debris in the tidal tails of major mergers with baryonic mass ratios of $1 \leq q < 4$. The merger-induced clumps are short-lived, and merge within a dynamical time of several tens million years. The number density of the clumpy systems is estimated to be $\sim 10^{-5}\, {\rm cMpc^{-3}}$, which is large enough to be detected in recent JWST surveys. △ Less

Submitted 26 June, 2024; v1 submitted 13 February, 2024; originally announced February 2024.

Comments: revised, 18 pages, 9 figures

Three Subtyping Algorithms for Binary Session Types and their Complexity Analyses (full version)

Authors: Thien Udomsrirungruang, Nobuko Yoshida

Abstract: Session types are a type discipline for describing and specifying communication behaviours of concurrent processes. Session subtyping, firstly introduced by Gay and Hole, is widely used for enlarging typability of session programs. This paper gives the complexity analysis of three algorithms for subtyping of synchronous binary session types. First, we analyse the complexity of the algorithm from t… ▽ More Session types are a type discipline for describing and specifying communication behaviours of concurrent processes. Session subtyping, firstly introduced by Gay and Hole, is widely used for enlarging typability of session programs. This paper gives the complexity analysis of three algorithms for subtyping of synchronous binary session types. First, we analyse the complexity of the algorithm from the original paper, which is based on an inductive tree search. We then introduce its optimised version, which improves the complexity, but is still exponential against the size of the two types. Finally, we propose a new quadratic algorithm based on a graph search using the concept of $\mathcal{XYZW}$-simulation, recently introduced by Silva et al. △ Less

Submitted 15 May, 2024; v1 submitted 10 February, 2024; originally announced February 2024.

Comments: 14 pages, 5 figures. Full version of a paper submitted to PLACES 2024

Modeling the core-halo mass relation in FDM halos

Authors: Hiroki Kawai, Ayuki Kamada, Kohei Kamada, Naoki Yoshida

Abstract: Fuzzy dark matter (FDM) is an intriguing candidate alternative to the standard cold dark matter (CDM). The FDM model predicts that dark halos have characteristic core structures generated by the effect of quantum pressure, which is different from the structure of CDM halos. We devise a semi-analytic model of a FDM halo density profile by assuming that the density distribution results from the redi… ▽ More Fuzzy dark matter (FDM) is an intriguing candidate alternative to the standard cold dark matter (CDM). The FDM model predicts that dark halos have characteristic core structures generated by the effect of quantum pressure, which is different from the structure of CDM halos. We devise a semi-analytic model of a FDM halo density profile by assuming that the density distribution results from the redistribution of mass in a halo with the Navarro-Frenk-White profile. We calculate the mass redistribution radius by considering dynamical relaxation within the FDM halo. We adopt a concentration-halo mass relation with lower concentration compared to that in the CDM model below the half mode mass, which originates from the suppressed matter density fluctuations at small length scales. Our model reproduces the core-halo mass relation (CHMR) found in the numerical simulation of \citet{2014NatPh..10..496S} at $z<1$. We show that the CHMR is well described by a double power law, unlike previous studies that approximate it by a single power law. Our model predictions are in reasonable agreement with the results of the largest FDM simulation of \citet{2021MNRAS.506.2603M} at $z=3$. We find that the core mass for a given halo mass follows the log-normal distribution, both in our model and in the simulation results for the first time, and quantitatively compare the variance of the distribution among them. Although our model does not fully explain the scatter of the CHMR, we show the scatter of the concentration-halo mass relation sizably contributes to them. △ Less

Submitted 13 June, 2024; v1 submitted 17 December, 2023; originally announced December 2023.

Comments: 22 pages, 17 figures, accepted by PRD

Report number: RESCEU-20/23

Checkpoint-based rollback recovery in session programming

Authors: Claudio Antares Mezzina, Francesco Tiezzi, Nobuko Yoshida

Abstract: To react to unforeseen circumstances or amend abnormal situations in communication-centric systems, programmers are in charge of "undoing" the interactions which led to an undesired state. To assist this task, session-based languages can be endowed with reversibility mechanisms. In this paper we propose a language enriched with programming facilities to commit session interactions, to roll back th… ▽ More To react to unforeseen circumstances or amend abnormal situations in communication-centric systems, programmers are in charge of "undoing" the interactions which led to an undesired state. To assist this task, session-based languages can be endowed with reversibility mechanisms. In this paper we propose a language enriched with programming facilities to commit session interactions, to roll back the computation to a previous commit point, and to abort the session. Rollbacks in our language always bring the system to previous visited states and a rollback cannot bring the system back to a point prior to the last commit. Programmers are relieved from the burden of ensuring that a rollback never restores a checkpoint imposed by a session participant different from the rollback requester. Such undesired situations are prevented at design-time (statically) by relying on a decidable compliance check at the type level, implemented in MAUDE. We show that the language satisfies error-freedom and progress of a session. △ Less

Submitted 3 December, 2023; originally announced December 2023.

Three-Dimensional Reconstruction of Weak-Lensing Mass Maps with a Sparsity Prior. II. Weighing Triaxial Cluster Halos

Authors: Shouzhuo Yang, Xiangchong Li, Naoki Yoshida

Abstract: Continuing work presented in Li et al. (2021), we performed a series of tests to our high-resolution three-dimensional mass map reconstruction algorithm \splinv{}. We test the mass reconstruction accuracy against realistic mock catalogs generated using shear field produced by triaxial halos with the inner density profile of $ρ\propto r^{-1}$ and of $ρ\propto r^{-1.5}$. The galaxy shape noise is mo… ▽ More Continuing work presented in Li et al. (2021), we performed a series of tests to our high-resolution three-dimensional mass map reconstruction algorithm \splinv{}. We test the mass reconstruction accuracy against realistic mock catalogs generated using shear field produced by triaxial halos with the inner density profile of $ρ\propto r^{-1}$ and of $ρ\propto r^{-1.5}$. The galaxy shape noise is modeled based on the Year-1 Subaru Hyper Suprime-Cam (HSC) Survey. After reviewing mathematical details of our algorithm and dark matter halo models, we determine an optimal value of the coefficient of the adaptive LASSO regression penalty term for single halo reconstruction. We successfully measure halo masses for massive triaxial halos; the mass determination accuracy is 5 percent for halos with $M = 10^{14.6}~M_\odot$ at $0.0625\leq z \leq 0.2425$, and 5 percent for those with $10^{14.8}~M_\odot$ at $0.0625\leq z \leq 0.4675$, and 20 percent for $M= 10^{15.0} ~M_\odot$ and $M=10^{15.2}~M_\odot$ in the redshift range $0.0625\leq z \leq 0.4675$. The redshift estimate accuracy is consistently below $Δz /z \leq 0.05$ for the above halo masses in the range $0.1525\leq z \leq 0.4675$. We also demonstrate that the orientation of triaxial halos and systematic error in our halo model do not affect reconstruction result significantly. Finally, we present results from reconstruction of mass distribution using shear catalogs produced by multiple halos, to show \splinv{}'s capability using realistic shear maps from ongoing and future galaxy surveys. △ Less

Submitted 30 November, 2023; originally announced December 2023.

Comments: 20 pages, 20 figures

Crash-Stop Failures in Asynchronous Multiparty Session Types

Authors: Adam D. Barwell, Ping Hou, Nobuko Yoshida, Fangyi Zhou

Abstract: Session types provide a typing discipline for message-passing systems. However, their theory often assumes an ideal world: one in which everything is reliable and without failures. Yet this is in stark contrast with distributed systems in the real world. To address this limitation, we introduce a new asynchronous multiparty session types (MPST) theory with crash-stop failures, where processes may… ▽ More Session types provide a typing discipline for message-passing systems. However, their theory often assumes an ideal world: one in which everything is reliable and without failures. Yet this is in stark contrast with distributed systems in the real world. To address this limitation, we introduce a new asynchronous multiparty session types (MPST) theory with crash-stop failures, where processes may crash arbitrarily and cease to interact after crashing. We augment asynchronous MPST and processes with crash handling branches, and integrate crash-stop failure semantics into types and processes. Our approach requires no user-level syntax extensions for global types, and features a formalisation of global semantics, which captures complex behaviours induced by crashed/crash handling processes. Our new theory covers the entire spectrum, ranging from the ideal world of total reliability to entirely unreliable scenarios where any process may crash, using optional reliability assumptions. Under these assumptions, we demonstrate the sound and complete correspondence between global and local type semantics, which guarantee deadlock-freedom, protocol conformance, and liveness of well-typed processes by construction, even in the presence of crashes. △ Less

Submitted 24 November, 2023; v1 submitted 20 November, 2023; originally announced November 2023.

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

Appearance of similar triangles by certain operations on triangles

Authors: Hiroki Naka, Takahiko Fujita, Naohiro Yoshida

Abstract: In this paper, a theorem about similar triangles is proved. It shows that two small and four large triangles similar to the original triangle can appear if we choose well among several intersections of the perpendicular bisectors of the sides with perpendicular lines of sides passing through the vertices of the triangle. In this paper, a theorem about similar triangles is proved. It shows that two small and four large triangles similar to the original triangle can appear if we choose well among several intersections of the perpendicular bisectors of the sides with perpendicular lines of sides passing through the vertices of the triangle. △ Less

Submitted 9 November, 2023; originally announced November 2023.

MSC Class: 51M04; 51M25

The Supersonic Project: Lighting up the faint end of the JWST UV luminosity function

Authors: Claire E. Williams, William Lake, Smadar Naoz, Blakesley Burkhart, Tommaso Treu, Federico Marinacci, Yurina Nakazato, Mark Vogelsberger, Naoki Yoshida, Gen Chiaki, Yeou S. Chiou, Avi Chen

Abstract: The James Webb Space Telescope (JWST) is capable of probing extremely early eras of our Universe when the supersonic relative motions between dark matter and baryonic overdensities modulate structure formation ($z>\sim 10$). We study low-mass galaxy formation including this "stream velocity" using high resolution AREPO hydrodynamics simulations, and present theoretical predictions of the UV lumino… ▽ More The James Webb Space Telescope (JWST) is capable of probing extremely early eras of our Universe when the supersonic relative motions between dark matter and baryonic overdensities modulate structure formation ($z>\sim 10$). We study low-mass galaxy formation including this "stream velocity" using high resolution AREPO hydrodynamics simulations, and present theoretical predictions of the UV luminosity function (UVLF) and galaxy stellar mass function (GSMF) down to extremely faint and low mass galaxies ($M_{UV}>\sim-15$, $10^4M_\odot<=M_*<=10^8 M_\odot)$. We show that, although the stream velocity suppresses early star formation overall, it induces a short period of rapid star formation in some larger dwarfs, leading to an enhancement in the faint-end of the UVLF at $z=12$. We demonstrate that JWST observations are close to this enhanced regime, and propose that the UVLF may constitute an important probe of the stream velocity at high redshift for JWST and future observatories. △ Less

Submitted 15 December, 2023; v1 submitted 5 October, 2023; originally announced October 2023.

Comments: 13 pages, 7 figures. Accepted ApJL

Quantum algorithm for the Vlasov simulation of the large-scale structure formation with massive neutrinos

Authors: Koichi Miyamoto, Soichiro Yamazaki, Fumio Uchida, Kotaro Fujisawa, Naoki Yoshida

Abstract: Investigating the cosmological implication of the fact that neutrino has finite mass is of importance for fundamental physics. In particular, massive neutrino affects the formation of the large-scale structure (LSS) of the universe, and, conversely, observations of the LSS can give constraints on the neutrino mass. Numerical simulations of the LSS formation including massive neutrino along with co… ▽ More Investigating the cosmological implication of the fact that neutrino has finite mass is of importance for fundamental physics. In particular, massive neutrino affects the formation of the large-scale structure (LSS) of the universe, and, conversely, observations of the LSS can give constraints on the neutrino mass. Numerical simulations of the LSS formation including massive neutrino along with conventional cold dark matter is thus an important task. For this, calculating the neutrino distribution in the phase space by solving the Vlasov equation is a suitable approach, but it requires solving the PDE in the $(6+1)$-dimensional space and is thus computationally demanding: Configuring $n_\mathrm{gr}$ grid points in each coordinate and $n_t$ time grid points leads to $O(n_\mathrm{gr}^6)$ memory space and $O(n_tn_\mathrm{gr}^6)$ queries to the coefficients in the discretized PDE. We propose a quantum algorithm for this task. Linearizing the Vlasov equation by neglecting the relatively weak self-gravity of the neutrino, we perform the Hamiltonian simulation to produce quantum states that encode the phase space distribution of neutrino. We also propose a way to extract the power spectrum of the neutrino density perturbations as classical data from the quantum state by quantum amplitude estimation with accuracy $ε$ and query complexity of order $\widetilde{O}((n_\mathrm{gr} + n_t)/ε)$. Our method also reduces the space complexity to $O(\mathrm{polylog}(n_\mathrm{gr}/ε))$ in terms of the qubit number, while using quantum random access memories with $O(n_\mathrm{gr}^3)$ entries. As far as we know, this is the first quantum algorithm for the LSS simulation that outputs the quantity of practical interest with guaranteed accuracy. △ Less

Submitted 26 February, 2024; v1 submitted 3 October, 2023; originally announced October 2023.

Report number: RESCEU-22/23

Journal ref: Phys. Rev. Research 6, 013200 (2024)

Early Structure Formation from Primordial Density Fluctuations with a Blue, Tilted Power Spectrum: High-Redshift Galaxies

Authors: Shingo Hirano, Naoki Yoshida

Abstract: Recent observations by the James Webb Space Telescope (JWST) discovered unexpectedly abundant luminous galaxies at high redshift, posing possibly a severe challenge to popular galaxy formation models. We study early structure formation in a cosmological model with a blue, tilted power spectrum (BTPS) given by $P(k) \propto k^{m_{\rm s}}$ with $m_{\rm s} > 1$ at small length scales. We run a set of… ▽ More Recent observations by the James Webb Space Telescope (JWST) discovered unexpectedly abundant luminous galaxies at high redshift, posing possibly a severe challenge to popular galaxy formation models. We study early structure formation in a cosmological model with a blue, tilted power spectrum (BTPS) given by $P(k) \propto k^{m_{\rm s}}$ with $m_{\rm s} > 1$ at small length scales. We run a set of cosmological $N$-body simulations and derive the abundance of dark matter halos and galaxies under simplified assumptions on star formation efficiency. The enhanced small-scale power allows rapid nonlinear structure formation at $z>7$, and galaxies with stellar mass exceeding $10^{10}\,M_\odot$ can be formed by $z=9$. Because of frequent mergers, the structure of galaxies and galaxy groups appears clumpy. The BTPS model reproduces the observed stellar mass density at $z=7-9$, and thus eases the claimed tension between galaxy formation theory and recent JWST observations. The large-scale structure of the present-day Universe is largely unaffected by the modification of the small-scale power spectrum. We conduct a systematic study by varying the slope of the small-scale power spectrum to derive constraints on the BTPS model from a set of observations of high-redshift galaxies. △ Less

Submitted 24 January, 2024; v1 submitted 20 June, 2023; originally announced June 2023.

Comments: 10 pages, 7 figures, 1 table, accepted for publication in ApJ

The Supersonic Project: Star Formation in Early Star Clusters without Dark Matter

Authors: William Lake, Smadar Naoz, Federico Marinacci, Blakesley Burkhart, Mark Vogelsberger, Claire E. Williams, Yeou S. Chiou, Gen Chiaki, Yurina Nakazato, Naoki Yoshida

Abstract: The formation mechanism of globular clusters (GCs) has long been debated by astronomers. It was recently proposed that Supersonically Induced Gas Objects (SIGOs), which formed in the early Universe due to the supersonic relative motion of baryons and dark matter at recombination, could be the progenitors of early globular clusters. In order to become GCs, SIGOs must form stars relatively efficient… ▽ More The formation mechanism of globular clusters (GCs) has long been debated by astronomers. It was recently proposed that Supersonically Induced Gas Objects (SIGOs), which formed in the early Universe due to the supersonic relative motion of baryons and dark matter at recombination, could be the progenitors of early globular clusters. In order to become GCs, SIGOs must form stars relatively efficiently despite forming outside of dark matter halos. We investigate the potential for star formation in SIGOs using cosmological hydrodynamic simulations, including the aforementioned relative motions of baryons and dark matter, molecular hydrogen cooling in primordial gas clouds, and including explicit star formation. We find that SIGOs do form stars and that the nascent star clusters formed through this process are accreted by dark matter halos on short timescales (a few hundreds of Myr). Thus, SIGOs may be found as intact substructures within these halos, analogous to many present-day GCs. From this result, we conclude that SIGOs are capable of forming star clusters with similar properties to globular clusters in the early Universe and we discuss their detectablity by upcoming JWST surveys. △ Less

Submitted 18 September, 2023; v1 submitted 1 June, 2023; originally announced June 2023.

Comments: 12 pages, 5 figures

Designing Asynchronous Multiparty Protocols with Crash-Stop Failures

Authors: Adam D. Barwell, Ping Hou, Nobuko Yoshida, Fangyi Zhou

Abstract: Session types provide a typing discipline for message-passing systems. However, most session type approaches assume an ideal world: one in which everything is reliable and without failures. Yet this is in stark contrast with distributed systems in the real world. To address this limitation, we introduce Teatrino, a code generation toolchain that utilises asynchronous multiparty session types (MPST… ▽ More Session types provide a typing discipline for message-passing systems. However, most session type approaches assume an ideal world: one in which everything is reliable and without failures. Yet this is in stark contrast with distributed systems in the real world. To address this limitation, we introduce Teatrino, a code generation toolchain that utilises asynchronous multiparty session types (MPST) with crash-stop semantics to support failure handling protocols. We augment asynchronous MPST and processes with crash handling branches. Our approach requires no user-level syntax extensions for global types and features a formalisation of global semantics, which captures complex behaviours induced by crashed/crash handling processes. The sound and complete correspondence between global and local type semantics guarantees deadlock-freedom, protocol conformance, and liveness of typed processes in the presence of crashes. Our theory is implemented in the toolchain Teatrino, which provides correctness by construction. Teatrino extends the Scribble multiparty protocol language to generate protocol-conforming Scala code, using the Effpi concurrent programming library. We extend both Scribble and Effpi to support crash-stop behaviour. We demonstrate the feasibility of our methodology and evaluate Teatrino with examples extended from both session type and distributed systems literature. △ Less

Submitted 15 May, 2023; v1 submitted 10 May, 2023; originally announced May 2023.

Comments: ECOOP 2023

Reionization and the ISM/Stellar Origins with JWST and ALMA (RIOJA): The core of the highest redshift galaxy overdensity at $z = 7.88$ confirmed by NIRSpec/JWST

Authors: Takuya Hashimoto, Javier Álvarez-Márquez, Yoshinobu Fudamoto, Luis Colina, Akio K. Inoue, Yurina Nakazato, Daniel Ceverino, Naoki Yoshida, Luca Costantin, Yuma Sugahara, Alejandro Crespo Gómez, Carmen Blanco-Prieto, Ken Mawatari, Santiago Arribas, Rui Marques-Chaves, Miguel Pereira-Santaella, Tom J. L. C. Bakx, Masato Hagimoto, Takeshi Hashigaya, Hiroshi Matsuo, Yoichi Tamura, Mitsutaka Usui, Yi W. Ren

Abstract: The protoclusters in the epoch of reionization, traced by galaxies overdensity regions, are ideal laboratories for studying the process of stellar assembly and cosmic reionization. We present the spectroscopic confirmation of the core of the most distant protocluster at $z = 7.88$, A2744-z7p9OD, with the James Webb Space Telescope NIRSpec integral field unit spectroscopy. The core region includes… ▽ More The protoclusters in the epoch of reionization, traced by galaxies overdensity regions, are ideal laboratories for studying the process of stellar assembly and cosmic reionization. We present the spectroscopic confirmation of the core of the most distant protocluster at $z = 7.88$, A2744-z7p9OD, with the James Webb Space Telescope NIRSpec integral field unit spectroscopy. The core region includes as many as 4 galaxies detected in [OIII] 4960 Å and 5008 Å in a small area of $\sim 3\arcsec \times 3\arcsec$, corresponding to $\sim$ 11 kpc $\times$ 11 kpc, after the lensing magnification correction. Three member galaxies are also tentatively detected in dust continuum in Atacama Large Millimeter/submillimeter Array Band 6, which is consistent with their red ultraviolet continuum slopes, $β\sim -1.3$. The member galaxies have stellar masses in the range of log($M_{*}/M_{\rm \odot}$) $\sim 7.6-9.2$ and star formation rates of $\sim 3-50$ $M_{\rm \odot}$ yr$^{-1}$, showing a diversity in their properties. FirstLight cosmological simulations reproduce the physical properties of the member galaxies including the stellar mass, [OIII] luminosity, and dust-to-stellar mass ratio, and predict that the member galaxies are on the verge of merging in a few to several tens Myr to become a large galaxy with $M_{\rm *}\sim 6\times10^{9} M_{\rm \odot}$. The presence of a multiple merger and evolved galaxies in the core region of A2744-z7p9OD indicates that environmental effects are already at work 650 Myr after the Big Bang. △ Less

Submitted 15 September, 2023; v1 submitted 8 May, 2023; originally announced May 2023.

Comments: 12 pages, 4 figures, 1 table, published in The Astrophysical Journal Letters

Simulations of Protoplanetary Disk Dispersal: Stellar Mass Dependence of the Disk Lifetime

Authors: A. Komaki, S. Fukuhara, T. K. Suzuki, N. Yoshida

Abstract: Recent infrared and submillimeter observations suggest that the protoplanetary disk lifetime depends on the central stellar mass. The disk dispersal is thought to be driven by viscous accretion, magneto-hydrodynamics (MHD) winds, and photoevaporation by the central star. We perform a set of one-dimensional simulations of long-term disk evolution that include all the three processes. We vary the st… ▽ More Recent infrared and submillimeter observations suggest that the protoplanetary disk lifetime depends on the central stellar mass. The disk dispersal is thought to be driven by viscous accretion, magneto-hydrodynamics (MHD) winds, and photoevaporation by the central star. We perform a set of one-dimensional simulations of long-term disk evolution that include all the three processes. We vary the stellar mass in the range of 0.5-7M$_{\odot}$, and study the mass dependence of the disk evolution. We show that a significant fraction of the disk gas is lost by MHD winds in the early stage, but the later disk evolution is mainly governed by photoevaporation. The disk radius decreases as photoevaporation clears out the gas in the outer disk efficiently. The qualitative evolutionary trends of the disk mass are remarkably similar for the wide range of the central stellar mass we consider, and the time evolution of the disk mass can be well fitted by a simple function. The dispersal time is approximately ten million years for low mass stars with weak mass dependence, but gets as short as two million years around a 7M$_{\odot}$ star. In the latter case, a prominent inner hole is formed by the combined effect of accretion and MHD winds within about one million years. The strength of the MHD wind and viscous accretion controls the overall mass-loss rate, but does not alter the dependence of the dispersal timescale on the central stellar mass. △ Less

Submitted 26 April, 2023; originally announced April 2023.

Comments: 14 pages, 10 figures, 1 table

Fabrication of a 64-Pixel TES Microcalorimeter Array with Iron Absorbers Uniquely Designed for 14.4-keV Solar Axion Search

Authors: Yuta Yagi, Tasuku Hayashi, Keita Tanaka, Rikuta Miyagawa, Ryo Ota, Noriko Y. Yamasaki, Kazuhisa Mitsuda, Nao Yoshida, Mikiko Saito, Takayuki Homma

Abstract: If a hypothetical elementary particle called an axion exists, to solve the strong CP problem, a 57Fe nucleus in the solar core could emit a 14.4-keV monochromatic axion through the M1 transition. If such axions are once more transformed into photons by a 57Fe absorber, a transition edge sensor (TES) X-ray microcalorimeter should be able to detect them efficiently. We have designed and fabricated a… ▽ More If a hypothetical elementary particle called an axion exists, to solve the strong CP problem, a 57Fe nucleus in the solar core could emit a 14.4-keV monochromatic axion through the M1 transition. If such axions are once more transformed into photons by a 57Fe absorber, a transition edge sensor (TES) X-ray microcalorimeter should be able to detect them efficiently. We have designed and fabricated a dedicated 64-pixel TES array with iron absorbers for the solar axion search. In order to decrease the effect of iron magnetization on spectroscopic performance, the iron absorber is placed next to the TES while maintaining a certain distance. A gold thermal transfer strap connects them. We have accomplished the electroplating of gold straps with high thermal conductivity. The residual resistivity ratio (RRR) was over 23, more than eight times higher than a previous evaporated strap. In addition, we successfully electroplated pure-iron films of more than a few micrometers in thickness for absorbers and a fabricated 64-pixel TES calorimeter structure. △ Less

Submitted 19 April, 2023; originally announced April 2023.

Comments: 5 pages, 5 figures, published in IEEE Transactions on Applied Superconductivity on 8 March 2023

Asymptotic expansion for batched bandits

Authors: Yechan Park, Nakahiro Yoshida

Abstract: In bandit algorithms, the randomly time-varying adaptive experimental design makes it difficult to apply traditional limit theorems to off-policy evaluation of the treatment effect. Moreover, the normal approximation by the central limit theorem becomes unsatisfactory for lack of information due to the small sample size of the inferior arm. To resolve this issue, we introduce a backwards asymptoti… ▽ More In bandit algorithms, the randomly time-varying adaptive experimental design makes it difficult to apply traditional limit theorems to off-policy evaluation of the treatment effect. Moreover, the normal approximation by the central limit theorem becomes unsatisfactory for lack of information due to the small sample size of the inferior arm. To resolve this issue, we introduce a backwards asymptotic expansion method and prove the validity of this scheme based on the partial mixing, that was originally introduced for the expansion of the distribution of a functional of a jump-diffusion process in a random environment. The theory is generalized in this paper to incorporate the backward propagation of random functions in the bandit algorithm. Besides the analytical validation, the simulation studies also support the new method. Our formulation is general and applicable to nonlinearly parametrized differentiable statistical models having an adaptive design. △ Less

Submitted 9 April, 2023; originally announced April 2023.

Quantum algorithm for collisionless Boltzmann simulation of self-gravitating systems

Authors: Soichiro Yamazaki, Fumio Uchida, Kotaro Fujisawa, Koichi Miyamoto, Naoki Yoshida

Abstract: The collisionless Boltzmann equation (CBE) is a fundamental equation that governs the dynamics of a broad range of astrophysical systems from space plasma to star clusters and galaxies. It is computationally expensive to integrate the CBE directly in a multi-dimensional phase space, and thus the applications to realistic astrophysical problems have been limited so far. Recently, Todorova & Steijl… ▽ More The collisionless Boltzmann equation (CBE) is a fundamental equation that governs the dynamics of a broad range of astrophysical systems from space plasma to star clusters and galaxies. It is computationally expensive to integrate the CBE directly in a multi-dimensional phase space, and thus the applications to realistic astrophysical problems have been limited so far. Recently, Todorova & Steijl (2020) proposed an efficient quantum algorithm to solve the CBE with significantly reduced computational complexity. We extend the algorithm to perform quantum simulations of self-gravitating systems, incorporating the method to calculate gravity with the major Fourier modes of the density distribution extracted from the solution-encoding quantum state. Our method improves the dependency of time and space complexities on Nv , the number of grid points in each velocity coordinate, compared to the classical simulation methods. We then conduct some numerical demonstrations of our method. We first run a 1+1 dimensional test calculation of free streaming motion on 64*64 grids using 13 simulated qubits and validate our method. We then perform simulations of Jeans collapse, and compare the result with analytic and linear theory calculations. It will thus allow us to perform large-scale CBE simulations on future quantum computers. △ Less

Submitted 23 May, 2024; v1 submitted 29 March, 2023; originally announced March 2023.

Comments: 12 pages, 10 figures

Machine Learning for Observational Cosmology

Authors: Kana Moriwaki, Takahiro Nishimichi, Naoki Yoshida

Abstract: An array of large observational programs using ground-based and space-borne telescopes is planned in the next decade. The forthcoming wide-field sky surveys are expected to deliver a sheer volume of data exceeding an exabyte. Processing the large amount of multiplex astronomical data is technically challenging, and fully automated technologies based on machine learning and artificial intelligence… ▽ More An array of large observational programs using ground-based and space-borne telescopes is planned in the next decade. The forthcoming wide-field sky surveys are expected to deliver a sheer volume of data exceeding an exabyte. Processing the large amount of multiplex astronomical data is technically challenging, and fully automated technologies based on machine learning and artificial intelligence are urgently needed. Maximizing scientific returns from the big data requires community-wide efforts. We summarize recent progress in machine learning applications in observational cosmology. We also address crucial issues in high-performance computing that are needed for the data processing and statistical analysis. △ Less

Submitted 25 May, 2023; v1 submitted 28 March, 2023; originally announced March 2023.

Comments: 55 pages, 8 figures, accepted for publication in Reports on Progress in Physics

Report number: YITP-23-40

An investigation of licensing of datasets for machine learning based on the GQM model

Authors: Junyu Chen, Norihiro Yoshida, Hiroaki Takada

Abstract: Dataset licensing is currently an issue in the development of machine learning systems. And in the development of machine learning systems, the most widely used are publicly available datasets. However, since the images in the publicly available dataset are mainly obtained from the Internet, some images are not commercially available. Furthermore, developers of machine learning systems do not ofte… ▽ More Dataset licensing is currently an issue in the development of machine learning systems. And in the development of machine learning systems, the most widely used are publicly available datasets. However, since the images in the publicly available dataset are mainly obtained from the Internet, some images are not commercially available. Furthermore, developers of machine learning systems do not often care about the license of the dataset when training machine learning models with it. In summary, the licensing of datasets for machine learning systems is in a state of incompleteness in all aspects at this stage. Our investigation of two collection datasets revealed that most of the current datasets lacked licenses, and the lack of licenses made it impossible to determine the commercial availability of the datasets. Therefore, we decided to take a more scientific and systematic approach to investigate the licensing of datasets and the licensing of machine learning systems that use the dataset to make it easier and more compliant for future developers of machine learning systems. △ Less

Submitted 23 March, 2023; originally announced March 2023.

The 300 pc resolution imaging of a z = 8.31 galaxy: Turbulent ionized gas and potential stellar feedback 600 million years after the Big Bang

Authors: Yoichi Tamura, Tom J. L. C. Bakx, Akio K. Inoue, Takuya Hashimoto, Tsuyoshi Tokuoka, Chihiro Imamura, Bunyo Hatsukade, Minju M. Lee, Kana Moriwaki, Takashi Okamoto, Kazuaki Ota, Hideki Umehata, Naoki Yoshida, Erik Zackrisson, Masato Hagimoto, Hiroshi Matsuo, Ikkoh Shimizu, Yuma Sugahara, Tsutomu T. Takeuchi

Abstract: We present the results of 300 pc resolution ALMA imaging of the [OIII] 88 $μ$m line and dust continuum emission from a $z = 8.312$ Lyman break galaxy MACS0416_Y1. The velocity-integrated [OIII] emission has three peaks which are likely associated with three young stellar clumps of MACS0416_Y1, while the channel map shows a complicated velocity structure with little indication of a global velocity… ▽ More We present the results of 300 pc resolution ALMA imaging of the [OIII] 88 $μ$m line and dust continuum emission from a $z = 8.312$ Lyman break galaxy MACS0416_Y1. The velocity-integrated [OIII] emission has three peaks which are likely associated with three young stellar clumps of MACS0416_Y1, while the channel map shows a complicated velocity structure with little indication of a global velocity gradient unlike what was found in [CII] 158 $μ$m at a larger scale, suggesting random bulk motion of ionized gas clouds inside the galaxy. In contrast, dust emission appears as two individual clumps apparently separating or bridging the [OIII]/stellar clumps. The cross correlation coefficient between dust and ultraviolet-related emission (i.e., [OIII] and ultraviolet continuum) is unity on a galactic scale, while it drops at < 1 kpc, suggesting well mixed geometry of multi-phase interstellar media on sub-kpc scales. If the cutoff scale characterizes different stages of star formation, the cutoff scale can be explained by gravitational instability of turbulent gas. We also report on a kpc-scale off-center cavity embedded in the dust continuum image. This could be a superbubble producing galactic-scale outflows, since the energy injection from the 4 Myr starburst suggested by a spectral energy distribution analysis is large enough to push the surrounding media creating a kpc-scale cavity. △ Less

Submitted 16 July, 2023; v1 submitted 20 March, 2023; originally announced March 2023.

Comments: 18 pages, 6 figures, 1 table, published in ApJ

Detections of [C II] 158 $μ$m and [O III] 88 $μ$m in a Local Lyman Continuum Emitter, Mrk 54, and its Implications to High-redshift ALMA Studies

Authors: Ryota Ura, Takuya Hashimoto, Akio K. Inoue, Dario Fadda, Matthew Hayes, Johannes Puschnig, Erik Zackrisson, Yoichi Tamura, Hiroshi Matsuo, Ken Mawatari, Yoshinobu Fudamoto, Masato Hagimoto, Nario Kuno, Yuma Sugahara, Satoshi Yamanaka, Tom J. L. C. Bakx, Yurina Nakazato, Mitsutaka Usui, Hidenobu Yajima, Naoki Yoshida

Abstract: We present integral field, far-infrared (FIR) spectroscopy of Mrk 54, a local Lyman Continuum Emitter (LCE), obtained with FIFI-LS on the Stratospheric Observatory for Infrared Astronomy. This is only the second time, after Haro 11, that [C II] 158 $μ$m and [O III] 88 $μ$m spectroscopy of the known LCEs have been obtained. We find that Mrk 54 has a strong [C II] emission that accounts for $\sim1$%… ▽ More We present integral field, far-infrared (FIR) spectroscopy of Mrk 54, a local Lyman Continuum Emitter (LCE), obtained with FIFI-LS on the Stratospheric Observatory for Infrared Astronomy. This is only the second time, after Haro 11, that [C II] 158 $μ$m and [O III] 88 $μ$m spectroscopy of the known LCEs have been obtained. We find that Mrk 54 has a strong [C II] emission that accounts for $\sim1$% of the total FIR luminosity, whereas it has only moderate [O III] emission, resulting in the low [O III]/[C II] luminosity ratio of $0.22\pm0.06$. In order to investigate whether [O III]/[C II] is a useful tracer of $f_{\rm esc}$ (LyC escape fraction), we examine the correlations of [O III]/[C II] and (i) the optical line ratio of $\rm O_{32} \equiv$ [O III] 5007 Å/[O II] 3727 Å, (ii) specific star formation rate, (iii) [O III] 88 $μ$m/[O I] 63 $μ$m ratio, (iv) gas phase metallicity, and (v) dust temperature based on a combined sample of Mrk 54 and the literature data from the Herschel Dwarf Galaxy Survey and the LITTLE THINGS Survey. We find that galaxies with high [O III]/[C II] luminosity ratios could be the result of high ionization (traced by $\rm O_{32}$), bursty star formation, high ionized-to-neutral gas volume filling factors (traced by [O III] 88 $μ$m/[O I] 63 $μ$m), and low gas-phase metallicities, which is in agreement with theoretical predictions. We present an empirical relation between the [O III]/[C II] ratio and $f_{\rm esc}$ based on the combination of the [O III]/[C II] and $\rm O_{32}$ correlation, and the known relation between $\rm O_{32}$ and $f_{\rm esc}$. The relation implies that high-redshift galaxies with high [O III]/[C II] ratios revealed by ALMA may have $f_{\rm esc}\gtrsim0.1$, significantly contributing to the cosmic reionization. △ Less

Submitted 20 March, 2023; originally announced March 2023.

Comments: 14 pages, 5 figures, Accepted for publication in ApJ

Upper Bound of Real Log Canonical Threshold of Tensor Decomposition and its Application to Bayesian Inference

Authors: Naoki Yoshida, Sumio Watanabe

Abstract: Tensor decomposition is now being used for data analysis, information compression, and knowledge recovery. However, the mathematical property of tensor decomposition is not yet fully clarified because it is one of singular learning machines. In this paper, we give the upper bound of its real log canonical threshold (RLCT) of the tensor decomposition by using an algebraic geometrical method and der… ▽ More Tensor decomposition is now being used for data analysis, information compression, and knowledge recovery. However, the mathematical property of tensor decomposition is not yet fully clarified because it is one of singular learning machines. In this paper, we give the upper bound of its real log canonical threshold (RLCT) of the tensor decomposition by using an algebraic geometrical method and derive its Bayesian generalization error theoretically. We also give considerations about its mathematical property through numerical experiments. △ Less

Submitted 3 April, 2023; v1 submitted 10 March, 2023; originally announced March 2023.

Strong depletion of $^{13}$C in CO induced by photolysis of CO$_{2}$ in the Martian atmosphere calculated by a photochemical model

Authors: Tatsuya Yoshida, Shohei Aoki, Yuichiro Ueno, Naoki Terada, Yuki Nakamura, Kimie Shiobara, Nao Yoshida, Hiromu Nakagawa, Shotaro Sakai, Shungo Koyama

Abstract: The isotopic signature of atmospheric carbon offers a unique tracer for the history of the Martian atmosphere and the origin of organic matter on Mars. Photolysis of CO$_{2}$ is known to induce strong isotopic fractionation of carbon between CO$_{2}$ and CO. However, its effect on the carbon isotopic compositions in the Martian atmosphere remains uncertain. Here we develop a 1-D photochemical mode… ▽ More The isotopic signature of atmospheric carbon offers a unique tracer for the history of the Martian atmosphere and the origin of organic matter on Mars. Photolysis of CO$_{2}$ is known to induce strong isotopic fractionation of carbon between CO$_{2}$ and CO. However, its effect on the carbon isotopic compositions in the Martian atmosphere remains uncertain. Here we develop a 1-D photochemical model considering isotopic fractionation via photolysis of CO$_{2}$ to estimate the vertical profiles of the carbon isotopic compositions of CO and CO$_{2}$ in the Martian atmosphere. We find that CO is depleted in $^{13}$C compared with CO$_{2}$ at each altitude due to the fractionation via CO$_{2}$ photolysis: the minimum value of $δ^{13}$C in CO is about $-170$ per mil under the standard eddy diffusion setting. This result supports the hypothesis that fractionated atmospheric CO is responsible for the production of the $^{13}$C-depleted organic carbon in Martian sediments detected by Curiosity Rover through the conversion of CO into organic materials and their deposition on the surface. The photolysis and transport-induced fractionation of CO we report here leads to a $\sim 15$ % decrease in the amount of inferred atmospheric loss when combined with the present-day fractionation of the atmosphere and previous studies of carbon escape to space. The fractionated isotopic composition of CO in the Martian atmosphere may be observed by ExoMars Trace Gas Orbiter (TGO) and ground-based telescopes, and escaping ion species produced by the fractionated carbon-bearing species may be detected by Martian Moons eXploration (MMX) in the future. △ Less

Submitted 24 February, 2023; originally announced February 2023.

Comments: 13 pages, 6 figures

Hybrid Multiparty Session Types -- Full Version

Authors: Lorenzo Gheri, Nobuko Yoshida

Abstract: Multiparty session types (MPST) are a specification and verification framework for distributed message-passing systems. The communication protocol of the system is specified as a global type, from which a collection of local types (local process implementations) is obtained by endpoint projection. A global type is a single disciplining entity for the whole system, specified by one designer that ha… ▽ More Multiparty session types (MPST) are a specification and verification framework for distributed message-passing systems. The communication protocol of the system is specified as a global type, from which a collection of local types (local process implementations) is obtained by endpoint projection. A global type is a single disciplining entity for the whole system, specified by one designer that has full knowledge of the communication protocol. On the other hand, distributed systems are often described in terms of their components: a different designer is in charge of providing a subprotocol for each component. The problem of modular specification of global protocols has been addressed in the literature, but the state of the art focuses only on dual input/output compatibility. Our work overcomes this limitation. We propose the first MPST theory of multiparty compositionality for distributed protocol specification that is semantics-preserving, allows the composition of two or more components, and retains full MPST expressiveness. We introduce hybrid types for describing subprotocols interacting with each other, define a novel compatibility relation, explicitly describe an algorithm for composing multiple subprotocols into a well-formed global type, and prove that compositionality preserves projection, thus retaining semantic guarantees, such as liveness and deadlock freedom. Finally, we test our work against real-world case studies and we smoothly extend our novel compatibility to MPST with delegation and explicit connections. △ Less

Submitted 20 March, 2023; v1 submitted 3 February, 2023; originally announced February 2023.

Existence of exact solution of the Susceptible-Exposed-Infectious-Recovered (SEIR) epidemic model

Authors: Norio Yoshida

Abstract: Exact solutions of the SEIR epidemic model are derived, and various properties of solutions are obtained directly from the exact solution. In this paper Abel differential equations play an important role in establishing the exact solution of SEIR differential system, in particular the number of infected individuals can be represented in a simple form by using a positive solution of an Abel differe… ▽ More Exact solutions of the SEIR epidemic model are derived, and various properties of solutions are obtained directly from the exact solution. In this paper Abel differential equations play an important role in establishing the exact solution of SEIR differential system, in particular the number of infected individuals can be represented in a simple form by using a positive solution of an Abel differential equation. It is shown that the parametric form of the exact solution satisfies some linear differential system including a positive solution of an Abel differential equation. △ Less

Submitted 1 February, 2023; originally announced February 2023.

Comments: 40 pages, 9 figures

MSC Class: 34A34

Journal ref: J. Differential Equations 355(2023), 103-143

Penalized estimation for non-identifiable models

Authors: Junichiro Yoshida, Nakahiro Yoshida

Abstract: We derive asymptotic properties of penalized estimators for singular models for which identifiability may break and the true parameter values can lie on the boundary of the parameter space. Selection consistency of the estimators is also validated. The problem that the true values lie on the boundary is dealt with by our previous results that are applicable to singular models, besides, penalized e… ▽ More We derive asymptotic properties of penalized estimators for singular models for which identifiability may break and the true parameter values can lie on the boundary of the parameter space. Selection consistency of the estimators is also validated. The problem that the true values lie on the boundary is dealt with by our previous results that are applicable to singular models, besides, penalized estimation and non-ergodic statistics. In order to overcome non-identifiability, we consider a suitable penalty such as the non-convex Bridge and the adaptive Lasso that stabilizes the asymptotic behavior of the estimator and shrinks inactive parameters. Then the estimator converges to one of the most parsimonious values among all the true values. In particular, the oracle property can also be obtained even if parametric structure of the singular model is so complex that likelihood ratio tests for model selection are labor intensive to perform. Among many potential applications, the examples handled in the paper are: (i) the superposition of parametric proportional hazard models and (ii) a counting process having intensity with multicollinear covariates. △ Less

Submitted 22 January, 2023; originally announced January 2023.

MSC Class: Primary 62F12; Secondary 62E20

Simulations of high-redshift [OIII] emitters: Chemical evolution and multi-line diagnostics

Authors: Yurina Nakazato, Naoki Yoshida, Daniel Ceverino

Abstract: Recent observations by James Webb Space Telescope discovered a number of high-redshift galaxies with strong emission lines from doubly ionized oxygen. Combined with ALMA observations of far-infrared lines, multi-line diagnostics can be applied to the high-redshift galaxies in order to probe the physical conditions of the inter-stellar medium. We study the formation and evolution of galaxies using… ▽ More Recent observations by James Webb Space Telescope discovered a number of high-redshift galaxies with strong emission lines from doubly ionized oxygen. Combined with ALMA observations of far-infrared lines, multi-line diagnostics can be applied to the high-redshift galaxies in order to probe the physical conditions of the inter-stellar medium. We study the formation and evolution of galaxies using the FirstLight simulation suite, which provides outputs of 62 high-resolution, zoom-in galaxy simulations. We devise a physical model of HII regions and calculate spatially resolved [OIII] line emission. We show that massive galaxies with stellar masses of $M_* > 10^9 M_\odot$ chemically evolve rapidly to $z=9$. Young stellar populations in the star-forming galaxies boost the [OIII] line emission, rendering the ratio of line luminosity to star formation rate larger than that for low-redshift galaxies, which is consistent with recent observations. Measuring the flux ratios of rest-frame optical and far-infrared lines allows us to estimate the physical conditions such as density and metallicity of the star-forming gas in high-redshift [OIII] emitters. △ Less

Submitted 29 June, 2023; v1 submitted 6 January, 2023; originally announced January 2023.

Comments: 10 pages, 5 figures, 1 table, ApJ in press

Quasi-maximum likelihood estimation and penalized estimation under non-standard conditions

Authors: Junichiro Yoshida, Nakahiro Yoshida

Abstract: The purpose of this article is to develop a general parametric estimation theory that allows the derivation of the limit distribution of estimators in non-regular models where the true parameter value may lie on the boundary of the parameter space or where even identifiability fails. For that, we propose a more general local approximation of the parameter space (at the true value) than previous st… ▽ More The purpose of this article is to develop a general parametric estimation theory that allows the derivation of the limit distribution of estimators in non-regular models where the true parameter value may lie on the boundary of the parameter space or where even identifiability fails. For that, we propose a more general local approximation of the parameter space (at the true value) than previous studies. This estimation theory is comprehensive in that it can handle penalized estimation as well as quasi-maximum likelihood estimation under such non-regular models. Besides, our results can apply to the so-called non-ergodic statistics, where the Fisher information is random in the limit, including the regular experiment that is locally asymptotically mixed normal. In penalized estimation, depending on the boundary constraint, even the Bridge estimator with $q<1$ does not necessarily give selection consistency. Therefore, some sufficient condition for selection consistency is described, precisely evaluating the balance between the boundary constraint and the form of the penalty. Examples handled in the paper are: (i) ML estimation of the generalized inverse Gaussian distribution, (ii) quasi-ML estimation of the diffusion parameter in a non-ergodic Itô process whose parameter space consists of positive semi-definite symmetric matrices, while the drift parameter is treated as nuisance and (iii) penalized ML estimation of variance components of random effects in linear mixed models. △ Less

Submitted 24 November, 2022; originally announced November 2022.

Comments: 35 pages

MSC Class: 62F12; 62E20

Some Martingale Properties of Simple Random Walk and Its Maximum Process

Authors: Takahiko Fujita, Shotaro Yagishita, Naohiro Yoshida

Abstract: In this paper, martingales related to simple random walks and their maximum process are investigated. First, a sufficient condition under which a function with three arguments, time, the random walk, and its maximum process becomes a martingale is presented, and as an application, an alternative way of deriving the Kennedy martingale is provided. Then, a complete characterization of a function wit… ▽ More In this paper, martingales related to simple random walks and their maximum process are investigated. First, a sufficient condition under which a function with three arguments, time, the random walk, and its maximum process becomes a martingale is presented, and as an application, an alternative way of deriving the Kennedy martingale is provided. Then, a complete characterization of a function with two arguments, the random walk and its maximum, being a martingale is presented. This martingale can be regarded as a discrete version of the Azéma--Yor martingale. As applications of discrete Azéma--Yor martingale, a proof of the Doob's inequalities is provided and a discrete Azéma--Yor solution for the Skorokhod embedding problem for the simple random walk is formulated and examined in detail. △ Less

Submitted 10 November, 2022; originally announced November 2022.

The Supersonic Project: The eccentricity and rotational support of SIGOs and DM GHOSts

Authors: Claire E. Williams, Smadar Naoz, William Lake, Yeou S. Chiou, Blakesley Burkhart, Federico Marinacci, Mark Vogelsberger, Gen Chiaki, Yurina Nakazato, Naoki Yoshida

Abstract: A supersonic relative velocity between dark matter (DM) and baryons (the stream velocity) at the time of recombination induces the formation of low mass objects with anomalous properties in the early Universe. We widen the scope of the `Supersonic Project' paper series to include objects we term Dark Matter + Gas Halos Offset by Streaming (DM GHOSts)--diffuse, DM-enriched structures formed because… ▽ More A supersonic relative velocity between dark matter (DM) and baryons (the stream velocity) at the time of recombination induces the formation of low mass objects with anomalous properties in the early Universe. We widen the scope of the `Supersonic Project' paper series to include objects we term Dark Matter + Gas Halos Offset by Streaming (DM GHOSts)--diffuse, DM-enriched structures formed because of a physical offset between the centers of mass of DM and baryonic overdensities. We present an updated numerical investigation of DM GHOSts and Supersonically Induced Gas Objects (SIGOs), including the effects of molecular cooling, in high resolution hydrodynamic simulations using the AREPO code. Supplemented by an analytical understanding of their ellipsoidal gravitational potentials, we study the population-level properties of these objects, characterizing their morphology, spin, radial mass, and velocity distributions in comparison to classical structures in non-streaming regions. The stream velocity causes deviations from sphericity in both the gas and DM components and lends greater rotational support to the gas. Low mass ($<\sim 10^{5.5}$ M$_\odot$) objects in regions of streaming demonstrate core-like rotation and mass profiles. Anomalies in the rotation and morphology of DM GHOSts could represent an early Universe analogue to observed ultra-faint dwarf galaxies with variations in DM content and unusual rotation curves. △ Less

Submitted 13 February, 2023; v1 submitted 3 November, 2022; originally announced November 2022.

Comments: 27 pages, 20 figures. Accepted ApJ

Exact solution of the Susceptible-Infectious-Recovered (SIR) epidemic model

Authors: Norio Yoshida

Abstract: Exact solution of the Susceptible-Infectious-Recovered (SIR) epidemic model is derived, and various properties of solution are obtained directly from the exact solution. It is shown that there exists an exact solution of an initial value problem for SIR differential system and that the parametric form of the exact solution satisfies some linear differential system. Exact solution of the Susceptible-Infectious-Recovered (SIR) epidemic model is derived, and various properties of solution are obtained directly from the exact solution. It is shown that there exists an exact solution of an initial value problem for SIR differential system and that the parametric form of the exact solution satisfies some linear differential system. △ Less

Submitted 2 October, 2022; originally announced October 2022.

Comments: 24 pages, 6 Postscript figures

MSC Class: 34A34

On the Expressiveness of Mixed Choice Sessions

Authors: Kirstin Peters, Nobuko Yoshida

Abstract: Session types provide a flexible programming style for structuring interaction, and are used to guarantee a safe and consistent composition of distributed processes. Traditional session types include only one-directional input (external) and output (internal) guarded choices. This prevents the session-processes to explore the full expressive power of the pi-calculus where the mixed choices are pro… ▽ More Session types provide a flexible programming style for structuring interaction, and are used to guarantee a safe and consistent composition of distributed processes. Traditional session types include only one-directional input (external) and output (internal) guarded choices. This prevents the session-processes to explore the full expressive power of the pi-calculus where the mixed choices are proved more expressive than the (non-mixed) guarded choices. To account this issue, recently Casal, Mordido, and Vasconcelos proposed the binary session types with mixed choices (CMV+). This paper carries a surprising, unfortunate result on CMV+: in spite of an inclusion of unrestricted channels with mixed choice, CMV+'s mixed choice is rather separate and not mixed. We prove this negative result using two methodologies (using either the leader election problem or a synchronisation pattern as distinguishing feature), showing that there exists no good encoding from the pi-calculus into CMV+, preserving distribution. We then close their open problem on the encoding from CMV+ into CMV (without mixed choice), proving its soundness and thereby that the encoding is good up to coupled similarity. △ Less

Submitted 6 September, 2022; originally announced September 2022.

Comments: In Proceedings EXPRESS/SOS 2022, arXiv:2208.14777. arXiv admin note: substantial text overlap with arXiv:2208.07041

ACM Class: F.1.2; F.3.2

Journal ref: EPTCS 368, 2022, pp. 113-130

Asymptotic expansion of an estimator for the Hurst coefficient

Authors: Yuliya Mishura, Hayate Yamagishi, Nakahiro Yoshida

Abstract: Asymptotic expansion is presented for an estimator of the Hurst coefficient of a fractional Brownian motion. For this, a recently developed theory of asymptotic expansion of the distribution of Wiener functionals is applied. The effects of the asymptotic expansion are demonstrated by numerical studies. Asymptotic expansion is presented for an estimator of the Hurst coefficient of a fractional Brownian motion. For this, a recently developed theory of asymptotic expansion of the distribution of Wiener functionals is applied. The effects of the asymptotic expansion are demonstrated by numerical studies. △ Less

Submitted 7 September, 2022; originally announced September 2022.

Comments: 34 pages, 10 figures