title: On the Minimax Risk of Dictionary Learning

abstract: We consider the problem of learning a dictionary matrix from a number of observed signals, which are assumed to be generated via a linear model with a common underlying dictionary. In particular, we derive lower bounds on the minimum achievable worst case mean squared error (MSE), regardless of computational complexity of the dictionary learning (DL) schemes. By casting DL as a classical (or frequentist) estimation problem, the lower bounds on the worst case MSE are derived by following an established information-theoretic approach to minimax estimation. The main conceptual contribution of this paper is the adaption of the information-theoretic approach to minimax estimation for the DL problem in order to derive lower bounds on the worst case MSE of any DL scheme. We derive three different lower bounds applying to different generative models for the observed signals. The first bound applies to a wide range of models, it only requires the existence of a covariance matrix of the (unknown) underlying coefficient vector. By specializing this bound to the case of sparse coefficient distributions, and assuming the true dictionary satisfies the restricted isometry property, we obtain a lower bound on the worst case MSE of DL schemes in terms of a signal to noise ratio (SNR). The third bound applies to a more restrictive subclass of coefficient distributions by requiring the non-zero coefficients to be Gaussian. While, compared with the previous two bounds, the applicability of this final bound is the most limited it is the tightest of the three bounds in the low SNR regime.


title: Toroidal vortices and the conglomeration of dust into rings in protoplanetary discs

abstract: We identify a new hydrodynamical instability in protoplanetary discs that may arise due to variations in the dust-to-gas ratio and may lead to concentration of dust grains within a disc. The instability can arise due to dust settling, which produces a vertical compositional entropy gradient. The entropy gradient drives a baroclinic instability that is capable of creating toroidal gas vortices that gather dust into rings. Such dust rings are potentially observable via continuum emission of the dust or scattered light. Indeed, this instability may offer an explanation for the rings recently observed in the discs around the young stars HL Tau and TW Hya that does not rely on clearing by protoplanets. The instability may also have wider ramifications, potentially aiding dust agglomeration, altering the radial migration of larger planetesimals, and modifying angular momentum transport within a disc.


title: The Complexity of Non-Iterated Probabilistic Justification Logic

abstract: The logic PJ is a probabilistic logic defined by adding (non-iterated) probability operators to the basic justification logic J. In this paper we establish upper and lower bounds for the complexity of the derivability problem in the logic PJ. The main result of the paper is that the complex- ity of the derivability problem in PJ remains the same as the complexity of the derivability problem in the underlying logic J, which is {\Pi}p2-complete. This implies hat the probability operators do not increase the complex- ity of the logic, although they arguably enrich the expressiveness of the language.


title: A survey of P2P multidimensional indexing structures

abstract: Traditional databases have long since reaped the benefits of multidimensional indexes. Numerous proposals in the literature describe multidimensional index designs for P2P systems. However, none of these designs have had real world implementations. Several proposals for P2P multidimensional indexes are reviewed and analyzed. Znet and VBI-tree are the most promising from a technical standpoint. All of the proposed designs assume honest nodes and are thus open to abuse. This is a critical flaw that must be solved before any of the proposed systems can be used.


title: Role of correlations and off-diagonal terms in binary disordered one dimensional systems

abstract: We investigate one dimensional tight binding model in the presence of a correlated binary disorder. The disorder is due to the interaction of particles with heavy immobile other species. Off-diagonal disorder is created by means of a fast periodic modulation of interspecies interaction. The method based on transfer matrix techniques allows us to calculate the energies of extended modes in the correlated binary disorder. We focus on $N$-mer correlations and regain known results for the case of purely diagonal disorder. For off-diagonal disorder we find resonant energies. We discuss ambiguous properties of those states and compare analytical results with numerical calculations. Separately we describe a special case of the dual random dimer model.


title: Analytic nuclear forces and molecular properties from full configuration interaction quantum Monte Carlo

abstract: Unbiased stochastic sampling of the one- and two-body reduced density matrices is achieved in full configuration interaction quantum Monte Carlo with the introduction of a second, "replica" ensemble of walkers, whose population evolves in imaginary time independently from the first, and which entails only modest additional computational overheads. The matrices obtained from this approach are shown to be representative of full configuration-interaction quality, and hence provide a realistic opportunity to achieve high-quality results for a range of properties whose operators do not necessarily commute with the hamiltonian. A density-matrix formulated quasi-variational energy estimator having been already proposed and investigated, the present work extends the scope of the theory to take in studies of analytic nuclear forces, molecular dipole moments and polarisabilities, with extensive comparison to exact results where possible. These new results confirm the suitability of the sampling technique and, where sufficiently large basis sets are available, achieve close agreement with experimental values, expanding the scope of the method to new areas of investigation.


title: Non-equilibrium spin-crossover in copper phthalocyanine

abstract: We demonstrate the tip induced control of the spin state of copper phthalocyanine (CuPc) on an insulator coated substrate. Accounting for electronic correlations, we find that, under the condition of energetic proximity of neutral excited states to the anionic groundstate, the system can undergo a population inversion towards these excited states. The resulting state of the system is accompanied by a change in the total spin quantum number. Experimental signatures of the crossover are the appearance of additional nodal planes in the topographical STM images as well as a strong suppression of the current near the center of the molecule. The robustness of the effect against moderate charge conserving relaxation processes has also been tested.


title: Kinetic Interpretation of Resonance Phenomena in Low Pressure Capacitively Coupled Radio Frequency Plasmas

abstract: The kinetic origin of resonance phenomena in capacitively coupled radio frequency plasmas is discovered based on particle-based numerical simulations. The analysis of the spatio-temporal distributions of plasma parameters such as the densities of hot and cold electrons, as well as the conduction and displacement currents reveals the mechanism of the formation of multiple electron beams during sheath expansion. The interplay between highly energetic beam electrons and low energetic bulk electrons is identified as the physical origin of the excitation of harmonics in the current.


title: Cyclic codes from the second class two-prime Whiteman's generalized cyclotomic sequence with order 6

abstract: Let $n_1=ef+1$ and $n_2=ef'+1$ be two distinct odd primes with positive integers $e,\ f,\ f'.$ In this paper, the two-prime Whiteman's generalized cyclotomic sequence of order $e=6$ is employed to construct several classes of cyclic codes over $\mathrm{GF}(q)$ with length $n_1n_2$. The lower bounds on the minimum distance of these cyclic codes are obtained.


title: Existence of weak solutions to a class of fourth order partial differential equations with Wasserstein gradient structure

abstract: We prove the global-in-time existence of nonnegative weak solutions to a class of fourth order partial differential equations on a convex bounded domain in arbitrary spatial dimensions. Our proof relies on the formal gradient flow structure of the equation with respect to the $L^2$-Wasserstein distance on the space of probability measures. We construct a weak solution by approximation via the time-discrete minimizing movement scheme; necessary compactness estimates are derived by entropy-dissipation methods. Our theory essentially comprises the thin film and Derrida-Lebowitz-Speer-Spohn equations.


title: A Sparse Stochastic Collocation Technique for High-Frequency Wave Propagation with Uncertainty

abstract: We consider the wave equation with highly oscillatory initial data, where there is uncertainty in the wave speed, initial phase and/or initial amplitude. To estimate quantities of interest related to the solution and their statistics, we combine a high-frequency method based on Gaussian beams with sparse stochastic collocation. Although the wave solution, $u^\varepsilon$, is highly oscillatory in both physical and stochastic spaces, we provide theoretical arguments and numerical evidence that quantities of interest based on local averages of $|u^\varepsilon|^2$ are smooth, with derivatives in the stochastic space uniformly bounded in $\varepsilon$, where $\varepsilon$ denotes the short wavelength. This observable related regularity makes the sparse stochastic collocation approach more efficient than Monte Carlo methods. We present numerical tests that demonstrate this advantage.


title: Uniqueness of dg enhancements for the derived category of a Grothendieck category

abstract: We prove that the derived category of a Grothendieck abelian category has a unique dg enhancement. Under some additional assumptions, we show that the same result holds true for its subcategory of compact objects. As a consequence, we deduce that the unbounded derived category of quasi-coherent sheaves on an algebraic stack and the category of perfect complexes on a noetherian concentrated algebraic stack with quasi-finite affine diagonal and enough perfect coherent sheaves have a unique dg enhancement. In particular, the category of perfect complexes on a noetherian scheme with enough locally free sheaves has a unique dg enhancement.


title: On the existence of a non-relativistic hermitian time operator

abstract: The existence of a hermitian time operator is proposed in the framework of non-relativistic quantum mechanics.The Heisenberg equation of motion is shown to yield constant rate of flow of time.It is shown to yield results consistent with classical expectations,when Bohr's correspondence principle is invoked.Further,it is shown that a massless free particle with non-zero momentum experiences timelessness in its own frame of reference whereas a negatively massive particle perceives forwardness of time.Moreover,application to quantum harmonic oscillators shows instantaneous transition of electrons from one eigenstate to another.


title: The Furstenberg Poisson Boundary and CAT(0) Cube Complexes

abstract: We show under weak hypotheses that $\partial X$, the Roller boundary of a finite dimensional CAT(0) cube complex $X$ is the Furstenberg-Poisson boundary of a sufficiently nice random walk on an acting group $\Gamma$. In particular, we show that if $\Gamma$ admits a nonelementary proper action on $X$, and $\mu$ is a generating probability measure of finite entropy and finite first logarithmic moment, then there is a $\mu$-stationary measure on $\partial X$ making it the Furstenberg-Poisson boundary for the $\mu$-random walk on $\Gamma$. We also show that the support is contained in the closure of the regular points. Regular points exhibit strong contracting properties.


title: Ab initio phase diagram of BaTiO$_3$ under epitaxial strain revisited

abstract: We revisit the phase diagram of BaTiO$_3$ under biaxial strain using a first principles-based effective Hamiltonian approach. We show that, in addition to the tetragonal ($c$), quasi-rhombohedral ($r$), and quasi-orthorhombic ($aa$) ferroelectric phases, that have been discussed previously, there are temperature and strain regions, in particular under tensile strain, where the system decomposes into multi-domain structures. In such cases, the strained system, at least on a local level, recovers the same phase sequence as the unclamped bulk material. Furthermore, we extend these results from the case of "uniform" biaxial strain to the situation where the two in-plane lattice constants are strained differently and show that similar considerations apply in this case.


title: An Improved Decision Procedure for Linear Time Mu-Calculus

abstract: An improved Present Future form (PF form) for linear time $\mu$-calculus ($\nu$TL) is presented in this paper. In particular, the future part of the new version turns into the conjunction of elements in the closure of a formula. We show that every closed $\nu$TL formula can be transformed into the new PF form. Additionally, based on the PF form, an algorithm for constructing Present Future form Graph (PFG), which can be utilized to describe models of a formula, is given. Further, an intuitive and efficient decision procedure for checking satisfiability of the guarded fragment of $\nu$TL formulas based on PFG is proposed and implemented in C++. The new decision procedure has the best time complexity over the existing ones despite the cost of exponential space. Finally, a PFG-based model checking approach for $\nu$TL is discussed where a counterexample can be obtained visually when a model violates a property.


title: JUNO and Neutrinoless Double Beta Decay

abstract: We study the impact of the precision determination of oscillation parameters in the JUNO experiment on half-life predictions for neutrinoless double beta decay. We show that the solar neutrino mixing angle can be measured by JUNO with below 1% uncertainty. This implies in particular that the minimal value of the effective mass in the inverted mass ordering will be known essentially without uncertainty. We demonstrate that this reduces the range of half-life predictions in order to test this value by a factor of two. The remaining uncertainty is caused by nuclear matrix elements. This has important consequences for future double beta decay experiments that aim at ruling out the inverted mass ordering or the Majorana nature of neutrinos.


title: Tunneling conductance in half-metal/conical magnet/superconductor junctions in the adiabatic and non-adiabatic regime: self-consistent calculations

abstract: The tunneling conductance in the half-metal/conical magnet/superconductor (HM/CM/SC) is investigated by the use of the combined Blonder-Tinkham-Klapwijk (BTK) formalism and the Bogoliubov-de Gennes (BdG) equations. We show that the conductance calculated self-consistently differs significantly from the one calculated in the non-self-consistent framework. The use of the self-consistent procedure ensures that the charge conservation is satisfied. Due to the spin band separation in the HM, the conductance in the subgap region is mainly determined by the anomalous Andreev reflection the probability of which strongly depends on the spin transmission in the CM layer. We show that the spin of electron injected from the HM can be transmitted through the CM to the SC adiabatically or non-adiabatically depending on the period of the exchange field modulation. We find that the conductance in the subgap region oscillates as a function of the CM layer thickness wherein the oscillations transform from irregular, in the non-adiabatic regime, to regular in the adiabatic case. In the non-adiabatic regime the decrease of the exchange field amplitude in the CM leads to the emergence of the conductance peak for one particular CM thickness in agreement with experiment [J.W.A Robinson, J. D. S Witt and M. G. Blamire, Science 329, 5987]. For both transport regimes the conductance is analyzed over a broad range of parameters determining the spiral magnetization in the CM.


title: Higher Order Statistics in Switched Diversity Systems

abstract: We analyze the level crossing rate (LCR) and the average fade duration of the output signal-to-noise-ratio (SNR) in generalized switched diversity systems. By using a common approach, we study these higher order statistics for two different kinds of configurations: (1) Colocated diversity, i.e. receiver equipped with multiple antennas, and (2) Distributed diversity, i.e. relaying link with multiple single-antenna threshold-based decode-and-forward (DF) relays. In both cases, we consider the switched diversity combining strategies Selection Combining and Switch \& Stay Combining (SSC). Whenever using threshold-based techniques such as DF or SSC, the output SNR is a discontinuous random process and hence classic Rice approach to calculate the LCR is not applicable. Thus, we use an alternative formulation in terms of the one and two-dimensional cumulative distribution functions of the output SNR. Our results are general, and hold for any arbitrary distribution of fading at the different diversity branches. Moreover, we develop a general asymptotic framework to calculate these higher order statistics in high mean SNR environments which only needs of the univariate probability density function.


title: An entropy based clustering order parameter for finite ensembles of oscillators

abstract: Based on the entropy concept, we define a new clustering order parameter $c$ feasible for finite systems of interacting oscillators. Unlike the generalized synchronization order parameters of the Kuramoto type, this new order parameter singles out the splay state as the unique state with $c = 0$, thus yielding a positive value whenever there is some kind of cluster formation in the system. It is therefore proposed to be monitored alongside the Kuramoto order parameters as a means to quantify the overall amount of clustering in the system.


title: On stochastic conservation laws and Malliavin calculus

abstract: For stochastic conservation laws driven by a semilinear noise term, we propose a generalization of the Kru\v{z}kov entropy condition by allowing the Kru\v{z}kov constants to be Malliavin differentiable random variables. Existence and uniqueness results are provided. Our approach sheds some new light on the stochastic entropy conditions put forth by Feng and Nualart [J. Funct. Anal., 2008] and Bauzet, Vallet, and Wittbold [J. Hyperbolic Differ. Equ., 2012].


title: Gamma ray tests of Minimal Dark Matter

abstract: We reconsider the model of Minimal Dark Matter (a fermionic, hypercharge-less quintuplet of the EW interactions) and compute its gamma ray signatures. We compare them with a number of gamma ray probes: the galactic halo diffuse measurements, the galactic center line searches and recent dwarf galaxies observations. We find that the original minimal model, whose mass is fixed at 9.4 TeV by the relic abundance requirement, is constrained by the line searches from the Galactic Center: it is ruled out if the Milky Way possesses a cuspy profile such as NFW but it is still allowed if it has a cored one. Observations of dwarf spheroidal galaxies are also relevant (in particular searches for lines), and ongoing astrophysical progresses on these systems have the potential to eventually rule out the model. We also explore a wider mass range, which applies to the case in which the relic abundance requirement is relaxed. Most of our results can be safely extended to the larger class of multi-TeV WIMP DM annihilating into massive gauge bosons.


title: Hairy Black Hole Stability in AdS, Quantum Mechanics on the Half-Line and Holography

abstract: We consider the linear stability of $4$-dimensional hairy black holes with mixed boundary conditions in Anti-de Sitter spacetime. We focus on the mass of scalar fields around the maximally supersymmetric vacuum of the gauged $\mathcal{N}=8$ supergravity in four dimensions, $m^{2}=-2l^{-2}$. It is shown that the Schr\"{o}dinger operator on the half-line, governing the $S^{2}$, $H^{2}$ or $\mathbb{R}^{2}$ invariant mode around the hairy black hole, allows for non-trivial self-adjoint extensions and each of them correspons to a class of mixed boundary conditions in the gravitational theory. Discarding the self-adjoint extensions with a negative mode impose a restriction on these boundary conditions. The restriction is given in terms of an integral of the potential in the Schr\"{o}dinger operator resembling the estimate of Simon for Schr\"{o}dinger operators on the real line. In the context of AdS/CFT duality, our result has a natural interpretation in terms of the field theory dual effective potential.


title: An infinite family of excluded minors for strong base-orderability

abstract: We discuss a conjecture of Ingleton on excluded minors for base-orderability, and, extending a result he stated, we prove that infinitely many of the matroids that he identified are excluded minors for base-orderability, as well as for the class of gammoids. We prove that a paving matroid is base-orderable if and only if it has no minor that is isomorphic to the cycle matroid of the complete graph on four vertices. For each k that is at least 2, we define the property of k-base-orderability, which lies strictly between base-orderability and strong base-orderability, and we show that k-base-orderable matroids form what Ingleton called a complete class. By generalizing an example of Ingleton, we construct a set of matroids, each of which is an excluded minor for k-base-orderability, but is (k-1)-base-orderable; the union of these sets, over all k, is an infinite set of base-orderable excluded minors for strong base-orderability.


title: The tail of the crossing probability in near-critical percolation --- an appendix to Ahlberg & Steif [arXiv:1405.7144]

abstract: We answer a question of Ahlberg and Steif (2014) by finding the tail behaviour of the crossing probability in near-critical planar percolation. Interestingly, this superexponentially small behaviour is different from the case of dynamical percolation, where the analogous tail probability was proved to be at least exponential and at most superpolynomial by Hammond, Mossel and Pete (2012). The proof is simple, given the scale covariance established by Garban, Pete and Schramm (2013).


title: How to Generate a Good Word Embedding?

abstract: We analyze three critical components of word embedding training: the model, the corpus, and the training parameters. We systematize existing neural-network-based word embedding algorithms and compare them using the same corpus. We evaluate each word embedding in three ways: analyzing its semantic properties, using it as a feature for supervised tasks and using it to initialize neural networks. We also provide several simple guidelines for training word embeddings. First, we discover that corpus domain is more important than corpus size. We recommend choosing a corpus in a suitable domain for the desired task, after that, using a larger corpus yields better results. Second, we find that faster models provide sufficient performance in most cases, and more complex models can be used if the training corpus is sufficiently large. Third, the early stopping metric for iterating should rely on the development set of the desired task rather than the validation loss of training embedding.


title: Statistics of time delay and scattering correlation functions in chaotic systems I. Random Matrix Theory

abstract: We consider the statistics of time delay in a chaotic cavity having $M$ open channels, in the absence of time-reversal invariance. In the random matrix theory approach, we compute the average value of polynomial functions of the time delay matrix $Q=-i\hbar S^\dag dS/dE$, where $S$ is the scattering matrix. Our results do not assume $M$ to be large. In a companion paper, we develop a semiclassical approximation to $S$-matrix correlation functions, from which the statistics of $Q$ can also be derived. Together, these papers contribute to establishing the conjectured equivalence between the random matrix and the semiclassical approaches.

categories: nlin.CD cond-mat.mes-hall math-ph math.MP


title: Summary of the searches for squarks and gluinos using $\sqrt{s}$ = 8 TeV pp collisions with the ATLAS experiment at the LHC

abstract: A summary is presented of ATLAS searches for gluinos and first- and second-generation squarks in final states containing jets and missing transverse momentum, with or without leptons or b-jets, in the $\sqrt{s}$ = 8 TeV data set collected at the Large Hadron Collider in 2012. This paper reports the results of new interpretations and statistical combinations of previously published analyses, as well as a new analysis. Since no significant excess of events over the Standard Model expectation is observed, the data are used to set limits in a variety of models. In all the considered simplified models that assume R-parity conservation, the limit on the gluino mass exceeds 1150 GeV at 95% confidence level, for an LSP mass smaller than 100 GeV. Furthermore, exclusion limits are set for left-handed squarks in a phenomenological MSSM model, a minimal Supergravity/Constrained MSSM model, R-parity-violation scenarios, a minimal gauge-mediated supersymmetry breaking model, a natural gauge mediation model, a non-universal Higgs mass model with gaugino mediation and a minimal model of universal extra dimensions.


title: Statistics of time delay and scattering correlation functions in chaotic systems II. Semiclassical Approximation

abstract: We consider $S$-matrix correlation functions for a chaotic cavity having $M$ open channels, in the absence of time-reversal invariance. Relying on a semiclassical approximation, we compute the average over $E$ of the quantities ${\rm Tr}[S^\dag(E-\epsilon)S(E+\epsilon)]^n$, for general positive integer $n$. Our result is an infinite series in $\epsilon$, whose coefficients are rational functions of $M$. From this we extract moments of the time delay matrix $Q=-i\hbar S^\dag dS/dE$, and check that the first 8 of them agree with the random matrix theory prediction from our previous paper.


title: Synthesis of Recursive ADT Transformations from Reusable Templates

abstract: Recent work has proposed a promising approach to improving scalability of program synthesis by allowing the user to supply a syntactic template that constrains the space of potential programs. Unfortunately, creating templates often requires nontrivial effort from the user, which impedes the usability of the synthesizer. We present a solution to this problem in the context of recursive transformations on algebraic data-types. Our approach relies on polymorphic synthesis constructs: a small but powerful extension to the language of syntactic templates, which makes it possible to define a program space in a concise and highly reusable manner, while at the same time retains the scalability benefits of conventional templates. This approach enables end-users to reuse predefined templates from a library for a wide variety of problems with little effort. The paper also describes a novel optimization that further improves the performance and scalability of the system. We evaluated the approach on a set of benchmarks that most notably includes desugaring functions for lambda calculus, which force the synthesizer to discover Church encodings for pairs and boolean operations.


title: Anomalies in the GRB spatial distribution

abstract: Swift's remarkable ability to quickly localize gamma-ray bursts has led to the accumulation of a sizable burst sample for which both angular locations and redshifts are measured. This sample has become large enough that it can potentially be used to probe angular anisotropies indicative of large-scale universal structure. In a previous work, a large clustering of gamma-ray bursts at redshift z about 2 was reported in the general direction of the constellations of Hercules and Corona Borealis. Since that report, a 42 per cent increase in the number of z about 2 gamma-ray bursts has been observed, warranting an updated analysis. Surprisingly, the cluster is more pronounced now than it was when it was first reported.


title: Generating Partially Synthetic Geocoded Public Use Data with Decreased Disclosure Risk Using Differential Smoothing

abstract: When collecting geocoded confidential data with the intent to disseminate, agencies often resort to altering the geographies prior to making data publicly available due to data privacy obligations. An alternative to releasing aggregated and/or perturbed data is to release multiply-imputed synthetic data, where sensitive values are replaced with draws from statistical models designed to capture important distributional features in the collected data. One issue that has received relatively little attention, however, is how to handle spatially outlying observations in the collected data, as common spatial models often have a tendency to overfit these observations. The goal of this work is to bring this issue to the forefront and propose a solution, which we refer to as "differential smoothing." After implementing our method on simulated data, highlighting the effectiveness of our approach under various scenarios, we illustrate the framework using data consisting of sale prices of homes in San Francisco.


title: Convergence, stability and robustness of multidimensional opinion dynamics in continuous time

abstract: We analyze a continuous time multidimensional opinion model where agents have heterogeneous but symmetric and compactly supported interaction functions. We consider Filippov solutions of the resulting dynamics and show strong Lyapunov stability of all equilibria in the relative interior of the set of equilibria. We investigate robustness of equilibria when a new agent with arbitrarily small weight is introduced to the system in equilibrium. Assuming the interaction functions to be indicators, we provide a necessary condition and a sufficient condition for robustness of the equilibria. Our necessary condition coincides with the necessary and sufficient condition obtained by Blondel et al. for one dimensional opinions.


title: On the red-giant luminosity bump

abstract: The increase in luminosity as a star evolves on the red-giant branch is interrupted briefly when the hydrogen-burning shell reaches the vicinity of the composition discontinuity left behind from the first convective dredge-up. The non-monotonic variation of luminosity causes an accumulation of stars, known as the `bump', in the distribution of stars in the colour-magnitude diagrams of stellar clusters, which has substantial diagnostic potential. Here I present numerical results on this behaviour and discuss the physical reason for the luminosity variation, with the goal of strengthening the understanding of origin of the phenomenon and hence of its diagnostic potential.


title: Clustering Tree-structured Data on Manifold

abstract: Tree-structured data usually contain both topological and geometrical information, and are necessarily considered on manifold instead of Euclidean space for appropriate data parameterization and analysis. In this study, we propose a novel tree-structured data parameterization, called Topology-Attribute matrix (T-A matrix), so the data clustering task can be conducted on matrix manifold. We incorporate the structure constraints embedded in data into the negative matrix factorization method to determine meta-trees from the T-A matrix, and the signature vector of each single tree can then be extracted by meta-tree decomposition. The meta-tree space turns out to be a cone space, in which we explore the distance metric and implement the clustering algorithm based on the concepts like Fr\'echet mean. Finally, the T-A matrix based clustering (TAMBAC) framework is evaluated and compared using both simulated data and real retinal images to illustrate its efficiency and accuracy.


title: Secure Partial Repair in Wireless Caching Networks with Broadcast Channels

abstract: We study security in partial repair in wireless caching networks where parts of the stored packets in the caching nodes are susceptible to be erased. Let us denote a caching node that has lost parts of its stored packets as a sick caching node and a caching node that has not lost any packet as a healthy caching node. In partial repair, a set of caching nodes (among sick and healthy caching nodes) broadcast information to other sick caching nodes to recover the erased packets. The broadcast information from a caching node is assumed to be received without any error by all other caching nodes. All the sick caching nodes then are able to recover their erased packets, while using the broadcast information and the nonerased packets in their storage as side information. In this setting, if an eavesdropper overhears the broadcast channels, it might obtain some information about the stored file. We thus study secure partial repair in the senses of information-theoretically strong and weak security. In both senses, we investigate the secrecy caching capacity, namely, the maximum amount of information which can be stored in the caching network such that there is no leakage of information during a partial repair process. We then deduce the strong and weak secrecy caching capacities, and also derive the sufficient finite field sizes for achieving the capacities. Finally, we propose optimal secure codes for exact partial repair, in which the recovered packets are exactly the same as erased packets.


title: Melvin Magnetic Fluxtube/Cosmology Correspondence

abstract: We explore a correspondence between Melvin magnetic fluxtubes and anisotropic cosmological solutions, which we call `Melvin cosmologies'. The correspondence via analytic continuation provides useful information in both directions. Solution generating techniques known on the fluxtube side can also be used for generating cosmological backgrounds. Melvin cosmologies interpolate between different limiting Kasner behaviors at early and late times. This has an analogue on the fluxtube side between limiting Levi-Civita behavior at small and large radii. We construct generalized Melvin fluxtubes and cosmologies in both Einstein-Maxwell theory and dilaton gravity and show that similar properties hold.


title: Identification of Stochastic Wiener Systems using Indirect Inference

abstract: We study identification of stochastic Wiener dynamic systems using so-called indirect inference. The main idea is to first fit an auxiliary model to the observed data and then in a second step, often by simulation, fit a more structured model to the estimated auxiliary model. This two-step procedure can be used when the direct maximum-likelihood estimate is difficult or intractable to compute. One such example is the identification of stochastic Wiener systems, i.e.,~linear dynamic systems with process noise where the output is measured using a non-linear sensor with additive measurement noise. It is in principle possible to evaluate the log-likelihood cost function using numerical integration, but the corresponding optimization problem can be quite intricate. This motivates studying consistent, but sub-optimal, identification methods for stochastic Wiener systems. We will consider indirect inference using the best linear approximation as an auxiliary model. We show that the key to obtain a reliable estimate is to use uncertainty weighting when fitting the stochastic Wiener model to the auxiliary model estimate. The main technical contribution of this paper is the corresponding asymptotic variance analysis. A numerical evaluation is presented based on a first-order finite impulse response system with a cubic non-linearity, for which certain illustrative analytic properties are derived.


title: Gamma-rays from Heavy Minimal Dark Matter

abstract: We consider the annihilation into gamma rays of Minimal Dark Matter candidates in the fermionic 5-plet and scalar 7-plet representations of $SU(2)_L$, taking into account both the Sommerfeld effect and the internal bremsstrahlung. Assuming the Einasto profile, we show that present measurements of the Galactic Center by the H.E.S.S. instrument exclude the 5-plet and 7-plet as the dominant form of dark matter for masses between 1 TeV and 20 TeV, in particular, the 5-plet mass leading to the observed dark matter density via thermal freeze-out. We also discuss prospects for the upcoming Cherenkov Telescope Array, which will be able to probe even heavier dark matter masses, including the scenario where the scalar 7-plet is thermally produced.


title: Order-disorder transition in the two dimensional interacting monomer-dimer model: Ising criticality

abstract: We study the order-disorder transition of the two dimensional interacting monomer-dimer model (IMD) which has two symmetric absorbing states. To be self-contained, we first estimate numerically the dynamic exponent $z$ of the two dimensional Ising model. From the relaxation dynamics of the magnetization at the critical point, we obtain $\beta/(\nu z) = 0.057~650(12)$, or $z = 2.168 26(45)$, where $\beta = \frac{1}{8}$ and $\nu=1$ are exactly known exponents. We, then, compare the critical relaxation of the order parameter at the transition point of the IMD with that of the Ising model. We found that the critical relaxation exponent $\beta/(\nu z)$ is in good agreement with the Ising model, unlike the recent claim by Nam et al [JSTAT {\bf (2014)}, P08011]. We also claim that the Binder cumulant is not an efficient quantity to locate the order-disorder transition point of the model with two symmetric absorbing states.


title: Light localisation and magneto-optic enhancement in Ni anti-dot arrays

abstract: The excitation of surface plasmons in magnetic nano-structures can strongly influence their magneto-optical properties. Here, we use photoemission electron microscopy to map the spatial distribution of the electric near-field on a nano-patterned magnetic surface that supports plasmon polaritons. By using different photon energies and polarization states of the incident light we reveal that the electric near-field is either concentrated in spots forming a hexagonal lattice with the same symmetry as the Ni nano-pattern or in stripes oriented along the $\Gamma$-K direction of the lattice and perpendicular to the polarization direction. We show that the polarization-dependent near-field enhancement on the patterned surface is directly correlated to both the excitation of surface plasmon polaritons on the patterned surface as well as the enhancement of the polar magneto-optical Kerr effect.


title: Putting Logic-Based Distributed Systems on Stable Grounds

abstract: In the Declarative Networking paradigm, Datalog-like languages are used to express distributed computations. Whereas recently formal operational semantics for these languages have been developed, a corresponding declarative semantics has been lacking so far. The challenge is to capture precisely the amount of nondeterminism that is inherent to distributed computations due to concurrency, networking delays, and asynchronous communication. This paper shows how a declarative, model-based semantics can be obtained by simply using the well-known stable model semantics for Datalog with negation. We show that the model-based semantics matches previously proposed formal operational semantics.


title: Initial state fluctuations from mid-peripheral to ultra-central collisions in a event-by-event transport approach

abstract: We have developed a relativistic kinetic transport approach that incorporates initial state fluctuations allowing to study the build up of elliptic flow $v_2$ and high order harmonics $v_3$, $v_4$ and $v_5$ for a fluid at fixed $\eta/s(T)$. We study the effect of the $\eta/s$ ratio and its T dependence on the build up of the $v_n(p_T)$ for two different beam energies: RHIC for Au+Au at $\sqrt{s}=200 \,GeV$ and LHC for $Pb+Pb$ at $\sqrt{s}=2.76 \,TeV$. We find that for the two different beam energies considered the suppression of the $v_n(p_T)$ due to the viscosity of the medium have different contributions coming from the cross over or QGP phase. Our study reveals that only in ultra-central collisions ($0 - 0.2 \%$) the $v_n(p_T)$ have a stronger sensitivity to the T dependence of $\eta/s$ in the QGP phase and this sensitivity increases with the order of the harmonic n. Moreover, the study of the correlations between the initial spatial anisotropies $\epsilon_n$ and the final flow coefficients $v_n$ shows that at LHC energies there is more correlation than at RHIC energies. The degree of correlation increases from peripheral to central collisions, but only in ultra-central collisions at LHC, we find that the linear correlation coefficient $C(n,n) \approx 1$ for $n=2,3,4$ and $5$. This suggests that the final correlations in the ($v_n$,$v_m$) space reflect the initial correlations in the ($\epsilon_n$,$\epsilon_m$) space.


title: Maximizing electrical power supply using FACTS devices

abstract: Modern society critically depends on the services electric power provides. Power systems rely on a network of power lines and transformers to deliver power from sources of power (generators) to the consumers (loads). However, when power lines fail (for example, through lightning or natural disasters) or when the system is heavily used, the network is often unable to fulfill all of the demand for power. While systems are vulnerable to these failures, increasingly, sophisticated control devices are being deployed to improve the efficiency of power systems. Such devices can also be used to improve the resiliency of power systems to failures. In this paper, we focus on using FACTS devices in this context. A FACTS device allows power grid operators to adjust the impedance parameters of power lines, thereby redistributing flow in the network and potentially increasing the amount of power that is supplied. Here we develop new approaches for determining the optimal parameter settings for FACTS devices in order to supply the maximal amount of power when networks are stressed, e.g. power line failures and heavy utilization.


title: TMD factorization and evolution at large $b_T$

abstract: In using transverse-momentum-dependent (TMD) parton densities and fragmentation functions, important non-perturbative information is at large transverse position $b_T$. This concerns both the TMD functions and their evolution. Fits to high energy data tend to predict too rapid evolution when extrapolated to low energies where larger values of $b_T$ dominate. I summarize a new analysis of the issues. It results in a proposal for much weaker $b_T$ dependence at large $b_T$ for the evolution kernel, while preserving the accuracy of the existing fits. The results are particularly important for using transverse-spin-dependent functions like the Sivers function.


title: Perturbation theory of non-perturbative QCD

abstract: Perturbation theory is shown to be working in the IR limit of pure SU(3) Yang-Mills theory in Landau gauge by an unconventional setting of the perturbative expansion. A dynamical mass is predicted for the gluon and the lattice data are reproduced fairly well by a second-order expansion, without any free parameter. The effective running coupling is small in the IR and the approximation can be improved by inclusion of higher order terms.


title: Meta-Kernelization using Well-Structured Modulators

abstract: Kernelization investigates exact preprocessing algorithms with performance guarantees. The most prevalent type of parameters used in kernelization is the solution size for optimization problems; however, also structural parameters have been successfully used to obtain polynomial kernels for a wide range of problems. Many of these parameters can be defined as the size of a smallest modulator of the given graph into a fixed graph class (i.e., a set of vertices whose deletion puts the graph into the graph class). Such parameters admit the construction of polynomial kernels even when the solution size is large or not applicable. This work follows up on the research on meta-kernelization frameworks in terms of structural parameters. We develop a class of parameters which are based on a more general view on modulators: instead of size, the parameters employ a combination of rank-width and split decompositions to measure structure inside the modulator. This allows us to lift kernelization results from modulator-size to more general parameters, hence providing smaller kernels. We show (i) how such large but well-structured modulators can be efficiently approximated, (ii) how they can be used to obtain polynomial kernels for any graph problem expressible in Monadic Second Order logic, and (iii) how they allow the extension of previous results in the area of structural meta-kernelization.


title: Explicit cost bounds of stochastic Galerkin approximations for parameterized PDEs with random coefficients

abstract: This work analyzes the overall computational complexity of the stochastic Galerkin finite element method (SGFEM) for approximating the solution of parameterized elliptic partial differential equations with both affine and non-affine random coefficients. To compute the fully discrete solution, such approaches employ a Galerkin projection in both the deterministic and stochastic domains, produced here by a combination of finite elements and a global orthogonal basis, defined on an isotopic total degree index set, respectively. To account for the sparsity of the resulting system, we present a rigorous cost analysis that considers the total number of coupled finite element systems that must be simultaneously solved in the SGFEM. However, to maintain sparsity as the coefficient becomes increasingly nonlinear in the parameterization, it is necessary to also approximate the coefficient by an additional orthogonal expansion. In this case we prove a rigorous complexity estimate for the number of floating point operations (FLOPs) required per matrix-vector multiplication of the coupled system. Based on such complexity estimates we also develop explicit cost bounds in terms of FLOPs to solve the stochastic Galerkin (SG) systems to a prescribed tolerance, which are used to compare with the minimal complexity estimates of a stochastic collocation finite element method (SCFEM), shown in our previous work [16]. Finally, computational evidence complements the theoretical estimates and supports our conclusion that, in the case that the coefficient is affine, the coupled SG system can be solved more efficiently than the decoupled SC systems. However, as the coefficient becomes more nonlinear, it becomes prohibitively expensive to obtain an approximation with the SGFEM.


title: Automatic Identification of Animal Breeds and Species Using Bioacoustics and Artificial Neural Networks

abstract: In this research endeavor, it was hypothesized that the sound produced by animals during their vocalizations can be used as identifiers of the animal breed or species even if they sound the same to unaided human ear. To test this hypothesis, three artificial neural networks (ANNs) were developed using bioacoustics properties as inputs for the respective automatic identification of 13 bird species, eight dog breeds, and 11 frog species. Recorded vocalizations of these animals were collected and processed using several known signal processing techniques to convert the respective sounds into computable bioacoustics values. The converted values of the vocalizations, together with the breed or species identifications, were used to train the ANNs following a ten-fold cross validation technique. Tests show that the respective ANNs can correctly identify 71.43\% of the birds, 94.44\% of the dogs, and 90.91\% of the frogs. This result show that bioacoustics and ANN can be used to automatically determine animal breeds and species, which together could be a promising automated tool for animal identification, biodiversity determination, animal conservation, and other animal welfare efforts.


title: A free boundary problem for fluid flow through porous media

abstract: This brief note addresses the free boundary problem arising from the steady two-dimensional seepage flow through a rectangular dam. The flow problem consists in finding the free boundary location, and the velocity and pressure fields. The problem solution is obtained by using an approximate model which provides the analytical expression for the free boundary profile. Within the range of validity of the model assumptions, the computational results show agreement with data in literature.


title: On additive shifts of multiplicative almost-subgroups in finite fields

abstract: We prove that for sets $A, B, C \subset \mathbb{F}_p$ with $|A|=|B|=|C| \leq \sqrt{p}$ and a fixed $0 \neq d \in \mathbb{F}_p$ holds $$ \max(|AB|, |(A+d)C|) \gg|A|^{1+1/26}. $$ In particular, $$ |A(A+1)| \gg |A|^{1 + 1/26} $$ and $$ \max(|AA|, |(A+1)(A+1)|) \gg |A|^{1 + 1/26}. $$ The first estimate improves the bound by Roche-Newton and Jones. In the general case of a field of order $q = p^m$ we obtain similar estimates with the exponent $1+1/559 + o(1)$ under the condition that $AB$ does not have large intersection with any subfield coset, answering a question of Shparlinski. Finally, we prove the estimate $$ \left| \sum_{x \in \mathbb{F}_q} \psi(x^n) \right| \ll q^{\frac{7 - 2\delta_2}{8}}n^{\frac{2+2\delta_2}{8}} $$ for Gauss sums over $\mathbb{F}_q$, where $\psi$ is a non-trivial additive character and $\delta_2 = 1/56 + o(1)$. The estimate gives an improvement over the classical Weil bound when $q^{1/2} \ll n = o\left( q^{29/57 + o(1)} \right)$.


title: Inequalities for operator space numerical radius of $2\times 2$ block matrices

abstract: In this paper, we study the relationship between operator space norm and operator space numerical radius on the matrix space $\mathcal{M}_n(X)$, when $X$ is a numerical radius operator space. Moreover, we establish several inequalities for operator space numerical radius and the maximal numerical radius norm of $2\times 2$ operator matrices and their off-diagonal parts. One of our main results states that if $(X, (O_n))$ is an operator space, then \begin{align*} \frac12\max\big(W_{\max}(x_1+x_2)&, W_{\max}(x_1-x_2) \big)\\ &\le W_{\max}\Big(\begin{bmatrix} 0 & x_1 \\ x_2 & 0 \end{bmatrix}\Big)\\ &\hspace{1.5cm}\le \frac12\left(W_{\max}(x_1+x_2)+ W_{\max}(x_1-x_2) \right) \end{align*} for all $x_1, x_2\in \mathcal{M}_n(X)$.


title: Bulk Local Operators, Conformal Descendants and Radial Quantization

abstract: We establish a construction of the bulk local operators in AdS by considering CFT at finite energy scale. Without assuming any prior knowledge about the bulk, the normalizable solution to the bulk free field equation automatically appears in the field theory arguments. In the radial quantization formalism, we find a properly regularized version of our initial construction. Possible generalizations beyond pure AdS are also discussed.


title: CMB lensing tomography with the DES Science Verification galaxies

abstract: We measure the cross-correlation between the galaxy density in the Dark Energy Survey (DES) Science Verification data and the lensing of the cosmic microwave background (CMB) as reconstructed with the Planck satellite and the South Pole Telescope (SPT). When using the DES main galaxy sample over the full redshift range $0.2 < z < 1.2$, a cross-correlation signal is detected at $6 \sigma$ and $4\sigma$ with SPT and Planck respectively. We then divide the DES galaxies into five photometric redshift bins, finding significant ($>$$2 \sigma$) detections in all bins. Comparing to the fiducial Planck cosmology, we find the redshift evolution of the signal matches expectations, although the amplitude is consistently lower than predicted across redshift bins. We test for possible systematics that could affect our result and find no evidence for significant contamination. Finally, we demonstrate how these measurements can be used to constrain the growth of structure across cosmic time. We find the data are fit by a model in which the amplitude of structure in the $z<1.2$ universe is $0.73 \pm 0.16$ times as large as predicted in the LCDM Planck cosmology, a $1.7\sigma$ deviation.


title: Cosmology from Cosmic Shear with DES Science Verification Data

abstract: We present the first constraints on cosmology from the Dark Energy Survey (DES), using weak lensing measurements from the preliminary Science Verification (SV) data. We use 139 square degrees of SV data, which is less than 3\% of the full DES survey area. Using cosmic shear 2-point measurements over three redshift bins we find $\sigma_8 (\Omega_{\rm m}/0.3)^{0.5} = 0.81 \pm 0.06$ (68\% confidence), after marginalising over 7 systematics parameters and 3 other cosmological parameters. We examine the robustness of our results to the choice of data vector and systematics assumed, and find them to be stable. About $20$\% of our error bar comes from marginalising over shear and photometric redshift calibration uncertainties. The current state-of-the-art cosmic shear measurements from CFHTLenS are mildly discrepant with the cosmological constraints from Planck CMB data; our results are consistent with both datasets. Our uncertainties are $\sim$30\% larger than those from CFHTLenS when we carry out a comparable analysis of the two datasets, which we attribute largely to the lower number density of our shear catalogue. We investigate constraints on dark energy and find that, with this small fraction of the full survey, the DES SV constraints make negligible impact on the Planck constraints. The moderate disagreement between the CFHTLenS and Planck values of $\sigma_8 (\Omega_{\rm m}/0.3)^{0.5}$ is present regardless of the value of $w$.


title: BPS black holes in a non-homogeneous deformation of the stu model of $N=2$, $D=4$ gauged supergravity

abstract: We consider a deformation of the well-known stu model of $N=2$, $D=4$ supergravity, characterized by a non-homogeneous special K\"{a}hler manifold, and by the smallest electric-magnetic duality Lie algebra consistent with its upliftability to five dimensions. We explicitly solve the BPS attractor equations and construct static supersymmetric black holes with radial symmetry, in the context of $\text{U}(1)$ dyonic Fayet-Iliopoulos gauging, focussing on axion-free solutions. Due to non-homogeneity of the scalar manifold, the model evades the analysis recently given in the literature. The relevant physical properties of the resulting black hole solution are discussed.


title: Sub-Riemannian distance on the Lie group $SO_0(2,1)$

abstract: A left-invariant sub-Riemannian metric $d$ on the shortened Lorentz group $SO_0(2,1)$ under the condition that $d$ is right-invariant relative to the orthogonal Lie subgroup $1\otimes SO(2)$ is studied. The distance between arbitrary two elements, the cut locus (as the union of the subgroup $1\otimes SO(2)$ with the antipodal set to the submanifold of symmetric matrices in the open solid torus $SO_0(2,1)$), and the conjugate set for the unit are found for $(SO_0(2,1),d)$.


title: Asymptotic scaling and continuum limit of pure SU(3) lattice gauge theory

abstract: Recently the Yang-Mills gradient flow of pure SU(3) lattice gauge theory has been calculated in the range from $\beta=6/g_0^2=6.3$ to~7.5 (Asakawa et al.), where $g_0^2$ is the bare coupling constant of the SU(3) Wilson action. Estimates of the deconfining phase transition are available from $\beta =5.7$ to~6.8 (Francis et al.). Here it is shown that the entire range from 5.7 to 7.5 is well described by a power series of the lattice spacing $a$ times the lambda lattice mass scale $\Lambda_L$, using asymptotic scaling in the 2-loop and 3-loop approximations for $a\Lambda_L$. In both cases identical ratios for gradient flows versus deconfinement observables are obtained. Differences in the normalization constants with respect to $\Lambda_L$ give a handle on their systematic errors.


title: Calculation of STOs electron repulsion integrals by ellipsoidal expansion and large-order approximations

abstract: For general two-electron two-centre integrals over Slater-type orbitals (STOs), the use of the Neumann expansion for the Coulomb interaction potential yields infinite series in terms of few basic functions. In many important cases the number of terms necessary to achieve convergence by a straightforward summation is large and one is forced to calculate the basic integrals of high order. We present a systematic approach to calculation of the higher-order terms in the Neumann series by large-order expansions of the basic integrals. The final expressions are shown to be transparent and straightforward to implement, and all auxiliary quantities can be calculated analytically. Moreover, numerical stability and computational efficiency are also discussed. Results of the present work can be used to speed up calculations of the STOs integral files, but also to study convergence of the Neumann expansion and develop appropriate convergence accelerators.


title: Effects of strong magnetic fields and rotation on white dwarf structure

abstract: In this work we compute models for relativistic white dwarfs in the presence of strong magnetic fields. These models possibly contribute to super-luminous SNIa. With an assumed axi-symmetric and poloidal magnetic field, we study the possibility of existence of super-Chandrasekhar magnetized white dwarfs by solving numerically the Einstein-Maxwell equations, by means of a pseudo-spectral method. We obtain a self-consistent rotating and non-rotating magnetized white dwarf models. According to our results, a maximum mass for a static magnetized white dwarf is 2.13 $\rm{M_{\odot}}$ in the Newtonian case and 2.09 $\rm{M_{\odot}}$ while taking into account general relativistic effects. Furthermore, we present results for rotating magnetized white dwarfs. The maximum magnetic field strength reached at the center of white dwarfs is of the order of $10^{15}\,$G in the static case, whereas for magnetized white dwarfs, rotating with the Keplerian angular velocity, is of the order of $10^{14}\,$G.


title: Integrated AlGaAs source of highly indistinguishable and energy-time entangled photons

abstract: The generation of nonclassical states of light in miniature chips is a crucial step towards practical implementations of future quantum technologies. Semiconductor materials are ideal to achieve extremely compact and massively parallel systems and several platforms are currently under development. In this context, spontaneous parametric down conversion in AlGaAs devices combines the advantages of room temperature operation, possibility of electrical injection and emission in the telecom band. Here we report on a chip-based AlGaAs source, producing indistinguishable and energy-time entangled photons with a brightness of $7.2\times10^6$ pairs/s and a signal-to-noise ratio of $141\pm12$. Indistinguishability between the photons is demonstrated via a Hong-Ou-Mandel experiment with a visibility of $89\pm3\%$, while energy-time entanglement is tested via a Franson interferometer leading to a value for the Bell parameter $ S=2.70\pm0.10$.


title: Lambda Bootis stars in the SuperWASP survey

abstract: We have analysed around 170 000 individual photometric WASP measurements of 15 well established lambda Bootis stars to search for variability. The lambda Bootis stars are a small group of late-B to early-F, Pop. I, stars that show moderate to extreme (surface) underabundances (up to a factor 100) of most Fe-peak elements, but solar abundances of lighter elements (C, N, O and S). They are excellent laboratories for the study of fundamental astrophysical processes such as diffusion, meridional circulation, stellar winds, and accretion in the presence of pulsation. From the 15 targets, eight are variable and seven are apparently constant with upper limits between 0.8 and 3.0 mmag. We present a detailed time series analysis and a comparison with previously published results. From an asteroseismologic study we conclude that the found chemical peculiarities are most probably restricted to the surface.


title: Buchdahl's inequality in five dimensional Gauss-Bonnet gravity

abstract: The Buchdahl limit for static spherically symmetric isotropic stars is generalised to the case of five dimensional Gauss-Bonnet gravity. Our result depends on the sign of the Gauss-Bonnet coupling constant $\alpha$. When $\alpha>0$, we find, unlike in general relativity, that the bound is dependent on the stellar structure, in particular the central energy density. We find that stable stellar structures can exist arbitrarily close to the event horizon. Thus stable stars can exist with extra mass in this theory compared to five dimensional general relativity. For $\alpha<0$ it is found that the Buchdahl bound is more restrictive than the general relativistic case.


title: Quantum Monte Carlo calculations of neutron matter with chiral three-body forces

abstract: Chiral effective field theory (EFT) enables a systematic description of low-energy hadronic interactions with controlled theoretical uncertainties. For strongly interacting systems, quantum Monte Carlo (QMC) methods provide some of the most accurate solutions, but they require as input local potentials. We have recently constructed local chiral nucleon-nucleon (NN) interactions up to next-to-next-to-leading order (N$^2$LO). Chiral EFT naturally predicts consistent many-body forces. In this paper, we consider the leading chiral three-nucleon (3N) interactions in local form. These are included in auxiliary field diffusion Monte Carlo (AFDMC) simulations. We present results for the equation of state of neutron matter and for the energies and radii of neutron drops. In particular, we study the regulator dependence at the Hartree-Fock level and in AFDMC and find that present local regulators lead to less repulsion from 3N forces compared to the usual nonlocal regulators.


title: Rigid Unit Modes in $sp$-$sp^2$ Hybridized Carbon Systems: Origin of Negative Thermal Expansion

abstract: Using density functional theory combined with quasi-harmonic approximation, we investigate the thermal expansion behaviors of three different types ({\alpha}, {\beta}, and {\gamma}) of graphyne which is a two-dimensional carbon allotrope composed of sp and $sp^2$ bonds. For each type of graphyne, we obtain the temperature dependent area variation by minimizing its free energy as a function of temperature, which is calculated by considering all the phonon modes in the whole Brillouin zone. We find that all three types of graphyne exhibit negative in-plane thermal expansion up to $T\lesssim1000$ K. The observed in-plane thermal contraction can be attributed partially to the ripple effect, similarly in graphene. The ripple effect itself, however, is not sufficient to explain anomalously larger thermal contraction found in graphyne than in graphene. Our deliberate analysis on the phonon modes observed in graphyne enables us to discover another source causing such thermal expansion anomaly. We find that there are particular phonon modes with frequencies around a few hundreds of cm$^{-1}$ existing exclusively in graphyne that may fill empty spaces resulting in area reduction. These modes are identified as "rigid unit modes" corresponding to the libration of each rigid unit composed of $sp^2$ bonds.


title: De Finetti theorems for a Boolean analogue of easy quantum groups

abstract: We show an organized form of quantum de Finetti theorem for Boolean independence. We define a Boolean analogue of easy quantum groups for the categories of interval partitions, which is a family of sequences of quantum semigroups. We construct the Haar states on those quantum semigroups. The proof of our de Finetti theorem is based on the analysis of the Haar states. [Modified]Definition of the Boolean quantum semigroups on categories of interval partitions [Delete]Classification of categories of interval partitions [Add]Proof of the positiveness of the Haar functionals (in particular they are Haar states)


title: Minimum X-ray source size for a lamp-post corona in light-bending models for AGN

abstract: The lamppost model is often used to describe the X-ray source geometry in AGN, where an infinitesimal point source is located on the black hole spin axis. This is especially invoked for Narrow Line Seyfert 1 (NLS1) galaxies, where an extremely broad iron line seen in episodes of low X-ray flux can both be explained by extremely strong relativistic effects as the source approaches the black hole horizon. The most extreme spectrum seen from the NLS1 1H0707-495 requires that the source is less than 1Rg above the event horizon in this geometry. However, the source must also be large enough to intercept sufficient seed photons from the disk to make the hard X-ray Compton continuum which produces the observed iron line/reflected spectrum. We use a fully relativistic ray tracing code to show that this implies that the source must be substantially larger than 1Rg in 1H0707-495 if the disk is the source of seed photons. Hence the source cannot fit as close as 1Rg to the horizon, so the observed spectrum and variability are not formed purely by effects of strong gravity but probably also by changes in corona and inner accretion flow geometry.


title: Constraints on modified Gauss-Bonnet gravity during big bang nucleosynthesis

abstract: The modified gravity is considered to be one of possible explanations of the accelerated expansions of the present and the early universe. We study effects of the modified gravity on big bang nucleosynthesis (BBN). If effects of the modified gravity are significant during the BBN epoch, they should be observed as changes of primordial light element abundances. We assume a $f(G)$ term with the Gauss-Bonnet term $G$, during the BBN epoch. A power-law relation of $df/dG \propto t^p$ where $t$ is the cosmic time was assumed for the function $f(G)$ as an example case. We solve time evolutions of physical variables during BBN in the $f(G)$ gravity model numerically, and analyzed calculated results. It is found that a proper solution for the cosmic expansion rate can be lost in some parameter region. In addition, we show that calculated results of primordial light element abundances can be significantly different from observational data. Especially, observational limits on primordial D abundance leads to the strongest constraint on the $f(G)$ gravity. We then derive constraints on parameters of the $f(G)$ gravity taking into account the existence of the solution of expansion rate and final light element abundances.


title: Incompressible Euler Equations and the Effect of Changes at a Distance

abstract: Because pressure is determined globally for the incompressible Euler equations, a localized change to the initial velocity will have an immediate effect throughout space. For solutions to be physically meaningful, one would expect such effects to decrease with distance from the localized change, giving the solutions a type of stability. Indeed, this is the case for solutions having spatial decay, as can be easily shown. We consider the more difficult case of solutions lacking spatial decay, and show that such stability still holds, albeit in a somewhat weaker form.


title: On quasi-periodicity properties of fractional integrals and fractional derivatives of periodic functions

abstract: This paper is devoted to the study of quasi-periodic properties of fractional order integrals and derivatives of periodic functions. Considering Riemann-Liouville and Caputo definitions, we discuss when the fractional derivative and when the fractional integral of a certain class of periodic functions satisfies particular properties. We study concepts close to the well known idea of periodic function, such as S-asymptotically periodic, asymptotically periodic or almost periodic function. Boundedness of fractional derivative and fractional integral of a periodic function is also studied.


title: A remark on observability of the wave equation with moving boundary

abstract: We deal with the wave equation with assigned moving boundary ($0<x<a(t)$) upon which Dirichlet or mixed boundary conditions are specified, here $a(t)$ is assumed to move slower than the light and periodically. Moreover $a$ is continuous, piecewise linear with two independent parameters. Our major concern will be an observation problem which is based measuring, at each $t>0$ of the transverse velocity at $a(t)$. The key to the results is the use of a reduction theorem.


title: Fluid dynamic propagation of initial baryon number perturbations on a Bjorken flow background

abstract: Baryon number density perturbations offer a possible route to experimentally measure baryon number susceptibilities and heat conductivity of the quark gluon plasma. We study the fluid dynamical evolution of local and event-by-event fluctuations of baryon number density, flow velocity and energy density on top of a (generalized) Bjorken expansion. To that end we use a background-fluctuation splitting and a Bessel-Fourier decomposition for the fluctuating part of the fluid dynamical fields with respect to the azimuthal angle, the radius in the transverse plane and rapidity. We examine how the time evolution of linear perturbations depends on the equation of state as well as on shear viscosity, bulk viscosity and heat conductivity for modes with different azimuthal, radial and rapidity wave numbers. Finally we discuss how this information is accessible to experiments in terms of the transverse and rapidity dependence of correlation functions for baryonic particles in high energy nuclear collisions.


title: Two weight $L^p$ estimates for paraproducts in non-homogeneous settings

abstract: We give a necessary and sufficient condition for the two weight $L^p$-estimates for paraproducts in non-homogeneous settings, $1<p<\infty$. We are mainly interested in the case $p\ne 2$, since the case $p=2$ is a well-known and easy corollary of the Carleson embedding theorem. The necessary and sufficient condition is given in terms of testing conditions of Sawyer type: for $p\le 2$ only one ("direct") testing condition is required, but for $p>2$ both "direct" and "adjoint" testing conditions are needed. An interesting feature is that the "adjoint" testing condition is that it is a testing condition not for the adjoint of the paraproduct, but for some auxiliary operator.


title: Exact evaluations of some Meijer G-functions and probability of all eigenvalues real for product of two Gaussian matrices

abstract: We provide a proof to a recent conjecture by Forrester [2014 J. Phys. A: Math. Theor. 47, 065202] regarding the algebraic and arithmetic structure of Meijer G-functions which appear in the expression for probability of all eigenvalues real for product of two real Gaussian matrices. In the process we come across several interesting identities involving Meijer G-functions.


title: Detailed analysis of the effects of stencil spatial variations with arbitrary high-order finite-difference Maxwell solver

abstract: Due to discretization effects and truncation to finite domains, many electromagnetic simulations present non-physical modifications of Maxwell's equations in space that may generate spurious signals affecting the overall accuracy of the result. Such modifications for instance occur when Perfectly Matched Layers (PMLs) are used at simulation domain boundaries to simulate open media. Another example is the use of arbitrary order Maxwell solver with domain decomposition technique that may under some condition involve stencil truncations at subdomain boundaries, resulting in small spurious errors that do eventually build up. In each case, a careful evaluation of the characteristics and magnitude of the errors resulting from these approximations, and their impact at any frequency and angle, requires detailed analytical and numerical studies. To this end, we present a general analytical approach that enables the evaluation of numerical discretization errors of fully three-dimensional arbitrary order finite-difference Maxwell solver, with arbitrary modification of the local stencil in the simulation domain. The analytical model is validated against simulations of domain decomposition technique and PMLs, when these are used with very high-order Maxwell solver, as well as in the infinite order limit of pseudo-spectral solvers. Results confirm that the new analytical approach enables exact predictions in each case. It also confirms that the domain decomposition technique can be used with very high-order Maxwell solver and a reasonably low number of guard cells with negligible effects on the whole accuracy of the simulation.


title: Similarity structures and de Rham decomposition

abstract: A similarity structure on a connected manifold M is a Riemannian metric on its universal cover such that the fundamental group of M acts by similarities. If the manifold M is compact, we show that the universal cover admits a de Rham decomposition with at most two factors, one of which is Euclidean. Very recently, after Belgun and Moroianu conjectured that the number of factors was at most one, Matveev and Nikolayevsky found an example with two factors. When the non-flat factor has dimension 2, we give a complete classification of the examples with two factors. In greater dimensions, we make the first steps towards such a classification by showing that M is a fibration (with singularities) by flat Riemannian manifolds; up to a finite covering of M, we may assume that these manifolds are flat tori. We also prove a version of the de Rham decomposition theorem for the universal covers of manifolds with locally metric connections. During the proof, we define a notion of transverse (not necessarily flat) similarity structure on foliations, and show that foliations endowed with such a structure are either transversally flat or transversally Riemannian. None of these results assumes analyticity.


title: A Combined Spitzer and Herschel Infrared Study of Gas and Dust in the Circumbinary Disk Orbiting V4046 Sgr

abstract: We present results from a spectroscopic Spitzer and Herschel mid-to-far-infrared study of the circumbinary disk orbiting the evolved (age ~12-23 Myr) close binary T Tauri system V4046 Sgr. Spitzer IRS spectra show emission lines of [Ne II], H_2 S(1), CO_2 and HCN, while Herschel PACS and SPIRE spectra reveal emission from [O I], OH, and tentative detections of H_2O and high-J transitions of CO. We measure [Ne III]/[Ne II] < 0.13, which is comparable to other X-ray/EUV luminous T Tauri stars that lack jets. We use the H_2 S(1) line luminosity to estimate the gas mass in the relatively warm surface layers of the inner disk. The presence of [O I] emission suggests that CO, H_2O, and/or OH is being photodissociated, and the lack of [C I] emission suggests any excess C may be locked up in HCN, CN and other organic molecules. Modeling of silicate dust grain emission features in the mid-infrared indicates that the inner disk is composed mainly of large (r~5 um) amorphous pyroxene and olivine grains (~86% by mass) with a relatively large proportion of crystalline silicates. These results are consistent with other lines of evidence indicating that planet building is ongoing in regions of the disk within ~30 AU of the central, close binary.


title: Untangling the Recombination Line Emission from HII Regions with Multiple Velocity Components

abstract: HII regions are the ionized spheres surrounding high-mass stars. They are ideal targets for tracing Galactic structure because they are predominantly found in spiral arms and have high luminosities at infrared and radio wavelengths. In the Green Bank Telescope HII Region Discovery Survey (GBT HRDS) we found that >30% of first Galactic quadrant HII regions have multiple hydrogen radio recombination line (RRL) velocities, which makes determining their Galactic locations and physical properties impossible. Here we make additional GBT RRL observations to determine the discrete HII region velocity for all 117 multiple-velocity sources within 18deg. < l < 65deg. The multiple-velocity sources are concentrated in the zone 22deg. < l < 32deg., coinciding with the largest regions of massive star formation, which implies that the diffuse emission is caused by leaked ionizing photons. We combine our observations with analyses of the electron temperature, molecular gas, and carbon recombination lines to determine the source velocities for 103 discrete H II regions (88% of the sample). With the source velocities known, we resolve the kinematic distance ambiguity for 47 regions, and thus determine their heliocentric distances.


title: Joint Optimal Number of RF chains and Power Allocation for Downlink Massive MIMO Systems

abstract: This paper investigates the downlink of massive multiple-input multiple-output (MIMO) systems that include a single cell Base Station (BS) equipped with large number of antennas serving multiple users. As the number of RF chains is getting large, the system model considered in this paper assumes a non negligible circuit power consumption. Hence, the aim of this work is to find the optimal balance between the power consumed by the RF chains and the transmitted power. First, assuming an equal power allocation among users, the optimal number of RF chains to be activated is analytically found. Then, for a given number of RF chains we derive analytically the optimal power allocation among users. Based on these analysis, we propose an iterative algorithm that computes jointly the optimal number of RF chains and the optimal power allocation vector. Simulations validate the analytical results and show the high performance provided by the proposed algorithm.


title: Computing transition rates for the 1-D stochastic Ginzburg--Landau--Allen--Cahn equation for finite-amplitude noise with a rare event algorithm

abstract: In this paper we compute and analyse the transition rates and duration of reactive trajectories of the stochastic 1-D Allen-Cahn equations for both the Freidlin-Wentzell regime (weak noise or temperature limit) and finite-amplitude white noise, as well as for small and large domain. We demonstrate that extremely rare reactive trajectories corresponding to direct transitions between two metastable states are efficiently computed using an algorithm called adaptive multilevel splitting. This algorithm is dedicated to the computation of rare events and is able to provide ensembles of reactive trajectories in a very efficient way. In the small noise limit, our numerical results are in agreement with large-deviation predictions such as instanton-like solutions, mean first passages and escape probabilities. We show that the duration of reactive trajectories follows a Gumbel distribution like for one degree of freedom systems. Moreover, the mean duration growths logarithmically with the inverse temperature. The prefactor given by the potential curvature grows exponentially with size. The main novelty of our work is that we also perform an analysis of reactive trajectories for large noises and large domains. In this case, we show that the position of the reactive front is essentially a random walk. This time, the mean duration grows linearly with the inverse temperature and quadratically with the size. Using a phenomenological description of the system, we are able to calculate the transition rate, although the dynamics is described by neither Freidlin--Wentzell or Eyring--Kramers type of results. Numerical results confirm our analysis.

categories: physics.flu-dyn cond-mat.stat-mech math.NA physics.class-ph


title: Subspace Alignment Based Domain Adaptation for RCNN Detector

abstract: In this paper, we propose subspace alignment based domain adaptation of the state of the art RCNN based object detector. The aim is to be able to achieve high quality object detection in novel, real world target scenarios without requiring labels from the target domain. While, unsupervised domain adaptation has been studied in the case of object classification, for object detection it has been relatively unexplored. In subspace based domain adaptation for objects, we need access to source and target subspaces for the bounding box features. The absence of supervision (labels and bounding boxes are absent) makes the task challenging. In this paper, we show that we can still adapt sub- spaces that are localized to the object by obtaining detections from the RCNN detector trained on source and applied on target. Then we form localized subspaces from the detections and show that subspace alignment based adaptation between these subspaces yields improved object detection. This evaluation is done by considering challenging real world datasets of PASCAL VOC as source and validation set of Microsoft COCO dataset as target for various categories.


title: Stochastic regularization effects of semi-martingales on random functions

abstract: In this paper we address an open question formulated in [17]. That is, we extend the It{\^o}-Tanaka trick, which links the time-average of a deterministic function f depending on a stochastic process X and F the solution of the Fokker-Planck equation associated to X, to random mappings f. To this end we provide new results on a class of adpated and non-adapted Fokker-Planck SPDEs and BSPDEs.


title: Absence of a four-body Efimov effect in the 2 + 2 fermionic problem

abstract: In the free three-dimensional space, we consider a pair of identical $\uparrow$ fermions of some species or in some internal state, and a pair of identical $\downarrow$ fermions of another species or in another state. There is a resonant $s$-wave interaction (that is of zero range and infinite scattering length) between fermions in different pairs, and no interaction within the same pair. We study whether this $2+2$ fermionic system can exhibit (as the $3+1$ fermionic system) a four-body Efimov effect in the absence of three-body Efimov effect, that is the mass ratio $\alpha$ between $\uparrow$ and $\downarrow$ fermions and its inverse are both smaller than 13.6069{\ldots}. For this purpose, we investigate scale invariant zero-energy solutions of the four-body Schr\''odinger equation, that is positively homogeneous functions of the coordinates of degree {$s-7/2$}, where $s$ is a generalized Efimov exponent {that becomes purely imaginary in the presence of a four-body Efimov effect.} Using rotational invariance in momentum space, it is found that the allowed values of $s$ are such that $M(s)$ has a zero eigenvalue; here the operator $M(s)$, that depends on the total angular momentum $\ell$, acts on functions of two real variables (the cosine of the angle between two wave vectors and the logarithm of the ratio of their moduli), and we write it explicitly in terms of an integral matrix kernel. We have performed a spectral analysis of $M(s)$, analytical and for an arbitrary imaginary $s$ for the continuous spectrum, numerical and limited to $s = 0$ and $\ell \le 12$ for the discrete spectrum. We conclude that no eigenvalue of $M(0)$ crosses zero over the mass ratio interval $\alpha \in [1, 13.6069\ldots]$, even if, in the parity sector $(-1)^{\ell}$, the continuous spectrum of $M(s)$ has everywhere a zero lower border. As a consequence, there is no possibility of a four-body Efimov effect for the 2+2 fermions. We also enunciated a conjecture for the fourth virial coefficient of the unitary spin-$1/2$ Fermi gas,inspired from the known analytical form of the third cluster coefficient and involving the integral over the imaginary $s$-axis of $s$ times the logarithmic derivative of the determinant of $M(s)$ summed over all angular momenta.The conjectured value is in contradiction with the experimental results.


title: On Automated Lemma Generation for Separation Logic with Inductive Definitions

abstract: Separation Logic with inductive definitions is a well-known approach for deductive verification of programs that manipulate dynamic data structures. Deciding verification conditions in this context is usually based on user-provided lemmas relating the inductive definitions. We propose a novel approach for generating these lemmas automatically which is based on simple syntactic criteria and deterministic strategies for applying them. Our approach focuses on iterative programs, although it can be applied to recursive programs as well, and specifications that describe not only the shape of the data structures, but also their content or their size. Empirically, we find that our approach is powerful enough to deal with sophisticated benchmarks, e.g., iterative procedures for searching, inserting, or deleting elements in sorted lists, binary search tress, red-black trees, and AVL trees, in a very efficient way.


title: Flatness for a Strongly Degenerate 1-D Parabolic Equation

abstract: We consider the degenerate equation $$\partial\_t f(t,x) - \partial\_x \left( x^{\alpha} \partial\_x f \right)(t,x) =0,$$ on the unit interval $x\in(0,1)$, in the strongly degenerate case $\alpha \in [1,2)$ with adapted boundary conditions at $x=0$ and boundary control at $x=1$. We use the flatness approach to construct explicit controls in some Gevrey classes steering the solution from any initial datum $f\_0 \in L^2(0,1)$ to zero in any time $T\textgreater{}0$.


title: Parity Biquandle Invariants of Virtual Knots

abstract: We define counting and cocycle enhancement invariants of virtual knots using parity biquandles. These invariants are determined by pairs consisting of a biquandle 2-cocycle \phi^0 and a map \phi^1 with certain compatibility conditions leading to one-variable or two-variable polynomial invariants of virtual knots. We provide examples to show that the parity cocycle invariants can distinguish virtual knots which are not distinguished by the corresponding non-parity invariants.


title: Gravitino Dark Matter and low-scale Baryogenesis

abstract: A very simple way to obtain comparable baryon and DM densities in the early Universe is through their contemporary production from the out-of-equilibrium decay of a mother particle, if both populations are suppressed by comparably small numbers, i.e. the CP violation in the decay and the branching fraction respectively. We present a detailed study of this kind of scenario in the context of a R-parity violating realization of the MSSM in which the baryon asymmetry and the gravitino Dark Matter are produced by the decay of a Bino. The implementation of this simple picture in a realistic particle framework results, however, quite involving, due to the non trivial determination of the abundance of the decaying Bino, as well as due to the impact of wash-out processes and of additional sources both for the baryon asymmetry and the DM relic density. In order to achieve a quantitative determination of the baryon and Dark Matter abundances, we have implemented and solved a system of coupled Boltzmann equations for the particle species involved in their generation, including all the relevant processes. In the most simple, but still general, limit, in which the processes determining the abundance and the decay rate of the Bino are mediated by degenerate right-handed squarks, the correct values of the DM and baryon relic densities are achieved for a Bino mass between 50 and 100 TeV, Gluino NLSP mass in the range 15-60 TeV and a gravitino mass between 100 GeV and few TeV. These high masses are unfortunately beyond the kinematical reach of LHC. On the contrary, an antiproton signal from the decays of the gravitino LSP might be within the sensitivity of AMS-02 and gamma-ray telescopes.


title: On Fej\'er monotone sequences and nonexpansive mappings

abstract: The notion of Fej\'er monotonicity has proven to be a fruitful concept in fixed point theory and optimization. In this paper, we present new conditions sufficient for convergence of Fej\'er monotone sequences and we also provide applications to the study of nonexpansive mappings. Various examples illustrate our results.


title: Entangling two transportable neutral atoms via local spin exchange

abstract: To advance quantum information science a constant pursuit is the search for physical systems that meet the stringent requirements for creating and preserving quantum entanglement. In atomic physics, robust two-qubit entanglement is typically achieved by strong, long-range interactions in the form of Coulomb interactions between ions or dipolar interactions between Rydberg atoms. While these interactions allow fast gates, atoms subject to these interactions must overcome the associated coupling to the environment and cross-talk among qubits. Local interactions, such as those requiring significant wavefunction overlap, can alleviate these detrimental effects yet present a new challenge: To distribute entanglement, qubits must be transported, merged for interaction, and then isolated for storage and subsequent operations. Here we show how, via a mobile optical tweezer, it is possible to prepare and locally entangle two ultracold neutral atoms, and then separate them while preserving their entanglement. While ground-state neutral atom experiments have measured dynamics consistent with spin entanglement, and detected entanglement with macroscopic observables, we are now able to demonstrate position-resolved two-particle coherence via application of a local gradient and parity measurements; this new entanglement-verification protocol could be applied to arbitrary spin-entangled states of spatially-separated atoms. The local entangling operation is achieved via ultracold spin-exchange interactions, and quantum tunneling is used to combine and separate atoms. Our toolset provides a framework for dynamically entangling remote qubits via local operations within a large-scale quantum register.


title: Distinguishing black-hole spin-orbit resonances by their gravitational wave signatures. II: Full parameter estimation

abstract: Gravitational waves from coalescing binary black holes encode the evolution of their spins prior to merger. In the post-Newtonian regime and on the precession timescale, this evolution has one of three morphologies, with the spins either librating around one of two fixed points ("resonances") or circulating freely. In this work we perform full parameter estimation on resonant binaries with fixed masses and spin magnitudes, changing three parameters: a conserved "projected effective spin" $\xi$ and resonant family $\Delta\Phi=0,\pi$ (which uniquely label the source), the inclination $\theta_{JN}$ of the binary's total angular momentum with respect to the line of sight (which determines the strength of precessional effects in the waveform), and the signal amplitude. We demonstrate that resonances can be distinguished for a wide range of binaries, except for highly symmetric configurations where precessional effects are suppressed. Motivated by new insight into double-spin evolution, we introduce new variables to characterize precessing black hole binaries which naturally reflects the timescale separation of the system and therefore better encode the dynamical information carried by gravitational waves.


title: Orientational properties of nematic disclinations

abstract: Topological defects play a pivotal role in the physics of liquid crystals and represent one of the most prominent and well studied aspects of mesophases. While in two-dimensional nematics, disclinations are traditionally treated as point-like objects, recent experimental studies on active nematics have suggested that half-strength disclinations might in fact possess a polar structure. In this article, we provide a precise definition of polarity for half-strength nematic disclinations, we introduce a simple and robust method to calculate this quantity from experimental and numerical data and we investigate how the orientational properties of half-strength disclinations affect their relaxational dynamics.


title: The sharpness of gamma-ray burst prompt emission spectra

abstract: We aim to obtain a measure of the curvature of time-resolved spectra that can be compared directly to theory. This tests the ability of models such as synchrotron emission to explain the peaks or breaks of GBM prompt emission spectra. We take the burst sample from the official Fermi GBM GRB time-resolved spectral catalog. We re-fit all spectra with a measured peak or break energy in the catalog best-fit models in various energy ranges, which cover the curvature around the spectral peak or break, resulting in a total of 1,113 spectra being analysed. We compute the sharpness angles under the peak or break of the triangle constructed under the model fit curves and compare to the values obtained from various representative emission models: blackbody, single-electron synchrotron, synchrotron emission from a Maxwellian or power-law electron distribution. We find that 35% of the time-resolved spectra are inconsistent with the single-electron synchrotron function, and 91% are inconsistent with the Maxwellian synchrotron function. The single temperature, single emission time and location blackbody function is found to be sharper than all the spectra. No general evolutionary trend of the sharpness angle is observed, neither per burst nor for the whole population. It is found that the limiting case, a single temperature Maxwellian synchrotron function, can only contribute up to $58^{+23}_{-18}$% of the peak flux. Our results show that even the sharpest but non-realistic case, the single-electron synchrotron function, cannot explain a large fraction of the observed GRB prompt spectra. Because of the fact that any combination of physically possible synchrotron spectra added together will always further broaden the spectrum, emission mechanisms other than optically thin synchrotron radiation are likely required in a full explanation of the spectral peaks or breaks of the GRB prompt emission phase.


title: First-principles nonequilibrium Green's function approach to transient photoabsorption: Application to atoms

abstract: We put forward a first-principle NonEquilibrium Green's Function (NEGF) approach to calculate the transient photoabsorption spectrum of optically thin samples. The method can deal with pump fields of arbitrary strength, frequency and duration as well as for overlapping and nonoverlapping pump and probe pulses. The electron-electron repulsion is accounted for by the correlation self-energy, and the resulting numerical scheme deals with matrices that scale quadratically with the system size. Two recent experiments, the first on helium and the second on krypton, are addressed. For the first experiment we explain the bending of the Autler-Townes absorption peaks with increasing the pump-probe delay $\t$, and relate the bending to the thickness and density of the gas. For the second experiment we find that sizable spectral structures of the pump-generated admixture of Kr ions are fingerprints of {\em dynamical correlation} effects, and hence they cannot be reproduced by time-local self-energy approximations. Remarkably, the NEGF approach also captures the retardation of the absorption onset of Kr$^{2+}$ with respect to Kr$^{1+}$ as a function of $\t$.


title: Estimating the Impact of Unknown Unknowns on Aggregate Query Results

abstract: It is common practice for data scientists to acquire and integrate disparate data sources to achieve higher quality results. But even with a perfectly cleaned and merged data set, two fundamental questions remain: (1) is the integrated data set complete and (2) what is the impact of any unknown (i.e., unobserved) data on query results? In this work, we develop and analyze techniques to estimate the impact of the unknown data (a.k.a., unknown unknowns) on simple aggregate queries. The key idea is that the overlap between different data sources enables us to estimate the number and values of the missing data items. Our main techniques are parameter-free and do not assume prior knowledge about the distribution. Through a series of experiments, we show that estimating the impact of unknown unknowns is invaluable to better assess the results of aggregate queries over integrated data sources.


title: The Power of Quantum Fourier Sampling

abstract: A line of work initiated by Terhal and DiVincenzo and Bremner, Jozsa, and Shepherd, shows that quantum computers can efficiently sample from probability distributions that cannot be exactly sampled efficiently on a classical computer, unless the PH collapses. Aaronson and Arkhipov take this further by considering a distribution that can be sampled efficiently by linear optical quantum computation, that under two feasible conjectures, cannot even be approximately sampled classically within bounded total variation distance, unless the PH collapses. In this work we use Quantum Fourier Sampling to construct a class of distributions that can be sampled by a quantum computer. We then argue that these distributions cannot be approximately sampled classically, unless the PH collapses, under variants of the Aaronson and Arkhipov conjectures. In particular, we show a general class of quantumly sampleable distributions each of which is based on an "Efficiently Specifiable" polynomial, for which a classical approximate sampler implies an average-case approximation. This class of polynomials contains the Permanent but also includes, for example, the Hamiltonian Cycle polynomial, and many other familiar #P-hard polynomials. Although our construction, unlike that proposed by Aaronson and Arkhipov, likely requires a universal quantum computer, we are able to use this additional power to weaken the conjectures needed to prove approximate sampling hardness results.


title: An Efficient Solver for Sparse Linear Systems Based on Rank-Structured Cholesky Factorization

abstract: Direct factorization methods for the solution of large, sparse linear systems that arise from PDE discretizations are robust, but typically show poor time and memory scalability for large systems. In this paper, we describe an efficient sparse, rank-structured Cholesky algorithm for solution of the positive definite linear system $A x = b$ when $A$ comes from a discretized partial-differential equation. Our approach combines the efficient memory access patterns of conventional supernodal Cholesky algorithms with the memory efficiency of rank-structured direct solvers. For several test problems arising from PDE discretizations, our method takes less memory than standard sparse Cholesky solvers and less wall-clock time than standard preconditioned iterations.