Kinetic energy choice in Hamiltonian/hybrid Monte Carlo

We consider how different choices of kinetic energy in Hamiltonian Monte Carlo affect algorithm performance. To this end, we introduce two quantities which can be easily evaluated, the composite gradient and the implicit noise. Results are established on integrator stability and geometric convergence, and we show that choices of kinetic energy that result in heavy-tailed momentum distributions can exhibit an undesirable negligible moves property, which we define. A general efficiency-robustness trade off is outlined, and implementations which rely on approximate gradients are also discussed. Two numerical studies illustrate our theoretical findings, showing that the standard choice which results in a Gaussian momentum distribution is not always optimal in terms of either robustness or efficiency.Comment: 15 pages (+7 page supplement, included here as an appendix), 2 figures (+1 in supplement

Following the relaxation dynamics of photoexcited aniline in the 273-266 nm region using time-resolved photoelectron imaging

Scanning tunnelling microscopy (STM) and Auger electron spectroscopy (AES) have been used to investigate the growth morphology of ultra-thin Pb films on the Ni3Al(111) face at room temperature. A previous study [K. Miśków and A. Krupski Appl Surf Sci 273, 2013, 554] using low-energy electron diffraction (LEED) and real time Auger intensity recording has demonstrated that an initial two-dimensional growth of the first Pb monolayer thick ‘wetting layer’ takes place. With further deposition and for T = 300 K, flat three atomic-layer-high islands are grown. Above 350 K, the Stranski–Krastanov growth mode was observed. In the current study, the analysis of STM measurements indicate and confirm that for coverage θ = 1.0 ML two-dimensional growth of the first Pb monolayer took place. Above θ > 1.0 ML, a three-dimensional growth of the Pb islands was observed with a strongly preferred atomic-scale ‘magic height (N),’ hexagonal shape and flat-tops. At coverage θ = 3.5 ML, only islands containing N = 3, 5, 7 and 11 atomic layers of Pb are observed. At the higher coverage θ = 5.5 ML, three types of regular hexagonal islands with side lengths of 25, 30 and 45 nm are observed. Furthermore, three different island adsorption configurations rotated by 10° ± 1° and 30° ± 6° with respect to each other were observed. After an annealing at T = 400 K of 5.5 ML of lead deposited at RT on the Ni3Al(111) the morphology of the surface changes. Post-anneal, islands of Pb are observed above the ‘wetting layer’ with an estimated average size and diameter of 768 ± 291 nm2 and 38.17 ± 6.56 nm and constant uniform height of two atomic layers (N = 2)

On-Line Fatigue Crack Growth Monitoring in Externally Pressurised Vessels Using the Alternating Current Potential Drop (ACPD) Technique

At ARE (Dunfermline) fatigue tests are carried out on internally stiffened, welded steel cylindrical pressure vessels. These vessels are constructed from high yield strength quenched and tempered steels with full penetration butt and T-butt welds

Post-outburst X-ray flux and timing evolution of Swift J1822.3-1606

Swift J1822.3-1606 was discovered on 2011 July 14 by the Swift Burst Alert Telescope following the detection of several bursts. The source was found to have a period of 8.4377 s and was identified as a magnetar. Here we present a phase-connected timing analysis and the evolution of the flux and spectral properties using RXTE, Swift, and Chandra observations. We measure a spin frequency of 0.1185154343(8) s$^{-1}$ and a frequency derivative of $-4.3\pm0.3\times10^{-15}$ at MJD 55761.0, in a timing analysis that include significant non-zero second and third frequency derivatives that we attribute to timing noise. This corresponds to an estimated spin-down inferred dipole magnetic field of $B\sim5\times10^{13}$ G, consistent with previous estimates though still possibly affected by unmodelled noise. We find that the post-outburst 1--10 keV flux evolution can be characterized by a double-exponential decay with decay timescales of $15.5\pm0.5$ and $177\pm14$ days. We also fit the light curve with a crustal cooling model which suggests that the cooling results from heat injection into the outer crust. We find that the hardness-flux correlation observed in magnetar outbursts also characterizes the outburst of Swift J1822.3-1606. We compare the properties of Swift J1822.3-1606 with those of other magnetars and their outbursts.Comment: 12 pages, 8 figures, accepted for publication in Ap

A Survey of 56 Mid-latitude EGRET Error Boxes for Radio Pulsars

We have conducted a radio pulsar survey of 56 unidentified gamma-ray sources from the 3rd EGRET catalog which are at intermediate Galactic latitudes (5 deg. < |b| < 73 deg.). For each source, four interleaved 35-minute pointings were made with the 13-beam, 1400-MHz multibeam receiver on the Parkes 64-m radio telescope. This covered the 95% error box of each source at a limiting sensitivity of about 0.2 mJy to pulsed radio emission for periods P > 10 ms and dispersion measures < 50 pc cm-3. Roughly half of the unidentified gamma-ray sources at |b| > 5 deg. with no proposed active galactic nucleus counterpart were covered in this survey. We detected nine isolated pulsars and four recycled binary pulsars, with three from each class being new. Timing observations suggest that only one of the pulsars has a spin-down luminosity which is even marginally consistent with the inferred luminosity of its coincident EGRET source. Our results suggest that population models, which include the Gould belt as a component, overestimate the number of isolated pulsars among the mid-latitude Galactic gamma-ray sources and that it is unlikely that Gould belt pulsars make up the majority of these sources. However, the possibility of steep pulsar radio spectra and the confusion of terrestrial radio interference with long-period pulsars (P > 200 ms) having very low dispersion measures (< 10 pc cm-3, expected for sources at a distance of less than about 1 kpc) prevent us from strongly ruling out this hypothesis. Our results also do not support the hypothesis that millisecond pulsars make up the majority of these sources. Non-pulsar source classes should therefore be further investigated as possible counterparts to the unidentified EGRET sources at intermediate Galactic latitudes.Comment: 24 pages, including 4 figures and 3 tables. Accepted for publication in Ap

Magnetar-like Emission from the Young Pulsar in Kes 75

We report detection of magnetar-like X-ray bursts from the young pulsar PSR J1846-0258, at the center of the supernova remnant Kes 75. This pulsar, long thought to be rotation-powered, has an inferred surface dipolar magnetic field of 4.9x10^13 G, higher than those of the vast majority of rotation-powered pulsars, but lower than those of the ~12 previously identified magnetars. The bursts were accompanied by a sudden flux increase and an unprecedented change in timing behavior. These phenomena lower the magnetic and rotational thresholds associated with magnetar-like behavior, and suggest that in neutron stars there exists a continuum of magnetic activity that increases with inferred magnetic field strength.Comment: 17 pages, 2 figures, accepted for publication in Science. Note: The content of this paper is embargoed until February 21, 200

An automated method for mapping geomorphological expressions of former subglacial meltwater pathways (hummock corridors) from high resolution digital elevation data

Elongated tracts of hummocks or ‘hummock corridors’, exposed on palaeo-ice sheet beds, are believed to represent former subglacial meltwater pathways. Here, we present a method, coded in MATLAB, for automatically detecting and mapping hummock corridors from high-resolution digital elevation models (DEMs). Initially the DEM is filtered to remove bed roughness outside the size range of hummocks. A Fast Fourier Transform is then performed to determine the dominant orientation of hummock corridors and remove misaligned features. Finally, image segmentation is used to isolate and extract the hummock corridors as a binary mask. We tested this automated approach visually and statistically against detailed manual mapping in three areas of Canada and northern Scandinavia. Results show that while the automated method does not perfectly reproduce the manual mapping, it successfully captures the general configuration, morphometry (length, width) and location of hummock corridors, despite variation in expression across and between sites. This technique is ideally suited to take advantage of newly available high-resolution digital elevation data (e.g. the ArcticDEM), whose enormous volume makes large-scale manual mapping prohibitively time consuming. Its application will enable efficient and comprehensive mapping of the spatial distribution of hummock corridors across palaeo-beds that is necessary for deriving insights into their formation and the organisation of subglacial meltwater flow beneath ice sheets