-
Long term variability of light-pollution in Bisei Town
Authors:
Ryosuke Itoh,
Syota Maeno
Abstract:
Bisei town, located in the west part of Japan, is known as a place where the local community protects its beautiful night sky from light pollution through its unique ordinances and the efforts of the local residents. It is also important to monitor in the quantity and quality of light pollution for precise measurement of astronomical observations. The fluorescent lamps in the city were gradually r…
▽ More
Bisei town, located in the west part of Japan, is known as a place where the local community protects its beautiful night sky from light pollution through its unique ordinances and the efforts of the local residents. It is also important to monitor in the quantity and quality of light pollution for precise measurement of astronomical observations. The fluorescent lamps in the city were gradually replaced with light emitting diode (LED) lamps. In order to investigate how much light pollution is affecting astronomical observation, we analyzed the archival photometric and spectroscopic data taken by the 101cm telescope that has been installed at Bisei Astronomical Observatory (BAO) since 2006. As a result, we found that there is no significant variability in sky brightness in optical bands, but from spectroscopic observation, we observed a blue humps around 4500 Åoriginating from LED lights from 2017 to 2023. The brightness of light pollution observed at BAO is not varied but the origin of light has gradually changed from fluorescent lamps to LED lamps.
△ Less
Submitted 20 December, 2023;
originally announced December 2023.
-
A measure-on-graph-valued diffusion: a particle system with collisions, and their applications
Authors:
Shuhei Mano
Abstract:
A diffusion taking value in probability measures on a graph with a vertex set $V$, $\sum_{i\in V}x_iδ_i$, is studied. The masses on each vertices satisfy the stochastic differential equation of the form $dx_i=\sum_{j\in N(i)}\sqrt{x_ix_j}dB_{ij}$ on the simplex, where $\{B_{ij}\}$ are independent standard Brownian motions with skew symmetry and $N(i)$ is the neighbour of the vertex $i$. A dual Mar…
▽ More
A diffusion taking value in probability measures on a graph with a vertex set $V$, $\sum_{i\in V}x_iδ_i$, is studied. The masses on each vertices satisfy the stochastic differential equation of the form $dx_i=\sum_{j\in N(i)}\sqrt{x_ix_j}dB_{ij}$ on the simplex, where $\{B_{ij}\}$ are independent standard Brownian motions with skew symmetry and $N(i)$ is the neighbour of the vertex $i$. A dual Markov chain on integer partitions to the Markov semigroup associated with the diffusion is used to show that the support of an extremal stationary state of the adjoint semigroup is an independent set of the graph. We also investigate the diffusion with a linear drift, which gives a killing of the dual Markov chain on a finite integer lattice. The Markov chain is used to study the unique stationary state of the diffusion, which generalizes the Dirichlet distribution. Two applications of the diffusions are discussed: analysis of an algorithm to find an independent set of a graph, and a Bayesian graph selection based on computation of probability of a sample by using coupling from the past.
△ Less
Submitted 10 March, 2023; v1 submitted 9 September, 2022;
originally announced September 2022.
-
Algebraic algorithm for direct sampling from toric models
Authors:
Shuhei Mano,
Nobuki Takayama
Abstract:
We show that Pfaffians or contiguity relations of hypergeometric functions of several variables give a direct sampling algorithm from toric models in statistics, which is a Markov chain on a lattice generated by a matrix $A$. A correspondence among graphical toric models and $A$-hypergeometric system is discussed and we give a sum formula of special values of $A$-hypergeometric polynomials. Some h…
▽ More
We show that Pfaffians or contiguity relations of hypergeometric functions of several variables give a direct sampling algorithm from toric models in statistics, which is a Markov chain on a lattice generated by a matrix $A$. A correspondence among graphical toric models and $A$-hypergeometric system is discussed and we give a sum formula of special values of $A$-hypergeometric polynomials. Some hypergeometric series which are interesting in view of statistics are presented.
△ Less
Submitted 28 October, 2021;
originally announced October 2021.
-
Asymptotic Moments Matching to Uniformly Minimum Variance Unbiased Estimation under Ewens Sampling Formula
Authors:
Masayo Y. Hirose,
Shuhei Mano
Abstract:
The Ewens sampling formula is a distribution related to the random partition of a positive integer. In this study, we investigate the issue of non-existence solutions in parameter estimation under the distribution. As a result, the first and second moments matching estimators to the uniformly minimum variance unbiased estimator are derived using the Ewens sampling formula in asymptotic sense. A Mo…
▽ More
The Ewens sampling formula is a distribution related to the random partition of a positive integer. In this study, we investigate the issue of non-existence solutions in parameter estimation under the distribution. As a result, the first and second moments matching estimators to the uniformly minimum variance unbiased estimator are derived using the Ewens sampling formula in asymptotic sense. A Monte Carlo simulation study is performed to evaluate the efficiency of the resulting estimators.
△ Less
Submitted 23 May, 2021;
originally announced May 2021.
-
Asymptotic bias reduction of maximum likelihood estimates via penalized likelihoods with differential geometry
Authors:
Masayo Y. Hirose,
Shuhei Mano
Abstract:
A procedure for asymptotic bias reduction of maximum likelihood estimates of generic estimands is developed. The estimator is realized as a plug-in estimator, where the parameter maximizes the penalized likelihood with a penalty function that satisfies a quasi-linear partial differential equation of the first order. The integration of the partial differential equation with the aid of differential…
▽ More
A procedure for asymptotic bias reduction of maximum likelihood estimates of generic estimands is developed. The estimator is realized as a plug-in estimator, where the parameter maximizes the penalized likelihood with a penalty function that satisfies a quasi-linear partial differential equation of the first order. The integration of the partial differential equation with the aid of differential geometry is discussed. Applications to generalized linear models, linear mixed-effects models, and a location-scale family are presented.
△ Less
Submitted 25 March, 2024; v1 submitted 30 November, 2020;
originally announced November 2020.
-
Characterization of the Nucleus, Morphology and Activity of Interstellar Comet 2I/Borisov by Optical and Near-Infrared GROWTH, Apache Point, IRTF, ZTF and Keck Observations
Authors:
Bryce T. Bolin,
Carey M. Lisse,
Mansi M. Kasliwal,
Robert Quimby,
Hanjie Tan,
Chris Copperwheat,
Zhong-Yi Lin,
Alessandro Morbidelli,
Lyu Abe,
Philippe Bendjoya,
James Bauer,
Kevin B. Burdge,
Michael Coughlin,
Christoffer Fremling,
Ryosuke Itoh,
Michael Koss,
Frank J. Masci,
Syota Maeno,
Eric E. Mamajek,
Federico Marocco,
Katsuhiro Murata,
Jean-Pierre Rivet,
Michael L. Sitko,
Daniel Stern,
David Vernet
, et al. (30 additional authors not shown)
Abstract:
We present visible and near-infrared photometric and spectroscopic observations of interstellar object 2I/Borisov taken from 2019 September 10 to 2019 November 29 using the GROWTH, the APO ARC 3.5 m and the NASA/IRTF 3.0 m combined with post and pre-discovery observations of 2I obtained by ZTF from 2019 March 17 to 2019 May 5. Comparison with imaging of distant Solar System comets shows an object…
▽ More
We present visible and near-infrared photometric and spectroscopic observations of interstellar object 2I/Borisov taken from 2019 September 10 to 2019 November 29 using the GROWTH, the APO ARC 3.5 m and the NASA/IRTF 3.0 m combined with post and pre-discovery observations of 2I obtained by ZTF from 2019 March 17 to 2019 May 5. Comparison with imaging of distant Solar System comets shows an object very similar to mildly active Solar System comets with an out-gassing rate of $\sim$10$^{27}$ mol/sec. The photometry, taken in filters spanning the visible and NIR range shows a gradual brightening trend of $\sim0.03$ mags/day since 2019 September 10 UTC for a reddish object becoming neutral in the NIR. The lightcurve from recent and pre-discovery data reveals a brightness trend suggesting the recent onset of significant H$_2$O sublimation with the comet being active with super volatiles such as CO at heliocentric distances $>$6 au consistent with its extended morphology. Using the advanced capability to significantly reduce the scattered light from the coma enabled by high-resolution NIR images from Keck adaptive optics taken on 2019 October 04, we estimate a diameter of 2I's nucleus of $\lesssim$1.4 km. We use the size estimates of 1I/'Oumuamua and 2I/Borisov to roughly estimate the slope of the ISO size-distribution resulting in a slope of $\sim$3.4$\pm$1.2, similar to Solar System comets and bodies produced from collisional equilibrium.
△ Less
Submitted 12 May, 2020; v1 submitted 30 October, 2019;
originally announced October 2019.
-
Application of the independent component analysis to the iKAGRA data
Authors:
KAGRA Collaboration,
T. Akutsu,
M. Ando,
K. Arai,
Y. Arai,
S. Araki,
A. Araya,
N. Aritomi,
H. Asada,
Y. Aso,
S. Atsuta,
K. Awai,
S. Bae,
Y. Bae,
L. Baiotti,
R. Bajpai,
M. A. Barton,
K. Cannon,
E. Capocasa,
M. Chan,
C. Chen,
K. Chen,
Y. Chen,
H. Chu,
Y-K. Chu
, et al. (227 additional authors not shown)
Abstract:
We apply the independent component analysis (ICA) to the real data from a gravitational wave detector for the first time. Specifically we use the iKAGRA data taken in April 2016, and calculate the correlations between the gravitational wave strain channel and 35 physical environmental channels. Using a couple of seismic channels which are found to be strongly correlated with the strain, we perform…
▽ More
We apply the independent component analysis (ICA) to the real data from a gravitational wave detector for the first time. Specifically we use the iKAGRA data taken in April 2016, and calculate the correlations between the gravitational wave strain channel and 35 physical environmental channels. Using a couple of seismic channels which are found to be strongly correlated with the strain, we perform ICA. Injecting a sinusoidal continuous signal in the strain channel, we find that ICA recovers correct parameters with enhanced signal-to-noise ratio, which demonstrates usefulness of this method. Among the two implementations of ICA used here, we find the correlation method yields the optimal result for the case environmental noises act on the strain channel linearly.
△ Less
Submitted 1 June, 2020; v1 submitted 8 August, 2019;
originally announced August 2019.
-
A nonparametric method to assess significance of events in search for gravitational waves with false discovery rate
Authors:
Hirotaka Yuzurihara,
Shuhei Mano,
Hideyuki Tagoshi
Abstract:
In this paper, we present a consistent procedure to assess the significance of gravitational wave events observed by laser interferometric gravitational wave detectors based on the background distribution of detection statistic. We propose a non-parametric method to estimate $p$-value. Based on the estimated $p$-values, we propose a new procedure to assess the significance of a particular event wi…
▽ More
In this paper, we present a consistent procedure to assess the significance of gravitational wave events observed by laser interferometric gravitational wave detectors based on the background distribution of detection statistic. We propose a non-parametric method to estimate $p$-value. Based on the estimated $p$-values, we propose a new procedure to assess the significance of a particular event with $q$-value which is the minimum false discovery rate that can be attained when calling the event significant. The $q$-value gives us a criterion on the significance of events which is different from $P_{\rm astro}$ which is used in the LIGO-Virgo analysis and in other analysis. The proposed procedure is applied to the 1-OGC and 2-OGC catalogs [2][3]. For most of the events which were claimed significant in [2] and [3], we also obtain the same results. However, there are differences in the significance for several marginal events. Since the proposed procedure does not require any assumptions on signal and noise, it is very simple and straightforward. The procedure is also applicable to other searches for gravitational waves whose background distribution of detection statistic is difficult to know.
△ Less
Submitted 12 October, 2021; v1 submitted 30 June, 2019;
originally announced July 2019.
-
First cryogenic test operation of underground km-scale gravitational-wave observatory KAGRA
Authors:
KAGRA Collaboration,
T. Akutsu,
M. Ando,
K. Arai,
Y. Arai,
S. Araki,
A. Araya,
N. Aritomi,
H. Asada,
Y. Aso,
S. Atsuta,
K. Awai,
S. Bae,
L. Baiotti,
M. A. Barton,
K. Cannon,
E. Capocasa,
C-S. Chen,
T-W. Chiu,
K. Cho,
Y-K. Chu,
K. Craig,
W. Creus,
K. Doi,
K. Eda
, et al. (179 additional authors not shown)
Abstract:
KAGRA is a second-generation interferometric gravitational-wave detector with 3-km arms constructed at Kamioka, Gifu in Japan. It is now in its final installation phase, which we call bKAGRA (baseline KAGRA), with scientific observations expected to begin in late 2019. One of the advantages of KAGRA is its underground location of at least 200 m below the ground surface, which brings small seismic…
▽ More
KAGRA is a second-generation interferometric gravitational-wave detector with 3-km arms constructed at Kamioka, Gifu in Japan. It is now in its final installation phase, which we call bKAGRA (baseline KAGRA), with scientific observations expected to begin in late 2019. One of the advantages of KAGRA is its underground location of at least 200 m below the ground surface, which brings small seismic motion at low frequencies and high stability of the detector. Another advantage is that it cools down the sapphire test mass mirrors to cryogenic temperatures to reduce thermal noise. In April-May 2018, we have operated a 3-km Michelson interferometer with a cryogenic test mass for 10 days, which was the first time that km-scale interferometer was operated at cryogenic temperatures. In this article, we report the results of this "bKAGRA Phase 1" operation. We have demonstrated the feasibility of 3-km interferometer alignment and control with cryogenic mirrors.
△ Less
Submitted 11 January, 2019;
originally announced January 2019.
-
KAGRA: 2.5 Generation Interferometric Gravitational Wave Detector
Authors:
T. Akutsu,
M. Ando,
K. Arai,
Y. Arai,
S. Araki,
A. Araya,
N. Aritomi,
H. Asada,
Y. Aso,
S. Atsuta,
K. Awai,
S. Bae,
L. Baiotti,
M. A. Barton,
K. Cannon,
E. Capocasa,
C-S. Chen,
T-W. Chiu,
K. Cho,
Y-K. Chu,
K. Craig,
W. Creus,
K. Doi,
K. Eda,
Y. Enomoto
, et al. (169 additional authors not shown)
Abstract:
The recent detections of gravitational waves (GWs) reported by LIGO/Virgo collaborations have made significant impact on physics and astronomy. A global network of GW detectors will play a key role to solve the unknown nature of the sources in coordinated observations with astronomical telescopes and detectors. Here we introduce KAGRA (former name LCGT; Large-scale Cryogenic Gravitational wave Tel…
▽ More
The recent detections of gravitational waves (GWs) reported by LIGO/Virgo collaborations have made significant impact on physics and astronomy. A global network of GW detectors will play a key role to solve the unknown nature of the sources in coordinated observations with astronomical telescopes and detectors. Here we introduce KAGRA (former name LCGT; Large-scale Cryogenic Gravitational wave Telescope), a new GW detector with two 3-km baseline arms arranged in the shape of an "L", located inside the Mt. Ikenoyama, Kamioka, Gifu, Japan. KAGRA's design is similar to those of the second generations such as Advanced LIGO/Virgo, but it will be operating at the cryogenic temperature with sapphire mirrors. This low temperature feature is advantageous for improving the sensitivity around 100 Hz and is considered as an important feature for the third generation GW detector concept (e.g. Einstein Telescope of Europe or Cosmic Explorer of USA). Hence, KAGRA is often called as a 2.5 generation GW detector based on laser interferometry. The installation and commissioning of KAGRA is underway and its cryogenic systems have been successfully tested in May, 2018. KAGRA's first observation run is scheduled in late 2019, aiming to join the third observation run (O3) of the advanced LIGO/Virgo network. In this work, we describe a brief history of KAGRA and highlights of main feature. We also discuss the prospects of GW observation with KAGRA in the era of O3. When operating along with the existing GW detectors, KAGRA will be helpful to locate a GW source more accurately and to determine the source parameters with higher precision, providing information for follow-up observations of a GW trigger candidate.
△ Less
Submitted 20 November, 2018;
originally announced November 2018.
-
A reversal phenomenon in estimation based on multiple samples from the Poisson--Dirichlet distribution
Authors:
Koji Tsukuda,
Shuhei Mano
Abstract:
Consider two forms of sampling from a population: (i) drawing $s$ samples of $n$ elements with replacement and (ii) drawing a single sample of $ns$ elements. In this paper, under the setting where the descending order population frequency follows the Poisson--Dirichlet distribution with parameter $θ$, we report that the magnitude relation of the Fisher information, which sample partitions converte…
▽ More
Consider two forms of sampling from a population: (i) drawing $s$ samples of $n$ elements with replacement and (ii) drawing a single sample of $ns$ elements. In this paper, under the setting where the descending order population frequency follows the Poisson--Dirichlet distribution with parameter $θ$, we report that the magnitude relation of the Fisher information, which sample partitions converted from samples (i) and (ii) possess, can change depending on the parameters, $n$, $s$, and $θ$. Roughly speaking, if $θ$ is small relative to $n$ and $s$, the Fisher information of (i) is larger than that of (ii); on the contrary, if $θ$ is large relative to $n$ and $s$, the Fisher information of (ii) is larger than that of (i). The result represents one aspect of random distributions.
△ Less
Submitted 2 February, 2018;
originally announced February 2018.
-
Construction of KAGRA: an Underground Gravitational Wave Observatory
Authors:
T. Akutsu,
M. Ando,
S. Araki,
A. Araya,
T. Arima,
N. Aritomi,
H. Asada,
Y. Aso,
S. Atsuta,
K. Awai,
L. Baiotti,
M. A. Barton,
D. Chen,
K. Cho,
K. Craig,
R. DeSalvo,
K. Doi,
K. Eda,
Y. Enomoto,
R. Flaminio,
S. Fujibayashi,
Y. Fujii,
M. -K. Fujimoto,
M. Fukushima,
T. Furuhata
, et al. (202 additional authors not shown)
Abstract:
Major construction and initial-phase operation of a second-generation gravitational-wave detector KAGRA has been completed. The entire 3-km detector is installed underground in a mine in order to be isolated from background seismic vibrations on the surface. This allows us to achieve a good sensitivity at low frequencies and high stability of the detector. Bare-bones equipment for the interferomet…
▽ More
Major construction and initial-phase operation of a second-generation gravitational-wave detector KAGRA has been completed. The entire 3-km detector is installed underground in a mine in order to be isolated from background seismic vibrations on the surface. This allows us to achieve a good sensitivity at low frequencies and high stability of the detector. Bare-bones equipment for the interferometer operation has been installed and the first test run was accomplished in March and April of 2016 with a rather simple configuration. The initial configuration of KAGRA is named {\it iKAGRA}. In this paper, we summarize the construction of KAGRA, including the study of the advantages and challenges of building an underground detector and the operation of the iKAGRA interferometer together with the geophysics interferometer that has been constructed in the same tunnel.
△ Less
Submitted 11 December, 2017; v1 submitted 30 November, 2017;
originally announced December 2017.
-
The status of KAGRA underground cryogenic gravitational wave telescope
Authors:
KAGRA Collaboration,
T. Akutsu,
M. Ando,
A. Araya,
N. Aritomi,
H. Asada,
Y. Aso,
S. Atsuta,
K. Awai,
M. A. Barton,
K. Cannon,
K. Craig,
W. Creus,
K. Doi,
K. Eda,
Y. Enomoto,
R. Flaminio,
Y. Fujii,
M. -K. Fujimoto,
T. Furuhata,
S. Haino,
K. Hasegawa,
K. Hashino,
K. Hayama,
S. Hirobayashi
, et al. (126 additional authors not shown)
Abstract:
KAGRA is a 3-km interferometric gravitational wave telescope located in the Kamioka mine in Japan. It is the first km-class gravitational wave telescope constructed underground to reduce seismic noise, and the first km-class telescope to use cryogenic cooling of test masses to reduce thermal noise. The construction of the infrastructure to house the interferometer in the tunnel, and the initial ph…
▽ More
KAGRA is a 3-km interferometric gravitational wave telescope located in the Kamioka mine in Japan. It is the first km-class gravitational wave telescope constructed underground to reduce seismic noise, and the first km-class telescope to use cryogenic cooling of test masses to reduce thermal noise. The construction of the infrastructure to house the interferometer in the tunnel, and the initial phase operation of the interferometer with a simple 3-km Michelson configuration have been completed. The first cryogenic operation is expected in 2018, and the observing runs with a full interferometer are expected in 2020s. The basic interferometer configuration and the current status of KAGRA are described.
△ Less
Submitted 13 October, 2017;
originally announced October 2017.
-
Partition structure and the A-hypergeometric distribution associated with the rational normal curve
Authors:
Shuhei Mano
Abstract:
A distribution whose normalization constant is an A-hypergeometric polynomial is called an A-hypergeometric distribution. Such a distribution is in turn a generalization of the generalized hypergeometric distribution on the contingency tables with fixed marginal sums. In this paper, we will see that an A-hypergeometric distribution with a homogeneous matrix of two rows, especially, that associated…
▽ More
A distribution whose normalization constant is an A-hypergeometric polynomial is called an A-hypergeometric distribution. Such a distribution is in turn a generalization of the generalized hypergeometric distribution on the contingency tables with fixed marginal sums. In this paper, we will see that an A-hypergeometric distribution with a homogeneous matrix of two rows, especially, that associated with the rational normal curve, appears in inferences involving exchangeable partition structures. An exact sampling algorithm is presented for the general (any number of rows) A-hypergeometric distributions. Then, the maximum likelihood estimation of the A-hypergeometric distribution associated with the rational normal curve, which is an algebraic exponential family, is discussed. The information geometry of the Newton polytope is useful for analyzing the full and the curved exponential family. Algebraic methods are provided for evaluating the A-hypergeometric polynomials.
△ Less
Submitted 26 October, 2017; v1 submitted 12 July, 2016;
originally announced July 2016.
-
The star-shaped Lambda-coalescent and Fleming-Viot process
Authors:
Robert Griffiths,
Shuhei Mano
Abstract:
The star-shaped $Λ$-coalescent and corresponding $Λ$-Fleming-Viot process where the $Λ$ measure has a single atom at unity are studied in this paper. The transition functions and stationary distribution of the $Λ$-Fleming-Viot process are derived in a two-type model with mutation. The distribution of the number of non-mutant lines back in time in the star-shaped $Λ$-coalescent is found. Extensions…
▽ More
The star-shaped $Λ$-coalescent and corresponding $Λ$-Fleming-Viot process where the $Λ$ measure has a single atom at unity are studied in this paper. The transition functions and stationary distribution of the $Λ$-Fleming-Viot process are derived in a two-type model with mutation. The distribution of the number of non-mutant lines back in time in the star-shaped $Λ$-coalescent is found. Extensions are made to a model with $d$ types, either with parent independent mutation or general Markov mutation, and an infinitely-many-types model when $d\to \infty$. An eigenfunction expansion for the transition functions is found which has polynomial right eigenfunctions and left eigenfunctions described by hyperfunctions. A further star-shaped model with general frequency dependent change is considered and the stationary distribution in the Fleming-Viot process derived. This model includes a star-shaped $Λ$-Fleming-Viot process with mutation and selection. In a general $Λ$-coalescent explicit formulae for the transition functions and stationary distribution when there is mutation are unknown, however in this paper explicit formulae are derived in the star-shaped coalescent.
△ Less
Submitted 8 May, 2016; v1 submitted 24 June, 2015;
originally announced June 2015.
-
Extreme sizes in the Gibbs-type exchangeable random partitions
Authors:
Shuhei Mano
Abstract:
Gibbs-type exchangeable random partitions, which is a class of multiplicative measures on the set of positive integer partitions, appear in various contexts, including Bayesian statistics, random combinatorial structures, and stochastic models of diversity in various phenomena. Some distributional results on ordered sizes in the Gibbs partition are established by introducing associated partial Bel…
▽ More
Gibbs-type exchangeable random partitions, which is a class of multiplicative measures on the set of positive integer partitions, appear in various contexts, including Bayesian statistics, random combinatorial structures, and stochastic models of diversity in various phenomena. Some distributional results on ordered sizes in the Gibbs partition are established by introducing associated partial Bell polynomials and analysis of the generating functions. The combinatorial approach is applied to derive explicit results on asymptotic behavior of the extreme sizes in the Gibbs partition. Especially, Ewens-Pitman partition, which is the sample from the Poisson-Dirichlet process and has been discussed from rather model-specific viewpoints, and a random partition which was recently introduced by Gnedin, are discussed in the details. As by-products, some formulas for the associated partial Bell polynomials are presented.
△ Less
Submitted 29 June, 2015; v1 submitted 9 June, 2013;
originally announced June 2013.
-
Kernel Approximate Bayesian Computation for Population Genetic Inferences
Authors:
Shigeki Nakagome,
Kenji Fukumizu,
Shuhei Mano
Abstract:
Approximate Bayesian computation (ABC) is a likelihood-free approach for Bayesian inferences based on a rejection algorithm method that applies a tolerance of dissimilarity between summary statistics from observed and simulated data. Although several improvements to the algorithm have been proposed, none of these improvements avoid the following two sources of approximation: 1) lack of sufficient…
▽ More
Approximate Bayesian computation (ABC) is a likelihood-free approach for Bayesian inferences based on a rejection algorithm method that applies a tolerance of dissimilarity between summary statistics from observed and simulated data. Although several improvements to the algorithm have been proposed, none of these improvements avoid the following two sources of approximation: 1) lack of sufficient statistics: sampling is not from the true posterior density given data but from an approximate posterior density given summary statistics; and 2) non-zero tolerance: sampling from the posterior density given summary statistics is achieved only in the limit of zero tolerance. The first source of approximation can be improved by adding a summary statistic, but an increase in the number of summary statistics could introduce additional variance caused by the low acceptance rate. Consequently, many researchers have attempted to develop techniques to choose informative summary statistics. The present study evaluated the utility of a kernel-based ABC method (Fukumizu et al. 2010, arXiv:1009.5736 and 2011, NIPS 24: 1549-1557) for complex problems that demand many summary statistics. Specifically, kernel ABC was applied to population genetic inference. We demonstrate that, in contrast to conventional ABCs, kernel ABC can incorporate a large number of summary statistics while maintaining high performance of the inference.
△ Less
Submitted 24 September, 2013; v1 submitted 14 May, 2012;
originally announced May 2012.
-
Duality between the two-locus Wright-Fisher Diffusion Model and the Ancestral Process with Recombination
Authors:
Shuhei Mano
Abstract:
Known results on the moments of the distribution generated by the two-locus Wright-Fisher diffusion model and a duality between the diffusion process and the ancestral process with recombination are briefly summarized. A numerical methods for computing moments by a Markov chain Monte Carlo and a method to compute closed-form expressions of the moments are presented. By using the duality argument p…
▽ More
Known results on the moments of the distribution generated by the two-locus Wright-Fisher diffusion model and a duality between the diffusion process and the ancestral process with recombination are briefly summarized. A numerical methods for computing moments by a Markov chain Monte Carlo and a method to compute closed-form expressions of the moments are presented. By using the duality argument properties of the ancestral recombination graph are studied in terms of the moments.
△ Less
Submitted 17 February, 2012; v1 submitted 26 January, 2012;
originally announced January 2012.
-
Ancestral Graph with Bias in Gene Conversion
Authors:
Shuhei Mano
Abstract:
Gene conversion is a mechanism by which a double-strand break in a DNA molecule is repaired using a homologous DNA molecule as a template. As a result, one gene is 'copied and pasted' onto the other gene. It was recently reported that the direction of gene conversion appears to be biased towards G and C nucleotides. In this paper a stochastic model of the dynamics of the bias in gene conversion is…
▽ More
Gene conversion is a mechanism by which a double-strand break in a DNA molecule is repaired using a homologous DNA molecule as a template. As a result, one gene is 'copied and pasted' onto the other gene. It was recently reported that the direction of gene conversion appears to be biased towards G and C nucleotides. In this paper a stochastic model of the dynamics of the bias in gene conversion is developed for a finite population of members in a multigene family. The dual process is the biased voter model, which generates an ancestral random graph for a given sample. An importance-sampling algorithm for computing the likelihood of the sample is also given.
△ Less
Submitted 20 January, 2012; v1 submitted 6 July, 2009;
originally announced July 2009.
-
Duality, Ancestral and Diffusion Processes in Models with Selection
Authors:
Shuhei Mano
Abstract:
The ancestral selection graph in population genetics was introduced by KroneNeuhauser (1997) as an analogue of the coalescent genealogy of a sample of genes from a neutrally evolving population. The number of particles in this graph, followed backwards in time, is a birth and death process with quadratic death and linear birth rates. In this paper an explicit form of the probability distribution…
▽ More
The ancestral selection graph in population genetics was introduced by KroneNeuhauser (1997) as an analogue of the coalescent genealogy of a sample of genes from a neutrally evolving population. The number of particles in this graph, followed backwards in time, is a birth and death process with quadratic death and linear birth rates. In this paper an explicit form of the probability distribution of the number of particles is obtained by using the density of the allele frequency in the corresponding diffusion model obtained by Kimura (1955). It is shown that the process of fixation of the allele in the diffusion model corresponds to convergence of the ancestral process to its stationary measure. The time to fixation of the allele conditional on fixation is studied in terms of the ancestral process.
△ Less
Submitted 2 February, 2009; v1 submitted 16 April, 2008;
originally announced April 2008.
-
Status of GRB Observations with the Suzaku Wideband All-sky Monitor
Authors:
M. S. Tashiro,
Y. Terada,
Y. Urata,
K. Onda,
N. Kodaka,
A. Endo,
M. Suzuki,
K. Morigami,
K. Yamaoka,
Y. E. Nakagawa,
S. Sugita,
Y. Fukazawa,
M. Ohno,
T. Takahashi,
C. Kira,
T. Uehara,
T. Tamagawa,
T. Enoto,
R. Miyawaki,
K. Nakazawa,
K. Makishima,
E. Sonoda,
M. Yamauchi,
S. Maeno,
H. Tanaka
, et al. (8 additional authors not shown)
Abstract:
The Wide-band All-sky Monitor (WAM) is a function of the large lateral BGO shield of the Hard X-ray Detector (HXD) onboard Suzaku. Its large geometrical area of 800 cm^2 per side, the large stopping power for the hard X-rays and the wide-field of view make the WAM an ideal detector for gamma-ray bursts (GRBs) observations in the energy range of 50-5000 keV. In fact, the WAM has observed 288 GRBs…
▽ More
The Wide-band All-sky Monitor (WAM) is a function of the large lateral BGO shield of the Hard X-ray Detector (HXD) onboard Suzaku. Its large geometrical area of 800 cm^2 per side, the large stopping power for the hard X-rays and the wide-field of view make the WAM an ideal detector for gamma-ray bursts (GRBs) observations in the energy range of 50-5000 keV. In fact, the WAM has observed 288 GRBs confirmed by other satellites, till the end of May 2007.
△ Less
Submitted 29 February, 2008;
originally announced February 2008.
-
Spectral evolution of GRB 060904A observed with Swift and Suzaku -- Possibility of Inefficient Electron Acceleration
Authors:
Daisuke Yonetoku,
Sachiko Tanabe,
Toshio Murakami,
Naomi Emura,
Yuka Aoyama,
Takashi Kidamura,
Hironobu Kodaira,
Yoshiki Kodama,
Ryota Kozaka,
Takuro Nashimoto,
Shinya Okuno,
Satoshi Yokota,
Satoru Yoshinari,
Keiichi Abe,
Kaori Onda,
Makoto S. Tashiro,
Yuji Urata,
Yujin E. Nakagawa,
Satoshi Sugita,
Kazutaka Yamaoka,
Atsumasa Yoshida,
Takuto Ishimura,
Nobuyuki Kawai,
Takashi Shimokawabe,
Kenzo Kinugasa
, et al. (17 additional authors not shown)
Abstract:
We observed an X-ray afterglow of GRB 060904A with the Swift and Suzaku satellites. We found rapid spectral softening during both the prompt tail phase and the decline phase of an X-ray flare in the BAT and XRT data. The observed spectra were fit by power-law photon indices which rapidly changed from $Γ= 1.51^{+0.04}_{-0.03}$ to $Γ= 5.30^{+0.69}_{-0.59}$ within a few hundred seconds in the promp…
▽ More
We observed an X-ray afterglow of GRB 060904A with the Swift and Suzaku satellites. We found rapid spectral softening during both the prompt tail phase and the decline phase of an X-ray flare in the BAT and XRT data. The observed spectra were fit by power-law photon indices which rapidly changed from $Γ= 1.51^{+0.04}_{-0.03}$ to $Γ= 5.30^{+0.69}_{-0.59}$ within a few hundred seconds in the prompt tail. This is one of the steepest X-ray spectra ever observed, making it quite difficult to explain by simple electron acceleration and synchrotron radiation. Then, we applied an alternative spectral fitting using a broken power-law with exponential cutoff (BPEC) model. It is valid to consider the situation that the cutoff energy is equivalent to the synchrotron frequency of the maximum energy electrons in their energy distribution. Since the spectral cutoff appears in the soft X-ray band, we conclude the electron acceleration has been inefficient in the internal shocks of GRB 060904A. These cutoff spectra suddenly disappeared at the transition time from the prompt tail phase to the shallow decay one. After that, typical afterglow spectra with the photon indices of 2.0 are continuously and preciously monitored by both XRT and Suzaku/XIS up to 1 day since the burst trigger time. We could successfully trace the temporal history of two characteristic break energies (peak energy and cutoff energy) and they show the time dependence of $\propto t^{-3} \sim t^{-4}$ while the following afterglow spectra are quite stable. This fact indicates that the emitting material of prompt tail is due to completely different dynamics from the shallow decay component. Therefore we conclude the emission sites of two distinct phenomena obviously differ from each other.
△ Less
Submitted 29 August, 2007;
originally announced August 2007.
-
Swift and Suzaku Observations of the X-Ray Afterglow from the GRB 060105
Authors:
M. S. Tashiro,
K. Abe,
L. Angelini,
S. Barthelmy,
N. Gehrels,
N. Ishikawa,
L. J. Kaluzienski,
N. Kawai,
R. L. Kelley,
K. Kinugasa,
H. Kodaira,
T. Kohmura,
K. Kubota,
Y. Maeda,
S. Maeno,
H. Murakami,
T. Murakami,
Y. E. Nakagawa,
K. Nakazawa,
J. Nousek,
S. Okuno,
K. Onda,
J. N. Reeves,
G. Ricker,
G. Sato
, et al. (12 additional authors not shown)
Abstract:
Results are presented of early X-ray afterglow observations of GRB 060105 by Swift and Suzaku. The bright, long gamma-ray burst GRB 060105 triggered the Swift Burst Alert Telescope (BAT) at 06:49:28 on 5 January 2006. The Suzaku team commenced a pre-planned target of opportunity observation at 19 ks (5.3 hr) after the Swift trigger. Following the prompt emission and successive very steep decay,…
▽ More
Results are presented of early X-ray afterglow observations of GRB 060105 by Swift and Suzaku. The bright, long gamma-ray burst GRB 060105 triggered the Swift Burst Alert Telescope (BAT) at 06:49:28 on 5 January 2006. The Suzaku team commenced a pre-planned target of opportunity observation at 19 ks (5.3 hr) after the Swift trigger. Following the prompt emission and successive very steep decay, a shallow decay was observed from T_0+187 s to T_0+1287 s. After an observation gap during T_0 +(1.5-3) ks, an extremely early steep decay was observed in T_0+(4-30) ks. The lightcurve flattened again at T_0+30 ks, and another steep decay followed from T_0+50 ks to the end of observations. Both steep decays exhibited decay indices of 2.3 - 2.4. This very early break, if it is a jet break, is the earliest case among X-ray afterglow observations, suggesting a very narrow jet whose opening angle is well below 1 degree. The unique Suzaku/XIS data allow us to set very tight upper limits on line emission or absorption in this GRB. For the reported pseudo-redshift of z=4.0+/-1.3 the upper limit on the iron line equivalent width is 50 eV.
△ Less
Submitted 15 September, 2006;
originally announced September 2006.
-
AdS3 Gravitational Instantons from Conformal Field Theory
Authors:
S. Mano
Abstract:
A conformal field theory on the boundary of three-dimensional asymptotic anti-de Sitter spaces which appear as near horizon geometry of D-brane bound states is discussed. It is shown that partition functions of gravitational instantons appear as high and low temperature limits of the partition function of the conformal field theory. The result reproduces phase transition between the anti-de Sitt…
▽ More
A conformal field theory on the boundary of three-dimensional asymptotic anti-de Sitter spaces which appear as near horizon geometry of D-brane bound states is discussed. It is shown that partition functions of gravitational instantons appear as high and low temperature limits of the partition function of the conformal field theory. The result reproduces phase transition between the anti-de Sitter space and the BTZ black hole in the bulk gravity.
△ Less
Submitted 1 July, 1999; v1 submitted 8 March, 1999;
originally announced March 1999.
-
Black Hole Entropy from BPS Closed String
Authors:
Shuhei Mano,
Yuhsuke Yoshida
Abstract:
By using BPS closed string, the entropy is calculated of the extremal five dimensional black hole consisting of Dirichlet onebranes, Dirichlet fivebranes and Kaluza-Klein momentum in the flat background approximation. In our formulations we consider two kinds of BPS closed strings with or without a winding number. In the former case heavy excitation modes of closed strings are used to derive the…
▽ More
By using BPS closed string, the entropy is calculated of the extremal five dimensional black hole consisting of Dirichlet onebranes, Dirichlet fivebranes and Kaluza-Klein momentum in the flat background approximation. In our formulations we consider two kinds of BPS closed strings with or without a winding number. In the former case heavy excitation modes of closed strings are used to derive the entropy. In the latter case we have no oscillator modes and consider collective motion of such massless closed strings. The entropy is given by the number of the ways how we divide the Kaluza-Klein momentum among the massless closed strings. In both cases the black hole entropy is the same as the Bekenstein-Hawking entropy. We argue that the collective modes of closed strings without winding is equivalent to a single closed string with winding. We propose that the two closed string pictures are connected with the open string pictures by the modular transformation.
△ Less
Submitted 11 June, 1999; v1 submitted 15 July, 1998;
originally announced July 1998.
-
Post-Newtonian Expansion of Gravitational Waves from a Particle in Circular Orbits around a Rotating Black Hole :Effects of Black Hole Absorption
Authors:
Hideyuki Tagoshi,
Shuhei Mano,
Eiichi Takasugi
Abstract:
When a particle moves around a Kerr black hole, it radiates gravitational waves.Some of these waves are absorbed by the black hole. We calculate such absorption of gravitational waves induced by a particle of mass mu in a circular orbit on an equatorial plane around a Kerr black hole of mass M. We assume that the velocity of the particle v is much smaller than the speed of light c and calculate…
▽ More
When a particle moves around a Kerr black hole, it radiates gravitational waves.Some of these waves are absorbed by the black hole. We calculate such absorption of gravitational waves induced by a particle of mass mu in a circular orbit on an equatorial plane around a Kerr black hole of mass M. We assume that the velocity of the particle v is much smaller than the speed of light c and calculate the energy absorption rate analytically. We adopt an analytic technique for the Teukolsky equation developed by Mano, Suzuki and Takasugi. We obtain the energy absorption rate to O((v/c)^8) compared to the lowest order. We find that the black hole absorption occurs at O((v/c)^5) beyond the Newtonian-quadrapole luminosity at infinity in the case when the black hole is rotating, which is O((v/c)^3) lower than the non-rotating case. Using the energy absorption rate, we investigate its effects on the orbital evolution of coalescing compact binaries.
△ Less
Submitted 25 November, 1997;
originally announced November 1997.
-
Analytic Solutions of the Teukolsky Equation and their Properties
Authors:
Shuhei Mano,
Eiichi Takasugi
Abstract:
The analytical solutions reported in our previous paper are given as series of hypergeometric or Coulomb wave functions. By using them, we can get the Teukolsky functions analytically in a desired accuracy. For the computation, the deep understanding of their properties is necessary. We summarize the main result: The relative normalization between the solutions with a spin weight s and -s is giv…
▽ More
The analytical solutions reported in our previous paper are given as series of hypergeometric or Coulomb wave functions. By using them, we can get the Teukolsky functions analytically in a desired accuracy. For the computation, the deep understanding of their properties is necessary. We summarize the main result: The relative normalization between the solutions with a spin weight s and -s is given analytically by using the Teukolsky-Starobinsky (T-S) identities. By examining the asymptotic behaviors of our solution and combined with the T-S identities and the Wronskian, we found nontrivial identities between the sums of coefficients of the series. These identities will serve to make various expression in simpler forms and also become a powerful tool to test the accuracy of the computation. As an application, we investigated the absorption rate and the evaporation rate of black hole and obtain interesting analytic results.
△ Less
Submitted 6 November, 1996;
originally announced November 1996.
-
Analytic Solutions of the Regge-Wheeler Equation and the Post-Minkowskian Expansion
Authors:
Shuhei Mano,
Hisao Suzuki,
Eiichi Takasugi
Abstract:
Analytic solutions of the Regge-Wheeler equation are presented in the form of series of hypergeometric functions and Coulomb wave functions which have different regions of convergence. Relations between these solutions are established. The series solutions are given as the Post-Minkowskian expansion with respect to a parameter $ε\equiv 2Mω$, $M$ being the mass of black hole. This expansion corre…
▽ More
Analytic solutions of the Regge-Wheeler equation are presented in the form of series of hypergeometric functions and Coulomb wave functions which have different regions of convergence. Relations between these solutions are established. The series solutions are given as the Post-Minkowskian expansion with respect to a parameter $ε\equiv 2Mω$, $M$ being the mass of black hole. This expansion corresponds to the post-Newtonian expansion when they are applied to the gravitational radiation from a particle in circular orbit around a black hole. These solutions can also be useful for numerical computations.
△ Less
Submitted 27 May, 1996;
originally announced May 1996.
-
Analytic Solutions of the Teukolsky Equation and their Low Frequency Expansions
Authors:
Shuhei Mano,
Hisao Suzuki,
Eiichi Takasugi
Abstract:
Analytic solutions of the Teukolsky equation in Kerr geometries are presented in the form of series of hypergeometric functions and Coulomb wave functions. Relations between these solutions are established. The solutions provide a very powerful method not only for examining the general properties of solutions and physical quantities when they are applied to, but also for numerical computations.…
▽ More
Analytic solutions of the Teukolsky equation in Kerr geometries are presented in the form of series of hypergeometric functions and Coulomb wave functions. Relations between these solutions are established. The solutions provide a very powerful method not only for examining the general properties of solutions and physical quantities when they are applied to, but also for numerical computations. The solutions are given in the expansion of a small parameter $ε\equiv 2Mω$, $M$ being the mass of black hole, which corresponds to Post-Minkowski expansion by $G$ and to post-Newtonian expansion when they are applied to the gravitational radiation from a particle in circular orbit around a black hole. It is expected that these solutions will become a powerful weapon to construct the theoretical template towards LIGO and VIRGO projects.
△ Less
Submitted 19 March, 1996; v1 submitted 13 March, 1996;
originally announced March 1996.