title: Robust Agent Teams via Socially-Attentive Monitoring

abstract: Agents in dynamic multi-agent environments must monitor their peers to execute individual and group plans. A key open question is how much monitoring of other agents' states is required to be effective: The Monitoring Selectivity Problem. We investigate this question in the context of detecting failures in teams of cooperating agents, via Socially-Attentive Monitoring, which focuses on monitoring for failures in the social relationships between the agents. We empirically and analytically explore a family of socially-attentive teamwork monitoring algorithms in two dynamic, complex, multi-agent domains, under varying conditions of task distribution and uncertainty. We show that a centralized scheme using a complex algorithm trades correctness for completeness and requires monitoring all teammates. In contrast, a simple distributed teamwork monitoring algorithm results in correct and complete detection of teamwork failures, despite relying on limited, uncertain knowledge, and monitoring only key agents in a team. In addition, we report on the design of a socially-attentive monitoring system and demonstrate its generality in monitoring several coordination relationships, diagnosing detected failures, and both on-line and off-line applications.


title: What you see is not what you get: how sampling affects macroscopic features of biological networks

abstract: We use mathematical methods from the theory of tailored random graphs to study systematically the effects of sampling on topological features of large biological signalling networks. Our aim in doing so is to increase our quantitative understanding of the relation between true biological networks and the imperfect and often biased samples of these networks that are reported in public data repositories and used by biomedical scientists. We derive exact explicit formulae for degree distributions and degree correlation kernels of sampled networks, in terms of the degree distributions and degree correlation kernels of the underlying true network, for a broad family of sampling protocols that include (un-)biased node and/or link undersampling as well as (un-)biased link oversampling. Our predictions are in excellent agreement with numerical simulations.


title: On Deducing Conditional Independence from d-Separation in Causal Graphs with Feedback (Research Note)

abstract: Pearl and Dechter (1996) claimed that the d-separation criterion for conditional independence in acyclic causal networks also applies to networks of discrete variables that have feedback cycles, provided that the variables of the system are uniquely determined by the random disturbances. I show by example that this is not true in general. Some condition stronger than uniqueness is needed, such as the existence of a causal dynamics guaranteed to lead to the unique solution.


title: What's in an Attribute? Consequences for the Least Common Subsumer

abstract: Functional relationships between objects, called `attributes', are of considerable importance in knowledge representation languages, including Description Logics (DLs). A study of the literature indicates that papers have made, often implicitly, different assumptions about the nature of attributes: whether they are always required to have a value, or whether they can be partial functions. The work presented here is the first explicit study of this difference for subclasses of the CLASSIC DL, involving the same-as concept constructor. It is shown that although determining subsumption between concept descriptions has the same complexity (though requiring different algorithms), the story is different in the case of determining the least common subsumer (lcs). For attributes interpreted as partial functions, the lcs exists and can be computed relatively easily; even in this case our results correct and extend three previous papers about the lcs of DLs. In the case where attributes must have a value, the lcs may not exist, and even if it exists it may be of exponential size. Interestingly, it is possible to decide in polynomial time if the lcs exists.


title: The Complexity of Reasoning with Cardinality Restrictions and Nominals in Expressive Description Logics

abstract: We study the complexity of the combination of the Description Logics ALCQ and ALCQI with a terminological formalism based on cardinality restrictions on concepts. These combinations can naturally be embedded into C^2, the two variable fragment of predicate logic with counting quantifiers, which yields decidability in NExpTime. We show that this approach leads to an optimal solution for ALCQI, as ALCQI with cardinality restrictions has the same complexity as C^2 (NExpTime-complete). In contrast, we show that for ALCQ, the problem can be solved in ExpTime. This result is obtained by a reduction of reasoning with cardinality restrictions to reasoning with the (in general weaker) terminological formalism of general axioms for ALCQ extended with nominals. Using the same reduction, we show that, for the extension of ALCQI with nominals, reasoning with general axioms is a NExpTime-complete problem. Finally, we sharpen this result and show that pure concept satisfiability for ALCQI with nominals is NExpTime-complete. Without nominals, this problem is known to be PSpace-complete.


title: Backbone Fragility and the Local Search Cost Peak

abstract: The local search algorithm WSat is one of the most successful algorithms for solving the satisfiability (SAT) problem. It is notably effective at solving hard Random 3-SAT instances near the so-called `satisfiability threshold', but still shows a peak in search cost near the threshold and large variations in cost over different instances. We make a number of significant contributions to the analysis of WSat on high-cost random instances, using the recently-introduced concept of the backbone of a SAT instance. The backbone is the set of literals which are entailed by an instance. We find that the number of solutions predicts the cost well for small-backbone instances but is much less relevant for the large-backbone instances which appear near the threshold and dominate in the overconstrained region. We show a very strong correlation between search cost and the Hamming distance to the nearest solution early in WSat's search. This pattern leads us to introduce a measure of the backbone fragility of an instance, which indicates how persistent the backbone is as clauses are removed. We propose that high-cost random instances for local search are those with very large backbones which are also backbone-fragile. We suggest that the decay in cost beyond the satisfiability threshold is due to increasing backbone robustness (the opposite of backbone fragility). Our hypothesis makes three correct predictions. First, that the backbone robustness of an instance is negatively correlated with the local search cost when other factors are controlled for. Second, that backbone-minimal instances (which are 3-SAT instances altered so as to be more backbone-fragile) are unusually hard for WSat. Third, that the clauses most often unsatisfied during search are those whose deletion has the most effect on the backbone. In understanding the pathologies of local search methods, we hope to contribute to the development of new and better techniques.


title: An Application of Reinforcement Learning to Dialogue Strategy Selection in a Spoken Dialogue System for Email

abstract: This paper describes a novel method by which a spoken dialogue system can learn to choose an optimal dialogue strategy from its experience interacting with human users. The method is based on a combination of reinforcement learning and performance modeling of spoken dialogue systems. The reinforcement learning component applies Q-learning (Watkins, 1989), while the performance modeling component applies the PARADISE evaluation framework (Walker et al., 1997) to learn the performance function (reward) used in reinforcement learning. We illustrate the method with a spoken dialogue system named ELVIS (EmaiL Voice Interactive System), that supports access to email over the phone. We conduct a set of experiments for training an optimal dialogue strategy on a corpus of 219 dialogues in which human users interact with ELVIS over the phone. We then test that strategy on a corpus of 18 dialogues. We show that ELVIS can learn to optimize its strategy selection for agent initiative, for reading messages, and for summarizing email folders.


title: Nonapproximability Results for Partially Observable Markov Decision Processes

abstract: We show that for several variations of partially observable Markov decision processes, polynomial-time algorithms for finding control policies are unlikely to or simply don't have guarantees of finding policies within a constant factor or a constant summand of optimal. Here "unlikely" means "unless some complexity classes collapse," where the collapses considered are P=NP, P=PSPACE, or P=EXP. Until or unless these collapses are shown to hold, any control-policy designer must choose between such performance guarantees and efficient computation.


title: On Reasonable and Forced Goal Orderings and their Use in an Agenda-Driven Planning Algorithm

abstract: The paper addresses the problem of computing goal orderings, which is one of the longstanding issues in AI planning. It makes two new contributions. First, it formally defines and discusses two different goal orderings, which are called the reasonable and the forced ordering. Both orderings are defined for simple STRIPS operators as well as for more complex ADL operators supporting negation and conditional effects. The complexity of these orderings is investigated and their practical relevance is discussed. Secondly, two different methods to compute reasonable goal orderings are developed. One of them is based on planning graphs, while the other investigates the set of actions directly. Finally, it is shown how the ordering relations, which have been derived for a given set of goals G, can be used to compute a so-called goal agenda that divides G into an ordered set of subgoals. Any planner can then, in principle, use the goal agenda to plan for increasing sets of subgoals. This can lead to an exponential complexity reduction, as the solution to a complex planning problem is found by solving easier subproblems. Since only a polynomial overhead is caused by the goal agenda computation, a potential exists to dramatically speed up planning algorithms as we demonstrate in the empirical evaluation, where we use this method in the IPP planner.


title: Asimovian Adaptive Agents

abstract: The goal of this research is to develop agents that are adaptive and predictable and timely. At first blush, these three requirements seem contradictory. For example, adaptation risks introducing undesirable side effects, thereby making agents' behavior less predictable. Furthermore, although formal verification can assist in ensuring behavioral predictability, it is known to be time-consuming. Our solution to the challenge of satisfying all three requirements is the following. Agents have finite-state automaton plans, which are adapted online via evolutionary learning (perturbation) operators. To ensure that critical behavioral constraints are always satisfied, agents' plans are first formally verified. They are then reverified after every adaptation. If reverification concludes that constraints are violated, the plans are repaired. The main objective of this paper is to improve the efficiency of reverification after learning, so that agents have a sufficiently rapid response time. We present two solutions: positive results that certain learning operators are a priori guaranteed to preserve useful classes of behavioral assurance constraints (which implies that no reverification is needed for these operators), and efficient incremental reverification algorithms for those learning operators that have negative a priori results.


title: A Model of Inductive Bias Learning

abstract: A major problem in machine learning is that of inductive bias: how to choose a learner's hypothesis space so that it is large enough to contain a solution to the problem being learnt, yet small enough to ensure reliable generalization from reasonably-sized training sets. Typically such bias is supplied by hand through the skill and insights of experts. In this paper a model for automatically learning bias is investigated. The central assumption of the model is that the learner is embedded within an environment of related learning tasks. Within such an environment the learner can sample from multiple tasks, and hence it can search for a hypothesis space that contains good solutions to many of the problems in the environment. Under certain restrictions on the set of all hypothesis spaces available to the learner, we show that a hypothesis space that performs well on a sufficiently large number of training tasks will also perform well when learning novel tasks in the same environment. Explicit bounds are also derived demonstrating that learning multiple tasks within an environment of related tasks can potentially give much better generalization than learning a single task.


title: Mean Field Methods for a Special Class of Belief Networks

abstract: The chief aim of this paper is to propose mean-field approximations for a broad class of Belief networks, of which sigmoid and noisy-or networks can be seen as special cases. The approximations are based on a powerful mean-field theory suggested by Plefka. We show that Saul, Jaakkola and Jordan' s approach is the first order approximation in Plefka's approach, via a variational derivation. The application of Plefka's theory to belief networks is not computationally tractable. To tackle this problem we propose new approximations based on Taylor series. Small scale experiments show that the proposed schemes are attractive.


title: On the Compilability and Expressive Power of Propositional Planning Formalisms

abstract: The recent approaches of extending the GRAPHPLAN algorithm to handle more expressive planning formalisms raise the question of what the formal meaning of "expressive power" is. We formalize the intuition that expressive power is a measure of how concisely planning domains and plans can be expressed in a particular formalism by introducing the notion of "compilation schemes" between planning formalisms. Using this notion, we analyze the expressiveness of a large family of propositional planning formalisms, ranging from basic STRIPS to a formalism with conditional effects, partial state specifications, and propositional formulae in the preconditions. One of the results is that conditional effects cannot be compiled away if plan size should grow only linearly but can be compiled away if we allow for polynomial growth of the resulting plans. This result confirms that the recently proposed extensions to the GRAPHPLAN algorithm concerning conditional effects are optimal with respect to the "compilability" framework. Another result is that general propositional formulae cannot be compiled into conditional effects if the plan size should be preserved linearly. This implies that allowing general propositional formulae in preconditions and effect conditions adds another level of difficulty in generating a plan.


title: Technical Paper Recommendation: A Study in Combining Multiple Information Sources

abstract: The growing need to manage and exploit the proliferation of online data sources is opening up new opportunities for bringing people closer to the resources they need. For instance, consider a recommendation service through which researchers can receive daily pointers to journal papers in their fields of interest. We survey some of the known approaches to the problem of technical paper recommendation and ask how they can be extended to deal with multiple information sources. More specifically, we focus on a variant of this problem - recommending conference paper submissions to reviewing committee members - which offers us a testbed to try different approaches. Using WHIRL - an information integration system - we are able to implement different recommendation algorithms derived from information retrieval principles. We also use a novel autonomous procedure for gathering reviewer interest information from the Web. We evaluate our approach and compare it to other methods using preference data provided by members of the AAAI-98 conference reviewing committee along with data about the actual submissions.


title: Partial-Order Planning with Concurrent Interacting Actions

abstract: In order to generate plans for agents with multiple actuators, agent teams, or distributed controllers, we must be able to represent and plan using concurrent actions with interacting effects. This has historically been considered a challenging task requiring a temporal planner with the ability to reason explicitly about time. We show that with simple modifications, the STRIPS action representation language can be used to represent interacting actions. Moreover, algorithms for partial-order planning require only small modifications in order to be applied in such multiagent domains. We demonstrate this fact by developing a sound and complete partial-order planner for planning with concurrent interacting actions, POMP, that extends existing partial-order planners in a straightforward way. These results open the way to the use of partial-order planners for the centralized control of cooperative multiagent systems.


title: Planning by Rewriting

abstract: Domain-independent planning is a hard combinatorial problem. Taking into account plan quality makes the task even more difficult. This article introduces Planning by Rewriting (PbR), a new paradigm for efficient high-quality domain-independent planning. PbR exploits declarative plan-rewriting rules and efficient local search techniques to transform an easy-to-generate, but possibly suboptimal, initial plan into a high-quality plan. In addition to addressing the issues of planning efficiency and plan quality, this framework offers a new anytime planning algorithm. We have implemented this planner and applied it to several existing domains. The experimental results show that the PbR approach provides significant savings in planning effort while generating high-quality plans.


title: Speeding Up the Convergence of Value Iteration in Partially Observable Markov Decision Processes

abstract: Partially observable Markov decision processes (POMDPs) have recently become popular among many AI researchers because they serve as a natural model for planning under uncertainty. Value iteration is a well-known algorithm for finding optimal policies for POMDPs. It typically takes a large number of iterations to converge. This paper proposes a method for accelerating the convergence of value iteration. The method has been evaluated on an array of benchmark problems and was found to be very effective: It enabled value iteration to converge after only a few iterations on all the test problems.


title: Conformant Planning via Symbolic Model Checking

abstract: We tackle the problem of planning in nondeterministic domains, by presenting a new approach to conformant planning. Conformant planning is the problem of finding a sequence of actions that is guaranteed to achieve the goal despite the nondeterminism of the domain. Our approach is based on the representation of the planning domain as a finite state automaton. We use Symbolic Model Checking techniques, in particular Binary Decision Diagrams, to compactly represent and efficiently search the automaton. In this paper we make the following contributions. First, we present a general planning algorithm for conformant planning, which applies to fully nondeterministic domains, with uncertainty in the initial condition and in action effects. The algorithm is based on a breadth-first, backward search, and returns conformant plans of minimal length, if a solution to the planning problem exists, otherwise it terminates concluding that the problem admits no conformant solution. Second, we provide a symbolic representation of the search space based on Binary Decision Diagrams (BDDs), which is the basis for search techniques derived from symbolic model checking. The symbolic representation makes it possible to analyze potentially large sets of states and transitions in a single computation step, thus providing for an efficient implementation. Third, we present CMBP (Conformant Model Based Planner), an efficient implementation of the data structures and algorithm described above, directly based on BDD manipulations, which allows for a compact representation of the search layers and an efficient implementation of the search steps. Finally, we present an experimental comparison of our approach with the state-of-the-art conformant planners CGP, QBFPLAN and GPT. Our analysis includes all the planning problems from the distribution packages of these systems, plus other problems defined to stress a number of specific factors. Our approach appears to be the most effective: CMBP is strictly more expressive than QBFPLAN and CGP and, in all the problems where a comparison is possible, CMBP outperforms its competitors, sometimes by orders of magnitude.


title: AIS-BN: An Adaptive Importance Sampling Algorithm for Evidential Reasoning in Large Bayesian Networks

abstract: Stochastic sampling algorithms, while an attractive alternative to exact algorithms in very large Bayesian network models, have been observed to perform poorly in evidential reasoning with extremely unlikely evidence. To address this problem, we propose an adaptive importance sampling algorithm, AIS-BN, that shows promising convergence rates even under extreme conditions and seems to outperform the existing sampling algorithms consistently. Three sources of this performance improvement are (1) two heuristics for initialization of the importance function that are based on the theoretical properties of importance sampling in finite-dimensional integrals and the structural advantages of Bayesian networks, (2) a smooth learning method for the importance function, and (3) a dynamic weighting function for combining samples from different stages of the algorithm. We tested the performance of the AIS-BN algorithm along with two state of the art general purpose sampling algorithms, likelihood weighting (Fung and Chang, 1989; Shachter and Peot, 1989) and self-importance sampling (Shachter and Peot, 1989). We used in our tests three large real Bayesian network models available to the scientific community: the CPCS network (Pradhan et al., 1994), the PathFinder network (Heckerman, Horvitz, and Nathwani, 1990), and the ANDES network (Conati, Gertner, VanLehn, and Druzdzel, 1997), with evidence as unlikely as 10^-41. While the AIS-BN algorithm always performed better than the other two algorithms, in the majority of the test cases it achieved orders of magnitude improvement in precision of the results. Improvement in speed given a desired precision is even more dramatic, although we are unable to report numerical results here, as the other algorithms almost never achieved the precision reached even by the first few iterations of the AIS-BN algorithm.


title: Conflict-Directed Backjumping Revisited

abstract: In recent years, many improvements to backtracking algorithms for solving constraint satisfaction problems have been proposed. The techniques for improving backtracking algorithms can be conveniently classified as look-ahead schemes and look-back schemes. Unfortunately, look-ahead and look-back schemes are not entirely orthogonal as it has been observed empirically that the enhancement of look-ahead techniques is sometimes counterproductive to the effects of look-back techniques. In this paper, we focus on the relationship between the two most important look-ahead techniques---using a variable ordering heuristic and maintaining a level of local consistency during the backtracking search---and the look-back technique of conflict-directed backjumping (CBJ). We show that there exists a "perfect" dynamic variable ordering such that CBJ becomes redundant. We also show theoretically that as the level of local consistency that is maintained in the backtracking search is increased, the less that backjumping will be an improvement. Our theoretical results partially explain why a backtracking algorithm doing more in the look-ahead phase cannot benefit more from the backjumping look-back scheme. Finally, we show empirically that adding CBJ to a backtracking algorithm that maintains generalized arc consistency (GAC), an algorithm that we refer to as GAC-CBJ, can still provide orders of magnitude speedups. Our empirical results contrast with Bessiere and Regin's conclusion (1996) that CBJ is useless to an algorithm that maintains arc consistency.


title: Deligne pairing and determinant bundle

abstract: Let $X \rightarrow S$ be a smooth projective surjective morphism, where $X$ and $S$ are integral schemes over complex numbers. Let L_0, L_1, .... L_{n-1}, L_{n} be line bundles over $X$. There is a natural isomorphism of the Deligne pairing $<L_{0},...,L_{n}>$ with the determinant line bundle ${\rm Det}(\otimes_{i=0}^{n} (L_i- {\mathcal O}_{X}))$.


title: Grounding the Lexical Semantics of Verbs in Visual Perception using Force Dynamics and Event Logic

abstract: This paper presents an implemented system for recognizing the occurrence of events described by simple spatial-motion verbs in short image sequences. The semantics of these verbs is specified with event-logic expressions that describe changes in the state of force-dynamic relations between the participants of the event. An efficient finite representation is introduced for the infinite sets of intervals that occur when describing liquid and semi-liquid events. Additionally, an efficient procedure using this representation is presented for inferring occurrences of compound events, described with event-logic expressions, from occurrences of primitive events. Using force dynamics and event logic to specify the lexical semantics of events allows the system to be more robust than prior systems based on motion profile.


title: Popular Ensemble Methods: An Empirical Study

abstract: An ensemble consists of a set of individually trained classifiers (such as neural networks or decision trees) whose predictions are combined when classifying novel instances. Previous research has shown that an ensemble is often more accurate than any of the single classifiers in the ensemble. Bagging (Breiman, 1996c) and Boosting (Freund and Shapire, 1996; Shapire, 1990) are two relatively new but popular methods for producing ensembles. In this paper we evaluate these methods on 23 data sets using both neural networks and decision trees as our classification algorithm. Our results clearly indicate a number of conclusions. First, while Bagging is almost always more accurate than a single classifier, it is sometimes much less accurate than Boosting. On the other hand, Boosting can create ensembles that are less accurate than a single classifier -- especially when using neural networks. Analysis indicates that the performance of the Boosting methods is dependent on the characteristics of the data set being examined. In fact, further results show that Boosting ensembles may overfit noisy data sets, thus decreasing its performance. Finally, consistent with previous studies, our work suggests that most of the gain in an ensemble's performance comes in the first few classifiers combined; however, relatively large gains can be seen up to 25 classifiers when Boosting decision trees.


title: Asymmetric Berry-Phase Interference Patterns in a Single-Molecule Magnet

abstract: A Mn4 single-molecule magnet displays asymmetric Berry-phase interference patterns in the transverse-field (HT) dependence of the magnetization tunneling probability when a longitudinal field (HL) is present, contrary to symmetric patterns observed for HL=0. Reversal of HL results in a reflection of the transverse-field asymmetry about HT=0, as expected on the basis of the time-reversal invariance of the spin-orbit Hamiltonian which is responsible for the tunneling oscillations. A fascinating motion of Berry-phase minima within the transverse-field magnitude-direction phase space results from a competition between noncollinear magnetoanisotropy tensors at the two distinct Mn sites.


title: Coset enumeration for certain infinitely presented groups

abstract: We describe an algorithm that computes the index of a finitely generated subgroup in a finitely $L$-presented group provided that this index is finite. This algorithm shows that the subgroup membership problem for finite index subgroups in a finitely $L$-presented group is decidable. As an application, we consider the low-index subgroups of some self-similar groups including the Grigorchuk group, the twisted twin of the Grigorchuk group, the Grigorchuk super-group and the Hanoi 3-group


title: Comments on the paper `Unifying Boxy Bulge and Planar Long Bar in the Milky Way' by Martinez-Valpuesta & Gerhard [arXiv:1105.0928]

abstract: We comment on the recent paper by Martinez-Valpuesta & Gerhard (arXiv:1105.0928, 2011), who suggest, as an alternative to the bulge + long bar hypothesis in the inner 4 kpc of our Galaxy, a single boxy-bulge structure with a twisted major axis. In principle, we find this proposal acceptable; indeed, from a purely morphological point of view, this is more a question of semantics than science, and possibly all of us are talking about the same thing. However, we think that the particular features of this new proposal of a "single twisted bulge/bar" scenario leaves certain observational facts unexplained, whereas the model of a misaligned bulge + long bar successfully explains them.


title: Minimal length in quantum space and integrations of the line element in Noncommutative Geometry

abstract: We question the emergence of a minimal length in quantum spacetime, comparing two notions that appeared at various points in the literature: on the one side, the quantum length as the spectrum of an operator L in the Doplicher Fredenhagen Roberts (DFR) quantum spacetime, as well as in the canonical noncommutative spacetime; on the other side, Connes' spectral distance in noncommutative geometry. Although on the Euclidean space the two notions merge into the one of geodesic distance, they yield distinct results in the noncommutative framework. In particular on the Moyal plane, the quantum length is bounded above from zero while the spectral distance can take any real positive value, including infinity. We show how to solve this discrepancy by doubling the spectral triple. This leads us to introduce a modified quantum length d'_L, which coincides exactly with the spectral distance d_D on the set of states of optimal localization. On the set of eigenstates of the quantum harmonic oscillator - together with their translations - d'_L and d_D coincide asymptotically, both in the high energy and large translation limits. At small energy, we interpret the discrepancy between d'_L and d_D as two distinct ways of integrating the line element on a quantum space. This leads us to propose an equation for a geodesic on the Moyal plane.


title: 3D Structure of Microwave Sources from Solar Rotation Stereoscopy vs Magnetic Extrapolations

abstract: We use rotation stereoscopy to estimate the height of a steady-state solar feature relative to the photosphere, based on its apparent motion in the image plane recorded over several days of observation. The stereoscopy algorithm is adapted to work with either one- or two-dimensional data (i.e. from images or from observations that record the projected position of the source along an arbitrary axis). The accuracy of the algorithm is tested on simulated data, and then the algorithm is used to estimate the coronal radio source heights associated with the active region NOAA 10956, based on multifrequency imaging data over 7 days from the Siberian Solar Radio Telescope near 5.7 GHz, the Nobeyama Radio Heliograph at 17 GHz, as well as one-dimensional scans at multiple frequencies spanning the 5.98--15.95 GHz frequency range from the RATAN-600 instrument. The gyroresonance emission mechanism, which is sensitive to the coronal magnetic field strength, is applied to convert the estimated radio source heights at various frequencies, h(f), to information about magnetic field vs. height B(h), and the results are compared to a magnetic field extrapolation derived from photospheric magnetic field observations obtained by Hinode and MDI. We found that the gyroresonant emission comes from the heights exceeding location of the third gyrolayer irrespectively on the magnetic extrapolation method; implications of this finding for the coronal magnetography and coronal plasma physics are discussed.


title: Indirect stabilization of weakly coupled systems with hybrid boundary conditions

abstract: We investigate stability properties of indirectly damped systems of evolution equations in Hilbert spaces, under new compatibility assumptions. We prove polynomial decay for the energy of solutions and optimize our results by interpolation techniques, obtaining a full range of power-like decay rates. In particular, we give explicit estimates with respect to the initial data. We discuss several applications to hyperbolic systems with {\em hybrid} boundary conditions, including the coupling of two wave equations subject to Dirichlet and Robin type boundary conditions, respectively.


title: Achievable Degrees of Freedom of the K-user Interference Channel with Partial Cooperation

abstract: In this paper, we consider the K-user interference channel with partial cooperation, where a strict subset of the K users cooperate. For the K-user interference channel with cooperating subsets of length M, the outer bound of the total degrees of freedom is KM/(M+1). In this paper, we propose a signal space-based interference alignment scheme that proves the achievability of these degrees of freedom for the case K=M+2. The proposed scheme consists of a design for the transmit precoding matrices and a processing algorithm which we call the Successive Interference Alignment (SIA) algorithm. The decoder of each message uses the SIA algorithm to process the signals received by the M cooperating receivers in order to get the maximum available degrees of freedom.


title: Coherent strong-field control of multiple states by a single chirped femtosecond laser pulse

abstract: We present a joint experimental and theoretical study on strong-field photo-ionization of sodium atoms using chirped femtosecond laser pulses. By tuning the chirp parameter, selectivity among the population in the highly excited states 5p, 6p, 7p and 5f, 6f is achieved. Different excitation pathways enabling control are identified by simultaneous ionization and measurement of photoelectron angular distributions employing the velocity map imaging technique. Free electron wave packets at an energy of around 1 eV are observed. These photoelectrons originate from two channels. The predominant 2+1+1 Resonance Enhanced Multi-Photon Ionization (REMPI) proceeds via the strongly driven two-photon transition $4s\leftarrow\leftarrow3s$, and subsequent ionization from the states 5p, 6p and 7p whereas the second pathway involves 3+1 REMPI via the states 5f and 6f. In addition, electron wave packets from two-photon ionization of the non-resonant transiently populated state 3p are observed close to the ionization threshold. A mainly qualitative five-state model for the predominant excitation channel is studied theoretically to provide insights into the physical mechanisms at play. Our analysis shows that by tuning the chirp parameter the dynamics is effectively controlled by dynamic Stark-shifts and level crossings. In particular, we show that under the experimental conditions the passage through an uncommon three-state "bow-tie" level crossing allows the preparation of coherent superposition states.


title: QCD Resummation for Single Spin Asymmetries

abstract: We study the transverse momentum dependent factorization for single spin asymmetries in Drell-Yan and semi-inclusive deep inelastic scattering processes at one-loop order. The next-to-leading order hard factors are calculated in the Ji-Ma-Yuan factorization scheme. We further derive the QCD resummation formalisms for these observables following the Collins-Soper-Sterman method. The results are expressed in terms of the collinear correlation functions from initial and/or final state hadrons coupled with the Sudakov form factor containing all order soft-gluon resummation effects. The scheme independent coefficients are calculated up to one-loop order.


title: About Thinning Invariant Partition Structures

abstract: Bernoulli-$p$ thinning has been well-studied for point processes. Here we consider three other cases: (1) sequences $(X_1,X_2,...)$; (2) gaps of such sequences $(X_{n+1}-X_1)_{n\in\mathbb{N}}$; (3) partition structures. For the first case we characterize the distributions which are simultaneously invariant under Bernoulli-$p$ thinning for all $p \in (0,1]$. Based on this, we make conjectures for the latter two cases, and provide a potential approach for proof. We explain the relation to spin glasses, which is complementary to important previous work of Aizenman and Ruzmaikina, Arguin, and Shkolnikov.


title: A Maass Lifting of $\Theta^3$ and Class Numbers of Real and imaginary Quadratic Fields

abstract: We give an explicit construct of a harmonic weak Maass form $F_{\Theta}$ that is a "lift" of $\Theta^3$, where $\Theta$ is the classical Jacobi theta function. Just as the Fourier coefficients of $\Theta^3$ are related to class numbers of imaginary quadratic fields, the Fourier coefficients of the "holomorphic part" of $F_{\Theta}$ are associated to class numbers of real quadratic fields.


title: Harmonic analysis of spherical sampling in diffusion MRI

abstract: In the last decade diffusion MRI has become a powerful tool to non-invasively study white-matter integrity in the brain. Recently many research groups have focused their attention on multi-shell spherical acquisitions with the aim of effectively mapping the diffusion signal with a lower number of q-space samples, hence enabling a crucial reduction of acquisition time. One of the quantities commonly studied in this context is the so-called orientation distribution function (ODF). In this setting, the spherical harmonic (SH) transform has gained a great deal of popularity thanks to its ability to perform convolution operations efficiently and accurately, such as the Funk-Radon transform notably required for ODF computation from q-space data. However, if the q-space signal is described with an unsuitable angular resolution at any b-value probed, aliasing (or interpolation) artifacts are unavoidably created. So far this aspect has been tackled empirically and, to our knowledge, no study has addressed this problem in a quantitative approach. The aim of the present work is to study more theoretically the efficiency of multi-shell spherical sampling in diffusion MRI, in order to gain understanding in HYDI-like approaches, possibly paving the way to further optimization strategies.


title: Scalable uniform construction of highly-conditional quantum gates

abstract: We present a scalable uniform technique for construction of highly conditional C$^n$-NOT quantum gates of trapped ion qubits, such as the Toffoli gate, without using ancilla states and circuits of an exorbitant number of concatenated one- and two-qubit gates. Apart from the initial dressing of the internal qubit states with vibrational phonons and the final restoration of the phonon ground state, our technique requires the application of just a single composite pulse on the target qubit and is applicable both in and outside the Lamb-Dicke regime. We design special narrowband composite pulses, which suppress all transitions but the conditional transition of the target qubit; moreover, these composite pulses significantly improve the spatial addressing selectivity.


title: Random stable laminations of the disk

abstract: We study large random dissections of polygons. We consider random dissections of a regular polygon with $n$ sides, which are chosen according to Boltzmann weights in the domain of attraction of a stable law of index $\theta\in(1,2]$. As $n$ goes to infinity, we prove that these random dissections converge in distribution toward a random compact set, called the random stable lamination. If $\theta=2$, we recover Aldous' Brownian triangulation. However, if $\theta\in(1,2)$, large faces remain in the limit and a different random compact set appears. We show that the random stable lamination can be coded by the continuous-time height function associated to the normalized excursion of a strictly stable spectrally positive L\'{e}vy process of index $\theta$. Using this coding, we establish that the Hausdorff dimension of the stable random lamination is almost surely $2-1/\theta$.


title: High-fidelity local addressing of trapped ions and atoms by composite sequences of laser pulses

abstract: A vital requirement for a quantum computer is the ability to locally address, with high fidelity, any of its qubits without affecting their neighbors. We propose an addressing method using composite sequences of laser pulses, which reduces dramatically the addressing error in a lattice of closely spaced atoms or ions, and at the same time significantly enhances the robustness of qubit manipulations. To this end, we design novel high-fidelity composite pulses for the most important single-qubit operations. In principle, this method allows one to beat the diffraction limit, for only atoms situated in a small spatial region around the center of the laser beam are excited, well within the laser beam waist.


title: Laminations in the language of leaves

abstract: Thurston defined invariant laminations, i.e. collections of chords of the unit circle $S^1$ (called \emph{leaves}) that are pairwise disjoint inside the open unit disk and satisfy a few dynamical properties. To be directly associated to a polynomial, a lamination has to be generated by an equivalence relation with specific properties on $S^1$; then it is called a \emph{q-lamination}. Since not all laminations are q-laminations, then from the point of view of studying polynomials the most interesting are those of them which are limits of q-laminations. In this paper we introduce an alternative definition of an invariant lamination, which involves only conditions on the leaves (and avoids gap invariance). The new class of laminations is slightly smaller than that defined by Thurston and is closed. We use this notion to elucidate the connection between invariant laminations and invariant equivalence relations on $S^1$.


title: Electron-cylotron maser radiation from electron holes: Upward current region

abstract: Electron holes are suggested to be an important source for generation of electron-cyclotron maser radiation. We demonstrate that electron holes generated in a ring-horseshoe distribution in the auroral-kilometric radiation source region have the capacity to emit band-limited radiation. The radiation is calculated in the proper frame of a circular model hole and shown to be strictly perpendicular in this frame. Its bandwidth under auroral conditions is of the order of $\sim1$ kHz, which is a reasonable value. {It is also shown that much of the drift of fine structure in the radiation can be interpreted as Doppler shift. Estimates based on data are in good agreement with theory. Growth and absorption rates have been obtained for the emitted radiation. However, the growth rate of a single hole obtained under conservative conditions is small, too small for reproducing the observed fine structure flux. Trapping of radiation inside the hole for the hole's lifetime helps amplifying the radiation additionally but introduces other problems. This entire set of questions is discussed at length and compared to radiation from the global horseshoe distribution.} The interior of the hole produces a weak absorption at slightly higher frequency than emission. The absorptivity is roughly two orders of magnitude below the growth rate of the radiation thus being weak even when the emission and absorption bands overlap. Transforming to the stationary observer's frame it is found that the radiation becomes oblique against the magnetic field. For approaching holes the radiated frequencies may even exceed the local electron cyclotron frequency.


title: Probing the high momentum component of the deuteron at high Q^2

abstract: The d(e,e'p) cross section at a momentum transfer of 3.5 (GeV/c)^2 was measured over a kinematical range that made it possible to study this reaction for a set of fixed missing momenta as a function of the neutron recoil angle theta_nq and to extract missing momentum distributions for fixed values of theta_nq up to 0.55 GeV/c. In the region of 35 (deg) <= theta_nq <= 45 (deg) recent calculations, which predict that final state interactions are small, agree reasonably well with the experimental data. Therefore these experimental reduced cross sections provide direct access to the high momentum component of the deuteron momentum distribution in exclusive deuteron electro-disintegration.


title: Full counting statistics of energy fluctuations in a driven quantum resonator

abstract: We consider the statistics of time-integrated energy fluctuations of a driven bosonic resonator (as measured by a QND detector), using the standard Keldysh prescription to define higher moments. We find that due to an effective cascading of fluctuations, these statistics are surprisingly non-classical: the low-temperature, quantum probability distribution is not equivalent to the high-temperature classical distribution evaluated at some effective temperature. Moreover, for a sufficiently large drive detuning and low temperatures, the Keldysh-ordered quasi-probability distribution characterizing these fluctuations fails to be positive-definite; this is similar to the full counting statistics of charge in superconducting systems. We argue that this indicates a kind of non-classical behaviour akin to that tested by Leggett-Garg inequalities.


title: Anomaly Induced Transport in Arbitrary Dimensions

abstract: Motivated by the consistency of a global anomaly with the second law of thermodynamics, we propose a form for the anomaly induced charge/energy transport in arbitrary even dimensions. In a given dimension, this form exhausts all second law constraints on anomaly induced transport at any given order in hydrodynamic derivative expansion. This is achieved by solving the second law constraints off-shell without resorting to hydrodynamic equations at lower orders. We also study various possible finite temperature corrections to such anomaly induced transport coefficients.


title: Rabi-vibronic resonance with large number of vibrational quanta

abstract: We study theoretically the Rabi oscillations of a resonantly driven two-level system linearly coupled to a harmonic oscillator (vibrational mode) with frequency, \omega_0. We show that for weak coupling, \omega_p \ll \omega_0, where \omega_p is the polaronic shift, Rabi oscillations are strongly modified in the vicinity of the Rabi-vibronic resonance \Omega_R = \omega_0, where \Omega_R is the Rabi frequency. The width of the resonance is (\Omega_R-\omega_0) \sim \omega_p^{2/3} \omega_0^{1/3} \gg \omega_p. Within this domain of \Omega_R the actual frequency of the Rabi oscillations exhibits a bistable behavior as a function of \Omega_R. Importantly, within the resonant domain, the oscillator is highly excited, which allows to treat it classically.


title: String theory of the Omega deformation

abstract: In this article, we construct a supersymmetric real mass deformation for the adjoint chiral multiplets in the gauge theory describing the dynamics of a stack of D2-branes in type II string theory. We do so by placing the D2-branes into the T-dual of a supersymmetric NS fluxbrane background. We furthermore note that this background is the string theoretic realization of an Omega-deformation of flat space in the directions transverse to the branes where the deformation parameters satisfy \epsilon_1 = - \epsilon_2. This \Omega-deformation therefore serves to give supersymmetric real masses to the chiral multiplets of the 3D gauge theory. To illustrate the physical effect of the real mass term, we derive BPS-saturated classical solutions for the branes rotating in this background. Finally, we reproduce all the same structure in the presence of NS fivebranes and comment on the relationship to the gauge theory/spin-chain correspondence of Nekrasov and Shatashvili.


title: Relativistic Quantum Information: developments in Quantum Information in general relativistic scenarios

abstract: Recently, there has been increased interest in understanding entanglement and quantum communication in black hole spacetimes and in using quantum information techniques to address questions in gravity. Studies on relativistic entanglement show the emergence of conceptually important qualitative differences to a non-relativistic treatment. For instance, entanglement was found to be an observer-dependent property that changes from the perspective of accelerated observers moving in flat spacetime. Relativisitic quantum information theory uses well-known tools coming from quantum information and quantum optics to study quantum effects provoked by gravity to learn information about the spacetime. We can take advantage of our knowledge about quantum correlations and effects produced by the gravitational interaction to set the basis for experimental proposals ultimately aiming at finding corrections due to quantum gravity effects, too mild to be directly observed. This doctoral thesis dissertation summarises most of the research carried out during my PhD on this topic.


title: The media effect in Axelrod's model explained

abstract: We revisit the problem of introducing an external global field -- the mass media -- in Axelrod's model of social dynamics, where in addition to their nearest neighbors, the agents can interact with a virtual neighbor whose cultural features are fixed from the outset. The finding that this apparently homogenizing field actually increases the cultural diversity has been considered a puzzle since the phenomenon was first reported more than a decade ago. Here we offer a simple explanation for it, which is based on the pedestrian observation that Axelrod's model exhibits more cultural diversity, i.e., more distinct cultural domains, when the agents are allowed to interact solely with the media field than when they can interact with their neighbors as well. In this perspective, it is the local homogenizing interactions that work towards making the absorbing configurations less fragmented as compared with the extreme situation in which the agents interact with the media only.


title: Random Latin square graphs

abstract: In this paper we introduce new models of random graphs, arising from Latin squares which include random Cayley graphs as a special case. We investigate some properties of these graphs including their clique, independence and chromatic numbers, their expansion properties as well as their connectivity and Hamiltonicity. The results obtained are compared with other models of random graphs and several similarities and differences are pointed out. For many properties our results for the general case are as strong as the known results for random Cayley graphs and sometimes improve the previously best results for the Cayley case.


title: Global simulations of the magnetic field evolution under the influence of the cosmic-ray-driven dynamo

abstract: We present three-dimensional global numerical simulations of the cosmic-ray (CR) driven dynamo in barred galaxies. We study the evolution of the interstellar medium of the barred galaxy in the presence of non-axisymmetric component of the potential, i.e., the bar. The magnetohydrodynamical dynamo is driven by CRs, which are continuously supplied to the disk by supernova (SN) remnants. No magnetic field is present at the beginning of simulations but one-tenth of SN explosions is a source of a small-scale randomly oriented dipolar magnetic field. In all models we assume that 10% of 10^51 erg SN kinetic energy output is converted into CR energy. To compare our results directly with the observed properties of galaxies, we construct realistic maps of polarized radio emission. The main result is that the CR-driven dynamo can amplify weak magnetic fields up to a few {\mu}G within a few Gyr in barred galaxies. The obtained e-folding time is equal to 300 Myr and the magnetic field reaches equipartition at time t ~ 4.0 Gyr. Initially, the completely random magnetic field evolves into large-scale structures. An even (quadrupole-type) configuration of the magnetic field with respect to the galactic plane can be observed. Additionally, the modeled magnetic field configuration resembles maps of the polarized intensity observed in barred galaxies. Polarization vectors are distributed along the bar and between spiral arms. Moreover, the drift of magnetic arms with respect to the spiral pattern in the gas density distribution is observed during the entire simulation time.


title: An Evolutionary Algorithm with Advanced Goal and Priority Specification for Multi-objective Optimization

abstract: This paper presents an evolutionary algorithm with a new goal-sequence domination scheme for better decision support in multi-objective optimization. The approach allows the inclusion of advanced hard/soft priority and constraint information on each objective component, and is capable of incorporating multiple specifications with overlapping or non-overlapping objective functions via logical 'OR' and 'AND' connectives to drive the search towards multiple regions of trade-off. In addition, we propose a dynamic sharing scheme that is simple and adaptively estimated according to the on-line population distribution without needing any a priori parameter setting. Each feature in the proposed algorithm is examined to show its respective contribution, and the performance of the algorithm is compared with other evolutionary optimization methods. It is shown that the proposed algorithm has performed well in the diversity of evolutionary search and uniform distribution of non-dominated individuals along the final trade-offs, without significant computational effort. The algorithm is also applied to the design optimization of a practical servo control system for hard disk drives with a single voice-coil-motor actuator. Results of the evolutionary designed servo control system show a superior closed-loop performance compared to classical PID or RPT approaches.


title: The GRT Planning System: Backward Heuristic Construction in Forward State-Space Planning

abstract: This paper presents GRT, a domain-independent heuristic planning system for STRIPS worlds. GRT solves problems in two phases. In the pre-processing phase, it estimates the distance between each fact and the goals of the problem, in a backward direction. Then, in the search phase, these estimates are used in order to further estimate the distance between each intermediate state and the goals, guiding so the search process in a forward direction and on a best-first basis. The paper presents the benefits from the adoption of opposite directions between the preprocessing and the search phases, discusses some difficulties that arise in the pre-processing phase and introduces techniques to cope with them. Moreover, it presents several methods of improving the efficiency of the heuristic, by enriching the representation and by reducing the size of the problem. Finally, a method of overcoming local optimal states, based on domain axioms, is proposed. According to it, difficult problems are decomposed into easier sub-problems that have to be solved sequentially. The performance results from various domains, including those of the recent planning competitions, show that GRT is among the fastest planners.


title: Popularity versus Similarity in Growing Networks

abstract: Popularity is attractive -- this is the formula underlying preferential attachment, a popular explanation for the emergence of scaling in growing networks. If new connections are made preferentially to more popular nodes, then the resulting distribution of the number of connections that nodes have follows power laws observed in many real networks. Preferential attachment has been directly validated for some real networks, including the Internet. Preferential attachment can also be a consequence of different underlying processes based on node fitness, ranking, optimization, random walks, or duplication. Here we show that popularity is just one dimension of attractiveness. Another dimension is similarity. We develop a framework where new connections, instead of preferring popular nodes, optimize certain trade-offs between popularity and similarity. The framework admits a geometric interpretation, in which popularity preference emerges from local optimization. As opposed to preferential attachment, the optimization framework accurately describes large-scale evolution of technological (Internet), social (web of trust), and biological (E.coli metabolic) networks, predicting the probability of new links in them with a remarkable precision. The developed framework can thus be used for predicting new links in evolving networks, and provides a different perspective on preferential attachment as an emergent phenomenon.

categories: physics.soc-ph cond-mat.stat-mech cs.NI cs.SI


title: Decomposition of operator semigroups on W*-algebras

abstract: We consider semigroups of operators on a W$^*$-algebra and prove, under appropriate assumptions, the existence of a Jacobs-DeLeeuw-Glicksberg type decomposition. This decomposition splits the algebra into a "stable" and "reversible" part with respect to the semigroup and yields, among others, a structural approach to the Perron-Frobenius spectral theory for completely positive operators on W$^*$-algebras.


title: Emergence of Bursts and Communities in Evolving Weighted Networks

abstract: Understanding the patterns of human dynamics and social interaction, and the way they lead to the formation of an organized and functional society are important issues especially for techno-social development. Addressing these issues of social networks has recently become possible through large scale data analysis of e.g. mobile phone call records, which has revealed the existence of modular or community structure with many links between nodes of the same community and relatively few links between nodes of different communities. The weights of links, e.g. the number of calls between two users, and the network topology are found correlated such that intra-community links are stronger compared to the weak inter-community links. This is known as Granovetter's "The strength of weak ties" hypothesis. In addition to this inhomogeneous community structure, the temporal patterns of human dynamics turn out to be inhomogeneous or bursty, characterized by the heavy tailed distribution of inter-event time between two consecutive events. In this paper, we study how the community structure and the bursty dynamics emerge together in an evolving weighted network model. The principal mechanisms behind these patterns are social interaction by cyclic closure, i.e. links to friends of friends and the focal closure, i.e. links to individuals sharing similar attributes or interests, and human dynamics by task handling process. These three mechanisms have been implemented as a network model with local attachment, global attachment, and priority-based queuing processes. By comprehensive numerical simulations we show that the interplay of these mechanisms leads to the emergence of heavy tailed inter-event time distribution and the evolution of Granovetter-type community structure. Moreover, the numerical results are found to be in qualitative agreement with empirical results from mobile phone call dataset.


title: Locally Inaccessible Information as a Fundamental Ingredient to Quantum Information

abstract: Quantum discord (QD) measures the fraction of the pairwise mutual information that is locally inaccessible, in a multipartite system. Fundamental aspects related to two important measures in quantum information theory the Entanglement of Formation (EOF) and the conditional entropy, can be understood in terms of the distribution of this form of Local Inaccessible Information (LII). As such, the EOF for an arbitrarily mixed bipartite system AB can be related to the gain or loss of LII due to the extra knowledge that a purifying ancillary system E has on the pair AB. Similarly, a clear meaning of the negativity of the conditional entropy for AB is given. We exemplify by showing that these relations elucidate important and yet not well understood quantum features, such as the bipartite entanglement sudden death and the distinction between EOF and QD for quantifying quantum correlation. For that we introduce the concept of LII flow which quantifies the LII shared in multipartite system when a sequential local measurements are performed.


title: On a problem of Erd\H{o}s and Rothschild on edges in triangles

abstract: Erd\H{o}s and Rothschild asked to estimate the maximum number, denoted by H(N,C), such that every N-vertex graph with at least CN^2 edges, each of which is contained in at least one triangle, must contain an edge that is in at least H(N,C) triangles. In particular, Erd\H{o}s asked in 1987 to determine whether for every C>0 there is \epsilon >0 such that H(N,C) > N^\epsilon, for all sufficiently large N. We prove that H(N,C) = N^{O(1/log log N)} for every fixed C < 1/4. This gives a negative answer to the question of Erd\H{o}s, and is best possible in terms of the range for C, as it is known that every N-vertex graph with more than (N^2)/4 edges contains an edge that is in at least N/6 triangles.


title: IVOA Architecture

abstract: This note describes the technical architecture of the IVOA. The description is decomposed into three levels. Level 0 is a general, high level summary of the IVOA Architecture. Level 1 provides more details about components and functionalities, still without being overly technical. Finally, Level 2 displays how the IVOA standards fit into the IVOA Architecture.


title: Conserved cosmological perturbation in Galileon models

abstract: We prove the existence of a fully nonlinear conserved curvature perturbation on large scales in Galileon-type scalar field models in two approaches. The first approach is based on the conservation of energy-momentum tensor of the Galileon field, which is also the familiar approach in understanding the conservation in $k$-essence or perfect fluid models. We show that the fluid corresponding to the Galileon field becomes perfect and barotropic on large scales, which is responsible to the conservation. The difference from $k$-essence model is that, besides the energy-momentum conservation, the Einstein equation must be employed in order to complete the proof of barotropy. In the second approach, we derive the fully non-perturbative action for the curvature perturbation $\zeta$ in Galileon model on large scales, and argue that $\zeta=const$ is indeed an exact solution on large scales. This conservation of curvature perturbation is important since it relates the later and the primordial universe.


title: Monopoles in the Plaquette Formulation of the 3D SU(2) Lattice Gauge Theory

abstract: Using a plaquette formulation for lattice gauge models we describe monopoles of the three dimensional SU(2) theory which appear as configurations in the complete axial gauge and violate the continuum Bianchi identity. Furthemore we derive a dual formulation for the Wilson loop in arbitrary representation and calculate the form of the interaction between generated electric flux and monopoles in the region of a weak coupling relevant for the continuum limit. The effective theory which controls the interaction is of the sine-Gordon type model. The string tension is calculated within the semiclassical approximation.


title: Dissecting the Wjj Anomaly: Diagnostic Tests of a Leptophobic Z'

abstract: We examine the scenario where a leptophobic Z' boson accounts for the excess of events in the Wjj channel as observed by CDF. We assume generation independent couplings for the Z' and obtain allowed regions for the four hadronic couplings using the cross section range quoted by CDF as well as constraints from dijet production at UA2. These coupling regions translate into well-determined rates for the associated production of Z/\gamma+Z' at the Tevatron and LHC, as well as W+Z' at the LHC,that are directly correlated with the Wjj rate observed at the Tevatron. The Wjj rate at the LHC is large and this channel should be observed soon once the SM backgrounds are under control. The rates for Z/\gamma+Z' associated production are smaller, and these processes should not yet have been observed at the Tevatron given the expected SM backgrounds. In addition, we also show that valuable coupling information is obtainable from the distributions of other kinematic variables, e.g., M_{WZ'}, p_T^W, and \cos \theta_W^*. Once detected, these associated production processes and the corresponding kinematic distributions examined here will provide further valuable information on the Z' boson couplings.


title: Discrete to continuum transition in multifractal spacetimes

abstract: We outline a field theory on a multifractal spacetime. The measure in the action is characterized by a varying Hausdorff dimension and logarithmic oscillations governed by a fundamental physical length. A fine hierarchy of length scales identifies different regimes, from a microscopic structure with discrete symmetries to an effectively continuum spacetime. Thanks to general arguments from fractal geometry, this scenario explicitly realizes two indirect or conjectured features of most quantum gravity models: a change of effective spacetime dimensionality with the probed scale, and the transition from a fundamentally discrete quantum spacetime to the continuum. It also allows us to probe ultramicroscopic scales where spectral methods based on ordinary geometry typically fail. Consequences for noncommutative field theories are discussed.


title: The Emergence of Leadership in Social Networks

abstract: We study a networked version of the minority game in which agents can choose to follow the choices made by a neighbouring agent in a social network. We show that for a wide variety of networks a leadership structure always emerges, with most agents following the choice made by a few agents. We find a suitable parameterisation which highlights the universal aspects of the behaviour and which also indicates where results depend on the type of social network.


title: Hypergraphs and City Street Networks

abstract: The map of a city's streets constitutes a particular case of spatial complex network. However a city is not limited to its topology: it is above all a geometrical object whose particularity is to organize into short and long axes called streets. In this article we present and discuss two algorithms aiming at recovering the notion of street from a graph representation of a city. Then we show that the length of the so-called streets scales logarithmically. This phenomenon leads to assume that a city is shaped into a logic of extension and division of space.


title: The Physics of the 'Heartbeat' State of GRS 1915+105

abstract: We present the first detailed phase-resolved spectral analysis of a joint Chandra High Energy Transmission Grating Spectrometer and Rossi X-ray Timing Explorer observation of the rho variability class in the microquasar GRS 1915+105. The rho cycle displays a high-amplitude, double-peaked flare that recurs roughly every 50 s, and is sometimes referred to as the "heartbeat" oscillation. The spectral and timing properties of the oscillation are consistent with the radiation pressure instability and the evolution of a local Eddington limit in the inner disk. We exploit strong variations in the X-ray continuum, iron emission lines, and the accretion disk wind to probe the accretion geometry over nearly six orders of magnitude in distance from the black hole. At small scales (1-10 R_g), we detect a burst of bremsstrahlung emission that appears to occur when a portion of the inner accretion disk evaporates due to radiation pressure. Jet activity, as inferred from the appearance of a short X-ray hard state, seems to be limited to times near minimum luminosity, with a duty cycle of ~10%. On larger scales (1e5-1e6 R_g) we use detailed photoionization arguments to track the relationship between the fast X-ray variability and the accretion disk wind. For the first time, we are able to show that changes in the broadband X-ray spectrum produce changes in the structure and density of the accretion disk wind on timescales as short as 5 seconds. These results clearly establish a causal link between the X-ray oscillations and the disk wind and therefore support the existence of a disk-jet-wind connection. Furthermore, our analysis shows that the mass loss rate in the wind may be sufficient to cause long-term oscillations in the accretion rate, leading to state transitions in GRS 1915+105.


title: Dark Energy and Neutrino Masses from Future Measurements of the Expansion History and Growth of Structure

abstract: We forecast the expected cosmological constraints from a combination of probes of both the universal expansion rate and matter perturbation growth, in the form of weak lensing tomography, galaxy tomography, supernovae, and the cosmic microwave background incorporating all cross-correlations between the observables for an extensive cosmological parameter set. We allow for non-zero curvature and parameterize our ignorance of the early universe by allowing for a non-negligible fraction of dark energy (DE) at high redshifts. We find that early DE density can be constrained to 0.2% of the critical density of the universe with Planck combined with a ground-based LSST-like survey, while curvature can be constrained to 0.06%. However, these additional degrees of freedom degrade our ability to measure late-time dark energy and the sum of neutrino masses. We find that the combination of cosmological probes can break degeneracies and constrain the sum of neutrino masses to 0.04 eV, present DE density also to 0.2% of the critical density, and the equation of state to 0.01 - roughly a factor of two degradation in the constraints overall compared to the case without allowing for early DE. The constraints for a space-based mission are similar. Even a modest 1% dark energy fraction of the critical density at high redshift, if not accounted for in future analyses, biases the cosmological parameters by up to 2 sigma. Our analysis suggests that throwing out nonlinear scales (multipoles > 1000) may not result in significant degradation in future parameter measurements when multiple cosmological probes are combined. We find that including cross-correlations between the different probes can result in improved constraints by up to a factor of 2 for the sum of neutrino masses and early dark energy density.


title: Extended narrow-line emission in the bright Seyfert 1.5 galaxy HE 2211-3903

abstract: Extended narrow-line regions (ENLRs) and extended emission-line regions (EELRs) have been the focus of integral field spectroscopy aiming at the inner kiloparsecs of nearby Seyfert galaxies as well as the larger environment of high redshift QSOs. Based on observations with the Wide Field Spectrograph WiFeS at the 2.3 m telescope of the Australian National University, we present spatially resolved emission-line diagnostics of the bright Seyfert 1.5 galaxy HE 2211-3903 which is drawn from a sample of the brightest Seyfert galaxies at z<0.06 with luminosities around the classical Seyfert/QSO demarcation. In addition to the previously known spiral arms of HE 2211-3903, the emission-line maps reveal a large scale ring with a radius of about 6 kpc which is connected to the active galactic nucleus (AGN) through a bar-like structure. The overall gas kinematics indicates a disk rotation pattern. The emission-line ratios show Seyfert-type, HII region-type, and composite classifications, while there is no strong evidence of LINER-type ratios. Shock ionization is likely to be negligible throughout the galaxy. The composite line ratios are explained via a mixing line between AGN and HII region photoionization. Composite line ratios are predominantly found in between the HII regions in the circum-nuclear region, the bar-like structure to the east of the nucleus, and the eastern half of the ring, suggesting AGN photoionization of the low-density interstellar medium in an ENLR on galaxy scales. The line ratios in the nucleus indicate N-enrichment, which is discussed in terms of chemical enrichment by Wolf-Rayet and Asymptotic Giant Branch stars during past and ongoing nuclear starburst activity.


title: Metallicity-Dependent quenching of Star Formation at High Redshift in Small Galaxies

abstract: [abridged] The star formation rates (SFR) of low-metallicity galaxies depend sensitively on the gas metallicity, because metals are crucial to mediating the transition from intermediate-temperature atomic gas to cold molecular gas, a necessary precursor to star formation. We study the impact of this effect on the star formation history of galaxies. We incorporate metallicity-dependent star formation and metal enrichment in a simple model that follows the evolution of a halo main progenitor. Our model shows that including the effect of metallicity leads to suppression of star formation at redshift z>2 in dark halos with masses <~ 10^11 Msun, with the suppression becoming near total for halos below ~10^9.5-10 Msun. We find that at high redshift the SFR cannot catch up with the gas inflow rate (IR), because the SFR is limited by the free-fall time, and because it is suppressed further by a lack of metals. As a result, in each galaxy the SFR is growing in time faster than the IR, and the integrated cosmic SFR density is rising with time. The suppressed in situ SFR at high z makes the growth of stellar mass dominated by ex situ SFR which implies that the specific SFR (sSFR) remains constant with time. The intensely accreted gas at high z is accumulating as an atomic gas reservoir. This provides additional fuel for star formation in 10^10 - 10^12 Msun halos at z ~ 1-3, which allows the SFR to exceed the instantaneous IR, and may enable an even higher outflow rate. At z<1, following the natural decline in IR with time due to the universal expansion, the SFR and sSFR are expected to drop. We specify the expected dependence of sSFR and metallicity on stellar mass and redshift. At a given z, and below a critical mass, these relations are predicted to be flat and rising respectively. Our model predictions qualitatively match some of the puzzling features in the observed star formation history.


title: Floquet Spectrum and Transport Through an Irradiated Graphene Ribbon

abstract: Graphene subject to a spatially uniform, circularly-polarized electric field supports a Floquet spectrum with properties akin to those of a topological insulator, including non-vanishing Chern numbers associated with bulk bands and current-carrying edge states. Transport properties of this system however are complicated by the non-equilibrium occupations of the Floquet states. We address this by considering transport in a two-terminal ribbon geometry for which the leads have well-defined chemical potentials, with an irradiated central scattering region. We demonstrate the presence of edge states, which for infinite mass boundary conditions may be associated with only one of the two valleys. At low frequencies, the bulk DC conductivity near zero energy is shown to be dominated by a series of states with very narrow anticrossings, leading to super-diffusive behavior. For very long ribbons, a ballistic regime emerges in which edge state transport dominates.


title: Temperature Equilibration Behind the Shock Front: an Optical and X-ray Study of RCW 86

abstract: We study the electron-proton temperature equilibration behind several shocks of the RCW 86 supernova remnant. To measure the proton temperature, we use published and new optical spectra, all from different locations on the remnant. For each location, we determine the electron temperature from X-ray spectra, and correct for temperature equilibration between the shock front and the location of the X-ray spectrum. We confirm the result of previous studies that the electron and proton temperatures behind shock fronts are consistent with equilibration for slow shocks and deviate for faster shocks. However, we can not confirm the previously reported trend of the electron temperature to proton temperature ratio of 1/v^2.


title: Using Ontologies for the Design of Data Warehouses

abstract: Obtaining an implementation of a data warehouse is a complex task that forces designers to acquire wide knowledge of the domain, thus requiring a high level of expertise and becoming it a prone-to-fail task. Based on our experience, we have detected a set of situations we have faced up with in real-world projects in which we believe that the use of ontologies will improve several aspects of the design of data warehouses. The aim of this article is to describe several shortcomings of current data warehouse design approaches and discuss the benefit of using ontologies to overcome them. This work is a starting point for discussing the convenience of using ontologies in data warehouse design.


title: Slime mould computes planar shapes

abstract: Computing a polygon defining a set of planar points is a classical problem of modern computational geometry. In laboratory experiments we demonstrate that a concave hull, a connected alpha-shape without holes, of a finite planar set is approximated by slime mould Physarum polycephalum. We represent planar points with sources of long-distance attractants and short-distance repellents and inoculate a piece of plasmodium outside the data set. The plasmodium moves towards the data and envelops it by pronounced protoplasmic tubes.


title: Numerical overcooling in shocks

abstract: We present a study of cooling in radiative shocks simulated with smoothed particle hydrodynamics (SPH) and adaptive mesh refinement codes. We obtain a similarity solution for a shock-tube problem in the presence of radiative cooling, and test how well the solution is reproduced in Gadget and Flash. Shock broadening governed by the details of the numerical scheme (artificial viscosity or Riemann solvers) leads to potentially significant overcooling in both codes. We interpret our findings in terms of a resolution criterion, and apply it to realistic simulations of cosmological accretion shocks onto galaxy haloes, cold accretion and thermal feedback from supernovae or active galactic nuclei. To avoid numerical overcooling of accretion shocks onto haloes that should develop a hot corona requires a particle or cell mass resolution of 10^6 M_sun, which is within reach of current state-of-the-art simulations. At this mass resolution, thermal feedback in the interstellar medium of a galaxy requires temperatures of supernova or AGN driven bubbles to be in excess of 10^7 K at densities of n_H=1.0 cm^-3, in order to avoid spurious suppression of the feedback by numerical overcooling.


title: The discovery and classification of 16 supernovae at high redshifts in ELAIS-S1 - the Stockholm VIMOS Supernova Survey II

abstract: Supernova surveys can be used to study a variety of subjects such as: (i) cosmology through type Ia supernovae (SNe), (ii) star-formation rates through core-collapse SNe, and (iii) supernova properties and their connection to host galaxy characteristics. The Stockholm VIMOS Supernova Survey (SVISS) is a multi-band imaging survey aiming to detect supernovae at redshift ~0.5 and derive thermonuclear and core-collapse supernova rates at high redshift. In this paper we present the supernovae discovered in the survey along with light curves and a photometric classification into thermonuclear and core-collapse types. To detect the supernovae in the VLT/VIMOS multi-epoch images, we used difference imaging and a combination of automatic and manual source detection to minimise the number of spurious detections. Photometry for the found variable sources was obtained and careful simulations were made to estimate correct errors. The light curves were typed using a Bayesian probability method and Monte Carlo simulations were used to study misclassification. We detected 16 supernovae, nine of which had a core-collapse origin and seven had a thermonuclear origin. The estimated misclassification errors are quite small, in the order of 5%, but vary with both redshift and type. The mean redshift of the supernovae is 0.58. Additionally, we found a variable source with a very extended light curve that could possibly be a pair instability supernova.


title: Gravitino Dark Matter and Neutrino Masses in Partial Split Supersymmetry

abstract: Partial Split Supersymmetry with bilinear R-parity violation allows to reproduce all neutrino mass and mixing parameters. The viable dark matter candidate in this model is the gravitino. We study the hypothesis that both possibilities are true: Partial Split Supersymmetry explains neutrino physics and that dark matter is actually composed of gravitinos. Since the gravitino has a small but non-zero decay probability, its decay products could be observed in astrophysical experiments. Combining bounds from astrophysical photon spectra with the bounds coming from the mass matrix in the neutrino sector we derive a stringent upper limit for the allowed gravitino mass. This mass limit is in good agreement with the results of direct dark matter searches.


title: $L_{p}[0,1] \setminus \bigcup\limits_{q>p} L_{q}[0,1]$ is spaceable for every $p>0$

abstract: In this short note we prove the result stated in the title; that is, for every $p>0$ there exists an infinite dimensional closed linear subspace of $L_{p}[0,1]$ every nonzero element of which does not belong to $\bigcup\limits_{q>p} L_{q}[0,1]$. This answers in the positive a question raised in 2010 by R. M. Aron on the spaceability of the above sets (for both, the Banach and quasi-Banach cases). We also complete some recent results from \cite{BDFP} for subsets of sequence spaces.


title: Uncovering the Spectral Energy Distribution in Active Galaxies Using High Ionization Mid-infrared Emission Lines

abstract: The shape of the spectral energy distribution of active galaxies in the EUV--soft X-ray band (13.6 eV to 1 keV) is uncertain because obscuration by dust and gas can hamper our view of the continuum. To investigate the shape of the spectral energy distribution in this energy band, we have generated a set of photoionization models which reproduce the small dispersion found in correlations between high-ionization mid-infrared emission lines in a sample of hard X-ray selected AGN. Our calculations show that a broken power-law continuum model is sufficient to reproduce the [Ne V]14.32 mm/[NeIII], [Ne V]24.32mm/[O IV]25.89mm and [O IV] 25.89mm/[Ne III] ratios, and does not require the addition of a "big bump" EUV model component. We constrain the EUV--soft X-ray slope, alpha_i, to be between 1.5 -- 2.0 and derive a best fit of alpha_i ~ 1.9 for Seyfert 1 galaxies, consistent with previous studies of intermediate redshift quasars. If we assume a blue bump model, most sources in our sample have derived temperatures between T_{BB}=10^{5.18} K to 10^{5.7} K, suggesting that the peak of this component spans a large range of energies extending from ~ lambda 600A to lambda 1900A. In this case, the best fitting peak energy that matches the mid-infrared line ratios of Seyfert 1 galaxies occurs between ~ lambda 700--1000A. Despite the fact that our results do not rule out the presence of an EUV bump, we conclude that our power-law model produces enough photons with energies > 4 Ry to generate the observed amount of mid-infrared emission in our sample of BAT AGN.


title: Resonant elastic soft x-ray scattering in oxygen-ordered YBa_2Cu_3O_{6+delta}

abstract: Static charge-density wave (CDW) and spin-density wave (SDW) order has been convincingly observed in La-based cuprates for some time. However, more recently it has been suggested by quantum oscillation, transport and thermodynamic measurements that density wave order is generic to underdoped cuprates and plays a significant role in YBa_2Cu_3O_{6+delta} (YBCO). We use resonant soft x-ray scattering at the Cu L and O K edges to search for evidence of density wave order in Ortho-II and Ortho-VIII oxygen-ordered YBCO. We report a null result -- no evidence for static CDW order -- in both Ortho-II and Ortho-VIII ordered YBCO. While this does not rule out static CDW order in the CuO_2 planes of YBCO, these measurements place limits on the parameter space (temperature, magnetic field, scattering vector) in which static CDW order may exist. In addition, we present a detailed analysis of the energy and polarization dependence of the Ortho-II superstructure Bragg reflection [0.5 0 0] at the Cu L edge. The intensity of this peak, which is due to the valence modulations of Cu in the chain layer, is compared with calculations using atomic scattering form factors deduced from x-ray absorption measurements. The calculated energy and polarization dependence of the scattering intensity is shown to agree very well with the measurement, validating the approach and providing a framework for analyzing future resonant soft x-ray scattering measurements.


title: Molecular Dynamics Simulation Study of Interaction between Model Rough Hydrophobic Surfaces

abstract: We study some aspects of hydrophobic interaction between molecular rough and flexible model surfaces. The model we use in this work is based on a model we used previously (Eun, C.; Berkowitz, M. L. J. Phys. Chem. B 2009, 113, 13222-13228), when we studied the interaction between model patches of lipid membranes. Our original model consisted of two graphene plates with attached polar headgroups; the plates were immersed in a water bath. The interaction between such plates can be considered as an example of a hydrophilic interaction. In the present work we modify our previous model by removing the charge from the zwitterionic headgroups. As a result of this procedure, the plate character changes; it becomes hydrophobic. By separating the total interaction (or potential of mean force, PMF) between plates into the direct and the water-mediated interactions we observe that the latter changes from repulsive to attractive, clearly emphasizing the important role of water as a medium. We also investigate the effect of roughness and flexibility of the headgroups on the interaction between plates and observe that roughness enhances the character of the hydrophobic interaction. The presence of a dewetting transition in a confined space between charge-removed plates confirms that the interaction between plates is strongly hydrophobic. In addition, we notice that there is a shallow local minimum in the PMF in case of charge-removed plates. We find that this minimum is associated with the configurational changes that flexible headgroups undergo, as the two plates are brought together.


title: Relative locality: A deepening of the relativity principle

abstract: We describe a recently introduced principle of relative locality which we propose governs a regime of quantum gravitational phenomena accessible to experimental investigation. This regime comprises phenomena in which $\hbar$ and $G_N$ may be neglected, while their ratio, the Planck mass $M_p =\sqrt{\hbar / G_N}$, is important. We propose that $M_p$ governs the scale at which momentum space may have a curved geometry. We find that there are striking consequences for the concept of locality. The description of events in spacetime now depends on the energy used to probe it. But there remains an invariant description of physics in phase space. There is furthermore a reasonable expectation that the geometry of momentum space can be measured experimentally using astrophysical observations.


title: The nonlinear power spectrum in clustering quintessence cosmologies

abstract: We study the nonlinear evolution of density perturbations in cosmologies where the late-time accelerated expansion is driven by a quintessence field with vanishing speed of sound. For these models matter and quintessence perturbations are comoving and it is possible to write a single continuity equation for the total density fluctuations given by a weighted sum of the two components. Including the Euler equation for the common velocity field we solve the evolution equations for the nonlinear, total density power spectrum in the Time-Renormalization Group approach. In fact any cosmological observable is directly related by gravity only to the total density perturbations, with the two components being individually unobservable. We show that the clustering of quintessence perturbations induces small corrections with respect to the nonlinear evolution of power spectrum in smooth quintessence models described by the same equation of state. Such small corrections, however, contrast with the large effect of a vanishing speed of sound on the linear growth function at low redshift. For this reason, models with the same normalization of the linear density power spectrum can present significantly different nonlinear corrections depending on the value of the sound speed. Although such differences vanish in the w \rightarrow -1, {\Lambda}CDM limit, we argue that the relation between linear and nonlinear growth of structures should be properly taken into account in constraining models with inhomogeneous dark energy.


title: Compact Optical Counterparts of Ultraluminous X-ray Sources

abstract: Using archival Hubble Space Telescope (HST) imaging data, we report the multiband photometric properties of 13 ultraluminous X-ray sources (ULXs) that have a unique compact optical counterpart. Both magnitude and color variation are detected at time scales of days to years. The optical color, variability, and X-ray to optical flux ratio indicate that the optical emission of most ULXs is dominated by X-ray reprocessing on the disk, similar to that of low mass X-ray binaries. For most sources, the optical spectrum is a power-law, $F_{\nu} \propto \nu^{\alpha}$ with $\alpha$ in the range 1.0 to 2.0 and the optically emitting region has a size on the order of 1e12 cm. Exceptions are NGC 2403 X-1 and M83 IXO 82, which show optical spectra consistent with direct emission from a standard thin disk, M101 ULX-1 and M81 ULS1, which have X-ray to optical flux ratios more similar to high-mass X-ray binaries, and IC 342 X-1, in which the optical light may be dominated by the companion star. Inconsistent extinction between the optical counterpart of NGC 5204 X-1 and the nearby optical nebulae suggests that they may be unrelated.


title: Non-geometric Fluxes, Asymmetric Strings and Nonassociative Geometry

abstract: We study closed bosonic strings propagating both in a flat background with constant H-flux and in its T-dual configurations. We define a conformal field theory capturing linear effects in the flux and compute scattering amplitudes of tachyons, where the Rogers dilogarithm plays a prominent role. For the scattering of four tachyons, a fluxed version of the Virasoro-Shapiro amplitude is derived and its pole structure is analyzed. In the case of an R-flux background obtained after three T-dualities, we find indications for a nonassociative target-space structure which can be described in terms of a deformed tri-product. Remarkably, this product is compatible with crossing symmetry of conformal correlation functions. We finally argue that the R-flux background flows to an asymmetric CFT.


title: Modelling stellar populations at high redshift

abstract: Stellar populations carry information about the formation of galaxies and their evolution up to the present epoch. A wealth of observational data are available nowadays, which are analysed with stellar population models in order to obtain key properties such as ages, star formation histories, stellar masses. Differences in the models and/or in the assumptions regarding the star formation history affect the derived properties as much as differences in the data. I shall review the interpretation of high-redshift galaxy data from a model perspective. While data quality dominates galaxy analysis at the highest possible redshifts (z>5), population modelling effects play the major part at lower redshifts. In particular, I discuss the cases of both star-forming galaxies at the peak of the cosmic star formation history as well as passive galaxies at redshift below 1 that are often used as cosmological probes. Remarks on the bridge between low and high-z massive galaxies conclude the contribution.


title: Non-autonomous 2-periodic Gumovski-Mira difference equations

abstract: We consider two types of non-autonomous 2-periodic Gumovski-Mira difference equations. We show that while the corresponding autonomous recurrences are conjugated, the behavior of the sequences generated by the 2-periodic ones differ dramatically: in one case the behavior of the sequences is simple (integrable) and in the other case it is much more complicated (chaotic). We also present a global study of the integrable case that includes which periods appear for the recurrence.


title: Instabilities in magnetized spherical Couette flow

abstract: We report 3D numerical simulations of the flow of an electrically conducting fluid in a spherical shell when a magnetic field is applied. Different spherical Couette configurations are investigated, by varying the rotation ratio between the inner and the outer sphere, the geometry of the imposed field, and the magnetic boundary conditions on the inner sphere. Either a Stewartson layer or a Shercliff layer, accompanied by a radial jet, can be generated depending on the rotation speeds and the magnetic field strength, and various non-axisymmetric destabilizations of the flow are observed. We show that instabilities arising from the presence of boundaries present striking similarities with the magnetorotational instability (MRI). To this end, we compare our numerical results to experimental observations of the Maryland experiment, who claimed to observe MRI in a similar setup.


title: Fluctuations of Matrix Entries of Regular Functions of Sample Covariance Random Matrices

abstract: We extend the results about the fluctuations of the matrix entries of regular functions of Wigner matrices to the case of sample covariance random matrices.


title: Sparse Non Gaussian Component Analysis by Semidefinite Programming

abstract: Sparse non-Gaussian component analysis (SNGCA) is an unsupervised method of extracting a linear structure from a high dimensional data based on estimating a low-dimensional non-Gaussian data component. In this paper we discuss a new approach to direct estimation of the projector on the target space based on semidefinite programming which improves the method sensitivity to a broad variety of deviations from normality. We also discuss the procedures which allows to recover the structure when its effective dimension is unknown.


title: Bayesian Sparsity-Path-Analysis of Genetic Association Signal using Generalized t Priors

abstract: We explore the use of generalized t priors on regression coefficients to help understand the nature of association signal within "hit regions" of genome-wide association studies. The particular generalized t distribution we adopt is a Student distribution on the absolute value of its argument. For low degrees of freedom we show that the generalized t exhibits 'sparsity-prior' properties with some attractive features over other common forms of sparse priors and includes the well known double-exponential distribution as the degrees of freedom tends to infinity. We pay particular attention to graphical representations of posterior statistics obtained from sparsity-path-analysis (SPA) where we sweep over the setting of the scale (shrinkage / precision) parameter in the prior to explore the space of posterior models obtained over a range of complexities, from very sparse models with all coefficient distributions heavily concentrated around zero, to models with diffuse priors and coefficients distributed around their maximum likelihood estimates. The SPA plots are akin to LASSO plots of maximum a posteriori (MAP) estimates but they characterise the complete marginal posterior distributions of the coefficients plotted as a function of the precision of the prior. Generating posterior distributions over a range of prior precisions is computationally challenging but naturally amenable to sequential Monte Carlo (SMC) algorithms indexed on the scale parameter. We show how SMC simulation on graphic-processing-units (GPUs) provides very efficient inference for SPA. We also present a scale-mixture representation of the generalized t prior that leads to an EM algorithm to obtain MAP estimates should only these be required.


title: Detection and Imaging of the Crab Nebula with the Nuclear Compton Telescope

abstract: The Nuclear Compton Telescope (NCT) is a balloon-borne Compton telescope designed for the study of astrophysical sources in the soft gamma-ray regime (200 keV--20 MeV). NCT's ten high-purity germanium crossed-strip detectors measure the deposited energies and three-dimensional positions of gamma-ray interactions in the sensitive volume, and this information is used to restrict the initial photon to a circle on the sky using the Compton scatter technique. Thus NCT is able to perform spectroscopy, imaging, and polarization analysis on soft gamma-ray sources. NCT is one of the next generation of Compton telescopes --- so-called compact Compton telescopes (CCTs) --- which can achieve effective areas comparable to COMPTEL's with an instrument that is a fraction of the size. The Crab Nebula was the primary target for the second flight of the NCT instrument, which occurred on 17--18 May 2009 in Fort Sumner, New Mexico. Analysis of 29.3 ks of data from the flight reveals an image of the Crab at a significance of 4-sigma. This is the first reported detection of an astrophysical source by a CCT.