A introduction to particle filtering is discussed starting with an overview of Bayesian inference from batch to sequential processors. Once the evolving Bayesian paradigm is established, simulation-based methods using sampling theory and Monte Carlo realizations are discussed. Here the usual limitations of nonlinear approximations and non-gaussian processes prevalent in classical nonlinear processing algorithms (e.g. Kalman filters) are no longer a restriction to perform Bayesian inference. It is shown how the underlying hidden or state variables are easily assimilated into this Bayesian construct. Importance sampling methods are then discussed and shown how they can be extended to sequential solutions implemented using Markovian state-space models as a natural evolution. With this in mind, the idea of a particle filter, which is a discrete representation of a probability distribution, is developed and shown how it can be implemented using sequential importance sampling/resampling methods. Finally, an application is briefly discussed comparing the performance of the particle filter designs with classical nonlinear filter implementations.

None

Direct imaging of exoplanets is limited by bright quasi-static speckles in the point spread function (PSF) of the central star. This limitation can be reduced by subtraction of reference PSF images. We have developed an algorithm to construct an optimal reference PSF image from an arbitrary set of reference images. This image is built as a linear combination of all available images and is optimized independently inside multiple subsections of the image to ensure that the absolute minimum residual noise is achieved within each subsection. The algorithm developed is completely general and can be used with many high contrast imaging observing strategies, such as angular differential imaging (ADI), roll subtraction, spectral differential imaging, reference star observations, etc. The performance of the algorithm is demonstrated for ADI data. It is shown that for this type of data the new algorithm provides a gain in sensitivity by up 22 to a factor 3 at small separation over the algorithm previously used.

The purpose of crucible scale testing with actual radioactive Tank 48H material was to duplicate the test results that had been previously performed on simulant Tank 48H material. The earlier crucible scale testing using simulants was successful in demonstrating that bench scale crucible tests produce results that are indicative of actual Fluidized Bed Steam Reforming (FBSR) pilot scale tests. Thus, comparison of the results using radioactive Tank 48H feed to those reported earlier with simulants would then provide proof that the radioactive tank waste behaves in a similar manner to the simulant. Demonstration of similar behavior for the actual radioactive Tank 48H slurry to the simulant is important as a preliminary or preparation step for the more complex bench-scale steam reformer unit that is planned for radioactive application in the Savannah River National Laboratory (SRNL) Shielded Cells Facility (SCF) later in 2008. The goals of this crucible-scale testing were to show 99% destruction of tetraphenylborate and to demonstrate that the final solid product produced is sodium carbonate. Testing protocol was repeated using the specifications of earlier simulant crucible scale testing, that is sealed high purity alumina crucibles containing a pre-carbonated and evaporated Tank 48H material. Sealing of the crucibles was …

The evolution of dense aluminum and carbon plasmas produced by laser irradiation of 500 {micro}m diameter semi-cylindrical targets was studied using soft x-ray laser interferometry. Plasmas created heating the cavity walls with 120 picosecond duration optical laser pulses of {approx} 1 x 10{sup 12} W cm{sup -2} peak intensity were observed to expand and converge on axis to form a localized high density plasma region. Electron density maps were measured using a 46.9 nm wavelength tabletop capillary discharge soft x-ray laser probe in combination with an amplitude division interferometer based on diffraction gratings. The measurements show that the plasma density on axis exceeds 1 x 10{sup 20} cm{sup -3}. The electron density profiles are compared with simulations conducted using the hydrodynamic code HYDRA, which show that the abrupt density increase near the axis is dominantly caused by the convergence of plasma generated at the bottom of the groove during laser irradiation.

Soft x-ray interferometry was used to measure the evolution of dense converging plasmas created by laser irradiation of 500 {micro}m diameter semi-cylindrical carbon targets. Optical laser pulses with an intensity of {approx} 1 x 10{sup 12} W cm{sup -2} and 120 ps duration were used to heat the surface of the cavities. The dense plasma formed expands from the walls converging slightly off the semi-cylinder's axis, giving rise to a bright localized high density plasma region. A sequence of electron density maps were measured using a 46.9 nm wavelength tabletop capillary discharge soft x-ray laser probe and a amplitude division interferometer based on diffraction gratings. The measured density profiles are compared with simulations conducted using the multi-dimensional hydrodynamic code HYDRA. The benchmarked model was then used to simulate particle trajectories which reveal that the increase in electron density near the axis is mainly the result of the convergence of plasma that originated at the bottom of the groove during laser irradiation.

The 100-D-24 Sample Building Drywell waste site was a drywell that received drainage from a floor drain in the 119-D Sample Building. Confirmatory sampling was conducted on November 3, 2005. The waste site meets the remedial action objectives specified in the Remaining Sites ROD. The results of confirmatory sampling show that residual contaminant concentrations do not preclude any future uses and allow for unrestricted use of shallow zone soils. The results also demonstrate that residual contaminant concentrations are protective of groundwater and the Columbia River.

Early detection of cancer is a key element in successful treatment of the disease. Understanding the particular type of cancer involved, its origins and probable course, is also important. PhIP (2-amino-1-methyl-6 phenylimidazo [4,5-b]pyridine), a heterocyclic amine produced during the cooking of meat at elevated temperatures, has been shown to induce mammary cancer in female, Sprague-Dawley rats. Tumors induced by PhIP have been shown to contain discreet cytogenetic signature patterns of gains and losses using comparative genomic hybridization (CGH). To determine if a protein signature exists for these tumors, we are analyzing expression levels of the protein products of the above-mentioned tumors in combination with a new bulk protein subtractive assay. This assay produces a panel of antibodies against proteins that are either on or off in the tumor. Hybridization of the antibody panel onto a 2-D gel of tumor or control protein will allow for identification of a distinct protein signature in the tumor. Analysis of several gene databases has identified a number of rat homologs of human cancer genes located in these regions of gain and loss. These genes include the oncogenes c-MYK, ERBB2/NEU, THRA and tumor suppressor genes EGR1 and HDAC3. The listed genes have been shown to …

Self-doping in graphene has been studied by examining single-layer epitaxially grown graphene samples with differing characteristic lateral terrace widths. Low energy electron microscopy was used to gain real-space information about the graphene surface morphology, which was compared with data obtained by angle-resolved photoemission spectroscopy to study the effect of the monolayer graphene terrace width on the low energy dispersions. By altering the graphene terrace width, we report significant changes in the electronic structure and quasiparticle relaxation time of the material, in addition to a terrace width-dependent doping effect.

The current conceptual model of radionuclide transport in unsaturated fractured rock includes water movement in fractures, with migration of the entrained radionuclides being retarded by diffusion into and sorption within the rock matrix. Water infiltration and radionuclide transport through low-permeability unsaturated fractured rock are episodic and intermittent in nature, at least at local scales. Under episodic flow conditions, the matrix is constantly imbibing or draining, and this fluctuating wetness both drives two-way advective movement of radionuclides, and forces changes in the matrix diffusivity. This work is intended to examine, both experimentally and numerically, how radionuclide transport under episodic flow conditions is affected by the interacting processes of imbibition and drainage, diffusion, and matrix sorption. Using Topopah Spring welded volcanic tuff, collected from the potential repository geologic unit at Yucca Mountain for storing high-level nuclear waste, we prepared a saw-cut fracture core (length 10.2 cm, diameter 4.4 cm, and fracture aperture 100 {micro}m). The dry core was packed into a flow reactor, flushed with CO{sub 2}, then saturated via slow pumping (0.01 mL/min) of synthetic groundwater. The fractured core was then flushed with air at >97% relative humidity (to simulate in situ unsaturated fractured rock conditions at Yucca Mountain), then the …

The objective of the Gypsy Project was to properly calculate seismic attributes and integrate these into a reservoir characterization project. Significant progress was made on the project in four areas. (1) Attenuation: In order for seismic inversion for rock properties or calculation of seismic attributes used to estimate rock properties to be performed validly, it is necessary to deal with seismic data that has had true amplitude and frequency content restored to account for earth filtering effects that are generally not included in seismic reservoir characterization methodologies. This requires the accurate measurement of seismic attenuation, something that is rarely achieved in practice. It is hoped that such measurements may also provide additional independent seismic attributes for use in reservoir characterization studies. In 2000, we were concerned with the ground truthing of attenuation measurements in the vicinity of wells. Our approach to the problem is one of extracting as time varying wavelet and relating temporal variations in the wavelet to an attenuation model of the earth. This method has the advantage of correcting for temporal variations in the reflectivity spectrum of the earth which confound the spectral ratio methodology which is the most commonly applied means of measuring attenuation from surface …

Time-dependent density-functional theory (TDDFT) provides an efficient, elegant, and formally exact way of describing the dynamics of interacting many-body quantum systems, circumventing the need for solving the full time-dependent Schroedinger equation. In the 20 years since it was first rigorously established in 1984, the field of TDDFT has made rapid and significant advances both formally as well as in terms of successful applications in chemistry, physics and materials science. Today, TDDFT has become the method of choice for calculating excitation energies of complex molecules, and is becoming increasingly popular for describing optical and spectroscopic properties of a variety of materials such as bulk solids, clusters and nanostructures. Other growing areas of applications of TDDFT are nonlinear dynamics of strongly excited electronic systems and molecular electronics. The purpose and scope of this Gordon Research Conference is to provide a platform for discussing the current state of the art of the rapidly progressing, highly interdisciplinary field of TDDFT, to identify and debate open questions, and to point out new promising research directions. The conference will bring together experts with a diverse background in chemistry, physics, and materials science.

Amine nitrogen inversion, difficult to observe in aqueous solution, is followed in a chiral, supramolecular host molecule with purely-rotational T-symmetry that reduces the local symmetry of encapsulated monoprotonated diamines and enables the observation and quantification of {Delta}G{double_dagger} for the combined hydrogen-bond breaking and nitrogen inversion rotation (NIR) process. Free energies of activation for the combined hydrogen-bond breaking and NIR process inside of the chiral assembly were determined by the NMR coalescence method. Activation parameters for ejection of the protonated amines from the assembly confirm that the NIR process responsible for the coalescence behavior occurs inside of the assembly rather than by a guest ejection/NIR/re-encapsulation mechanism. For one of the diamines, N,N,N{prime},N{prime}-tetramethylethylenediamine (TMEDA), the relative energy barriers for the hydrogen-bond breaking and NIR process were calculated at the G3(MP2)//B3LYP/6-31++G(d,p) level of theory, and these agreed well with the experimental data.

The electrochemical performance of LiFePO4/C composites in lithium cells is closely correlated to pressed pellet conductivities measured by AC impedance methods. These composite conductivities are a strong function not only of the amount of carbon but of its structure and distribution. Ideally, the amount of carbon in composites should be minimal (less than about 2 wtpercent) so as not to decrease the energy density unduly. This is particularly important for plug-in hybrid electric vehicle applications (PHEVs) where both high power and moderate energy density are required. Optimization of the carbon structure, particularly the sp2/sp3 and disordered/graphene (D/G) ratios, improves the electronic conductivity while minimizing the carbon amount. Manipulation of the carbon structure can be achieved via the use of synthetic additives including ironcontaining graphitization catalysts. Additionally, combustion synthesis techniques allow co-synthesis of LiFePO4 and carbon fibers or nanotubes, which can act as"nanowires" for the conduction of current during cell operation.

The definition of what makes an article interesting varies from user to user and continually evolves even for a single user. As a result, for news recommendation systems, useless document features can not be determined a priori and all features are usually considered for interestingness classification. Consequently, the presence of currently useless features degrades classification performance [1], particularly over the initial set of news articles being classified. The initial set of document is critical for a user when considering which particular news recommendation system to adopt. To address these problems, we introduce an improved version of the naive Bayes classifier with online feature selection. We use correlation to determine the utility of each feature and take advantage of the conditional independence assumption used by naive Bayes for online feature selection and classification. The augmented naive Bayes classifier performs 28% better than the traditional naive Bayes classifier in recommending news articles from the Yahoo! RSS feeds.

New corrosion-resistant, iron-based amorphous metals have been identified from published data or developed through combinatorial synthesis, and tested to determine their relative corrosion resistance. Many of these materials can be applied as coatings with advanced thermal spray technology. Two compositions have corrosion resistance superior to wrought nickel-based Alloy C-22 (UNS No. N06022) in some very aggressive environments, including concentrated calcium-chloride brines at elevated temperature. Two Fe-based amorphous metal formulations have been found that appear to have corrosion resistance comparable to, or better than that of Ni-based Alloy C-22, based on breakdown potential and corrosion rate. Both Cr and Mo provide corrosion resistance, B enables glass formation, and Y lowers critical cooling rate (CCR). SAM1651 has yttrium added, and has a nominal critical cooling rate of only 80 Kelvin per second, while SAM2X7 (similar to SAM2X5) has no yttrium, and a relatively high critical cooling rate of 610 Kelvin per second. Both amorphous metal formulations have strengths and weaknesses. SAM1651 (yttrium added) has a low critical cooling rate (CCR), which enables it to be rendered as a completely amorphous thermal spray coating. Unfortunately, it is relatively difficult to atomize, with powders being irregular in shape. This causes the powder to be …

The multi-institutional High Performance Corrosion Resistant Materials (HPCRM) Team is cosponsored by the Defense Advanced Projects Agency (DARPA) Defense Science Office (DSO) and the Department of Energy (DOE) Office of Civilian Radioactive Waste Management (OCRWM), and has developed new corrosion-resistant, iron-based amorphous metals that can be applied as coatings with advanced thermal spray technology. Two compositions have corrosion resistance superior to wrought nickel-based Alloy C-22 (UNS No. N06022) in very aggressive environments, including concentrated calcium-chloride brines at elevated temperature. Corrosion costs the Department of Defense billions of dollars every year, with an immense quantity of material in various structures undergoing corrosion. For example, in addition to fluid and seawater piping, ballast tanks, and propulsions systems, approximately 345 million square feet of structure aboard naval ships and crafts require costly corrosion control measures. The use of advanced corrosion-resistant materials to prevent the continuous degradation of this massive surface area would be extremely beneficial. The Fe-based corrosion-resistant, amorphous-metal coatings under development may prove of importance for applications on ships. Such coatings could be used as an 'integral drip shield' on spent fuel containers, as well as protective coatings that could be applied over welds, thereby preventing exposure to environments that might cause …

The topography, geology, hydrology, vegetation, and wildlife of the Oak Ridge Reservation (ORR) provide a complex and intricate array of resources that directly impact land stewardship and use decisions (Fig. 1). The purpose of this document is to consolidate general information regarding the natural resources and physical characteristics of the ORR. The ORR, encompassing 33,114 acres (13,401 ha) of federally owned land and three Department of Energy (DOE) installations, is located in Roane and Anderson Counties in east Tennessee, mostly within the corporate limits of the city of Oak Ridge and southwest of the population center of Oak Ridge. The ORR is bordered on the north and east by the population center of the city of Oak Ridge and on the south and west by the Clinch River/Melton Hill Lake impoundment. All areas of the ORR are relatively pristine when compared with the surrounding region, especially in the Valley and Ridge Physiographic Province (Fig. 2). From the air, the ORR is clearly a large and nearly continuous island of forest within a landscape that is fragmented by urban development and agriculture. Satellite imagery from 2006 was used to develop a land-use/land-cover cover map of the ORR and surrounding lands (Fig. …

A carbon credit market is developing in the United States. Information is needed by buyers and sellers of carbon credits so that the market functions equitably and efficiently. Analyses have been conducted to determine the optimal forest management regime to employ for each of the major commercial tree species so that profitability of timber production only or the combination of timber production and carbon sequestration is maximized. Because the potential of a forest ecosystem to sequester carbon depends on the tree species, site quality and management regimes utilized, analyses have determined how to optimize carbon sequestration by determining how to optimally manage each species, given a range of site qualities, discount rates, prices of carbon credits and other economic variables. The effects of a carbon credit market on the method and profitability of forest management, the cost of sequestering carbon, the amount of carbon that can be sequestered, and the amount of timber products produced has been determined.

The Hanford Seismic Assessment Program (HSAP) provides an uninterrupted collection of high-quality raw and processed seismic data from the Hanford Seismic Network for the U.S. Department of Energy and its contractors. The Hanford Seismic Assessment Team locates and identifies sources of seismic activity and monitors changes in the historical pattern of seismic activity at the Hanford Site. The data are compiled, archived, and published for use by the Hanford Site for waste management, Natural Phenomena Hazards assessments, and engineering design and construction. In addition, the seismic monitoring organization works with the Hanford Site Emergency Services Organization to provide assistance in the event of a significant earthquake on the Hanford Site. The Hanford Seismic Network and the Eastern Washington Regional Network consist of 41 individual sensor sites and 15 radio relay sites maintained by the Hanford Seismic Assessment Team. For the Hanford Seismic Network, 16 local earthquakes were recorded during the third quarter of fiscal year 2007. The largest event (magnitude 2.0) occurred on April 16, 2007 and was located 4 km southwest of the 400 Area in the Columbia River basalts at a depth of approximately 3 km. Stratigraphically, 7 earthquakes occurred in the Columbia River basalts (approximately 0-5 km …

In this work we present a NIDS cluster as a scalable solution for realizing high-performance, stateful network intrusion detection on commodity hardware. The design addresses three challenges: (i) distributing traffic evenly across an extensible set of analysis nodes in a fashion that minimizes the communication required for coordination, (ii) adapting the NIDS's operation to support coordinating its low-level analysis rather than just aggregating alerts; and (iii) validating that the cluster produces sound results. Prototypes of our NIDS cluster now operate at the Lawrence Berkeley National Laboratory and the University of California at Berkeley. In both environments the clusters greatly enhance the power of the network security monitoring.

The thesis of this extended abstract is simple. High productivity comes from high level infrastructures. To measure this, we introduce a methodology that goes beyond the tradition of timing software in serial and tuned parallel modes. We perform a classroom productivity study involving 29 students who have written a homework exercise in a low level language (MPI message passing) and a high level language (Star-P with MATLAB client). Our conclusions indicate what perhaps should be of little surprise: (1) the high level language is always far easier on the students than the low level language. (2) The early versions of the high level language perform inadequately compared to the tuned low level language, but later versions substantially catch up. Asymptotically, the analogy must hold that message passing is to high level language parallel programming as assembler is to high level environments such as MATLAB, Mathematica, Maple, or even Python. We follow the Kepner method that correctly realizes that traditional speedup numbers without some discussion of the human cost of reaching these numbers can fail to reflect the true human productivity cost of high performance computing. Traditional data compares low level message passing with serial computation. With the benefit of a …

The branching fractions for the decays B {yields} P{ell}{nu}{sub {ell}}, where P are the pseudoscalar charmless mesons {pi}{sup {+-}}, {pi}{sup 0}, {eta} and {eta}{prime} and {ell} is an electron or muon, are measured with B{sup 0} and B{sup {+-}} mesons found in the recoil of a second B meson decaying as B {yields} D{ell}{nu}{sub {ell}} or B {yields} D*{ell}{nu}{sub {ell}}. The measurements are based on a data set of 348 fb{sup -1} of e{sup +}e{sup -} collisions at {radical}s = 10.58 GeV recorded with the BABAR detector. Assuming isospin symmetry, measured pionic branching fractions are combined into {Beta}(B{sup 0} {yields} {pi}{sup -}{ell}{sup +}{nu}{sub {ell}}) = (1.54 {+-} 0.17{sub (stat)} {+-} 0.09{sub (syst)}) x 10{sup -4}. First evidence of the B{sup +} {yields} {eta}{ell}{sup +}{nu}{sub {ell}} decay is seen; its branching fraction is measured to be {Beta}(B{sup +} {yields} {eta}{ell}{sup +}{nu}{sub {ell}}) = (0.64 {+-} 0.20{sub (stat)} {+-} 0.03{sub (syst)}) x 10{sup -4}. It is determined that {Beta}(B{sup +} {yields} {eta}{prime}{ell}{sup +}{nu}{sub {ell}}) < 0.47 x 10{sup -4} to 90% confidence. Partial branching fractions for the pionic decays in ranges of the momentum transfer and various published calculations of the B {yields} {pi} hadronic form factor are used to obtain values …

The author assumes that the reader is familiar with the Spherical harmonics, Pn, method and the discrete ordinates, S{sub n}, method; for a derivation of the equations used in these methods. I will only discuss the Boltzmann equation in one dimension, and the Sn method using Gaussian quadrature. I will do this merely to simplify the following discussion; once you understand the concepts presented here you can easily extend the conclusions to more general situations. Why are the spherical harmonics P{sub n} and discrete ordinate S{sub n} methods, or more correctly the P{sub n} and S{sub n+1} methods, equivalent, e.g., P{sub 3} is equivalent to S{sub 4}? When the S{sub n} method uses a Gaussian quadrature most textbooks will tell you that both methods are equivalent to assuming that the angular flux can be represented by a Legendre polynomial expansion of order n. Most textbooks are wrong [1]. We know that the S{sub n} method constrains the ''particles'' to travel in discrete directions; when Gaussian quadrature is used these discrete directions correspond to the zeros of the Legendre polynomial P{sub n+1}({mu}). What is not immediately obvious is that the P{sub n} method constrains the ''particles'' in exactly the same way. …