Pseudo-High-Order Symplectic Integrators

Symplectic N-body integrators are widely used to study problems in celestial mechanics. The most popular algorithms are of 2nd and 4th order, requiring 2 and 6 substeps per timestep, respectively. The number of substeps increases rapidly with order in timestep, rendering higher-order methods impractical. However, symplectic integrators are often applied to systems in which perturbations between bodies are a small factor of the force due to a dominant central mass. In this case, it is possible to create optimized symplectic algorithms that require fewer substeps per timestep. This is achieved by only considering error terms of order epsilon, and neglecting those of order epsilon^2, epsilon^3 etc. Here we devise symplectic algorithms with 4 and 6 substeps per step which effectively behave as 4th and 6th-order integrators when epsilon is small. These algorithms are more efficient than the usual 2nd and 4th-order methods when applied to planetary systems.Comment: 14 pages, 5 figures. Accepted for publication in the Astronomical Journa

Predicting the movements of permanently installed electrodes on an active landslide using time-lapse geoelectrical resistivity data only

If electrodes move during geoelectrical resistivity monitoring and their new positions are not incorporated in the inversion, then the resulting tomographic images exhibit artefacts that can obscure genuine time-lapse resistivity changes in the subsurface. The effects of electrode movements on time-lapse resistivity tomography are investigated using a simple analytical model and real data. The correspondence between the model and the data is sufficiently good to be able to predict the effects of electrode movements with reasonable accuracy. For the linear electrode arrays and 2D inversions under consideration, the data are much more sensitive to longitudinal than transverse or vertical movements. Consequently the model can be used to invert the longitudinal offsets of the electrodes from their known baseline positions using only the time-lapse ratios of the apparent resistivity data. The example datasets are taken from a permanently installed electrode array on an active lobe of a landslide. Using two sets with different levels of noise and subsurface resistivity changes, it is found that the electrode positions can be recovered to an accuracy of 4 % of the baseline electrode spacing. This is sufficient to correct the artefacts in the resistivity images, and provides for the possibility of monitoring the movement of the landslide and its internal hydraulic processes simultaneously using electrical resistivity tomography only

A framework for initialising a dynamic clustering algorithm: ART2-A

Algorithms in the Adaptive Resonance Theory (ART) family adapt to structural changes in data as new information presents, making it an exciting candidate for dynamic online clustering of big health data. Its use however has largely been restricted to the signal processing field. In this paper we introduce an refinement of the ART2-A method within an adapted separation and concordance (SeCo) framework which has been shown to identify stable and reproducible solutions from repeated initialisations that also provides evidence for an appropriate number of initial clusters that best calibrates the algorithm with the data presented. The results show stable, reproducible solutions for a mix of real-world heath related datasets and well known benchmark datasets, selecting solutions which better represent the underlying structure of the data than using a single measure of separation. The scalability of the method and it's facility for dynamic online clustering makes it suitable for finding structure in big data

Temperature and Pressure Sensors Based on Spin-Allowed Broadband Luminescence of Doped Orthorhombic Perovskite Structures

Systems and methods that are capable of measuring pressure or temperature based on luminescence are discussed herein. These systems and methods are based on spin-allowed broadband luminescence of sensors with orthorhombic perovskite structures of rare earth aluminates doped with chromium or similar transition metals, such as chromium-doped gadolinium aluminate. Luminescence from these sensors can be measured to determine at least one of temperature or pressure, based on either the intense luminescence of these sensors, even at high temperatures, or low temperature techniques discussed herein

IgG anti-apolipoprotein A-1 antibodies in patients with systemic lupus erythematosus are associated with disease activity and corticosteroid therapy: an observational study.

IgG anti-apolipoprotein A-1 (IgG anti-apoA-1) antibodies are present in patients with systemic lupus erythematosus (SLE) and may link inflammatory disease activity and the increased risk of developing atherosclerosis and cardiovascular disease (CVD) in these patients. We carried out a rigorous analysis of the associations between IgG anti-apoA-1 levels and disease activity, drug therapy, serology, damage, mortality and CVD events in a large British SLE cohort

Temperature Sensing Above 1000 C Using Cr-Doped GdAlO3 Spin-Allowed Broadband Luminescence

Cr-doped GdAlO3 (Cr:GdAlO3) is shown to produce remarkably high-intensity spin-allowed broadband luminescence with sufficiently long decay times to make effective luminescence-decay-time based temperature measurements above 1000 C. This phosphor is therefore an attractive alternative to the much lower luminescence intensity rare-earth-doped thermographic phosphors that are typically utilized at these elevated temperatures. In particular, Cr:GdAlO3 will be preferred over rare-earth-doped phosphors, such as Dy:YAG, at temperatures up to 1200 C for intensity-starved situations when the much lower emission intensity from rare-earth-doped phosphors is insufficient for accurate temperature measurements in the presence of significant radiation background. While transition-metal-doped phosphors such as Cr:Al2O3 (ruby) are known to exhibit high luminescence intensity at low dopant concentrations, quenching due to nonradiative decay pathways competing with the (sup 2)E to (sup 4)A(sub 2) radiative transition (R line) has typically restricted their use for temperature sensing to below 600 C. Thermal quenching of the broadband (sup 4)T(sub 2) to (sup 4)A(sub 2) radiative transition from Cr:GdAlO3, however, is delayed until much higher temperatures (above 1000 C). This spin-allowed broadband emission persists to high temperatures because the lower-lying (sup 2)E energy level acts as a reservoir to thermally populate the higher shorter-lived (sup 4)T(sub 2) energy level and because the activation energy for nonradiative crossover relaxation from the (sup 4)T(sub 2) level to the (sup 4)A(sub 2) ground state is high. The strong crystal field associated with the tight bonding of the AlO6 octahedra in the GdAlO3 perovskite structure is responsible for this behavior

Addendum to "Coherent radio pulses from GEANT generated electromagnetic showers in ice"

We reevaluate our published calculations of electromagnetic showers generated by GEANT 3.21 and the radio frequency pulses they produce in ice. We are prompted by a recent report showing that GEANT 3.21-modeled showers are sensitive to internal settings in the electron tracking subroutine. We report the shower and pulse characteristics obtained with different settings of GEANT 3.21 and with GEANT 4. The default setting of electron tracking in GEANT 3.21 we used in previous work speeds up the shower simulation at the cost of information near the end of the tracks. We find that settings tracking electron and positron to lower energy yield a more accurate calculation, a more intense shower, and proportionately stronger radio pulses at low frequencies. At high frequencies the relation between shower tracking algorithm and pulse spectrum is more complex. We obtain radial distributions of shower particles and phase distributions of pulses from 100 GeV showers that are consistent with our published results.Comment: 4 pages, 3 figure