title: Quantum-state transfer from an ion to a photon

abstract: A quantum network requires information transfer between distant quantum computers, which would enable distributed quantum information processing and quantum communication. One model for such a network is based on the probabilistic measurement of two photons, each entangled with a distant atom or atomic ensemble, where the atoms represent quantum computing nodes. A second, deterministic model transfers information directly from a first atom onto a cavity photon, which carries it over an optical channel to a second atom; a prototype with neutral atoms has recently been demonstrated. In both cases, the central challenge is to find an efficient transfer process that preserves the coherence of the quantum state. Here, following the second scheme, we map the quantum state of a single ion onto a single photon within an optical cavity. Using an ion allows us to prepare the initial quantum state in a deterministic way, while the cavity enables high-efficiency photon generation. The mapping process is time-independent, allowing us to characterize the interplay between efficiency and fidelity. As the techniques for coherent manipulation and storage of multiple ions at a single quantum node are well established, this process offers a promising route toward networks between ion-based quantum computers.


title: Molecular Hydrogen and [Fe II] in Active Galactic Nuclei III: LINERS and Star Forming Galaxies

abstract: We study the kinematics and excitation mechanisms of H2 and [Fe II] lines in a sample of 67 emission-line galaxies with Infrared Telescope Facility SpeX near-infrared (NIR, 0.8-2.4 micrometers) spectroscopy together with new photoionisation models. H2 emission lines are systematically narrower than narrow-line region (NLR) lines, suggesting that the two are, very likely, kinematically disconnected. The new models and emission-line ratios show that the thermal excitation plays an important role not only in active galactic nuclei but also in star forming galaxies. The importance of the thermal excitation in star forming galaxies may be associated with the presence of supernova remnants close to the region emitting H2 lines. This hypothesis is further supported by the similarity between the vibrational and rotational temperatures of H2. We confirm that the diagram involving the line ratios H2 2.121/Br_gamma and [Fe II] 1.257/Pa_beta is an efficient tool for separating emission-line objects according to their dominant type of activity. We suggest new limits to the line ratios in order to discriminate between the different types of nuclear activity.


title: Strong bounds on Onsager coefficients and efficiency for three terminal thermoelectric transport in a magnetic field

abstract: For thermoelectric transport in the presence of a magnetic field that breaks time-reversal symmetry, a strong bound on the Onsager coefficients is derived within a general set-up using three terminals. Asymmetric Onsager coefficients lead to a maximum efficiency substantially smaller than the Carnot efficiency reaching only \eta_C/4 in the limit of strong asymmetry. Related bounds are derived for efficiency at maximum power, which can become larger than the Curzon-Ahlborn value \eta_C/2, and for a cooling device. Our approach reveals that in the presence of reversible currents the standard analysis based on the positivity of entropy production is incomplete without considering the role of current conservation explicitly.


title: NLO QCD corrections to the production of Higgs plus two jets at the LHC

abstract: We present the calculation of the NLO QCD corrections to the associated production of a Higgs boson and two jets, in the infinite top-mass limit. We discuss the technical details of the computation and we show the numerical impact of the radiative corrections on several observables at the LHC. The results are obtained by using a fully automated framework for fixed order NLO QCD calculations based on the interplay of the packages GoSam and Sherpa. The evaluation of the virtual corrections constitutes an application of the d-dimensional integrand-level reduction to theories with higher dimensional operators. We also present first results for the one-loop matrix elements of the partonic processes with a quark-pair in the final state, which enter the hadronic production of a Higgs boson together with three jets in the infinite top-mass approximation.


title: Towards quantum gravity measurement by cold atoms

abstract: We propose an experiment for the measurement of gravitational effect on cold atoms by applying a one-dimensional vertically sinusoidal oscillation to the magneto-optical trap; and observe the signature of low quantum energy shift of quantum bound states as a consequence of gravitational fluctuation. To this end, we present brief details of the experiment on a BEC, and a simplistic calculation of the Gross-Pitaevskii solution using Thomas-Fermi approximation with focus on the density of the BEC, for the time-dependent perturbation. Moreover, we calculate the power induced by quantum gravity on a generic atomic ensemble. We also address the possible challenges for the measurement of the expected results. And finally, we discuss the prospect of further developing this experiment towards measuring the effect of quantum spacetime fluctuations on cold atoms.


title: Discrete symmetries, roots of unity, and lepton mixing

abstract: We investigate the possibility that the first column of the lepton mixing matrix U is given by u_1 = (2,-1,-1)^T/sqrt{6}. In a purely group-theoretical approach, based on residual symmetries in the charged-lepton and neutrino sectors and on a theorem on vanishing sums of roots of unity, we discuss the finite groups which can enforce this. Assuming that there is only one residual symmetry in the Majorana neutrino mass matrix, we find the almost unique solution Z_q x S_4 where the cyclic factor Z_q with q = 1,2,3,... is irrelevant for obtaining u_1 in U. Our discussion also provides a natural mechanism for achieving this goal. Finally, barring vacuum alignment, we realize this mechanism in a class of renormalizable models.


title: Circumcenter of Mass and generalized Euler line

abstract: We define and study a variant of the center of mass of a polygon and, more generally, of a simplicial polytope which we call the Circumcenter of Mass (CCM). The Circumcenter of Mass is an affine combination of the circumcenters of the simplices in a triangulation of a polytope, weighted by their volumes. For an inscribed polytope, CCM coincides with the circumcenter. Our motivation comes from the study of completely integrable discrete dynamical systems, where the CCM is an invariant of the discrete bicycle (Darboux) transformation and of recuttings of polygons. We show that the CCM satisfies an analog of Archimedes' Lemma, a familiar property of the center of mass. We define and study a generalized Euler line associated to any simplicial polytope, extending the previously studied Euler line associated to the quadrilateral. We show that the generalized Euler line for polygons consists of all centers satisfying natural continuity and homogeneity assumptions and Archimedes' Lemma. Finally, we show that CCM can also be defined in the spherical and hyperbolic settings.


title: Parahoric induction and chamber homology for SL2

abstract: We consider the special linear group G=SL2 over a p-adic field, and its diagonal subgroup M=GL1. Parabolic induction of representations from M to G induces a map in equivariant homology, from the Bruhat-Tits building of M to that of G. We compute this map at the level of chain complexes, and show that it is given by parahoric induction (as defined by J.-F. Dat).


title: Hierarchical quantum communication

abstract: A general approach to study the hierarchical quantum information splitting (HQIS) is proposed and the same is used to systematically investigate the possibility of realizing HQIS using different classes of 4-qubit entangled states that are not connected by SLOCC. Explicit examples of HQIS using 4-qubit cluster state and 4-qubit |\Omega> state are provided. Further, the proposed HQIS scheme is generalized to introduce two new aspects of hierarchical quantum communication. To be precise, schemes of probabilistic hierarchical quantum information splitting and hierarchical quantum secret sharing are obtained by modifying the proposed HQIS scheme. A number of practical situations where hierarchical quantum communication would be of use are also presented.


title: Scalable time reversal of Raman echo quantum memory and quantum waveform conversion of light pulse

abstract: We have found a new hidden symmetry of time reversal light-atom interaction in the photon echo quantum memory with Raman atomic transition. The time-reversed quantum memory creates generalized conditions for ideal compression/decompression of time duration of the input light pulses and its wavelength. Based on a general analytical approach to this scheme, we have studied the optimal conditions for the light field compression/decompression in resonant atomic systems characterized by realistic spectral properties. The demonstrated necessary conditions for the effective quantum conversion of the light waveform and wavelength are also discussed for various possible realizations of the quantum memory scheme. The performed study promises new capabilities for fundamental study of the light-atom interaction and deterministic quantum manipulation of the light field, significant for quantum communication and quantum computing.


title: A spectroscopic investigation of the O-type star population in four Cygnus OB associations. I. Determination of the binary fraction

abstract: Context. Establishing the multiplicity of O-type stars is the first step towards accurately determining their stellar parameters. Moreover, the distribution of the orbital parameters provides observational clues to the way that O-type stars form and to the interactions during their evolution. Aims. Our objective is to constrain the multiplicity of a sample of O-type stars belonging to poorly investigated OB associations in the Cygnus complex and for the first time to provide orbital parameters for binaries identified in our sample. Such information is relevant to addressing the issue of the binarity in the context of O-type star formation scenarios. Methods. We performed a long-term spectroscopic survey of nineteen O-type stars. We searched for radial velocity variations to unveil binaries on timescales from a few days up to a few years, on the basis of a large set of optical spectra. Results. We confirm the binarity for four objects: HD193443, HD228989, HD229234 and HD194649. We derive for the first time the orbital solutions of three systems, and we confirm the values of the fourth, showing that these four systems all have orbital periods shorter than 10 days. Besides these results, we also detect several objects that show non-periodic line profile variations in some of their spectral lines. These variations mainly occur in the spectral lines, that are generally affected by the stellar wind and are not likely to be related to binarity. Conclusions. The minimal binary fraction in our sample is estimated to be 21%, but it varies from one OB association to the next. Indeed, 33% of O stars belonging to CygOB1 are binary systems, 0% in CygOB3, 0% in CygOB8, and 33% in CygOB9. Our spectroscopic investigation also stresses the absence of long-period systems among the stars in our sample...


title: The H\"older continuity of spectral measures of an extended CMV matrix

abstract: We prove results about the H\"older continuity of the spectral measures of the extended CMV matrix, given power law bounds of the solution of the eigenvalue equation. We thus arrive at a unitary analogue of the results of Damanik, Killip and Lenz about the spectral measure of the discrete Schr\"odinger operator.


title: Hamiltonian Dynamics for an alternative action describing Maxwell's equations

abstract: We develop a complete Dirac's canonical analysis for an alternative action that yields Maxwell's four-dimensional equations of motion. We study in detail the full symmetries of the action by following all steps of Dirac's method in order to obtain a detailed description of symmetries. Our results indicate that such an action does not have the same symmetries than Maxwell theory, namely, the model is not a gauge theory and the number of physical degrees of freedom are different.


title: Word Storms: Multiples of Word Clouds for Visual Comparison of Documents

abstract: Word clouds are a popular tool for visualizing documents, but they are not a good tool for comparing documents, because identical words are not presented consistently across different clouds. We introduce the concept of word storms, a visualization tool for analysing corpora of documents. A word storm is a group of word clouds, in which each cloud represents a single document, juxtaposed to allow the viewer to compare and contrast the documents. We present a novel algorithm that creates a coordinated word storm, in which words that appear in multiple documents are placed in the same location, using the same color and orientation, in all of the corresponding clouds. In this way, similar documents are represented by similar-looking word clouds, making them easier to compare and contrast visually. We evaluate the algorithm in two ways: first, an automatic evaluation based on document classification; and second, a user study. The results confirm that unlike standard word clouds, a coordinated word storm better allows for visual comparison of documents.


title: MuLan Measurement of the Positive Muon Lifetime and Determination of the Fermi Constant

abstract: We report results from the MuLan measurement of the positive muon lifetime. The experiment was conducted at the Paul Scherrer Institute using a time-structured surface muon beam and a segmented plastic scintillator array. Two different in-vacuum muon stopping targets were used: a ferromagnetic foil with a large internal magnetic field and a quartz crystal in a moderate external magnetic field. From a total of 1.6 x 10^{12} decays, we obtained the muon lifetime tau_mu = 2196980.3(2.2) ps (1.0 ppm) and Fermi constant G_F = 1.1663787(6) x 10^{-5} GeV^{-2} (0.5 ppm).


title: Distributions, their primitives and integrals with applications to differential equations

abstract: In this paper we will study integrability of distributions whose primitives are left regulated functions and locally or globally integrable in the Henstock--Kurzweil, Lebesgue or Riemann sense. Corresponding spaces of distributions and their primitives are defined and their properties are studied. Basic properties of primitive integrals are derived and applications to systems of first order nonlinear distributional differential equations and to an $m$th order distributional differential equation are presented. The domain of solutions can be unbounded, as shown by concrete examples.


title: Large Radio Telescopes for Anomalous Microwave Emission Observations

abstract: We discuss in this paper the problem of the Anomalous Microwave Emission (AME) in the light of ongoing or future observations to be performed with the largest fully steerable radio telescope in the world. High angular resolution observations of the AME will enable astronomers to drastically improve the knowledge of the AME mechanisms as well as the interplay between the different constituents of the interstellar medium in our galaxy. Extragalactic observations of the AME have started as well, and high resolution is even more important in this kind of observations. When cross-correlating with IR-dust emission, high angular resolution is also of fundamental importance in order to obtain unbiased results. The choice of the observational frequency is also of key importance in continuum observation. We calculate a merit function that accounts for the signal-to-noise ratio (SNR) in AME observation given the current state-of-the-art knowledge and technology. We also include in our merit functions the frequency dependence in the case of multifrequency observations. We briefly mention and compare the performance of four of the largest radiotelescopes in the world and hope the observational programs in each of them will be as intense as possible.


title: The interplay between the electromagnetic and the muonic longitudinal profile at production

abstract: The electromagnetic and the muonic longitudinal profile at production enclosure important information about the primary particle and the hadronic interactions that rule the shower development. In fact, these two profiles provide two different insights of the shower: the electromagnetic component gives a measurement of the energy and the strength of the neutral pion channel; while the muonic profile, being intimately related with the charged mesons decays, can be used as a direct probe for the high energy hadronic interactions. In this work we explore the interplay between the electromagnetic and muonic profiles, by analysing their phenomenologic behaviour for different primary masses and energies, zenith angles, and also different high energy hadronic interaction models. We have found that the muonic longitudinal profile at production displays universal features similar to what is known for the electromagnetic one. Moreover, we show that both profiles have new primary mass composition variables which are fairly independent of the high energy hadronic interaction model. Finally we discuss how the information in the electromagnetic and the muonic longitudinal profile can be combined to break the degeneracy between the primary mass composition and the high energy hadronic physics.


title: Time and a Temporally Statistical Quantum Geometrodynamics

abstract: This paper is an exposition of the author's recent work (arXiv:1001.3382 [hep-th], arXiv:1005.5430 [cond-mat.dis-nn], arXiv:1212.4956 [quant-ph]) on modeling M-theory vacua and quantum mechanical observers in the framework of a temporally statistical description of quantum geometrodynamics, including measurement processes based on the canonical theory of quantum gravity. In this paper we deal with several fundamental issues of time: the time-less problem in canonical quantum gravity; the physical origin of state reductions; and time-reversal symmetry breaking. We first model the observers and consider the time-less problem by invoking the time reparametrization symmetry breaking in the quantum mechanical world as seen by the observers. We next construct the hidden time variable theory, using a model of the gauged and affinized S-duality symmetry in type IIB string theory, as the statistical theory of time and explain the physical origin of state reductions using it. Finally, by the extension of the time reparametrization symmetry to all of the temporal hidden variables, we treat the issue of time reversal symmetry breaking as the spontaneous breaking of this extended time reparametrization symmetry. The classification of unitary time-dependent processes and the geometrizations of unitary and non-unitary time evolutions using the language of the derived category are also investigated.


title: Crystal Structure of BaFe2Se3 as a Function of Temperature and Pressure: Phase Transition Phenomena and High-Order Expansion of Landau Potential

abstract: BaFe2Se3 (Pnma, CsAg2I3-type structure), recently assumed to show superconductivity at ~ 11 K, exhibits a pressure-dependent structural transition to the CsCu2Cl3-type structure (Cmcm space group) around 60 kbar, as evidenced from pressure-dependent synchrotron powder diffraction data. Temperature-dependent synchrotron powder diffraction data indicate an evolution of the room-temperature BaFe2Se3 structure towards a high symmetry CsCu2Cl3 form upon heating. Around 425 K BaFe2Se3 undergoes a reversible, first order isostructural transition, that is supported by the differential scanning calorimetry data. The temperature-dependent structural changes occur in two stages, as determined by the alignment of the FeSe4 tetrahedra and corresponding adjustments of the positions of Ba atoms. On further heating, a second order phase transformation into the Cmcm structure is observed at 660 K. A rather unusual combination of isostructural and second-order phase transformations is parameterized within phenomenological theory assuming high-order expansion of Landau potential. A generic phase diagram mapping observed structures is proposed on the basis of the parameterization.


title: Partial classification of the Baumslag-Solitar group von Neumann algebras

abstract: We prove that the rational number |n/m| is an invariant of the group von Neumann algebra of the Baumslag-Solitar group BS(n,m). More precisely, if L(BS(n,m)) is isomorphic with L(\BS(n',m')), then |n'/m'| = |n/m| or |m/n|. We obtain this result by associating to abelian, but not maximal abelian, subalgebras of a II_1 factor, an equivalence relation that can be of type III. In particular, we associate to L(BS(n,m)) a canonical equivalence relation of type III_|n/m|.


title: Plancherel representations of $U(\infty)$ and correlated Gaussian Free Fields

abstract: We study asymptotics of traces of (noncommutative) monomials formed by images of certain elements of the universal enveloping algebra of the infinite-dimensional unitary group in its Plancherel representations. We prove that they converge to (commutative) moments of a Gaussian process that can be viewed as a collection of simply yet nontrivially correlated two-dimensional Gaussian Free Fields. The limiting process has previously arisen via the global scaling limit of spectra for submatrices of Wigner Hermitian random matrices. The results of the present work were announced in arXiv:1203.3010.


title: Giant Magnetized Outflows from the Centre of the Milky Way

abstract: The nucleus of the Milky Way is known to harbour regions of intense star formation activity as well as a super-massive black hole. Recent Fermi space telescope observations have revealed regions of \gamma-ray emission reaching far above and below the Galactic Centre, the so-called Fermi bubbles. It is uncertain whether these were generated by nuclear star formation or by quasar-like outbursts of the central black hole and no information on the structures' magnetic field has been reported. Here we report on the detection of two giant, linearly-polarized radio Lobes, containing three ridge-like sub-structures, emanating from the Galactic Centre. The Lobes each extend ~60 deg, bear a close correspondence to the Fermi bubbles, are located in the Galactic bulge, and are permeated by strong magnetic fields of up to 15 \mu G. Our data signal that the radio Lobes originate in a bi-conical, star-formation (rather than black hole) driven outflow from the Galaxy's central 200 pc that transports a massive magnetic energy of ~10^55 erg into the Galactic halo. The ridges wind around this outflow and, we suggest, constitute a `phonographic' record of nuclear star formation activity over at least 10 Myr.


title: Minimal Bounds on Nonlinearity in Auditory Processing

abstract: Time-reversal symmetry breaking is a key feature of nearly all natural sounds, caused by the physics of sound production. While attention has been paid to the response of the auditory system to "natural stimuli," very few psychophysical tests have been performed. We conduct psychophysical measurements of time-frequency acuity for both "natural" notes (sharp attack, long decay) and time-reversed ones. Our results demonstrate significantly greater precision, arising from enhanced temporal acuity, for such "natural" sounds over both their time-reversed versions and theoretically optimal gaussian pulses, without a corresponding decrease in frequency acuity. These data rule out models of auditory processing that obey a modified "uncertainty principle" between temporal and frequency acuity and suggest the existence of statistical priors for naturalistic stimuli, in the form of sharp-attack, long-decay notes. We are additionally able to calculate a minimal theoretical bound on the order of the nonlinearity present in auditory processing. We find that only matching pursuit, spectral derivatives, and reassigned spectrograms are able to satisfy this criterion.


title: Simulation of optical interstellar scintillation

abstract: Stars twinkle because their light propagates through the atmosphere. The same phenomenon is expected on a longer time scale when the light of remote stars crosses an interstellar turbulent molecular cloud, but it has never been observed at optical wavelengths. The aim of the study described in this paper is to fully simulate the scintillation process, starting from the molecular cloud description as a fractal object, ending with the simulations of fluctuating stellar light curves. Fast Fourier transforms are first used to simulate fractal clouds. Then, the illumination pattern resulting from the crossing of background star light through these refractive clouds is calculated from a Fresnel integral that also uses fast Fourier transform techniques. Regularisation procedure and computing limitations are discussed, along with the effect of spatial and temporal coherency (source size and wavelength passband). We quantify the expected modulation index of stellar light curves as a function of the turbulence strength --characterised by the diffraction radius $R_{diff}$-- and the projected source size, introduce the timing aspects, and establish connections between the light curve observables and the refractive cloud. We extend our discussion to clouds with different structure functions from Kolmogorov-type turbulence. Our study confirms that current telescopes of ~4m with fast-readout, wide-field detectors have the capability of discovering the first interstellar optical scintillation effects. We also show that this effect should be unambiguously distinguished from any other type of variability through the observation of desynchronised light curves, simultaneously measured by two distant telescopes.


title: Charged particle detection performances of CMOS pixel sensors produced in a 0.18 um process with a high resistivity epitaxial layer

abstract: The apparatus of the ALICE experiment at CERN will be upgraded in 2017/18 during the second long shutdown of the LHC (LS2). A major motivation for this upgrade is to extend the physics reach for charmed and beauty particles down to low transverse momenta. This requires a substantial improvement of the spatial resolution and the data rate capability of the ALICE Inner Tracking System (ITS). To achieve this goal, the new ITS will be equipped with 50 um thin CMOS Pixel Sensors (CPS) covering either the 3 innermost layers or all the 7 layers of the detector. The CPS being developed for the ITS upgrade at IPHC (Strasbourg) is derived from the MIMOSA 28 sensor realised for the STAR-PXL at RHIC in a 0.35 um CMOS process. In order to satisfy the ITS upgrade requirements in terms of readout speed and radiation tolerance, a CMOS process with a reduced feature size and a high resistivity epitaxial layer should be exploited. In this respect, the charged particle detection performance and radiation hardness of the TowerJazz 0.18 um CMOS process were studied with the help of the first prototype chip MIMOSA 32. The beam tests performed with negative pions of 120 GeV/c at the CERN-SPS allowed to measure a signal-to-noise ratio (SNR) for the non-irradiated chip in the range between 22 and 32 depending on the pixel design. The chip irradiated with the combined dose of 1 MRad and 10^13 n_eq/cm^2 was observed to yield a SNR ranging between 11 and 23 for coolant temperatures varying from 15 C to 30 C. These SNR values were measured to result in particle detection efficiencies above 99.5% and 98% before and after irradiation respectively. These satisfactory results allow to validate the TowerJazz 0.18 um CMOS process for the ALICE ITS upgrade.


title: Hochschild cohomology of triangular string algebras and its ring structure

abstract: We compute the Hochschild cohomology groups $\HH^*(A)$ in case $A$ is a triangular string algebra, and show that its ring structure is trivial.


title: Some Interesting Integer Polynomial Maps

abstract: We introduce three simple polynomial maps with integer coefficients that have interesting dynamical properties modulo primes.


title: Abelian Projections of the Mass-deformed ABJM theory and Weakly Curved Dual Geometry

abstract: We construct ${\cal N}=2,4$ supersymmetric abelian projections of the ${\cal N}=6$ mass-deformed ABJM theory. There are well-defined dual background geometries for the ${\cal N}=2$ abelian theory, while those geometries are unclear for the ${\cal N}=4$ abelian theory. The ${\cal N}=2$ theory is built on the supersymmetric vacua of the mass-deformed ABJM theory, which are proven to have one-to-one correspondence with the ${\mathbb Z}_k$ quotient of Lin-Lunin-Maldacena geometries. We select one special vacuum of the mass-deformed ABJM theory and show that the corresponding geometry is weakly curved at every point of the entire space transverse to the M2-branes in the large $N$ limit.


title: A pure Dirac's method for Yang-Mills expressed as a constrained BF-like theory

abstract: A pure Dirac's method of Yang-Mills expressed as a constrained BF-like theory is performed. In this paper we study an action principle composed by the coupling of two topological BF-like theories, which at the Lagrangian level reproduces Yang-Mills equations. By a pure Dirac's method we mean that we consider all the variables that occur in the Lagrangian density as dynamical variables and not only those ones that involve temporal derivatives. The analysis in the complete phase space enable us to calculate the extended Hamiltonian, the extended action, the constraint algebra, the gauge transformations and then we carry out the counting of degrees of freedom. We show that the constrained BF-like theory correspond at classical level to Yang-Mills theory. From the results obtained, we discuss briefly the quantization of the theory. In addition we compare our results with alternatives models that have been reported in the literature


title: Low-Ionization Outflows in High Eddington Ratio Quasars

abstract: The broad MgII doublet has been much studied in connection with its potentially important role as a virial estimator of black hole mass in high redshift quasars. An important task is therefore identification of any line components likely related to broadening by non-virial motions. High s/n median composite spectra (binned in the "4D eigenvector 1" context of Sulentic et al. 2007) were constructed for the brightest 680 SDSS DR7 quasars in the 0.4 < z < 0.75 range where both MgII 2800 and Hbeta are recorded in the same spectra. Composite spectra representing 90% of the quasars confirm previous findings that FWHM(MgII 2800) is about 20% narrower than FWHM(Hbeta). The situation is clearly different for the most extreme (Population A) sources which are the highest Eddington radiators in the sample. In the median spectra of these sources FWHM MgII 2800 is equal to or greater than FWHM(Hbeta) and shows a significant blueshift relative to Hbeta. We interpret the MgII 2800 blueshift as the signature of a radiation-driven wind or outflow in the highest accreting quasars. In this interpretation the MgII 2800 line width -- affected by blueshifted emission -- is unsuitable for virial mass estimation in ~ 10% of quasars.


title: Pair Breaking Caused by Magnetic Impurities in the High-T$_\text{C}$ Superconductor Bi$_{2.1}$Sr$_{1.9}$Ca(Cu$_{1-x}$Fe$_{x}$)$_{2}$O$_{y}$

abstract: Conventional superconductivity is robust against the addition of impurities unless the impurities are magnetic in which case superconductivity is quickly suppressed. Here we present a study of the cuprate superconductor Bi$_2$Sr$_2$Ca$_1$Cu$_2$O$_{8+\delta}$ that is intentionally doped with the magnetic impurity, Fe. Through the use of our Tomographic Density of States (TDoS) technique, we find that while the superconducting gap magnitude is essentially unaffected by the inclusion of iron, the onset of superconductivity, T$_{C}$, and the pair-breaking rate are strongly dependent and correlated. These findings suggest that, in the cuprates, the pair-breaking rate is critical to the determination of T$_{C}$ and that magnetic impurities do not disrupt the strength of pairing but rather the lifetime of the pairs.


title: A cold-atom random laser

abstract: Conventional lasers make use of optical cavities to provide feedback to gain media. Conversely, mirrorless lasers can be built by using disordered structures to induce multiple scattering, which increases the effective path length in the gain medium and thus provides the necessary feedback. These so-called random lasers potentially offer a new and simple mean to address applications such as lighting. To date, they are all based on condensed-matter media. Interestingly, light or microwave amplification by stimulated emission occurs also naturally in stellar gases and planetary atmospheres. The possibility of additional scattering-induced feedback (that is, random lasing) has been discussed and could explain unusual properties of some space masers. Here, we report the experimental observation of random lasing in a controlled, cold atomic vapour, taking advantage of Raman gain. By tuning the gain frequency in the vicinity of a scattering resonance, we observe an enhancement of the light emission of the cloud due to random lasing. The unique possibility to both control the experimental parameters and to model the microscopic response of our system provides an ideal test bench for better understanding natural lasing sources, in particular the role of resonant scattering feedback in astrophysical lasers.


title: EELA Operations: A standalone regional dashboard implementation

abstract: Grid operators in EGEE use a dedicated dashboard as their central operational tool, stable and scalable for the last 5 years despite continuous upgrade from specifications by users, monitoring tools or data providers. In EGEE-III, regionalisation of operations led the tool developers to conceive a standalone instance of this tool. Hereby, we will present the concept and the EELA-II implementation. Indeed, there-engineering of this tool led to an easily deployable package that canconnect to EELA-II specific information sources such as EVENTUM, EELAGOCDB like or SAM EELA instance through the three components of thepackage: the generic and scalable data access mechanism, Lavoisier; thewidely spread php framework Symfony, for configuration flexibility and a Mysql database.


title: Oxygen and light sensitive field-effect transistors based on ZnO nanoparticles attached to individual double-wall carbon nanotubes

abstract: The attachment of semiconducting nanoparticles to carbon nanotubes is one of the most challenging subjects in nanotechnology. Successful high coverage attachment and control over the charge transfer mechanism and photo-current generation opens a wide field of new applications such as highly effective solar cells and fibre-enhanced polymers. In this work we study the charge transfer in individual double-wall carbon nanotubes highly covered with uniform ZnO nanoparticles. The synthetic colloidal procedure was chosen to avoid long-chained ligands at the nanoparticle-nanotube interface. The resulting composite material was used as conductive channel in a field effect transistor device and the electrical photo-response was analysed under various conditions. By means of the transfer characteristics we could elucidate the mechanism of charge transfer from non-covalently attached semiconducting nanoparticles to carbon nanotubes. The role of positive charges remaining on the nanoparticles is discussed in terms of a gating effect.


title: FRW in cosmological self-creation theory: Hamiltonian approach

abstract: We use the Brans-Dicke theory from the framework of General Relativity (Einstein frame), but now the total energy momentum tensor fulfills the following condition $\rm[\frac{1}{\phi}(8\pi T^{\mu \nu (M)}+T^{\mu\nu (\phi)})]_{;\nu}=0$. We take as a first model the flat FRW metric in the Hamilton-Jacobi scheme and we present the Lagrange-Charpit approach in order to find classical solutions. In the quantum scheme, once we determine the characteristic surfaces, the quantum solution is obtained. These two classes of solutions are found for all values of the barotropic parameter $\gamma$.


title: Tensor product weight modules over the Virasoro algebra

abstract: The tensor product of highest weight modules with intermediate series modules over the Virasoro algebra was discussed by Zhang [Z] in 1997. Since then the irreducibility problem for the tensor products has been open. In this paper, we determine the necessary and sufficient conditions for these tensor products to be simple. From non-simple tensor products, we can get other interesting simple Virasoro modules. We also obtain that any two such tensor products are isomorphic if and only if the corresponding highest weight modules and intermediate series modules are isomorphic respectively. Our method is to develop a "shifting technique" and to widely use Feigin-Fuchs' Theorem on singular vectors of Verma modules over the Virasoro algebra.


title: Wicking in a powder

abstract: We investigate the wicking in granular media by considering layers of grains at the surface of a liquid, and discuss the critical contact angle below which spontaneous impregnation takes place. This angle is found to be on the order of 55{\deg} for monodisperse layers, significantly smaller than 90{\deg}, the threshold value for penetrating assemblies of tubes: owing to geometry, impregnating grains is more demanding than impregnating tubes. We also consider the additional effects of polydispersity and pressure on this wetting transition and discuss the corresponding shift observed on the critical contact angle.


title: Adaptive Electricity Scheduling in Microgrids

abstract: Microgrid (MG) is a promising component for future smart grid (SG) deployment. The balance of supply and demand of electric energy is one of the most important requirements of MG management. In this paper, we present a novel framework for smart energy management based on the concept of quality-of-service in electricity (QoSE). Specifically, the resident electricity demand is classified into basic usage and quality usage. The basic usage is always guaranteed by the MG, while the quality usage is controlled based on the MG state. The microgrid control center (MGCC) aims to minimize the MG operation cost and maintain the outage probability of quality usage, i.e., QoSE, below a target value, by scheduling electricity among renewable energy resources, energy storage systems, and macrogrid. The problem is formulated as a constrained stochastic programming problem. The Lyapunov optimization technique is then applied to derive an adaptive electricity scheduling algorithm by introducing the QoSE virtual queues and energy storage virtual queues. The proposed algorithm is an online algorithm since it does not require any statistics and future knowledge of the electricity supply, demand and price processes. We derive several "hard" performance bounds for the proposed algorithm, and evaluate its performance with trace-driven simulations. The simulation results demonstrate the efficacy of the proposed electricity scheduling algorithm.


title: Log-integrability of Rademacher Fourier series, with applications to random analytic functions

abstract: We prove that any power of the logarithm of Fourier series with random signs is integrable. This result has applications to the distribution of values of random Taylor series, one of which answers a long-standing question by J.-P. Kahane.


title: The Origin and Optical Depth of Ionizing Radiation in the "Green Pea" Galaxies

abstract: Although Lyman continuum (LyC) radiation from star-forming galaxies likely drove the reionization of the Universe, observations of star-forming galaxies at low redshift generally indicate low LyC escape fractions. However, the extreme [O III]/[O II] ratios of the z=0.1-0.3 Green Pea galaxies may be due to high escape fractions of ionizing radiation. To analyze the LyC optical depths and ionizing sources of these rare, compact starbursts, we compare nebular photoionization and stellar population models with observed emission lines in the Peas' SDSS spectra. We focus on the six most extreme Green Peas, the galaxies with the highest [O III]/[O II] ratios and the best candidates for escaping ionizing radiation. The Balmer line equivalent widths and He I {\lambda}3819 emission in the extreme Peas support young ages of 3-5 Myr, and He II {\lambda}4686 emission in five extreme Peas signals the presence of hard ionizing sources. Ionization by active galactic nuclei or high-mass X-ray binaries is inconsistent with the Peas' line ratios and ages. Although stacked spectra reveal no Wolf-Rayet (WR) features, we tentatively detect WR features in the SDSS spectra of three extreme Peas. Based on the Peas' ages and line ratios, we find that WR stars, chemically homogeneous O stars, or shocks could produce the observed He II emission. If hot stars are responsible, the Peas' optical depths are ambiguous. However, accounting for emission from shocks lowers the inferred optical depth and suggests that the Peas may be optically thin. The Peas' ages likely optimize the escape of Lyman-continuum radiation; they are old enough for supernovae and stellar winds to reshape the interstellar medium, but young enough to possess large numbers of UV-luminous O or WR stars.


title: Some remarks on the Pigola-Rigoli-Setti version of the Omori-Yau maximum principle

abstract: We prove that the hypotheses in the version of the Omori-Yau maximum principle that was given by Pigola-Rigoli-Setti are logically equivalent to the assumption that the manifold carries a $C^2$ proper function whose gradient and Hessian (Laplacian) are bounded. In particular, this result extends the scope of the original Omori-Yau principle, formulated in terms of lower bounds for curvature.


title: Shadow of a Kaluza-Klein rotating dilaton black hole

abstract: We study the shadow produced by a spinning Kaluza-Klein black hole in Einstein gravity coupled to a Maxwell field and a dilaton. The size and the shape of the shadow depend on the mass, the charge, and the angular momentum. We find that, for fixed values of these parameters, the shadow is slightly larger and less deformed than for its Kerr-Newman counterpart.


title: A robust Bayesian approach to modelling epistemic uncertainty in common-cause failure models

abstract: In a standard Bayesian approach to the alpha-factor model for common-cause failure, a precise Dirichlet prior distribution models epistemic uncertainty in the alpha-factors. This Dirichlet prior is then updated with observed data to obtain a posterior distribution, which forms the basis for further inferences. In this paper, we adapt the imprecise Dirichlet model of Walley to represent epistemic uncertainty in the alpha-factors. In this approach, epistemic uncertainty is expressed more cautiously via lower and upper expectations for each alpha-factor, along with a learning parameter which determines how quickly the model learns from observed data. For this application, we focus on elicitation of the learning parameter, and find that values in the range of 1 to 10 seem reasonable. The approach is compared with Kelly and Atwood's minimally informative Dirichlet prior for the alpha-factor model, which incorporated precise mean values for the alpha-factors, but which was otherwise quite diffuse. Next, we explore the use of a set of Gamma priors to model epistemic uncertainty in the marginal failure rate, expressed via a lower and upper expectation for this rate, again along with a learning parameter. As zero counts are generally less of an issue here, we find that the choice of this learning parameter is less crucial. Finally, we demonstrate how both epistemic uncertainty models can be combined to arrive at lower and upper expectations for all common-cause failure rates. Thereby, we effectively provide a full sensitivity analysis of common-cause failure rates, properly reflecting epistemic uncertainty of the analyst on all levels of the common-cause failure model.


title: Follow the Leader If You Can, Hedge If You Must

abstract: Follow-the-Leader (FTL) is an intuitive sequential prediction strategy that guarantees constant regret in the stochastic setting, but has terrible performance for worst-case data. Other hedging strategies have better worst-case guarantees but may perform much worse than FTL if the data are not maximally adversarial. We introduce the FlipFlop algorithm, which is the first method that provably combines the best of both worlds. As part of our construction, we develop AdaHedge, which is a new way of dynamically tuning the learning rate in Hedge without using the doubling trick. AdaHedge refines a method by Cesa-Bianchi, Mansour and Stoltz (2007), yielding slightly improved worst-case guarantees. By interleaving AdaHedge and FTL, the FlipFlop algorithm achieves regret within a constant factor of the FTL regret, without sacrificing AdaHedge's worst-case guarantees. AdaHedge and FlipFlop do not need to know the range of the losses in advance; moreover, unlike earlier methods, both have the intuitive property that the issued weights are invariant under rescaling and translation of the losses. The losses are also allowed to be negative, in which case they may be interpreted as gains.


title: First-order sideband transitions with flux-driven asymmetric transmon qubits

abstract: We demonstrate rapid, first-order sideband transitions between a superconducting resonator and a frequency-modulated transmon qubit. The qubit contains a substantial asymmetry between its Josephson junctions leading to a linear portion of the energy band near the resonator frequency. The sideband transitions are driven with a magnetic flux signal of a few hundred MHz coupled to the qubit. This modulates the qubit splitting at a frequency near the detuning between the dressed qubit and resonator frequencies, leading to rates up to 85 MHz for exchanging quanta between the qubit and resonator.


title: Nonlinear instability in simulations of Large Plasma Device turbulence

abstract: Several simulations of turbulence in the Large Plasma Device (LAPD) [W. Gekelman et al., Rev. Sci. Inst. 62, 2875 (1991)] are energetically analyzed and compared with each other and with the experiment. The simulations use the same model, but different axial boundary conditions. They employ either periodic, zero-value, zero-derivative, or sheath axial boundaries. The linear stability physics is different between the scenarios because the various boundary conditions allow the drift wave instability to access different axial structures, and the sheath boundary simulation contains a conducting wall mode instability which is just as unstable as the drift waves. Nevertheless, the turbulence in all the simulations is relatively similar because it is primarily driven by a robust nonlinear instability that is the same for all cases. The nonlinear instability preferentially drives $k_\parallel = 0$ potential energy fluctuations, which then three-wave couple to $k_\parallel \ne 0$ potential energy fluctuations in order to access the adiabatic response to transfer their energy to kinetic energy fluctuations. The turbulence self-organizes to drive this nonlinear instability, which destroys the linear eigenmode structures, making the linear instabilities ineffective.


title: Holography without translational symmetry

abstract: We propose massive gravity as a holographic framework for describing a class of strongly interacting quantum field theories with broken translational symmetry. Bulk gravitons are assumed to have a Lorentz-breaking mass term as a substitute for spatial inhomogeneities. This breaks momentum-conservation in the boundary field theory. At finite chemical potential, the gravity duals are charged black holes in asymptotically anti-de Sitter spacetime. The conductivity in these systems generally exhibits a Drude peak that approaches a delta function in the massless gravity limit. Furthermore, the optical conductivity shows an emergent scaling law: $|\sigma(\omega)| \approx {A \over \omega^{\alpha}} + B$. This result is consistent with that found earlier by Horowitz, Santos, and Tong who introduced an explicit inhomogeneous lattice into the system.


title: Signatures of the topological spin of Josephson vortices in topological superconductors

abstract: We consider a modified setup for measuring the Aharonov-Casher phase which consists of a Josephson vortex trapped in an annular topological superconducting junction. The junction encloses both electric charge and magnetic flux. We discover a deviation from the Aharonov-Casher prediction whose origin we identify in an additive universal topological phase that remarkably depends only on the parity of the number of vortices enclosed by the junction. We show that this phase is $\pm 2\pi$ times the topological spin of the Josephson vortex and is proportional to the Chern number. The presence of this phase can be measured through its effect on the junction's voltage characteristics, thus revealing the topological properties of the Josephson vortex and the superconducting state.


title: Stellar Populations in the Central 0.5 pc of the Galaxy I: A New Method for Constructing Luminosity Functions and Surface-Density Profiles

abstract: We present new high angular resolution near-infrared spectroscopic observations of the nuclear star cluster surrounding the Milky Way's central supermassive black hole. Using the integral-field spectrograph OSIRIS on Keck II behind the laser-guide-star adaptive optics system, this spectroscopic survey enables us to separate early-type (young, 4-6 Myr) and late-type (old, >1 Gyr) stars with a completeness of 50% down to K' = 15.5 mag, which corresponds to ~10 \msun for the early-type stars. This work increases the radial extent of reported OSIRIS/Keck measurements by more than a factor of 3 from 4" to 14" (0.16 pc to 0.56 pc), along the projected disk of young stars. For our analysis, we implement a new method of completeness correction using a combination of star-planting simulations and Bayesian inference. We assign probabilities for the spectral type of every source detected in deep imaging down to K' = 15.5 mag using information from spectra, simulations, number counts, and the distribution of stars. The inferred radial surface-density profiles, $\Sigma(R) \propto R^{-\Gamma}$, for the young stars and late-type giants are consistent with earlier results ($\Gamma_{early} = 0.93 \pm 0.09$, $\Gamma_{late} = 0.16 \pm 0.07$). The late-type surface-density profile is approximately flat out to the edge of the survey. While the late-type stellar luminosity function is consistent with the Galactic bulge, the completeness-corrected luminosity function of the early-type stars has significantly more young stars at faint magnitudes compared to previous surveys with similar depth. This luminosity function indicates that the corresponding mass function of the young stars is likely less top-heavy than that inferred from previous surveys.


title: Stellar Populations in the Central 0.5 pc of the Galaxy II: The Initial Mass Function

abstract: The supermassive black hole at the center of the Milky Way plays host to a massive, young cluster that may have formed in one of the most inhospitable environments in the Galaxy. We present new measurements of the global properties of this cluster, including the initial mass function (IMF), age, and cluster mass. These results are based on Keck laser-guide-star adaptive optics observations used to identify the young stars and measure their Kp-band luminosity function as presented in Do et al. 2013. A Bayesian inference methodology is developed to simultaneously fit the global properties of the cluster utilizing the observations and extensive simulations of synthetic star clusters. We find that the slope of the mass function for this cluster is \alpha = 1.7 +/- 0.2, which is steeper than previously reported, but still flatter than the traditional Salpeter slope of 2.35. The age of the cluster is between 2.5-5.8 Myr with 95% confidence, which is a younger age than typically adopted but consistent within the uncertainties of past measurements. The exact age of the cluster is difficult to determine since our results show two distinct age solutions (3.9 Myr and 2.8 Myr) due to model degeneracies in the relative number of Wolf-Rayet and OB stars. The total cluster mass is between 14,000 - 37,000 \msun above 1 \msun and it is necessary to include multiple star systems in order to fit the observed luminosity function and the number of observed Wolf-Rayet stars. The new IMF slope measurement is now consistent with X-ray observations indicating a factor of 10 fewer X-ray emitting pre-main-sequence stars than expected when compared with a Salpeter IMF. The young cluster at the Galactic center is one of the few definitive examples of an IMF that deviates significantly from the near-universal IMFs found in the solar neighborhood.


title: Coloring triangle-free rectangle overlap graphs with $O(\log\log n)$ colors

abstract: Recently, it was proved that triangle-free intersection graphs of $n$ line segments in the plane can have chromatic number as large as $\Theta(\log\log n)$. Essentially the same construction produces $\Theta(\log\log n)$-chromatic triangle-free intersection graphs of a variety of other geometric shapes---those belonging to any class of compact arc-connected sets in $\mathbb{R}^2$ closed under horizontal scaling, vertical scaling, and translation, except for axis-parallel rectangles. We show that this construction is asymptotically optimal for intersection graphs of boundaries of axis-parallel rectangles, which can be alternatively described as overlap graphs of axis-parallel rectangles. That is, we prove that triangle-free rectangle overlap graphs have chromatic number $O(\log\log n)$, improving on the previous bound of $O(\log n)$. To this end, we exploit a relationship between off-line coloring of rectangle overlap graphs and on-line coloring of interval overlap graphs. Our coloring method decomposes the graph into a bounded number of subgraphs with a tree-like structure that "encodes" strategies of the adversary in the on-line coloring problem. Then, these subgraphs are colored with $O(\log\log n)$ colors using a combination of techniques from on-line algorithms (first-fit) and data structure design (heavy-light decomposition).


title: Eventual linear convergence of the Douglas Rachford iteration for basis pursuit

abstract: We provide a simple analysis of the Douglas-Rachford splitting algorithm in the context of $\ell^1$ minimization with linear constraints, and quantify the asymptotic linear convergence rate in terms of principal angles between relevant vector spaces. In the compressed sensing setting, we show how to bound this rate in terms of the restricted isometry constant. More general iterative schemes obtained by $\ell^2$-regularization and over-relaxation including the dual split Bregman method are also treated, which answers the question how to choose the relaxation and soft-thresholding parameters to accelerate the asymptotic convergence rate. We make no attempt at characterizing the transient regime preceding the onset of linear convergence.


title: Timelike U-dualities in Generalised Geometry

abstract: We study timelike U-dualities acting in three and four directions of 11-dimensional supergravity, which form the groups $SL(2)\times SL(3)$ and SL(5). Using generalised geometry, we find that timelike U-dualities, despite previous conjectures, do not change the signature of the spacetime. Furthermore, we prove that the spacetime signature must be $(-,+,...,+)$ when the U-duality modular group is either $\frac{SL(2)\times SL(3)}{SO(1,1)\times SO(2,1)}$ or $\frac{SL(5)}{SO(3,2)}$. We find that for some dual solutions it is necessary to include a trivector field which is related to the existence of non-geometric fluxes in lower dimensions. In the second part of the paper, we explicitly study the action of the dualities on supergravity solutions corresponding to M2-branes. For a finite range of the transformation, the action of $SL(2)\otimes SL(3)$ on the worldvolume of uncharged M2-branes charges them while it changes the charge of extreme M2-branes. It thus acts as a Harrison transformation. At the limits of the range, we obtain the "subtracted geometries" which correspond to an infinite Harrison boost. Outside this range the trivector field becomes non-zero and we obtain a dual solution that cannot be uniquely written in terms of a metric, 3-form and trivector. Instead it corresponds to a family of solutions linked by a local SO(1,1) rotation. The SL(5) duality is used to act on a smeared extreme M2-brane giving a brane-like solution carrying momentum in the transverse direction that the brane was delocalised along.


title: Finite temperature effects on anisotropic pressure and equation of state of dense neutron matter in an ultrastrong magnetic field

abstract: Spin polarized states in dense neutron matter with recently developed Skyrme effective interaction (BSk20 parametrization) are considered in the magnetic fields $H$ up to $10^{20}$ G at finite temperature. In a strong magnetic field, the total pressure in neutron matter is anisotropic, and the difference between the pressures parallel and perpendicular to the field direction becomes significant at $H>H_{th}\sim10^{18}$ G. The longitudinal pressure decreases with the magnetic field and vanishes in the critical field $10^{18}<H_c\lesssim10^{19}$ G, resulting in the longitudinal instability of neutron matter. With increasing the temperature, the threshold $H_{th}$ and critical $H_c$ magnetic fields also increase. The appearance of the longitudinal instability prevents the formation of a fully spin polarized state in neutron matter and only the states with moderate spin polarization are accessible. The anisotropic equation of state is determined at densities and temperatures relevant for the interiors of magnetars. The entropy of strongly magnetized neutron matter turns out to be larger than the entropy of the nonpolarized matter. This is caused by some specific details in the dependence of the entropy on the effective masses of neutrons with spin up and spin down in a polarized state.


title: A Search for Photometric Variability in L and T type Brown Dwarf Atmospheres

abstract: Using the Gemini infrared camera on the 3-meter Shane telescope at Lick Observatory, we have searched for broad-band J and K' photometric variability for a sample of 15 L and T-type brown dwarfs, including 7 suspected spectral binaries. Four of the dwarfs - 2MASS J0939-2448, 2MASS J1416+1348A, 2MASS J1711+2232, and 2MASS J2139+0220 - exhibit statistically significant variations over timescales ranging from ~0.5 hr to 6 days. Our detection of variability in 2MASS J2139+0220 confirms that reported by Radigan et al., and periodogram and phase dispersion minimization analysis also confirms a variability period of approximately 7.6 +/- 0.2 hours. Remarkably, two of the four variables are known or candidate binary systems, including 2MASS J2139+0220, for which we find only marginal evidence of radial velocity variation over the course of a year. This result suggests that some spectral binary candidates may appear as such due to the blending of cloudy and non-cloudy regions in a single "patchy" atmosphere. Our results are consistent with an overall variability fraction of 35 +/- 5 %, with no clear evidence of greater variability among brown dwarfs at the L dwarf/T dwarf transition.


title: A three-phase free boundary problem with melting ice and dissolving gas

abstract: We develop a mathematical model for a three-phase free boundary problem in one dimension that involves the interactions between gas, water and ice. The dynamics are driven by melting of the ice layer, while the pressurized gas also dissolves within the meltwater. The model incorporates a Stefan condition at the water-ice interface along with Henry's law for dissolution of gas at the gas-water interface. We employ a quasi-steady approximation for the phase temperatures and then derive a series solution for the interface positions. A non-standard feature of the model is an integral free boundary condition that arises from mass conservation owing to changes in gas density at the gas-water interface, which makes the problem non-self-adjoint. We derive a two-scale asymptotic series solution for the dissolved gas concentration, which because of the non-self-adjointness gives rise to a Fourier series expansion in eigenfunctions that do not satisfy the usual orthogonality conditions. Numerical simulations of the original governing equations are used to validate the series approximations.


title: A Comparison of Bimolecular Reaction Models for Stochastic Reaction Diffusion Systems

abstract: Stochastic reaction-diffusion models have become an important tool in studying how both noise in the chemical reaction process and the spatial movement of molecules influences the behavior of biological systems. There are two primary spatially-continuous models that have been used in recent studies: the diffusion limited reaction model of Smoluchowski, and a second approach popularized by Doi. Both models treat molecules as points undergoing Brownian motion. The former represents chemical reactions between two reactants through the use of reactive boundary conditions, with two molecules reacting instantly upon reaching a fixed separation (called the reaction-radius). The Doi model uses reaction potentials, whereby two molecules react with a fixed probability per unit time, $\lambda$, when separated by less than the reaction radius. In this work we study the rigorous relationship between the two models. For the special case of a protein diffusing to a fixed DNA binding site, we prove that the solution to the Doi model converges to the solution of the Smoluchowski model as $\lambda \to \infty$, with a rigorous $O(\lambda^{-1/2 + \epsilon})$ error bound (for any fixed $\epsilon > 0$). We investigate by numerical simulation, for biologically relevant parameter values, the difference between the solutions and associated reaction time statistics of the two models. As the reaction-radius is decreased, for sufficiently large but fixed values of $\lambda$, these differences are found to increase like the inverse of the binding radius.


title: Enhancement of level-crossing resonances in rubidium atoms by frequency control of the exciting radiation field

abstract: We studied magneto-optical resonances caused by excited-state level crossings in a nonzero magnetic field. Experimental measurements were performed on the transitions of the $D_2$ line of rubidium. These measured signals were described by a theoretical model that takes into account all neighboring hyperfine transitions, the mixing of magnetic sublevels in an external magnetic field, the coherence properties of the exciting laser radiation, and the Doppler effect. Good agreement between the experimental measurements and the theoretical model could be achieved over a wide range of laser power densities. We further showed that the contrasts of the level-crossing peaks can be sensitive to changes in the frequency of the exciting laser radiation as small as several tens of megahertz when the hyperfine splitting of the exciting state is larger than the Doppler broadening.


title: Quantum correlations of two-mode Gaussian systems in a thermal environment

abstract: In the framework of the theory of open systems based on completely positive quantum dynamical semigroups, we give a description of continuous variable quantum entanglement and quantum discord for a system consisting of two non-interacting non-resonant bosonic modes embedded in a thermal environment. We study the time evolution of logarithmic negativity, which characterizes the degree of entanglement, and show that in the case of an entangled initial squeezed thermal state, entanglement suppression takes place for all temperatures of the environment, including zero temperature. We analyze the time evolution of the Gaussian quantum discord, which is a measure of all quantum correlations in the bipartite state, including entanglement, and show that discord decays asymptotically in time under the effect of the thermal bath. We describe also the time evolution of classical correlations and quantum mutual information, which measures the total correlations of the quantum system.


title: Markov Equivalence Classes for Maximal Ancestral Graphs

abstract: Ancestral graphs are a class of graphs that encode conditional independence relations arising in DAG models with latent and selection variables, corresponding to marginalization and conditioning. However, for any ancestral graph, there may be several other graphs to which it is Markov equivalent. We introduce a simple representation of a Markov equivalence class of ancestral graphs, thereby facilitating model search. \ More specifically, we define a join operation on ancestral graphs which will associate a unique graph with a Markov equivalence class. We also extend the separation criterion for ancestral graphs (which is an extension of d-separation) and provide a proof of the pairwise Markov property for joined ancestral graphs.


title: Learning Hierarchical Object Maps Of Non-Stationary Environments with mobile robots

abstract: Building models, or maps, of robot environments is a highly active research area; however, most existing techniques construct unstructured maps and assume static environments. In this paper, we present an algorithm for learning object models of non-stationary objects found in office-type environments. Our algorithm exploits the fact that many objects found in office environments look alike (e.g., chairs, recycling bins). It does so through a two-level hierarchical representation, which links individual objects with generic shape templates of object classes. We derive an approximate EM algorithm for learning shape parameters at both levels of the hierarchy, using local occupancy grid maps for representing shape. Additionally, we develop a Bayesian model selection algorithm that enables the robot to estimate the total number of objects and object templates in the environment. Experimental results using a real robot equipped with a laser range finder indicate that our approach performs well at learning object-based maps of simple office environments. The approach outperforms a previously developed non-hierarchical algorithm that models objects but lacks class templates.


title: A constraint satisfaction approach to the robust spanning tree problem with interval data

abstract: Robust optimization is one of the fundamental approaches to deal with uncertainty in combinatorial optimization. This paper considers the robust spanning tree problem with interval data, which arises in a variety of telecommunication applications. It proposes a constraint satisfaction approach using a combinatorial lower bound, a pruning component that removes infeasible and suboptimal edges, as well as a search strategy exploring the most uncertain edges first. The resulting algorithm is shown to produce very dramatic improvements over the mathematical programming approach of Yaman et al. and to enlarge considerably the class of problems amenable to effective solutions


title: On the Construction of the Inclusion Boundary Neighbourhood for Markov Equivalence Classes of Bayesian Network Structures

abstract: The problem of learning Markov equivalence classes of Bayesian network structures may be solved by searching for the maximum of a scoring metric in a space of these classes. This paper deals with the definition and analysis of one such search space. We use a theoretically motivated neighbourhood, the inclusion boundary, and represent equivalence classes by essential graphs. We show that this search space is connected and that the score of the neighbours can be evaluated incrementally. We devise a practical way of building this neighbourhood for an essential graph that is purely graphical and does not explicitely refer to the underlying independences. We find that its size can be intractable, depending on the complexity of the essential graph of the equivalence class. The emphasis is put on the potential use of this space with greedy hill -climbing search


title: Tree-dependent Component Analysis

abstract: We present a generalization of independent component analysis (ICA), where instead of looking for a linear transform that makes the data components independent, we look for a transform that makes the data components well fit by a tree-structured graphical model. Treating the problem as a semiparametric statistical problem, we show that the optimal transform is found by minimizing a contrast function based on mutual information, a function that directly extends the contrast function used for classical ICA. We provide two approximations of this contrast function, one using kernel density estimation, and another using kernel generalized variance. This tree-dependent component analysis framework leads naturally to an efficient general multivariate density estimation technique where only bivariate density estimation needs to be performed.


title: Bipolar Possibilistic Representations

abstract: Recently, it has been emphasized that the possibility theory framework allows us to distinguish between i) what is possible because it is not ruled out by the available knowledge, and ii) what is possible for sure. This distinction may be useful when representing knowledge, for modelling values which are not impossible because they are consistent with the available knowledge on the one hand, and values guaranteed to be possible because reported from observations on the other hand. It is also of interest when expressing preferences, to point out values which are positively desired among those which are not rejected. This distinction can be encoded by two types of constraints expressed in terms of necessity measures and in terms of guaranteed possibility functions, which induce a pair of possibility distributions at the semantic level. A consistency condition should ensure that what is claimed to be guaranteed as possible is indeed not impossible. The present paper investigates the representation of this bipolar view, including the case when it is stated by means of conditional measures, or by means of comparative context-dependent constraints. The interest of this bipolar framework, which has been recently stressed for expressing preferences, is also pointed out in the representation of diagnostic knowledge.


title: Learning with Scope, with Application to Information Extraction and Classification

abstract: In probabilistic approaches to classification and information extraction, one typically builds a statistical model of words under the assumption that future data will exhibit the same regularities as the training data. In many data sets, however, there are scope-limited features whose predictive power is only applicable to a certain subset of the data. For example, in information extraction from web pages, word formatting may be indicative of extraction category in different ways on different web pages. The difficulty with using such features is capturing and exploiting the new regularities encountered in previously unseen data. In this paper, we propose a hierarchical probabilistic model that uses both local/scope-limited features, such as word formatting, and global features, such as word content. The local regularities are modeled as an unobserved random parameter which is drawn once for each local data set. This random parameter is estimated during the inference process and then used to perform classification with both the local and global features--- a procedure which is akin to automatically retuning the classifier to the local regularities on each newly encountered web page. Exact inference is intractable and we present approximations via point estimates and variational methods. Empirical results on large collections of web data demonstrate that this method significantly improves performance from traditional models of global features alone.


title: Qualitative MDPs and POMDPs: An Order-Of-Magnitude Approximation

abstract: We develop a qualitative theory of Markov Decision Processes (MDPs) and Partially Observable MDPs that can be used to model sequential decision making tasks when only qualitative information is available. Our approach is based upon an order-of-magnitude approximation of both probabilities and utilities, similar to epsilon-semantics. The result is a qualitative theory that has close ties with the standard maximum-expected-utility theory and is amenable to general planning techniques.


title: Introducing Variable Importance Tradeoffs into CP-Nets

abstract: The ability to make decisions and to assess potential courses of action is a corner-stone of many AI applications, and usually this requires explicit information about the decision-maker s preferences. IN many applications, preference elicitation IS a serious bottleneck.The USER either does NOT have the time, the knowledge, OR the expert support required TO specify complex multi - attribute utility functions. IN such cases, a method that IS based ON intuitive, yet expressive, preference statements IS required. IN this paper we suggest the USE OF TCP - nets, an enhancement OF CP - nets, AS a tool FOR representing, AND reasoning about qualitative preference statements.We present AND motivate this framework, define its semantics, AND show how it can be used TO perform constrained optimization.


title: Planning under Continuous Time and Resource Uncertainty: A Challenge for AI

abstract: We outline a class of problems, typical of Mars rover operations, that are problematic for current methods of planning under uncertainty. The existing methods fail because they suffer from one or more of the following limitations: 1) they rely on very simple models of actions and time, 2) they assume that uncertainty is manifested in discrete action outcomes, 3) they are only practical for very small problems. For many real world problems, these assumptions fail to hold. In particular, when planning the activities for a Mars rover, none of the above assumptions is valid: 1) actions can be concurrent and have differing durations, 2) there is uncertainty concerning action durations and consumption of continuous resources like power, and 3) typical daily plans involve on the order of a hundred actions. This class of problems may be of particular interest to the UAI community because both classical and decision-theoretic planning techniques may be useful in solving it. We describe the rover problem, discuss previous work on planning under uncertainty, and present a detailed, but very small, example illustrating some of the difficulties of finding good plans.


title: Generalized Instrumental Variables

abstract: This paper concerns the assessment of direct causal effects from a combination of: (i) non-experimental data, and (ii) qualitative domain knowledge. Domain knowledge is encoded in the form of a directed acyclic graph (DAG), in which all interactions are assumed linear, and some variables are presumed to be unobserved. We provide a generalization of the well-known method of Instrumental Variables, which allows its application to models with few conditional independeces.


title: Finding Optimal Bayesian Networks

abstract: In this paper, we derive optimality results for greedy Bayesian-network search algorithms that perform single-edge modifications at each step and use asymptotically consistent scoring criteria. Our results extend those of Meek (1997) and Chickering (2002), who demonstrate that in the limit of large datasets, if the generative distribution is perfect with respect to a DAG defined over the observable variables, such search algorithms will identify this optimal (i.e. generative) DAG model. We relax their assumption about the generative distribution, and assume only that this distribution satisfies the {em composition property} over the observable variables, which is a more realistic assumption for real domains. Under this assumption, we guarantee that the search algorithms identify an {em inclusion-optimal} model; that is, a model that (1) contains the generative distribution and (2) has no sub-model that contains this distribution. In addition, we show that the composition property is guaranteed to hold whenever the dependence relationships in the generative distribution can be characterized by paths between singleton elements in some generative graphical model (e.g. a DAG, a chain graph, or a Markov network) even when the generative model includes unobserved variables, and even when the observed data is subject to selection bias.


title: Continuation Methods for Mixing Heterogenous Sources

abstract: A number of modern learning tasks involve estimation from heterogeneous information sources. This includes classification with labeled and unlabeled data as well as other problems with analogous structure such as competitive (game theoretic) problems. The associated estimation problems can be typically reduced to solving a set of fixed point equations (consistency conditions). We introduce a general method for combining a preferred information source with another in this setting by evolving continuous paths of fixed points at intermediate allocations. We explicitly identify critical points along the unique paths to either increase the stability of estimation or to ensure a significant departure from the initial source. The homotopy continuation approach is guaranteed to terminate at the second source, and involves no combinatorial effort. We illustrate the power of these ideas both in classification tasks with labeled and unlabeled data, as well as in the context of a competitive (min-max) formulation of DNA sequence motif discovery.


title: Interpolating Conditional Density Trees

abstract: Joint distributions over many variables are frequently modeled by decomposing them into products of simpler, lower-dimensional conditional distributions, such as in sparsely connected Bayesian networks. However, automatically learning such models can be very computationally expensive when there are many datapoints and many continuous variables with complex nonlinear relationships, particularly when no good ways of decomposing the joint distribution are known a priori. In such situations, previous research has generally focused on the use of discretization techniques in which each continuous variable has a single discretization that is used throughout the entire network. \ In this paper, we present and compare a wide variety of tree-based algorithms for learning and evaluating conditional density estimates over continuous variables. These trees can be thought of as discretizations that vary according to the particular interactions being modeled; however, the density within a given leaf of the tree need not be assumed constant, and we show that such nonuniform leaf densities lead to more accurate density estimation. We have developed Bayesian network structure-learning algorithms that employ these tree-based conditional density representations, and we show that they can be used to practically learn complex joint probability models over dozens of continuous variables from thousands of datapoints. We focus on finding models that are simultaneously accurate, fast to learn, and fast to evaluate once they are learned.


title: Iterative Join-Graph Propagation

abstract: The paper presents an iterative version of join-tree clustering that applies the message passing of join-tree clustering algorithm to join-graphs rather than to join-trees, iteratively. It is inspired by the success of Pearl's belief propagation algorithm as an iterative approximation scheme on one hand, and by a recently introduced mini-clustering i. success as an anytime approximation method, on the other. The proposed Iterative Join-graph Propagation IJGP belongs to the class of generalized belief propagation methods, recently proposed using analogy with algorithms in statistical physics. Empirical evaluation of this approach on a number of problem classes demonstrates that even the most time-efficient variant is almost always superior to IBP and MC i, and is sometimes more accurate by as much as several orders of magnitude.


title: An Information-Theoretic External Cluster-Validity Measure

abstract: In this paper we propose a measure of clustering quality or accuracy that is appropriate in situations where it is desirable to evaluate a clustering algorithm by somehow comparing the clusters it produces with ``ground truth' consisting of classes assigned to the patterns by manual means or some other means in whose veracity there is confidence. Such measures are refered to as ``external'. Our measure also has the characteristic of allowing clusterings with different numbers of clusters to be compared in a quantitative and principled way. Our evaluation scheme quantitatively measures how useful the cluster labels of the patterns are as predictors of their class labels. In cases where all clusterings to be compared have the same number of clusters, the measure is equivalent to the mutual information between the cluster labels and the class labels. In cases where the numbers of clusters are different, however, it computes the reduction in the number of bits that would be required to encode (compress) the class labels if both the encoder and decoder have free acccess to the cluster labels. To achieve this encoding the estimated conditional probabilities of the class labels given the cluster labels must also be encoded. These estimated probabilities can be seen as a model for the class labels and their associated code length as a model cost.


title: Causes and Explanations in the Structural-Model Approach: Tractable Cases

abstract: In this paper, we continue our research on the algorithmic aspects of Halpern and Pearl's causes and explanations in the structural-model approach. To this end, we present new characterizations of weak causes for certain classes of causal models, which show that under suitable restrictions deciding causes and explanations is tractable. To our knowledge, these are the first explicit tractability results for the structural-model approach.


title: The Thing That We Tried Didn't Work Very Well : Deictic Representation in Reinforcement Learning

abstract: Most reinforcement learning methods operate on propositional representations of the world state. Such representations are often intractably large and generalize poorly. Using a deictic representation is believed to be a viable alternative: they promise generalization while allowing the use of existing reinforcement-learning methods. Yet, there are few experiments on learning with deictic representations reported in the literature. In this paper we explore the effectiveness of two forms of deictic representation and a na\"{i}ve propositional representation in a simple blocks-world domain. We find, empirically, that the deictic representations actually worsen learning performance. We conclude with a discussion of possible causes of these results and strategies for more effective learning in domains with objects.


title: Factorization of Discrete Probability Distributions

abstract: We formulate necessary and sufficient conditions for an arbitrary discrete probability distribution to factor according to an undirected graphical model, or a log-linear model, or other more general exponential models. This result generalizes the well known Hammersley-Clifford Theorem.


title: Statistical Decisions Using Likelihood Information Without Prior Probabilities

abstract: This paper presents a decision-theoretic approach to statistical inference that satisfies the likelihood principle (LP) without using prior information. Unlike the Bayesian approach, which also satisfies LP, we do not assume knowledge of the prior distribution of the unknown parameter. With respect to information that can be obtained from an experiment, our solution is more efficient than Wald's minimax solution.However, with respect to information assumed to be known before the experiment, our solution demands less input than the Bayesian solution.


title: Reduction of Maximum Entropy Models to Hidden Markov Models

abstract: We show that maximum entropy (maxent) models can be modeled with certain kinds of HMMs, allowing us to construct maxent models with hidden variables, hidden state sequences, or other characteristics. The models can be trained using the forward-backward algorithm. While the results are primarily of theoretical interest, unifying apparently unrelated concepts, we also give experimental results for a maxent model with a hidden variable on a word disambiguation task; the model outperforms standard techniques.


title: Distributed Planning in Hierarchical Factored MDPs

abstract: We present a principled and efficient planning algorithm for collaborative multiagent dynamical systems. All computation, during both the planning and the execution phases, is distributed among the agents; each agent only needs to model and plan for a small part of the system. Each of these local subsystems is small, but once they are combined they can represent an exponentially larger problem. The subsystems are connected through a subsystem hierarchy. Coordination and communication between the agents is not imposed, but derived directly from the structure of this hierarchy. A globally consistent plan is achieved by a message passing algorithm, where messages correspond to natural local reward functions and are computed by local linear programs; another message passing algorithm allows us to execute the resulting policy. When two portions of the hierarchy share the same structure, our algorithm can reuse plans and messages to speed up computation.


title: Expectation Propogation for approximate inference in dynamic Bayesian networks

abstract: We describe expectation propagation for approximate inference in dynamic Bayesian networks as a natural extension of Pearl s exact belief propagation.Expectation propagation IS a greedy algorithm, converges IN many practical cases, but NOT always.We derive a DOUBLE - loop algorithm, guaranteed TO converge TO a local minimum OF a Bethe free energy.Furthermore, we show that stable fixed points OF (damped) expectation propagation correspond TO local minima OF this free energy, but that the converse need NOT be the CASE .We illustrate the algorithms BY applying them TO switching linear dynamical systems AND discuss implications FOR approximate inference IN general Bayesian networks.


title: Coordinates: Probabilistic Forecasting of Presence and Availability

abstract: We present methods employed in Coordinate, a prototype service that supports collaboration and communication by learning predictive models that provide forecasts of users s AND availability.We describe how data IS collected about USER activity AND proximity FROM multiple devices, IN addition TO analysis OF the content OF users, the time of day, and day of week. We review applications of presence forecasting embedded in the Priorities application and then present details of the Coordinate service that was informed by the earlier efforts.


title: Unconstrained Influence Diagrams

abstract: We extend the language of influence diagrams to cope with decision scenarios where the order of decisions and observations is not determined. As the ordering of decisions is dependent on the evidence, a step-strategy of such a scenario is a sequence of dependent choices of the next action. A strategy is a step-strategy together with selection functions for decision actions. The structure of a step-strategy can be represented as a DAG with nodes labeled with action variables. We introduce the concept of GS-DAG: a DAG incorporating an optimal step-strategy for any instantiation. We give a method for constructing GS-DAGs, and we show how to use a GS-DAG for determining an optimal strategy. Finally we discuss how analysis of relevant past can be used to reduce the size of the GS-DAG.


title: CFW: A Collaborative Filtering System Using Posteriors Over Weights Of Evidence

abstract: We describe CFW, a computationally efficient algorithm for collaborative filtering that uses posteriors over weights of evidence. In experiments on real data, we show that this method predicts as well or better than other methods in situations where the size of the user query is small. The new approach works particularly well when the user s query CONTAINS low frequency(unpopular) items.The approach complements that OF dependency networks which perform well WHEN the size OF the query IS large.Also IN this paper, we argue that the USE OF posteriors OVER weights OF evidence IS a natural way TO recommend similar items collaborative - filtering task.


title: A Bayesian Network Scoring Metric That Is Based On Globally Uniform Parameter Priors

abstract: We introduce a new Bayesian network (BN) scoring metric called the Global Uniform (GU) metric. This metric is based on a particular type of default parameter prior. Such priors may be useful when a BN developer is not willing or able to specify domain-specific parameter priors. The GU parameter prior specifies that every prior joint probability distribution P consistent with a BN structure S is considered to be equally likely. Distribution P is consistent with S if P includes just the set of independence relations defined by S. We show that the GU metric addresses some undesirable behavior of the BDeu and K2 Bayesian network scoring metrics, which also use particular forms of default parameter priors. A closed form formula for computing GU for special classes of BNs is derived. Efficiently computing GU for an arbitrary BN remains an open problem.


title: Efficient Nash Computation in Large Population Games with Bounded Influence

abstract: We introduce a general representation of large-population games in which each player s influence ON the others IS centralized AND limited, but may otherwise be arbitrary.This representation significantly generalizes the class known AS congestion games IN a natural way.Our main results are provably correct AND efficient algorithms FOR computing AND learning approximate Nash equilibria IN this general framework.


title: Dimension Correction for Hierarchical Latent Class Models

abstract: Model complexity is an important factor to consider when selecting among graphical models. When all variables are observed, the complexity of a model can be measured by its standard dimension, i.e. the number of independent parameters. When hidden variables are present, however, standard dimension might no longer be appropriate. One should instead use effective dimension (Geiger et al. 1996). This paper is concerned with the computation of effective dimension. First we present an upper bound on the effective dimension of a latent class (LC) model. This bound is tight and its computation is easy. We then consider a generalization of LC models called hierarchical latent class (HLC) models (Zhang 2002). We show that the effective dimension of an HLC model can be obtained from the effective dimensions of some related LC models. We also demonstrate empirically that using effective dimension in place of standard dimension improves the quality of models learned from data.


title: Almost-everywhere algorithmic stability and generalization error

abstract: We explore in some detail the notion of algorithmic stability as a viable framework for analyzing the generalization error of learning algorithms. We introduce the new notion of training stability of a learning algorithm and show that, in a general setting, it is sufficient for good bounds on generalization error. In the PAC setting, training stability is both necessary and sufficient for learnability.\ The approach based on training stability makes no reference to VC dimension or VC entropy. There is no need to prove uniform convergence, and generalization error is bounded directly via an extended McDiarmid inequality. As a result it potentially allows us to deal with a broader class of learning algorithms than Empirical Risk Minimization. \ We also explore the relationships among VC dimension, generalization error, and various notions of stability. Several examples of learning algorithms are considered.


title: Value Function Approximation in Zero-Sum Markov Games

abstract: This paper investigates value function approximation in the context of zero-sum Markov games, which can be viewed as a generalization of the Markov decision process (MDP) framework to the two-agent case. We generalize error bounds from MDPs to Markov games and describe generalizations of reinforcement learning algorithms to Markov games. We present a generalization of the optimal stopping problem to a two-player simultaneous move Markov game. For this special problem, we provide stronger bounds and can guarantee convergence for LSTD and temporal difference learning with linear value function approximation. We demonstrate the viability of value function approximation for Markov games by using the Least squares policy iteration (LSPI) algorithm to learn good policies for a soccer domain and a flow control problem.


title: General Lower Bounds based on Computer Generated Higher Order Expansions

abstract: In this article we show the rough outline of a computer algorithm to generate lower bounds on the exponential function of (in principle) arbitrary precision. We implemented this to generate all necessary analytic terms for the Boltzmann machine partition function thus leading to lower bounds of any order. It turns out that the extra variational parameters can be optimized analytically. We show that bounds upto nineth order are still reasonably calculable in practical situations. The generated terms can also be used as extra correction terms (beyond TAP) in mean field expansions.


title: Monitoring a Complez Physical System using a Hybrid Dynamic Bayes Net

abstract: The Reverse Water Gas Shift system (RWGS) is a complex physical system designed to produce oxygen from the carbon dioxide atmosphere on Mars. If sent to Mars, it would operate without human supervision, thus requiring a reliable automated system for monitoring and control. The RWGS presents many challenges typical of real-world systems, including: noisy and biased sensors, nonlinear behavior, effects that are manifested over different time granularities, and unobservability of many important quantities. In this paper we model the RWGS using a hybrid (discrete/continuous) Dynamic Bayesian Network (DBN), where the state at each time slice contains 33 discrete and 184 continuous variables. We show how the system state can be tracked using probabilistic inference over the model. We discuss how to deal with the various challenges presented by the RWGS, providing a suite of techniques that are likely to be useful in a wide range of applications. In particular, we describe a general framework for dealing with nonlinear behavior using numerical integration techniques, extending the successful Unscented Filter. We also show how to use a fixed-point computation to deal with effects that develop at different time scales, specifically rapid changes occurring during slowly changing processes. We test our model using real data collected from the RWGS, demonstrating the feasibility of hybrid DBNs for monitoring complex real-world physical systems.


title: Polynomial Value Iteration Algorithms for Detrerminstic MDPs

abstract: Value iteration is a commonly used and empirically competitive method in solving many Markov decision process problems. However, it is known that value iteration has only pseudo-polynomial complexity in general. We establish a somewhat surprising polynomial bound for value iteration on deterministic Markov decision (DMDP) problems. We show that the basic value iteration procedure converges to the highest average reward cycle on a DMDP problem in heta(n^2) iterations, or heta(mn^2) total time, where n denotes the number of states, and m the number of edges. We give two extensions of value iteration that solve the DMDP in heta(mn) time. We explore the analysis of policy iteration algorithms and report on an empirical study of value iteration showing that its convergence is much faster on random sparse graphs.


title: Decayed MCMC Filtering

abstract: Filtering---estimating the state of a partially observable Markov process from a sequence of observations---is one of the most widely studied problems in control theory, AI, and computational statistics. Exact computation of the posterior distribution is generally intractable for large discrete systems and for nonlinear continuous systems, so a good deal of effort has gone into developing robust approximation algorithms. This paper describes a simple stochastic approximation algorithm for filtering called {em decayed MCMC}. The algorithm applies Markov chain Monte Carlo sampling to the space of state trajectories using a proposal distribution that favours flips of more recent state variables. The formal analysis of the algorithm involves a generalization of standard coupling arguments for MCMC convergence. We prove that for any ergodic underlying Markov process, the convergence time of decayed MCMC with inverse-polynomial decay remains bounded as the length of the observation sequence grows. We show experimentally that decayed MCMC is at least competitive with other approximation algorithms such as particle filtering.


title: Formalizing Scenario Analysis

abstract: We propose a formal treatment of scenarios in the context of a dialectical argumentation formalism for qualitative reasoning about uncertain propositions. Our formalism extends prior work in which arguments for and against uncertain propositions were presented and compared in interaction spaces called Agoras. We now define the notion of a scenario in this framework and use it to define a set of qualitative uncertainty labels for propositions across a collection of scenarios. This work is intended to lead to a formal theory of scenarios and scenario analysis.