title: Critical solutions in topologically gauged N=8 CFTs in three dimensions

abstract: In this paper we discuss some special (critical) background solutions that arise in topological gauged ${\mathcal N}=8$ three-dimensional CFTs with SO(N) gauge group. These solutions solve the TMG equations (containing the parameters $\mu$ and $l$) for a certain set of values of $\mu l$ obtained by varying the number of scalar fields with a VEV. Apart from Minkowski, chiral round $AdS_3$ and null-warped $AdS_3$ (or Schr\"odinger(z=2)) we identify also a more exotic solution recently found in $TMG$ by Ertl, Grumiller and Johansson. We also discuss the spectrum, symmetry breaking pattern and the supermultiplet structure in the various backgrounds and argue that some properties are due to their common origin in a conformal phase. Some of the scalar fields, including all higgsed ones, turn out to satisfy three-dimensional singleton field equations. Finally, we note that topologically gauged ${\mathcal N}=6$ ABJ(M) theories have a similar, but more restricted, set of background solutions.


title: Sobolev and isoperimetric inequalities for submanifolds in weighted ambient spaces

abstract: In this paper, we prove a Sobolev and isoperimetric inequalities for submanifold in weighted manifold. Our results generalize the Hoffman-Spruck's inequalities.


title: Algorithms for Large-scale Whole Genome Association Analysis

abstract: In order to associate complex traits with genetic polymorphisms, genome-wide association studies process huge datasets involving tens of thousands of individuals genotyped for millions of polymorphisms. When handling these datasets, which exceed the main memory of contemporary computers, one faces two distinct challenges: 1) Millions of polymorphisms come at the cost of hundreds of Gigabytes of genotype data, which can only be kept in secondary storage; 2) the relatedness of the test population is represented by a covariance matrix, which, for large populations, can only fit in the combined main memory of a distributed architecture. In this paper, we present solutions for both challenges: The genotype data is streamed from and to secondary storage using a double buffering technique, while the covariance matrix is kept across the main memory of a distributed memory system. We show that these methods sustain high-performance and allow the analysis of enormous dataset


title: Colored bosons on top FBA and angular cross section for $t \bar t$ production

abstract: With full data set that corresponds to an integrated luminosity of 9.4 fb$^{-1}$, CDF has updated the top quark forward-backward asymmetry (FBA) as functions of rapidity difference $|\Delta y|$ and $t\bar t$ invariant mass $M_{t\bar t}$. Beside the sustained inconsistency between experiments and standard model (SM) predictions at large $|\Delta y|$ and $M_{t\bar t}$, an unexpected large first Legendre moment with $a_1= 0.39\pm 0.108 $ is found. In order to solve the large top FBA, we study the contributions of color triplet scalar and color octet vector boson. We find that the top FBA at $|\Delta y| >1$ and $M_{t\bar t} > 450$ GeV in triplet and octet model could be enhanced to be around 30% and 20%, whereas the first Legendre moment is $a^{\bf Di}_1= 0.38$ and $a^{\bf Axi}_1= 0.23$, respectively.


title: The parabolic Anderson model in a dynamic random environment: space-time ergodicity for the quenched Lyapunov exponent

abstract: We continue our study of the parabolic Anderson equation $\partial u(x,t)/\partial t = \kappa\Delta u(x,t) + \xi(x,t)u(x,t)$, $x\in\Z^d$, $t\geq 0$, where $\kappa \in [0,\infty)$ is the diffusion constant, $\Delta$ is the discrete Laplacian, and $\xi$ plays the role of a \emph{dynamic random environment} that drives the equation. The initial condition $u(x,0)=u_0(x)$, $x\in\Z^d$, is taken to be non-negative and bounded. The solution of the parabolic Anderson equation describes the evolution of a field of particles performing independent simple random walks with binary branching: particles jump at rate $2d\kappa$, split into two at rate $\xi \vee 0$, and die at rate $(-\xi) \vee 0$. We assume that $\xi$ is stationary and ergodic under translations in space and time, is not constant and satisfies $\E(|\xi(0,0)|)<\infty$, where $\E$ denotes expectation w.r.t.\ $\xi$. Our main object of interest is the quenched Lyapunov exponent $\lambda_0 (\kappa) = \lim_{t\to\infty} \frac{1}{t}\log u(0,t)$. In earlier work we showed that under certain mild space-time mixing assumptions the limit exists $\xi$-a.s., is finite and continuous on $[0,\infty)$, is globally Lipschitz on $(0,\infty)$, is not Lipschitz at 0, and satisfies $\lambda_0(0) = \E(\xi(0,0))$ and $\lambda_0(\kappa) > \E(\xi(0,0))$ for $\kappa \in (0,\infty)$.In the present paper we show that $\lim_{\kappa\to\infty} \lambda_0(\kappa) =\E(\xi(0,0))$ under an additional space-time mixing condition on $\xi$. This result shows that the parabolic Anderson model exhibits space-time ergodicity in the limit of large diffusivity. This fact is interesting because there are choices of $\xi$ that fulfill our assumption for which the annealed Lyapunov exponent $\lambda_1(\kappa) = \lim_{t\to\infty} \frac{1}{t}\log \E(u(0,t))$ is infinite on $[0,\infty)$, a situation that is referred to as strongly catalytic behavior.


title: Elaboration of the alpha-Model Derived from the BCS Theory of Superconductivity

abstract: The single-band alpha-model of superconductivity [H. Padamsee, et al., J. Low Temp. Phys. {\bf 12}, 387 (1973)] is a popular model that was adapted from the single-band Bardeen-Cooper-Schrieffer (BCS) theory of superconductivity mainly to allow fits to electronic heat capacity versus temperature T data that deviate from the BCS prediction. The model assumes that the normalized superconducting order parameter Delta(T)/Delta(0) and therefore the normalized London penetration depth lambda_L(T)/lambda_L(0) are the same as in BCS theory, calculated using the BCS value alpha_BCS = 1.764 of alpha = Delta(0)/(kB Tc), where kB is Boltzmann's constant and Tc is the superconducting transition temperature. On the other hand, to calculate the electronic free energy, entropy, heat capacity and thermodynamic critical field versus T, the alpha-model takes alpha to be an adjustable parameter. Here we write the BCS equations and limiting behaviors for the superconducting state thermodynamic properties explicitly in terms of alpha, as needed for calculations within the alpha-model, and present plots of the results versus T and alpha that are compared with the respective BCS predictions. Mechanisms such as gap anisotropy and strong coupling that can cause deviations of the thermodynamics from the BCS predictions, especially the heat capacity jump at Tc, are considered. Extensions of the alpha-model that have appeared in the literature such as the two-band model are also discussed. Tables of values of Delta(T)/Delta(0) and lambda_L(T)/lambda_L(0) calculated from the BCS theory using alpha = alpha_BCS are provided, which are the same in the alpha-model by assumption. Tables of values of the entropy, heat capacity and thermodynamic critical field versus T for seven values of alpha including alpha_BCS are also presented.


title: Extrapolation-based implicit-explicit general linear methods

abstract: For many systems of differential equations modeling problems in science and engineering, there are natural splittings of the right hand side into two parts, one non-stiff or mildly stiff, and the other one stiff. For such systems implicit-explicit (IMEX) integration combines an explicit scheme for the non-stiff part with an implicit scheme for the stiff part. In a recent series of papers two of the authors (Sandu and Zhang) have developed IMEX GLMs, a family of implicit-explicit schemes based on general linear methods. It has been shown that, due to their high stage order, IMEX GLMs require no additional coupling order conditions, and are not marred by order reduction. This work develops a new extrapolation-based approach to construct practical IMEX GLM pairs of high order. We look for methods with large absolute stability region, assuming that the implicit part of the method is A- or L-stable. We provide examples of IMEX GLMs with optimal stability properties. Their application to a two dimensional test problem confirms the theoretical findings.


title: A Participatory Universe of J. A. Wheeler as an Intentional Correlate of Embodied Subjects and an Example of Purposiveness in Physics

abstract: This paper investigates the role of human subjectivity and its delimiters in articulating the universe in physics and cosmology. As a case study, we reflect upon the complex of ideas of the so called Participatory Universe by later J. A. Wheeler. The objective of the paper is to explicate the role of the human agency as a centre of disclosure and manifestation of the universe as well the as teleology of scientific representation of the world implied by the intrinsic purposiveness of human actions.


title: Sensitivity of measured fission yields on prompt-neutron corrections

abstract: The amount of emitted prompt neutrons from the fission fragments increases as a function of excitation energy. Yet it is not fully understood whether the increase in \nu(A) as a function of E_{n} is mass dependent. The share of excitation energies among the fragments is still under debate, but there are reasons to believe that the excess in neutron emission originates only from the heavy fragments, leaving \nu_{light}(A) almost unchanged. In this work we investigated the consequences of a mass-dependent increase in \nu(A) on the final mass and energy distributions. The assumptions on \nu(A) are essential when analysing measurements based on the 2E-technique. This choice showed to be significant on the measured observables. For example, the post-neutron emission mass yield distribution revealed changes up to 10-30%. The outcome of this work pinpoint the urgent need to determine \nu(A) experimentally, and in particular, how \nu(A) changes as a function of incident-neutron energy. Until then, many fission yields in the data libraries could be largely affected, since they were analysed based on another assumption on the neutron emission.


title: Local characterization of 1d topologically ordered states

abstract: We consider 1d Hamiltonian systems whose ground states display symmetry protected topological order. We show that ground states within the topological phase cannot be connected with each other through LOCC between a bipartition of the system. Our claim is demonstrated by analyzing the entanglement spectrum and Renyi entropies of different physical systems providing examples for symmetry protected topological phases. Specifically, we consider spin-1/2 Cluster-Ising model and a class of spin-1 models undergoing quantum phase transitions to the Haldane phase. Our results provide a probe for simmetry-protected topological order, that holds true even at the system's local scale. Therefore our analysis can serve as as local experimental test for topological order.


title: Binary Frequencies in a Sample of Globular Clusters. I. Methodology and Initial Results

abstract: Binary stars are thought to be a controlling factor in globular cluster evolution, since they can heat the environmental stars by converting their binding energy to kinetic energy during dynamical interactions. Through such interaction, the binaries determine the time until core collapse. To test predictions of this model, we have determined binary fractions for 35 clusters. Here we present our methodology with a representative globular cluster NGC 4590. We use HST archival ACS data in the F606W and F814W bands and apply PSF-fitting photometry to obtain high quality color-magnitude diagrams. We formulate the star superposition effect as a Poisson probability distribution function, with parameters optimized through Monte-Carlo simulations. A model-independent binary fraction of (6.2 +- 0.3)% is obtained by counting stars that extend to the red side of the residual color distribution after accounting for the photometric errors and the star superposition effect. A model-dependent binary fraction is obtained by constructing models with a known binary fraction and an assumed binary mass-ratio distribution function. This leads to a binary fraction range of 6.8% to 10.8%, depending on the assumed shape to the binary mass ratio distribution, with the best fit occurring for a binary distribution that favors low mass ratios (and higher binary fractions). We also represent the method for radial analysis of the binary fraction in the representative case of NGC 6981, which shows a decreasing trend for the binary fraction towards the outside, consistent with theoretical predictions for the dynamical effect on the binary fraction.


title: Ultra Low Power Associative Computing with Spin Neurons and Resistive Crossbar Memory

abstract: Emerging resistive-crossbar memory (RCM) technology can be promising for computationally-expensive analog pattern-matching tasks. However, the use of CMOS analog-circuits with RCM would result in large power-consumption and poor scalability, thereby eschewing the benefits of RCM-based computation. We propose the use of low-voltage, fast-switching, magneto-metallic spin-neurons for ultra low-power non-Boolean computing with RCM. We present the design of analog associative memory for face recognition using RCM, where, substituting conventional analog circuits with spin-neurons can achieve ~100x lower power. This makes the proposed design ~1000x more energy-efficient than a 45nm-CMOS digital ASIC, thereby significantly enhancing the prospects of RCM based computational hardware.


title: In search of superluminal quantum communications: recent experiments and possible improvements

abstract: As shown in the famous \emph{EPR} paper (Einstein, Podolsky e Rosen,1935), Quantum Mechanics is non-local. The Bell theorem and the experiments by Aspect and many others, ruled out the possibility of explaining quantum correlations between entangled particles using local hidden variables models (except for implausible combinations of loopholes). Some authors (Bell, Eberhard, Bohm and Hiley) suggested that quantum correlations could be due to superluminal communications (tachyons) that propagate isotropically with velocity \emph{$v_{t}>c$} in a preferred reference frame. For finite values of \emph{$v_{t}$}, Quantum Mechanics and superluminal models lead to different predictions. Some years ago a Geneva group and our group did experiments on entangled photons to evidence possible discrepancies between experimental results and quantum predictions. Since no discrepancy was found, these experiments established only lower bounds for the possible tachyon velocities \emph{$v_{t}$}. Here we propose an improved experiment that should lead us to explore a much larger range of possible tachyon velocities \emph{$v_{t}$} for any possible direction of velocity $\vec{V}$ of the tachyons preferred frame.


title: The temperature evolution of the shear viscosity in a unitary Fermi gas

abstract: We present an ab initio determination of the shear viscosity for the unitary Fermi gas based on finite temperature quantum Monte Carlo (QMC) calculations and the Kubo linear-response formalism. The results are confronted with the bound for the shear viscosity originating from hydrodynamic fluctuations. Assuming smoothness of the frequency dependent shear viscosity eta(omega), we show that the bound is violated in the low temperature regime and the violation occurs simultaneously with the onset of the Cooper paring in the system. In order to preserve the hydrodynamic bound in QMC eta(omega) has to possess a sharp structure located in the vicinity of zero frequency which is not resolved by an analytic continuation procedure.


title: Ordered $*$-Semigroups and a $C^*$-Correspondence for a Partial Isometry

abstract: Certain $*$-semigroups are associated with the universal $C^*$-algebra generated by a partial isometry, which is itself the universal $C^*$-algebra of a $*$-semigroup. A fundamental role for a $*$-structure on a semigroup is emphasized, and ordered and matricially ordered $*$-semigroups are introduced, along with their universal $C^*$-algebras. The universal $C^*$-algebra generated by a partial isometry is isomorphic to a relative Cuntz-Pimsner $C^*$-algebra of a $C^*$-correspondence over the $C^*$-algebra of a matricially ordered $*$-semigroup. One may view the $C^*$-algebra of a partial isometry as the crossed product algebra associated with a dynamical system defined by a complete order map modelled by a partial isometry acting on a matricially ordered $*$-semigroup.


title: Cumulative Probability of Blast Fragmentation Effect

abstract: This paper presents formulae for calculation of cumulative probability of effect made by blast fragments. Analysis with Mott distribution, discrete fragment enumeration, spatial non-uniformity, numerical issues, and a generalisation for a set of effects are also discussed.


title: Wave-particle duality: an information-based approach

abstract: Recently, Bohr's complementarity principle was assessed in setups involving delayed choices. These works argued in favor of a reformulation of the aforementioned principle so as to account for situations in which a quantum system would simultaneously behave as wave and particle. Here we defend a framework that, supported by well-known experimental results and consistent with the decoherence paradigm, allows us to interpret complementarity in terms of correlations between the system and an {\em informer}. Our proposal offers formal definition and operational interpretation for the dual behavior in terms of both nonlocal resources and the couple work-information. Most importantly, our results provide a generalized information-based trade-off for the wave-particle duality and a causal interpretation for delayed-choice experiments.


title: Trapping Electromagnetic Solitons in Cylinders

abstract: Electromagnetic waves, in vacuum or dielectrics, can be confined in unbounded cylinders in such a way that they turn around the main axis. For particular choices of the cylinder's section, interesting stationary configurations may be assumed. By refining some results obtained in previous papers, additional more complex situations are examined here. For such peculiar guided waves an explicit expression is given in terms of Bessel's functions. Possible applications are in the development of whispering gallery resonators.


title: Near Heisenberg limited atomic clocks in the presence of decoherence

abstract: The ultimate stability of atomic clocks is limited by the quantum noise of the atoms. To reduce this noise it has been suggested to use entangled atomic ensembles with reduced atomic noise. Potentially this can push the stability all the way to the limit allowed by the Heisenberg uncertainty relation, which is denoted the Heisenberg limit. In practice, however, entangled states are often more prone to decoherence, which may prevent reaching this performance. Here we present an adaptive measurement protocol that in the presence of a realistic source of decoherence enables us to get near Heisenberg limited stability of atomic clocks using entangled atoms. The protocol may thus realize the full potential of entanglement for quantum metrology despite the detrimental influence of decoherence.


title: The effect of phenotypic selection on stochastic gene expression

abstract: Genetically identical cells in the same population can take on phenotypically variable states, leading to differentiated responses to external signals, such as nutrients and drug-induced stress. Many models and experiments have focused on a description based on discrete phenotypic states. Here we consider the effects of selection acting on a single trait, which we explicitly link to the variable number of proteins expressed by a gene. Considering different regulatory models for the gene under selection, we calculate the steady-state distribution of expression levels and show how the population adapts its expression to enhance its fitness. We quantitatively relate the overall fitness of the population to the heritability of expression levels, and their diversity within the population. We show how selection can increase or decrease the variability in the population, alter the stability of bimodal states, and impact the switching rates between metastable attractors.


title: Macro Lambda Calculus

abstract: The goal of our Macro Lambda Calculus project (MLC) is to encode lambda terms into interaction nets. Its software implementation will accept input in the notation similar to lambda calculus allowing macro definitions. Output is similar to interaction calculus and is suitable for our Interaction Nets Compiler program (INC). In this paper, we describe the interaction system for call-by-need evaluation and the mechanism of encoding lambda terms into this system which MLC is based on.


title: Gauge invariant backreaction in general single field models of inflation

abstract: In a general single field inflationary model we consider the effects of long wavelength scalar fluctuations on the effective expansion rate and equation of state seen by a class of free falling observers, using a physical gauge invariant formulation. In a previous work we showed that for a free massive inflaton no backreaction is observed within some constraints. In this paper we extend the validity of our previous results to the case of an arbitrary self-interacting inflation potential, working to second order in cosmological perturbation theory and to all order in slow-roll approximation. For these general inflationary models, we also show the equivalence of the free falling observers to the ones comoving with the inflaton field.


title: Phantom Instability of Viscous Dark Energy in Anisotropic Space-Time

abstract: Phantom dark energy is a proposal that explains the current observations that mildly favor the equation of state of dark energy $\omega^{de}$ crossing -1 at 68% confidence level. However, phantom fields are generally ruled out by ultraviolet quantum instabilities. To overcome this discrepancy, in this paper we propose a mechanism to show that how the presence of bulk viscosity in the cosmic fluid can temporarily drive the fluid into the phantom region ($\omega<-1$). As time is going on, phantom decays and ultimately $\omega^{de}$ approaches to -1. Then we show these quintessence and phantom descriptions of non-viscous and viscous dark energy and reconstruct the potential of these two scalar fields. Also a diagnostic for these models are performed by using the statefinder pairs $\{s,r\}$. All results are obtained in an anisotropic space-time which is a generalization of FLRW universe.


title: A competing risks approach for nonparametric estimation of transition probabilities in a non-Markov illness-death model

abstract: Competing risks model time to first event and type of first event. An example from hospital epidemiology is the incidence of hospital-acquired infection, which has to account for hospital discharge of non-infected patients as a competing risk. An illness-death model would allow to further study hospital outcomes of infected patients. Such a model typically relies on a Markov assumption. However, it is conceivable that the future course of an infected patient does not only depend on the time since hospital admission and current infection status but also on the time since infection. We demonstrate how a modified competing risks model can be used for nonparametric estimation of transition probabilities when the Markov assumption is violated.


title: Prediction of the dynamic oscillation threshold in a clarinet model with a linearly increasing blowing pressure : influence of noise

abstract: This paper presents an analysis of the effects of noise and precision on a simplified model of the clarinet driven by a variable control parameter. When the control parameter is varied the clarinet model undergoes a dynamic bifurcation. A consequence of this is the phenomenon of bifurcation delay: the bifurcation point is shifted from the static oscillation threshold to an higher value called dynamic oscillation threshold. In a previous work [8], the dynamic oscillation threshold is obtained analytically. In the present article, the sensitivity of the dynamic threshold on precision is analyzed as a stochastic variable introduced in the model. A new theoretical expression is given for the dynamic thresholds in presence of the stochastic variable, providing a fair prediction of the thresholds found in finite-precision simulations. These dynamic thresholds are found to depend on the increase rate and are independent on the initial value of the parameter, both in simulations and in theory.


title: The finiteness problem for automaton semigroups is undecidable

abstract: The finiteness problem for automaton groups and semigroups has been widely studied, several partial positive results are known. However we prove that, in the most general case, the problem is undecidable. We study the case of automaton semigroups. Given a NW-deterministic Wang tile set, we construct an Mealy automaton, such that the plane admit a valid Wang tiling if and only if the Mealy automaton generates a finite semigroup. The construction is similar to a construction by Kari for proving that the nilpotency problem for cellular automata is unsolvable. Moreover Kari proves that the tiling of the plane is undecidable for NW-deterministic Wang tile set. It follows that the finiteness problem for automaton semigroup is undecidable.


title: A fourth-order model for MEMS with clamped boundary conditions

abstract: The dynamical and stationary behaviors of a fourth-order equation in the unit ball with clamped boundary conditions and a singular reaction term are investigated. The equation arises in the modeling of microelectromechanical systems (MEMS) and includes a positive voltage parameter $\lambda$. It is shown that there is a threshold value $\lambda_*>0$ of the voltage parameter such that no radially symmetric stationary solution exists for $\lambda>\lambda_*$, while at least two such solutions exist for $\lambda\in (0,\lambda_*)$. Local and global well-posedness results are obtained for the corresponding hyperbolic and parabolic evolution problems as well as the occurrence of finite time singularities when $\lambda>\lambda_*$.


title: A solution to the Pompeiu problem

abstract: Let $f \in L_{loc}^1 (\R^n)\cap \mathcal{S}$, where $\mathcal{S}$ is the Schwartz class of distributions, and $$\int_{\sigma (D)} f(x) dx = 0 \quad \forall \sigma \in G, \qquad (*)$$ where $D\subset \R^n$ is a bounded domain, the closure $\bar{D}$ of which is diffeomorphic to a closed ball, and $S$ is its boundary. Then the comp$ is connected and path connected. By $G$ the group of all rigid motions of $\R^n$ is denoted. This group consists of all translations and rotations. A proof of the following theorem is given. Theorem 1. {\it Assume that $n=2$, $f\not\equiv 0$, and (*) holds. Then $D$ is a ball.} Corollary. {\it If the problem $(\nabla^2+k^2)u=0$ in $D$, $u_N|_S=0$, $u|_S=const\neq 0$ has a solution, then $D$ is a ball.} Here $N$ is the outer unit normal to $S$.


title: A discreteness criterion for groups containing parabolic isometries

abstract: This note will prove a discreteness criterion for groups of orientation-preserving isometries of the hyperbolic space which contain a parabolic element. It can be viewed as a generalization of the well-known results of Shimizu-Leutbecher and Jorgensen in dimensions 2 and 3, and is closely related to Waterman's inequality in higher dimensions. Unlike his algebraic method, the argument presented here is geometric and yields an improved asymptotic bound.


title: Fossil Evidence for the Two-phase Formation of Elliptical Galaxies

abstract: Massive early-type galaxies have undergone dramatic structural evolution over the last 10 Gyr. A companion paper shows that nearby elliptical galaxies with M*>1.3x10^{11} M_sun generically contain three photometric subcomponents: a compact inner component with effective radius Re<1 kpc, an intermediate-scale middle component with Re~2.5 kpc, and an extended outer envelope with $R_e \approx 10$ kpc. Here we attempt to relate these substructures with the properties of early-type galaxies observed at higher redshifts. We find that a hypothetical structure formed from combining the inner plus the middle components of local ellipticals follows a strikingly tight stellar mass-size relation, one that resembles the distribution of early-type galaxies at z~1.5. Outside of the central kpc, the median stellar mass surface density profiles of this composite structure agree closest with those of massive galaxies that have similar cumulative number density at 1.5<z<2.0 within the uncertainty. We propose that the central substructures in nearby ellipticals are the evolutionary descendants of the "red nuggets" formed under highly dissipative ("wet") conditions at high redshifts, as envisioned in the initial stages of the two-phase formation scenario recently advocated for massive galaxies. Subsequent accretion, plausibly through dissipationless ("dry") minor mergers, builds the outer regions of the galaxy identified as the outer envelope in our decomposition. The large scatter exhibited by this component on the stellar mass-size plane testifies to the stochastic nature of the accretion events.


title: Incremental Computation of Pseudo-Inverse of Laplacian: Theory and Applications

abstract: A divide-and-conquer based approach for computing the Moore-Penrose pseudo-inverse of the combinatorial Laplacian matrix $(\bb L^+)$ of a simple, undirected graph is proposed. % The nature of the underlying sub-problems is studied in detail by means of an elegant interplay between $\bb L^+$ and the effective resistance distance $(\Omega)$. Closed forms are provided for a novel {\em two-stage} process that helps compute the pseudo-inverse incrementally. Analogous scalar forms are obtained for the converse case, that of structural regress, which entails the breaking up of a graph into disjoint components through successive edge deletions. The scalar forms in both cases, show absolute element-wise independence at all stages, thus suggesting potential parallelizability. Analytical and experimental results are presented for dynamic (time-evolving) graphs as well as large graphs in general (representing real-world networks). An order of magnitude reduction in computational time is achieved for dynamic graphs; while in the general case, our approach performs better in practice than the standard methods, even though the worst case theoretical complexities may remain the same: an important contribution with consequences to the study of online social networks.


title: Mechanisms of Multi-strain Coexistence in Host-Phage Systems with Nested Infection Networks

abstract: Bacteria and their viruses ("bacteriophages") coexist in natural environments forming complex infection networks. Recent empirical findings suggest that phage-bacteria infection networks often possess a nested structure such that there is a hierarchical relationship among who can infect whom. Here we consider how nested infection networks may affect phage and bacteria dynamics using a multi-type Lotka-Volterra framework with cross-infection. Analysis of similar models have, in the past, assumed simpler interaction structures as a first step towards tractability. We solve the proposed model, finding trade-off conditions on the life-history traits of both bacteria and viruses that allow coexistence in communities with nested infection networks. First, we find that bacterial growth rate should decrease with increasing defense against infection. Second, we find that the efficiency of viral infection should decrease with host range. Next, we establish a relationship between relative densities and the curvature of life history trade-offs. We compare and contrast the current findings to the "Kill-the-Winner" model of multi-species phage-bacteria communities. Finally, we discuss a suite of testable hypotheses stemming from the current model concerning relationships between infection range, life history traits and coexistence in complex phage-bacteria communities.


title: ClusterCluster: Parallel Markov Chain Monte Carlo for Dirichlet Process Mixtures

abstract: The Dirichlet process (DP) is a fundamental mathematical tool for Bayesian nonparametric modeling, and is widely used in tasks such as density estimation, natural language processing, and time series modeling. Although MCMC inference methods for the DP often provide a gold standard in terms asymptotic accuracy, they can be computationally expensive and are not obviously parallelizable. We propose a reparameterization of the Dirichlet process that induces conditional independencies between the atoms that form the random measure. This conditional independence enables many of the Markov chain transition operators for DP inference to be simulated in parallel across multiple cores. Applied to mixture modeling, our approach enables the Dirichlet process to simultaneously learn clusters that describe the data and superclusters that define the granularity of parallelization. Unlike previous approaches, our technique does not require alteration of the model and leaves the true posterior distribution invariant. It also naturally lends itself to a distributed software implementation in terms of Map-Reduce, which we test in cluster configurations of over 50 machines and 100 cores. We present experiments exploring the parallel efficiency and convergence properties of our approach on both synthetic and real-world data, including runs on 1MM data vectors in 256 dimensions.


title: Altered Maxwell equations in the length gauge

abstract: The length gauge uses a scalar potential to describe a laser field, thus treating it as a longitudinal field rather than as a transverse field. This distinction is revealed in the fact that the Maxwell equations that relate to the length gauge are not the same as those for transverse fields. In particular, a source term is necessary in the length-gauge Maxwell equations, whereas the Coulomb-gauge description of plane waves possesses the basic property of transverse fields that they propagate with no source terms at all. This difference is shown to be importantly consequential in some previously unremarked circumstances; and it explains why the G\"oppert-Mayer gauge transformation does not provide the security that might be expected of full gauge equivalence.


title: Interplay of three kinds of motion in the disk counterpart of type II spicules: up-flow, transversal and torsional motions

abstract: Recently, it was shown that the complex dynamical behaviour of spicules has to be interpreted as the result of simultaneous action of three kinds of motion: (1) field aligned flows, (2) swaying motions, and (3) torsional motions. We use high-quality observations from CRISP at the SST to investigate signs of these different kinetic modes in spicules on the disk. Earlier, rapid blue-shifted excursions (RBEs), short-lived absorption features in the blue wing of chromospheric spectral lines, were identified as the disk counterpart of type II spicules. Here we report the existence of similar absorption features in the red wing of the Ca II 8542 and Halpha lines: rapid red-shifted excursions (RREs). RREs are found over the whole solar disk and are located in the same regions as RBEs: in the vicinity of magnetic field concentrations. RREs have similar characteristics as RBEs: they have similar lengths, widths, lifetimes, and average Doppler velocity. The striking similarity of RREs to RBEs implies that RREs are a manifestation of the same physical phenomenon and that spicules harbour motions that can result in a net red-shift when observed on-disk. We find that RREs are less abundant than RBEs and we interpret the higher number of RBEs and the decreased imbalance towards the limb as an indication that field-aligned up-flows have a significant contribution to the net Dopplershift. Most RREs and RBEs are observed in isolation but we find many examples of parallel and touching RRE/RBE pairs which appear to be part of the same spicule. We interpret the existence of these RRE/RBE pairs as signs that torsional motion is an important characteristic of spicules. The fact that most RBEs and RREs are observed in isolation agrees with the idea that transversal swaying motion is another important kinetic mode. We find examples of transitions from RRE to RBE and vice versa.


title: Black-hole solution without curvature singularity

abstract: An exact solution of the vacuum Einstein field equations over a nonsimply-connected manifold is presented. This solution is spherically symmetric and has no curvature singularity. It can be considered to be a regularization of the Schwarzschild solution over a simply-connected manifold, which has a curvature singularity at the center. Spherically symmetric collapse of matter in R^4 may result in this nonsingular black-hole solution, if quantum-gravity effects allow for topology change near the center.


title: Nonthermal emission model of isolated X-ray pulsar RX J0420.0-5022

abstract: In the present paper an alternative theoretical interpretation to the generally assumed thermal emission models of the observed X-ray spectrum of isolated pulsar RX J0420.0-5022 is presented. It is well known that the distribution function of relativistic particles is one-dimensional at the pulsar surface. However, cyclotron instability causes an appearance of transverse momenta of relativistic electrons, which as a result, start to radiate in the synchrotron regime. On the basis of the Vlasov's kinetic equation we study the process of the quasi-linear diffusion (QLD) developed by means of the cyclotron instability. This mechanism provides generation of optical and X-ray emission on the light cylinder lengthscales. The analysis of the three archival XMM-Newton observations of RX J0420.0-5022 is performed. Considering a different approach of the synchrotron emission theory, the spectral energy distribution is obtained that is in a good agreement with the observational data. A fit to the X-ray spectrum is performed using both the present synchrotron emission model spectrum absorbed by cold interstellar matter and generally assumed absorbed black-body model.


title: Elliptic Flow of Protons and Antiprotons in Au+Au Collisions at $\sqrt{s_{NN}}$ = 7.7--62.4 GeV within Alternative Scenarios of Three-Fluid Dynamics

abstract: Analysis of elliptic flow of protons and antiprotons in Au+Au collisions is performed in a wide range of incident energies $\sqrt{s_{NN}}$ = 7.7--62.4 GeV. Simulations has been done within the three-fluid model employing a purely hadronic equation of state (EoS) and two versions of the EoS involving deconfinement transition: an EoS with the first-order phase transition and that with a smooth crossover transition. It is found that the proton data are reproduced approximately to the same extent within all of the scenarios, including the hadronic one, while the deconfinement scenarios look certainly preferable for the antiproton elliptic flow. The fact that difference between elliptic flows of protons and antiprotons decreases with the incident energy rise is a consequence of reducing baryon stopping rather than an onset of deconfinement.


title: Hunting for GeV emission from the gamma-ray binary HESS J0632+057

abstract: In the last decade Cherenkov telescopes on the ground and space-based gamma-ray instruments have identified a new sub-class of high mass X-ray binaries (HMXB), whose emission is dominated by gamma rays. To date only five of these systems have been definitively identified. However at GeV energies there is still one, HESS J0632+057, that has no reported detection by the Fermi LAT. A deep search for gamma-ray emission of HESS J0632+057 has been performed using more than 3.5 years of Fermi-LAT data. We present the results of this search in the context of the other known gamma-ray binary systems.


title: Dilepton Spectroscopy of QCD Matter at Collider Energies

abstract: Low-mass dilepton spectra as measured in high-energy heavy-ion collisions are a unique tool to obtain spectroscopic information about the strongly interacting medium produced in these reactions. Specifically, in-medium modifications of the vector spectral function, which is well known in the vacuum, can be deduced from the thermal radiation off the expanding QCD fireball. This, in particular, allows to investigate the fate of the $\rho$ resonance in the dense medium, and possibly infer from it signatures of the (partial) restoration of chiral symmetry, which is spontaneously broken in the QCD vacuum. After briefly reviewing calculations of thermal dilepton emission rates from hot QCD matter, utilizing effective hadronic theory, lattice QCD or resummed perturbative QCD, we focus on applications to dilepton spectra at heavy-ion collider experiments at RHIC and LHC. This includes invariant-mass spectra at full RHIC energy with transverse-momentum dependencies and azimuthal asymmetries, as well as a systematic investigation of the excitation function down to fixed-target energies, thus making contact to previous precision measurements at the SPS. Furthermore, predictions for the energy frontier at the LHC are presented in both dielectron and dimuon channels.


title: Embedding of Blink Frequency in Electrooculography Signal using Difference Expansion based Reversible Watermarking Technique

abstract: In the past few years, like other fields, rapid expansion of digitization and globalization has influenced the medical field as well. For progress of diagnostic results most of the reputed hospitals and diagnostic centres all over the world have started exchanging medical information. In this proposed method, the calculated diagnostic parametric values of the original Electrooculography (EOG) signal are embedded as a watermark by using Difference Expansion (DE) algorithm based reversible watermarking technique. The extracted watermark provides the required parametric values at the recipient end without any post computation of the recovered EOG signal. By computing the parametric values from the recovered signal, the integrity of the extracted watermark can be validated. The time domain features of EOG signal are calculated for the generation of watermark. In the current work, various features are studied and two major features related to blink frequency are used to generate the watermark. The high Signal to Noise Ratio (SNR) and the Bit Error Rate (BER) claim the robustness of the proposed method.


title: Nonlinear Charged Black Holes in Anti-de Sitter Quasi-Topological Gravity

abstract: In this paper, we present the static charged solutions of quartic quasitopological gravity in the presence of a nonlinear electromagnetic field. Two branches of these solutions present black holes with one or two horizons or a naked singularity depending on the charge and mass of the black hole. The entropy of the charged black holes of fourth order quasitopological gravity through the use of Wald formula is computed and the mass, temperature and the charge of these black holes are found as well. We show that black holes with spherical, flat and hyperbolical horizon in quasitopological gravity are stable for any allowed quasitopological parameters.


title: Parametrization of white matter manifold-like structures using principal surfaces

abstract: In this manuscript, we are concerned with data generated from a diffusion tensor imaging (DTI) experiment. The goal is to parameterize manifold-like white matter tracts, such as the corpus callosum, using principal surfaces. We approach the problem by finding a geometrically motivated surface-based representation of the corpus callosum and visualize the fractional anisotropy (FA) values projected onto the surface; the method applies to any other diffusion summary as well as to other white matter tracts. We provide an algorithm that 1) constructs the principal surface of a corpus callosum; 2) flattens the surface into a parametric 2D map; 3) projects associated FA values on the map. The algorithm was applied to a longitudinal study containing 466 diffusion tensor images of 176 multiple sclerosis (MS) patients observed at multiple visits. For each subject and visit the study contains a registered DTI scan of the corpus callosum at roughly 20,000 voxels. Extensive simulation studies demonstrate fast convergence and robust performance of the algorithm under a variety of challenging scenarios.


title: On Differentially Private Filtering for Event Streams

abstract: Rigorous privacy mechanisms that can cope with dynamic data are required to encourage a wider adoption of large-scale monitoring and decision systems relying on end-user information. A promising approach to develop these mechanisms is to specify quantitative privacy requirements at design time rather than as an afterthought, and to rely on signal processing techniques to achieve satisfying trade-offs between privacy and performance specifications. This paper discusses, from the signal processing point of view, an event stream analysis problem introduced in the database and cryptography literature. A discrete-valued input signal describes the occurrence of events contributed by end-users, and a system is supposed to provide some output signal based on this information, while preserving the privacy of the participants. The notion of privacy adopted here is that of event-level differential privacy, which provides strong privacy guarantees and has important operational advantages. Several mechanisms are described to provide differentially private output signals while minimizing the impact on performance. These mechanisms demonstrate the benefits of leveraging system theoretic techniques to provide privacy guarantees for dynamic systems.


title: A Modification of LLR

abstract: The Lucas-Lehmer (LL) primality test for Mersenne numbers is the fastest known primality test. In 1969, Hans Riesel published a modification of LL to test numbers of the form $N = h \cdot 2^n - 1$, where $h < 2^n$ is an odd integer and $n \ge 2$ \cite{Riesel}. This test is now known as the Lucas-Lehmer-Riesel (LLR) primality test. In Algorithm \ref{PrimalityAlgorithm}, we present a modification of LLR which works for any odd integer $N$. A probabilistic version of our algorithm runs in expected time $\tilde{O}(\log^3 N)$, and a deterministic version in expected $\tilde{O}(\log^4 N)$. We conclude with a conjecture which, if true, would imply that there exists a polynomial time algorithm for factoring integers.


title: Generalized Elliptic Integrals and Applications

abstract: We use some general properties, presented in previous work, to evaluate special cases of integrals relating Rogers-Ramanujan continued fraction, eta function and elliptic integrals.


title: Experimental Neutron-Induced Fission Fragment Mass Yields of 232Th and 238U at Energies from 10 to 33 MeV

abstract: Development of nuclear energy applications requires data for neutron-induced reactions for actinides in a wide neutron energy range. Here we describe measurements of pre-neutron emission fission fragment mass yields of 232Th and 238U at incident neutron energies from 10 to 33 MeV. The measurements were done at the quasi-monoenergetic neutron beam of the Louvain-la-Neuve cyclotron facility CYCLONE; a multi-section twin Frisch-gridded ionization chamber was used to detect fission fragments. For the peak neutron energies at 33, 45 and 60 MeV, the details of the data analysis and the experimental results have been published before and in this work we present data analysis in the low-energy tail of the neutron energy spectra. The preliminary measurement results are compared with available experimental data and theoretical predictions.


title: Local Radiation-Driven Instabilities in Post-Main Sequence Massive Stars

abstract: Late in their evolution, massive stars may undergo periods of violent instability and mass loss, but the mechanism responsible for these episodes has not been identified. We study one potential contributor: the development of local radiation-driven instabilities in the outer layers of main sequence (MS) and post-MS massive stars. We construct a sequence of massive stellar evolution models and investigate where they are subject to local radiative instabilities, both in the presence of magnetic fields and without them,and at a range of metallicities. We find that these types of instabilities do not occur in solar-metallicity MS stars up to 100\,M$_\odot$, but they set in immediately post-MS for stars heavier than $\sim 25$\,M$_\odot$. Once an instability appears, it involves a significant amount of mass in the star's upper layers (up to $\sim 1$ per cent of the initial stellar mass), suggesting that radiation-driven instabilities are a potentially viable mechanism for dynamic mass loss. We find that the presence of magnetic fields at strengths low enough not to disturb the hydrostatic balance of the star does not alter these results. Stars with sub-Solar metallicity also show instability, but their instabilities involve less mass and appear later in the star's evolution.


title: Examples of backreaction of small scale inhomogeneities in cosmology

abstract: In previous work, we introduced a new framework to treat large scale backreaction effects due to small scale inhomogeneities in general relativity. We considered one-parameter families of spacetimes for which such backreaction effects can occur, and we proved that, provided the weak energy condition on matter is satisfied, the leading effect of small scale inhomogeneities on large scale dynamics is to produce a traceless effective stress-energy tensor that itself satisfies the weak energy condition. In this work, we illustrate the nature of our framework by providing two explicit examples of one-parameter families with backreaction. The first, based on previous work of Berger, is a family of polarized vacuum Gowdy spacetimes on a torus, which satisfies all of the assumptions of our framework. As the parameter approaches its limiting value, the metric uniformly approaches a smooth background metric, but spacetime derivatives of the deviation of the metric from the background metric do not converge uniformly to zero. The limiting metric has nontrivial backreaction from the small scale inhomogeneities, with an effective stress-energy that is traceless and satisfies the weak energy condition, in accord with our theorems. Our second one-parameter family consists of metrics which have a uniform Friedmann-Lemaitre-Robertson-Walker limit. This family satisfies all of our assumptions with the exception of the weak energy condition for matter. In this case, the limiting metric has an effective stress-energy tensor which is not traceless. We emphasize the importance of imposing energy conditions on matter in studies of backreaction.


title: Vortex matter in low dimensional systems with proximity induced superconductivity

abstract: We study theoretically the vortex matter structure in low dimensional (LD) systems with superconducting order induced by proximity to a bulk superconductor. We analyze the effects of microscopic coupling mechanisms between the two systems and the effects of possible mismatch in the band structures of these materials on the energy spectrum of vortex-core electrons. The unusual structure of vortex cores is discussed in the context of recent tunneling microscopy/spectroscopy experiments.


title: One-Dimensional Navier-Stokes Finite Element Flow Model

abstract: This technical report documents the theoretical, computational, and practical aspects of the one-dimensional Navier-Stokes finite element flow model. The document is particularly useful to those who are interested in implementing, validating and utilizing this relatively-simple and widely-used model.


title: Dark energy properties from large future galaxy surveys

abstract: We perform a detailed forecast on how well a {\sc Euclid}-like survey will be able to constrain dark energy and neutrino parameters from a combination of its cosmic shear power spectrum, galaxy power spectrum, and cluster mass function measurements. We find that the combination of these three probes vastly improves the survey's potential to measure the time evolution of dark energy. In terms of a dark energy figure-of-merit defined as $(\sigma(w_{\mathrm p}) \sigma(w_a))^{-1}$, we find a value of 690 for {\sc Euclid}-like data combined with {\sc Planck}-like measurements of the cosmic microwave background (CMB) anisotropies in a 10-dimensional cosmological parameter space, assuming a $\Lambda$CDM fiducial cosmology. For the more commonly used 7-parameter model, we find a figure-of-merit of 1900 for the same data combination. We consider also the survey's potential to measure dark energy perturbations in models wherein the dark energy is parameterised as a fluid with a nonstandard non-adiabatic sound speed, and find that in an \emph{optimistic} scenario in which $w_0$ deviates by as much as is currently observationally allowed from $-1$, models with $\hat{c}_\mathrm{s}^2 = 10^{-6}$ and $\hat{c}_\mathrm{s}^2 = 1$ can be distinguished at more than $2\sigma$ significance. We emphasise that constraints on the dark energy sound speed from cluster measurements are strongly dependent on the modelling of the cluster mass function; significantly weaker sensitivities ensue if we modify our model to include fewer features of nonlinear dark energy clustering. Finally, we find that the sum of neutrino masses can be measured with a $1 \sigma$ precision of 0.015~eV, (abridged)


title: Coherent Josephson qubit suitable for scalable quantum integrated circuits

abstract: We demonstrate a planar, tunable superconducting qubit with energy relaxation times up to 44 microseconds. This is achieved by using a geometry designed to both minimize radiative loss and reduce coupling to materials-related defects. At these levels of coherence, we find a fine structure in the qubit energy lifetime as a function of frequency, indicating the presence of a sparse population of incoherent, weakly coupled two-level defects. This is supported by a model analysis as well as experimental variations in the geometry. Our `Xmon' qubit combines facile fabrication, straightforward connectivity, fast control, and long coherence, opening a viable route to constructing a chip-based quantum computer.


title: Determinacy and J\'onsson cardinals in $L(\mathbb{R})$

abstract: Assume $\mathsf{ZF}+\mathsf{AD}+V=L(\mathbb{R})$ and let $\kappa<\Theta$ be an uncountable cardinal. We show that $\kappa$ is J\'onsson, and that if $\mathrm{cof}(\kappa)=\omega$ then $\kappa$ is Rowbottom. We also establish some other partition properties.


title: Phase Space Formulation of Population Dynamics in Ecology

abstract: A phase space theory for population dynamics in Ecology is presented. This theory applies for a certain class of dynamical systems, that will be called M-systems, for which a conserved quantity, the M-function, can be defined in phase space. This M-function is the generator of time displacements and contains all the dynamical information of the system. In this sense the M-function plays the role of the hamiltonian function for mechanical systems. In analogy with Hamilton theory we derive equations of motion as derivatives over the resource function in phase space. A M-bracket is defined which allows one to perform a geometrical approach in analogy to Poisson bracket of hamiltonian systems. We show that the equations of motion can be derived from a variational principle over a functional J of the trajectories. This functional plays for M-systems the same role than the action S for hamiltonian systems. Finally, three important systems in population dynamics, namely, Lotka-Volterra, self-feeding and logistic evolution, are shown to be M-systems.


title: Terahertz Conductivity of the Heavy-Fermion State in CeCoIn5

abstract: The optical properties of thin films of the heavy-fermion compound CeCoIn5, which were deposited by molecular beam epitaxy onto MgF2 substrates, have been studied at frequencies 7 to 45 cm^{-1} (corresponding to 0.2 to 1.3 THz) and temperatures 2 to 300 K. We observe an electrodynamic behavior which is typical for heavy fermions, namely Drude-like conductivity with a relaxation rate at rather low frequencies. This relaxation rate increases almost linearly with temperature up to at least 30 K. The coherent heavy-fermion state, characterized by an increase of the effective mass, continuously evolves upon cooling and is not fully developed for temperatures as low as 5 K.


title: Sharing of Semantically Enhanced Information for the Adaptive Execution of Business Processes

abstract: Motivated from the Context Aware Computing, and more particularly from the Data-Driven Process Adaptation approach, we propose the Semantic Context Space (SCS) Engine which aims to facilitate the provision of adaptable business processes. The SCS Engine provides a space which stores semantically annotated data and it is open to other processes, systems, and external sources for information exchange. The specified implementation is inspired from the Semantic TupleSpace and uses the JavaSpace Service of the Jini Framework (changed to Apache River lately) as an underlying basis. The SCS Engine supplies an interface where a client can execute the following operations: (i) write: which inserts in the space available information along with its respective meta-information, (ii) read: which retrieves from the space information which meets specific meta-information constrains, and (iii) take: which retrieves and simultaneously deletes from the space information which meets specific meta-information constrains. In terms of this thesis the available types of meta-information are based on ontologies described in RDFS or WSML. The applicability of the SCS Engine implementation in the context of data-driven process adaptation has been ensured by an experimental evaluation of the provided operations. Ultimately, we refer to open issues which could be addressed to enrich the proposed Engine with additional features.


title: Controlling the position of traveling waves in reaction-diffusion systems

abstract: We present a method to control the position as a function of time of one-dimensional traveling wave solutions to reaction-diffusion systems according to a pre-specified protocol of motion. Given this protocol, the control function is found as the solution of a perturbatively derived integral equation. Two cases are considered. First, we derive an analytical expression for the space ($x$) and time ($t$) dependent control function $f\left(x,t\right)$ that is valid for arbitrary protocols and many reaction-diffusion systems. These results are close to numerically computed optimal controls. Second, for stationary control of traveling waves in one-component systems, the integral equation reduces to a Fredholm integral equation of the first kind. In both cases, the control can be expressed in terms of the uncontrolled wave profile and its propagation velocity, rendering detailed knowledge of the reaction kinetics unnecessary.


title: Duality of Entanglement Norms

abstract: We consider four norms on tensor product spaces that have appeared in quantum information theory and demonstrate duality relationships between them. We show that the product numerical radius is dual to the robustness of entanglement, and we similarly show that the S(k)-norm is dual to the projective tensor norm. We show that, analogous to how the product numerical radius and the S(k)-norm characterize k-block positivity of operators, there is a natural version of the projective tensor norm that characterizes Schmidt number. In this way we obtain an elementary new proof of the cross norm criterion for separability, and we also generalize both the cross norm and realignment criteria to the case of arbitrary Schmidt number.


title: Socially optimal charging strategies for electric vehicles

abstract: Electric vehicles represent a promising technology for reducing emissions and dependence on fossil fuels and have started entering different automotive markets. In order to bolster their adoption by consumers and hence enhance their penetration rate, a charging station infrastructure needs to be deployed. This paper studies decentralized policies that assign electric vehicles to a network of charging stations with the goal to achieve little to no queueing. This objective is especially important for electric vehicles, whose charging times are fairly long. The social optimality of the proposed policies is established in the many-server regime, where each station is equipped with multiple charging slots. Further, convergence issues of the algorithm that achieves the optimal policy are examined. Finally, the results provide insight on how to address questions related to the optimal location deployment of the infrastructure.


title: Determining neutron star masses and radii using energy-resolved waveforms of X-ray burst oscillations

abstract: Simultaneous, precise measurements of the mass $M$ and radius $R$ of neutron stars can yield uniquely valuable information about the still uncertain properties of cold matter at several times the density of nuclear matter. One method that could be used to measure $M$ and $R$ is to analyze the energy-dependent waveforms of the X-ray flux oscillations seen during some thermonuclear bursts from some neutron stars. These oscillations are thought to be produced by X-ray emission from hotter regions on the surface of the star that are rotating at or near the spin frequency of the star. Here we explore how well $M$ and $R$ could be determined by generating, and analyzing using Bayesian techniques, synthetic energy-resolved X-ray data that we produce assuming a future space mission having 2--30 keV energy coverage and an effective area of 10 m$^2$, such as the proposed \textit{LOFT} or \textit{AXTAR} missions. We find that if the hot spot is within 10$^\circ$ of the rotation equator, both $M$ and $R$ can usually be determined with an uncertainty of about 10% if there are $10^6$ total counts from the spot, whereas waveforms from spots within 20$^\circ$ of the rotation pole provide no useful constraints. These constraints can usually be achieved even if the burst oscillations vary with time and data from multiple bursts must be used to obtain 10$^6$ counts from the hot spot. This is therefore a promising method to constrain $M$ and $R$ tightly enough to discriminate strongly between competing models of cold, high-density matter.


title: The PAV algorithm optimizes binary proper scoring rules

abstract: There has been much recent interest in application of the pool-adjacent-violators (PAV) algorithm for the purpose of calibrating the probabilistic outputs of automatic pattern recognition and machine learning algorithms. Special cost functions, known as proper scoring rules form natural objective functions to judge the goodness of such calibration. We show that for binary pattern classifiers, the non-parametric optimization of calibration, subject to a monotonicity constraint, can be solved by PAV and that this solution is optimal for all regular binary proper scoring rules. This extends previous results which were limited to convex binary proper scoring rules. We further show that this result holds not only for calibration of probabilities, but also for calibration of log-likelihood-ratios, in which case optimality holds independently of the prior probabilities of the pattern classes.


title: Semiclassical description of collapses and revivals of quantum wave packets in bounded domains

abstract: We study a special kind of semiclassical limit of quantum dynamics on a circle and in a box (infinite potential well with hard walls) as the Planck constant tends to zero and time tends to infinity. The results give detailed information about all stages of evolution of quantum wave packets: semiclassical motion, collapses, revivals, as well as intermediate stages. In particular, we rigorously justify the fact that the spatial distribution of a wave packet is most of the time close to uniform. This fact was previously known only from numerical calculations. We apply the obtained results to a problem of classical mechanics: deciding whether recently suggested functional classical mechanics is preferable to traditional Newtonian one from the quantum-mechanical point of view. To do this, we study the semiclassical limit of the Husimi functions of quantum states. We show that functional mechanics remains valid at larger time scales than Newtonian one and, therefore, is preferable. Finally, we analyse the quantum dynamics in a box in case when the size of the box is known with a random error. We show that, in this case, the probability distribution of the position of a quantum particle is not almost periodic, but tends to a limit distribution as time indefinitely increases.


title: A primer on information theory, with applications to neuroscience

abstract: Given the constant rise in quantity and quality of data obtained from neural systems on many scales ranging from molecular to systems', information-theoretic analyses became increasingly necessary during the past few decades in the neurosciences. Such analyses can provide deep insights into the functionality of such systems, as well as a rigid mathematical theory and quantitative measures of information processing in both healthy and diseased states of neural systems. This chapter will present a short introduction to the fundamentals of information theory, especially suited for people having a less firm background in mathematics and probability theory. To begin, the fundamentals of probability theory such as the notion of probability, probability distributions, and random variables will be reviewed. Then, the concepts of information and entropy (in the sense of Shannon), mutual information, and transfer entropy (sometimes also referred to as conditional mutual information) will be outlined. As these quantities cannot be computed exactly from measured data in practice, estimation techniques for information-theoretic quantities will be presented. The chapter will conclude with the applications of information theory in the field of neuroscience, including questions of possible medical applications and a short review of software packages that can be used for information-theoretic analyses of neural data.


title: Hybrid Textures of Majorana Neutrino Mass Matrix and Current Experimental Tests

abstract: Motivated by recent measurements of a relatively large \theta_{13} in the Daya Bay and RENO reactor neutrino experiments, we carry out a systematic analysis of the hybrid textures of Majorana neutrino mass matrix $M_\nu$, which contain one texture zero and two equal nonzero matrix elements. We show that three neutrino masses (m_1, m_2, m_3) and three leptonic CP-violating phases (\delta, \rho, \sigma) can fully be determined from two neutrino mass-squared differences (\delta m^2, \Delta m^2) and three flavor mixing angles (\theta_{12}, \theta_{23}, \theta_{13}). Out of sixty logically possible patterns of $M_\nu$, thirty-nine are found to be compatible with current experimental data at the 3\sigma level. We demonstrate that the texture zero of $M_\nu$ is stable against one-loop quantum corrections, while the equality between two independent elements not. Phenomenological implications of $M_\nu$ for the neutrinoless double-beta decay and leptonic CP violation are discussed, and a realization of the texture zero and equality by means of discrete flavor symmetries is illustrated.


title: On weights for relative motives with integral coefficients

abstract: The goal of this paper is to define a certain Chow weight structure for the category of Voevodsky's motivic complexes with integral coefficients (as described by Cisinski and Deglise) over any excellent finite-dimensional separated scheme $S$. Our results are parallel to (though substantially weaker than) the corresponding 'rational coefficient' statements proved by D. Hebert and the author. As an immediate consequence of the existence of 'weights', we obtain certain (Chow)-weight spectral sequences and filtrations for any (co)homology of $S$-motives.


title: One-shot lossy quantum data compression

abstract: We provide a framework for one-shot quantum rate distortion coding, in which the goal is to determine the minimum number of qubits required to compress quantum information as a function of the probability that the distortion incurred upon decompression exceeds some specified level. We obtain a one-shot characterization of the minimum qubit compression size for an entanglement-assisted quantum rate-distortion code in terms of the smooth max-information, a quantity previously employed in the one-shot quantum reverse Shannon theorem. Next, we show how this characterization converges to the known expression for the entanglement-assisted quantum rate distortion function for asymptotically many copies of a memoryless quantum information source. Finally, we give a tight, finite blocklength characterization for the entanglement-assisted minimum qubit compression size of a memoryless isotropic qubit source subject to an average symbol-wise distortion constraint.


title: On inferring extinction laws in z~6 quasars as signatures of supernova dust

abstract: Unusual extinction curves of high-redshift QSOs have been taken as evidence that dust is primarily produced by supernovae at high redshift. In particular, the 3000 A Todini-Ferrara-Maiolino kink in the extinction curve of the z = 6.20 SDSS J1048+4637 has been attributed to supernova dust. Here we discuss the challenges in inferring robust extinction curves of high-redshift QSOs and critically assess previous claims of detection of supernova dust. In particular, we address the sensitivity to the choice of intrinsic QSO spectrum, the need for a long wavelength baseline, and the drawbacks in fitting theoretical extinction curves. In a sample of 21 QSOs at z ~ 6 we detect significant ultraviolet extinction using existing broad-band optical, near-infrared, and Spitzer photometry. The median extinction curve is consistent with a Small Magellanic Cloud curve with A_1450 ~ 0.7 mag and does not exhibit any conspicuous (restframe) 2175 A or 3000 A features. For two QSOs, SDSS J1044-0125 at z = 5.78 and SDSS J1030+0524 at z = 6.31, we further present X-shooter spectra covering the wavelength range 0.9-2.5 um. The resulting non-parametric extinction curves do not exhibit the 3000 A kink. Finally, in a re-analysis of literature spectra of SDSS J1048+4637, we do not find evidence for a conspicuous kink. We conclude that the existing evidence for a 3000 A feature is weak and that the overall dust properties at high and low redshift show no significant differences. This, however, does not preclude supernovae from dominating the dust budget at high redshift.


title: An Almost-Linear-Time Algorithm for Approximate Max Flow in Undirected Graphs, and its Multicommodity Generalizations

abstract: In this paper, we introduce a new framework for approximately solving flow problems in capacitated, undirected graphs and apply it to provide asymptotically faster algorithms for the maximum $s$-$t$ flow and maximum concurrent multicommodity flow problems. For graphs with $n$ vertices and $m$ edges, it allows us to find an $\epsilon$-approximate maximum $s$-$t$ flow in time $O(m^{1+o(1)}\epsilon^{-2})$, improving on the previous best bound of $\tilde{O}(mn^{1/3} poly(1/\epsilon))$. Applying the same framework in the multicommodity setting solves a maximum concurrent multicommodity flow problem with $k$ commodities in $O(m^{1+o(1)}\epsilon^{-2}k^2)$ time, improving on the existing bound of $\tilde{O}(m^{4/3} poly(k,\epsilon^{-1})$. Our algorithms utilize several new technical tools that we believe may be of independent interest: - We give a non-Euclidean generalization of gradient descent and provide bounds on its performance. Using this, we show how to reduce approximate maximum flow and maximum concurrent flow to the efficient construction of oblivious routings with a low competitive ratio. - We define and provide an efficient construction of a new type of flow sparsifier. In addition to providing the standard properties of a cut sparsifier our construction allows for flows in the sparse graph to be routed (very efficiently) in the original graph with low congestion. - We give the first almost-linear-time construction of an $O(m^{o(1)})$-competitive oblivious routing scheme. No previous such algorithm ran in time better than $\tilde{{\Omega}}(mn)$. We also note that independently Jonah Sherman produced an almost linear time algorithm for maximum flow and we thank him for coordinating submissions.


title: The structure of Bayes nets for vision recognition

abstract: This paper is part of a study whose goal is to show the effciency of using Bayes networks to carry out model based vision calculations. [Binford et al. 1987] Recognition proceeds by drawing up a network model from the object's geometric and functional description that predicts the appearance of an object. Then this network is used to find the object within a photographic image. Many existing and proposed techniques for vision recognition resemble the uncertainty calculations of a Bayes net. In contrast, though, they lack a derivation from first principles, and tend to rely on arbitrary parameters that we hope to avoid by a network model. The connectedness of the network depends on what independence considerations can be identified in the vision problem. Greater independence leads to easier calculations, at the expense of the net's expressiveness. Once this trade-off is made and the structure of the network is determined, it should be possible to tailor a solution technique for it. This paper explores the use of a network with multiply connected paths, drawing on both techniques of belief networks [Pearl 86] and influence diagrams. We then demonstrate how one formulation of a multiply connected network can be solved.


title: Summary of A New Normative Theory of Probabilistic Logic

abstract: By probabilistic logic I mean a normative theory of belief that explains how a body of evidence affects one's degree of belief in a possible hypothesis. A new axiomatization of such a theory is presented which avoids a finite additivity axiom, yet which retains many useful inference rules. Many of the examples of this theory--its models do not use numerical probabilities. Put another way, this article gives sharper answers to the two questions: 1.What kinds of sets can used as the range of a probability function? 2.Under what conditions is the range set of a probability function isomorphic to the set of real numbers in the interval 10,1/ with the usual arithmetical operations?


title: Probability Distributions Over Possible Worlds

abstract: In Probabilistic Logic Nilsson uses the device of a probability distribution over a set of possible worlds to assign probabilities to the sentences of a logical language. In his paper Nilsson concentrated on inference and associated computational issues. This paper, on the other hand, examines the probabilistic semantics in more detail, particularly for the case of first-order languages, and attempts to explain some of the features and limitations of this form of probability logic. It is pointed out that the device of assigning probabilities to logical sentences has certain expressive limitations. In particular, statistical assertions are not easily expressed by such a device. This leads to certain difficulties with attempts to give probabilistic semantics to default reasoning using probabilities assigned to logical sentences.


title: Hierarchical Evidence and Belief Functions

abstract: Dempster/Shafer (D/S) theory has been advocated as a way of representing incompleteness of evidence in a system's knowledge base. Methods now exist for propagating beliefs through chains of inference. This paper discusses how rules with attached beliefs, a common representation for knowledge in automated reasoning systems, can be transformed into the joint belief functions required by propagation algorithms. A rule is taken as defining a conditional belief function on the consequent given the antecedents. It is demonstrated by example that different joint belief functions may be consistent with a given set of rules. Moreover, different representations of the same rules may yield different beliefs on the consequent hypotheses.


title: Decision-Theoretic Control of Problem Solving: Principles and Architecture

abstract: This paper presents an approach to the design of autonomous, real-time systems operating in uncertain environments. We address issues of problem solving and reflective control of reasoning under uncertainty in terms of two fundamental elements: l) a set of decision-theoretic models for selecting among alternative problem-solving methods and 2) a general computational architecture for resource-bounded problem solving. The decisiontheoretic models provide a set of principles for choosing among alternative problem-solving methods based on their relative costs and benefits, where benefits are characterized in terms of the value of information provided by the output of a reasoning activity. The output may be an estimate of some uncertain quantity or a recommendation for action. The computational architecture, called Schemer-ll, provides for interleaving of and communication among various problem-solving subsystems. These subsystems provide alternative approaches to information gathering, belief refinement, solution construction, and solution execution. In particular, the architecture provides a mechanism for interrupting the subsystems in response to critical events. We provide a decision theoretic account for scheduling problem-solving elements and for critical-event-driven interruption of activities in an architecture such as Schemer-II.


title: Induction and Uncertainty Management Techniques Applied to Veterinary Medical Diagnosis

abstract: This paper discusses a project undertaken between the Departments of Computing Science, Statistics, and the College of Veterinary Medicine to design a medical diagnostic system. On-line medical data has been collected in the hospital database system for several years. A number of induction methods are being used to extract knowledge from the data in an attempt to improve upon simple diagnostic charts used by the clinicians. They also enhance the results of classical statistical methods - finding many more significant variables. The second part of the paper describes an essentially Bayesian method of evidence combination using fuzzy events at an initial step. Results are presented and comparisons are made with other methods.


title: KNET: Integrating Hypermedia and Bayesian Modeling

abstract: KNET is a general-purpose shell for constructing expert systems based on belief networks and decision networks. Such networks serve as graphical representations for decision models, in which the knowledge engineer must define clearly the alternatives, states, preferences, and relationships that constitute a decision basis. KNET contains a knowledge-engineering core written in Object Pascal and an interface that tightly integrates HyperCard, a hypertext authoring tool for the Apple Macintosh computer, into a novel expert-system architecture. Hypertext and hypermedia have become increasingly important in the storage management, and retrieval of information. In broad terms, hypermedia deliver heterogeneous bits of information in dynamic, extensively cross-referenced packages. The resulting KNET system features a coherent probabilistic scheme for managing uncertainty, an objectoriented graphics editor for drawing and manipulating decision networks, and HyperCard's potential for quickly constructing flexible and friendly user interfaces. We envision KNET as a useful prototyping tool for our ongoing research on a variety of Bayesian reasoning problems, including tractable representation, inference, and explanation.


title: A Method for Using Belief Networks as Influence Diagrams

abstract: This paper demonstrates a method for using belief-network algorithms to solve influence diagram problems. In particular, both exact and approximation belief-network algorithms may be applied to solve influence-diagram problems. More generally, knowing the relationship between belief-network and influence-diagram problems may be useful in the design and development of more efficient influence diagram algorithms.


title: Process, Structure, and Modularity in Reasoning with Uncertainty

abstract: Computational mechanisms for uncertainty management must support interactive and incremental problem formulation, inference, hypothesis testing, and decision making. However, most current uncertainty inference systems concentrate primarily on inference, and provide no support for the larger issues. We present a computational approach to uncertainty management which provides direct support for the dynamic, incremental aspect of this task, while at the same time permitting direct representation of the structure of evidential relationships. At the same time, we show that this approach responds to the modularity concerns of Heckerman and Horvitz [Heck87]. This paper emphasizes examples of the capabilities of this approach. Another paper [D'Am89] details the representations and algorithms involved.


title: Probabilistic Causal Reasoning

abstract: Predicting the future is an important component of decision making. In most situations, however, there is not enough information to make accurate predictions. In this paper, we develop a theory of causal reasoning for predictive inference under uncertainty. We emphasize a common type of prediction that involves reasoning about persistence: whether or not a proposition once made true remains true at some later time. We provide a decision procedure with a polynomial-time algorithm for determining the probability of the possible consequences of a set events and initial conditions. The integration of simple probability theory with temporal projection enables us to circumvent problems that nonmonotonic temporal reasoning schemes have in dealing with persistence. The ideas in this paper have been implemented in a prototype system that refines a database of causal rules in the course of applying those rules to construct and carry out plans in a manufacturing domain.


title: Modeling uncertain and vague knowledge in possibility and evidence theories

abstract: This paper advocates the usefulness of new theories of uncertainty for the purpose of modeling some facets of uncertain knowledge, especially vagueness, in AI. It can be viewed as a partial reply to Cheeseman's (among others) defense of probability.


title: A Temporal Logic for Uncertain Events and An Outline of A Possible Implementation in An Extension of PROLOG

abstract: There is uncertainty associated with the occurrence of many events in real life. In this paper we develop a temporal logic to deal with such uncertain events and outline a possible implementation in an extension of PROLOG. Events are represented as fuzzy sets with the membership function giving the possibility of occurrence of the event in a given interval of time. The developed temporal logic is simple but powerful. It can determine effectively the various temporal relations between uncertain events or their combinations. PROLOG provides a uniform substrate on which to effectively implement such a temporal logic for uncertain events


title: Uncertainty Management for Fuzzy Decision Support Systems

abstract: A new approach for uncertainty management for fuzzy, rule based decision support systems is proposed: The domain expert's knowledge is expressed by a set of rules that frequently refer to vague and uncertain propositions. The certainty of propositions is represented using intervals [a, b] expressing that the proposition's probability is at least a and at most b. Methods and techniques for computing the overall certainty of fuzzy compound propositions that have been defined by using logical connectives 'and', 'or' and 'not' are introduced. Different inference schemas for applying fuzzy rules by using modus ponens are discussed. Different algorithms for combining evidence that has been received from different rules for the same proposition are provided. The relationship of the approach to other approaches is analyzed and its problems of knowledge acquisition and knowledge representation are discussed in some detail. The basic concepts of a rule-based programming language called PICASSO, for which the approach is a theoretical foundation, are outlined.


title: Probability as a Modal Operator

abstract: This paper argues for a modal view of probability. The syntax and semantics of one particularly strong probability logic are discussed and some examples of the use of the logic are provided. We show that it is both natural and useful to think of probability as a modal operator. Contrary to popular belief in AI, a probability ranging between 0 and 1 represents a continuum between impossibility and necessity, not between simple falsity and truth. The present work provides a clear semantics for quantification into the scope of the probability operator and for higher-order probabilities. Probability logic is a language for expressing both probabilistic and logical concepts.


title: Truth Maintenance Under Uncertainty

abstract: This paper addresses the problem of resolving errors under uncertainty in a rule-based system. A new approach has been developed that reformulates this problem as a neural-network learning problem. The strength and the fundamental limitations of this approach are explored and discussed. The main result is that neural heuristics can be applied to solve some but not all problems in rule-based systems.


title: Bayesian Assessment of a Connectionist Model for Fault Detection

abstract: A previous paper [2] showed how to generate a linear discriminant network (LDN) that computes likely faults for a noisy fault detection problem by using a modification of the perceptron learning algorithm called the pocket algorithm. Here we compare the performance of this connectionist model with performance of the optimal Bayesian decision rule for the example that was previously described. We find that for this particular problem the connectionist model performs about 97% as well as the optimal Bayesian procedure. We then define a more general class of noisy single-pattern boolean (NSB) fault detection problems where each fault corresponds to a single :pattern of boolean instrument readings and instruments are independently noisy. This is equivalent to specifying that instrument readings are probabilistic but conditionally independent given any particular fault. We prove: 1. The optimal Bayesian decision rule for every NSB fault detection problem is representable by an LDN containing no intermediate nodes. (This slightly extends a result first published by Minsky & Selfridge.) 2. Given an NSB fault detection problem, then with arbitrarily high probability after sufficient iterations the pocket algorithm will generate an LDN that computes an optimal Bayesian decision rule for that problem. In practice we find that a reasonable number of iterations of the pocket algorithm produces a network with good, but not optimal, performance.


title: On the Logic of Causal Models

abstract: This paper explores the role of Directed Acyclic Graphs (DAGs) as a representation of conditional independence relationships. We show that DAGs offer polynomially sound and complete inference mechanisms for inferring conditional independence relationships from a given causal set of such relationships. As a consequence, d-separation, a graphical criterion for identifying independencies in a DAG, is shown to uncover more valid independencies then any other criterion. In addition, we employ the Armstrong property of conditional independence to show that the dependence relationships displayed by a DAG are inherently consistent, i.e. for every DAG D there exists some probability distribution P that embodies all the conditional independencies displayed in D and none other.


title: The Optimality of Satisficing Solutions

abstract: This paper addresses a prevailing assumption in single-agent heuristic search theory- that problem-solving algorithms should guarantee shortest-path solutions, which are typically called optimal. Optimality implies a metric for judging solution quality, where the optimal solution is the solution with the highest quality. When path-length is the metric, we will distinguish such solutions as p-optimal.


title: An Empirical Comparison of Three Inference Methods

abstract: In this paper, an empirical evaluation of three inference methods for uncertain reasoning is presented in the context of Pathfinder, a large expert system for the diagnosis of lymph-node pathology. The inference procedures evaluated are (1) Bayes' theorem, assuming evidence is conditionally independent given each hypothesis; (2) odds-likelihood updating, assuming evidence is conditionally independent given each hypothesis and given the negation of each hypothesis; and (3) a inference method related to the Dempster-Shafer theory of belief. Both expert-rating and decision-theoretic metrics are used to compare the diagnostic accuracy of the inference methods.


title: Parallel Belief Revision

abstract: This paper describes a formal system of belief revision developed by Wolfgang Spohn and shows that this system has a parallel implementation that can be derived from an influence diagram in a manner similar to that in which Bayesian networks are derived. The proof rests upon completeness results for an axiomatization of the notion of conditional independence, with the Spohn system being used as a semantics for the relation of conditional independence.


title: Stochastic Sensitivity Analysis Using Fuzzy Influence Diagrams

abstract: The practice of stochastic sensitivity analysis described in the decision analysis literature is a testimonial to the need for considering deviations from precise point estimates of uncertainty. We propose the use of Bayesian fuzzy probabilities within an influence diagram computational scheme for performing sensitivity analysis during the solution of probabilistic inference and decision problems. Unlike other parametric approaches, the proposed scheme does not require resolving the problem for the varying probability point estimates. We claim that the solution to fuzzy influence diagrams provides as much information as the classical point estimate approach plus additional information concerning stochastic sensitivity. An example based on diagnostic decision making in microcomputer assembly is used to illustrate this idea. We claim that the solution to fuzzy influence diagrams provides as much information as the classical point estimate approach plus additional interval information that is useful for stochastic sensitivity analysis.


title: A Representation of Uncertainty to Aid Insight into Decision Models

abstract: Many real world models can be characterized as weak, meaning that there is significant uncertainty in both the data input and inferences. This lack of determinism makes it especially difficult for users of computer decision aids to understand and have confidence in the models. This paper presents a representation for uncertainty and utilities that serves as a framework for graphical summary and computer-generated explanation of decision models. The application described that tests the methodology is a computer decision aid designed to enhance the clinician-patient consultation process for patients with angina (chest pain due to lack of blood flow to the heart muscle). The angina model is represented as a Bayesian decision network. Additionally, the probabilities and utilities are treated as random variables with probability distributions on their range of possible values. The initial distributions represent information on all patients with anginal symptoms, and the approach allows for rapid tailoring to more patientspecific distributions. This framework provides a metric for judging the importance of each variable in the model dynamically.


title: Rational Nonmonotonic Reasoning

abstract: Nonmonotonic reasoning is a pattern of reasoning that allows an agent to make and retract (tentative) conclusions from inconclusive evidence. This paper gives a possible-worlds interpretation of the nonmonotonic reasoning problem based on standard decision theory and the emerging probability logic. The system's central principle is that a tentative conclusion is a decision to make a bet, not an assertion of fact. The system is rational, and as sound as the proof theory of its underlying probability log.


title: A Comparison of Decision Analysis and Expert Rules for Sequential Diagnosis

abstract: There has long been debate about the relative merits of decision theoretic methods and heuristic rule-based approaches for reasoning under uncertainty. We report an experimental comparison of the performance of the two approaches to troubleshooting, specifically to test selection for fault diagnosis. We use as experimental testbed the problem of diagnosing motorcycle engines. The first approach employs heuristic test selection rules obtained from expert mechanics. We compare it with the optimal decision analytic algorithm for test selection which employs estimated component failure probabilities and test costs. The decision analytic algorithm was found to reduce the expected cost (i.e. time) to arrive at a diagnosis by an average of 14% relative to the expert rules. Sensitivity analysis shows the results are quite robust to inaccuracy in the probability and cost estimates. This difference suggests some interesting implications for knowledge acquisition.


title: Multiple decision trees

abstract: This paper describes experiments, on two domains, to investigate the effect of averaging over predictions of multiple decision trees, instead of using a single tree. Other authors have pointed out theoretical and commonsense reasons for preferring the multiple tree approach. Ideally, we would like to consider predictions from all trees, weighted by their probability. However, there is a vast number of different trees, and it is difficult to estimate the probability of each tree. We sidestep the estimation problem by using a modified version of the ID3 algorithm to build good trees, and average over only these trees. Our results are encouraging. For each domain, we managed to produce a small number of good trees. We find that it is best to average across sets of trees with different structure; this usually gives better performance than any of the constituent trees, including the ID3 tree.


title: Probabilistic Inference and Probabilistic Reasoning

abstract: Uncertainty enters into human reasoning and inference in at least two ways. It is reasonable to suppose that there will be roles for these distinct uses of uncertainty also in automated reasoning.


title: Probabilistic and Non-Monotonic Inference

abstract: (l) I have enough evidence to render the sentence S probable. (la) So, relative to what I know, it is rational of me to believe S. (2) Now that I have more evidence, S may no longer be probable. (2a) So now, relative to what I know, it is not rational of me to believe S. These seem a perfectly ordinary, common sense, pair of situations. Generally and vaguely, I take them to embody what I shall call probabilistic inference. This form of inference is clearly non-monotonic. Relatively few people have taken this form of inference, based on high probability, to serve as a foundation for non-monotonic logic or for a logical or defeasible inference. There are exceptions: Jane Nutter [16] thinks that sometimes probability has something to do with non-monotonic reasoning. Judea Pearl [ 17] has recently been exploring the possibility. There are any number of people whom one might call probability enthusiasts who feel that probability provides all the answers by itself, with no need of help from logic. Cheeseman [1], Henrion [5] and others think it useful to look at a distribution of probabilities over a whole algebra of statements, to update that distribution in the light of new evidence, and to use the latest updated distribution of probability over the algebra as a basis for planning and decision making. A slightly weaker form of this approach is captured by Nilsson [15], where one assumes certain probabilities for certain statements, and infers the probabilities, or constraints on the probabilities of other statement. None of this corresponds to what I call probabilistic inference. All of the inference that is taking place, either in Bayesian updating, or in probabilistic logic, is strictly deductive. Deductive inference, particularly that concerned with the distribution of classical probabilities or chances, is of great importance. But this is not to say that there is no important role for what earlier logicians have called "ampliative" or "inductive" or "scientific" inference, in which the conclusion goes beyond the premises, asserts more than do the premises. This depends on what David Israel [6] has called "real rules of inference". It is characteristic of any such logic or inference procedure that it can go wrong: that statements accepted at one point may be rejected at a later point. Research underlying the results reported here has been partially supported by the Signals Warfare Center of the United States Army.


title: Epistemological Relevance and Statistical Knowledge

abstract: For many years, at least since McCarthy and Hayes (1969), writers have lamented, and attempted to compensate for, the alleged fact that we often do not have adequate statistical knowledge for governing the uncertainty of belief, for making uncertain inferences, and the like. It is hardly ever spelled out what "adequate statistical knowledge" would be, if we had it, and how adequate statistical knowledge could be used to control and regulate epistemic uncertainty.