The Diffusion Kernel Filter

A particle filter method is presented for the discrete-time filtering problem with nonlinear ItA ` stochastic ordinary differential equations (SODE) with additive noise supposed to be analytically integrable as a function of the underlying vector-Wiener process and time. The Diffusion Kernel Filter is arrived at by a parametrization of small noise-driven state fluctuations within branches of prediction and a local use of this parametrization in the Bootstrap Filter. The method applies for small noise and short prediction steps. With explicit numerical integrators, the operations count in the Diffusion Kernel Filter is shown to be smaller than in the Bootstrap Filter whenever the initial state for the prediction step has sufficiently few moments. The established parametrization is a dual-formula for the analysis of sensitivity to gaussian-initial perturbations and the analysis of sensitivity to noise-perturbations, in deterministic models, showing in particular how the stability of a deterministic dynamics is modeled by noise on short times and how the diffusion matrix of an SODE should be modeled (i.e. defined) for a gaussian-initial deterministic problem to be cast into an SODE problem. From it, a novel definition of prediction may be proposed that coincides with the deterministic path within the branch of prediction whose information entropy at the end of the prediction step is closest to the average information entropy over all branches. Tests are made with the Lorenz-63 equations, showing good results both for the filter and the definition of prediction.

Partitions selection strategy for set of clustering solutions

Clustering is a difficult task: there is no single cluster definition and the data can have more than one underlying structure. Pareto-based multi-objective genetic algorithms (e.g., MOCK Multi-Objective Clustering with automatic K-determination and MOCLE-Multi-Objective Clustering Ensemble) were proposed to tackle these problems. However, the output of such algorithms can often contains a high number of partitions, becoming difficult for an expert to manually analyze all of them. In order to deal with this problem, we present two selection strategies, which are based on the corrected Rand, to choose a subset of solutions. To test them, they are applied to the set of solutions produced by MOCK and MOCLE in the context of several datasets. The study was also extended to select a reduced set of partitions from the initial population of MOCLE. These analysis show that both versions of selection strategy proposed are very effective. They can significantly reduce the number of solutions and, at the same time, keep the quality and the diversity of the partitions in the original set of solutions. (C) 2010 Elsevier B.V. All rights reserved.

C(k)-solvability near the characteristic set for a class of planar complex vector fields of infinite type

This paper deals with semi-global C(k)-solvability of complex vector fields of the form L = partial derivative/partial derivative t + x(r) (a(x) + ib(x))partial derivative/partial derivative x, r >= 1, defined on Omega(epsilon) = (-epsilon, epsilon) x S(1), epsilon > 0, where a and b are C(infinity) real-valued functions in (-epsilon, epsilon). It is shown that the interplay between the order of vanishing of the functions a and b at x = 0 influences the C(k)-solvability at Sigma = {0} x S(1). When r = 1, it is permitted that the functions a and b of L depend on the x and t variables, that is, L = partial derivative/partial derivative t + x(a(x, t) + ib(x, t))partial derivative/partial derivative x, where (x, t) is an element of Omega(epsilon).

Gevrey solvability near the characteristic set for a class of planar complex vector fields of infinite type

We study the Gevrey solvability of a class of complex vector fields, defined on Omega(epsilon) = (-epsilon, epsilon) x S(1), given by L = partial derivative/partial derivative t + (a(x) + ib(x))partial derivative/partial derivative x, b not equivalent to 0, near the characteristic set Sigma = {0} x S(1). We show that the interplay between the order of vanishing of the functions a and b at x = 0 plays a role in the Gevrey solvability. (C) 2008 Elsevier Inc. All rights reserved.

A geometric mass-preserving redistancing scheme for the level set function

In this paper we describe and evaluate a geometric mass-preserving redistancing procedure for the level set function on general structured grids. The proposed algorithm is adapted from a recent finite element-based method and preserves the mass by means of a localized mass correction. A salient feature of the scheme is the absence of adjustable parameters. The algorithm is tested in two and three spatial dimensions and compared with the widely used partial differential equation (PDE)-based redistancing method using structured Cartesian grids. Through the use of quantitative error measures of interest in level set methods, we show that the overall performance of the proposed geometric procedure is better than PDE-based reinitialization schemes, since it is more robust with comparable accuracy. We also show that the algorithm is well-suited for the highly stretched curvilinear grids used in CFD simulations. Copyright (C) 2010 John Wiley & Sons, Ltd.

Partial spectral projected gradient method with active-set strategy for linearly constrained optimization

A method for linearly constrained optimization which modifies and generalizes recent box-constraint optimization algorithms is introduced. The new algorithm is based on a relaxed form of Spectral Projected Gradient iterations. Intercalated with these projected steps, internal iterations restricted to faces of the polytope are performed, which enhance the efficiency of the algorithm. Convergence proofs are given and numerical experiments are included and commented. Software supporting this paper is available through the Tango Project web page: http://www.ime.usp.br/similar to egbirgin/tango/.

A Floating-point Extended Kalman Filter Implementation for Autonomous Mobile Robots

Localization and Mapping are two of the most important capabilities for autonomous mobile robots and have been receiving considerable attention from the scientific computing community over the last 10 years. One of the most efficient methods to address these problems is based on the use of the Extended Kalman Filter (EKF). The EKF simultaneously estimates a model of the environment (map) and the position of the robot based on odometric and exteroceptive sensor information. As this algorithm demands a considerable amount of computation, it is usually executed on high end PCs coupled to the robot. In this work we present an FPGA-based architecture for the EKF algorithm that is capable of processing two-dimensional maps containing up to 1.8 k features at real time (14 Hz), a three-fold improvement over a Pentium M 1.6 GHz, and a 13-fold improvement over an ARM920T 200 MHz. The proposed architecture also consumes only 1.3% of the Pentium and 12.3% of the ARM energy per feature.

All-optical switching device for infrared based on PbTe quantum dots

Multilayers of PbTe quantum dots embedded in SiO2 were fabricated by alternate use of Pulsed Laser Deposition (PLD) and Plasma Enhanced Chemical Vapor Deposition (PECVD) techniques. The morphological properties of the nanostructured material were studied by means of High Resolution Transmission Electron Microscopy (HRTEM), Grazing-Incidence Small-Angle X-ray scattering (GISAXS) and X-ray Reflectometry (XRR) techniques. A preliminary analysis of the GISAXS spectra provided information about the multilayer periodicity and its relationship to the size of the deposited PbTe nanoparticles. Finally multilayers were fabricated inside a Fabry-Perot cavity. The device was characterized by means of Scanning Electron Microscopy (SEM). Transmittance measurements show the device functionality in the infrared region. (C) 2007 Elsevier Ltd. All rights reserved.

Aerosol inorganic composition at a tropical site: Discrepancies between filter-based sampling and a semi-continuous method

The concentrations of the water-soluble inorganic aerosol species, ammonium (NH4+), nitrate (NO3-), chloride (Cl-), and sulfate (SO42-), were measured from September to November 2002 at a pasture site in the Amazon Basin (Rondnia, Brazil) (LBA-SMOCC). Measurements were conducted using a semi-continuous technique (Wet-annular denuder/Steam-Jet Aerosol Collector: WAD/SJAC) and three integrating filter-based methods, namely (1) a denuder-filter pack (DFP: Teflon and impregnated Whatman filters), (2) a stacked-filter unit (SFU: polycarbonate filters), and (3) a High Volume dichotomous sampler (HiVol: quartz fiber filters). Measurements covered the late dry season (biomass burning), a transition period, and the onset of the wet season (clean conditions). Analyses of the particles collected on filters were performed using ion chromatography (IC) and Particle-Induced X-ray Emission spectrometry (PIXE). Season-dependent discrepancies were observed between the WAD/SJAC system and the filter-based samplers. During the dry season, when PM2.5 (D-p <= 2.5 mu m) concentrations were similar to 100 mu g m(-3), aerosol NH4+ and SO42- measured by the filter-based samplers were on average two times higher than those determined by the WAD/SJAC. Concentrations of aerosol NO3- and Cl- measured with the HiVol during daytime, and with the DFP during day- and nighttime also exceeded those of the WAD/SJAC by a factor of two. In contrast, aerosol NO3- and Cl- measured with the SFU during the dry season were nearly two times lower than those measured by the WAD/SJAC. These differences declined markedly during the transition period and towards the cleaner conditions during the onset of the wet season (PM2.5 similar to 5 mu g m(-3)); when filter-based samplers measured on average 40-90% less than the WAD/SJAC. The differences were not due to consistent systematic biases of the analytical techniques, but were apparently a result of prevailing environmental conditions and different sampling procedures. For the transition period and wet season, the significance of our results is reduced by a low number of data points. We argue that the observed differences are mainly attributable to (a) positive and negative filter sampling artifacts, (b) presence of organic compounds and organosulfates on filter substrates, and (c) a SJAC sampling efficiency of less than 100%.

Infrared chiral anomaly at finite temperature

We study the Schwinger model at finite temperature and show that a temperature dependent chiral anomaly may arise from the long distance behavior of the electric field. At high temperature this anomaly depends linearly on the temperature T and is present not only in the two point function, but also in all even point amplitudes. (C) 2011 Elsevier B.V. All rights reserved.

Infrared emissions, visible up-conversion, thermoluminescence and defect centres in Er(3)Al(5)O(12) phosphor obtained by solution combustion reaction

The Er(3)Al(5)O(12) phosphor powders were prepared using the solution combustion method. Formation and homogeneity of the Er(3)Al(5)O(12) phosphor powders have been verified by X-ray diffraction and energy-dispersive X-ray analysis respectively. The frequency up-conversion from Er(3)Al(5)O(12) phosphor powder corresponding to the (2)H(9/2) -> (4)I(15/2), (2)H(11/2) -> (4)I(15/2), (4)S(3/2) -> (4)I(15/2), (4)F(9/2) -> (4)I(15/2) and the infrared emission (IR) due to the (4)I(13/2) -> (4)I(15/2) transitions lying at similar to 410, similar to 524, similar to 556, 645-680 nm and at similar to 1.53 mu m respectively upon excitation with a Ti-Sapphire pulsed/CW laser have been reported. The mechanism responsible for the frequency up-conversion and IR emission is discussed in detail. Defect centres induced by radiation were studied using the techniques of thermoluminescence and electron spin resonance. A single glow peak at 430A degrees C is observed and the thermoluminescence results show the presence of a defect center which decays at high temperature. Electron spin resonance studies indicate a center characterized by a g-factor equal to 2.0056 and it is observed that this center is not related to the thermoluminescence peak. A negligibly small concentration of cation and anion vacancies appears to be present in the phosphor in accordance with the earlier theoretical predictions.

Infrared emission and defect centres in Er and Yb codoped Y(3)Al(5)O(12) phosphors

YAG phosphor powders doped/codoped with Er(3+)/(Er(3+) + Yb(3+)) have been synthesised by using the solution combustion method. The effect of direct pumping into the (4)I(11/2) level under 980 nm excitation of doped/codoped Er(3+)/Yb(3+)-Er(3+) in Y(3)Al(5)O(12) (YAG) phosphor responsible for an infrared (IR) emission peaking at similar to 1.53 mu m corresponding to the (4)I(13/2)->(4)I(15/2) transition has been studied. YAG exhibits three thermally-stimulated luminescence (TSL) peaks at around 140A degrees C, 210A degrees C and 445A degrees C. Electron spin resonance (ESR) studies were carried out to identify the centres responsible for the TSL peaks. The room temperature ESR spectrum of irradiated phosphor appears to be a superposition of two distinct centres. One of the centres (centre I) with principal g-value 2.0176 is identified as O(-) ion, while centre II with an isotropic g-factor 2.0020 is assigned to an F(+) centre (singly ionised oxygen vacancy). An additional defect centre is observed during thermal-annealing experiments and this centre (assigned to F(+) centre) seems to originate from an F-centre (oxygen vacancy with two electrons) and these two centres appear to correlate with the observed high-temperature TSL peak in YAG phosphor.

Infrared luminescence, thermoluminescence and defect centres in Er and Yb co-doped ZnAl(2)O(4) phosphor

Er and Yb co-doped ZnAl(2)O(4) phosphors were prepared by solution combustion synthesis and the identification of Er and Yb were done by energy-dispersive X-ray analysis (EDX) studies. A luminescence at 1.5 mu m, due to the (4)I(13/2) ->(4)I(15/2) transition, has been studied in the NIR region in Er and Yb co-doped ZnAl(2)O(4) phosphors upon 980 nm CW pumping. Er-doped ZnAl(2)O(4) exhibits two thermally stimulated luminescence (TSL) peaks around 174A degrees C and 483A degrees C, while Yb co-doped ZnAl(2)O(4) exhibits TSL peaks around 170A degrees C and 423A degrees C. Electron spin resonance (ESR) studies were carried out to identify defect centres responsible for TSL peaks observed in the phosphors. Room temperature ESR spectrum appears to be a superposition of two distinct centres. These centres are assigned to an O(-) ion and F(+) centre. O(-) ion appears to correlate with the 174A degrees C TSL peak and F(+) centre appears to relate with the high temperature TSL peak at 483A degrees C in ZnAl(2)O(4):Er phosphor.

NIR to visible up-conversion, infrared luminescence, thermoluminescence and defect centres in Y(2)O(3) : Er phosphor

Er(3+) doped Y(2)O(3) phosphor was prepared by the solution combustion method and characterized using powder x-ray diffraction and energy-dispersive analysis of x-ray mapping studies. Room temperature near infrared (NIR) to green up-conversion (UC) emissions in the region 520-580 nm {((2)H(11/2), (4)S(3/2)) -> (4)I(15/2)} and red UC emissions in the region 650-700 nm ((4)F(9/2) -> (4)I(15/2)) of Er(3+) ions have been observed upon direct excitation to the (4)I(11/2) level using similar to 972 nm laser radiation of nanosecond pulses. The possible mechanisms for the UC processes have been discussed on the basis of the energy level scheme, the pump power dependence as well as based on the temporal evolution. The excited state absorption is observed to be the dominant mechanism for the UC process. Y(2)O(3) : Er exhibits one thermally stimulated luminescence (TSL) peak around 367 degrees C. Electron spin resonance (ESR) studies were carried out to study the defect centres induced in the phosphor by gamma irradiation and also to identify the centres responsible for the TSL peak. Room temperature ESR spectrum of irradiated phosphor appears to be a superposition of at least three distinct centres. One of them (centre I) with principal g-values g(parallel to) = 2.0415 and g(perpendicular to) = 2.0056 is identified as O(2)(-) centre while centre II with an isotropic g-factor 2.0096 is assigned to an F(+)-centre (singly ionized oxygen vacancy). Centre III is also assigned to an F(+)-centre with a small g-factor anisotropy (g(parallel to) = 1.974 and g(perpendicular to) = 1.967). Additional defect centres are observed during thermal annealing experiments and one of them appearing around 330 degrees C grows with the annealing temperature. This centre (assigned to an F(+)-centre) seems to originate from an F-centre (oxygen vacancy with two electrons) and the F-centre appears to correlate with the observed TSL peak in Y2O3 : Er phosphor. The trap depth for this peak has been determined to be 0.97 eV from TSL data.

Thermoluminescence, infrared reflectivity and electron paramagnetic resonance properties of hemimorphite

Silicate mineral hemimorphite has been investigated concerning its TL, IR and EPR properties. A broad TL peak around 180 degrees C and a weaker and narrower peak around 360 degrees C were found in a sample annealed at 600 degrees C for I h and then irradiated. The deconvolution using the CGCD method revealed peaks around 132, 169, 222 and 367 degrees C. The reflectivity measurements showed several bands in the NIR region due to H(2)O, OH and Al-OH complexes. No band was observed in the visible region. The thermal treatments were carried out from similar to 110 to 940 degrees C and dehydration was observed, first causing a diminishing optical absorption in general and the disappearance of water and hydroxyl absorption bands. The EPR spectrum of natural hemimorphite, presented Cu(2+) signals at g = 2.4 and g = 2.1 plus E(1)` signal superposed to Fe(3+) signal around g = 2.0. (C) 2008 Elsevier Ltd. All rights reserved.