Bayesian Bounds on Parameter Estimation Accuracy for Compact Coalescing Binary Gravitational Wave Signals

A global network of laser interferometric gravitational wave detectors is projected to be in operation by around the turn of the century. Here, the noisy output of a single instrument is examined. A gravitational wave is assumed to have been detected in the data and we deal with the subsequent problem of parameter estimation. Specifically, we investigate theoretical lower bounds on the minimum mean-square errors associated with measuring the parameters of the inspiral waveform generated by an orbiting system of neutron stars/black holes. Three theoretical lower bounds on parameter estimation accuracy are considered: the Cramer-Rao bound (CRB); the Weiss-Weinstein bound (WWB); and the Ziv-Zakai bound (ZZB). We obtain the WWB and ZZB for the Newtonian-form of the coalescing binary waveform, and compare them with published CRB and numerical Monte-Carlo results. At large SNR, we find that the theoretical bounds are all identical and are attained by the Monte-Carlo results. As SNR gradually drops below 10, the WWB and ZZB are both found to provide increasingly tighter lower bounds than the CRB. However, at these levels of moderate SNR, there is a significant departure between all the bounds and the numerical Monte-Carlo results.Comment: 17 pages (LaTeX), 4 figures. Submitted to Physical Review

Searching for t->c(u)h with dipole moments

A discovery of flavour-changing Higgs-boson decays would constitute an undeniable signal of new physics. We derive model-independent constraints on the tch and tuh couplings that arise from the bounds on hadronic electric dipole moments. Comparisons of the present and future sensitivities with both the direct LHC constraints and the indirect limits from D-meson physics are also presented.Comment: 18 pages, 2 figures; v2: references added and typos corrected; coincides with version to be published in JHE

Numerical studies of the fluid and optical fields associated with complex cavity flows

Numerical solutions for the flowfield about several cavity configurations have been computed using the Reynolds averaged Navier-Stokes equations. Comparisons between numerical and experimental results are made in two dimensions for free shear layers and a rectangular cavity, and in three dimensions for the transonic aero-window problem of the Stratospheric Observatory for Infrared Astronomy (SOFIA). Results show that dominant acoustic frequencies and magnitudes of the self excited resonant cavity flows compare well with the experiment. In addition, solution sensitivity to artificial dissipation and grid resolution levels are determined. Optical path distortion due to the flow field is modelled geometrically and is found to match the experiment. The fluid field was computed using a diagonalized scheme within an overset mesh framework. An existing code, OVERFLOW, was utilized with the additions of characteristic boundary condition and output routines required for reduction of the unsteady data. The newly developed code is directly applicable to a generalized three dimensional structured grid zone. Details are provided in a paper included in Appendix A

Nonintrusive approximation of parametrized limits of matrix power algorithms -- application to matrix inverses and log-determinants

We consider in this work quantities that can be obtained as limits of powers of parametrized matrices, for instance the inverse matrix or the logarithm of the determinant. Under the assumption of affine dependence in the parameters, we use the Empirical Interpolation Method (EIM) to derive an approximation for powers of these matrices, from which we derive a nonintrusive approximation for the aforementioned limits. We derive upper bounds of the error made by the obtained formula. Finally, numerical comparisons with classical intrusive and nonintrusive approximation techniques are provided: in the considered test-cases, our algorithm performs well compared to the nonintrusive ones

Resummation of heavy jet mass and comparison to LEP data

The heavy jet mass distribution in e+e- collisions is computed to next-to-next-to-next-to leading logarithmic (NNNLL) and next-to-next-to leading fixed order accuracy (NNLO). The singular terms predicted from the resummed distribution are confirmed by the fixed order distributions allowing a precise extraction of the unknown soft function coefficients. A number of quantitative and qualitative comparisons of heavy jet mass and the related thrust distribution are made. From fitting to ALEPH data, a value of alpha_s is extracted, alpha_s(m_Z)=0.1220 +/- 0.0031, which is larger than, but not in conflict with, the corresponding value for thrust. A weighted average of the two produces alpha_s(m_Z) = 0.1193 +/- 0.0027, consistent with the world average. A study of the non-perturbative corrections shows that the flat direction observed for thrust between alpha_s and a simple non-perturbative shape parameter is not lifted in combining with heavy jet mass. The Monte Carlo treatment of hadronization gives qualitatively different results for thrust and heavy jet mass, and we conclude that it cannot be trusted to add power corrections to the event shape distributions at this accuracy. Whether a more sophisticated effective field theory approach to power corrections can reconcile the thrust and heavy jet mass distributions remains an open question.Comment: 33 pages, 14 figures. v2 added effect of lower numerical cutoff with improved extraction of the soft function constants; power correction discussion clarified. v3 small typos correcte

Internal states of model isotropic granular packings. III. Elastic properties

In this third and final paper of a series, elastic properties of numerically simulated isotropic packings of spherical beads assembled by different procedures and subjected to a varying confining pressure P are investigated. In addition P, which determines the stiffness of contacts by Hertz's law, elastic moduli are chiefly sensitive to the coordination number, the possible values of which are not necessarily correlated with the density. Comparisons of numerical and experimental results for glass beads in the 10kPa-10MPa range reveal similar differences between dry samples compacted by vibrations and lubricated packings. The greater stiffness of the latter, in spite of their lower density, can hence be attributed to a larger coordination number. Voigt and Reuss bounds bracket bulk modulus B accurately, but simple estimation schemes fail for shear modulus G, especially in poorly coordinated configurations under low P. Tenuous, fragile networks respond differently to changes in load direction, as compared to load intensity. The shear modulus, in poorly coordinated packings, tends to vary proportionally to the degree of force indeterminacy per unit volume. The elastic range extends to small strain intervals, in agreement with experimental observations. The origins of nonelastic response are discussed. We conclude that elastic moduli provide access to mechanically important information about coordination numbers, which escape direct measurement techniques, and indicate further perspectives.Comment: Published in Physical Review E 25 page

Predicting response bounds for friction-damped gas turbine blades with uncertain friction coupling

Friction dampers are often used to reduce high amplitude vibration within gas turbines: they are a robust solution that are able to withstand extreme operating environments. Although the turbine blades are manufactured to tight tolerances, there can be significant variability in the overall response of the assembly. Uncertainties associated with the frictional contact properties are a major factor contributing to this variability. This paper applies a recently developed method for predicting response bounds to friction-damped gas turbines when the characteristics of the friction dampers are unknown, including uncertainty regarding the underlying functional form of the friction law. The approach taken is to represent the frictional contact using a describing function, and formulate an optimisation problem to seek upper and lower bounds on a chosen response metric, such as displacement amplitude. Constraints are chosen that describe known properties of the frictional nonlinearity, without needing to specify a particular constitutive law. The method was validated by comparison with numerical and experimental results from an idealised test system. The experimental test rig consisted of an array of eight beams coupled by pin-contact friction dampers. A modal description of this test rig formed the basis of a numerical model, which uses the Harmonic Balance Method (HBM) for nonlinear simulations. A set of Monte Carlo tests was carried out numerically and experimentally for both a two-beam sub-assembly as well as for the full eight-beam assembly. Comparisons with numerical results showed excellent agreement providing confident verification of the implementation, and comparisons with experimental results revealed that the bounds became less conservative as the system complexity increased. Overall the results are promising: upper and lower response bounds for an array of friction-damped systems can be computed at similar cost to a single HBM simulation, giving reliable bounds that are valid for both parametric and model uncertainties associated with the friction couplings.Mitsubishi Heavy Industrie