The Bifurcation Interpreter: A Step Towards the Automatic Analysis of Dynamical Systems

The Bifurcation Interpreter is a computer program that autonomously explores the steady-state orbits of one-parameter families of periodically- driven oscillators. To report its findings, the Interpreter generates schematic diagrams and English text descriptions similar to those appearing in the science and engineering research literature. Given a system of equations as input, the Interpreter uses symbolic algebra to automatically generate numerical procedures that simulate the system. The Interpreter incorporates knowledge about dynamical systems theory, which it uses to guide the simulations, to interpret the results, and to minimize the effects of numerical error.

PROTRACTED HYPOFRACTIONATED RADIOTHERAPY for GRAVES' OPHTHALMOPATHY: A PILOT STUDY of CLINICAL and RADIOLOGIC RESPONSE

Purpose: To evaluate the clinical and radiologic response of patients with Graves' ophthalmopathy given low-dose orbital radiotherapy (RT) with a protracted fractionation.Methods and Materials: Eighteen patients (36 orbits) received orbital RT with a total dose of 10 Gy, fractionated in 1 Gy once a week over 10 weeks. of these, 9 patients received steroid therapy as well. Patients were evaluated clinically and radiologically at 6 months after treatment. Clinical response assessment was carried out using three criteria: by physical examination, by a modified clinical activity score, and by a verbal questionnaire considering the 10 most common signs and symptoms of the disease. Radiologic response was assessed by magnetic resonance imaging.Results: Improvement in ocular pain, palpebral edema, visual acuity, and ocular motility was observed in all patients. Significant decrease in symptoms such as tearing (p < 0.001) diplopia (p = 0.008), conjunctival hyperemia (p = 0.002), and ocular grittiness (p = 0.031) also occurred. Magnetic resonance imaging showed decrease in ocular muscle thickness and in the intensity of the T2 sequence signal in the majority of patients. Treatments were well tolerated, and to date no complications from treatment have been observed. There was no statistical difference in clinical and radiologic response between patients receiving RT alone and those receiving RT plus steroid therapy.Conclusion: RT delivered in at a low dose and in a protracted scheme should be considered as a useful therapeutic option for patients with Graves' ophthalmopathy. (C) 2012 Elsevier Inc.

The Depiction of Coronal Structure in White-Light Images

Morgan, Huw, Habbal, S.R., Woo, R., (2006) 'The Depiction of Coronal Structure in White-Light Images', Solar Physics 236(2) pp.263-272 RAE2008

Ocena wiarygodności metod identyfikacji strumieni meteoroidowych

Wydział Fizyki:Instytut Obserwatorium Astronomiczne

The Construction of Arbitrary Stable Dynamics in Non-Linear Neural Networks

In this paper, two methods for constructing systems of ordinary differential equations realizing any fixed finite set of equilibria in any fixed finite dimension are introduced; no spurious equilibria are possible for either method. By using the first method, one can construct a system with the fewest number of equilibria, given a fixed set of attractors. Using a strict Lyapunov function for each of these differential equations, a large class of systems with the same set of equilibria is constructed. A method of fitting these nonlinear systems to trajectories is proposed. In addition, a general method which will produce an arbitrary number of periodic orbits of shapes of arbitrary complexity is also discussed. A more general second method is given to construct a differential equation which converges to a fixed given finite set of equilibria. This technique is much more general in that it allows this set of equilibria to have any of a large class of indices which are consistent with the Morse Inequalities. It is clear that this class is not universal, because there is a large class of additional vector fields with convergent dynamics which cannot be constructed by the above method. The easiest way to see this is to enumerate the set of Morse indices which can be obtained by the above method and compare this class with the class of Morse indices of arbitrary differential equations with convergent dynamics. The former set of indices are a proper subclass of the latter, therefore, the above construction cannot be universal. In general, it is a difficult open problem to construct a specific example of a differential equation with a given fixed set of equilibria, permissible Morse indices, and permissible connections between stable and unstable manifolds. A strict Lyapunov function is given for this second case as well. This strict Lyapunov function as above enables construction of a large class of examples consistent with these more complicated dynamics and indices. The determination of all the basins of attraction in the general case for these systems is also difficult and open.

The Synthesis of Arbitrary Stable Dynamics in Non-linear Neural Networks II: Feedback and Universality

We wish to construct a realization theory of stable neural networks and use this theory to model the variety of stable dynamics apparent in natural data. Such a theory should have numerous applications to constructing specific artificial neural networks with desired dynamical behavior. The networks used in this theory should have well understood dynamics yet be as diverse as possible to capture natural diversity. In this article, I describe a parameterized family of higher order, gradient-like neural networks which have known arbitrary equilibria with unstable manifolds of known specified dimension. Moreover, any system with hyperbolic dynamics is conjugate to one of these systems in a neighborhood of the equilibrium points. Prior work on how to synthesize attractors using dynamical systems theory, optimization, or direct parametric. fits to known stable systems, is either non-constructive, lacks generality, or has unspecified attracting equilibria. More specifically, We construct a parameterized family of gradient-like neural networks with a simple feedback rule which will generate equilibrium points with a set of unstable manifolds of specified dimension. Strict Lyapunov functions and nested periodic orbits are obtained for these systems and used as a method of synthesis to generate a large family of systems with the same local dynamics. This work is applied to show how one can interpolate finite sets of data, on nested periodic orbits.

The synchronization of superparamagnetic beads driven by a micro-magnetic ratchet.

We present theoretical, numerical, and experimental analyses on the non-linear dynamic behavior of superparamagnetic beads exposed to a periodic array of micro-magnets and an external rotating field. The agreement between theoretical and experimental results revealed that non-linear magnetic forcing dynamics are responsible for transitions between phase-locked orbits, sub-harmonic orbits, and closed orbits, representing different mobility regimes of colloidal beads. These results suggest that the non-linear behavior can be exploited to construct a novel colloidal separation device that can achieve effectively infinite separation resolution for different types of beads, by exploiting minor differences in their bead's properties. We also identify a unique set of initial conditions, which we denote the "devil's gate" which can be used to expeditiously identify the full range of mobility for a given bead type.

Contrasting the efficiency of radiation belt losses caused by ducted and nonducted whistler-mode waves from ground-based transmitters

It has long been recognized that whistler-mode waves can be trapped in plasmaspheric whistler ducts which guide the waves. For nonguided cases these waves are said to be "nonducted", which is dominant for L < 1.6. Wave-particle interactions are affected by the wave being ducted or nonducted. In the field-aligned ducted case, first-order cyclotron resonance is dominant, whereas nonducted interactions open up a much wider range of energies through equatorial and off-equatorial resonance. There is conflicting information as to whether the most significant particle loss processes are driven by ducted or nonducted waves. In this study we use loss cone observations from the DEMETER and POES low-altitude satellites to focus on electron losses driven by powerful VLF communications transmitters. Both satellites confirm that there are well-defined enhancements in the flux of electrons in the drift loss cone due to ducted transmissions from the powerful transmitter with call sign NWC. Typically, ∼80% of DEMETER nighttime orbits to the east of NWC show electron flux enhancements in the drift loss cone, spanning a L range consistent with first-order cyclotron theory, and inconsistent with nonducted resonances. In contrast, ∼1% or less of nonducted transmissions originate from NPM-generated electron flux enhancements. While the waves originating from these two transmitters have been predicted to lead to similar levels of pitch angle scattering, we find that the enhancements from NPM are at least 50 times smaller than those from NWC. This suggests that lower-latitude, nonducted VLF waves are much less effective in driving radiation belt pitch angle scattering. Copyright 2010 by the American Geophysical Union.

Resonant and Secular Orbital Interations

In stable solar systems, planets remain in nearly elliptical orbits around their stars. Over longer timescales, however, their orbital shapes and sizes change due to mutual gravitational perturbations. Orbits of satellites around a planet vary for the same reason. Because of their interactions, the orbits of planets and satellites today are different from what they were earlier. In order to determine their original orbits, which are critical constraints on formation theories, it is crucial to understand how orbits evolve over the age of the Solar System. Depending on their timescale, we classify orbital interactions as either short-term (orbital resonances) or long-term (secular evolution). My work involves examples of both interaction types. Resonant history of the small Neptunian satellites In satellite systems, tidal migration brings satellite orbits in and out of resonances. During a resonance passage, satellite orbits change dramatically in a very short period of time. We investigate the resonant history of the six small Neptunian moons. In this unique system, the exotic orbit of the large captured Triton (with a circular, retrograde, and highly tilted orbit) influences the resonances among the small satellites very strongly. We derive an analytical framework which can be applied to Neptune's satellites and to similar systems. Our numerical simulations explain the current orbital tilts of the small satellites as well as constrain key physical parameters of both Neptune and its moons. Secular orbital interactions during eccentricity damping Long-term periodic changes of orbital shape and orientation occur when two or more planets orbit the same star. The variations of orbital elements are superpositions of the same number of fundamental modes as the number of planets in the system. We investigate how this effect interacts with other perturbations imposed by external disturbances, such as the tides and relativistic effects. Through analytical studies of a system consisting of two planets, we find that an external perturbation exerted on one planet affects the other indirectly. We formulate a general theory for how both orbits evolve in response to an arbitrary externally-imposed slow change in eccentricity.

Investigation and Development of a Fully 3D Tilt Capable Hybrid SPECT - CT System for Dedicated Breast Imaging

X-ray mammography has been the gold standard for breast imaging for decades, despite the significant limitations posed by the two dimensional (2D) image acquisitions. Difficulty in diagnosing lesions close to the chest wall and axilla, high amount of structural overlap and patient discomfort due to compression are only some of these limitations. To overcome these drawbacks, three dimensional (3D) breast imaging modalities have been developed including dual modality single photon emission computed tomography (SPECT) and computed tomography (CT) systems. This thesis focuses on the development and integration of the next generation of such a device for dedicated breast imaging. The goals of this dissertation work are to: [1] understand and characterize any effects of fully 3-D trajectories on reconstructed image scatter correction, absorbed dose and Hounsifeld Unit accuracy, and [2] design, develop and implement the fully flexible, third generation hybrid SPECT-CT system capable of traversing complex 3D orbits about a pendant breast volume, without interference from the other. Such a system would overcome artifacts resulting from incompletely sampled divergent cone beam imaging schemes and allow imaging closer to the chest wall, which other systems currently under research and development elsewhere cannot achieve.

The dependence of x-ray scatter radiation on object shape, size, material composition and the CT acquisition trajectory, was investigated with a well-established beam stop array (BSA) scatter correction method. While the 2D scatter to primary ratio (SPR) was the main metric used to characterize total system scatter, a new metric called ‘normalized scatter contribution’ was developed to compare the results of scatter correction on 3D reconstructed volumes. Scatter estimation studies were undertaken with a sinusoidal saddle (±15° polar tilt) orbit and a traditional circular (AZOR) orbit. Clinical studies to acquire data for scatter correction were used to evaluate the 2D SPR on a small set of patients scanned with the AZOR orbit. Clinical SPR results showed clear dependence of scatter on breast composition and glandular tissue distribution, otherwise consistent with the overall phantom-based size and density measurements. Additionally, SPR dependence was also observed on the acquisition trajectory where 2D scatter increased with an increase in the polar tilt angle of the system.

The dose delivered by any imaging system is of primary importance from the patient’s point of view, and therefore trajectory related differences in the dose distribution in a target volume were evaluated. Monte Carlo simulations as well as physical measurements using radiochromic film were undertaken using saddle and AZOR orbits. Results illustrated that both orbits deliver comparable dose to the target volume, and only slightly differ in distribution within the volume. Simulations and measurements showed similar results, and all measured dose values were within the standard screening mammography-specific, 6 mGy dose limit, which is used as a benchmark for dose comparisons.

Hounsfield Units (HU) are used clinically in differentiating tissue types in a reconstructed CT image, and therefore the HU accuracy of a system is very important, especially when using non-traditional trajectories. Uniform phantoms filled with various uniform density fluids were used to investigate differences in HU accuracy between saddle and AZOR orbits. Results illustrate the considerably better performance of the saddle orbit, especially close to the chest and nipple region of what would clinically be a pedant breast volume. The AZOR orbit causes shading artifacts near the nipple, due to insufficient sampling, rendering a major portion of the scanned phantom unusable, whereas the saddle orbit performs exceptionally well and provides a tighter distribution of HU values in reconstructed volumes.

Finally, the third generation, fully-suspended SPECT-CT system was designed in and developed in our lab. A novel mechanical method using a linear motor was developed for tilting the CT system. A new x-ray source and a custom made 40 x 30 cm2 detector were integrated on to this system. The SPECT system was nested, in the center of the gantry, orthogonal to the CT source-detector pair. The SPECT system tilts on a goniometer, and the newly developed CT tilting mechanism allows ±15° maximum polar tilting of the CT system. The entire gantry is mounted on a rotation stage, allowing complex arbitrary trajectories for each system, without interference from the other, while having a common field of view. This hybrid system shows potential to be used clinically as a diagnostic tool for dedicated breast imaging.

Proper-motion binaries in the Hipparcos catalogue: comparison with radial velocity data

Context. This paper is the last in a series devoted to the analysis of the binary content of the Hipparcos Catalogue. Aims. The comparison of the proper motions constructed from positions spanning a short (Hipparcos) or long time (Tycho-2) makes it possible to uncover binaries with periods of the order of or somewhat larger than the short time span (in this case, the 3 yr duration of the Hipparcos mission), since the unrecognised orbital motion will then add to the proper motion. Methods. A list of candidate proper motion binaries is constructed from a carefully designed χ2 test evaluating the statistical significance of the difference between the Tycho-2 and Hipparcos proper motions for 103 134 stars in common between the two catalogues (excluding components of visual systems). Since similar lists of proper-motion binaries have already been constructed, the present paper focuses on the evaluation of the detection efficiency of proper-motion binaries, using different kinds of control data (mostly radial velocities). The detection rate for entries from the Ninth Catalogue of Spectroscopic Binary Orbits (SB9) is evaluated, as well as for stars like barium stars, which are known to be all binaries, and finally for spectroscopic binaries identified from radial velocity data in the Geneva-Copenhagen survey of F and G dwarfs in the solar neighbourhood. Results. Proper motion binaries are efficiently detected for systems with parallaxes in excess of ∼20 mas, and periods in the range 1000-30 000 d. The shortest periods in this range (1000-2000 d, i.e. once to twice the duration of the Hipparcos mission) may appear only as DMSA/G binaries (accelerated proper motion in the Hipparcos Double and Multiple System Annex). Proper motion binaries detected among SB9 systems having periods shorter than about 400 d hint at triple systems, the proper-motion binary involving a component with a longer orbital period. A list of 19 candidate triple systems is provided. Binaries suspected of having low-mass (brown-dwarf-like) companions are listed as well. Among the 37 barium stars with parallaxes larger than 5 mas, only 7 exhibit no evidence for duplicity whatsoever (be it spectroscopic or astrometric). Finally, the fraction of proper-motion binaries shows no significant variation among the various (regular) spectral classes, when due account is taken for the detection biases. © ESO 2007.

Chemically tagging the Hyades stream: does it partly originate from the Hyades cluster?

The Hyades stream has long been thought to be a dispersed vestige of the Hyades cluster. However, recent analyses of the parallax distribution, of the mass function, and of the action-space distribution of stream stars have shown it to be rather composed of orbits trapped at a resonance of a density disturbance. This resonant scenario should leave a clearly different signature in the element abundances of stream stars than the dispersed cluster scenario, since the Hyades cluster is chemically homogeneous. Here, we study the metallicity as well as the element abundances of Li, Na, Mg, Fe, Zr, Ba, La, Ce, Nd and Eu for a random sample of stars belonging to the Hyades stream, and compare them with those of stars from the Hyades cluster. From this analysis: (i) we independently confirm that the Hyades stream cannot be solely composed of stars originating in the Hyades cluster; (ii) we show thatsomestars (namely 2/21) from the Hyades stream nevertheless have abundances compatible with an origin in the cluster; (iii) we emphasize that the use of Li as a chemical tag of the cluster origin of main-sequence stars is very efficient in the range 5500K ≤Teff≤ 6200K, since the Li sequence in the Hyades cluster is very tight, while at the same time spanning a large abundance range; (iv) we show that, while this evaporated population has a metallicity excess of ~0.2 dex with respect to the local thin-disc population, identical to that of the Hyades cluster, the remainder of the Hyades stream population has still a metallicity excess of ~0.06-0.15 dex, consistent with an origin in the inner Galaxy and (v) we show that the Hyades stream can be interpreted as an inner 4:1 resonance of the spiral pattern: this then also reproduces an orbital family compatible with the Sirius stream, and places the origin of the Hyades stream up to 1kpc inwards from the solar radius, which might explain the observed metallicity excess of the stream population. © 2011 The Authors Monthly Notices of the Royal Astronomical Society © 2011 RAS.

Positronium-atom collisions

New results are presented for Ps(1s) scattering by H(1s), He(1(1)S) and Li(2s). Calculations have been performed in a coupled state framework, usually employing pseudostates, and allowing for excitation of both the Ps and the atom. In the Ps(1s)-H(1s) calculations the H- formation channel has also been included using a highly accurate H- wave function. Resonances resulting from unstable states in which the positron orbits H- have been calculated and analysed. The new Ps(1s)-He(1(1)S) calculations still fail to resolve existing discrepancies between theory and experiment at very low energies. The possible importance of the Ps(-) formation channel in all three collision systems is discussed. (C) 2004 Elsevier B.V. All rights reserved.

Optical colors of 56 near-Earth objects: trends with size and orbit

We present the results of BVRIZ photometry of 56 near-Earth objects (NEOs) obtained with the 1-m Jacobus Kapteyn telescope on La Palma during 2000 and 2001. Our sample includes many NEOs with particularly deep 1 - mum pyroxene/olivine absorption bands, similar to Q-type asteroids. We also classify three NEOs with particularly blue colors. No D-type asteroids were found, placing an upper limit of similar to2% on the fraction of the NEO population originating in the outer main belt or the Trojan clouds. The ratio of dark to bright objects in our sample was found to be 0.40, significantly higher than current theoretical predictions. As well as classifying the NEOs, we have investigated color trends with size and orbit. We see a general trend for larger silicate objects to have shallower absorption bands but find no significant difference in the distribution of taxonomic classes at small and large sizes. Our data clearly show that different taxonomic classes tend to occupy different regions of (a, e) space. By comparing our data with current model predictions for NEO dynamical evolution we see that Q- R-, and V-type NEOs tend to have orbits associated with "fast track" delivery from the main belt, whereas S-type NEOs tend to have orbits associated with "slow track" delivery. This outcome would be expected if space weathering occurs on time scales of >10(6) years. (C) 2003 Elsevier Science (USA). All rights reserved.

CCD photometry of distant comets II

R band CCD photometric observations of short period Jupiter family comets in the heliocentric region of 2.11 AU less than or equal to R-h less than or equal to 5.63 AU were performed using the 4.2 m William Herschel Telescope on La Palma in December 1998. 22 comets were targeted, including the comet- asteroid transition object 49P/Arend-Rigaux. Out of a total of ten detected comets. six were seen to display substantial outgassing (48P, 65P, 74P, 103Pt 128P, and 139P), with the remaining four comets (7P. 9P, 22P. and 49P) being stellar in appearance. Nuclear radius measurements and relative dust production rates in terms of Af rho were measured for these comets, along with upper limits for the remaining twelve undetected comets (6P, 44P, 51P, 54P, 57P: 63P, 71P, 73P, 79P, 86P, 87P, and 100P). The inactive comets had nuclear radii in the range 1.8 km less than or equal to r(N) less than or equal to 4.4 km, while upper limits for the active and undetected comets (assuming they all lay within the held of view) were between 0.6 km and 12.7 km, for an assumed albedo of 0.04. Even if one applies the previously measured maximum axis ratio of 2.6:1 and the minimum measured albedo of 0.02 to the undetected comets. their projected semi-major axes are all constrained to below 8 km. For the active comets, photometric profiles of their dust comae were measured and are consistent with those of steady state coma models.