title: New results on the search for spin-exotic mesons with COMPASS

abstract: The COMPASS fixed-target experiment at the CERN-SPS studies the structure and spectrum of hadrons. One important goal using hadron beams is the search for new states, in particular spin-exotic mesons and glueballs. As a first input to the puzzle, COMPASS observed a significant $J^{PC}$ spin-exotic signal in the 2004 pilot run data (190\,GeV/$c$ $\pi^{-}$ beam, Pb target) in three charged pion final states consistent with the disputed $\pi_1(1600)$. We started our hadron spectroscopy programme in 2008 by collecting very high statistics using a 190 GeV/$c$ negative pion beam scattered off a liquid hydrogen (proton) target. The current status and new results from the 2008 data on the search for the $\pi_1(1600)$ resonance with exotic $J^{PC}=1^{-+}$ quantum numbers obtained from partial-wave analyses of the $\rho\pi$ and $\eta'\pi$ decay channels are presented.


title: Electrolytic Synthesis and Characterizations of Silver Nanopowder

abstract: This work reports a simple, novel, cost effective and eco-friendly electrolytic synthesis of silver nanoparticles using AgNO3 as metal precursor. The synthesis rate is much faster than other methods and this approach is suitable for large scale production. They are characterized by XRD, SEM and FT-IR techniques to analyze size, morphology and functional groups. XRD studies reveal a high degree of crystallinity and monophasic Ag nanoparticles. Their particle size is found to be 24 nm and specific surface area (SSA) is 24 m2/g. Analysis of Ag nanoparticles SSA reports that increasing their SSA improves their antibacterial actions. Microbiology assay founds that Ag nanoparticles are effective against E.coli and B.megaterium bacteria. SSA of bacteria analysis reveals that it plays a major role while reacting with antimicrobial agents.


title: Extended Lagrangian free energy molecular dynamics

abstract: Extended free energy Lagrangians are proposed for first principles molecular dynamics simulations at finite electronic temperatures for plane-wave pseudopotential and local orbital density matrix based calculations. Thanks to the extended Lagrangian description the electronic degrees of freedom can be integrated by stable geometric schemes that conserve the free energy. For the local orbital representations both the nuclear and electronic forces have simple and numerically efficient expressions that are well suited for reduced complexity calculations. A rapidly converging recursive Fermi operator expansion method that does not require the calculation of eigenvalues and eigenfunctions for the construction of the fractionally occupied density matrix is discussed. An efficient expression for the Pulay force that is valid also for density matrices with fractional occupation occurring at finite electronic temperatures is also demonstrated.


title: Modern observations of Hubble's first-discovered Cepheid in M31

abstract: We present a modern ephemeris and modern light curve of the first-discovered Cepheid variable in M31, Edwin Hubble's M31-V1. Observers of the American Association of Variable Star Observers undertook these observations during the latter half of 2010. The observations were in support of an outreach program by the Space Telescope Science Institute's Hubble Heritage project, but the resulting data are the first concentrated observations of M31-V1 made in modern times. AAVSO observers obtained 214 V-band, Rc-band, and unfiltered observations from which a current ephemeris was derived. The ephemeris derived from these observations is JD(Max) = 2455430.5(+/-0.5) + 31.4 (+/-0.1) E. The period derived from the 2010 data are in agreement with the historic values of the period, but the single season of data precludes a more precise determination of the period or measurement of the period change using these data alone. However, using an ephemeris based upon the period derived by Baade and Swope we are able to fit all of the observed data acceptably well. Continued observations in the modern era will be very valuable in linking these modern data with data from the 1920s-30s and 1950s, and will enable us to measure period change in this historic Cepheid. In particular, we strongly encourage intensive observations of this star around predicted times of maximum to constrain the date of maximum to better than 0.5 days.


title: Scalar Hair from a Derivative Coupling of a Scalar Field to the Einstein Tensor

abstract: We consider a gravitating system of vanishing cosmological constant consisting of an electromagnetic field and a scalar field coupled to the Einstein tensor. A Reissner-Nordstr\"om black hole undergoes a second-order phase transition to a hairy black hole of generally anisotropic hair at a certain critical temperature which we compute. The no-hair theorem is evaded due to the coupling between the scalar field and the Einstein tensor. Within a first order perturbative approach we calculate explicitly the properties of a hairy black hole configuration near the critical temperature and show that it is energetically favorable over the corresponding Reissner-Nordstr\"om black hole.


title: Isoparametric hypersurfaces in Damek-Ricci spaces

abstract: We construct uncountably many isoparametric families of hypersurfaces in Damek-Ricci spaces. We characterize those of them that have constant principal curvatures by means of the new concept of generalized Kahler angle. It follows that, in general, these examples are inhomogeneous and have nonconstant principal curvatures. We also find new cohomogeneity one actions on quaternionic hyperbolic spaces, and an isoparametric family of inhomogeneous hypersurfaces with constant principal curvatures in the Cayley hyperbolic plane.


title: Environment and properties of emitting electrons in blazar jets: Mrk 421 as a laboratory

abstract: Here we report our recent study on the spectral energy distribution (SED) of the high frequency BL Lac object Mrk 421 in different luminosity states. We used a full-fledged chi2-minimization procedure instead of more commonly used "eyeball" fit to model the observed flux of the source (from optical to very high energy), with a Synchrotron-Self-Compton (SSC) emission mechanism. Our study shows that the synchrotron power and peak frequency remain constant with varying source activity, and the magnetic field decreases with the source activity while the break energy of electron spectrum and the Doppler factor increase. Since a lower magnetic field and higher density of electrons result in increased electron-photon scattering efficiency, the Compton power increases, so does the total emission.


title: Rugged free-energy landscapes in disordered spin systems

abstract: This thesis is an attempt to provide a new outlook on complex systems, as well as some physical answers for certain models, taking a computational approach. We have focused on disordered systems, addressing two traditional problems in three spatial dimensions: the Edwards-Anderson spin glass and the Diluted Antiferromagnet in a Field (the physical realisation of the random-field Ising model). These systems have been studied by means of large-scale Monte Carlo simulations, exploiting a variety of platforms, which include the Janus special-purpose supercomputer. Two main themes are explored throughout: a) the relationship between the (experimentally unreachable) equilibrium phase and the non-equilibrium evolution and b) the computation and efficient treatment of rugged free-energy landscapes. We perform a thorough study of the low-temperature phase of the D=3 Edwards-Anderson spin glass, where we establish a time-length dictionary and a finite-time scaling formalism to link, in a quantitative way, the experimental non-equilibrium regime and the finite-size equilibrium phase. At the experimentally relevant scales, the replica symmetry breaking theory emerges as the appropriate theoretical picture. We also introduce Tethered Monte Carlo, a general strategy for the study of systems with rugged free-energy landscapes. This formalism provides a general method to guide the exploration of the configuration space by constraining one or more reaction coordinates. From these tethered simulations, the Helmholtz potential associated to the reaction coordinates is reconstructed, yielding all the information about the system. We use this method to provide a comprehensive picture of the critical behaviour in the Diluted Antiferromagnet in a Field.


title: Phase diagrams of $\rm La_{1-x}Ca_xMnO_3$ in Double Exchange Model with added antiferromagnetic and Jahn-Teller interaction

abstract: The phase diagram of the multivalent manganites $\rm La_{1-x}Ca_xMnO_3$, in space of temperature and doping $x$, is a challenge for the theoretical physics. It is an important test for the model used to study these compounds and the method of calculation. To obtain theoretically this diagram for $x<0.5$, we consider the two-band Double Exchange Model for manganites with added Jahn-Teller coupling and antiferromagnetic Heisenberg term. In order to calculate Curie and N\'{e}el temperatures we derive an effective Heisenberg model for a vector which describes the local orientation of the total magnetization of the system. The exchange constants of this model are different for different space directions and depend on the density of $e_g$ electrons, antiferromagnetic constants and the Jahn-Teller energy. To reproduce the well known phase transitions from A-type antiferromagnetism to ferromagnetism at low $x$ and C-type antiferromagnetism to G-type antiferromagnetism at large $x$, we argue that the antiferromagnetic exchange constants should depend on the lattice direction. We show that ferromagnetic to A-type antiferromagnetic transition results from the Jahn-Teller distortion. Accounting adequately for the magnon-magnon interaction, Curie and N\'{e}el temperatures are calculated. The results are in very good agreement with the experiment and provide values for the model parameters, which best describe the behavior of the critical temperature for $x<0.5$.


abstract: The necessity of a theory of General Topology and, most of all, of Algebraic Topology on locally finite metric spaces comes from many areas of research in both Applied and Pure Mathematics: Molecular Biology, Mathematical Chemistry, Computer Science, Topological Graph Theory and Metric Geometry. In this paper we propose the basic notions of such a theory and some applications: we replace the classical notions of continuous function, homeomorphism and homotopic equivalence with the notions of NPP-function, NPP-local-isomorphism and NPP-homotopy (NPP stands for Nearest Point Preserving); we also introduce the notion of NPP-isomorphism. We construct three invariants under NPP-isomorphisms and, in particular, we define the fundamental group of a locally finite metric space. As first applications, we propose the following: motivated by the longstanding question whether there is a purely metric condition which extends the notion of amenability of a group to any metric space, we propose the property SN (Small Neighborhood); motivated by some applicative problems in Computer Science, we prove the analog of the Jordan curve theorem in $\mathbb Z^2$; motivated by a question asked during a lecture at Lausanne, we extend to any locally finite metric space a recent inequality of P.N.Jolissaint and Valette regarding the $\ell_p$-distortion.


title: Limiting distribution of maximal crossing and nesting of Poissonized random matchings

abstract: The notion of $r$-crossing and $r$-nesting of a complete matching was introduced and a symmetry property was proved by Chen et al. [Trans. Amer. Math. Soc. 359 (2007) 1555-1575]. We consider random matchings of large size and study their maximal crossing and their maximal nesting. It is known that the marginal distribution of each of them converges to the GOE Tracy-Widom distribution. We show that the maximal crossing and the maximal nesting becomes independent asymptotically, and we evaluate the joint distribution for the Poissonized random matchings explicitly to the first correction term. This leads to an evaluation of the asymptotic of the covariance. Furthermore, we compute the explicit second correction term in the distribution function of two objects: (a) the length of the longest increasing subsequence of Poissonized random permutation and (b) the maximal crossing, and hence also the maximal nesting, of Poissonized random matching.


title: "Thermodynamique cach\'ee des particules" and the quantum potential

abstract: According to de Broglie, temperature plays a basic role in quantum Hamilton-Jacobi theory. Here we show that a possible dependence on the temperature of the integration constants of the relativistic quantum Hamilton-Jacobi may lead to corrections to the dispersion relations. The change of the relativistic equations is simply described by means of a thermal coordinate.


title: Transfer and scattering of wave packets by a nonlinear trap

abstract: In the framework of a one-dimensional model with a tightly localized self-attractive nonlinearity, we study the formation and transfer (dragging) of a trapped mode by "nonlinear tweezers", as well as the scattering of coherent linear wave packets on the stationary localized nonlinearity. The use of the nonlinear trap for the dragging allows one to pick up and transfer the relevant structures without grabbing surrounding "garbage". A stability border for the dragged modes is identified by means of of analytical estimates and systematic simulations. In the framework of the scattering problem, the shares of trapped, reflected, and transmitted wave fields are found. Quasi-Airy stationary modes with a divergent norm, that may be dragged by the nonlinear trap moving at a constant acceleration, are briefly considered too.


title: Multiparton interactions and multiparton distributions in QCD

abstract: After a brief recapitulation of the general interest of parton densities, we discuss multiple hard interactions and multiparton distributions. We report on recent theoretical progress in their QCD description, on outstanding conceptual problems and on possibilities to use multiparton distributions as a laboratory to test and improve our understanding of hadron structure.


title: The Minimal Dimensionless Standard Model (MDSM) and its Cosmology

abstract: Consider the minimal renormalizable extension of the Standard Model with purely dimensionless couplings, successful electroweak symmetry breaking (via the Coleman-Weinberg mechanism) and a see-saw mechanism for neutrino mass: we will call this the Minimal Dimensionless Standard Model (MDSM). In fact, 3 closely related models fit the bill: MDSM_1, MDSM_2 and MDSM_3. We analyze the theoretical and observational constraints on these models. We argue that, when they are minimally coupled to gravity, they can accomplish several important cosmological tasks (inflation, dark matter, leptogenesis) in a way that is economical, predictive and tightly woven into the fabric of known physics. One of the models (MDSM_3), which includes an extra U(1)_{B-L} gauge symmetry, seems particularly promising.


title: Mathematics of Knowledge Refinement: Probabilistic Arithmetic, with no unknowns and no infinity. Part I. Generalized Probabilistic Arithmetic. Basic definitions and properties

abstract: An approach to build Probabilistic Arithmetic in which initial values of all correlated random variables are known, but with varying degrees of accuracy. As a result of the proposed Probabilistic Arithmetic operations, variable values, degrees of their accuracy and correlations are refined. Probabilistic Generalized Addition (PGA) and Probabilistic Generalized Multiplication (PGM) operations on correlated random variables are defined and their basic properties identified and described: \bullet Proposed PGA and PGM operations possess inverse operations - subtraction and division. \bullet There is no difference between direct and inverse operations: addition and subtraction, multiplication and division (this is why these operations are called "Generalized"). \bullet Division by approximately zero is possible and the result never equals to \infty, making this approach promising in computational and computer mathematics. \bullet Unlike the usual hyperbola, which, when the argument is changing from + \infty to 0, has a gap at 0, the Generalized Probabilistic Hyperbola, under certain combination of initial accuracies, is continuous at 0. First, as the argument changes from +\infty to 0, it behaves like a typical hyperbola, monotonically increasing. However, after reaching certain maximum value, it starts to decrease to 0 monotonically. Last property gives a hope that the proposed approach might be promising in Quantum Physics, as it would allow a more adequate description of macro and micro physics, as well as of the transition from one to the other.


title: Confinement: a real-time visualization

abstract: Due to the mechanism of confinement, as known from quantum chromodynamics, it is difficult to observe individual particles carrying fractional quantum number (e.g. quark with fractional electric charge). A condensed matter example of fractionalized particles is spinons in quasi-one-dimensional spin systems, which are domain walls in the background of Neel configurations carrying spin-1/2. Using the time-evolving block decimation algorithm, we visualize the nontrivial domain wall dynamics induced by the confine mechanism in a two-leg spin-1/2 ladder. It can be illustrated by a simple single-particle picture of Bloch oscillation, not only qualitatively but also quantitatively. We propose the experimental realization and the real time detection of the domain wall dynamics in the ultra-cold boson systems of $^{87}$Rb.


title: Toward Higher Spin dS3/CFT2

abstract: I take steps toward the construction of a CFT dual to Vasiliev's higher spin gravity in three dimensional de Sitter space. There are two main claims. The first is that higher spin de Sitter symmetries are related to extended Virasoro symmetries, as in AdS; this is verified explicitly for the case of W_3 asymptotic symmetry. The associated chiral algebra has imaginary central charge. The second (conjectural) claim, inspired by work of Gaberdiel and Gopakumar in AdS_3/CFT_2, is that an appropriate CFT can be identified as an exotic non-unitary WZW coset model at complex level.


title: Scattering for solutions of NLS in the exterior of a 2d star--shaped obstacle

abstract: We prove that solutions to non-linear Schr\"odinger equations in two dimensions and in the exterior of a bounded and smooth star-shaped obstacle scatter in the energy space. The non-linear potential is defocusing and grows at least as the quintic power.


title: Quasi-potentials and regularization of currents, and applications

abstract: Let $Y$ be a compact K\"ahler manifold. We show that the weak regularization $K_n$ of Dinh and Sibony for the diagonal $\Delta_Y$ (see Section 2 for more detail) is compatible with wedge product in the following sense: If $T$ is a positive $dd^c$-closed $(p,p)$ current and $\theta$ is a smooth $(q,q)$ form then there is a sequence of positive $dd^c$-closed $(p+q,p+q)$ currents $S_n$ whose masses converge to 0 so that $-S_n\leq K_n(T\wedge \theta)-K_n(T)\wedge \theta \leq S_n$ for all $n$. We also prove a result concerning the quasi-potentials of positive closed currents. We give two applications of these results. First, we prove a corresponding compatibility with wedge product for the pullback operator defined in our previous paper. Second, we define an intersection product for positive $dd^c$-closed currents. This intersection is symmetric and has a local nature.


title: Free Resolutions and Sparse Determinantal Ideals

abstract: A sparse generic matrix is a matrix whose entries are distinct variables and zeros. Such matrices were studied by Giusti and Merle who computed some invariants of their ideals of maximal minors. In this paper we extend these results by computing a minimal free resolution for all such sparse determinantal ideals. We do so by introducing a technique for pruning minimal free resolutions when a subset of the variables is set to zero. Our technique correctly computes a minimal free resolution in two cases of interest: resolutions of monomial ideals, and ideals resolved by the Eagon-Northcott Complex. As a consequence we can show that sparse determinantal ideals have a linear resolution over the integers, and that the projective dimension depends only on the number of columns of the matrix which are identically zero. Finally, we show that all such ideals have the property that regardless of the term order chosen, the Betti numbers of the ideal and its initial ideal are the same. In particular the nonzero generators of these ideals form a universal Gr\"obner basis.


title: Generating Functionals and Lagrangian PDEs

abstract: We introduce the concept of Type-I/II generating functionals defined on the space of boundary data of a Lagrangian field theory. On the Lagrangian side, we define an analogue of Jacobi's solution to the Hamilton-Jacobi equation for field theories, and we show that by taking variational derivatives of this functional, we obtain an isotropic submanifold of the space of Cauchy data, described by the so-called multisymplectic form formula. We also define a Hamiltonian analogue of Jacobi's solution, and we show that this functional is a Type-II generating functional. We finish the paper by defining a similar framework of generating functions for discrete field theories, and we show that for the linear wave equation, we recover the multisymplectic conservation law of Bridges.


title: The Dual Potential, the involution kernel and Transport in Ergodic Optimization

abstract: Consider the shift $\sigma$ acting on the Bernoulli space $\Sigma={1,2,...,n}^\mathbb{N}$. We denote $\hat{\Sigma}= {1,2,...,n}^\mathbb{Z}$. We analyze several properties of the maximizing probability $\mu_{\infty,A}$ of a Holder potential $A: \Sigma \to \mathbb{R}$. Associated to $A(x)$, via the involution kernel, $W: \hat{\Sigma} \to \mathbb{R}$, it is known that can we get the dual potential $A^*(y)$, where $(x,y)\in \hat{\Sigma}$. Consider $\mu_{\infty, A^*}$ a maximizing probability for $A^*$. We would like to consider the transport problem from $\mu_{\infty,A}$ to $\mu_{\infty,A^*}$. In this case, it is natural to consider the cost function $c(x,y) = I(x) - W(x,y) +\gamma $, where $I$ is the deviation function. The pair of functions for the Kantorovich Transport dual Problem are $(-V,-V^*$), where we denote the two calibrated sub-actions by $V$ and $V^*$, respectively, for $A$ and $A^*$ for $\mu_{\infty,A}$. We analyze the graph property for the optimal plan $\hat{\mu}$.

categories: math.DS cond-mat.stat-mech math-ph math.MP math.OC math.PR


title: Statistical model uncertainty and OPERA-like time-of-flight measurements

abstract: Time-of-flight measurements such as the OPERA and MINOS experiments rely crucially on statistical analysis (as well as many other ingredients) for their conclusions. The nature of these experiments leads to a simple class of statistical models for the results; however, which model in the class is appropriate is not known exactly, as this depends on information obtained experimentally, which is subject to noise and other errors. To obtain robust conclusions, this problem, known as "model uncertainty," needs to be addressed, with quantitative bounds on the effect such uncertainty may have on the final result. The OPERA (and MINOS) analysis appears to take steps to mitigate the effects of model uncertainty, though without quantifying any remaining effect. We describe one of the strategies used (averaging individual probability distributions), and point out a potential source of error if this is not carried out correctly. We then argue that the correct version of this strategy is not the most effective, and suggest possible alternatives. These alternatives may give more accurate statistical results using the same data, allowing, for example, more accurate determination of the dependence of the anomalous time shift on energy. Which strategies work and how well can only be evaluated with access to the full data. Whether or not the anomalous result from OPERA turns out to be confirmed, we believe that techniques such as those presented here may be appropriate for the analysis of other timing experiments of this type.


title: Panchromatic observations of the textbook GRB 110205A: constraining physical mechanisms of prompt emission and afterglow

abstract: We present a comprehensive analysis of a bright, long duration (T90 ~ 257 s) GRB 110205A at redshift z= 2.22. The optical prompt emission was detected by Swift/UVOT, ROTSE-IIIb and BOOTES telescopes when the GRB was still radiating in the gamma-ray band. Nearly 200 s of observations were obtained simultaneously from optical, X-ray to gamma-ray, which makes it one of the exceptional cases to study the broadband spectral energy distribution across 6 orders of magnitude in energy during the prompt emission phase. By fitting the time resolved prompt spectra, we clearly identify, for the first time, an interesting two-break energy spectrum, roughly consistent with the standard GRB synchrotron emission model in the fast cooling regime. Although the prompt optical emission is brighter than the extrapolation of the best fit X/gamma-ray spectra, it traces the gamma-ray light curve shape, suggesting a relation to the prompt high energy emission. The synchrotron + SSC scenario is disfavored by the data, but the models invoking a pair of internal shocks or having two emission regions can interpret the data well. Shortly after prompt emission (~ 1100 s), a bright (R = 14.0) optical emission hump with very steep rise (alpha ~ 5.5) was observed which we interpret as the emission from the reverse shock. It is the first time that the rising phase of a reverse shock component has been closely observed. The full optical and X-ray afterglow lightcurves can be interpreted within the standard reverse shock (RS) + forward shock (FS) model. In general, the high quality prompt emission and afterglow data allow us to apply the standard fireball shock model to extract valuable information about the GRB including the radiation mechanism, radius of prompt emission R, initial Lorentz factor of the outflow, the composition of the ejecta, as well as the collimation angle and the total energy budget.


title: A topological interpretation of the walk distances

abstract: The walk distances in graphs have no direct interpretation in terms of walk weights, since they are introduced via the \emph{logarithms} of walk weights. Only in the limiting cases where the logarithms vanish such representations follow straightforwardly. The interpretation proposed in this paper rests on the identity $\ln\det B=\tr\ln B$ applied to the cofactors of the matrix $I-tA,$ where $A$ is the weighted adjacency matrix of a weighted multigraph and $t$ is a sufficiently small positive parameter. In addition, this interpretation is based on the power series expansion of the logarithm of a matrix. Kasteleyn (1967) was probably the first to apply the foregoing approach to expanding the determinant of $I-A$. We show that using a certain linear transformation the same approach can be extended to the cofactors of $I-tA,$ which provides a topological interpretation of the walk distances.


title: The Local Field Factor and Microscopic Cascading: A Self-Consistent Method Applied to Confined Systems of Molecules

abstract: We use a simplified self-consistent method to address nonlinear-optical cascading phenomena, which shows added microscopic cascading contributions in high-ordered nonlinear susceptibilities through fifth order. These cascading terms in the microscopic regime encompass all possible scalar cascading configurations. The imposition of geometric constraints further influences the predicted cascading contributions and opens up additional design parameters for nonlinear-optical materials. These results are used in approximating the effective fifth-order susceptibility in thin films of C60 monomers of varying thickness and concentration. This paper contains the corrections to the original paper that appeared in the Journal of Physics B as reflected in the content of the corrigendum that followed.


title: Does the Superluminal Neutrino Uncover Torsion?

abstract: I investigate the possibility of the propagation of neutrino with superluminal speed through matter in the context of the relation between gravity, spin and torsion. Using a lemma of Penrose and earlier works on the relation between spin, torsion and gravity I glimpse on a frame work in which superluminal speed of the neutrinos moving through matter become possible. In presence of torsoin neutrinos are found to follow spacelike geodesics by tunneling through the light cone into the spacelike region and consequently appear to have superluminal speed in our timelike world, whereas photons always follow null geodesics. This frame work may set new frontiers for spacetime physics.


title: Partial quasi-morphisms and quasi-states on cotangent bundles, and symplectic homogenization

abstract: For a closed connected manifold N, we construct a family of functions on the Hamiltonian group G of the cotangent bundle T^*N, and a family of functions on the space of smooth functions with compact support on T^*N. These satisfy properties analogous to those of partial quasi-morphisms and quasi-states of Entov and Polterovich. The families are parametrized by the first real cohomology of N. In the case N=T^n the family of functions on G coincides with Viterbo's symplectic homogenization operator. These functions have applications to the algebraic and geometric structure of G, to Aubry-Mather theory, to restrictions on Poisson brackets, and to symplectic rigidity.


title: Electronic band structure of ferro-pnictide superconductors from ARPES experiment

abstract: ARPES experiments on iron based superconductors show that the differences between the measured and calculated electronic band structures look insignificant but can be crucial for understanding of the mechanism of high temperature superconductivity. Here we focus on those differences for 111 and 122 compounds and discuss the observed correlation of the experimental band structure with the superconductivity.


title: Detection of Outflowing and Extraplanar Gas in Disks in an Assembling Galaxy Cluster at z=0.37

abstract: We detect ionized gas characteristics indicative of winds in three disk-dominated galaxies that are members of a super-group at z=0.37 that will merge to form a Coma-mass cluster. All three galaxies are IR-luminous (L_IR > 4 x 10^{10} L_sun, SFR >8 M_sun per year) and lie outside the X-ray cores of the galaxy groups. We find that the most IR-luminous galaxy has strong blue and redshifted emission lines with velocities of ~ +/-200 km/s and a third, blueshifted (~ 900 km/s) component. This galaxy's line-widths (Hb, [OIII]5007, [NII], Ha) correspond to velocities of 100-1000 km/s. We detect extraplanar gas in two of three galaxies with SFR > 8 M_sun per year whose orientations are approximately edge-on and which have IFU spaxels off the stellar disk. IFU maps reveal that the extraplanar gas extends to r_h ~ 10 kpc; [NII] and Ha line-widths correspond to velocities of ~200-400 km/s in the disk and decrease to ~50-150 km/s above the disk. Multi-wavelength observations indicate that the emission is dominated by star formation. Including the most IR-luminous galaxy we find that 18% of supergroup members with SFR > 8 M_sun per year show ionized gas characteristics indicative of outflows. This is a lower limit as showing that gas is outflowing in the remaining, moderately inclined, galaxies requires a non-trivial decoupling of contributions to the emission lines from rotational and turbulent motion. Ionized gas mass loss in these winds is ~0.1 M_sun per year for each galaxy, although the winds are likely to entrain significantly larger amounts of mass in neutral and molecular gas.


title: Cosmological evolution of warm dark matter fluctuations I: Efficient computational framework with Volterra integral equations

abstract: We study the complete cosmological evolution of dark matter (DM) density fluctuations for DM particles that decoupled being ultrarelativistic during the radiation dominated era which is the case of keV scale warm DM (WDM). The new framework presented here can be applied to other types of DM and in particular we extend it to cold DM (CDM). The collisionless and linearized Boltzmann-Vlasov equations (B-V) for WDM and neutrinos in the presence of photons and coupled to the linearized Einstein equations are studied in detail in the presence of anisotropic stress with the Newtonian potential generically different from the spatial curvature perturbations. We recast this full system of B-V equations for DM and neutrinos into a system of coupled Volterra integral equations. These Volterra-type equations are valid both in the radiation dominated (RD) and matter dominated (MD) eras during which the WDM particles are ultrarelativistic and then nonrelativistic. This generalizes the so-called Gilbert integral equation only valid for nonrelativistic particles in the MD era. We succeed to reduce the system of four Volterra integral equations for the density and anisotropic stress fluctuations of DM and neutrinos into a system of only two coupled Volterra equations. The kernels and inhomogeneities in these equations are explicitly given functions. Combining the Boltzmann-Vlasov equations and the linearized Einstein equations constrain the initial conditions on the distribution functions and gravitational potentials. In the absence of neutrinos the anisotropic stress vanishes and the Volterra-type equations reduce to a single integral equation. These Volterra integral equations provide a useful and precise framework to compute the primordial WDM fluctuations over a wide range of scales including small scales up to k ~ 1/5 kpc.


title: Runaway massive stars from R136: VFTS 682 is very likely a "slow runaway"

abstract: We conduct a theoretical study on the ejection of runaway massive stars from R136 --- the central massive, star-burst cluster in the 30 Doradus complex of the Large Magellanic Cloud. Specifically, we investigate the possibility of the very massive star (VMS) VFTS 682 being a runaway member of R136. Recent observations of the above VMS, by virtue of its isolated location and its moderate peculiar motion, have raised the fundamental question whether isolated massive star formation is indeed possible. We perform the first realistic N-body computations of fully mass-segregated R136-type star clusters in which all the massive stars are in primordial binary systems. These calculations confirm that the dynamical ejection of a VMS from a R136-like cluster, with kinematic properties similar to those of VFTS 682, is common. Hence the conjecture of isolated massive star formation is unnecessary to account for this VMS. Our results are also quite consistent with the ejection of 30 Dor 016, another suspected runaway VMS from R136. We further note that during the clusters' evolution, mergers of massive binaries produce a few single stars per cluster with masses significantly exceeding the canonical upper-limit of 150 solar mass. The observations of such single super-canonical stars in R136, therefore, do not imply an IMF with an upper limit greatly exceeding the accepted canonical 150 solar mass limit, as has been suggested recently, and they are consistent with the canonical upper limit.


title: Resolving astrophysical uncertainties in dark matter direct detection

abstract: We study the impact of the assumed velocity distribution of galactic dark matter particles on the interpretation of results from nuclear recoil detectors. By converting experimental data to variables that make the astrophysical unknowns explicit, different experiments can be compared without implicit assumptions concerning the dark matter halo. We extend this framework to include the annual modulation signal, as well as multiple target elements. Recent results from DAMA, CoGeNT and CRESST-II can be brought into agreement if the velocity distribution is very anisotropic and thus allows a large modulation fraction. However constraints from CDMS and XENON cannot be evaded by appealing to such astrophysical uncertainties alone.


title: On Symmetric and Asymmetric Light Dark Matter

abstract: We examine cosmological, astrophysical and collider constraints on thermal dark matter (DM) with mass mX in the range 1 MeV to 10 GeV. Cosmic microwave background (CMB) observations, which severely constrain light symmetric DM, can be evaded if the DM relic density is sufficiently asymmetric. CMB constraints require the present anti-DM to DM ratio to be less than 2*10^{-6} (10^{-1}) for DM mass mX = 1 MeV (10 GeV) with ionizing efficiency factor f ~ 1. We determine the minimum annihilation cross section for achieving these asymmetries subject to the relic density constraint; these cross sections are larger than the usual thermal annihilation cross section. On account of collider constraints, such annihilation cross sections can only be obtained by invoking light mediators. These light mediators can give rise to significant DM self-interactions, and we derive a lower bound on the mediator mass from elliptical DM halo shape constraints. We find that halo shapes require a mediator with mass mphi > 4 * 10^{-2} MeV (40 MeV) for mX = 1 MeV (10 GeV). We map all of these constraints to the parameter space of DM-electron and DM-nucleon scattering cross sections for direct detection. For DM-electron scattering, a significant fraction of the parameter space is already ruled out by beam-dump and supernova cooling constraints.


title: Do all O stars form in star clusters?

abstract: The question whether or not massive stars can form in isolation or only in star clusters is of great importance for the theory of (massive) star-formation as well as for the stellar initial mass function of whole galaxies (IGIMF-theory). While a seemingly easy question it is rather difficult to answer. Several physical processes (e.g. star-loss due to stellar dynamics or gas expulsion) and observational limitations (e.g. dust obscuration of young clusters, resolution) pose severe challenges to answer this question. In this contribution we will present the current arguments in favour and against the idea that all O stars form in clusters.


title: Damping of tensor modes in inflation

abstract: We discuss the damping of tensor modes due to anisotropic stress in inflation. The effect is negligible in standard inflation and may be significantly large in inflation models that involve drastic production of free-streaming particles.


title: Modeling Collisional Cascades In Debris Disks: Steep Dust-Size Distributions

abstract: We explore the evolution of the mass distribution of dust in collision-dominated debris disks, using the collisional code introduced in our previous paper. We analyze the equilibrium distribution and its dependence on model parameters by evolving over 100 models to 10 Gyr. With our numerical models, we confirm that systems reach collisional equilibrium with a mass distribution that is steeper than the traditional solution by Dohnanyi (1969). Our model yields a quasi steady-state slope of n(m) ~ m^{-1.88} [n(a) ~ a^{-3.65}] as a robust solution for a wide range of possible model parameters. We also show that a simple power-law function can be an appropriate approximation for the mass distribution of particles in certain regimes. The steeper solution has observable effects in the submillimeter and millimeter wavelength regimes of the electromagnetic spectrum. We assemble data for nine debris disks that have been observed at these wavelengths and, using a simplified absorption efficiency model, show that the predicted slope of the particle mass distribution generates SEDs that are in agreement with the observed ones.


title: On the misalignment of the directly imaged planet {\beta} Pictoris b with the system's warped inner disk

abstract: The vertical warp in the debris disk Beta Pictoris -- an inclined inner disk extending into a flat outer disk -- has long been interpreted as the signpost of a planet on an inclined orbit. Direct images spanning 2004-2010 have revealed Beta Pictoris b, a planet with a mass and orbital distance consistent with this picture. However, it was recently reported that the orbit of planet b is aligned with the flat outer disk, not the inclined inner disk, and thus lacks the inclination to warp the disk. We explore three scenarios for reconciling the apparent misalignment of the directly imaged planet Beta Pictoris b with the warped inner disk of Beta Pictoris: observational uncertainty, an additional planet, and damping of planet b's inclination. We find that, at the extremes of the uncertainties, the orbit of Beta Pictoris b has the inclination necessary to produce the observed warp. We also find that if planet b were aligned with the flat outer disk, it would prevent another planet from creating a warp with the observed properties; therefore planet b itself must be responsible for the warp. Finally, planet b's inclination could have been damped by dynamical friction and still produce the observed disk morphology, but the feasibility of damping depends on disk properties and the presence of other planets. More precise observations of the orbit of planet b and the position angle of the outer disk will allow us to distinguish between the first and third scenario.


title: Phase separation in low-density neutron matter

abstract: Low-density neutron matter has been studied extensively for many decades, with a view to better understanding the properties of neutron-star crusts and neutron-rich nuclei. Neutron matter is beyond experimental control, but in the past decade it has become possible to create systems of fermionic ultracold atomic gases in a regime close to low-density neutron matter. In both these contexts pairing is significant, making simple perturbative approaches impossible to apply and necessitating ab initio microscopic simulations. Atomic experiments have also probed polarized matter. In this work, we study population-imbalanced neutron matter (possibly relevant to magnetars and to density functionals of nuclei) arriving at the lowest-energy configuration to date. For small to intermediate relative fractions the system turns out to be fully normal, while beyond a critical polarization we find phase coexistence between a partially polarized normal neutron gas and a balanced superfluid gas. As in cold atoms, a homogeneous polarized superfluid is close to stability but not stable with respect to phase separation. We also study the dependence of the critical polarization on the density.


title: Identifying Non-Resonant Kepler Planetary Systems

abstract: The Kepler mission has discovered a plethora of multiple transiting planet candidate exosystems, many of which feature putative pairs of planets near mean motion resonance commensurabilities. Identifying potentially resonant systems could help guide future observations and enhance our understanding of planetary formation scenarios. We develop and apply an algebraic method to determine which Kepler 2-planet systems cannot be in a 1st-4th order resonance, given the current, publicly available data. This method identifies when any potentially resonant angle of a system must circulate. We identify and list 70 near-resonant systems which cannot actually reside in resonance, assuming a widely-used formulation for deriving planetary masses from their observed radii and that these systems do not contain unseen bodies that affect the interactions of the observed planets. This work strengthens the argument that a high fraction of exoplanetary systems may be near resonance but not actually in resonance.


title: Cosmological evolution of warm dark matter fluctuations II: Solution from small to large scales and keV sterile neutrinos

abstract: We solve the cosmological evolution of warm dark matter (WDM) density fluctuations with the Volterra integral equations of paper I. In the absence of neutrinos, the anisotropic stress vanishes and the Volterra equations reduce to a single integral equation. We solve numerically this equation both for DM fermions decoupling at equilibrium and DM sterile neutrinos decoupling out of equilibrium. We give the exact analytic solution for the density fluctuations and gravitational potential at zero wavenumber. We compute the density contrast as a function of the scale factor a for a wide range of wavenumbers k. At fixed a, the density contrast grows with k for k < k_c while it decreases for k > k_c, (k_c ~ 1.6/Mpc). The density contrast depends on k and a mainly through the product k a exhibiting a self-similar behavior. Our numerical density contrast for small k gently approaches our analytic solution for k = 0. For fixed k < 1/(60 kpc), the density contrast generically grows with a while for k > 1/(60 kpc) it exhibits oscillations since the RD era which become stronger as k grows. We compute the transfer function of the density contrast for thermal fermions and for sterile neutrinos in: a) the Dodelson-Widrow (DW) model and b) in a model with sterile neutrinos produced by a scalar particle decay. The transfer function grows with k for small k and then decreases after reaching a maximum at k = k_c reflecting the time evolution of the density contrast. The integral kernels in the Volterra equations are nonlocal in time and their falloff determine the memory of the past evolution since decoupling. This falloff is faster when DM decouples at equilibrium than when it decouples out of equilibrium. Although neutrinos and photons can be neglected in the MD era, they contribute in the MD era through their memory from the RD era.


title: On 'Nothing'

abstract: Nothing---the absence of spacetime---can be either an endpoint of tunneling, as in the bubble of nothing, or a starting point for tunneling, as in the quantum creation of a universe. We argue that these two tunnelings can be treated within a unified framework, and that, in both cases, nothing should be thought of as the limit of anti-de Sitter space in which the curvature length approaches zero. To study nothing, we study decays in models with perturbatively stabilized extra dimensions, which admit not just bubbles of nothing---topology-changing transitions in which the extra dimensions pinch off and a hole forms in spacetime---but also a whole family of topology-preserving transitions that nonetheless smoothly hollow out and approach the bubble of nothing in one limit. The bubble solutions that are close to this limit, bubbles of next-to- nothing, give us a controlled setting in which to understand nothing. Armed with this understanding, we are able to embed proposed mechanisms for the reverse process, tunneling from nothing to something, within the relatively secure foundation of the Coleman-De Luccia formalism and show that the Hawking-Turok instanton does not mediate the quantum creation of a universe.


title: Numerical Simulations of Supernova Dust Destruction. II. Metal-Enriched Ejecta Knots

abstract: Following our previous work, we investigate through hydrodynamic simulations the destruction of newly-formed dust grains by sputtering in the reverse shocks of supernova remnants. Using an idealized setup of a planar shock impacting a dense, spherical clump, we implant a population of Lagrangian particles into the clump to represent a distribution of dust grains in size and composition. We vary the relative velocity between the reverse shock and ejecta clump to explore the effects of shock-heating and cloud compression. Because supernova ejecta will be metal-enriched, we consider gas metallicities from Z/Zsun = 1 to 100 and their influence on cooling properties of the cloud and the thermal sputtering rates of embedded dust grains. We post-process the simulation output to calculate grain sputtering for a variety of species and size distributions. In the metallicity regime considered in this paper, the balance between increased radiative cooling and increased grain erosion depends on the impact velocity of the reverse shock. For slow shocks (velocity less than or equal to 3000 km/s), the amount of dust destruction is comparable across metallicities, or in some cases is decreased with increased metallicity. For higher shock velocities (velocity greater than or equal to 5000 km/s), an increase in metallicity from Z/Zsun = 10 to 100 can lead to an additional 24% destruction of the initial dust mass. While the total dust destruction varies widely across grain species and simulation parameters, our most extreme cases result in complete destruction for some grain species and only 44% dust mass survival for the most robust species. These survival rates are important in understanding how early supernovae contribute to the observed dust masses in high-redshift galaxies.


title: Phase transitions in center-stabilized lattice gauge theories

abstract: We simulate four-dimensional center-stabilized lattice Yang-Mills theories on R^3 x S^1 with a newly developed pseudo-heatbath algorithm. We analyze the phase structure of such theories, namely the bulk transition and the spontaneous breaking of the center symmetry associated with the compact direction.


title: Linear series on curves: stability and Clifford index

abstract: We study concepts of stabilities associated to a smooth complex curve together with a linear series on it. In particular we investigate the relation between stability of the associated Dual Span Bundle and linear stability. Our result implies a stability condition related to the Clifford index of the curve. Furthermore, in some of the cases, we prove that a stronger stability holds: cohomological stability. Eventually using our results we obtain stable vector bundles of integral slope 3, and prove that they admit theta-divisors.


title: Construction of an NP Problem with an Exponential Lower Bound

abstract: In this paper we present a Hashed-Path Traveling Salesperson Problem (HPTSP), a new type of problem which has the interesting property of having no polynomial time solutions. Next we show that HPTSP is in the class NP by demonstrating that local information about sub-routes is insufficient to compute the complete value of each route. As a consequence, via Ladner's theorem, we show that the class NPI is non-empty.


title: Gap Symmetry in KFe_2As_2 and cos 4\theta gap component in LiFeAs

abstract: We revisit the issue of the gap symmetry in KFe2As2, which is an Fe-pnictide superconductor with only hole pockets. Previous theoretical studies mostly argued for a d-wave gap in KFe2As2 since transport and thermodynamic measurements point to the presence of the gap nodes. However, a d-wave gap is inconsistent with recent laser-based angle-resolved photoemission measurements. We propose the scenario for a nodal s-wave superconductivity induced by a non-magnetic intra-band and inter-band interactions between fermions near hole pockets. The superconducting gap changes sign between the hole pockets and has cos 4\theta angular dependence and accidental nodes on one or several hole pockets. We argue that strong angle dependence is the consequence of near-degeneracy between inter-pocket and intra-pocket interaction on the hole pockets. We also analyze cos 4\theta angular dependence of the gap in other Fe-pnictides and compare theoretical results with the photoemission experiments of LiFeAs.


title: Swayed by Friends or by the Crowd?

abstract: We have conducted three empirical studies of the effects of friend recommendations and general ratings on how online users make choices. These two components of social influence were investigated through user studies on Mechanical Turk. We find that for a user deciding between two choices an additional rating star has a much larger effect than an additional friend's recommendation on the probability of selecting an item. Equally important, negative opinions from friends are more influential than positive opinions, and people exhibit more random behavior in their choices when the decision involves less cost and risk. Our results can be generalized across different demographics, implying that individuals trade off recommendations from friends and ratings in a similar fashion.


title: Subphotospheric heating in GRBs: analysis and modeling of GRB090902B as observed by Fermi

abstract: We analyze the spectral evolution of GRB 090902B and show that subphotospheric dissipation can explain both the spectra and the spectral evolution. The emission from a GRB photosphere can give rise to a variety of spectral shapes. The spectrum can have a shape close to that of a Planck function (as is observed during the first half of GRB090902B) or be broadened, resembling a typical Band function (as is observed during the second half of GRB090902B). The shape mainly depends on the strength and location of the dissipation in the jet, the ratio of the energy densities of thermal photons and of the electrons at the dissipation site, as well as on the strength of the magnetic field. We further discuss numerical models of the dissipation and relate these to the observed spectra.


title: Optimal Region of Latching Activity in an Adaptive Potts Model for Networks of Neurons

abstract: In statistical mechanics, the Potts model is a model for interacting spins with more than two discrete states. Neural networks which exhibit features of learning and associative memory can also be modeled by a system of Potts spins. A spontaneous behavior of hopping from one discrete attractor state to another (referred to as latching) has been proposed to be associated with higher cognitive functions. Here we propose a model in which both the stochastic dynamics of Potts models and an adaptive potential function are present. A latching dynamics is observed in a limited region of the noise(temperature)-adaptation parameter space. We hence suggest noise as a fundamental factor in such alternations alongside adaptation. From a dynamical systems point of view, the noise-adaptation alternations may be the underlying mechanism for multi-stability in attractor based models. An optimality criterion for realistic models is finally inferred.


title: Out of Medium Fragmentation from Long-Lived Jet Showers

abstract: We study the time structure of vacuum jet evolution via a simple uncertainty principle estimate in the kinematic range explored by current heavy ion collisions at the LHC. We observe that a large fraction of the partonic splittings occur at large times, of the order of several fm. We compare the time distribution of vacuum splittings with the distribution of path lengths traversed by jets in a heavy ion collision. We find that if no medium induced modification of the jet dynamics were present, a very large fraction (larger than 80% for inclusive jets) of the jet splittings would occur outside of the medium. We confront this observation with current available data on jet properties in heavy ion collisions and discuss its implications for the dynamics of jet-medium interactions.


title: Operational Method for Finite Difference Equations

abstract: In this article I present a fast and direct method for solving several types of linear finite difference equations (FDE) with constant coefficients. The method is based on a polynomial form of the translation operator and its inverse, and can be used to find the particular solution of the FDE. This work raises the possibility of developing new ways to expand the scope of the operational methods.


title: X-ray Signatures of Circumnuclear Gas in AGN

abstract: X-ray spectra of AGN are complex. X-ray absorption and emission features trace gas covering a wide range of column densities and ionization states. High resolution spectra show the absorbing gas to be outflowing, perhaps in the form of an accretion disk wind. The absorbing complex shapes the form of the X-ray spectrum while X-ray reverberation and absorption changes explain the spectral and timing behaviour of AGN. We discuss recent progress, highlighting some new results and reviewing the implications that can be drawn from the data.


title: Discovery, classification, and scientific exploration of transient events from the Catalina Real-time Transient Survey

abstract: Exploration of the time domain - variable and transient objects and phenomena - is rapidly becoming a vibrant research frontier, touching on essentially every field of astronomy and astrophysics, from the Solar system to cosmology. Time domain astronomy is being enabled by the advent of the new generation of synoptic sky surveys that cover large areas on the sky repeatedly, and generating massive data streams. Their scientific exploration poses many challenges, driven mainly by the need for a real-time discovery, classification, and follow-up of the interesting events. Here we describe the Catalina Real-Time Transient Survey (CRTS), that discovers and publishes transient events at optical wavelengths in real time, thus benefiting the entire community. We describe some of the scientific results to date, and then focus on the challenges of the automated classification and prioritization of transient events. CRTS represents a scientific and a technological testbed and precursor for the larger surveys in the future, including the Large Synoptic Survey Telescope (LSST) and the Square Kilometer Array (SKA).


title: Heavy quarks in a magnetic field

abstract: The motion of a heavy charged quark in a magnetic field is analyzed in the vacuum of strongly coupled CFT. The motion of the quark is dissipative. It moves in spiral until it eventually comes to rest. The world-sheet geometry is locally AdS_2 but has a time dependent horizon. The string profile in the static gauge extends from the boundary till a point where an embedding singularity exists. Connections with other circular string motions are established.


title: Sources of FCNC in $SU(3)_C \otimes SU(3)_L \otimes U(1)_X$ models

abstract: There are different models based on the gauge symmetry $SU(3)_C \otimes SU(3)_L \otimes U(1)_X$ (331) some of them includes exotic particles and others are constructed without any exotic charges assigned to the fermionic spectrum. Each model build up on 331 symmetry has its own interesting properties according to the representations of the gauge group used for the fermionic spectrum; that is the main reason to explore and identify the possible sources of flavor changing neutral currents and lepton flavor violation at tree level.


title: Irregular Labellings of Circulant Graphs

abstract: We investigate the \textit{irregularity strength} ($s(G)$) and \textit{total vertex irregularity strength} ($tvs(G)$) of circulant graphs $Ci_n(1,2,...,k)$ and prove that $tvs(Ci_n(1,2,...,k))=\lceil\frac{n+2k}{2k+1}\rceil$, while $s(Ci_n(1,2,...,k))=\lceil\frac{n+2k-1}{2k}\rceil$ except the case when $(n \bmod 4k = 2k+1 \wedge k\bmod 2=1) \vee n=2k+1$ and $s(Ci_n(1,2,...,k))=\lceil\frac{n+2k-1}{2k}\rceil+1$.


title: Bayesian Gaussian Copula Factor Models for Mixed Data

abstract: Gaussian factor models have proven widely useful for parsimoniously characterizing dependence in multivariate data. There is a rich literature on their extension to mixed categorical and continuous variables, using latent Gaussian variables or through generalized latent trait models acommodating measurements in the exponential family. However, when generalizing to non-Gaussian measured variables the latent variables typically influence both the dependence structure and the form of the marginal distributions, complicating interpretation and introducing artifacts. To address this problem we propose a novel class of Bayesian Gaussian copula factor models which decouple the latent factors from the marginal distributions. A semiparametric specification for the marginals based on the extended rank likelihood yields straightforward implementation and substantial computational gains, critical for scaling to high-dimensional applications. We provide new theoretical and empirical justifications for using this likelihood in Bayesian inference. We propose new default priors for the factor loadings and develop efficient parameter-expanded Gibbs sampling for posterior computation. The methods are evaluated through simulations and applied to a dataset in political science. The methods in this paper are implemented in the R package bfa.


title: Radio Variability Studies of Gamma-Ray Blazars with the OVRO 40 m Telescope

abstract: Since late 2007, we have been monitoring a large sample of known and likely gamma-ray-loud blazars at 15 GHz twice per week with the Owens Valley Radio Observatory (OVRO) 40 m Telescope. Our initial sample included the 1158 sources above declination -20 degrees from the Candidate Gamma-Ray Blazar Survey (CGRaBS), and we have since added nearly 400 more sources, including all blazars associated with Fermi Large Area Telescope (LAT) detections in the First AGN Catalog (1LAC). Here, we describe the new sample and present results for 2008 through early 2011. Using statistical likelihood analyses, we compare the variability amplitude for various sub-populations within our sample. These include comparisons of gamma-ray-loud versus gamma-ray-quiet objects, BL Lac objects versus flat-spectrum radio quasars, and a study of the variability amplitude trend with redshift. We also describe KuPol, the new digital Ku-band receiver being constructed for the 40 m telescope. This new receiver will provide total intensity and linear polarization measurements over the 12-18 GHz band, with 16 MHz spectral resolution.


title: Algebraic Properties of Codimension Series of PI-Algebras

abstract: For a PI-algebra R over a field of characteristic 0 let T(R) be the T-ideal of the polynomial identities of R and let c(R,t) be the codimension series of R (i.e., the generating function of the codimension sequence of R). Let A, B and R be PI-algebras such that T(R)=T(A)T(B). We show that if c(A,t) and c(B,t) are rational functions, then c(R,t) is also rational. If c(A,t) is rational and c(B,t) is algebraic, then c(R,t) is also algebraic. The proof is based on the fact that the product of two exponential generating functions behaves as the exponential generating function of the sequence of the degrees of the outer tensor products of two sequences of representations of symmetric groups.


title: Constraining Dark Matter Signal from a Combined Analysis of Milky Way Satellites with the Fermi-LAT

abstract: Dwarf spheroidal galaxies have a large mass to light ratio and low astrophysical background, and are therefore considered one of the most promising targets for dark matter searches in the gamma-ray band. By applying a joint likelihood analysis, the power of resultant limits in case of no detection can be enhanced and robust constraints on the dark matter parameter space can be obtained. We present results from a combined analysis of 10 dwarf spheroidal galaxies using Fermi-LAT data. Different annihilation channels have been analyzed and uncertainties from astrophysical properties have been taken into account.


title: Anonymous Meeting in Networks

abstract: A team consisting of an unknown number of mobile agents, starting from different nodes of an unknown network, possibly at different times, have to meet at the same node. Agents are anonymous (identical), execute the same deterministic algorithm and move in synchronous rounds along links of the network. Which configurations are gatherable and how to gather all of them deterministically by the same algorithm? We give a complete solution of this gathering problem in arbitrary networks. We characterize all gatherable configurations and give two universal deterministic gathering algorithms, i.e., algorithms that gather all gatherable configurations. The first algorithm works under the assumption that an upper bound n on the size of the network is known. In this case our algorithm guarantees gathering with detection, i.e., the existence of a round for any gatherable configuration, such that all agents are at the same node and all declare that gathering is accomplished. If no upper bound on the size of the network is known, we show that a universal algorithm for gathering with detection does not exist. Hence, for this harder scenario, we construct a second universal gathering algorithm, which guarantees that, for any gatherable configuration, all agents eventually get to one node and stop, although they cannot tell if gathering is over. The time of the first algorithm is polynomial in the upper bound n on the size of the network, and the time of the second algorithm is polynomial in the (unknown) size itself. Our results have an important consequence for the leader election problem for anonymous agents in arbitrary graphs. For anonymous agents in graphs, leader election turns out to be equivalent to gathering with detection. Hence, as a by-product, we obtain a complete solution of the leader election problem for anonymous agents in arbitrary graphs.


title: Renormalization group invariance in the Pinch Technique

abstract: We show how to construct, using an elementary extension of the Pinch Technique, all off-shell Green's functions of a non-Abelian gauge theory so that they are locally gauge-invariant and renormalization-group invariant (RGI), as the S-matrix is, as well as being process-independent, coupling-constant independent (dimensional transmutation), and satisfying QED-like Ward identities. We call these PT-RGI Green's functions and outline how to construct an approximate three-gluon PT-RGI vertex with three physical scales and no dependence on the renormalization point $\mu$. Properties of the PT-RGI Schwinger-Dyson equations are discussed, mostly in the context of a modified form of $\phi^3_6$. The PT-RGI property of all off-shell Green's functions, plus other work of long ago, leads to a near-realization of the old dreams of S-matrix theorists.


title: A new interpretation of Giant radio pulses from the Crab pulsar

abstract: The Crab pulsar experienced a major flare in 2010 as observed by Fermi LAT. Observations by the Hubble Space Telescope indicate that the flare was accompanied by a structural change in the anvil region of the Crab Nebula. In the framework of a photometric analysis we reconstruct the energetics of this event. Reconnection zones near the light cylinder are expected to release energy by accelerating beams of electrons, leading to flares of varying amplitude. In this case the major flare would have reduced the magnetic energy stored in the reconnection zones, and would thus have had an impact on the properties of the giant radio flares presumably originating from these regions. We test this scenario by observing giant radio pulses with the Westerbork Synthesis Radio Telescope.


title: The joint graphical lasso for inverse covariance estimation across multiple classes

abstract: We consider the problem of estimating multiple related but distinct graphical models on the basis of a high-dimensional data set with observations that belong to distinct classes. A motivating example occurs in the analysis of gene expression data for tissue samples with and without cancer. In this case, we might wish to estimate a gene expression network for the normal tissue and a gene expression network for the tumor tissue. We expect the two gene expression networks to be similar but not identical to each other, and so more accurate estimation of these two networks may be possible using a joint approach. We propose the joint graphical lasso for this purpose. Rather than estimating a graphical model for each class separately, or a single graphical model across all classes, we borrow strength across the classes in order to estimate multiple graphical models that share certain characteristics, such as the locations or weights of nonzero edges. Our approach is based upon maximizing a penalized log likelihood. We employ fused lasso or group lasso penalties, and implement a very fast computational approach that solves the joint graphical lasso problem. In a simulation study we demonstrate that our proposed approach leads to more accurate estimation of networks and covariance structure than competing approaches. We further illustrate our proposal on a publicly-available lung cancer gene expression data set.


title: Pulsar Emission Geometry and Accelerating Field Strength

abstract: The high-quality Fermi LAT observations of gamma-ray pulsars have opened a new window to understanding the generation mechanisms of high-energy emission from these systems. The high statistics allow for careful modeling of the light curve features as well as for phase resolved spectral modeling. We modeled the LAT light curves of the Vela and CTA 1 pulsars with simulated high-energy light curves generated from geometrical representations of the outer gap and slot gap emission models, within the vacuum retarded dipole and force-free fields. A Markov Chain Monte Carlo maximum likelihood method was used to explore the phase space of the magnetic inclination angle, viewing angle, maximum emission radius, and gap width. We also used the measured spectral cutoff energies to estimate the accelerating parallel electric field dependence on radius, under the assumptions that the high-energy emission is dominated by curvature radiation and the geometry (radius of emission and minimum radius of curvature of the magnetic field lines) is determined by the best fitting light curves for each model. We find that light curves from the vacuum field more closely match the observed light curves and multiwavelength constraints, and that the calculated parallel electric field can place additional constraints on the emission geometry.


title: Symmetry Energy from Systematic of Isobaric Analog States

abstract: Excitation energies to isobaric states, that are analogs of ground states, are dominated by contributions from the symmetry energy. This opens up a possibility of investigating the symmetry energy on nucleus-by-nucleus basis. Upon correcting energies of measured nuclear levels for shell and pairing effects, we find that the lowest energies for a given isospin rise in proportion to the square of isospin, allowing for an interpretation of the coefficient of proportionality in terms of a symmetry coefficient for a given nucleus. In the (A,Z) regions where there are enough data, we demonstrate a Z-independence of that coefficient. We further concentrate on the A-dependence of the coefficient, in order to learn about the density dependence of symmetry energy in uniform matter, given the changes of the density in the surface region. In parallel to the analysis of data, we carry out an analysis of the coefficient for nuclei calculated within the Skyrme-Hartree-Fock (SHF) approach, with known symmetry energy for uniform matter. While the data from isobaric analog states suggest a simple interpretation for the A-dependent symmetry coefficient, in terms of the surface and volume symmetry coefficients, the SHF results point to a more complicated situation within the isovector sector than in the isoscalar, with much stronger curvature effects in the first. We exploit the SHF results in estimating the curvature contributions to the symmetry coefficient. That assessment is hampered by instabilities of common Skyrme parameterizations of nuclear interactions.


title: Peculiar anisotropic stationary spherically symmetric solution of Einstein equations

abstract: Motivated by studies on gravitational lenses, we present an exact solution of the field equations of general relativity, which is static and spherically symmetric, has no mass but has a non-vanishing spacelike components of the stress-energy-momentum tensor. In spite of its strange nature, this solution provides with non-trivial descriptions of gravitational effects. We show that the main aspects found in the \emph{dark matter phenomena} can be satisfactorily described by this geometry. We comment on the relevance it could have to consider non-vanishing spacelike components of the stress-energy-momentum tensor ascribed to dark matter.


title: The Average Likelihood Ratio for Large-scale Multiple Testing and Detecting Sparse Mixtures

abstract: Large-scale multiple testing problems require the simultaneous assessment of many p-values. This paper compares several methods to assess the evidence in multiple binomial counts of p-values: the maximum of the binomial counts after standardization (the `higher-criticism statistic'), the maximum of the binomial counts after a log-likelihood ratio transformation (the `Berk-Jones statistic'), and a newly introduced average of the binomial counts after a likelihood ratio transformation. Simulations show that the higher criticism statistic has a superior performance to the Berk-Jones statistic in the case of very sparse alternatives (sparsity coefficient $\beta \gtrapprox 0.75$), while the situation is reversed for $\beta \lessapprox 0.75$. The average likelihood ratio is found to combine the favorable performance of higher criticism in the very sparse case with that of the Berk-Jones statistic in the less sparse case and thus appears to dominate both statistics. Some asymptotic optimality theory is considered but found to set in too slowly to illuminate the above findings, at least for sample sizes up to one million. In contrast, asymptotic approximations to the critical values of the Berk-Jones statistic that have been developed by Wellner and Koltchinskii (2003) and Jager and Wellner (2007) are found to give surprisingly accurate approximations even for quite small sample sizes.


title: Non-classical Solution to Hessian Equation from Cartan Isoparametric Cubic

abstract: We show how to construct a non-smooth solution to Hessian fully nonlinear second-order uniformly elliptic equation using the Cartan isoparametric cubic in 5 dimensions.


title: Exact Sequences in Non-Exact Categories (An Application to Semimodules)

abstract: We consider a notion of exact sequences in any -not necessarily exact- pointed category relative to a given (E;M)-factorization structure. We apply this notion to introduce and investigate a new notion of exact sequences of semimodules over semirings relative to the canonical image factorization. Several homological results are proved using the new notion of exactness including some restricted versions of the Short Five Lemma and the Snake Lemma opening the door for introducing and investigating homology objects in such categories. Our results apply in particular to the variety of commutative monoids extending results in homological varieties to relative homological varieties.


title: Separating enhancement from loss: plasmonic nanocavities in the weak coupling regime

abstract: By modifying the density of optical states at the location of an emitter, weak cavity-emitter coupling can enable a host of potential applications in quantum optics, from the development of low- threshold lasers to brighter single-photon sources for quantum cryptography. Although some of the first demonstrations of spontaneous emission modification occurred in metallic structures, it was only after the recent demonstration of cavity quantum electrodynamics effects in dielectric optical cavities that metal-based optical cavities were considered for quantum optics applications. Advantages of metal-optical cavities include their compatibility with a large variety of emitters and their broadband cavity spectra, which enable enhancement of spectrally-broad emitters. Here, we demonstrate a metal- based nanocavity structure that achieves radiative emission rate enhancements of 1000, opening up the possibility of pursuing cavity electrodynamics investigations with intrinsically broad optical emitters, including organic dyes and colloidal quantum dots.


title: The Kauffman bracket skein module of two-bridge links

abstract: We calculate the Kauffman bracket skein module (KBSM) of the complement of all two-bridge links. For a two-bridge link, we show that the KBSM of its complement is free over the ring $\BC[t^{\pm 1}]$ and when reducing $t=-1$, it is isomorphic to the ring of regular functions on the character variety of the link group.


title: Recent Belle results from Y(5S) sample

abstract: The large data sample recorded with the Belle detector at the Y(5S) energy provides a unique opportunity to study the poorly-known $\bs$ meson. Two analyses, performed with a data sample representing an integrated luminosity of 121 /fb, are presented: the measurement of the Bs0->J/psi f0(980) and Bs0->J/psi f0(1370) branching fractions, and the 5sigma observation of the decay Bs0-bar->Lambda(c)+ pi- Lambda-bar which is the first observation of a baryonic Bs0 decay. In addition, we present new results of a measurement of the CKM angle phi_1(beta) with B pi tagged events.


title: Colour flux-tubes in static Pentaquark and Tetraquark systems

abstract: The colour fields created by the static tetraquark and pentaquark systems are computed in quenched SU(3) lattice QCD, with gauge invariant lattice operators, in a 24^3 x 48 lattice at beta=6.2. We generate our quenched configurations with GPUs, and detail the respective benchmanrks in different SU(N) groups. While at smaller distances the coulomb potential is expected to dominate, at larger distances it is expected that fundamental flux tubes, similar to the flux-tube between a quark and an antiquark, emerge and confine the quarks. In order to minimize the potential the fundamental flux tubes should connect at 120o angles. We compute the square of the colour fields utilizing plaquettes, and locate the static sources with generalized Wilson loops and with APE smearing. The tetraquark system is well described by a double-Y-shaped flux-tube, with two Steiner points, but when quark-antiquark pairs are close enough the two junctions collapse and we have an X-shaped flux-tube, with one Steiner point. The pentaquark system is well described by a three-Y-shaped flux-tube where the three flux the junctions are Steiner points.


title: Search for transient neutrino sources with IceCube

abstract: The IceCube detector, which is embedded in the glacial ice at the geographic South Pole, is the first neutrino telescope to comprise a volume of one cubic kilometer. The search for neutrinos of astrophysical origin is among the primary goals of IceCube. Point source candidates include Galactic objects such as supernova remnants (SNRs) as well as extragalactic objects such as Active Galactic Nuclei (AGN) and Gamma-Ray Bursts (GRBs). Offline and online searches for transient sources like GRBs and supernovae (SNe) are presented. Triggered searches use satellite measurements from Fermi, SWIFT and Konus. Complementary to the triggered offline search, an online neutrino multiplet selection allows IceCube to trigger a network of optical telescopes, which can then identify a possible electromagnetic counterpart. This allows to probe for mildly relativistic jets in SNe and hence to reveal the connection between GRBs, SNe and relativistic jets. Results from IceCube's triggered GRB search and a first limit on relativistic jets in SNe from the optical follow-up program are presented.


title: String tension at finite temperature Lattice QCD

abstract: The critical curve of string tension as a function of the temperature is computed in SU(3) Lattice QCD. We present the results for the string tension utilizing a pair of Polyakov loop and antiloop, with two different techniques. We compare the colour averaged loop-antiloop which is gauge invariant but is only adequate to study the string tension, and the colour singlet loop-antiloop using the Landau gauge fixing which also enables to compute the coulomb part of the free energy.


title: OpenWeather: a peer-to-peer weather data transmission protocol

abstract: The study of the weather is performed using instruments termed weather stations. These weather stations are distributed around the world, collecting the data from the different phenomena. Several weather organizations have been deploying thousands of these instruments, creating big networks to collect weather data. These instruments are collecting the weather data and delivering it for later processing in the collections points. Nevertheless, all the methodologies used to transmit the weather data are based in protocols non adapted for this purpose. Thus, the weather stations are limited by the data formats and protocols used in them, not taking advantage of the real-time data available on them. We research the weather instruments, their technology and their network capabilities, in order to provide a solution for the mentioned problem. OpenWeather is the protocol proposed to provide a more optimum and reliable way to transmit the weather data. We evaluate the environmental factors, such as location or bandwidth availability, in order to design a protocol adapted to the requirements established by the automatic weather stations. A peer to peer architecture is proposed, providing a functional implementation of OpenWeather protocol. The evaluation of the protocol is executed in a real scenario, providing the hints to adapt the protocol to a common automatic weather station.


title: Exploring Linkablility of Community Reviewing

abstract: Large numbers of people all over the world read and contribute to various review sites. Many contributors are understandably concerned about privacy; specifically, about linkability of reviews (and accounts) across review sites. In this paper, we study linkability of community reviewing and try to answer the question: to what extent are "anonymous" reviews linkable, i.e., likely authored by the same contributor? Based on a very large set of reviews from a popular site (Yelp), we show that a high percentage of ostensibly anonymous reviews can be linked with very high confidence. This is despite the fact that we use very simple models and equally simple features set. Our study suggests that contributors reliably expose their identities in reviews. This has important implications for cross-referencing accounts between different review sites. Also, techniques used in our study could be adopted by review sites to give contributors feedback about privacy of their reviews.


title: Hybridization of wave functions in one-dimensional localization

abstract: A quantum particle can be localized in a disordered potential, the effect known as Anderson localization. In such a system, correlations of wave functions at very close energies may be described, due to Mott, in terms of a hybridization of localized states. We revisit this hybridization description and show that it may be used to obtain quantitatively exact expressions for some asymptotic features of correlation functions, if the tails of the wave functions and the hybridization matrix elements are assumed to have log-normal distributions typical for localization effects. Specifically, we consider three types of one-dimensional systems: a strictly one-dimensional wire and two quasi-one-dimensional wires with unitary and orthogonal symmetries. In each of these models, we consider two types of correlation functions: the correlations of the density of states at close energies and the dynamic response function at low frequencies. For each of those correlation functions, within our method, we calculate three asymptotic features: the behavior at the logarithmically large "Mott length scale", the low-frequency limit at length scale between the localization length and the Mott length scale, and the leading correction in frequency to this limit. In the several cases, where exact results are available, our method reproduces them within the precision of the orders in frequency considered.


title: The Incompressible Euler Limit of the Boltzmann Equation with Accommodation Boundary Condition

abstract: The convergence of solutions of the incompressible Navier-Stokes equations set in a domain with boundary to solutions of the Euler equations in the large Reynolds number limit is a challenging open problem both in 2 and 3 space dimensions. In particular it is distinct from the question of existence in the large of a smooth solution of the initial-boundary value problem for the Euler equations. The present paper proposes three results in that direction. First, if the solutions of the Navier-Stokes equations satisfy a slip boundary condition with vanishing slip coefficient in the large Reynolds number limit, we show by an energy method that they converge to the classical solution of the Euler equations on its time interval of existence. Next we show that the incompressible Navier-Stokes limit of the Boltzmann equation with Maxwell's accommodation condition at the boundary is governed by the Navier-Stokes equations with slip boundary condition, and we express the slip coefficient at the fluid level in terms of the accommodation parameter at the kinetic level. This second result is formal, in the style of [Bardos-Golse-Levermore, J. Stat. Phys. 63 (1991), 323-344]. Finally, we establish the incompressible Euler limit of the Boltzmann equation set in a domain with boundary with Maxwell's accommodation condition assuming that the accommodation parameter is small enough in terms of the Knudsen number. Our proof uses the relative entropy method following closely the analysis in [L. Saint-Raymond, Arch. Ration. Mech. Anal. 166 (2003), 47-80] in the case of the 3-torus, except for the boundary terms, which require special treatment.


title: Optical Spectroscopic Survey of High Latitude WISE Selected Sources

abstract: We report on the results of an optical spectroscopic survey at high Galactic latitude (b\geq30{\deg}) of a sample of WISE-selected targets, grouped by WISE W1 ({\lambda}_eff = 3.4 \mum) flux, which we use to characterize the sources WISE detected. We observed 762 targets in 10 disjoint fields centered on ultra-luminous infrared galaxy (ULIRG) candidates using the DEIMOS spectrograph on Keck II. We find 0.30 \pm 0.02 galaxies arcmin^-2 with a median redshift of z=0.33 \pm 0.01 for the sample with W1 \geq 120 \muJy. The foreground stellar densities in our survey range from 0.23 \pm 0.07 arcmin^(-2) to 1.1 \pm 0.1 arcmin^(-2) for the same sample. We obtained spectra that produced science grade redshifts for \geq 90% of our targets for sources with W1 flux \geq 120 \muJy that also had i-band flux \geq 18 \muJy. We used for targeting very preliminary data reductions available to the team in August of 2010. Our results therefore present a conservative estimate of what is possible to achieve using WISE's Preliminary Data Release for the study of field galaxies.


title: Colour flux-tubes in static Pentaquark systems

abstract: The colour fields created by the static tetraquark and pentaquark systems are computed in quenched SU(3) lattice QCD, with gauge invariant lattice operators, in a 24^3 x 48 lattice at beta=6.2 . We generate our quenched configurations with GPUs, and detail the respective benchmanrks in different SU(N) groups. While at smaller distances the coulomb potential is expected to dominate, at larger distances it is expected that fundamental flux tubes, similar to the flux-tube between a quark and an antiquark, emerge and confine the quarks. In order to minimize the potential the fundamental flux tubes should connect at 120o angles. We compute the square of the colour fields utilizing plaquettes, and locate the static sources with generalized Wilson loops and with APE smearing. The tetraquark system is well described by a double-Y-shaped flux-tube, with two Steiner points, but when quark-antiquark pairs are close enough the two junctions collapse and we have an X-shaped flux-tube, with one Steiner point. The pentaquark system is well described by a three-Y-shaped flux-tube where the three flux the junctions are Steiner points.


title: Application of metasurface description for multilayered metamaterials and an alternative theory for metamaterial perfect absorber

abstract: We analyze single and multilayered metamaterials by modeling each layer as a metasurface with effective surface electric and magnetic susceptibility derived through a thin film approximation. Employing a transfer matrix method, these metasurfaces can be assembled into multilayered metamaterials to realize certain functionalities. We demonstrate numerically that this approach provides an alternative interpretation of metamaterial-based perfect absorption, showing that the underlying mechanism is a modified Fabry-Perot resonance. This method provides a general approach applicable for decoupled or weakly coupled multilayered metamaterials.


title: From the Boltzmann Equation to the Euler Equations in the Presence of Boundaries

abstract: The fluid dynamic limit of the Boltzmann equation leading to the Euler equations for an incompressible fluid with constant density in the presence of material boundaries shares some important features with the better known inviscid limit of the Navier-Stokes equations. The present paper slightly extends recent results from [C. Bardos, F. Golse, L. Paillard, Comm. Math. Sci., 10 (2012), 159--190] to the case of boundary conditions for the Boltzmann equation more general than Maxwell's accomodation condition.


title: Searching for Unmodeled Sources Using the Earth Occultation Data from the Fermi GBM

abstract: Employing the 12 NaI detectors in the Fermi GBM, the Earth Occultation Technique (EOT) can be used to measure the fluxes of x-ray and gamma-ray sources. Each time a source passes behind the Earth (or emerges from behind the Earth), a step-like feature is produced in the detector count rate. With a predefined catalog of source positions, the times of the occultation steps can be calculated, the individual steps fit, and the fluxes derived. However, in order to find new sources and generate a complete catalog, a method is needed for generating an image of the sky. An imaging algorithm has been developed to generate all-sky images using the GBM data. Here we present imaging results from ~2.5 years of data in the 12-25 keV and 100-300 keV energy bands.


title: The Photometric Periods of the Nova-Like Cataclysmic Variable LQ Pegasi (PG 2133+115)

abstract: We present a time-resolved differential photometric study and time series analysis of the nova-like cataclysmic variable star LQ Peg. We discover three periodicities in the photometry, one with a period of 3.42 +/- 0.03 hours, and another with a period of 56.8 +/- 0.01 hours. We interpret these to be the apsidal superhump and precessional periods of the accretion disk, respectively, and predict that the orbital period of LQ Peg is 3.22 +/- 0.03 hours. The third periodicity, with a period of 41.3 +/- 0.01 hours, we interpret to be the nodal precessional period of the accretion disk. We also report a flare that lasted four minutes and had an energy in visible light of (1.2 +/- 0.3) x 10^36 ergs, or 10^4-5 times more energetic than the largest solar flares, comparable to the most energetic visible-light stellar flares known. We calculate the absolute magnitude of LQ Peg to be M_J = 4.78 +/- 0.54, and its distance to be 800 +/- 200 pc.


title: Constraints on neutron star mass and radius in GS 1826-24 from sub-Eddington X-ray bursts

abstract: We investigate the constraints on neutron star mass and radius in GS 1826-24 from models of lightcurves and spectral evolution of type I X-ray bursts. This source shows remarkable agreement with theoretical calculations of burst energies, recurrence times, and lightcurves. We first exploit this agreement to set the overall luminosity scale of the observed bursts. When combined with a measured blackbody normalization, this leads to a distance and anisotropy independent measurement of the ratio between the redshift 1+z and color correction factor f_c. We find 1+z=1.19-1.28 for f_c=1.4-1.5. We then compare the evolution of the blackbody normalization with flux in the cooling tail of bursts with predictions from spectral models of Suleimanov et al. (2011b). The observations are well described by the models at luminosities greater than about one third of the peak luminosity, with deviations emerging at luminosities below that. We show that this comparison leads to distance independent upper limits on R_\infty and neutron star mass of R_\infty\lesssim 9.0-13.2 km and M<1.2-1.7 M_\odot, respectively, for solar abundance of hydrogen at the photosphere and a range of metallicity and surface gravity. The radius limits are low in comparison to previous measurements. This may be indicative of a subsolar hydrogen fraction in the GS 1826-24 photosphere, or of larger color corrections than that predicted by spetral models. Our analysis also gives an upper limit on the distance to GS 1826-24 of d<4.0-5.5 kpc \xi_b^{-1/2}, where \xi_b is the degree of anisotropy of the burst emission.


title: The Photometric Period of the Cataclysmic Variable HV Andromedae

abstract: We present four nights of time-resolved photometry of the cataclysmic variable star HV And. Our time series analysis has revealed a prominent period at 3.368 +/- 0.060 hours, as well as some low frequency power. We interpret this signal, from saw-tooth waves in the light curve, as evidence of superhumps in HV And.


title: A penalized likelihood approach to estimate within-household contact networks from egocentric data

abstract: Acute infectious diseases are transmitted over networks of social contacts. Epidemic models are used to predict the spread of emergent pathogens and compare intervention strategies. Many of these models assume equal probability of contact within mixing groups (homes, schools, etc.), but little work has inferred the actual contact network, which may influence epidemic estimates. We develop a penalized likelihood method to infer contact networks within households, a key area for disease transmission. Using egocentric surveys of contact behavior in Belgium, we estimate within-household contact networks for six different age compositions. Our estimates show dependency in contact behavior and vary substantively by age composition, with fewer contacts occurring in older households. Our results are relevant for epidemic models used to make policy recommendations.


title: Results from the binaries LS I +61{\deg}303 and LS 5039 after 2.5 years of Fermi monitoring

abstract: The Fermi Large Area Telescope (LAT) has made the first definitive GeV detections of the binaries LS I +61^{\circ}303 and LS5039 in the first year after its launch in August 2008. These detections were unambiguous because, apart from a reduced positional uncertainty, the gamma-ray emission in each case was orbitally modulated with the corresponding orbital period. The LAT results posed new questions about the nature of these objects, after the unexpected observation of an exponential cutoff in the GeV gamma-ray spectra of both LS I +61^{\circ}303 and LS5039, at least along part of their orbital motion. We present here the analysis of new data from the LAT, comprising 2.5 years of observations through which LS I +61^{\circ}303 continues to provide some surprises. We find an increase in flux in March 2009 and a steady decrease in the flux fraction modulation. The LAT now detects emission up to 30 GeV, where prior datasets led to upper limits only. At the same time, contemporaneous TeV observations either no longer detected the source, or found it -at least in some orbits- close to periastron, far from the usual phases in which the source usually appeared at TeV energies. The on-source exposure of LS 5039 has also drastically increased along the last years, and whilst our analysis shows no new behavior in comparison with our earlier report, the higher statistics of the current dataset allows for a deeper investigation of its orbital and spectral evolution.


title: A simple mathematical model for assessment of anti-toxin antibodies

abstract: The toxins associated with infectious diseases are potential targets for inhibitors which have the potential for prophylactic or therapeutic use. Many antibodies have been generated for this purpose, and the objective of this study was to develop a simple mathematical model that may be used to evaluate the potential protective effect of antibodies. This model was used to evaluate the contributions of antibody affinity and concentration to reducing antibody-receptor complex formation and internalization. The model also enables prediction of the antibody kinetic constants and concentration required to provide a specified degree of protection. We hope that this model, once validated experimentally, will be a useful tool for in vitro selection of potentially protective antibodies for progression to in vivo evaluation.


title: Revisiting k-means: New Algorithms via Bayesian Nonparametrics

abstract: Bayesian models offer great flexibility for clustering applications---Bayesian nonparametrics can be used for modeling infinite mixtures, and hierarchical Bayesian models can be utilized for sharing clusters across multiple data sets. For the most part, such flexibility is lacking in classical clustering methods such as k-means. In this paper, we revisit the k-means clustering algorithm from a Bayesian nonparametric viewpoint. Inspired by the asymptotic connection between k-means and mixtures of Gaussians, we show that a Gibbs sampling algorithm for the Dirichlet process mixture approaches a hard clustering algorithm in the limit, and further that the resulting algorithm monotonically minimizes an elegant underlying k-means-like clustering objective that includes a penalty for the number of clusters. We generalize this analysis to the case of clustering multiple data sets through a similar asymptotic argument with the hierarchical Dirichlet process. We also discuss further extensions that highlight the benefits of our analysis: i) a spectral relaxation involving thresholded eigenvectors, and ii) a normalized cut graph clustering algorithm that does not fix the number of clusters in the graph.