title: Static magnetic moments revealed by muon spin relaxation and thermodynamic measurements in quantum spin ice Yb$_2$Ti$_2$O$_7$

abstract: We present muon spin relaxation ($\mu$SR) and specific-heat versus temperature $C(T)$ measurements on polycrystalline and single-crystal samples of the pyrochlore magnet Yb$_2$Ti$_2$O$_7$. $C(T)$ exhibits a sharp peak at a $T_\mathrm{C}$ of 0.21 and 0.26~K for the single-crystal and polycrystalline samples respectively. For both samples, the magnetic entropy released between 50~mK and 30~K amounts to $R\ln2$ per Yb. At temperatures below $T_\mathrm{C}$ we observe a steep drop in the asymmetry of the zero-field $\mu$SR time spectra at short time scales, as well as a decoupling of the muon spins from the internal field in longitudinal magnetic fields of $\leq0.25$~T for both the polycrystalline and single-crystal samples. These muon data are indicative of static magnetic moments. Our results are consistent with the onset of long-range magnetic order in both forms of Yb$_2$Ti$_2$O$_7$.


title: ITLinQ: A New Approach for Spectrum Sharing in Device-to-Device Communication Systems

abstract: We consider the problem of spectrum sharing in device-to-device communication systems. Inspired by the recent optimality condition for treating interference as noise, we define a new concept of "information-theoretic independent sets" (ITIS), which indicates the sets of links for which simultaneous communication and treating the interference from each other as noise is information-theoretically optimal (to within a constant gap). Based on this concept, we develop a new spectrum sharing mechanism, called "information-theoretic link scheduling" (ITLinQ), which at each time schedules those links that form an ITIS. We first provide a performance guarantee for ITLinQ by characterizing the fraction of the capacity region that it can achieve in a network with sources and destinations located randomly within a fixed area. Furthermore, we demonstrate how ITLinQ can be implemented in a distributed manner, using an initial 2-phase signaling mechanism which provides the required channel state information at all the links. Through numerical analysis, we show that distributed ITLinQ can outperform similar state-of-the-art spectrum sharing mechanisms, such as FlashLinQ, by more than a 100% of sum-rate gain, while keeping the complexity at the same level. Finally, we discuss a variation of the distributed ITLinQ scheme which can also guarantee fairness among the links in the network and numerically evaluate its performance.


title: Electron Transport Properties of Composite Ferroelectrics

abstract: We study electron transport in composite ferroelectrics --- materials consisting of metallic grains embedded in a ferroelectric matrix. Due to its complex tunable morphology the thermodynamic properties of these materials can be essentially different from bulk or thin-film ferroelectrics. We calculate the conductivity of composite ferroelectrics by taking into account the interplay between charge localization, multiple grain boundaries, strong Coulomb repulsion, and ferroelectric order parameter. We show that the ferroelectricity plays a crucial role on the temperature behavior of the conductivity in the vicinity of the ferroelectric-paraelectric transition.


title: Intrinsic graphene field effect transistor on amorphous carbon films

abstract: Fabrication of graphene field effect transistor is described which uses an intrinsic graphene on the surface of as deposited hydrogenated amorphous carbon films. Ambipolar characteristic has been demonstrated typical for graphene devices, which changes to unipolar characteristic if the surface graphene was etched in oxygen plasma. Because amorphous carbon films can be growth easily, with unlimited dimensions and no transfer of graphene is necessary, this can open new perspective for graphene electronics.


title: Entanglement genesis by ancilla-based parity measurement in 2D circuit QED

abstract: We present an indirect two-qubit parity meter in planar circuit quantum electrodynamics, realized by discrete interaction with an ancilla and a subsequent projective ancilla measurement with a dedicated, dispersively coupled resonator. Quantum process tomography and successful entanglement by measurement demonstrate that the meter is intrinsically quantum non-demolition. Separate interaction and measurement steps allow commencing subsequent data qubit operations in parallel with ancilla measurement, offering time savings over continuous schemes.


title: A Balloon-borne Measurement of High Latitude Atmospheric Neutrons Using a LiCAF Neutron Detector

abstract: PoGOLino is a scintillator-based neutron detector. Its main purpose is to provide data on the neutron flux in the upper stratosphere at high latitudes at thermal and nonthermal energies for the PoGOLite instrument. PoGOLite is a balloon borne hard X-ray polarimeter for which the main source of background stems from high energy neutrons. No measurements of the neutron environment for the planned flight latitude and altitude exist. Furthermore this neutron environment changes with altitude, latitude and solar activity, three variables that will vary throughout the PoGOLite flight. PoGOLino was developed to study the neutron environment and the influences from these three variables upon it. PoGOLino consists of two Europium doped Lithium Calcium Aluminium Fluoride (Eu:LiCAF) scintillators, each of which is sandwiched between 2 Bismuth Germanium Oxide (BGO) scintillating crystals, which serve to veto signals produced by gamma-rays and charged particles. This allows the neutron flux to be measured even in high radiation environments. Measurements of neutrons in two separate energy bands are achieved by placing one LiCAF detector inside a moderating polyethylene shield while the second detector remains unshielded. The PoGOLino instrument was launched on March 20th 2013 from the Esrange Space Center in Northern Sweden to an altitude of 30.9 km. A description of the detector design and read-out system is presented. A detailed set of simulations of the atmospheric neutron environment performed using both PLANETOCOSMICS and Geant4 will also be described. The comparison of the neutron flux measured during flight to predictions based on these simulations will be presented and the consequences for the PoGOLite background will be discussed.


title: Quantum corrections to scalar field dynamics in a slow-roll space-time

abstract: We consider the dynamics of a quantum scalar field in the background of a slow-roll inflating Universe. We compute the one-loop quantum corrections to the field and Friedmann equation of motion, in both a 1PI and a 2PI expansion, to leading order in slow-roll. Generalizing the works of [1-3], we then solve these equations to compute the effect on the primordial power spectrum, for the case of a self-interacting inflaton and a self-interacting spectator field. We find that for the inflaton the corrections are negligible due to the smallness of the coupling constant despite the large IR enhancement of the loop contributions. For a curvaton scenario, on the other hand, we find tension in using the 1PI loop corrections, which may indicate that the quantum corrections could be non-perturbatively large in this case, thus requiring resummation.


title: The Steklov spectrum of surfaces: asymptotics and invariants

abstract: We obtain precise asymptotics for the Steklov eigenvalues on a compact Riemannian surface with boundary. It is shown that the number of connected components of the boundary, as well as their lengths, are invariants of the Steklov spectrum. The proofs are based on pseudodifferential techniques for the Dirichlet-to-Neumann operator and on a number-theoretic argument.


title: Extended Lyman-alpha emission from a high-z DLA at z=3.115

abstract: We searched for star formation activity associated with high-z Damped Lyman-alpha systems (DLAs) with Subaru telescope. We used a set of narrow-band (NB) filters whose central wavelengths correspond to the redshifted Lyman-alpha emission lines of targeted DLA absorbers at 3<z<4.5. We detected one apparent NB-excess object located 3.80 arcsec (~28kpc) away from the quasar SDSS J031036.84+005521.7. Follow-up spectroscopy revealed an asymmetric Lyman-alpha emission at z_em=3.115+/-0.003, which perfectly matches the sub-DLA trough at z_abs=3.1150 with logN(HI)/cm^-2=20.05. The Lyman-alpha luminosity is estimated to be L(LyA)=1.07x10^42 erg s^-1, which corresponds to a star formation rate of 0.97 M_\odot yr^-1. Interestingly, the detected Lyman-alpha emission is spatially extended with a sharp peak. The large extent of the Lyman-alpha emission is remarkably one-sided toward the quasar line-of-sight, and is redshifted. The observed spatially asymmetric surface brightness profile can be qualitatively explained by a model of a DLA host galaxy, assuming a galactic outflow and a clumpy distribution of HI clouds in the circumgalactic medium. This large Lyman-alpha extension, which is similar to those found in Rauch et al. (2008), could be the result of complicated anisotropic radiative transfer through the surrounding neutral gas embedded in the DLA.


title: The symplectic arc algebra is formal

abstract: We prove a formality theorem for the Fukaya categories of the symplectic manifolds underlying symplectic Khovanov cohomology, over fields of characteristic zero. The key ingredient is the construction of a degree one Hochschild cohomology class on a Floer A-infinity algebra associated to the (k,k)-nilpotent slice Y, obtained by counting holomorphic discs which satisfy a suitable conormal condition at infinity in a partial compactification of Y. The partial compactification is obtained as the Hilbert scheme of a partial compactification of a Milnor fibre. A sequel to this paper will prove formality of the symplectic cup and cap bimodules, and infer that symplectic Khovanov cohomology and Khovanov cohomology have the same total rank over characteristic zero fields.


title: Superdescendants of the D1D5 CFT and their dual 3-charge geometries

abstract: We describe how to obtain the gravity duals of semiclassical states in the D1-D5 CFT that are superdescendants of a class of RR ground states. On the gravity side, the configurations we construct are regular and asymptotically reproduce the 3-charge D1-D5-P black hole compactified on $S^1\times T^4$. The geometries depend trivially on the $T^4$ directions but non-trivially on the remaining 6D space. In the decoupling limit, they reduce to asymptotically AdS$_3 \times S^3 \times T^4$ spaces that are dual to CFT states obtained by acting with (exponentials of) the operators of the superconformal algebra. As explicit examples, we generalise the solution first constructed in arXiv:1306.1745 and discuss another class of states that have a more complicated dual geometry. By using the free orbifold description of the CFT we calculate the average values for momentum and the angular momenta of these configurations. Finally we compare the CFT results with those obtained in the bulk from the asymptotically $M^{1,4} \times S^1\times T^4$ region.


title: CKM benchmarks for electron EDM experiments

abstract: All current experiments searching for an electron EDM d_e are performed with atoms and diatomic molecules. Motivated by significant recent progress in searches for an EDM-type signal in diatomic molecules with an uncompensated electron spin, we provide an estimate for the expected signal in the Standard Model due to the CKM phase. We find that the main contribution originates from the effective electron-nucleon operator $\bar{e} i\gamma_5 e \bar{N}N$, induced by a combination of weak and electromagnetic interactions at $O(G_F^2\alpha^2)$, and not by the CKM-induced electron EDM itself. When the resulting atomic P,T-odd mixing is interpreted as an {\it equivalent} electron EDM, this estimate leads to the benchmark $d_e^{equiv}(CKM) \sim 10^{-38}$ ecm.


title: An approach to intersection theory on singular varieties using motivic complexes

abstract: We introduce techniques of Suslin, Voevodsky, and others into the study of singular varieties. Our approach is modeled after Goresky-MacPherson intersection homology. We provide a formulation of perversity cycle spaces leading to perversity homology theory and a companion perversity cohomology theory based upon generalized cocycle spaces. These theories lead to conditions on pairs of cycles which can be intersected and a suitable equivalence relation on cocycles/cycles enabling pairings on equivalence classes. We establish suspension and splitting theorems, as well as a localization property. Some examples of intersections on singular varieties are computed.


title: Universal symmetry-protected topological invariants for symmetry-protected topological states

abstract: Symmetry-protected topological (SPT) states are short-range entangled states with a symmetry G. They belong to a new class of quantum states of matter which are classified by the group cohomology $H^{d+1}(G,\mathbb{R}/\mathbb{Z})$ in d-dimensional space. In this paper, we propose a class of symmetry- protected topological invariants that may allow us to fully characterize SPT states with a symmetry group G (ie allow us to measure the cocycles in $H^{d+1}(G,\mathbb{R}/\mathbb{Z})$ that characterize the SPT states). We give an explicit and detailed construction of symmetry-protected topological invariants for 2+1D SPT states. Such a construction can be directly generalized to other dimensions.


title: Periodic solutions of semi-explicit differential-algebraic equations with time-dependent constraints

abstract: In this paper we investigate the properties of the set of T-periodic solutions of semi-explicit parametrized Differential-Algebraic Equations with non-autonomous constraints of a particular type. We provide simple, degree theoretic conditions for the existence of branches of T-periodic solutions of the considered equations. Our approach is based on topological arguments about differential equations on implicitly defined manifolds, combined with elementary facts of matrix analysis.


title: Water Adsorption at the Tetrahedral Titania Surface Layer of SrTiO$_3$(110)-(4$\times$1)

abstract: The interaction of water with oxide surfaces is of great interest for both fundamental science and applications. We present a combined theoretical [density functional theory (DFT)] and experimental [Scanning Tunneling Microscopy (STM), photoemission spectroscopy (PES)] study of water interaction with the two-dimensional titania overlayer that terminates the SrTiO$_3$(110)-(4$\times$1) surface and consists of TiO$_4$ tetrahedra. STM, core-level and valence band PES show that H$_2$O neither adsorbs nor dissociates on the stoichiometric surface at room temperature, while it dissociates at oxygen vacancies. This is in agreement with DFT calculations, which show that the energy barriers for water dissociation on the stoichiometric and reduced surfaces are 1.7 and 0.9 eV, respectively. We propose that water weakly adsorbs on two-dimensional, tetrahedrally coordinated overlayers.


title: QUADRO: A supervised dimension reduction method via Rayleigh quotient optimization

abstract: We propose a novel Rayleigh quotient based sparse quadratic dimension reduction method - named QUADRO (Quadratic Dimension Reduction via Rayleigh Optimization) - for analyzing high- dimensional data. Unlike in the linear setting where Rayleigh quotient optimization coincides with classification, these two problems are very different under nonlinear settings. In this paper, we clarify this difference and show that Rayleigh quotient optimization may be of independent scientific interests. One major challenge of Rayleigh quotient optimization is that the variance of quadratic statistics involves all fourth cross-moments of predictors, which are infeasible to compute for high-dimensional applications and may accumulate too many stochastic errors. This issue is resolved by considering a family of elliptical models. Moreover, for heavy-tail distributions, robust estimates of mean vectors and covariance matrices are employed to guarantee uniform convergence in estimating nonpolynomially many parameters, even though only the fourth moments are assumed. Methodologically, QUADRO is based on elliptical models which allow us to formulate the Rayleigh quotient maximization as a convex optimization problem. Computationally, we propose an efficient linearized augmented Lagrangian method to solve the constrained optimization problem. Theoretically, we provide explicit rates of convergence in terms of Rayleigh quotient under both Gaussian and general elliptical models. Thorough numerical results on both synthetic and real datasets are also provided to back up our theoretical results.


title: What is the probability that a random integral quadratic form in $n$ variables is isotropic?

abstract: We show that the density of quadratic forms in $n$ variables over ${\mathbb Z}_p$ that are isotropic is a rational function in $p$, where the rational function is independent of $p$, and we determine this rational function explicitly. As a consequence, for each $n$, we determine the probability that a random integral quadratic form in $n$ variables is isotropic. In particular, we show that the probability that a random integral quaternary quadratic form is isotropic is $\approx 97.0\%$, in the case where the coefficients of the quadratic form are independently and uniformly distributed in the range $[-X,X]$ with $X\to\infty$. When random integral quaternary quadratic forms are chosen with respect to the Gaussian Orthogonal Ensemble (GOE), the probability of isotropy increases to $\approx 98.3\%$.


title: Experimental investigation of the trachea oscillation and its role in the pitch production

abstract: Several experiments have been performed to investigate the mechanical vibrations associated with trachea and larynx when Italian vowels are emitted. The mechanical measurements have been made by using two laser Doppler vibrometers (based on the well-known not-invasive optical measurement technique) coupled with the acoustic field measured by high-quality certified microphones. The recorded signals are analyzed by using well-established methods in time and frequency domains. The signals of the mechanical vibrations along the trachea and the larynx are compared with those of the acoustic ones. Focusing the attention of the signals' onsets, we can observe an upward propagation of the mechanical vibrations for which it is possible to estimate a delay between the traces. We observe that the mechanical oscillations at the trachea start before the larynx and acoustic oscillations. Moreover these tracheal oscillations are self-oscillations in time and are associated with the pitch production, indicating a further hydrodynamic instability at trachea. This leads to new insights in the mechanism controlling the pitch in the speech.


title: A vanishing theorem for log canonical pairs after De Fernex-Ein

abstract: We extend the main vanishing theorem in a paper of de Fernex and Ein to singular varieties without assuming locally complete intersection.


title: A Mathematica Package for Calculation of One-Loop Penguins in FCNC Processes

abstract: In this work, we present a Mathematica package Peng4BSM@LO which calculates the contributions to the Wilson Coefficients of certain effective operators originating from the one-loop penguin Feynman diagrams. Both vector and scalar external legs are considered. The key feature of our package is the ability to find the corresponding expressions in almost any New Physics model which extends the SM and has no flavour changing neutral current (FCNC) transitions at the tree level.


title: Long division unites - long union divides, a model for social network evolution

abstract: A remarkable phenomenon in the time evolution of many networks such as cultural, political, national and economic systems, is the recurrent transition between the states of union and division of nodes. In this work, we propose a phenomenological modeling, inspired by the maxim "long union divides and long division unites", in order to investigate the evolutionary characters of these networks composed of the entities whose behaviors are dominated by these two events. The nodes are endowed with quantities such as identity, ingredient, richness (power), openness (connections), age, distance, interaction etc. which determine collectively the evolution in a probabilistic way. Depending on a tunable parameter, the time evolution of this model is mainly an alternative domination of union or division state, with a possible state of final union dominated by one single node.


title: Simplices in the Endomorphism Semiring of a Finite Chain

abstract: We establish new results concerning endomorphisms of a finite chain if the cardinality of the image of such endomorphism is no more than some fixed number. The semiring of all such endomorphisms can be seen as a simplex whose vertices are the constant endomorphisms. We explore the properties of these simplices.


title: On q-algebraic equations and their power series solutions

abstract: We study the existence of formal power series solutions to q-algebraic equations. When a solution exists, we give a sufficient condition on the equation for this solution to have a positive radius of convergence. We emphasize on the case where the solution is divergent, giving a sharp estimate on the growth of the coefficients. As a consequence, we obtain a bound on the q-Gevrey order of the formal solution, which is optimal in some cases. Various examples illustrate our main results.


title: Composable Languages for Bioinformatics: The NYoSh experiment

abstract: Language workbenches are software engineering tools that help domain experts develop solutions to various classes of problems. Some of these tools focus on non-technical users and provide languages to help organize knowledge while other workbenches provide means to create new programming languages. A key advantage of language workbenches is that they support the composition of independently developed languages. This capability is useful when developing programs that can benefit from different levels of abstraction. We reasoned that language workbenches could be useful to develop bioinformatics software solutions. In order to evaluate the potential of language workbenches in bioinformatics, we tested a prominent workbench by developing an alternative to shell scripting. While shell scripts are widely used in bioinformatics to automate computational analysis, existing scripting languages do not provide many of the features present in modern programming languages. We report on our design of NYoSh (Not Your ordinary Shell). NYoSh was implemented as a collection of languages that can be composed to write programs as expressive and concise as shell scripts. NYoSh offers a concrete illustration of the advantages that language workbench technologies can bring to bioinformatics. For instance, NYoSh scripts can be edited with an environment-aware editor that provides semantic error detection and can be compiled interactively with an automatic build and deployment system. In contrast to shell scripts, NYoSh scripts can be written in a modern development environment, supporting context dependent intentions and can be extended seamlessly with new abstractions and language constructs. We demonstrate language extension and composition by presenting a tight integration of NYoSh scripts with the GobyWeb system. The NYoSh Workbench prototype is distributed at http://nyosh.campagnelab.org


title: Modelling Gamma Ray Emission and Pair Production in High-Intensity Laser-Matter Interactions

abstract: In high-intensity (> 10^21W/cm^2) laser-matter interactions gamma-ray photon emission by the electrons can strongly affect the electron's dynamics and copious numbers of electron-positron pairs can be produced by the emitted photons. We show how these processes can be included in simulations by coupling a Monte-Carlo algorithm describing the emission to a particle-in-cell code. The Monte-Carlo algorithm includes quantum corrections to the photon emission, which we show must be included if the pair production rate is to be correctly determined. The accuracy, convergence and energy conservation properties of the Monte-Carlo algorithm are analysed in simple test problems.


title: Bayesian Discovery of Threat Networks

abstract: A novel unified Bayesian framework for network detection is developed, under which a detection algorithm is derived based on random walks on graphs. The algorithm detects threat networks using partial observations of their activity, and is proved to be optimum in the Neyman-Pearson sense. The algorithm is defined by a graph, at least one observation, and a diffusion model for threat. A link to well-known spectral detection methods is provided, and the equivalence of the random walk and harmonic solutions to the Bayesian formulation is proven. A general diffusion model is introduced that utilizes spatio-temporal relationships between vertices, and is used for a specific space-time formulation that leads to significant performance improvements on coordinated covert networks. This performance is demonstrated using a new hybrid mixed-membership blockmodel introduced to simulate random covert networks with realistic properties.

categories: cs.SI cs.LG math.ST physics.soc-ph stat.ML stat.TH


title: Random Matrix Spectra as a Time Series

abstract: Spectra of ordered eigenvalues of finite Random Matrices are interpreted as a time series. Dataadaptive techniques from signal analysis are applied to decompose the spectrum in clearly differentiated trend and fluctuation modes, avoiding possible artifacts introduced by standard unfolding techniques. The fluctuation modes are scale invariant and follow different power laws for Poisson and Gaussian ensembles, which already during the unfolding allows to distinguish the two cases.


title: Linking accretion flow and particle acceleration in jets - II. Self-similar jet models with full relativistic MHD gravitational mass

abstract: We present a new, semi-analytic formalism to model the acceleration and collimation of relativistic jets in a gravitational potential. The gravitational energy density includes the kinetic, thermal, and electromagnetic mass contributions. The solutions are close to self-similar throughout the integration, from very close to the black hole to the region where gravity is unimportant. The field lines are tied to the conditions very close to the central object and eventually overcollimate, possibly leading to a collimation shock. This collimation shock could provide the conditions for diffusive shock acceleration, leading to the observed electron populations with a power-law energy distribution in jets. We provide the derivation, a detailed analysis of a solution, and describe the effects the parameters have on the properties of the solutions, such as the Lorentz factor and location of the collimation shock. We also discuss the deviations from self-similarity. By comparing the new gravity term with the gravity term obtained from a non-relativistic formalism in a previous work, we show they are equivalent in the non-relativistic limit. This equivalence shows the approach taken in that work is valid and allows us to comment on its limitations.


title: Fusion: A general framework for hierarchical tilings

abstract: One well studied way to construct quasicrystalline tilings is via inflate-and-subdivide (a.k.a. substitution) rules. These produce self-similar tilings--the Penrose, octagonal, and pinwheel tilings are famous examples. We present a different model for generating hierarchical tilings we call "fusion rules". Inflate-and-subdivide rules are a special case of fusion rules, but general fusion rules are more flexible and allow for defects, changes in geometry, and even constrained randomness. A condition that produces homogeneous structures and a method for computing frequency for fusion tiling spaces are discussed.


title: Constraints on Jones transmission matrices from time-reversal invariance and discrete spatial symmetries

abstract: Optical spectroscopies are most often used to probe dynamical correlations in materials, but they are also a probe of symmetry. Polarization anisotropies are of course sensitive to structural anisotropies, but have been much less used as a probe of more exotic symmetry breakings in ordered states. In this paper a Jones transfer matrix formalism is discussed to infer the existence of exotic broken symmetry states of matter from their electrodynamic response for a full complement of possible broken symmetries including reflection, rotation, rotation-reflection, inversion, and time-reversal. A specific condition to distinguish the case of macroscopic time-reversal symmetry breaking is particularly important as in a dynamical experiment like optics, one must distinguish $reciprocity$ from time-reversal symmetry as dissipation violates strict time-reversal symmetry of an experiment. Different forms of reciprocity can be distinguished, but only one is a sufficient (but not necessary) condition for macroscopic time-reversal symmetry breaking. I show the constraints that a Jones matrix develops under the presence or absence of such symmetries. These constraints typically appear in the form of an algebra relating matrix elements or overall constraints (transposition, unitarity, hermiticity, normality, etc.) on the form of the Jones matrix. I work out a number of examples including the trivial case of a ferromagnet, and the less trivial cases of magnetoelectrics and vector and scalar spin "chiral" states. I show that the formalism can be used to demonstrate that Kerr rotation must be absent in time-reversal symmetric chiral materials. The formalism here is discussed with an eye towards its use in time-domain THz spectroscopy in transmission, but with small modifications it is more generally applicable.


title: EM for phylogenetic topology reconstruction on non-homogeneous data

abstract: Background: The reconstruction of the phylogenetic tree topology of four taxa is, still nowadays, one of the main challenges in phylogenetics. Its difficulties lie in considering not too restrictive evolutionary models, and correctly dealing with the long-branch attraction problem. The correct reconstruction of 4-taxon trees is crucial for making quartet-based methods work and being able to recover large phylogenies. Results: In this paper we consider an expectation-maximization method for maximizing the likelihood of (time nonhomogeneous) evolutionary Markov models on trees. We study its success on reconstructing 4-taxon topologies and its performance as input method in quartet-based phylogenetic reconstruction methods such as QFIT and QuartetSuite. Our results show that the method proposed here outperforms neighbor-joining and the usual (time-homogeneous continuous-time) maximum likelihood methods on 4-leaved trees with among-lineage instantaneous rate heterogeneity, and perform similarly to usual continuous-time maximum-likelihood when data satisfies the assumptions of both methods. Conclusions: The method presented in this paper is well suited for reconstructing the topology of any number of taxa via quartet-based methods and is highly accurate, specially regarding largely divergent trees and time nonhomogeneous data.


title: Fundamental stellar parameters and metallicities from Bayesian spectroscopy. Application to low- and high-resolution spectra

abstract: We present a unified framework to derive fundamental stellar parameters by combining all available observational and theoretical information for a star. The algorithm relies on the method of Bayesian inference, which for the first time directly integrates the spectroscopic analysis pipeline based on the global spectrum synthesis and allows for comprehensive and objective error calculations given the priors. Arbitrary input datasets can be included into our analysis and other stellar quantities, in addition to stellar age, effective temperature, surface gravity, and metallicity, can be computed on demand. We lay out the mathematical framework of the method and apply it to several observational datasets, including high- and low-resolution spectra (UVES, NARVAL, HARPS, SDSS/SEGUE). We find that simpler approximations for the spectroscopic PDF, which are inherent to past Bayesian approaches, lead to deviations of several standard deviations and unreliable errors on the same data. By its flexibility and the simultaneous analysis of multiple independent measurements for a star, it will be ideal to analyse and cross-calibrate the large ongoing and forthcoming surveys, like Gaia-ESO, SDSS, Gaia and LSST.


title: CANDELS+3D-HST: compact SFGs at z~2-3, the progenitors of the first quiescent galaxies

abstract: We analyze the star-forming and structural properties of 45 massive (log(M/Msun)>10) compact star-forming galaxies (SFGs) at 2<z<3 to explore whether they are progenitors of compact quiescent galaxies at z~2. The optical/NIR and far-IR Spitzer/Herschel colors indicate that most compact SFGs are heavily obscured. Nearly half (47%) host an X-ray bright AGN. In contrast, only about 10% of other massive galaxies at that time host AGNs. Compact SFGs have centrally-concentrated light profiles and spheroidal morphologies similar to quiescent galaxies, and are thus strikingly different from other SFGs. Most compact SFGs lie either within the SFR-M main sequence (65%) or below (30%), on the expected evolutionary path towards quiescent galaxies. These results show conclusively that galaxies become more compact before they lose their gas and dust, quenching star formation. Using extensive HST photometry from CANDELS and grism spectroscopy from the 3D-HST survey, we model their stellar populations with either exponentially declining (tau) star formation histories (SFHs) or physically-motivated SFHs drawn from semi-analytic models (SAMs). SAMs predict longer formation timescales and older ages ~2 Gyr, which are nearly twice as old as the estimates of the tau models. While both models yield good SED fits, SAM SFHs better match the observed slope and zero point of the SFR-M main sequence. Some low-mass compact SFGs (log(M/Msun)=10-10.6) have younger ages but lower sSFRs than that of more massive galaxies, suggesting that the low-mass galaxies reach the red sequence faster. If the progenitors of compact SFGs are extended SFGs, state-of-the-art SAMs show that mergers and disk instabilities are both able to shrink galaxies, but disk instabilities are more frequent (60% versus 40%) and form more concentrated galaxies. We confirm this result via high-resolution hydrodynamic simulations.


title: Scaling relations and mass bias in hydrodynamical f(R) gravity simulations of galaxy clusters

abstract: We investigate the impact of chameleon-type f(R) gravity models on the properties of galaxy clusters and groups. Our f(R) simulations follow for the first time also the hydrodynamics of the intracluster and intragroup medium. This allows us to assess how f(R) gravity alters the X-ray scaling relations of clusters and how hydrostatic and dynamical mass estimates are biased when modifications of gravity are ignored in their determination. We find that velocity dispersions and intracluster medium temperatures are both increased by up to 1/3 in f(R) gravity in low-mass halos, while the difference disappears in massive objects. The mass scale of the transition depends on the background value f_R0 of the scalar degree of freedom. These changes in temperature and velocity dispersion alter the mass-temperature and X-ray luminosity-temperature scaling relations and bias dynamical and hydrostatic mass estimates that do not explicitly account for modified gravity towards higher values. Recently, a relative enhancement of X-ray compared to weak lensing masses was found by the Planck Collaboration (2013). We demonstrate that an explanation for this offset may be provided by modified gravity and the associated bias effects, which interestingly are of the required size. Finally, we find that the abundance of subhalos at fixed cluster mass is only weakly affected by f(R) gravity.


title: A Merger Shock in Abell 2034

abstract: We present a $250\,$ks Chandra observation of the cluster merger A2034 with the aim of understanding the nature of a sharp edge previously characterized as a cold front. The new data reveal that the edge is coherent over a larger opening angle and is significantly more bow-shock-shaped than previously thought. Within $\sim 27\,$degrees about the axis of symmetry of the edge the density, temperature and pressure drop abruptly by factors of $1.83^{+0.09}_{-0.08}$, $1.85^{+0.41}_{-0.41}$ and $3.4^{+0.8}_{-0.7}$, respectively. This is inconsistent with the pressure equilibrium expected of a cold front and we conclude that the edge is a shock front. We measure a Mach number $M = 1.59^{+0.06}_{-0.07}$ and corresponding shock velocity $v_{\rm shock}\simeq 2057\,$km/s. Using spectra collected at the MMT with the Hectospec multi-object spectrograph we identify 328 spectroscopically confirmed cluster members. Significantly, we find a local peak in the projected galaxy density associated with a bright cluster galaxy which is located just ahead of the nose of the shock. The data are consistent with a merger viewed within $\sim 23\,$degrees of the plane of the sky. The merging subclusters are now moving apart along a north-south axis approximately $0.3\,$Gyr after a small impact parameter core passage. The gas core of the secondary subcluster, which was driving the shock, appears to have been disrupted by the merger. Without a driving 'piston' we speculate that the shock is dying. Finally, we propose that the diffuse radio emission near the shock is due to the revival of pre-existing radio plasma which has been overrun by the shock.


title: 3D simulations of the early stages of AGN jets: geometry, thermodynamics and backflow

abstract: We investigate the interplay between jets from Active Galactic Nuclei (AGNs) and the surrounding InterStellar Medium (ISM) through full 3D, high resolution, Adaptive Mesh Refinement simulations performed with the FLASH code. We follow the jet- ISM system for several Myr in its transition from an early, compact source to an extended one including a large cocoon. During the jet evolution, we identify three major evolutionary stages and we find that, contrary to the prediction of popular theoretical models, none of the simulations shows a self-similar behavior. We also follow the evolution of the energy budget, and find that the fraction of input power deposited into the ISM (the AGN coupling constant) is of order of a few percent during the first few Myr. This is in broad agreement with galaxy formation models employing AGN feedback. However, we find that in these early stages, this energy is deposited only in a small fraction (< 1%) of the total ISM volume. Finally we demonstrate the relevance of backflows arising within the extended cocoon generated by a relativistic AGN jet within the ISM of its host galaxy, previously proposed as a mechanism for self-regulating the gas accretion onto the central object. These backflows tend later to be destabilized by the 3D dynamics, rather than by hydrodynamic (Kelvin- Helmholtz) instabilities. Yet, in the first few hundred thousand years, backflows may create a central accretion region of significant extent, and convey there as much as a few millions of solar masses.


title: Geometric and dynamic distortions in anisotropic galaxy clustering

abstract: We examine the signature of dynamic (redshift-space) distortions and geometric distortions (including the Alcock-Paczynski effect) in the context of the galaxy power spectrum measured in upcoming galaxy redshift surveys. Information comes from both the baryon acoustic oscillation (BAO) feature and the broadband power spectrum shape. Accurate modeling is required to extract this information without introducing systematic biases in the result. We consider an analytic model for the power spectrum of dark matter halos in redshift space, based on the distribution function expansion, and compare with halo clustering measured in N-body simulations. We forecast that the distribution function model is sufficiently accurate to allow the inclusion of broadband information on scales down to k~0.2 h/Mpc, with somewhat better accuracy for higher bias halos. Compared with a BAO-only analysis with reconstruction, including broadband shape information can improve unbiased constraints on distance measures H(z) and D_A(z) by ~30% and 20%, respectively, for a galaxy sample similar to the DESI luminous red galaxies. The gains in precision are larger in the absence of BAO reconstruction. Furthermore, including broadband shape information allows the measurement of structure growth, through redshift-space distortions. For the same galaxy sample, the distribution function model is able to constrain f*sigma_8 to ~2%, when simultaneously fitting for H(z) and D_A(z). We discuss techniques to optimize the analysis of the power spectrum, including removing modes near the line-of-sight that are particularly challenging to model, and whether these approaches can improve parameter constraints. We find that such techniques are unlikely to significantly improve constraints on geometry, although they may allow higher precision measurements of redshift-space distortions.


title: Testing the No-Hair Theorem with Event Horizon Telescope Observations of Sagittarius A*

abstract: The advent of the Event Horizon Telescope (EHT), a millimeter-wave very-long baseline interferometric array, has enabled spatially-resolved studies of the sub-horizon-scale structure for a handful of supermassive black holes. Among these, the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), presents the largest angular cross section. Thus far, these studies have focused upon measurements of the black hole spin and the validation of low-luminosity accretion models. However, a critical input into the analysis of EHT data is the structure of the black hole spacetime, and thus these observations provide the novel opportunity to test the applicability of the Kerr metric to astrophysical black holes. Here we present the first simulated images of a radiatively inefficient accretion flow (RIAF) around Sgr A* employing a quasi-Kerr metric that contains an independent quadrupole moment in addition to the mass and spin that fully characterize a black hole in general relativity. We show that these images can be significantly different from the images of a RIAF around a Kerr black hole with the same spin and demonstrate the feasibility of testing the no-hair theorem by constraining the quadrupolar deviation from the Kerr metric with existing EHT data. Equally important, we find that the disk inclination and spin orientation angles are robust to the inclusion of additional parameters, providing confidence in previous estimations assuming the Kerr metric based upon EHT observations. However, at present the limits upon potential modifications of the Kerr metric remain weak.


title: Physical Properties of the Inner Shocks in Hot, Tilted Black Hole Accretion Flows

abstract: Simulations of hot, pressure supported, tilted black hole accretion flows, in which the angular momentum of the flow is misaligned with the black hole spin axis, can exhibit two non-axisymmetric shock structures in the inner regions of the flow. We analyze the strength and significance of these shock structures in simulations with tilt angles of 10 and 15 degrees. By integrating fluid trajectories in the simulations through the shocks, and tracking the variations of fluid quantities along these trajectories, we show that these shocks are strong, with substantial compression ratios, in contrast to earlier claims. However, they are only moderately relativistic. We also show that the two density enhancements resembling flow streams in their shape are in fact merely post-shock compressions, as fluid trajectories cut across, rather than flow along, them. The dissipation associated with the shocks is a substantial fraction ($\simeq3-12$ percent) of the rest mass energy advected into the hole, and therefore comparable to the dissipation expected from turbulence. The shocks should therefore make order unity changes in the observed properties of black hole accretion flows that are tilted.


title: MyGIsFOS: an automated code for parameter determination and detailed abundance analysis in cool stars

abstract: The current and planned high-resolution, high-multiplexity stellar spectroscopic surveys, as well as the swelling amount of under-utilized data present in public archives have led to an increasing number of efforts to automate the crucial but slow process to retrieve stellar parameters and chemical abundances from spectra. We present MyGIsFOS, a code designed to derive atmospheric parameters and detailed stellar abundances from medium - high resolution spectra of cool (FGK) stars. We describe the general structure and workings of the code, present analyses of a number of well studied stars representative of the parameter space MyGIsFOS is designed to cover, and examples of the exploitation of MyGIsFOS very fast analysis to assess uncertainties through Montecarlo tests. MyGIsFOS aims to reproduce a ``traditional'' manual analysis by fitting spectral features for different elements against a precomputed grid of synthetic spectra. Fe I and Fe II lines can be employed to determine temperature, gravity, microturbulence, and metallicity by iteratively minimizing the dependence of Fe I abundance from line lower energy and equivalent width, and imposing Fe I - Fe II ionization equilibrium. Once parameters are retrieved, detailed chemical abundances are measured from lines of other elements. MyGIsFOS replicates closely the results obtained in similar analyses on a set of well known stars. It is also quite fast, performing a full parameter determination and detailed abundance analysis in about two minutes per star on a mainstream desktop computer. Currently, its preferred field of application are high-resolution and/or large spectral coverage data (e.g UVES, X-Shooter, HARPS, Sophie).


title: On Constructing Constrained Tree Automata Recognizing Ground Instances of Constrained Terms

abstract: An inductive theorem proving method for constrained term rewriting systems, which is based on rewriting induction, needs a decision procedure for reduction-completeness of constrained terms. In addition, the sufficient complete property of constrained term rewriting systems enables us to relax the side conditions of some inference rules in the proving method. These two properties can be reduced to intersection emptiness problems related to sets of ground instances for constrained terms. This paper proposes a method to construct deterministic, complete, and constraint-complete constrained tree automata recognizing ground instances of constrained terms.


title: Random Generation of Nondeterministic Finite-State Tree Automata

abstract: Algorithms for (nondeterministic) finite-state tree automata (FTAs) are often tested on random FTAs, in which all internal transitions are equiprobable. The run-time results obtained in this manner are usually overly optimistic as most such generated random FTAs are trivial in the sense that the number of states of an equivalent minimal deterministic FTA is extremely small. It is demonstrated that nontrivial random FTAs are obtained only for a narrow band of transition probabilities. Moreover, an analytic analysis yields a formula to approximate the transition probability that yields the most complex random FTAs, which should be used in experiments.


title: Signatures of spin-triplet excitations in optical conductivity of valence bond solids

abstract: We show that the optical responses below the Mott gap can be used to probe the spin-triplet excitations in valence bond solid (VBS) phases in Mott insulators. The optical conductivity in this regime arises due to the electronic polarization mechanism via virtual electron hopping processes. We apply this mechanism to the Hubbard model with spin-orbit couplings and/or the corresponding spin model with significant Dzyaloshinskii-Moriya (DM) interactions, and compute the optical conductivity of VBS states on both ideal and deformed Kagome lattices. In case of the deformed Kagome lattice, we study the antiferromagnet, Rb$_2$Cu$_3$SnF$_{12}$ with the pinwheel VBS state. In case of the ideal Kagome lattice, we explore the optical conductivity signatures of the spin-triplet excitations for three VBS states with (1) a 12-site unit cell, (2) a 36-site unit cell with six-fold rotation symmetry, and (3) a 36-site unit cell with three-fold rotation symmetry, respectively. We find that increasing the DM interactions generally leads to broad and smooth features in the optical conductivity with interesting experimental consequences. The optical conductivity reflects the features of the spin-triplet excitations that can be measured in future experiments.


title: Relativistic Self-similar Dynamic Collapses of Black Holes in General Polytropic Spherical Clouds

abstract: We study the hydrodynamic self-similar mass collapses of general polytropic (GP) spherical clouds to central Schwarzschild black holes and void evolution with or without shocks. In order to grossly capture characteristic effects of general relativity (GR) outside yet close to the event horizon of a Schwarzschild black hole and to avoid mathematical complexity, we adopt the approximation of the Paczynski-Wiita gravity to replace the simple Newtonian gravity in our model formulation. A new dimensionless parameter s appears with the physical meaning of the square of the ratio of the sound speed to the speed of light $c$. Various self-similar dynamic solutions are constructed for a polytropic index $\gamma>4/3$. Two (for small enough $s<1$) or no (for large enough $s<1$) expansion-wave collapse solutions (EWCSs) with central event horizons exist when $\gamma>4/3$, representing the collapse of static singular GP spheres towards the central singularity of spacetime. Such GP spherical dynamic mass collapse is shown to be highly efficient for the rapid formation of supermassive black holes (SMBHs; mass range of $10^6-10^{10}M_{\odot}$) in the early universe or even hypermassive black holes (HMBHs; mass range of $10^{10}-10^{12}M_{\odot}$) if extremely massive mass reservoirs could be sustained for a sufficiently long time, which may evolve into hard X-ray/gamma ray sources or quasars according to their surroundings. Self-similar dynamic solutions of a GP gas are also proposed for the stellar mass black hole formation during the violent supernova explosion of a massive progenitor star, the timescale of which is estimated of $10^{-3}$ seconds. Rebound shocks travelling in supernovae are also discussed based on our self-similar shock expansion solutions.


title: Expressiveness of Visibly Pushdown Transducers

abstract: Visibly pushdown transducers (VPTs) are visibly pushdown automata extended with outputs. They have been introduced to model transformations of nested words, i.e. words with a call/return structure. As trees and more generally hedges can be linearized into (well) nested words, VPTs are a natural formalism to express tree transformations evaluated in streaming. This paper aims at characterizing precisely the expressive power of VPTs with respect to other tree transducer models.


title: Node Query Preservation for Deterministic Linear Top-Down Tree Transducers

abstract: This paper discusses the decidability of node query preservation problems for XML document transformations. We assume a transformation given by a deterministic linear top-down data tree transducer (abbreviated as DLT^V) and an n-ary query based on runs of a tree automaton. We say that a DLT^V Tr strongly preserves a query Q if there is a query Q' such that for every document t, the answer set of Q' for Tr(t) is equal to the answer set of Q for t. Also we say that Tr weakly preserves Q if there is a query Q' such that for every t_d in the range of Tr, the answer set of Q' for t_d is equal to the union of the answer set of Q for t such that t_d = Tr(t). We show that the weak preservation problem is coNP-complete and the strong preservation problem is in 2-EXPTIME.


title: XPath Node Selection over Grammar-Compressed Trees

abstract: XML document markup is highly repetitive and therefore well compressible using grammar-based compression. Downward, navigational XPath can be executed over grammar-compressed trees in PTIME: the query is translated into an automaton which is executed in one pass over the grammar. This result is well-known and has been mentioned before. Here we present precise bounds on the time complexity of this problem, in terms of big-O notation. For a given grammar and XPath query, we consider three different tasks: (1) to count the number of nodes selected by the query, (2) to materialize the pre-order numbers of the selected nodes, and (3) to serialize the subtrees at the selected nodes.


title: Mechanism of Ultrafast Relaxation of a Photo-Carrier in Antiferromagnetic Spin Background

abstract: We study the relaxation mechanism of a highly excited carrier propagating in the antiferromagnetic background modeled by the $t$-$J$ Hamiltonian on a square lattice. We show that the relaxation consists of two distinct stages. The initial ultrafast stage with the relaxation time $\tau\sim (\hbar/t_0)(J/t_0)^{-2/3}$ (where $t_0$ is the hopping integral and $J$ is the exchange interaction) is based on generation of string states in the close proximity of the carrier. This unusual scaling of $\tau$ is obtained by means of comparison of numerical results with a simplified $t$-$J_z$ model on a Bethe lattice. In the subsequent (much slower) stage local spin excitations are carried away by magnons. The relaxation time on the two-leg ladder system is an order of magnitude longer due to the lack of string excitations. This further reinforces the importance of string excitations for the ultrafast relaxation in the two-dimensional system.


title: PAMELA and AMS-02 positron and electron spectra are reproduced by three-dimensional cosmic-ray modeling

abstract: The PAMELA collaboration recently released the $e^+$ absolute spectrum between 1 and 300 GeV in addition to the positron fraction and the $e^-$ spectrum previously measured in the same period. We use the newly developed three-dimensional upgrade of the DRAGON package to model those data. This code allows us to consider a realistic spiral arm source distribution in the Galaxy, which impacts the high-energy shape of the propagated spectra. At low energy we treat solar modulation with the HelioProp code and compare its results with those obtained using the usual force-field approximation. We show that all PAMELA data sets can be consistently, and accurately, described in terms of a standard background on top of which a charge symmetric $e^+$ + $e^-$ extra component with harder injection spectrum is added; this extra contribution is peaked at about 1 - 10 TeV and may originate from a diffuse population of sources located in the Galactic arms. For the first time, we compute the energy required to sustain such a relevant positron flux in the Galaxy, finding that it is naturally compatible with an astrophysical origin. We considered several reference propagation setups; we find that models with a low (or null) reacceleration - tuned against light nuclei data - nicely describe both PAMELA leptonic and hadronic data with no need to introduce a low-energy break in the proton and helium spectra, as it would be required for high reacceleration models. We also compare our models with the preliminary $e^-$ and $e^+$ absolute spectra recently measured by AMS-02. We find that those data, differently from what is inferred from the positron fraction alone, favor a high energy cutoff (about 10 TeV) of the extra component if this is uniquely generated in the Galactic arms. A lower cutoff may be allowed if a relevant contribution from powerful $e^-$ + $e^+$ nearby accelerators is invoked.


title: Bayesian quantum frequency estimation in presence of collective dephasing

abstract: We advocate a Bayesian approach to optimal quantum frequency estimation - an important issue for future quantum enhanced atomic clock operation. The approach provides a clear insight into the interplay between decoherence and the extent of the prior knowledge in determining the optimal interrogation times and optimal estimation strategies. We propose a general framework capable of describing local oscillator noise as well as additional collective atomic dephasing effects. For a Gaussian noise the average Bayesian cost can be expressed using the quantum Fisher information and thus we establish a direct link between the two, often competing, approaches to quantum estimation theory


title: Effects of Fluid Velocity Gradients on Heavy Quark Energy Loss

abstract: We use holographic duality to analyze the drag force on, and consequent energy loss of, a heavy quark moving through a strongly coupled conformal fluid with non-vanishing gradients in its velocity and temperature. We derive the general expression for the drag force to first order in the fluid gradients. Using this general expression, we show that a quark that is instantaneously at rest, relative to the fluid, in a fluid whose velocity is changing with time feels a nonzero force. And, we show that for a quark that is moving ultra-relativistically, the first order gradient "corrections" become larger than the zeroth order drag force, suggesting that the gradient expansion may be unreliable in this regime. We illustrate the importance of the fluid gradients for heavy quark energy loss by considering a fluid with one-dimensional boost invariant Bjorken expansion as well as the strongly coupled plasma created by colliding sheets of energy.


title: The proportion of plane cubic curves over ${\mathbb Q}$ that everywhere locally have a point

abstract: We show that the proportion of plane cubic curves over ${\mathbb Q}_p$ that have a ${\mathbb Q}_p$-rational point is a rational function in $p$, where the rational function is independent of $p$, and we determine this rational function explicitly. As a consequence, we obtain the density of plane cubic curves over ${\mathbb Q}$ that have points everywhere locally; numerically, this density is shown to be $\approx 97.3\%$.


title: Monte Carlo studies on the expanding behavior of the early universe in the Lorentzian type IIB matrix model

abstract: The type IIB matrix model is a conjectured nonperturbative formulation of superstring theory. Recently the Lorentzian version of the model has been studied by Monte Carlo simulation, and it has been shown that only three out of nine spatial directions start to expand after a critical time. We extend this work by investigating the expanding behavior for much longer time. We find that the 3d space expands exponentially for some period of time, which may be interpreted as inflation. We also simulate a simplified model, which is expected to capture some qualitative features of the original model at much later times. We observe that the exponential expansion eventually changes into a power-law t^(1/2) behavior, which agrees with the expanding behavior of the Friedmann-Robertson-Walker (FRW) universe in the radiation dominated era.


title: Quenched disorder and vestigial nematicity in the pseudo-gap regime of the cuprates

abstract: We have carried out a theoretical analysis of the Landau-Ginzburg-Wilson effective field theory of a classical incommensurate CDW in the presence of weak quenched disorder. While the possibility a sharp phase transition and long-range CDW order are precluded in such systems, we show that any discrete symmetry breaking aspect of the charge order -- nematicity in the case of the unidirectional (stripe) CDW we consider explicitly -- generically survives up to a non-zero critical disorder strength. Such "vestigial order", which is subject to unambiguous macroscopic detection, can serve as an avatar of what would be CDW order in the ideal, zero disorder limit. Various recent experiments in the pseudo-gap regime of the hole-doped cuprate high-temperature superconductors are readily interpreted in light of these results.


title: The First Pulse of the Extremely Bright GRB 130427A: A Test Lab for Synchrotron Shocks

abstract: Gamma-ray burst (GRB) 130427A is one of the most energetic GRBs ever observed. The initial pulse up to 2.5 s is possibly the brightest well-isolated pulse observed to date. A fine time resolution spectral analysis shows power-law decays of the peak energy from the onset of the pulse, consistent with models of internal synchrotron shock pulses. However, a strongly correlated power-law behavior is observed between the luminosity and the spectral peak energy that is inconsistent with curvature effects arising in the relativistic outflow. It is difficult for any of the existing models to account for all of the observed spectral and temporal behaviors simultaneously.


title: Parallel Tempering Simulation of the three-dimensional Edwards-Anderson Model with Compact Asynchronous Multispin Coding on GPU

abstract: Monte Carlo simulations of the Ising model play an important role in the field of computational statistical physics, and they have revealed many properties of the model over the past few decades. However, the effect of frustration due to random disorder, in particular the possible spin glass phase, remains a crucial but poorly understood problem. One of the obstacles in the Monte Carlo simulation of random frustrated systems is their long relaxation time making an efficient parallel implementation on state-of-the-art computation platforms highly desirable. The Graphics Processing Unit (GPU) is such a platform that provides an opportunity to significantly enhance the computational performance and thus gain new insight into this problem. In this paper, we present optimization and tuning approaches for the CUDA implementation of the spin glass simulation on GPUs. We discuss the integration of various design alternatives, such as GPU kernel construction with minimal communication, memory tiling, and look-up tables. We present a binary data format, Compact Asynchronous Multispin Coding (CAMSC), which provides an additional $28.4\%$ speedup compared with the traditionally used Asynchronous Multispin Coding (AMSC). Our overall design sustains a performance of 33.5 picoseconds per spin flip attempt for simulating the three-dimensional Edwards-Anderson model with parallel tempering, which significantly improves the performance over existing GPU implementations.


title: NSMAXG: A new magnetic neutron star spectral model in XSPEC

abstract: The excellent sensitivity of X-ray telescopes, such as Chandra and XMM-Newton, is ideal for the study of cooling neutron stars, which can emit at these energies. In order to exploit the wealth of information contained in the high quality data, a thorough knowledge of the radiative properties of neutron star atmospheres is necessary. A key factor affecting photon emission is magnetic fields, and neutron stars are known to have strong surface magnetic fields. Here I briefly describe our latest work on constructing magnetic (B >= 10^10 G) atmosphere models of neutron stars and the NSMAXG implementation of these models in XSPEC. Our results allow for more robust extractions of neutron star parameters from observations.


title: On the analysis of a coupled kinetic-fluid model with local alignment forces

abstract: This paper studies global existence, hydrodynamic limit, and large-time behavior of weak solutions to a kinetic flocking model coupled to the incompressible Navier-Stokes equations. The model describes the motion of particles immersed in a Navier-Stokes fluid interacting through local alignment. We first prove the existence of weak solutions using energy and $L^p$ estimates together with the velocity averaging lemma. We also rigorously establish a hydrodynamic limit corresponding to strong noise and local alignment. In this limit, the dynamics can be totally described by a coupled compressible Euler - incompressible Navier-Stokes system. The proof is via relative entropy techniques. Finally, we show a conditional result on the large-time behavior of classical solutions. Specifically, if the mass-density satisfies a uniform in time integrability estimate, then particles align with the fluid velocity exponentially fast without any further assumption on the viscosity of the fluid.


title: On the $N_f$-dependence of gluonic observables

abstract: We compute $t_0$, $w_0$ and the topological susceptibility, defined at finite gradient flow time for two-flavour QCD. The use of three lattice spacings and pion masses between 192 and 500 MeV together with a careful error analysis allow to approach the continuum limit of the two-flavour theory despite significant auto-correlations. A comparison to $N_f=0$ results shows the size of sea quark effects in $t_0^2\chi$, with $\chi$ the topological susceptibility, and low energy observables such as $t_0/w_0^2$ and $t_0/r_0^2$.


title: Active Galactic Nuclei, Neutrinos, and Interacting Cosmic Rays in NGC 253 & NGC 1068

abstract: The galaxies M82, NGC 253, NGC 1068, and NGC 4945 have been detected in gamma-rays by Fermi. Previously, we developed and tested a model for cosmic ray interactions in the starburst galaxy M82. Now, we aim to explore the differences between starburst and active galactic nuclei (AGN) environments by applying our self-consistent model to the starburst galaxy NGC 253 and the Seyfert galaxy NGC 1068. Assuming constant cosmic-ray acceleration efficiency by supernova remnants with Milky-Way parameters, we calculate the cosmic-ray proton and primary and secondary electron/positron populations, predict the radio and gamma-ray spectra, and compare with published measurements. We find that our models easily fits the observed gamma-ray spectrum for NGC 253 while constraining the cosmic ray source spectral index and acceleration efficiency. However, we encountered difficultly modeling the observed radio data and constraining the speed of the galactic wind and the magnetic field strength, unless the gas mass is less than currently preferred values. Additionally, our starburst model consistently underestimates the observed gamma-ray flux and overestimates the radio flux for NGC 1068; these issues would be resolved if the AGN is the primary source of gamma-rays. We discuss the implications of these results and make predictions for the neutrino fluxes for both galaxies.


title: Dynamic Integration of ALM Tools for Agile Software Development

abstract: The paper describes the need for and goals of tool-integration within software development processes. In particular we focus on agile software development but are not limited to. The integration of tools and data between the different domains of the process is essential for an efficient, effective and customized software development. We describe what the next steps in the pursuit of integration are and how major goals can be achieved. Beyond theoretical and architectural considerations we describe the prototypical implementation of an open platform approach. The paper introduces platform apps and a functionality store as general concepts to make apps and their functionalities available to the community. We describe the implementation of the approach and how it can be practically utilized. The description is based on one major use case and further steps are motivated by various other examples.


title: The Dust Scattering Halo of Cygnus X-3

abstract: Dust grains scatter X-ray light through small angles, producing a diffuse halo image around bright X-ray point sources situated behind a large amount of interstellar material. We present analytic solutions to the integral for the dust scattering intensity, which allow for a Bayesian analysis of the scattering halo around Cygnus X-3. Fitting the optically thin 4-6 keV halo surface brightness profile yields the dust grain size and spatial distribution. We assume a power law distribution of grain sizes ($n \propto a^{-p}$) and fit for $p$, the grain radius cut-off amax, and dust mass column. We find that a $p \approx 3.5$ dust grain size distribution with amax $\approx 0.2$ um fits the halo profile relatively well, whether the dust is distributed uniformly along the line of sight or in clumps. We find that a model consisting of two dust screens, representative of foreground spiral arms, requires the foreground Perseus arm to contain 80\% of the total dust mass. The remaining 20\% of the dust, which may be associated with the outer spiral arm of the Milky Way, is located within 1 kpc of Cyg X-3. Regardless of which model was used, we found $\tau_{\rm sca} \sim 2 \ E_{\rm keV}^{-2}$. We examine the energy resolved halos of Cyg X-3 from 1 - 6 keV and find that there is a sharp drop in scattering halo intensity when $E < 2-3$ keV, which cannot be explained with multiple scattering effects. We hypothesize that this may be caused by large dust grains or material with unique dielectric properties, causing the scattering cross-section to depart from the Rayleigh-Gans approximation that is used most often in X-ray scattering studies. The foreground Cyg OB2 association, which contains several evolved stars with large extinction values, is a likely culprit for grains of unique size or composition.


title: Transmission electron microscopy and ferromagnetic resonance investigations of tunnel magnetic junctions using Co2MnGe Heusler alloy as magnetic electrodes

abstract: HRTEM, nano-beam electronic diffraction, energy dispersive X-rays scanning spectroscopy, Vibrating Sample Magnetometry (VSM) and FerroMagnetic Resonance (FMR) techniques are used in view of comparing (static and dynamic) magnetic and structural properties of Co2MnGe (13 nm)/Al2O3 (3 nm)/Co (13 nm) tunnel magnetic junctions (TMJ), deposited on various single crystalline substrates (a-plane sapphire, MgO(100) and Si(111)). They allow for providing a correlation between these magnetic properties and the fine structure investigated at atomic scale. The Al2O3 tunnel barrier is always amorphous and contains a large concentration of Co atoms, which, however, is significantly reduced when using a sapphire substrate. The Co layer is polycrystalline and shows larger grains for films grown on a sapphire substrate. The VSM investigation reveals in-plane anisotropy only for samples grown on a sapphire substrate. The FMR spectra of the TMJs are compared to the obtained ones with a single Co and Co2MnGe films of identical thickness deposited on a sapphire substrate. As expected, two distinct modes are detected in the TMJs while only one mode is observed in each single film. For the TMJ grown on a sapphire substrate the FMR behavior does not significantly differ from the superposition of the individual spectra of the single films, allowing for concluding that the exchange coupling between the two magnetic layers is too small to give rise to observable shifts. For TMJs grown on a Si or on a MgO substrate the resonance spectra reveal one mode which is nearly identical to the obtained one in the single Co film, while the other observed resonance shows a considerably smaller intensity and cannot be described using the magnetic parameters appropriate to the single Co2MnGe film.


title: Adaptive Learning of Region-based pLSA Model for Total Scene Annotation

abstract: In this paper, we present a region-based pLSA model to accomplish the task of total scene annotation. To be more specific, we not only properly generate a list of tags for each image, but also localizing each region with its corresponding tag. We integrate advantages of different existing region-based works: employ efficient and powerful JSEG algorithm for segmentation so that each region can easily express meaningful object information; the introduction of pLSA model can help better capturing semantic information behind the low-level features. Moreover, we also propose an adaptive padding mechanism to automatically choose the optimal padding strategy for each region, which directly increases the overall system performance. Finally we conduct 3 experiments to verify our ideas on Corel database and demonstrate the effectiveness and accuracy of our system.


title: PANDA: Pose Aligned Networks for Deep Attribute Modeling

abstract: We propose a method for inferring human attributes (such as gender, hair style, clothes style, expression, action) from images of people under large variation of viewpoint, pose, appearance, articulation and occlusion. Convolutional Neural Nets (CNN) have been shown to perform very well on large scale object recognition problems. In the context of attribute classification, however, the signal is often subtle and it may cover only a small part of the image, while the image is dominated by the effects of pose and viewpoint. Discounting for pose variation would require training on very large labeled datasets which are not presently available. Part-based models, such as poselets and DPM have been shown to perform well for this problem but they are limited by shallow low-level features. We propose a new method which combines part-based models and deep learning by training pose-normalized CNNs. We show substantial improvement vs. state-of-the-art methods on challenging attribute classification tasks in unconstrained settings. Experiments confirm that our method outperforms both the best part-based methods on this problem and conventional CNNs trained on the full bounding box of the person.


title: On multiple peaks and moderate deviations for supremum of Gaussian field

abstract: We prove two theorems concerning extreme values of general Gaussian fields. Our first theorem concerns with the concept of multiple peaks. A theorem of Chatterjee states that when a centered Gaussian field admits the so-called superconcentration property, it typically attains values near its maximum on multiple near-orthogonal sites, known as multiple peaks. We improve his theorem in two aspects: (i) the number of peaks attained by our bound is of the order $\exp(c / \sigma^2)$ (as opposed to Chatterjee's polynomial bound in $1/\sigma$), where $\sigma$ is the standard deviation of the supremum of the Gaussian field, which is assumed to have variance at most $1$ and (ii) our bound need not assume that the correlations are non-negative. We also prove a similar result based on the superconcentration of the free energy. As primary applications, we infer that for the S-K spin glass model on the $n$-hypercube and directed polymers on $\mathbb{Z}_n^2$, there are polynomially (in $n$) many near-orthogonal sites that achieve values near their respective maxima. Our second theorem gives an upper bound on moderate deviation for the supremum of a general Gaussian field. While the Gaussian isoperimetric inequality implies a sub-Gaussian concentration bound for the supremum, we show that the exponent in that bound can be improved under the assumption that the expectation of the supremum is of the same order as that of the independent case.


title: Hydroxyl, water, ammonia, carbon monoxide and neutral carbon towards the Sgr A complex

abstract: We observed OH, H$_2$O, HN$_3$, C$^{18}$O, and C$_I$ towards the +50 km/s cloud (M-0.02-0.07), the CND and the +20 km/s (M-0.13-0.08) cloud in the Sgr A complex with the VLA, Odin and SEST. Strong OH absorption, H$_2$O emission and absorption lines were seen at all three positions. Strong C$^{18}$O emissions were seen towards the +50 and +20 km/s clouds. The CND is rich in H$_2$O and OH, and these abundances are considerably higher than in the surrounding clouds, indicating that shocks, star formation and clump collisions prevail in those objects. A comparison with the literature reveals that it is likely that PDR chemistry including grain surface reactions, and perhaps also the influences of shocks has led to the observed abundances of the observed molecular species studied here. In the redward high-velocity line wings of both the +50 and +20 km/s clouds and the CND, the very high H$_2$O abundances are suggested to be caused by the combined action of shock desorption from icy grain mantles and high-temperature, gas-phase shock chemistry. Only three of the molecules are briefly discussed here. For OH and H$_2$O three of the nine observed positions are shown, while a map of the C$^{18}$O emission is provided. An extensive paper was recently published with Open Access (Karlsson et al. 2013; http://www.aanda.org/articles/aa/pdf/2013/06/aa20471-12.pdf ).


title: The virialization density of peaks with general density profiles under spherical collapse

abstract: We calculate the non-linear virialization density, $\Delta_c$, of halos under spherical collapse from peaks with an arbitrary initial and final density profile. This is in contrast to the standard calculation of $\Delta_c$ which assumes top-hat profiles. Given our formalism, the non-linear halo density can be calculated once the shape of the initial peak's density profile and the shape of the virialized halo's profile are provided. We solve for $\Delta_c$ for halos in an Einstein de-Sitter and $\Lambda$CDM universe. As examples, we consider power-law initial profiles as well as spherically averaged peak profiles calculated from the statistics of a Gaussian random field. We find that, depending on the profiles used, $\Delta_c$ is smaller by a factor of a few to as much as a factor of 10 as compared to the density given by the standard calculation ($\approx 200$). Using our results, we show that, for halo finding algorithms that identify halos through an over-density threshold, the halo mass function measured from cosmological simulations can be enhanced at all halo masses by a factor of a few. This difference could be important when using numerical simulations to assess the validity of analytic models of the halo mass function.


title: On Nonrigid Shape Similarity and Correspondence

abstract: An important operation in geometry processing is finding the correspondences between pairs of shapes. The Gromov-Hausdorff distance, a measure of dissimilarity between metric spaces, has been found to be highly useful for nonrigid shape comparison. Here, we explore the applicability of related shape similarity measures to the problem of shape correspondence, adopting spectral type distances. We propose to evaluate the spectral kernel distance, the spectral embedding distance and the novel spectral quasi-conformal distance, comparing the manifolds from different viewpoints. By matching the shapes in the spectral domain, important attributes of surface structure are being aligned. For the purpose of testing our ideas, we introduce a fully automatic framework for finding intrinsic correspondence between two shapes. The proposed method achieves state-of-the-art results on the Princeton isometric shape matching protocol applied, as usual, to the TOSCA and SCAPE benchmarks.


title: Shock Reflection-Diffraction, von Neumann's Conjectures and Nonlinear Equations of Mixed Type

abstract: Shock waves are fundamental in nature. One of the most fundamental problems in fluid mechanics is shock reflection-diffraction by wedges. The complexity of reflection-diffraction configurations was first reported by Ernst Mach in 1878. The problems remained dormant until the 1940s when John von Neumann, as well as other mathematical/experimental scientists, began extensive research into all aspects of shock reflection-diffraction phenomena. In this paper we start with shock reflection-diffraction phenomena and historic perspectives, their fundamental scientific issues and theoretical roles in the mathematical theory of hyperbolic systems of conservation laws. Then we present how the global shock reflection-diffraction problem can be formulated as a boundary value problem in an unbounded domain for nonlinear conservation laws of mixed hyperbolic-elliptic type, and describe the von Neumann conjectures: the sonic conjecture and the detachment conjecture. Finally we discuss some recent developments in solving the von Neumann conjectures and establishing a mathematical theory of shock reflection-diffraction, including the existence, regularity, and stability of global regular configurations of shock reflection-diffraction by wedges.


title: Large-volume results in SU(2) with adjoint fermions

abstract: Taming finite-volume effects is a crucial ingredient in order to identify the existence of IR fixed points. We present the latest results from our numerical simulations of SU(2) gauge theory with 2 Dirac fermions in the adjoint representation on large volumes. We compare with previous results, and extrapolate to thermodynamic limit when possible.


title: The Glauber-Sudarshan and Kirkwood-Rihaczec functions

abstract: It is shown how to write the Kirkwood-Rihaczec quasiprobability distribution as an expectation value of the vacuum state. We do this, by writing the position eigenstates as a "displacment" of the vacumm. We also give a relation between the Glauber-Sudarshan and Kirkwood-Rihaczec quasiprobability distributions.


title: Compressive Measurement Designs for Estimating Structured Signals in Structured Clutter: A Bayesian Experimental Design Approach

abstract: This work considers an estimation task in compressive sensing, where the goal is to estimate an unknown signal from compressive measurements that are corrupted by additive pre-measurement noise (interference, or clutter) as well as post-measurement noise, in the specific setting where some (perhaps limited) prior knowledge on the signal, interference, and noise is available. The specific aim here is to devise a strategy for incorporating this prior information into the design of an appropriate compressive measurement strategy. Here, the prior information is interpreted as statistics of a prior distribution on the relevant quantities, and an approach based on Bayesian Experimental Design is proposed. Experimental results on synthetic data demonstrate that the proposed approach outperforms traditional random compressive measurement designs, which are agnostic to the prior information, as well as several other knowledge-enhanced sensing matrix designs based on more heuristic notions.


title: Model for the Atomic Dielectric Response in Time Dependent Laser Fields

abstract: A nonlocal quantum model is presented for calculating the atomic dielectric response to a strong laser electric field. By replacing the Coulomb potential with a nonlocal potential in the Schrodinger equation, a 3+1D calculation of the time-dependent electric dipole moment can be replaced with a 0+1D integral equation, offering significant computational savings. The model is benchmarked against an established ionization model and \textit{ab initio} simulation of the time-dependent Schrodinger equation. The reduced computational overhead makes the model a promising candidate to incorporate full quantum mechanical time dynamics in laser pulse propagation simulations.

categories: physics.atom-ph physics.optics physics.plasm-ph quant-ph


title: Cu(Ir$_{1-x}$Cr$_x$)$_2$S$_4$: a model system for studying nanoscale phase coexistence at the metal-insulator transition

abstract: Increasingly, nanoscale phase coexistence and hidden broken symmetry states are being found in the vicinity of metal-insulator transitions (MIT), for example, in high temperature superconductors, heavy fermion and colossal magnetoresistive materials, but their importance and possible role in the MIT and related emergent behaviors is not understood. Despite their ubiquity, they are hard to study because they produce weak diffuse signals in most measurements. Here we propose Cu(Ir$_{1-x}$Cr$_x$)$_2$S$_4$ as a model system, where robust local structural signals lead to key new insights. We demonstrate a hitherto unobserved coexistence of a Ir$^{4+}$ charge-localized dimer phase and Cr-ferromagnetism. The resulting phase diagram that takes into account the short range dimer order, is highly reminiscent of a generic MIT phase diagram similar to the cuprates. We suggest that the presence of quenched strain from dopant ions acts as an arbiter deciding between the competing ground states.


title: General entropy-like uncertainty relations in finite dimensions

abstract: We revisit entropic formulations of the uncertainty principle for an arbitrary pair of positive operator-valued measures (POVM) $A$ and $B$, acting on finite dimensional Hilbert space. Salicr\'u generalized $(h,\phi)$-entropies, including R\'enyi and Tsallis ones among others, are used as uncertainty measures associated with the distribution probabilities corresponding to the outcomes of the observables. We obtain a nontrivial lower bound for the sum of generalized entropies for any pair of entropic functionals, which is valid for both pure and mixed states. The bound depends on the overlap triplet $(c_A,c_B,c_{A,B})$ with $c_A$ (resp. $c_B$) being the overlap between the elements of the POVM $A$ (resp. $B$) and $c_{A,B}$ the overlap between the pair of POVM. Our approach is inspired by that of de Vicente and S\'anchez-Ruiz [Phys.\ Rev.\ A \textbf{77}, 042110 (2008)] and consists in a minimization of the entropy sum subject to the Landau-Pollak inequality that links the maximum probabilities of both observables. We solve the constrained optimization problem in a geometrical way and furthermore, when dealing with R\'enyi or Tsallis entropic formulations of the uncertainty principle, we overcome the H\"older conjugacy constraint imposed on the entropic indices by the Riesz-Thorin theorem. In the case of nondegenerate observables, we show that for given $c_{A,B} > \frac{1}{\sqrt2}$, the bound obtained is optimal; and that, for R\'enyi entropies, our bound improves Deutsch one, but Maassen-Uffink bound prevails when $c_{A,B} \leq\frac12$. Finally, we illustrate by comparing our bound with known previous results in particular cases of R\'enyi and Tsallis entropies.


title: Submodular Welfare Maximization

abstract: An overview of different variants of the submodular welfare maximization problem in combinatorial auctions. In particular, I studied the existing algorithmic and game theoretic results for submodular welfare maximization problem and its applications in other areas such as social networks.


title: Estimation of Extreme Quantiles for Functions of Dependent Random Variables

abstract: We propose a new method for estimating the extreme quantiles for a function of several dependent random variables. In contrast to the conventional approach based on extreme value theory, we do not impose the condition that the tail of the underlying distribution admits an approximate parametric form, and, furthermore, our estimation makes use of the full observed data. The proposed method is semiparametric as no parametric forms are assumed on all the marginal distributions. But we select appropriate bivariate copulas to model the joint dependence structure by taking the advantage of the recent development in constructing large dimensional vine copulas. Consequently a sample quantile resulted from a large bootstrap sample drawn from the fitted joint distribution is taken as the estimator for the extreme quantile. This estimator is proved to be consistent. The reliable and robust performance of the proposed method is further illustrated by simulation.


title: Observing interferences between past and future quantum states in resonance fluorescence

abstract: In quantum physics, measurement results are random but their statistics can be predicted assuming some knowledge about the system in the past. Additional knowledge from a future measurement deeply changes the statistics in the present and leads to purely quantum features. In particular conditioned average outcomes of a weak measurement, so-called weak values, were shown to go beyond the conventional range, give a way to directly measure complex quantities, and can be used to enhance the sensitivity of quantum meters. Recently, these concepts have been considered in the general case of open quantum systems where decoherence occurs. Then, what are the properties of weak values for the unavoidable measurement associated to decoherence, the one performed by the environment? Here, we answer this question in the simplest open quantum system: a quantum bit in presence of a relaxation channel. We continuously monitor the fluorescence emitted by a superconducting qubit driven at resonance. Conditioned on initial preparation and final single shot measurement outcome of the qubit state, we probe weak values displaying all the above properties. The fluorescence signal exhibits interferences between oscillations associated to past and future quantum states. The measured data are in complete agreement with theory.


title: Power Spectrum and Non-Gaussianities in Anisotropic Inflation

abstract: We study the planar regime of curvature perturbations for single field inflationary models in an axially symmetric Bianchi I background. In a theory with standard scalar field action, the power spectrum for such modes has a pole as the planarity parameter goes to zero. We show that constraints from back reaction lead to a strong lower bound on the planarity parameter for high-momentum planar modes and use this bound to calculate the signal-to-noise ratio of the anisotropic power spectrum in the CMB, which in turn places an upper bound on the Hubble scale during inflation allowed in our model. We find that non-Gaussianities for these planar modes are enhanced for the flattened triangle and the squeezed triangle configurations, but show that the estimated values of the f_NL parameters remain well below the experimental bounds from the CMB for generic planar modes (other, more promising signatures are also discussed). For a standard action, f_NL from the squeezed configuration turns out to be larger compared to that from the flattened triangle configuration in the planar regime. However, in a theory with higher derivative operators, non-Gaussianities from the flattened triangle can become larger than the squeezed configuration in a certain limit of the planarity parameter.


title: Holographic description of quantum black hole on a computer

abstract: The discovery of the fact that black holes radiate particles and eventually evaporate led Hawking to pose the well-known information loss paradox. This paradox caused a long and serious debate since it claims that the fundamental laws of quantum mechanics may be violated. A possible cure appeared recently from superstring theory, a consistent theory of quantum gravity: if the holographic description of a quantum black hole based on the gauge/gravity duality is correct, the information is not lost and quantum mechanics remains valid. Here we test this gauge/gravity duality on a computer at the level of quantum gravity for the first time. The black hole mass obtained by Monte Carlo simulation of the dual gauge theory reproduces precisely the quantum gravity effects in an evaporating black hole. This result opens up totally new perspectives towards quantum gravity since one can simulate quantum black holes through dual gauge theories.


title: A nonlocal kinetic energy functional for an inhomogeneous two-dimensional Fermi gas

abstract: The average-density approximation is used to construct a nonlocal kinetic energy functional for an inhomogeneous two-dimensional Fermi gas. This functional is then used to formulate a Thomas-Fermi von Weizs\"acker-like theory for the description of the ground state properties of the system. The quality of the kinetic energy functional is tested by performing a fully self-consistent calculation for an ideal, harmonically confined, two-dimensional system. Good agreement with exact results are found, with the number and kinetic energy densities exhibiting oscillatory structure associated with the nonlocality of the energy functional. Most importantly, this functional shows a marked improvement over the two-dimensional Thomas-Fermi von Weizs\"acker theory, particularly in the vicinity of the classically forbidden region.


title: Higher-rank instanton cohomology and the quilted Atiyah-Floer conjecture

abstract: Given a closed, connected, oriented 3-manifold with positive first Betti number, one can define an instanton Floer group as well as a quilted Lagrangian Floer group. The quilted Atiyah-Floer conjecture states that these cohomology groups are isomorphic. We initiate a program for proving this conjecture.


title: Diffusion Parameters of Flows in Stable Multi-class Queueing Networks

abstract: We consider open multi-class queueing networks with general arrival processes, general processing time sequences and Bernoulli routing. The network is assumed to be operating under an arbitrary work-conserving scheduling policy that makes the system stable. We study the variability of flows within the network. Computable expressions for quantifying flow variability have previously been discussed in the literature. However, in this paper, we shed more light on such analysis to justify the use of these expressions in the asymptotic analysis of network flows. Towards that end, we find a simple diffusion limit for the inter-class flows and establish the relation to asymptotic (co-)variance rates.


title: On the components of the gauge group for PU(r)-bundles

abstract: We discuss a general procedure for using characteristic classes to study the components of the gauge group for a principal G-bundle. To illustrate this, we work out the case where G is the projective unitary group.


title: Tracking and Characterizing Botnets Using Automatically Generated Domains

abstract: Modern botnets rely on domain-generation algorithms (DGAs) to build resilient command-and-control infrastructures. Recent works focus on recognizing automatically generated domains (AGDs) from DNS traffic, which potentially allows to identify previously unknown AGDs to hinder or disrupt botnets' communication capabilities. The state-of-the-art approaches require to deploy low-level DNS sensors to access data whose collection poses practical and privacy issues, making their adoption problematic. We propose a mechanism that overcomes the above limitations by analyzing DNS traffic data through a combination of linguistic and IP-based features of suspicious domains. In this way, we are able to identify AGD names, characterize their DGAs and isolate logical groups of domains that represent the respective botnets. Moreover, our system enriches these groups with new, previously unknown AGD names, and produce novel knowledge about the evolving behavior of each tracked botnet. We used our system in real-world settings, to help researchers that requested intelligence on suspicious domains and were able to label them as belonging to the correct botnet automatically. Additionally, we ran an evaluation on 1,153,516 domains, including AGDs from both modern (e.g., Bamital) and traditional (e.g., Conficker, Torpig) botnets. Our approach correctly isolated families of AGDs that belonged to distinct DGAs, and set automatically generated from non-automatically generated domains apart in 94.8 percent of the cases.


title: Unitary S Matrices With Long-Range Correlations and the Quantum Black Hole

abstract: We propose an S matrix approach to the quantum black hole in which causality, unitarity and their interrelation play a prominent role. Assuming the 't Hooft S matrix ansatz for a gravitating region surrounded by an asymptotically flat space-time we find a non-local transformation which changes the standard causality requirement but is a symmetry of the unitarity condition of the S matrix. This new S matrix then implies correlations between the in and out states of the theory with the involvement of a third entity which in the case of a quantum black hole, we argue is the horizon S matrix. Such correlations are thus linked to preserving the unitarity of the S matrix and to the fact that entangling unitary operators are nonlocal. The analysis is performed within the Bogoliubov S matrix framework by considering a spacetime consisting of causal complements with a boundary in between. No particular metric or lagrangian dynamics need be invoked even to obtain an evolution equation for the full S matrix. Constraints imposed by the new causality requirement and implications for the effectiveness of field theoretical descriptions and for complementarity are also discussed. We find that the tension between information preservation and complementarity may be resolved provided the full quantum gravity theory either through symmetries or fine tuning forbids the occurrence of closed time like curves of information flow. Then, even if causality is violated near the horizon at any intermediate stage, a standard causal ordering may be preserved for the observer away from the horizon. In the context of the black hole, the novelty of our formulation is that it appears well suited to understand unitarity at any intermediate stage of black hole evaporation. Moreover, it is applicable generally to all theories with long range correlations including the final state projection models.


title: Geometric multigrid for an implicit-time immersed boundary method

abstract: The immersed boundary (IB) method is an approach to fluid-structure interaction that uses Lagrangian variables to describe the structure and Eulerian variables to describe the fluid. Explicit time stepping schemes for the IB method require solvers only for Eulerian equations, for which fast Cartesian grid solution methods are available. Such methods are relatively straightforward to develop and are widely used in practice but often require very small time steps to maintain stability. Implicit-time IB methods permit the stable use of large time steps, but efficient implementations of such methods require significantly more complex solvers that effectively treat both Lagrangian and Eulerian variables simultaneously. Several different approaches to solving the coupled Lagrangian-Eulerian equations have been proposed, but a complete understanding of this problem is still emerging. This paper presents a geometric multigrid method for an implicit-time discretization of the IB equations. This multigrid scheme uses a generalization of box relaxation that is shown to handle problems in which the physical stiffness of the structure is very large. Numerical examples are provided to illustrate the effectiveness and efficiency of the algorithms described herein. These tests show that using multigrid as a preconditioner for a Krylov method yields improvements in both robustness and efficiency as compared to using multigrid as a solver. They also demonstrate that with a time step 100--1000 times larger than that permitted by an explicit IB method, the multigrid-preconditioned implicit IB method is approximately 50--200 times more efficient than the explicit method.