Model selection and parameter estimation of dynamical systems using a novel variant of approximate Bayesian computation

Model selection is a challenging problem that is of importance in many branches of the sciences and engineering, particularly in structural dynamics. By definition, it is intended to select the most likely model among a set of competing models that best matches the dynamic behaviour of a real structure and better predicts the measured data. The Bayesian approach which is based essentially on the evaluation of a likelihood function is one of the most popular approach to deal with model selection and parameter estimation issues. However, in some circumstances, the likelihood function is either intractable or not available even in a closed form. To overcome this issue, the likelihood-free or approximate Bayesian computation (ABC) algorithm has been introduced in the literature, which relaxes the need for an explicit likelihood function to measure the level of agreement between model predictions and measurements. However, ABC algorithms suffer from a low acceptance rate of samples which is actually a common problem with the traditional Bayesian methods. To overcome this shortcoming and alleviate the computational burden, a new variant of the ABC algorithm based on an ellipsoidal Nested Sampling (NS) technique is introduced in this paper; it has been called ABC-NS. Through this paper, it will be shown how the new algorithm is a promising alternative to deal with parameter estimation and model selection issues. It promises drastic speedups and provides a good approximation of the posterior distributions. To demonstrate its robust computational efficiency, four illustrative examples are given. Firstly, the efficiency of the algorithm is demonstrated to deal with parameter estimation. Secondly, two examples based on simulated and real data are given to demonstrate the efficiency of the algorithm to deal with model selection in structural dynamics

Identification of piecewise-linear mechanical oscillators via Bayesian model selection and parameter estimation

The problem of identifying single degree-of-freedom (SDOF) nonlinear mechanical oscillators with piecewise-linear (PWL) restoring forces is considered. PWL nonlinear systems are a class of models that specify or approximate nonlinear systems via a set of locally-linear maps, each defined over different operating regions. They are useful in modelling hybrid phenomena common in practical situations, such as, systems with different modes of operation, or systems whose dynamics change because of physical limits or thresholds. However, identifying PWL models can be a challenging task when the number of operating regions and their partitions are unknown. This paper formulates the identification of oscillators with PWL restoring forces as a task of concurrent model selection and parameter estimation, where the selection of the number of linear regions is treated as a model selection task and identifying the associated system parameters as a task of parameter estimation. In this study, PWL maps in restoring forces with up to four regions are considered, and the task of model selection and parameter estimation task is addressed in a Bayesian framework. A likelihood-free Approximate Bayesian Computation (ABC) scheme is followed, which is easy to implement and provides a simplified way of doing model selection. The proposed approach has been demonstrated using two numerical examples and an experimental study, where ABC has been used to select models and identify parameters from among four SDOF PWL systems with different number of PWL regions. The results demonstrate the flexibility of using the proposed Bayesian approach for identifying the correct model and parameters of PWL systems, in addition to furnishing uncertainty estimates of the identified parameters

Exotic ρ±ρ0\rho^\pm\rho^0 state photoproduction

It is shown that the list of unusual mesons planned for a careful study in photoproduction can be extended by the exotic states $X^\pm(1600)$ with $I^G(J ^{PC})=2^+(2^{++})$ which should be looked for in the $\rho^\pm\rho^0$ decay channels in the reactions $\gamma N\to\rho^\pm\rho^0N$ and $\gamma N\to\rho^\pm \rho^0\Delta$. The full classification of the $\rho^\pm\rho^0$ states by their quantum numbers is presented. A simple model for the spin structure of the $\gamma p\to f_2(1270)p$, $\gamma p\to a^0_2(1320)p$, and $\gamma N\to X^\pm (N, \Delta)$ reaction amplitudes is formulated and the tentative estimates of the corresponding cross sections at the incident photon energy $E_\gamma\approx 6$ GeV are obtained: $\sigma(\gamma p\to f_2(1270)p)\approx0.12$ $\mu$b, $\sigma(\gamma p\to a^0_2(1320)p)\approx0.25$ $\mu$b, $\sigma(\gamma N\to X^\pm N\to\rho^\pm\rho^0N)\approx0.018$ $\mu$b, and $\sigma(\gamma p\to X^-\Delta^{++ }\to\rho^-\rho^0\Delta^{++})\approx0.031$ $\mu$b. The problem of the $X^\pm$ signal extraction from the natural background due to the other $\pi^\pm\pi^0 \pi^+\pi^-$ production channels is discussed. In particular the estimates are presented for the $\gamma p\to h_1(1170)\pi^+n$, $\gamma p\to\rho'^{+}n\to \pi^+\pi^0\pi^+\pi^-n$, and $\gamma p\to\omega\rho^0p$ reaction cross sections. Our main conclusion is that the search for the exotic $X^\pm(2^+(2^{++}))$ states is quite feasible at JEFLAB facility. The expected yield of the $\gamma N\to X^\pm N\to\rho^\pm\rho^0N$ events in a 30-day run at the 100% detection efficiency approximates $2.8\times10^6$ events.Comment: 19 pages, revtex, 1 figure in postscipt, some comments and references added, a few minor typos corrected, to be published in Phys. Rev.

All-sky search for long-duration gravitational wave transients with initial LIGO

We present the results of a search for long-duration gravitational wave transients in two sets of data collected by the LIGO Hanford and LIGO Livingston detectors between November 5, 2005 and September 30, 2007, and July 7, 2009 and October 20, 2010, with a total observational time of 283.0 days and 132.9 days, respectively. The search targets gravitational wave transients of duration 10-500 s in a frequency band of 40-1000 Hz, with minimal assumptions about the signal waveform, polarization, source direction, or time of occurrence. All candidate triggers were consistent with the expected background; as a result we set 90% confidence upper limits on the rate of long-duration gravitational wave transients for different types of gravitational wave signals. For signals from black hole accretion disk instabilities, we set upper limits on the source rate density between 3.4×10-5 and 9.4×10-4 Mpc-3 yr-1 at 90% confidence. These are the first results from an all-sky search for unmodeled long-duration transient gravitational waves. © 2016 American Physical Society

All-sky search for long-duration gravitational wave transients with initial LIGO

We present the results of a search for long-duration gravitational wave transients in two sets of data collected by the LIGO Hanford and LIGO Livingston detectors between November 5, 2005 and September 30, 2007, and July 7, 2009 and October 20, 2010, with a total observational time of 283.0 days and 132.9 days, respectively. The search targets gravitational wave transients of duration 10-500 s in a frequency band of 40-1000 Hz, with minimal assumptions about the signal waveform, polarization, source direction, or time of occurrence. All candidate triggers were consistent with the expected background; as a result we set 90% confidence upper limits on the rate of long-duration gravitational wave transients for different types of gravitational wave signals. For signals from black hole accretion disk instabilities, we set upper limits on the source rate density between 3.4×10-5 and 9.4×10-4 Mpc-3 yr-1 at 90% confidence. These are the first results from an all-sky search for unmodeled long-duration transient gravitational waves. © 2016 American Physical Society

Search for Tensor, Vector, and Scalar Polarizations in the Stochastic Gravitational-Wave Background

The detection of gravitational waves with Advanced LIGO and Advanced Virgo has enabled novel tests of general relativity, including direct study of the polarization of gravitational waves. While general relativity allows for only two tensor gravitational-wave polarizations, general metric theories can additionally predict two vector and two scalar polarizations. The polarization of gravitational waves is encoded in the spectral shape of the stochastic gravitational-wave background, formed by the superposition of cosmological and individually unresolved astrophysical sources. Using data recorded by Advanced LIGO during its first observing run, we search for a stochastic background of generically polarized gravitational waves. We find no evidence for a background of any polarization, and place the first direct bounds on the contributions of vector and scalar polarizations to the stochastic background. Under log-uniform priors for the energy in each polarization, we limit the energy densities of tensor, vector, and scalar modes at 95% credibility to Ω0T<5.58×10-8, Ω0V<6.35×10-8, and Ω0S<1.08×10-7 at a reference frequency f0=25 Hz. © 2018 American Physical Society

On the progenitor of binary neutron star merger GW170817

On 2017 August 17 the merger of two compact objects with masses consistent with two neutron stars was discovered through gravitational-wave (GW170817), gamma-ray (GRB 170817A), and optical (SSS17a/AT 2017gfo) observations. The optical source was associated with the early-type galaxy NGC 4993 at a distance of just ∼40 Mpc, consistent with the gravitational-wave measurement, and the merger was localized to be at a projected distance of ∼2 kpc away from the galaxy's center. We use this minimal set of facts and the mass posteriors of the two neutron stars to derive the first constraints on the progenitor of GW170817 at the time of the second supernova (SN). We generate simulated progenitor populations and follow the three-dimensional kinematic evolution from binary neutron star (BNS) birth to the merger time, accounting for pre-SN galactic motion, for considerably different input distributions of the progenitor mass, pre-SN semimajor axis, and SN-kick velocity. Though not considerably tight, we find these constraints to be comparable to those for Galactic BNS progenitors. The derived constraints are very strongly influenced by the requirement of keeping the binary bound after the second SN and having the merger occur relatively close to the center of the galaxy. These constraints are insensitive to the galaxy's star formation history, provided the stellar populations are older than 1 Gyr