title: Multidimensional counting grids: Inferring word order from disordered bags of words

abstract: Models of bags of words typically assume topic mixing so that the words in a single bag come from a limited number of topics. We show here that many sets of bag of words exhibit a very different pattern of variation than the patterns that are efficiently captured by topic mixing. In many cases, from one bag of words to the next, the words disappear and new ones appear as if the theme slowly and smoothly shifted across documents (providing that the documents are somehow ordered). Examples of latent structure that describe such ordering are easily imagined. For example, the advancement of the date of the news stories is reflected in a smooth change over the theme of the day as certain evolving news stories fall out of favor and new events create new stories. Overlaps among the stories of consecutive days can be modeled by using windows over linearly arranged tight distributions over words. We show here that such strategy can be extended to multiple dimensions and cases where the ordering of data is not readily obvious. We demonstrate that this way of modeling covariation in word occurrences outperforms standard topic models in classification and prediction tasks in applications in biology, text modeling and computer vision.


title: Partial Order MCMC for Structure Discovery in Bayesian Networks

abstract: We present a new Markov chain Monte Carlo method for estimating posterior probabilities of structural features in Bayesian networks. The method draws samples from the posterior distribution of partial orders on the nodes; for each sampled partial order, the conditional probabilities of interest are computed exactly. We give both analytical and empirical results that suggest the superiority of the new method compared to previous methods, which sample either directed acyclic graphs or linear orders on the nodes.


title: A Geometric Traversal Algorithm for Reward-Uncertain MDPs

abstract: Markov decision processes (MDPs) are widely used in modeling decision making problems in stochastic environments. However, precise specification of the reward functions in MDPs is often very difficult. Recent approaches have focused on computing an optimal policy based on the minimax regret criterion for obtaining a robust policy under uncertainty in the reward function. One of the core tasks in computing the minimax regret policy is to obtain the set of all policies that can be optimal for some candidate reward function. In this paper, we propose an efficient algorithm that exploits the geometric properties of the reward function associated with the policies. We also present an approximate version of the method for further speed up. We experimentally demonstrate that our algorithm improves the performance by orders of magnitude.


title: Iterated risk measures for risk-sensitive Markov decision processes with discounted cost

abstract: We demonstrate a limitation of discounted expected utility, a standard approach for representing the preference to risk when future cost is discounted. Specifically, we provide an example of the preference of a decision maker that appears to be rational but cannot be represented with any discounted expected utility. A straightforward modification to discounted expected utility leads to inconsistent decision making over time. We will show that an iterated risk measure can represent the preference that cannot be represented by any discounted expected utility and that the decisions based on the iterated risk measure are consistent over time.


title: Price Updating in Combinatorial Prediction Markets with Bayesian Networks

abstract: To overcome the #P-hardness of computing/updating prices in logarithm market scoring rule-based (LMSR-based) combinatorial prediction markets, Chen et al. [5] recently used a simple Bayesian network to represent the prices of securities in combinatorial predictionmarkets for tournaments, and showed that two types of popular securities are structure preserving. In this paper, we significantly extend this idea by employing Bayesian networks in general combinatorial prediction markets. We reveal a very natural connection between LMSR-based combinatorial prediction markets and probabilistic belief aggregation,which leads to a complete characterization of all structure preserving securities for decomposable network structures. Notably, the main results by Chen et al. [5] are corollaries of our characterization. We then prove that in order for a very basic set of securities to be structure preserving, the graph of the Bayesian network must be decomposable. We also discuss some approximation techniques for securities that are not structure preserving.


title: Identifiability of Causal Graphs using Functional Models

abstract: This work addresses the following question: Under what assumptions on the data generating process can one infer the causal graph from the joint distribution? The approach taken by conditional independence-based causal discovery methods is based on two assumptions: the Markov condition and faithfulness. It has been shown that under these assumptions the causal graph can be identified up to Markov equivalence (some arrows remain undirected) using methods like the PC algorithm. In this work we propose an alternative by defining Identifiable Functional Model Classes (IFMOCs). As our main theorem we prove that if the data generating process belongs to an IFMOC, one can identify the complete causal graph. To the best of our knowledge this is the first identifiability result of this kind that is not limited to linear functional relationships. We discuss how the IFMOC assumption and the Markov and faithfulness assumptions relate to each other and explain why we believe that the IFMOC assumption can be tested more easily on given data. We further provide a practical algorithm that recovers the causal graph from finitely many data; experiments on simulated data support the theoretical findings.


title: Nonparametric Divergence Estimation with Applications to Machine Learning on Distributions

abstract: Low-dimensional embedding, manifold learning, clustering, classification, and anomaly detection are among the most important problems in machine learning. The existing methods usually consider the case when each instance has a fixed, finite-dimensional feature representation. Here we consider a different setting. We assume that each instance corresponds to a continuous probability distribution. These distributions are unknown, but we are given some i.i.d. samples from each distribution. Our goal is to estimate the distances between these distributions and use these distances to perform low-dimensional embedding, clustering/classification, or anomaly detection for the distributions. We present estimation algorithms, describe how to apply them for machine learning tasks on distributions, and show empirical results on synthetic data, real word images, and astronomical data sets.


title: Compressed Inference for Probabilistic Sequential Models

abstract: Hidden Markov models (HMMs) and conditional random fields (CRFs) are two popular techniques for modeling sequential data. Inference algorithms designed over CRFs and HMMs allow estimation of the state sequence given the observations. In several applications, estimation of the state sequence is not the end goal; instead the goal is to compute some function of it. In such scenarios, estimating the state sequence by conventional inference techniques, followed by computing the functional mapping from the estimate is not necessarily optimal. A more formal approach is to directly infer the final outcome from the observations. In particular, we consider the specific instantiation of the problem where the goal is to find the state trajectories without exact transition points and derive a novel polynomial time inference algorithm that outperforms vanilla inference techniques. We show that this particular problem arises commonly in many disparate applications and present experiments on three of them: (1) Toy robot tracking; (2) Single stroke character recognition; (3) Handwritten word recognition.


title: Fast MCMC sampling for Markov jump processes and continuous time Bayesian networks

abstract: Markov jump processes and continuous time Bayesian networks are important classes of continuous time dynamical systems. In this paper, we tackle the problem of inferring unobserved paths in these models by introducing a fast auxiliary variable Gibbs sampler. Our approach is based on the idea of uniformization, and sets up a Markov chain over paths by sampling a finite set of virtual jump times and then running a standard hidden Markov model forward filtering-backward sampling algorithm over states at the set of extant and virtual jump times. We demonstrate significant computational benefits over a state-of-the-art Gibbs sampler on a number of continuous time Bayesian networks.


title: New Probabilistic Bounds on Eigenvalues and Eigenvectors of Random Kernel Matrices

abstract: Kernel methods are successful approaches for different machine learning problems. This success is mainly rooted in using feature maps and kernel matrices. Some methods rely on the eigenvalues/eigenvectors of the kernel matrix, while for other methods the spectral information can be used to estimate the excess risk. An important question remains on how close the sample eigenvalues/eigenvectors are to the population values. In this paper, we improve earlier results on concentration bounds for eigenvalues of general kernel matrices. For distance and inner product kernel functions, e.g. radial basis functions, we provide new concentration bounds, which are characterized by the eigenvalues of the sample covariance matrix. Meanwhile, the obstacles for sharper bounds are accounted for and partially addressed. As a case study, we derive a concentration inequality for sample kernel target-alignment.


title: Symbolic Dynamic Programming for Discrete and Continuous State MDPs

abstract: Many real-world decision-theoretic planning problems can be naturally modeled with discrete and continuous state Markov decision processes (DC-MDPs). While previous work has addressed automated decision-theoretic planning for DCMDPs, optimal solutions have only been defined so far for limited settings, e.g., DC-MDPs having hyper-rectangular piecewise linear value functions. In this work, we extend symbolic dynamic programming (SDP) techniques to provide optimal solutions for a vastly expanded class of DCMDPs. To address the inherent combinatorial aspects of SDP, we introduce the XADD - a continuous variable extension of the algebraic decision diagram (ADD) - that maintains compact representations of the exact value function. Empirically, we demonstrate an implementation of SDP with XADDs on various DC-MDPs, showing the first optimal automated solutions to DCMDPs with linear and nonlinear piecewise partitioned value functions and showing the advantages of constraint-based pruning for XADDs.


title: An Efficient Algorithm for Computing Interventional Distributions in Latent Variable Causal Models

abstract: Probabilistic inference in graphical models is the task of computing marginal and conditional densities of interest from a factorized representation of a joint probability distribution. Inference algorithms such as variable elimination and belief propagation take advantage of constraints embedded in this factorization to compute such densities efficiently. In this paper, we propose an algorithm which computes interventional distributions in latent variable causal models represented by acyclic directed mixed graphs(ADMGs). To compute these distributions efficiently, we take advantage of a recursive factorization which generalizes the usual Markov factorization for DAGs and the more recent factorization for ADMGs. Our algorithm can be viewed as a generalization of variable elimination to the mixed graph case. We show our algorithm is exponential in the mixed graph generalization of treewidth.


title: Adjustment Criteria in Causal Diagrams: An Algorithmic Perspective

abstract: Identifying and controlling bias is a key problem in empirical sciences. Causal diagram theory provides graphical criteria for deciding whether and how causal effects can be identified from observed (nonexperimental) data by covariate adjustment. Here we prove equivalences between existing as well as new criteria for adjustment and we provide a new simplified but still equivalent notion of d-separation. These lead to efficient algorithms for two important tasks in causal diagram analysis: (1) listing minimal covariate adjustments (with polynomial delay); and (2) identifying the subdiagram involved in biasing paths (in linear time). Our results improve upon existing exponential-time solutions for these problems, enabling users to assess the effects of covariate adjustment on diagrams with tens to hundreds of variables interactively in real time.


title: Learning mixed graphical models from data with p larger than n

abstract: Structure learning of Gaussian graphical models is an extensively studied problem in the classical multivariate setting where the sample size n is larger than the number of random variables p, as well as in the more challenging setting when p>>n. However, analogous approaches for learning the structure of graphical models with mixed discrete and continuous variables when p>>n remain largely unexplored. Here we describe a statistical learning procedure for this problem based on limited-order correlations and assess its performance with synthetic and real data.


title: Robust learning Bayesian networks for prior belief

abstract: Recent reports have described that learning Bayesian networks are highly sensitive to the chosen equivalent sample size (ESS) in the Bayesian Dirichlet equivalence uniform (BDeu). This sensitivity often engenders some unstable or undesirable results. This paper describes some asymptotic analyses of BDeu to explain the reasons for the sensitivity and its effects. Furthermore, this paper presents a proposal for a robust learning score for ESS by eliminating the sensitive factors from the approximation of log-BDeu.


title: Distributed Anytime MAP Inference

abstract: We present a distributed anytime algorithm for performing MAP inference in graphical models. The problem is formulated as a linear programming relaxation over the edges of a graph. The resulting program has a constraint structure that allows application of the Dantzig-Wolfe decomposition principle. Subprograms are defined over individual edges and can be computed in a distributed manner. This accommodates solutions to graphs whose state space does not fit in memory. The decomposition master program is guaranteed to compute the optimal solution in a finite number of iterations, while the solution converges monotonically with each iteration. Formulating the MAP inference problem as a linear program allows additional (global) constraints to be defined; something not possible with message passing algorithms. Experimental results show that our algorithm's solution quality outperforms most current algorithms and it scales well to large problems.


title: The Structure of Signals: Causal Interdependence Models for Games of Incomplete Information

abstract: Traditional economic models typically treat private information, or signals, as generated from some underlying state. Recent work has explicated alternative models, where signals correspond to interpretations of available information. We show that the difference between these formulations can be sharply cast in terms of causal dependence structure, and employ graphical models to illustrate the distinguishing characteristics. The graphical representation supports inferences about signal patterns in the interpreted framework, and suggests how results based on the generated model can be extended to more general situations. Specific insights about bidding games in classical auction mechanisms derive from qualitative graphical models.


title: Sparse matrix-variate Gaussian process blockmodels for network modeling

abstract: We face network data from various sources, such as protein interactions and online social networks. A critical problem is to model network interactions and identify latent groups of network nodes. This problem is challenging due to many reasons. For example, the network nodes are interdependent instead of independent of each other, and the data are known to be very noisy (e.g., missing edges). To address these challenges, we propose a new relational model for network data, Sparse Matrix-variate Gaussian process Blockmodel (SMGB). Our model generalizes popular bilinear generative models and captures nonlinear network interactions using a matrix-variate Gaussian process with latent membership variables. We also assign sparse prior distributions on the latent membership variables to learn sparse group assignments for individual network nodes. To estimate the latent variables efficiently from data, we develop an efficient variational expectation maximization method. We compared our approaches with several state-of-the-art network models on both synthetic and real-world network datasets. Experimental results demonstrate SMGBs outperform the alternative approaches in terms of discovering latent classes or predicting unknown interactions.


title: Hierarchical Maximum Margin Learning for Multi-Class Classification

abstract: Due to myriads of classes, designing accurate and efficient classifiers becomes very challenging for multi-class classification. Recent research has shown that class structure learning can greatly facilitate multi-class learning. In this paper, we propose a novel method to learn the class structure for multi-class classification problems. The class structure is assumed to be a binary hierarchical tree. To learn such a tree, we propose a maximum separating margin method to determine the child nodes of any internal node. The proposed method ensures that two classgroups represented by any two sibling nodes are most separable. In the experiments, we evaluate the accuracy and efficiency of the proposed method over other multi-class classification methods on real world large-scale problems. The results show that the proposed method outperforms benchmark methods in terms of accuracy for most datasets and performs comparably with other class structure learning methods in terms of efficiency for all datasets.


title: Tightening MRF Relaxations with Planar Subproblems

abstract: We describe a new technique for computing lower-bounds on the minimum energy configuration of a planar Markov Random Field (MRF). Our method successively adds large numbers of constraints and enforces consistency over binary projections of the original problem state space. These constraints are represented in terms of subproblems in a dual-decomposition framework that is optimized using subgradient techniques. The complete set of constraints we consider enforces cycle consistency over the original graph. In practice we find that the method converges quickly on most problems with the addition of a few subproblems and outperforms existing methods for some interesting classes of hard potentials.


title: Rank/Norm Regularization with Closed-Form Solutions: Application to Subspace Clustering

abstract: When data is sampled from an unknown subspace, principal component analysis (PCA) provides an effective way to estimate the subspace and hence reduce the dimension of the data. At the heart of PCA is the Eckart-Young-Mirsky theorem, which characterizes the best rank k approximation of a matrix. In this paper, we prove a generalization of the Eckart-Young-Mirsky theorem under all unitarily invariant norms. Using this result, we obtain closed-form solutions for a set of rank/norm regularized problems, and derive closed-form solutions for a general class of subspace clustering problems (where data is modelled by unions of unknown subspaces). From these results we obtain new theoretical insights and promising experimental results.


title: Measuring the Hardness of Stochastic Sampling on Bayesian Networks with Deterministic Causalities: the k-Test

abstract: Approximate Bayesian inference is NP-hard. Dagum and Luby defined the Local Variance Bound (LVB) to measure the approximation hardness of Bayesian inference on Bayesian networks, assuming the networks model strictly positive joint probability distributions, i.e. zero probabilities are not permitted. This paper introduces the k-test to measure the approximation hardness of inference on Bayesian networks with deterministic causalities in the probability distribution, i.e. when zero conditional probabilities are permitted. Approximation by stochastic sampling is a widely-used inference method that is known to suffer from inefficiencies due to sample rejection. The k-test predicts when rejection rates of stochastic sampling a Bayesian network will be low, modest, high, or when sampling is intractable.


title: Risk Bounds for Infinitely Divisible Distribution

abstract: In this paper, we study the risk bounds for samples independently drawn from an infinitely divisible (ID) distribution. In particular, based on a martingale method, we develop two deviation inequalities for a sequence of random variables of an ID distribution with zero Gaussian component. By applying the deviation inequalities, we obtain the risk bounds based on the covering number for the ID distribution. Finally, we analyze the asymptotic convergence of the risk bound derived from one of the two deviation inequalities and show that the convergence rate of the bound is faster than the result for the generic i.i.d. empirical process (Mendelson, 2003).


title: Kernel-based Conditional Independence Test and Application in Causal Discovery

abstract: Conditional independence testing is an important problem, especially in Bayesian network learning and causal discovery. Due to the curse of dimensionality, testing for conditional independence of continuous variables is particularly challenging. We propose a Kernel-based Conditional Independence test (KCI-test), by constructing an appropriate test statistic and deriving its asymptotic distribution under the null hypothesis of conditional independence. The proposed method is computationally efficient and easy to implement. Experimental results show that it outperforms other methods, especially when the conditioning set is large or the sample size is not very large, in which case other methods encounter difficulties.


title: Smoothing Multivariate Performance Measures

abstract: A Support Vector Method for multivariate performance measures was recently introduced by Joachims (2005). The underlying optimization problem is currently solved using cutting plane methods such as SVM-Perf and BMRM. One can show that these algorithms converge to an eta accurate solution in O(1/Lambda*e) iterations, where lambda is the trade-off parameter between the regularizer and the loss function. We present a smoothing strategy for multivariate performance scores, in particular precision/recall break-even point and ROCArea. When combined with Nesterov's accelerated gradient algorithm our smoothing strategy yields an optimization algorithm which converges to an eta accurate solution in O(min{1/e,1/sqrt(lambda*e)}) iterations. Furthermore, the cost per iteration of our scheme is the same as that of SVM-Perf and BMRM. Empirical evaluation on a number of publicly available datasets shows that our method converges significantly faster than cutting plane methods without sacrificing generalization ability.


title: Belief Propagation by Message Passing in Junction Trees: Computing Each Message Faster Using GPU Parallelization

abstract: Compiling Bayesian networks (BNs) to junction trees and performing belief propagation over them is among the most prominent approaches to computing posteriors in BNs. However, belief propagation over junction tree is known to be computationally intensive in the general case. Its complexity may increase dramatically with the connectivity and state space cardinality of Bayesian network nodes. In this paper, we address this computational challenge using GPU parallelization. We develop data structures and algorithms that extend existing junction tree techniques, and specifically develop a novel approach to computing each belief propagation message in parallel. We implement our approach on an NVIDIA GPU and test it using BNs from several applications. Experimentally, we study how junction tree parameters affect parallelization opportunities and hence the performance of our algorithm. We achieve speedups ranging from 0.68 to 9.18 for the BNs studied.


title: Sparse Topical Coding

abstract: We present sparse topical coding (STC), a non-probabilistic formulation of topic models for discovering latent representations of large collections of data. Unlike probabilistic topic models, STC relaxes the normalization constraint of admixture proportions and the constraint of defining a normalized likelihood function. Such relaxations make STC amenable to: 1) directly control the sparsity of inferred representations by using sparsity-inducing regularizers; 2) be seamlessly integrated with a convex error function (e.g., SVM hinge loss) for supervised learning; and 3) be efficiently learned with a simply structured coordinate descent algorithm. Our results demonstrate the advantages of STC and supervised MedSTC on identifying topical meanings of words and improving classification accuracy and time efficiency.


title: Testing whether linear equations are causal: A free probability theory approach

abstract: We propose a method that infers whether linear relations between two high-dimensional variables X and Y are due to a causal influence from X to Y or from Y to X. The earlier proposed so-called Trace Method is extended to the regime where the dimension of the observed variables exceeds the sample size. Based on previous work, we postulate conditions that characterize a causal relation between X and Y. Moreover, we describe a statistical test and argue that both causal directions are typically rejected if there is a common cause. A full theoretical analysis is presented for the deterministic case but our approach seems to be valid for the noisy case, too, for which we additionally present an approach based on a sparsity constraint. The discussed method yields promising results for both simulated and real world data.


title: A New Class of Four-Dimensional N=1 Supergravity with Non-minimal Derivative Couplings

abstract: In the N=1 four-dimensional new-minimal supergravity framework, we supersymmetrise the coupling of the scalar kinetic term to the Einstein tensor. This coupling, although introduces a non-minimal derivative interaction of curvature to matter, it does not introduce harmful higher-derivatives. For this construction, we employ off-shell chiral and real linear multiplets. Physical scalars are accommodated in the chiral multiplet whereas curvature resides in a linear one.


title: Phase-Space analysis of Teleparallel Dark Energy

abstract: We perform a detailed dynamical analysis of the teleparallel dark energy scenario, which is based on the teleparallel equivalent of General Relativity, in which one adds a canonical scalar field, allowing also for a nonminimal coupling with gravity. We find that the universe can result in the quintessence-like, dark-energy-dominated solution, or to the stiff dark-energy late-time attractor, similarly to standard quintessence. However, teleparallel dark energy possesses an additional late-time solution, in which dark energy behaves like a cosmological constant, independently of the specific values of the model parameters. Finally, during the evolution the dark energy equation-of-state parameter can be either above or below -1, offering a good description for its observed dynamical behavior and its stabilization close to the cosmological-constant value.


title: Graphical Models for Bandit Problems

abstract: We introduce a rich class of graphical models for multi-armed bandit problems that permit both the state or context space and the action space to be very large, yet succinctly specify the payoffs for any context-action pair. Our main result is an algorithm for such models whose regret is bounded by the number of parameters and whose running time depends only on the treewidth of the graph substructure induced by the action space.


title: Is Quasar Optical Variability a Damped Random Walk?

abstract: The damped random walk (DRW) model is increasingly used to model the variability in quasar optical light curves, but it is still uncertain whether the DRW model provides an adequate description of quasar optical variability across all time scales. Using a sample of OGLE quasar light curves, we consider four modifications to the DRW model by introducing additional parameters into the covariance function to search for deviations from the DRW model on both short and long time scales. We find good agreement with the DRW model on time scales that are well sampled by the data (from a month to a few years), possibly with some intrinsic scatter in the additional parameters, but this conclusion depends on the statistical test employed and is sensitive to whether the estimates of the photometric errors are correct to within ~10%. On very short time scales (below a few months), we see some evidence of the existence of a cutoff below which the correlation is stronger than the DRW model, echoing the recent finding of Mushotzky et al. (2011) using quasar light curves from Kepler. On very long time scales (> a few years), the light curves do not constrain models well, but are consistent with the DRW model.


title: Constraints on the Lifetimes of Disks Resulting from Tidally Destroyed Rocky Planetary Bodies

abstract: Spitzer IRAC observations of 15 metal-polluted white dwarfs reveal infrared excesses in the spectral energy distributions of HE 0110-5630, GD 61, and HE 1349-2305. All three of these stars have helium-dominated atmospheres, and their infrared emissions are consistent with warm dust produced by the tidal destruction of (minor) planetary bodies. This study brings the number of metal-polluted, helium and hydrogen atmosphere white dwarfs surveyed with IRAC to 53 and 38 respectively. It also nearly doubles the number of metal-polluted helium-rich white dwarfs found to have closely orbiting dust by Spitzer. From the increased statistics for both atmospheric types with circumstellar dust, we derive a typical disk lifetime of log[t_{disk} (yr)] = 5.6+-1.1 (ranging from 3*10^4 - 5*10^6 yr). This assumes a relatively constant rate of accretion over the timescale where dust persists, which is uncertain. We find that the fraction of highly metal-polluted helium-rich white dwarfs that have an infrared excess detected by Spitzer is only 23 per cent, compared to 48 per cent for metal-polluted hydrogen-rich white dwarfs, and we conclude from this difference that the typical lifetime of dusty disks is somewhat shorter than the diffusion time scales of helium-rich white dwarf. We also find evidence for higher time-averaged accretion rates onto helium-rich stars compared to the instantaneous accretion rates onto hydrogen-rich stars; this is an indication that our picture of evolved star-planetary system interactions is incomplete. We discuss some speculative scenarios that can explain the observations.


title: Using Kepler transit observations to measure stellar spot belt migration rates

abstract: Planetary transits provide a unique opportunity to investigate the surface distributions of star spots. Our aim is to determine if, with continuous observation (such as the data that will be provided by the Kepler mission), we can in addition measure the rate of drift of the spot belts. We begin by simulating magnetic cycles suitable for the Sun and more active stars, incorporating both flux emergence and surface transport. This provides the radial magnetic field distribution on the stellar surface as a function of time. We then model the transit of a planet whose orbital axis is misaligned with the stellar rotation axis. Such a planet could occult spots at a range of latitudes. This allows us to complete the forward modelling of the shape of the transit lightcurve. We then attempt the inverse problem of recovering spot locations from the transit alone. From this we determine if transit lightcurves can be used to measure spot belt locations as a function of time. We find that for low-activity stars such as the Sun, the 3.5 year Kepler window is insufficient to determine this drift rate. For more active stars, it may be difficult to distinguish subtle differences in the nature of flux emergence, such as the degree of overlap of the "butterfly wings". The rate and direction of drift of the spot belts can however be determined for these stars. This would provide a critical test of dynamo theory.


title: Measuring X-ray variability in faint/sparsely-sampled AGN

abstract: We study the statistical properties of the Normalized Excess Variance of variability process characterized by a red-noise power spectral density (PSD), as the case of Active Galactic Nuclei (AGN). We perform Monte Carlo simulations of lightcurves, assuming both a continuous and a sparse sampling pattern and various signal-to-noise (S/N) ratios. We show that the normalized excess variance is a biased estimate of the variance even in the case of continuously sampled lightcurves. The bias depends on the PSD slope and on the sampling pattern, but not on the S/N ratio. We provide a simple formula to account for the bias, which yields unbiased estimates with an accuracy better than 15%. We show that the normalized excess variance estimates based on single lightcurves (especially for sparse sampling and S/N less than 3) are highly uncertain (even if corrected for bias) and we propose instead the use of an ensemble estimate, based on multiple lightcurves of the same object, or on the use of lightcurves of many objects. These estimates have symmetric distributions, known errors, and can also be corrected for biases.We use our results to estimate the ability to measure the intrinsic source variability in current data, and show that they could also be useful in the planning of the observing strategy of future surveys such as those provided by X-ray missions studying distant and/or faint AGN populations and, more in general, in the estimation of the variability amplitude of sources that will result from future surveys such as Pan-STARRS, and LSST.


title: The XMM Cluster Survey: The interplay between the brightest cluster galaxy and the intra-cluster medium via AGN feedback

abstract: Using a sample of 123 X-ray clusters and groups drawn from the XMM-Cluster Survey first data release, we investigate the interplay between the brightest cluster galaxy (BCG), its black hole, and the intra-cluster/group medium (ICM). It appears that for groups and clusters with a BCG likely to host significant AGN feedback, gas cooling dominates in those with Tx > 2 keV while AGN feedback dominates below. This may be understood through the sub-unity exponent found in the scaling relation we derive between the BCG mass and cluster mass over the halo mass range 10^13 < M500 < 10^15Msol and the lack of correlation between radio luminosity and cluster mass, such that BCG AGN in groups can have relatively more energetic influence on the ICM. The Lx - Tx relation for systems with the most massive BCGs, or those with BCGs co-located with the peak of the ICM emission, is steeper than that for those with the least massive and most offset, which instead follows self-similarity. This is evidence that a combination of central gas cooling and powerful, well fuelled AGN causes the departure of the ICM from pure gravitational heating, with the steepened relation crossing self-similarity at Tx = 2 keV. Importantly, regardless of their black hole mass, BCGs are more likely to host radio-loud AGN if they are in a massive cluster (Tx > 2 keV) and again co-located with an effective fuel supply of dense, cooling gas. This demonstrates that the most massive black holes appear to know more about their host cluster than they do about their host galaxy. The results lead us to propose a physically motivated, empirical definition of 'cluster' and 'group', delineated at 2 keV.


title: Carbon Detection in Early-Time Optical Spectra of Type Ia Supernovae

abstract: While O is often seen in spectra of Type Ia supernovae (SNe Ia) as both unburned fuel and a product of C burning, C is only occasionally seen at the earliest times, and it represents the most direct way of investigating primordial white dwarf material and its relation to SN Ia explosion scenarios and mechanisms. In this paper, we search for C absorption features in 188 optical spectra of 144 low-redshift (z < 0.1) SNe Ia with ages <3.6 d after maximum brightness. These data were obtained as part of the Berkeley SN Ia Program (BSNIP; Silverman et al. 2012) and represent the largest set of SNe Ia in which C has ever been searched. We find that ~11 per cent of the SNe studied show definite C absorption features while ~25 per cent show some evidence for C II in their spectra. Also, if one obtains a spectrum at t < -5 d, then there is a better than 30 per cent chance of detecting a distinct absorption feature from C II. SNe Ia that show C are found to resemble those without C in many respects, but objects with C tend to have bluer optical colours than those without C. The typical expansion velocity of the C II {\lambda}6580 feature is measured to be 12,000-13,000 km/s, and the ratio of the C II {\lambda}6580 to Si II {\lambda}6355 velocities is remarkably constant with time and among different objects with a median value of ~1.05. While the pseudo-equivalent widths (pEWs) of the C II {\lambda}6580 and C II {\lambda}7234 features are found mostly to decrease with time, we see evidence of a significant increase in pEW between ~12 and 11 d before maximum brightness, which is actually predicted by some theoretical models. The range of pEWs measured from the BSNIP data implies a range of C mass in SN Ia ejecta of about (2-30) * 10^-3 M_Sun.


title: A comment on anti-brane singularities in warped throats

abstract: We compute the imaginary self-dual (ISD) and imaginary anti-self-dual (IASD) fluxes for the Klebanov-Strassler background perturbed by a stack of p anti-D3 branes. We show that, at linear order in p, they both have a singularity in the near-brane region. While one can argue that the IASD flux may disappear at full non-linear level, no such argument exists for the ISD mode. An analogy with anti-D6 backreaction suggests that such singularity may survive once full backreaction is taken into account and may be a universal feature of anti-brane solutions.


title: B2 0954+25A: a typical Fermi blazar or a gamma-loud Narrow Line Seyfert 1

abstract: B2 0954+25A, detected by the {\it Fermi} satellite, is a blazar with interesting observational properties: it has been observed to transit from a jet dominated to a disk dominated state; its radio spectrum appears flat at all observing frequencies (down to 74 MHz); optically, the H$\beta$ line profile is asymmetric. The flatness of radio spectrum suggests that the isotropic emission from radio lobes is very weak, despite the large size of its jet ($\gtrsim$ 500 kpc). Its broad--band spectral energy distribution is surprisingly similar to that of the prototypical $\gamma$--ray, radio loud, Narrow Line Seyfert 1 ($\gamma$--NLS1) galaxy PMN J0948+0022. In this work we revisit the mass estimates of B2 0954+25A considering only the symmetric component of the H$\beta$ line and find (1--3) $\times 10^8$ M$_{\sun}$. In light of our composite analysis, we propose to classify the source as a transition object between the class of Flat Spectrum Radio Quasar and $\gamma$--ray, radio loud NLS1. A comparison with two members of each class (3C 273 and PMN J0948+0022) is discussed.


title: On the average density profile of dark-matter halos in the inner regions of massive early-type galaxies

abstract: We study a sample of 39 massive early-type lens galaxies at redshift z < 0.3 to determine the slope of the average dark-matter density profile in the innermost regions. We keep the strong lensing and stellar population synthesis modeling as simple as possible to measure the galaxy total and luminous masses. By rescaling the values of the Einstein radius and dark-matter projected mass with the values of the luminous effective radius and mass, we combine all the data of the galaxies in the sample. We find that between 0.3 and 0.9 times the value of the effective radius the average logarithmic slope of the dark-matter projected density profile is -1.0 +/- 0.2 (i.e., approximately isothermal) or -0.7 +/- 0.5 (i.e., shallower than isothermal), if, respectively, a constant Chabrier or heavier, Salpeter-like stellar IMF is adopted. These results provide positive evidence of the influence of the baryonic component on the contraction of the galaxy dark-matter halos, compared to the predictions of dark matter-only cosmological simulations, and open a new way to test models of structure formation and evolution within the standard LCDM cosmological scenario.


title: Dissipativity of the delay semigroup

abstract: Under mild conditions a delay semigroup can be transformed into a (generalized) contraction semigroup by modifying the inner product on the (Hilbert) state space into an equivalent inner product. Applications to stability of differential equations with delay and stochastic differential equations with delay are given as examples.


title: Critical Temperature Associated to Symmetry Breaking of Klein--Gordon fields versus Condensation Temperature in a Weakly interacting Bose--Einstein Gas

abstract: We deduce the relation between the critical temperature associated to the symmetry breaking of scalar fields with one--loop correction potential immersed in a thermal bath and the condensation temperature of the aforementioned system, assuming a harmonic oscillator type potential. We show that these two temperatures are related through the \emph{scale} associated to the system. In this aim, we infer the order of magnitude for the \emph{scale} as a function of the corresponding healing length, in order to give a criterium to compare both temperatures. Additionally, we prove that the condensation temperature is independent of the thermal bath within the semiclassical approximation, for a positive coupling constant, assuming that the thermal bath contribution is the lowest energy associated to the system.


title: A QCD Analysis of Double Parton Scattering: Color Correlations, Interference Effects and Evolution

abstract: We derive a factorization formula for the double Drell-Yan cross section in terms of double parton distribution functions (dPDFs). Diparton flavor, spin and color correlations and parton-exchange interference terms contribute, even for unpolarized beams. Soft radiation effects are nontrivial for the color correlation and interference contributions, and are described by non-perturbative soft functions. We provide a field-theoretic definition of the quark dPDFs and study some of their basic properties, including discrete symmetries and their interpretation in a non-relativistic quark model. We calculate the renormalization group evolution of the quark dPDFs and of the soft functions. The evolution receives contributions from both ultraviolet and rapidity divergences. We find that color correlation and interference effects are Sudakov suppressed, greatly reducing the number of dPDFs needed to describe double parton scattering at high energy experiments.


title: Repercussions of Flavour Symmetry Breaking on CP Violation in D-Meson Decays

abstract: We investigate to what extent the recently measured value for a non-vanishing direct CP asymmetry in D0 -> K+ K- and D0 -> pi+ pi- decays can be accommodated in the Standard Model (SM) or extensions with a constrained flavour sector, for instance from a sequential 4th generation of quarks (4G). From the comparison with D0 -> K- pi+ branching ratios, we establish large U-spin symmetry (d <-> s) breaking effects with large strong phases between different interfering amplitudes. On the basis of conservative estimates on amplitude ratios - which are supported by an analysis of the breaking of a (c <-> u) symmetry in non-leptonic B0 decays - we find that, in the SM, direct CP asymmetries in the pi+ pi- or K+ K- modes (or in their difference) of the order of several per mille are still plausible. Due to the constraints on the new CP phases in the 4G model, only moderate effects compared to the SM estimates are possible. We suggest CP studies at LHCb as well as at (Super)B-factories of several distinctive modes, such as D+ -> Kbar(*)0 pi+, phi pi+ and Ds -> K(*)0 pi+, phi pi+ (K+) etc., which should shed more light on the short- and long-distance issues underlying CP violation in non-leptonic D-meson decays.


title: Understanding surface-adsorption of proteins: the Vroman effect

abstract: It is now well accepted that cellular responses to materials in a biological medium reflect greatly the adsorbed biomolecular layer, rather than the material itself. Here, we study by molecular dynamic simulations the competitive protein adsorption on a surface (Vroman effect), i.e. the non-monotonic behavior of the amount of protein adsorbed on a surface in contact with plasma as a function of contact time and plasma concentration. We find a complex behavior, with regimes during which small and large proteins are not necessarily competing between them, but are both competing with others in solution. We show how the effect can be understood, controlled and inverted.


title: Nonlinear, electrocatalytic swimming in the presence of salt

abstract: A small, bimetallic particle in a hydrogen peroxide solution can propel itself by means of an electrocatalytic reaction. The swimming is driven by a flux of ions around the particle. We model this process for the presence of a monovalent salt, where reaction-driven proton currents induce salt ion currents. A theory for thin diffuse layers is employed, which yields nonlinear, coupled transport equations. The boundary conditions include a compact Stern layer of adsorbed ions. Electrochemical processes on the particle surface are modeled with a first order reaction of the Butler-Volmer type. The equations are solved numerically for the swimming speed. An analytical approximation is derived under the assumption that the decomposition of hydrogen peroxide occurs mainly without inducing an electric current. We find that the swimming speed increases linearly with hydrogen peroxide concentration for small concentrations. The influence of ion diffusion on the reaction rate can lead to a concave shape of the function of speed vs. hydrogen peroxide concentration. The compact layer of ions on the particle diminishes the reaction rate and consequently reduces the speed. Our results are consistent with published experimental data.

categories: cond-mat.soft physics.bio-ph physics.chem-ph physics.flu-dyn


title: The BF theory as an electric Julia-Toulouse condensate

abstract: The Julia-Toulouse mechanism is used to show that the topological Abelian BF term may be induced by the condensation of electric charges. As an application we discuss the subtle question of including consistently magnetic defects into the Maxwell-BF theory in a way to avoid the usual problems of current conservation, charge quantization, Elitzur's theorem violation and the reality of the Dirac brane, produced by the non-minimal coupling. We also discuss a new way of obtaining the Dirac's veto, which is based on the so-called Dirac brane symmetry.


title: Studies of material properties under irradiation at BNL Linear Isotope Producer (BLIP)

abstract: Effects of proton beams irradiating materials considered for targets in high-power accelerator experiments have been under study using the Brookhaven National Laboratory's (BNL) 200 MeV Linac. The primary objectives of the study that includes a wide array of materials and alloys ranging between low and high-Z are to (a) observe changes in physio-mechanical properties which are important in maintaining high-power target functionality, (b) identify possible limits of proton flux or fluence above which certain material seize to maintain integrity, (c) study the role of material operating temperature in inducing or maintaining radiation damage reversal, and (d) correlate radiation damage effects of different species such as energetic protons and neutrons on materials by utilizing reactor and particle accelerator experience data. These objectives are specifically being addressed in the latest material irradiation study linked to the Long Baseline Neutrino Experiment (LBNE). Observations on irradiation effects on materials considered for high-power targets and other beam intercepting elements, such as collimators, from past studies and preliminary observations of the ongoing LBNE study are presented in this paper.


title: On the net reproduction rate of continuous structured populations with distributed states at birth

abstract: We consider a nonlinear structured population model with a distributed recruitment term. The question of the existence of non-trivial steady states can be treated (at least!) in three different ways. One approach is to study spectral properties of a parametrized family of unbounded operators. The alternative approach, on which we focus here, is based on the reformulation of the problem as an integral equation. In this context we introduce a density dependent net reproduction rate and discuss its relationship to a biologically meaningful quantity. Finally, we briefly discuss a third approach, which is based on the finite rank approximation of the recruitment operator.


title: Ideal-Modified Bosonic Gas Trapped in an Arbitrary Three Dimensional Power-Law Potential

abstract: We analyze the effects caused by an anomalous single-particle dispersion relation suggested in several quantum-gravity models, upon the thermodynamics of a Bose-Einstein condensate trapped in a generic 3-dimensional power-law potential. We prove that the shift in the condensation temperature, caused by a deformed dispersion relation, described as a non-trivial function of the number of particles and the shape associated to the corresponding trap, could provide bounds for the parameters associated to such deformation. Additionally, we calculate the fluctuations in the number of particles as a criterium of thermodynamic stability for these systems. We show that the apparent instability caused by the anomalous fluctuations in the thermodynamic limit can be suppressed considering the lowest energy associated to the system in question.


title: Sea Ice Brightness Temperature as a Function of Ice Thickness: Computed curves for AMSR-E and SMOS (frequencies from 1.4 to 89 GHz)

abstract: The relationship between sea ice thickness and microwave brightness temperature is explored. Parameterized ice profiles are fed to a radiative-transfer-based sea ice emissivity model (Microwave Emmission of Layered Snowpack, MEMLS). Complex permittivities, required as input for the simulation, are determined using a semi-empirical mixture model. Since the thickness-radiance relation is not fixed but can vary significantly depending upon past and current weather conditions, we determine a range of brightness temperature values for each thickness. This is done using a bootstrap model in which the salinities are varied based on variances supplied with the thickness-salinity curve and the complex permittivities are varied based on variance supplied by the mixture model. The results suggest that scattering is one of the most important parameters determining sea ice brightness temperature, especially for new and forming ice. Further work must be done to accurately model both scattering and complex permittivities in sea ice.


title: Muon Collider

abstract: Both e+e- and {\mu}+{\mu}- colliders have been proposed as possible candidates for a lepton collider to complement and extend the reach of the Large Hadron Collider (LHC) at CERN. The physics program that could be pursued by a new lepton collider (e+e- or {\mu}+{\mu}-) with sufficient luminosity would include understanding the mechanism behind mass generation and electroweak symmetry breaking; searching for, and possibly discovering, supersymmetric particles; and hunting for signs of extra spacetime dimensions and quantum gravity. However, the appropriate energy reach for such a collider is currently unknown, and will only be determined following initial physics results at the LHC. It is entirely possible that such results will indicate that a lepton collider with a collision energy well in excess of 1 TeV will be required to illuminate the physics uncovered at LHC. Such a requirement would require consideration of muons as the lepton of choice for such a collider.


title: A discrete model for the apparent viscosity of polydisperse suspensions including maximum packing fraction

abstract: Based on the notion of a construction process consisting of the stepwise addition of particles to the pure fluid, a discrete model for the apparent viscosity as well as for the maximum packing fraction of polydisperse suspensions of spherical, non-colloidal particles is derived. The model connects the approaches by Bruggeman and Farris and is valid for large size ratios of consecutive particle classes during the construction process, appearing to be the first model consistently describing polydisperse volume fractions and maximum packing fraction within a single approach. In that context, the consistent inclusion of the maximum packing fraction into effective medium models is discussed. Furthermore, new generalized forms of the well-known Quemada and Krieger equations allowing for the choice of a second-order Taylor coefficient for the volume fraction ($\phi^2$-coefficient), found by asymptotic matching, are proposed. The model for the maximum packing fraction as well as the complete viscosity model are compared to experimental data from the literature showing good agreement. As a result, the new model is shown to replace the empirical Sudduth model for large diameter ratios. The extension of the model to the case of small size ratios is left for future work.


title: On origin of 1/f noise in manganites: memoryless transport against mysterious slow fluctuators

abstract: An alternative explanation of 1/f-noise in manganites is suggested and discussed


title: Modified Bosonic Gas Trapped in a Generic 3-dim Power Law Potential

abstract: We analyze the consequences caused by an anomalous single-particle dispersion relation suggested in several quantum-gravity models, upon the thermodynamics of a Bose-Einstein condensate trapped in a generic 3-dimensional power-law potential. We prove that the condensation temperature is shifted as a consequence of such deformation and show that this fact could be used to provide bounds on the deformation parameters. Additionally, we show that the shift in the condensation temperature, described as a non-trivial function of the number of particles and the trap parameters, could be used as a criterion to analyze the effects caused by a deformed dispersion relation in weakly interacting systems and also in finite size systems.


title: An Adaptive Mechanism for Accurate Query Answering under Differential Privacy

abstract: We propose a novel mechanism for answering sets of count- ing queries under differential privacy. Given a workload of counting queries, the mechanism automatically selects a different set of "strategy" queries to answer privately, using those answers to derive answers to the workload. The main algorithm proposed in this paper approximates the optimal strategy for any workload of linear counting queries. With no cost to the privacy guarantee, the mechanism improves significantly on prior approaches and achieves near-optimal error for many workloads, when applied under (\epsilon, \delta)-differential privacy. The result is an adaptive mechanism which can help users achieve good utility without requiring that they reason carefully about the best formulation of their task.


title: Theorematum quorundam arithmeticorum demonstrationes

abstract: Euler proves that the sum of two 4th powers can't be a 4th power and that the difference of two distinct non-zero 4th powers can't be a 4th power and Fermat's theorem that the equation x(x+1)/2=y^4 can only be solved in integers if x=1 and the final theorem y^3+1=x^2 can only be solves for x=3 and y=2 in integers. The paper is translated from Euler's Latin original into German.


title: IDR muon capture front end and variations

abstract: The (International Design Report) IDR neutrino factory scenario for capture, bunching, phase-energy rotation and initial cooling of micro's produced from a proton source target is explored. It requires a drift section from the target, a bunching section and a -E rotation section leading into the cooling channel. The rf frequency changes along the bunching and rotation transport in order to form the 's into a train of equal-energy bunches suitable for cooling and acceleration. Optimization and variations are discussed. An important concern is rf limitations within the focusing magnetic fields, mitigation procedures are described. The method can be extended to provide muons for a micro+-micro < Collider, variations toward optimizing that extension are discussed.


title: Disorder promotes ferromagnetism: Rounding of the quantum phase transition in Sr_{1-x}Ca_xRuO_3

abstract: The subtle interplay of randomness and quantum fluctuations at low temperatures gives rise to a plethora of unconventional phenomena in systems ranging from quantum magnets and correlated electron materials to ultracold atomic gases. Particularly strong disorder effects have been predicted to occur at zero-temperature quantum phase transitions. Here, we demonstrate that the composition-driven ferromagnetic-to-paramagnetic quantum phase transition in Sr1-xCaxRuO3 is completely destroyed by the disorder introduced via the different ionic radii of the randomly distributed Sr and Ca ions. Using a magneto-optical technique, we map the magnetic phase diagram in the composition-temperature space. We find that the ferromagnetic phase is significantly extended by the disorder and develops a pronounced tail over a broad range of the composition x. These findings are explained by a microscopic model of smeared quantum phase transitions in itinerant magnets. Moreover, our theoretical study implies that correlated disorder is even more powerful in promoting ferromagnetism than random disorder.


title: Summary of session 3 on synchrotron radiation and beam dynamics

abstract: Below we summarize presentations, discussions and general conclusions of the Workshop session on "Beam Dynamics Issues". Major subjects include effects due to synchrotron radiation (SR), cryogenic loads, electron cloud, impedances, intra-beam scattering (IBS) and beam-beam interactions.


title: Flavour Condensate and the Dark Sector of the Universe

abstract: This thesis is devoted to the development of a nonperturbative quantum field theoretical approach to flavour physics, with special attention to cosmological applications. Neutrino flavour oscillation is nowadays a fairly well-established experimental fact. However, the formulation of flavour oscillations in a relativistic field theoretical framework presents non-trivial difficulties. A nonperturbative approach for building flavour states has been proposed by Blasone, Vitiello and coworkers. The formalism implies a non-trivial physical vacuum (called "flavour vacuum"), which might act as a source of Dark Energy. Furthermore, such a vacuum has been recognized as the effective vacuum state arising in the low energy limit of a string theoretical model, D-particle Foam Model. In the attempt of probing the observable phenomenology of the D-particle foam model, a simple toy model (two scalars with mixing \`a la Blasone & Vitiello on a adiabatically expanding background) has been studied, proving that the flavour vacuum might behave as Dark Energy under certain assumptions. The first work presented in this thesis represents a development of this approach. A more realistic model is considered, which includes two flavoured Dirac fermions on a generic Friedmann-Robertson-Walker universe. In this framework we show that the flavour vacuum presents different features, which are incompatible with Dark Energy. Motivated by this discrepancy, we next embark on the analysis of a simple supersymmetric model in flat spacetime (free Wess-Zumino), proving that the bosonic component of flavour vacuum acts as Dark Energy, whereas the fermionic as a source of Dark Matter. Finally we develop a new method of calculation that open the way to a nonperturbative extension of these results for interactive theories.


title: Solitons supported by singular spatial modulation of the Kerr nonlinearity

abstract: We introduce a setting based on the one-dimensional (1D) nonlinear Schroedinger equation (NLSE) with the self-focusing (SF) cubic term modulated by a singular function of the coordinate, |x|^{-a}. It may be additionally combined with the uniform self-defocusing (SDF) nonlinear background, and with a similar singular repulsive linear potential. The setting, which can be implemented in optics and BEC, aims to extend the general analysis of the existence and stability of solitons in NLSEs. Results for fundamental solitons are obtained analytically and verified numerically. The solitons feature a quasi-cuspon shape, with the second derivative diverging at the center, and are stable in the entire existence range, which is 0 < a < 1. Dipole (odd) solitons are found too. They are unstable in the infinite domain, but stable in the semi-infinite one. In the presence of the SDF background, there are two subfamilies of fundamental solitons, one stable and one unstable, which exist together above a threshold value of the norm (total power of the soliton). The system which additionally includes the singular repulsive linear potential emulates solitons in a uniform space of the fractional dimension, 0 < D < 1. A two-dimensional extension of the system, based on the quadratic nonlinearity, is formulated too.


title: High Frame-rate Imaging Based Photometry, Photometric Reduction of Data from Electron-multiplying Charge Coupled Devices (EMCCDs)

abstract: The EMCCD is a type of CCD that delivers fast readout times and negligible readout noise, making it an ideal detector for high frame rate applications which improve resolution, like lucky imaging or shift-and-add. This improvement in resolution can potentially improve the photometry of faint stars in extremely crowded fields significantly by alleviating crowding. Alleviating crowding is a prerequisite for observing gravitational microlensing in main sequence stars towards the galactic bulge. However, the photometric stability of this device has not been assessed. The EMCCD has sources of noise not found in conventional CCDs, and new methods for handling these must be developed. We aim to investigate how the normal photometric reduction steps from conventional CCDs should be adjusted to be applicable to EMCCD data. One complication is that a bias frame cannot be obtained conventionally, as the output from an EMCCD is not normally distributed. Also, the readout process generates spurious charges in any CCD, but in EMCCD data, these charges are visible as opposed to the conventional CCD. Furthermore we aim to eliminate the photon waste associated with lucky imaging by combining this method with shift-and-add. A simple probabilistic model for the dark output of an EMCCD is developed. Fitting this model with the expectation-maximization algorithm allows us to estimate the bias, readout noise, amplification, and spurious charge rate per pixel and thus correct for these phenomena. To investigate the stability of the photometry, corrected frames of a crowded field are reduced with a PSF fitting photometry package, where a lucky image is used as a reference. We find that it is possible to develop an algorithm that elegantly reduces EMCCD data and produces stable photometry at the 1% level in an extremely crowded field.


title: A Relic Star Cluster in the Sextans Dwarf Spheroidal Galaxy - Implications for Early Star and Galaxy Formation

abstract: We present tentative evidence for the existence of a dissolved star cluster in the Sextans dwarf spheroidal galaxy. In a sample of six stars, we identify three (possibly four) stars around [Fe/H] =-2.7 that are highly clustered in a multi-dimensional chemical abundance space. The estimated initial stellar mass of the cluster is M*,init = 1.9^+1.5_-0.9 (1.6^+1.2_-0.8)*10^5 Msol assuming a Salpeter (Kroupa) initial mass function (IMF). If corroborated by follow-up spectroscopy, this ancient star cluster at [Fe/H] =-2.7 is the most metal-poor system identified to date. Inspired by this finding, we also present a new way to interpret the cumulative metallicity functions of dwarf galaxies. From available observational data, we speculate that the ultra-faint dwarf galaxy population, or a significant fraction thereof, and the more luminous, classical dwarf spheroidal population were formed in different environments and would thus be distinct in origin.


title: Shielding experiments by the JASMIN collaboration at Fermilab (II) - Radioactivity measurement induced by secondary particles from the anti-proton production target

abstract: The JASMIN Collaboration has performed an experiment to conduct measurements of nuclear reaction rates around the anti-proton production (Pbar) target at the Fermi National Accelerator Laboratory (FNAL). At the Pbar target station, the target, consisting an Inconel 600 cylinder, was irradiated by a 120 GeV/c proton beam from the FNAL Main Injector. The beam intensity was 3.6 x 10**12 protons per second. Samples of Al, Nb, Cu, and Au were placed near the target to investigate the spatial and energy distribution of secondary particles emitted from it. After irradiation, the induced activities of the samples were measured by studying their gamma ray spectra using HPGe detectors. The production rates of 30 nuclides induced in Al, Nb, Cu, Au samples were obtained. These rates increase for samples placed in a forward (small angle) position relative to the target. The angular dependence of these reaction rates becomes larger for increasing threshold energy. These experimental results are compared with Monte Carlo calculations. The calculated results generally agree with the experimental results to within a factor of 2 to 3.


title: On infinite-dimensional state spaces

abstract: It is well-known that the canonical commutation relation $[x,p]=i$ can be realized only on an infinite-dimensional Hilbert space. While any finite set of experimental data can also be explained in terms of a finite-dimensional Hilbert space by approximating the commutation relation, Occam's razor prefers the infinite-dimensional model in which $[x,p]=i$ holds on the nose. This reasoning one will necessarily have to make in any approach which tries to detect the infinite-dimensionality. One drawback of using the canonical commutation relation for this purpose is that it has unclear operational meaning. Here, we identify an operationally well-defined context from which an analogous conclusion can be drawn: if two unitary transformations $U,V$ on a quantum system satisfy the relation $V^{-1}U^2V=U^3$, then finite-dimensionality entails the relation $UV^{-1}UV=V^{-1}UVU$; this implication strongly fails in some infinite-dimensional realizations. This is a result from combinatorial group theory for which we give a new proof. This proof adapts to the consideration of cases where the assumed relation $V^{-1}U^2V=U^3$ holds only up to $\eps$ and then yields a lower bound on the dimension.


title: Source-like solution for radial imbibition into a homogeneous semi-infinite porous medium

abstract: We describe the imbibition process from a point source into a homogeneous semi-infinite porous material. When body forces are negligible, the advance of the wetting front is driven by capillary pressure and resisted by viscous forces. With the assumption that the wetting front assumes a hemispherical shape, our analytical results show that the absorbed volume flow rate is approximately constant with respect to time, and that the radius of the wetting evolves in time as r \approx t^(1/3). This cube-root law for the long-time dynamics is confirmed by experiments using a packed cell of glass microspheres with average diameter of 42 {\mu}m. This result complements the classical one-dimensional imbibition result where the imbibition length l \approx t^(1/2), and studies in axisymmetric porous cones with small opening angles where l \approx t^(1/4) at long times.


title: A Comparative Review of Dimension Reduction Methods in Approximate Bayesian Computation

abstract: Approximate Bayesian computation (ABC) methods make use of comparisons between simulated and observed summary statistics to overcome the problem of computationally intractable likelihood functions. As the practical implementation of ABC requires computations based on vectors of summary statistics, rather than full data sets, a central question is how to derive low-dimensional summary statistics from the observed data with minimal loss of information. In this article we provide a comprehensive review and comparison of the performance of the principal methods of dimension reduction proposed in the ABC literature. The methods are split into three nonmutually exclusive classes consisting of best subset selection methods, projection techniques and regularization. In addition, we introduce two new methods of dimension reduction. The first is a best subset selection method based on Akaike and Bayesian information criteria, and the second uses ridge regression as a regularization procedure. We illustrate the performance of these dimension reduction techniques through the analysis of three challenging models and data sets.


title: Study of degenerate parabolic system modeling the hydrogen displacement in a nuclear waste repository

abstract: Our goal is the mathematical analysis of a two phase (liquid and gas) two components (water and hydrogen) system modeling the hydrogen displacement in a storage site for radioactive waste. We suppose that the water is only in the liquid phase and is incompressible. The hydrogen in the gas phase is supposed compressible and could be dissolved into the water with the Henry's law. The flow is described by the conservation of the mass of each components. The model is treated without simplified assumptions on the gas density. This model is degenerated due to vanishing terms. We establish an existence result for the nonlinear degenerate parabolic system based on new energy estimate on pressures.


title: A Hierachical Evolutionary Algorithm for Multiobjective Optimization in IMRT

abstract: Purpose: Current inverse planning methods for IMRT are limited because they are not designed to explore the trade-offs between the competing objectives between the tumor and normal tissues. Our goal was to develop an efficient multiobjective optimization algorithm that was flexible enough to handle any form of objective function and that resulted in a set of Pareto optimal plans. Methods: We developed a hierarchical evolutionary multiobjective algorithm designed to quickly generate a diverse Pareto optimal set of IMRT plans that meet all clinical constraints and reflect the trade-offs in the plans. The top level of the hierarchical algorithm is a multiobjective evolutionary algorithm (MOEA). The genes of the individuals generated in the MOEA are the parameters that define the penalty function minimized during an accelerated deterministic IMRT optimization that represents the bottom level of the hierarchy. The MOEA incorporates clinical criteria to restrict the search space through protocol objectives and then uses Pareto optimality among the fitness objectives to select individuals. Results: Acceleration techniques implemented on both levels of the hierarchical algorithm resulted in short, practical runtimes for optimizations. The MOEA improvements were evaluated for example prostate cases with one target and two OARs. The modified MOEA dominated 11.3% of plans using a standard genetic algorithm package. By implementing domination advantage and protocol objectives, small diverse populations of clinically acceptable plans that were only dominated 0.2% by the Pareto front could be generated in a fraction of an hour. Conclusions: Our MOEA produces a diverse Pareto optimal set of plans that meet all dosimetric protocol criteria in a feasible amount of time. It optimizes not only beamlet intensities but also objective function parameters on a patient-specific basis.


title: No-signaling Quantum Key Distribution: Solution by Linear Programming

abstract: We outline a straightforward approach for obtaining a secret key rate using only no-signaling constraints and linear programming. Assuming an individual attack, we consider all possible joint probabilities. Initially, we study only the case where Eve has binary outcomes, and we impose constraints due to the no-signaling principle and given measurement outcomes. Within the remaining space of joint probabilities, by using linear programming, we get bound on the probability of Eve correctly guessing Bob's bit. We then make use of an inequality that relates this guessing probability to the mutual information between Bob and a more general Eve, who is not binary-restricted. Putting our computed bound together with the Csisz\'ar-K\"orner formula, we obtain a positive key generation rate. The optimal value of this rate agrees with known results, but was calculated in a more straightforward way, offering the potential of generalization to different scenarios.


title: When is Better Best? A multiobjective perspective

abstract: Purpose: To identify the most informative methods for reporting results of treatment planning comparisons. Methods: Seven papers from the past year of International Journal of Radiation Oncology Biology Physics reported on comparisons of treatment plans for IMRT and IMAT. The papers were reviewed to identify methods of comparisons. Decision theoretical concepts were used to evaluate the study methods and highlight those that provide the most information. Results: None of the studies examined the correlation between objectives. Statistical comparisons provided some information but not enough to make provide support for a robust decision analysis. Conclusion: The increased use of treatment planning studies to evaluate different methods in radiation therapy requires improved standards for designing the studies and reporting the results.


title: Physical Layer Security for Two-Way Untrusted Relaying with Friendly Jammers

abstract: In this paper, we consider a two-way relay network where two sources can communicate only through an untrusted intermediate relay, and investigate the physical layer security issue of this two-way relay scenario. Specifically, we treat the intermediate relay as an eavesdropper from which the information transmitted by the sources needs to be kept secret, despite the fact that its cooperation in relaying this information is essential. We indicate that a non-zero secrecy rate is indeed achievable in this two-way relay network even without external friendly jammers. As for the system with friendly jammers, after further analysis, we can obtain that the secrecy rate of the sources can be effectively improved by utilizing proper jamming power from the friendly jammers. Then, we formulate a Stackelberg game model between the sources and the friendly jammers as a power control scheme to achieve the optimized secrecy rate of the sources, in which the sources are treated as the sole buyer and the friendly jammers are the sellers. In addition, the optimal solutions of the jamming power and the asking prices are given and a distributed updating algorithm to obtain the Stakelberg equilibrium is provided for the proposed game. Finally, the simulations results verify the properties and the efficiency of the proposed Stackelberg game based scheme.


title: Bonding trends within ternary Isocoordinate chalcogenide glasses GeAsSe

abstract: A structural study is presented of ab-initio molecular dynamics simulations of Ge-As-Se calcogenide glasses performed at the same mean coordination number but differing stoichiometry ranging between Se rich and Se poor glasses. Starting configurations are generated via Reverse Monte Carlo (RMC) simulations of Extended X-ray Absorption Fine Structure (EXAFS) measurements of experimental samples. Structural analysis is presented illustrating the bonding trends found with changing stoichiometry.


title: A new approach to asymptotics of Maclaurin coefficients of algebraic functions

abstract: We propose a general method for deriving asymptotics of the Maclaurin series coefficients of algebraic functions that is based on a procedure of K. V. Safonov and multivariate singularity analysis. We test the feasibility of this this approach by experimenting on several examples.


title: Search for large extra dimensions in dimuon and dielectron events in pp collisions at sqrt(s) = 7 TeV

abstract: Results are presented from a search for large, extra spatial dimensions in events with either two isolated muons or two isolated electrons. The data are from proton-proton interactions at sqrt(s) = 7 TeV collected with the CMS detector at the LHC. The size of the data sample corresponds to an integrated luminosity of approximately 2 inverse femtobarns. The observed dimuon and dielectron mass spectra are found to be consistent with standard-model expectations. Depending on the number of extra dimensions, the 95% confidence level limits from the combined dimuon and dielectron channels range from Ms > 2.4 TeV to Ms > 3.8 TeV, where Ms characterizes the scale for the onset of quantum gravity.


title: Electron self-injection for the acceleration in laser-pulse-wakes in the presence of a `strong' external magnetic field

abstract: An external static magnetic field with its strength B~10T may result in the laser wake wave-breaking upon changing the electron motion in the vicinity of maximal density ramp of a wave period. This, as shown by numerical simulations, can change the resonance character of the electron self-injection in the laser wake-field; a total charge loaded in the acceleration phase of laser pulse wake can be controlled by a proper choice of the magnetic field strength.


title: Towards consistent mapping of distant worlds: secondary-eclipse scanning of the exoplanet HD189733b

abstract: Mapping distant worlds is the next frontier for exoplanet infrared photometry studies. Ultimately, constraining spatial and temporal properties of an exoplanet atmosphere will provide further insight into its physics. For tidally-locked hot Jupiters that transit and are eclipsed by their host star, the first steps are now possible. Our aim is to constrain an exoplanet's shape, brightness distribution (BD) and system parameters from its light curve. Notably, we rely on the eclipse scanning. We use archived Spitzer 8-{\mu}m data of HD189733 (6 transits, 8 secondary eclipses, and a phase curve) in a global MCMC procedure for mitigating systematics. We also include HD189733's out-of-transit radial velocity measurements. We find a 6-{\sigma} deviation from the expected occultation of a uniformly-bright disk. This deviation emerges mainly from HD189733b's thermal pattern, not from its shape. We indicate that the correlation of the orbital eccentricity, e, and BD (uniform time offset) does also depend on the stellar density, \rho*, and the impact parameter, b (e-b-\rho*-BD correlation). For HD189733b, we find that relaxing the e-constraint and using more complex BDs lead to lower stellar/planetary densities and a more localized and latitudinally-shifted hot spot. We obtain an improved constraint on the upper limit of HD189733b's orbital eccentricity, e<0.011 (95%), when including the RV measurements. Our study provides new insights into the analysis of exoplanet light curves and a proper framework for future eclipse-scanning observations. Observations of the same exoplanet at different wavelengths will improve the constraints on its system parameters while ultimately yielding a large-scale time-dependent 3D map of its atmosphere. Finally, we discuss the perspective of extending our method to observations in the visible, in particular to better understand exoplanet albedos.


title: Four Fundamental Foreground Power Spectrum Shapes for 21 cm Cosmology Observations

abstract: Contamination from instrumental effects interacting with bright astrophysical sources is the primary impediment to measuring Epoch of Reionization and BAO 21 cm power spectra---an effect called mode-mixing. In this paper we identify four fundamental power spectrum shapes produced by mode-mixing that will affect all upcoming observations. We are able, for the first time, to explain the wedge-like structure seen in advanced simulations and to forecast the shape of an 'EoR window' that is mostly free of contamination. Understanding the origins of these contaminations also enables us to identify calibration and foreground subtraction errors below the imaging limit, providing a powerful new tool for precision observations.


title: Quantum Enhanced Magnetometer with Low-Frequency Squeezing

abstract: We report the demonstration of a magnetometer with noise-floor reduction below the shot-noise level. This magnetometer, based on a nonlinear magneto-optical rotation effect, is enhanced by the injection of a squeezed vacuum state into its input. The noise spectrum shows squeezed noise reduction of about 2 dB spanning from close to 100 Hz to several megahertz. We also report on the observation of two different regimes of operation of such a magnetometer: one in which the detection noise is limited by the quantum noise of the light probe only, and one in which we see additional noise originating from laser noise which is rotated into the vacuum polarization.


title: Cosmicflows-2: SNIa Calibration and H0

abstract: The construction of the Cosmicflows-2 compendium of distances involves the merging of distance measures contributed by the following methods: (Cepheid) Period-Luminosity, Tip of the Red Giant Branch (TRGB), Surface Brightness Fluctuation (SBF), Luminosity-Linewidth (TF), Fundamental Plane (FP), and Type Ia supernova (SNIa). The method involving SNIa is at the top of an interconnected ladder, providing accurate distances to well beyond the expected range of distortions to Hubble flow from peculiar motions. In this paper, the SNIa scale is anchored by 36 TF spirals with Cepheid or TRGB distances, 56 SNIa hosts with TF distances, and 61 groups or clusters hosting SNIa with Cepheid, SBF, TF, or FP distances. With the SNIa scale zero point set, a value of the Hubble Constant is evaluated over a range of redshifts 0.03 < z < 0.5, assuming a cosmological model with Omega_m = 0.27 and Omega_Lambda = 0.73. The value determined for the Hubble Constant is H0 = 75.9 \pm 3.8 km s-1 Mpc-1.


title: Pseudogap in underdoped Ba$_{1-x}$K$_{x}$Fe$_{2}$As$_{2}$ as seen via optical conductivity

abstract: We report the observation of a pseudogap in the \emph{ab}-plane optical conductivity of underdoped Ba$_{1-x}$K$_{x}$Fe$_{2}$As$_{2}$ ($x = 0.2$ and 0.12) single crystals. Both samples show prominent gaps opened by a spin density wave (SDW) order and superconductivity at the transition temperatures $T_{\it SDW}$ and $T_c$, respectively. In addition, we observe an evident pseudogap below $T^{\ast} \sim$ 75 K, a temperature much lower than $T_{\it SDW}$ but much higher than $T_{c}$. A spectral weight analysis shows that the pseudogap is closely connected to the superconducting gap, indicating the possibility of its being a precursor of superconductivity. The doping dependence of the gaps is also supportive of such a scenario.


title: Developing and applying heterogeneous phylogenetic models with XRate

abstract: Modeling sequence evolution on phylogenetic trees is a useful technique in computational biology. Especially powerful are models which take account of the heterogeneous nature of sequence evolution according to the "grammar" of the encoded gene features. However, beyond a modest level of model complexity, manual coding of models becomes prohibitively labor-intensive. We demonstrate, via a set of case studies, the new built-in model-prototyping capabilities of XRate (macros and Scheme extensions). These features allow rapid implementation of phylogenetic models which would have previously been far more labor-intensive. XRate's new capabilities for lineage-specific models, ancestral sequence reconstruction, and improved annotation output are also discussed. XRate's flexible model-specification capabilities and computational efficiency make it well-suited to developing and prototyping phylogenetic grammar models. XRate is available as part of the DART software package: http://biowiki.org/DART .


title: Effective embedding of residually hyperbolic groups into direct products of extensions of centralizers

abstract: For any torsion-free hyperbolic group $\Gamma$ and any group $G$ that is fully residually $\Gamma$, we construct algorithmically a finite collection of homomorphisms from $G$ to groups obtained from $\Gamma$ by extensions of centralizers, at least one of which is injective. When $G$ is residually $\Gamma$, this gives a effective embedding of $G$ into a direct product of such groups. We also give an algorithmic construction of a diagram encoding the set of homomorphisms from a given finitely presented group to $\Gamma$.


title: On curvature and hyperbolicity of monotone Hamiltonian systems

abstract: Assume that a Hamiltonian system is monotone. In this paper, we give several characterizations on when such a system is Anosov. Assuming that a monotone Hamiltonian system has no conjugate point, we show that there are two distributions which are invariant under the Hamiltonian flow. We show that a monotone Hamiltonian flow without conjugate point is Anosov if and only if these distributions are transversal. We also show that if the reduced curvature of the Hamiltonian system is non-positive, then the flow is Anosov if and only if the reduced curvature is negative somewhere along each trajectory.


title: One-step synthesis of KxFe2-ySe2 single crystals for high critical current density

abstract: We have established a simple process that allows for the one-step synthesis of KxFe2-ySe2 single crystals, which exhibit high critical current density Jc. The post annealing and quenching technique has improved the homogeneity of as-grown crystals, resulting in full shielding of the external magnetic field. The quenched crystals show a superconducting transition at Tconset = 32.9 K and Tczero = 32.1 K. The upper critical fields \mu_{0}Hc2(0) for H//ab and H//c are estimated to be ~206 and ~50 T, respectively. The critical current densities Jc for H//ab and H//c reach as high as 1.0\times10^{5} and 3.4\times10^{4} A/cm2 at 5 K. Furthermore, Jc exhibits a high field performance and a significantly weak temperature dependence up to 5 T, suggesting strong pinning. These results demonstrate that KxFe2-ySe2 would be a promising candidate material for practical applications.


title: Chandra pulsar survey (ChaPS)

abstract: Taking advantage of the high sensitivity of the Chandra Advanced CCD Imaging Spectrometer, we have conducted a snap-shot survey of pulsars previously undetected in X-rays. We detected 12 pulsars and established deep flux limits for 11 pulsars. Using these new results, we revisit the relationship between the X-ray luminosity, L_psr_x, and spin-down power, Edot. We find that the obtained limits further increase the extremely large spread in the non-thermal X-ray efficiencies, eta_psr_x=L_psr_x/Edot, with some of them being now below 1e-5. Such a spread cannot be explained by poorly known distances or by beaming of pulsar radiation. We also find evidence of a break in the dependence of L_psr_x on Edot, such that pulsars become more X-ray efficient at Edot<~ 1e34-1e35 erg/s. We examine the relationship between the gamma-ray luminosity, L_psr_g, and Edot, which exhibits a smaller scatter compared to that in X-rays. This confirms that the very large spread in the X-ray efficiencies cannot be explained just by beaming because the gamma-ray emission is generally expected to be beamed stronger than the X-ray emission. Intriguingly, there is also an indication of a break in the L_psr_g(Edot) dependence at Edot~1e35 erg/s, with lower-Edot pulsars becoming less gamma-ray efficient. We also examine the distance-independent L_psr_f/L_psr_x ratio as a function of Edot for a sample of gamma-ray pulsars observed by Chandra and find that it peaks at Edot~1e35 erg/s, showing that the breaks cannot originate from poorly measured distances. We discuss the implications of our findings for existing models of magnetospheric emission and venues for further exploration.


title: Honeycomb Lattice Potentials and Dirac Points

abstract: We prove that the two-dimensional Schroedinger operator with a potential having the symmetry of a honeycomb structure has dispersion surfaces with conical singularities (Dirac points) at the vertices of its Brillouin zone. No assumptions are made on the size of the potential. We then prove the robustness of such conical singularities to a restrictive class of perturbations, which break the honeycomb lattice symmetry. General small perturbations of potentials with Dirac points do not have Dirac points; their dispersion surfaces are smooth. The presence of Dirac points in honeycomb structures is associated with many novel electronic and optical properties of materials such as graphene.

categories: math-ph cond-mat.mtrl-sci math.AP math.MP quant-ph


title: Vector- and Pseudoscalar-baryon coupled channel systems

abstract: In this manuscript, I will report the details of our recent work on the vector meson-baryon (VB) interaction, which we studied with the motivation of finding dynamical generation of resonances in the corresponding systems. We started our study by building a formalism based on the hidden local symmetry and calculating the leading order contributions to the scattering equations by summing the diagrams with: (a) a vector meson exchange in the t-channel (b) an octet baryon exchange in the s-, u-channels and (c) a contact interaction arising from the part of the vector meson-baryon Lagrangian which is related to the anomalous magnetic moment of the baryons. We find the contribution from all these sources, except the s-channel, to be important. The amplitudes obtained by solving the coupled channel Bethe-Salpeter equations for the systems with total strangeness zero, show generation of one isospin 3/2, spin 1/2 resonance and three isospin 1/2 resonances: two with spin 3/2 and one with spin 1/2. We identify these resonances with $\Delta$ (1900) $S_{31}$, $N^*$(2080) $D_{13}$, $N^*$(1700) $D_{13}$, and $N^*$(2090) $S_{11}$, respectively. We have further extended our study by including pseudoscalar meson-baryon (PB) as the coupled channels of VB systems. For this, we obtain the PB $\rightarrow$ VB amplitudes by using the Kroll-Ruddermann term where, considering the vector meson dominance phenomena, the photon is replaced by a vector meson. The calculations done within this formalism reveal a very strong coupling of the VB channels to the low-lying resonances like $\Lambda$(1405) and $\Lambda$(1670), which can have important implications on certain reactions producing them. In addition to this, we find that the effect of coupling the higher mass states to the lighter channels is not restricted to increasing the width of those states, it can be far more strong.


title: A functional model for pure $\Gamma$-contractions

abstract: A pair of commuting operators $(S,P)$ defined on a Hilbert space $\mathcal H$ for which the closed symmetrized bidisc $$ \Gamma= \{(z_1+z_2,z_1z_2):: |z_1|\leq 1,\, |z_2|\leq 1 \}\subseteq \mathbb C^2, $$ is a spectral set is called a $\Gamma$-contraction in the literature. A $\Gamma$-contraction $(S,P)$ is said to be pure if $P$ is a pure contraction, i.e, ${P^*}^n \rightarrow 0$ strongly as $n \rightarrow \infty $. Here we construct a functional model and produce a complete unitary invariant for a pure $\Gamma$-contraction. The key ingredient in these constructions is an operator, which is the unique solution of the operator equation $$ S-S^*P=D_PXD_P, \textup{where} X\in \mathcal B(\mathcal D_P), $$and is called the fundamental operator of the $\Gamma$-contraction $(S,P)$. We also discuss some important properties of the fundamental operator.


title: Stabilized Lattice Boltzmann-Enskog method for compressible flows and its application to one and two-component fluids in nanochannels

abstract: A numerically stable method to solve the discretized Boltzmann-Enskog equation describing the behavior of non ideal fluids under inhomogeneous conditions is presented. The algorithm employed uses a Lagrangian finite-difference scheme for the treatment of the convective term and a forcing term to account for the molecular repulsion together with a Bhatnagar-Gross-Krook relaxation term. In order to eliminate the spurious currents induced by the numerical discretization procedure, we use a trapezoidal rule for the time integration together with a version of the two-distribution method of He et al. (J. Comp. Phys 152, 642 (1999)). Numerical tests show that, in the case of one component fluid in the presence of a spherical potential well, the proposed method reduces the numerical error by several orders of magnitude. We conduct another test by considering the flow of a two component fluid in a channel with a bottleneck and provide information about the density and velocity field in this structured geometry.


title: The internally 4-connected binary matroids with no M(K5\e)-minor

abstract: Let AG(3,2)xU(1,1) denote the binary matroid obtained from the direct sum of AG(3,2) and a coloop by completing the 3-point lines between every element in AG(3,2) and the coloop. We prove that every internally 4-connected binary matroid that does not have a minor isomorphic to M(K5\e) is isomorphic to a minor of (AG(3,2)xU(1,1))*.


title: Energy dependence of inclusive photon elliptic flow in heavy-ion collision models

abstract: We present a comparison of inclusive photon elliptic flow parameter (v_{2}) measured at RHIC and SPS high energy heavy-ion collision experiments to calculations done using the AMPT and UrQMD models. The new results discussed includes the comparison of the model calculations of photon v_{2} to corresponding measurements at the forward rapidities. We observe that the AMPT model which includes partonic interactions and quark coalescence as a mechanism of hadronization is in good agreement with the measurements even at forward rapidities (2.3 < \eta < 3.9) at RHIC as was previously observed for measurements at midrapidity. At the top SPS energy the contribution from partonic effects are smaller than that at RHIC energy, based on the comparison of the measured photon v_{2} with those from the AMPT default and UrQMD model calculations. We find that if the measurements in RHIC beam energy scan (BES) and LHC energies would require an energy dependent partonic cross section in the AMPT models, then the observed longitudinal scaling of v_{2} at top RHIC energies (19.6-200 GeV) will be violated. We also discuss the relation between the inclusive photon v_{2} and those of their parent \pi^{0}'s for the beam energies of 7.7 GeV to 2.76 TeV. The model results show that the transverse momentum (p_{\mathrm T}) integrated v_{2} of \pi^{0} is larger by about 44% relative to those of the inclusive photons. Finally we present the expectations of inclusive photon v_{2} for the RHIC beam energy scan (BES) program and LHC from the transport models, so that they can be compared to corresponding measurements using the data already collected at RHIC and LHC.


title: The Complex Energy Method Applied to the Nd Scattering with a Model Three-Body Force

abstract: Using the complex energy method, the problem of nucleon-deuteron scattering is solved with a simple three-body force having a separable form. Our results are compared with the results of modern direct two-variable calculations and a good agreement is found. This forms a firm base for other applications of the complex energy method.


title: Multiscale non-adiabatic dynamics with radiative decay, case study on the post-ionization fragmentation of rare-gas tetramers

abstract: In this supplementary material, we recollect, for reader's convenience, the general scheme of suggested multiscale model (Sec. 1), and basic informations about approaches used for pilot study: a detailed description of the interaction model (Sec. 2) and dynamical methods used for the dark dynamics step (Sec. 3) reported previously in two preceding studies [1, 2]. In addition, a detailed description of the treatment of radiative processes is also given (Sec. 4).


title: Gone with the Wind ON_Mars (GOWON): A Wind-Driven Networked System of Mobile Sensors on Mars

abstract: We propose a revolutionary way of studying the sur-face of Mars using a wind-driven network of mobile sensors- Gone with the Wind ON_Mars (GOWON). GOWON is envisioned to be a scalable, 100% self energy-generating and distributed system that allows in-situ mapping of a wide range of phenomena in a much larger portion of the surface of Mars compared to earlier missions. It could radically improve the possibility of finding rare phenomena like bio signatures through random wind-driven search. It could explore difficult terrains that were beyond the reach of previous missions, such as regions with very steep slopes, cluttered surfaces and/or sand dunes; GOWON is envisioned as an on going mission with a long life span. It could achieve any of NASA's scientific objectives on Mars in a cost-effective way, leaving a long lasting sensing and searching infrastructure on Mars. GOWON is a 2012 Step B invitee for NASA Innovative Advanced Concept (NIAC). It addresses the challenge area of the Mars Surface System Capabilities area. We believe the challenge to be near-term, i.e., 2018-2024.