A Bayesian method for pulsar template generation

Extracting Times of Arrival from pulsar radio signals depends on the knowledge of the pulsars pulse profile and how this template is generated. We examine pulsar template generation with Bayesian methods. We will contrast the classical generation mechanism of averaging intensity profiles with a new approach based on Bayesian inference. We introduce the Bayesian measurement model imposed and derive the algorithm to reconstruct a "statistical template" out of noisy data. The properties of these "statistical templates" are analysed with simulated and real measurement data from PSR B1133+16. We explain how to put this new form of template to use in analysing secondary parameters of interest and give various examples: We implement a nonlinear filter for determining ToAs of pulsars. Applying this method to data from PSR J1713+0747 we derive ToAs self consistently, meaning all epochs were timed and we used the same epochs for template generation. While the average template contains fluctuations and noise as unavoidable artifacts, we find that the "statistical template" derived by Bayesian inference quantifies fluctuations and remaining uncertainty. This is why the algorithm suggested turns out to reconstruct templates of statistical significance from ten to fifty single pulses. A moving data window of fifty pulses, taking out one single pulse at the beginning and adding one at the end of the window unravels the characteristics of the methods to be compared. It shows that the change induced in the classical reconstruction is dominated by random fluctuations for the average template, while statistically significant changes drive the dynamics of the proposed method's reconstruction. The analysis of phase shifts with simulated data reveals that the proposed nonlinear algorithm is able to reconstruct correct phase information along with an acceptable estimation of the remaining uncertainty.Comment: 21 pages, 16 figures, submitted to MNRA

Comparison of porcine thorax to gelatine blocks for wound

Published online first in International Journal of Legal Medicine. The support of EPSRC and The Home Office are recognised. Open Access, this article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http:/ /creativecommons.org/licenses/by/4.0/)Tissue simulants are typically used in ballistic testing as substitutes for biological tissues. Many simulants have been used, with gelatine amongst the most common. While two concentrations of gelatine (10 and 20 %) have been used extensively, no agreed standard exists for the preparation of either. Comparison of ballistic damage produced in both concentrations is lacking. The damage produced in gelatine is also questioned, with regards to what it would mean for specific areas of living tissue. The aim of the work discussed in this paper was to consider how damage caused by selected pistol and rifle ammunition varied in different simulants. Damage to gelatine blocks 10 and 20 % in concentration were tested with 9 mm Luger (9 × 19 full metal jacket; FMJ) rounds, while damage produced by .223 Remington (5.56 × 45 Federal Premium® Tactical® Bonded®) rounds to porcine thorax sections (skin, underlying tissue, ribs, lungs, ribs, underlying tissue, skin; backed by a block of 10 % gelatine) were compared to 10 and 20 % gelatine blocks. Results from the .223 Remington rifle round, which is one that typically expands on impact, revealed depths of penetration in the thorax arrangement were significantly different to 20 % gelatine, but not 10 % gelatine. The level of damage produced in the simulated thoraxes was smaller in scale to that witnessed in both gelatine concentrations,though greater debris was produced in the thoraxes.The support of EPSRC and The Home Office are recognised

Data-driven discovery of coordinates and governing equations

The discovery of governing equations from scientific data has the potential to transform data-rich fields that lack well-characterized quantitative descriptions. Advances in sparse regression are currently enabling the tractable identification of both the structure and parameters of a nonlinear dynamical system from data. The resulting models have the fewest terms necessary to describe the dynamics, balancing model complexity with descriptive ability, and thus promoting interpretability and generalizability. This provides an algorithmic approach to Occam's razor for model discovery. However, this approach fundamentally relies on an effective coordinate system in which the dynamics have a simple representation. In this work, we design a custom autoencoder to discover a coordinate transformation into a reduced space where the dynamics may be sparsely represented. Thus, we simultaneously learn the governing equations and the associated coordinate system. We demonstrate this approach on several example high-dimensional dynamical systems with low-dimensional behavior. The resulting modeling framework combines the strengths of deep neural networks for flexible representation and sparse identification of nonlinear dynamics (SINDy) for parsimonious models. It is the first method of its kind to place the discovery of coordinates and models on an equal footing.Comment: 25 pages, 6 figures; added acknowledgment

Trivalent scandium, yttrium and lanthanide complexes with thia-oxa and selena-oxa macrocycles and crown ether coordination

Complexes of the oxa-thia macrocycles [18]aneO4S2, [15]aneO3S2 and the oxa-selena macrocycle [18]aneO4Se2 (L) of types [MCl2(L)]FeCl4 (M = Sc or Y) were prepared from [ScCl3(thf)3] or [YCl2(THF)5][YCl4(THF)2] and the ligand in anhydrous MeCN, using FeCl3 as a chloride abstractor. The [MI2(L)]I, [LaI3(L)] and [LuI2(L)]I have been prepared from the ligands and the appropriate anhydrous metal triiodide in MeCN. Complexes of type [LaI3(crown)] and [LuI2(crown)]I (crown = 18-crown-6, 15-crown-5) were made for comparison. Use of the metal iodide results in complexes with high solubility compared to the corresponding chlorides, although also with increased sensitivity to moisture. All complexes were characterised by microanalysis, IR, (1)H, (45)Sc and (77)Se NMR spectroscopy as appropriate. X-ray crystal structures are reported for [ScCl2([18]aneO4S2)][FeCl4], [ScI2([18]aneO4S2)]I, [YCl2(18-crown-6)]3[Y2Cl9], [YCl2([18]aneO4S2)][FeCl4], [LaI3(15-crown-5)], [LaI2(18-crown-6)(MeCN)]I, [LuI(18-crown-6)(MeCN)2]I2, [Lu(15-crown-5)(MeCN)2(OH2)]I3, [LaI3([18]aneO4S2)], [LaI([18]aneO4S2)(OH2)]I2, [LaI3([18]aneO4Se2)] and [LuI2([18]aneO4Se2)]I. In each complex all the neutral donor atoms of the macrocycles are coordinated to the metal centre, showing very rare examples of these oxophilic metal centres coordinated to thioether groups, and the first examples of coordinated selenoether donors. In some cases MeCN or adventitious water displaces halide ligands, but not the S/Se donors from La or Lu complexes. A complex of the oxa-tellura macrocycle [18]aneO4Te2, [ScCl2([18]aneO4Te2)][FeCl4] was isolated, but is unstable in MeCN solution, depositing elemental Te. YCl3 and 18-crown-6 produced [YCl2(18-crown-6)]3[Y2Cl9], the asymmetric unit of which contains two cations with a trans-YCl2 arrangement and a third with a cis-YCl2 group

Deep Learning How to Fit an Intravoxel Incoherent Motion Model to Diffusion-Weighted MRI

Purpose: This prospective clinical study assesses the feasibility of training a deep neural network (DNN) for intravoxel incoherent motion (IVIM) model fitting to diffusion-weighted magnetic resonance imaging (DW-MRI) data and evaluates its performance. Methods: In May 2011, ten male volunteers (age range: 29 to 53 years, mean: 37 years) underwent DW-MRI of the upper abdomen on 1.5T and 3.0T magnetic resonance scanners. Regions of interest in the left and right liver lobe, pancreas, spleen, renal cortex, and renal medulla were delineated independently by two readers. DNNs were trained for IVIM model fitting using these data; results were compared to least-squares and Bayesian approaches to IVIM fitting. Intraclass Correlation Coefficients (ICC) were used to assess consistency of measurements between readers. Intersubject variability was evaluated using Coefficients of Variation (CV). The fitting error was calculated based on simulated data and the average fitting time of each method was recorded. Results: DNNs were trained successfully for IVIM parameter estimation. This approach was associated with high consistency between the two readers (ICCs between 50 and 97%), low intersubject variability of estimated parameter values (CVs between 9.2 and 28.4), and the lowest error when compared with least-squares and Bayesian approaches. Fitting by DNNs was several orders of magnitude quicker than the other methods but the networks may need to be re-trained for different acquisition protocols or imaged anatomical regions. Conclusion: DNNs are recommended for accurate and robust IVIM model fitting to DW-MRI data. Suitable software is available at (1)

Opportunities for weed manipulation using GMHT row crops

The herbicides and cultivation systems available in most non-GM crops allow farmers little flexibility as to when they control weeds. However, glyphosate and glufosinate-ammonium, as used in GM herbicide tolerant crops, offer the opportunity to control large weeds and weed control can be timed according to the agronomic and environmental aims of the user. This paper will use sugar beet as a model crop and report results where different approaches to weed control have been used and discuss their relevance in the wider agricultural and environmental contextNon peer reviewe

The Revamped FISA: Striking a Better Balance Between the Government\u27s Need to Protect Itself and the 4th Amendment

The investigations of the 9/11 terrorist attacks highlighted a series of lapses in intelligence-sharing within the federal government regarding terrorist operations. One area closely examined by Congress,\u27 the judiciary, and many legal and political commentators is the appropriate scope of intelligence collection within the United States concerning foreign threats to the nation\u27s security ( foreign intelligence ). Domestic intelligence collection is a particularly complex sphere of national security as gathering intelligence on American soil requires balancing the privacy rights of individuals guaranteed by the Fourth Amendment against the nation\u27s need to protect itself. The Foreign Intelligence Surveillance Act of 1978 ( FISA ) governs the conduct of electronic surveillance and physical searches carried out for foreign intelligence purposes within the United States. FISA establishes procedures for collecting foreign intelligence information, which are parallel to, and independent of, the conventional law enforcement channels used to secure judicial approval for searches and electronic surveillance. Under FISA, federal investigators submit applications for foreign intelligence surveillance to a secret court that exists for the sole purpose of reviewing government requests to gather information pursuant to the statute

PSR J1453+1902 and the radio luminosities of solitary versus binary millisecond pulsars

We present 3 yr of timing observations for PSR J1453+1902, a 5.79-ms pulsar discovered during a 430-MHz drift-scan survey with the Arecibo telescope. Our observations show that PSR J1453+1902 is solitary and has a proper motion of 8(2) mas/yr. At the nominal distance of 1.2 kpc estimated from the pulsar's dispersion measure, this corresponds to a transverse speed of 46(11) km/s, typical of the millisecond pulsar population. We analyse the current sample of 55 millisecond pulsars in the Galactic disk and revisit the question of whether the luminosities of isolated millisecond pulsars are different from their binary counterparts. We demonstrate that the apparent differences in the luminosity distributions seen in samples selected from 430-MHz surveys can be explained by small-number statistics and observational selection biases. An examination of the sample from 1400-MHz surveys shows no differences in the distributions. The simplest conclusion from the current data is that the spin, kinematic, spatial and luminosity distributions of isolated and binary millisecond pulsars are consistent with a single homogeneous population.Comment: 8 pages, 5 figures and 3 tables, accepted for publication by MNRA