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.

CCD photometry of distant comets. III - Ensemble properties of Jupiter-family comets

We describe the results of a ground-based observational "snapshot" study of Jupiter-family comets in the heliocentric range 2.29 AU less than or equal to R-h less than or equal to 5.72 AU. Results are presented based on observations from the 1m JKT on the island of La Palma. A total of 25 comets were targeted with 15 being positively detected. Broad-band VRI photometry was performed to determine dimensions, colour indices, and dust production rates in terms of the "A frho" formalism. The results for selected comets are compared with previous investigations. Ensemble properties of the Jupiter- family population have been investigated by combining the results presented here with those of Lowry et al. (1999), and Lowry & Fitzsimmons (2001). We find that the cumulative size distribution of the Jupiter-family comets can be described by a power law of the form; Sigma(> r) proportional to r(-1.6+/- 0.1). This size distribution is considerably shallower than that found for the observed Edgeworth-Kuiper belt objects, which may reflect either an intrinsic difference at small km- sizes in the belt, or the various processes affecting the nuclei of comets as their orbits evolve from the Edgeworth- Kuiper belt to the inner Solar system. Also, there would appear to be no correlation between nuclear absolute magnitude and perihelion distance. Finally, for the sample of active comets, there is a distinct correlation between absolute R band magnitude and perihelion distance, which can be explained by either a discovery bias towards brighter comets or in terms of "rubble" mantle formation.

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.

WASP-1b and WASP-2b: Two new transiting exoplanets detected with SuperWASP and SOPHIE

We have detected low-amplitude radial-velocity variations in two stars, USNO-B1.0 1219-0005465 (GSC 02265-00107 = WASP-1) and USNO-B1.0 0964-0543604 (GSC 00522-01199 = WASP-2). Both stars were identified as being likely host stars of transiting exoplanets in the 2004 SuperWASP wide-field transit survey. Using the newly commissioned radial-velocity spectrograph SOPHIE at the Observatoire de Haute-Provence, we found that both objects exhibit reflex orbital radial-velocity variations with amplitudes characteristic of planetary-mass companions and in-phase with the photometric orbits. Line-bisector studies rule out faint blended binaries as the cause of either the radial-velocity variations or the transits. We perform preliminary spectral analyses of the host stars, which together with their radial-velocity variations and fits to the transit light curves yield estimates of the planetary masses and radii. WASP-1b and WASP-2b have orbital periods of 2.52 and 2.15 d, respectively. Given mass estimates for their F7V and K1V primaries, we derive planet masses 0.80-0.98 and 0.81-0.95 times that of Jupiter, respectively. WASP-1b appears to have an inflated radius of at least 1.33 RJup, whereas WASP-2b has a radius in the range 0.65-1.26 RJup.

WASP-14b: 7.3 M-J transiting planet in an eccentric orbit

We report the discovery of a 7.3 M-J exoplanet WASP-14b, one of the most massive transiting exoplanets observed to date. The planet orbits the 10th-magnitude F5V star USNO-B1 11118-0262485 with a period of 2.243 752 d and orbital eccentricity e = 0.09. A simultaneous fit of the transit light curve and radial velocity measurements yields a planetary mass of 7.3 +/- 0.5 M-J and a radius of 1.28 +/- 0.08 R-J. This leads to a mean density of about 4.6 g cm(-3) making it the densest transiting exoplanets yet found at an orbital period less than 3 d. We estimate this system to be at a distance of 160 +/- 20 pc. Spectral analysis of the host star reveals a temperature of 6475 +/- 100 K, log g = 4.07 cm s(-2) and v sin i = 4.9 +/- 1.0 km s(-1), and also a high lithium abundance, log N(Li) = 2.84 +/- 0.05. The stellar density, effective temperature and rotation rate suggest an age for the system of about 0.5-1.0 Gyr.

A SuperWASP search for additional transiting planets in 24 known systems

We present results from a search for additional transiting planets in 24 systems already known to contain a transiting planet. We model the transits due to the known planet in each system and subtract these models from light curves obtained with the SuperWASP (Wide Angle Search for Planets) survey instruments. These residual light curves are then searched for evidence of additional periodic transit events. Although we do not find any evidence for additional planets in any of the planetary systems studied, we are able to characterize our ability to find such planets by means of Monte Carlo simulations. Artificially generated transit signals corresponding to planets with a range of sizes and orbital periods were injected into the SuperWASP photometry and the resulting light curves searched for planets. As a result, the detection efficiency as a function of both the radius and orbital period of any second planet is calculated. We determine that there is a good (>50 per cent) chance of detecting additional, Saturn-sized planets in P ~ 10 d orbits around planet-hosting stars that have several seasons of SuperWASP photometry. Additionally, we confirm previous evidence of the rotational stellar variability of WASP-10, and refine the period of rotation. We find that the period of the rotation is 11.91 +/- 0.05 d, and the false alarm probability for this period is extremely low (~10-13).

The Low Density Transiting Exoplanet WASP-15b

We report the discovery of a low-density exoplanet transiting an 11th magnitude star in the Southern hemisphere. WASP-15b, which orbits its host star with a period P = 3.7520656 ± 0.0000028 d, has a mass M p = 0.542 ± 0.050 M J and radius R p = 1.428 ± 0.077 R J, and is therefore one of the least dense transiting exoplanets so far discovered (?p = 0.247 ± 0.035 g cm-3). An analysis of the spectrum of the host star shows it to be of spectral type around F5, with an effective temperature T eff = 6300 ± 100 K and [Fe/H] = -0.17 ± 0.11.

Wasp-17b: An Ultra-Low Density Planet in a Probable Retrograde Orbit

We report the discovery of the transiting giant planet WASP-17b, the least-dense planet currently known. It is 1.6 Saturn masses, but 1.5-2 Jupiter radii, giving a density of 6%-14% that of Jupiter. WASP-17b is in a 3.7 day orbit around a sub-solar metallicity, V = 11.6, F6 star. Preliminary detection of the Rossiter-McLaughlin effect suggests that WASP-17b is in a retrograde orbit (? ˜ -150°), indicative of a violent history involving planet-planet or star-planet scattering. WASP-17b's bloated radius could be due to tidal heating resulting from recent or ongoing tidal circularization of an eccentric orbit, such as the highly eccentric orbits that typically result from scattering interactions. It will thus be important to determine more precisely the current orbital eccentricity by further high-precision radial velocity measurements or by timing the secondary eclipse, both to reduce the uncertainty on the planet's radius and to test tidal-heating models. Owing to its low surface gravity, WASP-17b's atmosphere has the largest scale height of any known planet, making it a good target for transmission spectroscopy.

A transit timing analysis of seven RISE light curves of the exoplanet system HAT-P-3

We present seven light curves of the exoplanet system HAT-P-3, taken as part of a transit timing programme using the rapid imager to search for exoplanets instrument on the Liverpool Telescope. The light curves are analysed using a Markov chain Monte Carlo algorithm to update the parameters of the system. The inclination is found to be i = 86.75+0.22-0.21°, the planet-star radius ratio to be Rp/R* = 0.1098+0.0010-0.0012 and the stellar radius to be R* = 0.834+0.018-0.026Rsolar, consistent with previous results but with a significant improvement in the precision. Central transit times and uncertainties for each light curve are also determined, and a residual permutation algorithm is used as an independent check on the errors. The transit times are found to be consistent with a linear ephemeris, and a new ephemeris is calculated as Tc(0) = 2454856.70118 +/- 0.00018 HJD and P = 2.899738 +/- 0.000007 d. Model timing residuals are fitted to the measured timing residuals to place upper mass limits for a hypothetical perturbing planet as a function of the period ratio. These show that we have probed for planets with masses as low as 0.33 and 1.81 M? in the interior and exterior 2:1 resonances, respectively, assuming the planets are initially in circular orbits.

Tight constraints on the existence of additional planets around HD 189733†

We report a transit timing study of the transiting exoplanetary system HD 189733. In total, we observed 10 transits in 2006 and 2008 with the 2.6-m Nordic Optical Telescope, and two transits in 2007 with the 4.2-m William Herschel Telescope. We used Markov Chain Monte Carlo simulations to derive the system parameters and their uncertainties, and our results are in a good agreement with previously published values. We performed two independent analyses of transit timing residuals to place upper mass limits on putative perturbing planets. The results show no evidence for the presence of planets down to 1 Earth mass near the 1:2 and 2:1 resonance orbits, and planets down to 2.2 Earth masses near the 3:5 and 5:3 resonance orbits with HD 189733b. These are the strongest limits to date on the presence of other planets in this system. Based on observations made with the Nordic Optical Telescope, operated on the island of La Palma jointly by Denmark, Finland, Iceland, Norway and Sweden, in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias. Based on observations made with the William Herschel Telescope operated on the island of La Palma by the Isaac Newton Group in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias. ‡

Spin-orbit angle measurements for six southern transiting planets. New insights into the dynamical origins of hot Jupiters

Context. Several competing scenarios for planetary-system formation and evolution seek to explain how hot Jupiters came to be so close to their parent stars. Most planetary parameters evolve with time, making it hard to distinguish between models. The obliquity of an orbit with respect to the stellar rotation axis is thought to be more stable than other parameters such as eccentricity. Most planets, to date, appear aligned with the stellar rotation axis; the few misaligned planets so far detected are massive (> 2 MJ). Aims: Our goal is to measure the degree of alignment between planetary orbits and stellar spin axes, to search for potential correlations with eccentricity or other planetary parameters and to measure long term radial velocity variability indicating the presence of other bodies in the system. Methods: For transiting planets, the Rossiter-McLaughlin effect allows the measurement of the sky-projected angle ß between the stellar rotation axis and a planet's orbital axis. Using the HARPS spectrograph, we observed the Rossiter-McLaughlin effect for six transiting hot Jupiters found by the WASP consortium. We combine these with long term radial velocity measurements obtained with CORALIE. We used a combined analysis of photometry and radial velocities, fitting model parameters with the Markov Chain Monte Carlo method. After obtaining ß we attempt to statistically determine the distribution of the real spin-orbit angle ?. Results: We found that three of our targets have ß above 90°: WASP-2b: ß = 153°+11-15, WASP-15b: ß = 139.6°+5.2-4.3 and WASP-17b: ß = 148.5°+5.1-4.2; the other three (WASP-4b, WASP-5b and WASP-18b) have angles compatible with 0°. We find no dependence between the misaligned angle and planet mass nor with any other planetary parameter. All six orbits are close to circular, with only one firm detection of eccentricity e = 0.00848+0.00085-0.00095 in WASP-18b. No long-term radial acceleration was detected for any of the targets. Combining all previous 20 measurements of ß and our six and transforming them into a distribution of ? we find that between about 45 and 85% of hot Jupiters have ? > 30°. Conclusions: Most hot Jupiters are misaligned, with a large variety of spin-orbit angles. We find observations and predictions using the Kozai mechanism match well. If these observational facts are confirmed in the future, we may then conclude that most hot Jupiters are formed from a dynamical and tidal origin without the necessity to use type I or II migration. At present, standard disc migration cannot explain the observations without invoking at least another additional process.

Estimating the masses of extra-solar planets

All extra-solar planet masses that have been derived spectroscopically are lower limits since the inclination of the orbit to our line-of-sight is unknown except for transiting systems. In theory, however, it is possible to determine the inclination angle, i, between the rotation axis of a star and an observer's line-of-sight from measurements of the projected equatorial velocity (v sin i), the stellar rotation period (P(rot)) and the stellar radius (R(*)). For stars which host planetary systems this allows the removal of the sin i dependency of extra-solar planet masses derived from spectroscopic observations under the assumption that the planetary orbits lie perpendicular to the stellar rotation axis.
We have carried out an extensive literature search and present a catalogue of v sin i, P(rot) and R(*) estimates for stars hosting extra-solar planets. In addition, we have used Hipparcos parallaxes and the Barnes-Evans relationship to further supplement the R(*) estimates obtained from the literature. Using this catalogue, we have obtained sin i estimates using a Markov-chain Monte Carlo analysis. This technique allows proper 1 Sigma two-tailed confidence limits to be placed on the derived sin i's along with the transit probability for each planet to be determined.
While we find that a small proportion of systems yield sin i's significantly greater than 1, most likely due to poor P(rot) estimations, the large majority are acceptable. We are further encouraged by the cases where we have data on transiting systems, as the technique indicates inclinations of similar to 90 degrees and high transit probabilities. In total, we are able to estimate the true masses of 133 extra-solar planets. Of these 133 extra-solar planets, only six have revised masses that place them above the 13M(J) deuterium burning limit; four of those six extra-solar planet candidates were already suspected to lie above the deuterium burning limit before correcting their masses for the sin i dependency. Our work reveals a population of high-mass extra-solar planets with low eccentricities, and we speculate that these extra-solar planets may represent the signature of different planetary formation mechanisms at work. Finally, we discuss future observations that should improve the robustness of this technique.