Three Hundred Years in Annapolis: Phase III Archaeological Investigations of the Anne Arundel County Courthouse Site (18AP63), Annapolis, Maryland

During the summer of 1994, Archaeology in Annapolis conducted archaeological investigations of the city block bounded by Franklin, South and Cathedral Streets in the city of Annapolis. This Phase III excavation was conducted as a means to identify subsurface cultural resources in the impact area associated with the proposed construction of the Anne Arundel County Courthouse addition. This impact area included both the upper and lower parking lots used by Courthouse employees. Investigations were conducted in the form of mechanical trenching and hand excavated units. Excavations in the upper lot area yielded significant information concerning the interior area of the block. Known as Bellis Court, this series of rowhouses was constructed in the late nineteenth century and was used as rental properties by African-Americans. The dwellings remained until the middle of the twentieth century when they were demolished in preparation for the construction of a Courthouse addition. Portions of the foundation of a house owned by William H. Bellis in the 1870s were also exposed in this area. Construction of this house was begun by William Nicholson around 1730 and completed by Daniel Dulany in 1732/33. It was demolished in 1896 by James Munroe, a Trustee for Bellis. Excavations in the upper lot also revealed the remains of a late seventeenth/early eighteenth century wood-lined cellar, believed to be part of the earliest known structure on Lot 58. After an initially rapid deposition of fill around 1828, this cellar was gradually covered with soil throughout the remainder of the nineteenth century. The fill deposit in the cellar feature yielded a mixed assemblage of artifacts that included sherds of early materials such as North Devon gravel-tempered earthenware, North Devon sgraffito and Northem Italian slipware, along with creamware, pearlware and whiteware. In the lower parking lot, numerous artifacts were recovered from yard scatter associated with the houses that at one time fronted along Cathedral Street and were occupied by African- Americans. An assemblage of late seventeenth century/early eighteenth century materials and several slag deposits from an early forge were recovered from this second area of study. The materials associated with the forge, including portions of a crucible, provided evidence of some of the earliest industry in Annapolis. Investigations in both the upper and lower parking lots added to the knowledge of the changing landscape within the project area, including a prevalence of open space in early periods, a surprising survival of impermanent structures, and a gradual regrading and filling of the block with houses and interior courts. Excavations at the Anne Arundel County Courthouse proved this to be a multi-component site, rich in cultural resources from Annapolis' Early Settlement Period through its Modern Period (as specified by Maryland's Comprehensive Historic Preservation Plan (Weissman 1986)). This report provides detailed interpretations of the archaeological findings of these Phase III investigations.

Computational Journalism: from Answering Question to Questioning Answers and Raising Good Questions

Our media is saturated with claims of ``facts'' made from data. Database research has in the past focused on how to answer queries, but has not devoted much attention to discerning more subtle qualities of the resulting claims, e.g., is a claim ``cherry-picking''? This paper proposes a Query Response Surface (QRS) based framework that models claims based on structured data as parameterized queries. A key insight is that we can learn a lot about a claim by perturbing its parameters and seeing how its conclusion changes. This framework lets us formulate and tackle practical fact-checking tasks --- reverse-engineering vague claims, and countering questionable claims --- as computational problems. Within the QRS based framework, we take one step further, and propose a problem along with efficient algorithms for finding high-quality claims of a given form from data, i.e. raising good questions, in the first place. This is achieved to using a limited number of high-valued claims to represent high-valued regions of the QRS. Besides the general purpose high-quality claim finding problem, lead-finding can be tailored towards specific claim quality measures, also defined within the QRS framework. An example of uniqueness-based lead-finding is presented for ``one-of-the-few'' claims, landing in interpretable high-quality claims, and an adjustable mechanism for ranking objects, e.g. NBA players, based on what claims can be made for them. Finally, we study the use of visualization as a powerful way of conveying results of a large number of claims. An efficient two stage sampling algorithm is proposed for generating input of 2d scatter plot with heatmap, evalutaing a limited amount of data, while preserving the two essential visual features, namely outliers and clusters. For all the problems, we present real-world examples and experiments that demonstrate the power of our model, efficiency of our algorithms, and usefulness of their results.

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.

Comparisons between “sand blast” and “centripetal effect accelerator” type erosion testers

There are two major types of erosion testing devices that are used throughout the world for quantifying particle impact erosion against a solid surface. The first of these uses pressurised air to accelerate abrasive particles through a nozzle so that they impinge upon a target specimen. The second adopts a rotating disc to accelerate abrasive particles using the centripetal effect so that they impinge upon a series of targets arranged around the periphery of the disc. This paper reports the findings of a collaborative project that was designed to compare the performance and results obtained from a rig of each of the two types mentioned above. The sand blast type rig was provided by The Department of Powder Science Technology (POSTEC) at The Telemark Technological Research and Development Centre (TEL-TEK), Porsgrunn, Norway while the centripetal effect accelerator was provided by The Wolfson Centre for Bulk Solids Handling Technology, University of Greenwich, London, UK. The test programme included tests against a wide range of materials that are commonly used in pneumatic handling facilities. (Pneumatic handling is a means of conveying and transporting powders and granular solid materials in bulk in industrial process plant, through pipelines using a gas as the carrier medium.) Olivine sand was used as the abrasive and it was projected against the test specimens at velocities and concentrations commensurate with those seen in pneumatic conveyors. In all instances the materials used in the test programme were taken from the same batch so that scatter of experimental results due to specimen variation was minimised. The paper contains a series of recommendations for erosion testing equipment. A discussion based on the results and their applicability to the prediction of wear in pneumatic conveyors concludes the paper.

Dimethylsulfide gas transfer coefficients from algal blooms in the Southern Ocean

Air–sea dimethylsulfide (DMS) fluxes and bulk air–sea gradients were measured over the Southern Ocean in February–March 2012 during the Surface Ocean Aerosol Production (SOAP) study. The cruise encountered three distinct phytoplankton bloom regions, consisting of two blooms with moderate DMS levels, and a high biomass, dinoflagellate-dominated bloom with high seawater DMS levels (> 15 nM). Gas transfer coefficients were considerably scattered at wind speeds above 5 m/s. Bin averaging the data resulted in a linear relationship between wind speed and mean gas transfer velocity consistent with that previously observed. However, the wind-speed-binned gas transfer data distribution at all wind speeds is positively skewed. The flux and seawater DMS distributions were also positively skewed, which suggests that eddy covariance-derived gas transfer velocities are consistently influenced by additional, log-normal noise. A flux footprint analysis was conducted during a transect into the prevailing wind and through elevated DMS levels in the dinoflagellate bloom. Accounting for the temporal/spatial separation between flux and seawater concentration significantly reduces the scatter in computed transfer velocity. The SOAP gas transfer velocity data show no obvious modification of the gas transfer–wind speed relationship by biological activity or waves. This study highlights the challenges associated with eddy covariance gas transfer measurements in biologically active and heterogeneous bloom environments.

Fault Diagnosis in Internal Combustion Engines Using Nonlinear Multivariate Statistics

This paper presents a statistical-based fault diagnosis scheme for application to internal combustion engines. The scheme relies on an identified model that describes the relationships between a set of recorded engine variables using principal component analysis (PCA). Since combustion cycles are complex in nature and produce nonlinear relationships between the recorded engine variables, the paper proposes the use of nonlinear PCA (NLPCA). The paper further justifies the use of NLPCA by comparing the model accuracy of the NLPCA model with that of a linear PCA model. A new nonlinear variable reconstruction algorithm and bivariate scatter plots are proposed for fault isolation, following the application of NLPCA. The proposed technique allows the diagnosis of different fault types under steady-state operating conditions. More precisely, nonlinear variable reconstruction can remove the fault signature from the recorded engine data, which allows the identification and isolation of the root cause of abnormal engine behaviour. The paper shows that this can lead to (i) an enhanced identification of potential root causes of abnormal events and (ii) the masking of faulty sensor readings. The effectiveness of the enhanced NLPCA based monitoring scheme is illustrated by its application to a sensor fault and a process fault. The sensor fault relates to a drift in the fuel flow reading, whilst the process fault relates to a partial blockage of the intercooler. These faults are introduced to a Volkswagen TDI 1.9 Litre diesel engine mounted on an experimental engine test bench facility.

A note on the prediction of liquid hold-up with the stratified roll wave regime for gas/liquidco-current flow in horizontal pipes.

This note presents a simple model for prediction of liquid hold-up in two-phase horizontal pipe flow for the stratified roll wave (St+RW) flow regime. Liquid hold-up data for horizontal two-phase pipe flow [1, 2, 3, 4, 5 and 6] exhibit a steady increase with liquid velocity and a more dramatic fall with increasing gas rate as shown by Hand et al. [7 and 8] for example. In addition the liquid hold-up is reported to show an additional variation with pipe diameter. Generally, if the initial liquid rate for the no-gas flow condition gives a liquid height below the pipe centre line, the flow patterns pass successively through the stratified (St), stratified ripple (St+R), stratified roll wave, film plus droplet (F+D) and finally the annular (A+D, A+RW, A+BTS) regimes as the gas rate is increased. Hand et al. [7 and 8] have given a detailed description of this progression in flow regime development and definitions of the patterns involved. Despite the fact that there are over one hundred models which have been developed to predict liquid hold-up, none have been shown to be universally useful, while only a handful have proven to be applicable to specific flow regimes [9, 10, 11 and 12]. One of the most intractable regimes to predict has been the stratified roll wave pattern where the liquid hold-up shows the most dramatic change with gas flow rate. It has been suggested that the momentum balance-type models, which give both hold-up and pressure drop prediction, can predict universally for all flow regimes but particularly in the case of the difficult stratified roll wave pattern. Donnelly [1] recently demonstrated that the momentum balance models experienced some difficulties in the prediction of this regime. Without going into lengthy details, these models differ in the assumed friction factor or shear stress on the surfaces within the pipe particularly at the liquid–gas interface. The Baker–Jardine model [13] when tested against the 0.0454 m i.d. data of Nguyen [2] exhibited a wide scatter for both liquid hold-up and pressure drop as shown in Fig. 1. The Andritsos–Hanratty model [14] gave better prediction of pressure drop but a wide scatter for liquid hold-up estimation (cf. Fig. 2) when tested against the 0.0935 m i.d. data of Hand [5]. The Spedding–Hand model [15], shown in Fig. 3 against the data of Hand [5], gave improved performance but was still unsatisfactory with the prediction of hold-up for stratified-type flows. The MARS model of Grolman [6] gave better prediction of hold-up (cf. Fig. 4) but deterioration in the estimation of pressure drop when tested against the data of Nguyen [2]. Thus no method is available that will accurately predict liquid hold-up across the whole range of flow patterns but particularly for the stratified plus roll wavy regime. The position is particularly unfortunate since the stratified-type regimes are perhaps the most predominant pattern found in multiphase lines.

Virtual state scattering with cold electrons: para-xylene and para-difluorobenzene

The scattering of electrons with kinetic energies down to a few meV by para-xylene and para-difluorobenzene has been observed experimentally with an electron beam energy resolution of 0.95 to 1.5 meV (full width half maximum). At low electron energies the collisions can be considered as cold scattering events because the de Broglie wavelength of the electron is considerably larger than the target dimensions. The scattering cross sections measured rise rapidly at low energy due to virtual state scattering. The nature of this scattering process is discussed using s- and p-wave phase shifts derived from the experimental data. Scattering lengths are derived of, respectively, -9.5+/-0.5 and -8.0+/-0.5 a.u. for para-xylene and para-difluorobenzene. The virtual state effect is interpreted in terms of nuclear diabatic and partially adiabatic models, involving the electronic and vibronic symmetries of the unoccupied orbitals in the target species. The concept of direct and indirect virtual state scattering is introduced, through which the present species, in common with carbon dioxide and benzene, scatter through an indirect virtual state process, whereas other species, such as perfluorobenzene, scatter through a direct process. (C) 2005 American Institute of Physics.

Clinical implementation of dynamic multileaf collimation for compensated bladder treatments.

BACKGROUND AND PURPOSE: To describe the clinical implementation of dynamic multileaf collimation (DMLC). Custom compensated four-field treatments of carcinoma of the bladder have been used as a simple test site for the introduction of intensity modulated radiotherapy.MATERIALS AND METHODS: Compensating intensity modulations are calculated from computed tomography (CT) data, accounting for scattered, as well as primary radiation. Modulations are converted to multileaf collimator (MLC) leaf and jaw settings for dynamic delivery on a linear accelerator. A full dose calculation is carried out, accounting for dynamic leaf and jaw motion and transmission through these components. Before treatment, a test run of the delivery is performed and an absolute dose measurement made in a water or solid water phantom. Treatments are verified by in vivo diode measurements and real-time electronic portal imaging. RESULTS: Seven patients have been treated using DMLC. The technique improves dose homogeneity within the target volume, reducing high dose areas and compensating for loss of scatter at the beam edge. A typical total treatment time is 20 min. CONCLUSIONS: Compensated bladder treatments have proven an effective test site for DMLC in an extremely busy clinic.

Direct Observation of strong coupling in a laser-shock compressed plasma

In this Letter we report on a near collective x-ray scattering experiment on shock-compressed targets. A highly coupled Al plasma was generated and probed by spectrally resolving an x-ray source forward scattered by the sample. A significant reduction in the intensity of the elastic scatter was observed, which we attribute to the formation of an incipient long-range order. This speculation is confirmed by x-ray scattering calculations accounting for both electron degeneracy and strong coupling effects. Measurements from rear side visible diagnostics are consistent with the plasma parameters inferred from x-ray scattering data. These results give the experimental evidence of the strongly coupled ionic dynamics in dense plasmas.

Wide angle crystal spectrometer for angularly and spectrally resolved x-ray scattering experiments

A novel wide angle spectrometer has been implemented with a highly oriented pyrolytic graphite crystal coupled to an image plate. This spectrometer has allowed us to look at the energy resolved spectrum of scattered x rays from a dense plasma over a wide range of angles ( ~ 30°) in a single shot. Using this spectrometer we were able to observe the temporal evolution of the angular scatter cross section from a laser shocked foil. A spectrometer of this type may also be useful in investigations of x-ray line transfer from laser-plasmas experiments.

Testing the reliability of the JEOL FEGSEM 6500F electron microprobe for quantitative major element analysis of glass shards from rhyolitic tephra.

Electronprobe microanalysis is now widely adopted in tephra studies as a technique for determining the major element geochemistry of individual glass shards. Accurate geochemical characterization is crucial for enabling robust tephra-based correlations; such information may also be used to link the tephra to a specific source and often to a particular eruption. In this article, we present major element analyses for rhyolitic natural glass standards analysed on three different microprobes and the new JEOL FEGSEM 6500F microprobe at Queen’s University Belfast. Despite the scatter in some elements, good comparability is demonstrated among data yielded from this new system, the previous Belfast JEOL-733 Superprobe, the JEOL-8200 Superprobe (Copenhagen) and the existing long-established microprobe facility in Edinburgh. Importantly, our results show that major elements analysed using different microprobes and variable operating conditions allow two high-silica glasses to be discriminated accurately.

Interaction of High Intensity Picosecond Laser Pulses with Large Pre-formed Plasmas

The interaction of short (1-2 ps) laser pulses with solid targets at irradiances of over 1016 Wcm~2 , in the presence of a substantial prepulse has been investigated. High absorption of laser energy is found even at high angles of incidence, with evidence for a resonance absorption peak being found for S, P, and circular polarizations. It is considered that this may be a result of refraction and beam filamentation, which causes loss of distinct polarization. Measurements of hard X-ray emission (~ 100 keV) confirm a resonance absorption type peak at 45-50°, again for all three cases. Typically, 5-15% of the incident light is back-reflected by stimulated Brillouin scatter, with spatially resolved spectra showing evidence of beam hot-spots at high intensity. The possibility that filamentation and refraction of the beam can explain the lack of polarization dependence in the absorption and hard X-ray emission data is discussed.

The main-sequence rotation-colour relation in the Coma Berenices open cluster

We present the results of a photometric survey of rotation rates in the Coma Berenices (Melotte 111) open cluster, using data obtained as part of the SuperWASP exoplanetary transit-search programme. The goal of the Coma survey was to measure precise rotation periods for main-sequence F, G and K dwarfs in this intermediate-age (~600 Myr) cluster, and to determine the extent to which magnetic braking has caused the stellar spin periods to converge. We find a tight, almost linear relationship between rotation period and J - K colour with an rms scatter of only 2 per cent. The relation is similar to that seen among F, G and K stars in the Hyades. Such strong convergence can only be explained if angular momentum is not at present being transferred from a reservoir in the deep stellar interiors to the surface layers. We conclude that the coupling time-scale for angular momentum transport from a rapidly spinning radiative core to the outer convective zone must be substantially shorter than the cluster age, and that from the age of Coma onwards stars rotate effectively as solid bodies. The existence of a tight relationship between stellar mass and rotation period at a given age supports the use of stellar rotation period as an age indicator in F, G and K stars of Hyades age and older. We demonstrate that individual stellar ages can be determined within the Coma population with an internal precision of the order of 9 per cent (rms), using a standard magnetic braking law in which rotation period increases with the square root of stellar age. We find that a slight modification to the magnetic-braking power law, P ~ t0.56, yields rotational and asteroseismological ages in good agreement for the Sun and other stars of solar age for which p-mode studies and photometric rotation periods have been published.

Quantum ground state of self-organized atomic crystals in optical resonators

Cold atoms, driven by a laser and simultaneously coupled to the quantum field of an optical resonator, may self-organize in periodic structures. These structures are supported by the optical lattice, which emerges from the laser light they scatter into the cavity mode and form when the laser intensity exceeds a threshold value. We study theoretically the quantum ground state of these structures above the pump threshold of self-organization by mapping the atomic dynamics of the self-organized crystal to a Bose-Hubbard model. We find that the quantum ground state of the self-organized structure can be the one of a Mott insulator, depending on the pump strength of the driving laser. For very large pump strengths, where the intracavity-field intensity is maximum and one would expect a Mott-insulator state, we find intervals of parameters where the phase is compressible. These states could be realized in existing experimental setups.