title: DarkBit: A GAMBIT module for computing dark matter observables and likelihoods

abstract: We introduce DarkBit, an advanced software code for computing dark matter constraints on various extensions to the Standard Model of particle physics, comprising both new native code and interfaces to external packages. This release includes a dedicated signal yield calculator for gamma-ray observations, which significantly extends current tools by implementing a cascade decay Monte Carlo, as well as a dedicated likelihood calculator for current and future experiments (gamlike). This provides a general solution for studying complex particle physics models that predict dark matter annihilation to a multitude of final states. We also supply a direct detection package that models a large range of direct detection experiments (DDcalc), and provides the corresponding likelihoods for arbitrary combinations of spin-independent and spin-dependent scattering processes. Finally, we provide custom relic density routines along with interfaces to DarkSUSY, micrOMEGAs, and the neutrino telescope likelihood package nuLike. DarkBit is written in the framework of the Global And Modular Beyond the Standard Model Inference Tool (GAMBIT), providing seamless integration into a comprehensive statistical fitting framework that allows users to explore new models with both particle and astrophysics constraints, and a consistent treatment of systematic uncertainties. In this paper we describe its main functionality, provide a guide to getting started quickly, and show illustrative examples for results obtained with DarkBit (both as a standalone tool and as a GAMBIT module). This includes a quantitative comparison between two of the main dark matter codes (DarkSUSY and micrOMEGAs), and application of DarkBit's advanced direct and indirect detection routines to a simple effective dark matter model.


title: CEPC Precision of Electroweak Oblique Parameters and Weakly Interacting Dark Matter: the Scalar Case

abstract: We investigate the sensitivity to weakly interacting scalar dark matter from future determination of electroweak oblique parameters in the Circular Electron-Positron Collider (CEPC) project. As illuminating examples, three dark matter models with scalar electroweak multiplets are studied. The multiplet couplings to the standard model Higgs doublet can break the mass degeneracy among the components, leading to nonzero contributions to oblique parameters. The dark matter candidate in these model is either a CP-even or CP-odd scalar, whose trilinear coupling to the Higgs boson could induce direct detection signals. For some moderate coupling values, we find that the CEPC sensitivity can be better than current direct detection experiments, exploring up to a TeV mass scale.


title: Characterizing the stellar population of NGC 1980

abstract: NGC 1980 is a young cluster that is located about 0.5 degrees south of the Orion Nebula Cluster (ONC). Recent studies by Bouy et al. and Pillitteri et al. have suggested that NGC 1980 contains an older population of stars compared to a much younger ONC, and that it belongs to a foreground population that may be located in front of the Orion A molecular gas by as much as 40 pc. In this work we present low-resolution spectra towards 148 young stars found towards the NGC 1980 region. We determine the spectral types of these stars, examine accretion signatures and measure the extinction towards them. We determine that based on these observations, the age of the population of NGC 1980 is indistinguishable from L1641, estimated to be ~3 Myr, comparable with the study by Fang et al.


title: Cavity-induced anti-correlated photon emission rates of a single ion

abstract: We report on the alteration of photon emission properties of a single trapped ion coupled to a high finesse optical fiber cavity. We show that the vacuum field of the cavity can simultaneously affect the emissions in both the infrared (IR) and ultraviolet (UV) branches of the $\Lambda-$type level system of $^{40}\mathrm{Ca}^+$ despite the cavity coupling only to the IR transition. The cavity induces strong emission in the IR transition through the Purcell effect resulting in a simultaneous suppression of the UV fluorescence. The measured suppression of this fluorescence is as large as 66% compared with the case without the cavity. Through analysis of the measurement results, we have obtained an ion-cavity coupling of $\bar{g}_0 = 2\pi\cdot (5.3 \pm 0.1)$ MHz, the largest ever reported so far for a single ion in the IR domain.


title: Numerical binary black hole mergers in dynamical Chern-Simons: I. Scalar field

abstract: Testing general relativity in the non-linear, dynamical, strong-field regime of gravity is one of the major goals of gravitational wave astrophysics. Performing precision tests of general relativity (GR) requires numerical inspiral, merger, and ringdown waveforms for binary black hole (BBH) systems in theories beyond GR. Currently, GR and scalar-tensor gravity are the only theories amenable to numerical simulations. In this article, we present a well-posed perturbation scheme for numerically integrating beyond-GR theories that have a continuous limit to GR. We demonstrate this scheme by simulating BBH mergers in dynamical Chern-Simons gravity (dCS), to linear order in the perturbation parameter. We present mode waveforms and energy fluxes of the dCS pseudoscalar field from our numerical simulations. We find good agreement with analytic predictions at early times, including the absence of pseudoscalar dipole radiation. We discover new phenomenology only accessible through numerics: a burst of dipole radiation during merger. We also quantify the self-consistency of the perturbation scheme. Finally, we estimate bounds that GR-consistent LIGO detections could place on the new dCS length scale, approximately $\ell \lesssim \mathcal{O}(10)~\mathrm{km}$.


title: Highlights from BNL and RHIC 2016

abstract: Highlights of news from Brookhaven National Laboratory (BNL) and results from the Relativistic Heavy Ion Collider (RHIC) in the period July 2015-2016 are presented. Transverse single spin asymmetries from polarized p+p collisions are presented for $\pi^0$ and jets as a function of Feynman $x_F$. An energy scan to study the $\sqrt{s}$ dependence of collectivity and flow for small systems was performed as well as a high luminosity Au$+$Au run. The failure of a quench protection diode resulted in a pause in the run to replace it, but otherwise performance of RHIC was the best ever. Experimental results discussed are an elegant measurement from STAR of the force between anti-protons using HBT correlations, flow in U+U collisions, an improved method of generating constituent quarks by PHENIX and new Number of Quark Participants (NQP) scaling of $E_T$ distributions in p$+$p,d$+$Au and Au$+$Au which worked well. New hard-scattering results as a function of $\sqrt{s}$ in Au$+$Au central collisions are presented. Also, measurements of the di-hadron acoplanarity for $\pi^0 +h$ and $\gamma+h$ in p+p collisions at $\sqrt{s}=500$ GeV are presented in terms of the out-of-plane transverse momentum $p_{\rm out}$ which differ from the prediction of the TMD framework of parton transverse momentum dynamics.


title: Upstream Causes of Downstream Effects

abstract: The United States Environmental Protection Agency considers nutrient pollution in stream ecosystems one of the U.S. most pressing environmental challenges. But limited independent replicates, lack of experimental randomization, and space- and time-varying confounding handicap causal inference on effects of nutrient pollution. In this paper the causal g-methods developed by Robins and colleagues are extended to allow for exposures to vary in time and space in order to assess the effects of nutrient pollution on chlorophyll a, a proxy for algal production. Publicly available data from the North Carolina Cape Fear River and a simulation study are used to show how causal effects of upstream nutrient concentrations on downstream chlorophyll a levels may be estimated from typical water quality monitoring data. Estimates obtained from the parametric g-formula, a marginal structural model, and a structural nested model indicate that chlorophyll a concentrations at Lock and Dam 1 were influenced by nitrate concentrations measured 86 to 109 km upstream, an area where four major industrial and municipal point sources discharge wastewater.


title: Modelling double charge exchange response function for tetraneutron system

abstract: This work is an attempt to model the $4n$ response function of a recent RIKEN experimental study of the double charge exchange $^4$He($^8$He,$^8$Be)$^4$n reaction in order to put in evidence an eventual enhancement mechanism of the zero energy cross section, including a near-threshold resonance. This resonance can indeed be reproduced only by adding to the standard nuclear Hamiltonian an unphysically large T=3/2 attractive 3n-force which destroys the neighboring nuclear chart. No other mechanisms like cusps or related structures were found.


title: Long Term Photometric Variability in Kepler Full Frame Images: Magnetic Cycles of Sun-Like Stars

abstract: Photometry from the Kepler mission is optimized to detect small, short duration signals like planet transits at the expense of long-term trends. This long-term variability can be recovered in photometry from the Full Frame Images (FFIs), a set of calibration data collected approximately monthly during the Kepler mission. Here, we present f3, an open-source package to perform photometry on the Kepler FFIs in order to detect changes in the brightness of stars in the Kepler field of view over long time baselines. We apply this package to a sample of 4,000 Sun-like stars with measured rotation periods. We find $\approx 10\%$ of these targets have long-term variability in their observed flux. For the majority of targets we identify the luminosity variations are either correlated or anticorrelated with the short-term variability due to starspots on the stellar surface. We find a transition between anticorrelated (starspot-dominated) variability and correlated (facula-dominated) variability between rotation periods of 15 and 25 days, suggesting the transition between the two modes is complete for stars at the age of the Sun. We also identify a sample of stars with apparently complete cycles, as well as a collection of short-period binaries with extreme photometric variation over the Kepler mission.


title: Visualizing the BEC-BCS crossover in the two-dimensional Fermi gas: pairing gaps and dynamical response functions from ab initio computations

abstract: Experiments with ultracold atoms provide a highly controllable laboratory setting with many unique opportunities for precision exploration of quantum many-body phenomena. The nature of such systems, with strong interaction and quantum entanglement, makes reliable theoretical calculations challenging. Especially difficult are excitation and dynamical properties, which are often the most directly relevant to experiment. We carry out exact numerical calculations, by Monte Carlo sampling of imaginary-time propagation of Slater determinants, to compute the pairing gap in the two-dimensional Fermi gas from first principles. Applying state-of-art analytic continuation techniques, we obtain the spectral function, and the density and spin structure factors providing unique tools to visualize the BEC-BCS crossover. These quantities will allow for a direct comparison with experiments.


title: Nonlinear saturation of the Weibel instability

abstract: The growth and saturation of magnetic fields due to the Weibel instability (WI) have important implications for laboratory and astrophysical plasmas, and this has drawn significant interest recently. Since the WI can generate a large magnetic field from no initial field, the maximum magnitudes achieved can have significant consequences for a number of applications. Hence, an understanding of the detailed dynamics driving the nonlinear saturation of the WI is important. This work considers the nonlinear saturation of the WI when counter-streaming populations of initially unmagnetized electrons are perturbed by a magnetic field oriented perpendicular to the direction of streaming. Previous works have found magnetic trapping to be important and connected electron skin depth spatial scales to the nonlinear saturation of the WI. 2 Results presented in this work are consistent with these findings for a high-temperature case. However, using a high-order continuum kinetic simulation tool, this work demonstrates that, when the electron populations are colder, a significant electrostatic potential develops that works with the magnetic field to create potential wells. The electrostatic field develops due to transverse flows induced by the WI, and in some cases is strengthened by a secondary instability. This field plays a key role in saturation of the WI for colder populations. The role of the electrostatic potential in Weibel instability saturation has not been studied in detail previously.


title: Status of the scalar singlet dark matter model

abstract: One of the simplest viable models for dark matter is an additional neutral scalar, stabilised by a $\mathbb{Z}_2$ symmetry. Using the GAMBIT package and combining results from four independent samplers, we present Bayesian and frequentist global fits of this model. We vary the singlet mass and coupling along with 13 nuisance parameters, including nuclear uncertainties relevant for direct detection, the local dark matter density, and selected quark masses and couplings. We include the dark matter relic density measured by Planck, direct searches with LUX, PandaX, SuperCDMS and XENON100, limits on invisible Higgs decays from the Large Hadron Collider, searches for high-energy neutrinos from dark matter annihilation in the Sun with IceCube, and searches for gamma rays from annihilation in dwarf galaxies with the Fermi-LAT. Viable solutions remain at couplings of order unity, for singlet masses between the Higgs mass and about 300 GeV, and at masses above $\sim$1 TeV. Only in the latter case can the scalar singlet constitute all of dark matter. Frequentist analysis shows that the low-mass resonance region, where the singlet is about half the mass of the Higgs, can also account for all of dark matter, and remains viable. However, Bayesian considerations show this region to be rather fine-tuned.


title: Metric Mahler measures over number fields

abstract: For an algebraic number $\alpha$, the metric Mahler measure $m_1(\alpha)$ was first studied by Dubickas and Smyth in 2001 and was later generalized to the $t$-metric Mahler measure $m_t(\alpha)$ by the author in 2010. The definition of $m_t(\alpha)$ involves taking an infimum over a certain collection $N$-tuples of points in $\overline{\mathbb Q}$, and from previous work of Jankauskas and the author, the infimum in the definition of $m_t(\alpha)$ is attained by rational points when $\alpha\in \mathbb Q$. As a consequence of our main theorem in this article, we obtain an analog of this result when $\mathbb Q$ is replaced with any imaginary quadratic number field of class number equal to $1$. Further, we study examples of other number fields to which our methods may be applied, and we establish various partial results in those cases.


title: FlavBit: A GAMBIT module for computing flavour observables and likelihoods

abstract: Flavour physics observables are excellent probes of new physics up to very high energy scales. Here we present FlavBit, the dedicated flavour physics module of the global-fitting package GAMBIT. FlavBit includes custom implementations of various likelihood routines for a wide range of flavour observables, including detailed uncertainties and correlations associated with LHCb measurements of rare, leptonic and semileptonic decays of B and D mesons, kaons and pions. It provides a generalised interface to external theory codes such as SuperIso, allowing users to calculate flavour observables in and beyond the Standard Model, and then test them in detail against all relevant experimental data. We describe FlavBit and its constituent physics in some detail, then give examples from supersymmetry and effective field theory illustrating how it can be used both as a standalone library for flavour physics, and within GAMBIT.


title: Radiation reaction for spinning bodies in effective field theory I: Spin-orbit effects

abstract: We compute the leading Post-Newtonian (PN) contributions at linear order in the spin to the radiation-reaction acceleration and spin evolution for binary systems, which enter at fourth PN order. The calculation is carried out, from first principles, using the effective field theory framework for spinning compact objects, in both the Newton-Wigner and covariant spin supplementary conditions. A non-trivial consistency check is performed on our results by showing that the energy loss induced by the resulting radiation-reaction force is equivalent to the total emitted power in the far zone, up to so-called "Schott terms." We also find that, at this order, the radiation reaction has no net effect on the evolution of the spins. The spin-spin contributions to radiation reaction are reported in a companion paper.


title: Global fits of GUT-scale SUSY models with GAMBIT

abstract: We present the most comprehensive global fits to date of three supersymmetric models motivated by grand unification: the Constrained Minimal Supersymmetric Standard Model (CMSSM), and its Non-Universal Higgs Mass generalisations NUHM1 and NUHM2. We include likelihoods from a number of direct and indirect dark matter searches, a large collection of electroweak precision and flavour observables, direct searches for supersymmetry at LEP and Runs I and II of the LHC, and constraints from Higgs observables. Our analysis improves on existing results not only in terms of the number of included observables, but also in the level of detail with which we treat them, our sampling techniques for scanning the parameter space, and our treatment of nuisance parameters. We show that stau co-annihilation is now ruled out in the CMSSM at more than 95\% confidence. Stop co-annihilation turns out to be one of the most promising mechanisms for achieving an appropriate relic density of dark matter in all three models, whilst avoiding all other constraints. We find high-likelihood regions of parameter space featuring light stops and charginos, making them potentially detectable in the near future at the LHC. We also show that tonne-scale direct detection will play a largely complementary role, probing large parts of the remaining viable parameter space, including essentially all models with multi-TeV neutralinos.


title: SpecBit, DecayBit and PrecisionBit: GAMBIT modules for computing mass spectra, particle decay rates and precision observables

abstract: We present the GAMBIT modules SpecBit, DecayBit and PrecisionBit. Together they provide a new framework for linking publicly available spectrum generators, decay codes and other precision observable calculations in a physically and statistically consistent manner. This allows users to automatically run various combinations of existing codes as if they are a single package. The modular design allows software packages fulfilling the same role to be exchanged freely at runtime, with the results presented in a common format that can be easily passed to downstream dark matter, collider and flavour codes. These modules constitute an essential part of the broader GAMBIT framework, a major new software package for performing global fits. In this paper we present the observable calculations, data, and likelihood functions implemented in the three modules, as well as the conventions and assumptions used in interfacing them with external codes. We also present 3-BIT-HIT, a command-line utility for computing mass spectra, couplings, decays and precision observables in the MSSM, which shows how the three modules can be easily used independently of GAMBIT.


title: On the cohomology of the mapping class group of the punctured projective plane

abstract: The mapping class group $\Gamma^k(N_g)$ of a non-orientable surface with punctures is studied via classical homotopy theory of configuration spaces. In particular, we obtain a non-orientable version of the Birman exact sequence. In the case of $\mathbb R {\rm P}^2$, we analize the Serre spectral sequence of a fiber bundle $F_k(\mathbb R {\rm P}^2)/\Sigma_k \to X_k \to BSO(3)$ where $X_k$ is a $K(\Gamma^k(\mathbb R {\rm P}^2),1)$ and $F_k(\mathbb R {\rm P}^2)/\Sigma_k$ denotes the configuration space of unordered $k$-tuples of distinct points in $\mathbb R {\rm P}^2$. As a consequence, we express the mod-2 cohomology of $\Gamma^k(\mathbb R {\rm P}^2)$ in terms of that of $F_k(\mathbb R {\rm P}^2)/\Sigma_k$.


title: Radiation reaction for spinning bodies in effective field theory II: Spin-spin effects

abstract: We compute the leading Post-Newtonian (PN) contributions at quadratic order in the spins to the radiation-reaction acceleration and spin evolution for binary systems, entering at four-and-a-half PN order. Our calculation includes the back-reaction from finite-size spin effects, which is presented for the first time. The computation is carried out, from first principles, using the effective field theory framework for spinning extended objects. At this order, nonconservative effects in the spin-spin sector are independent of the spin supplementary conditions. A non-trivial consistency check is performed by showing that the energy loss induced by the resulting radiation-reaction force is equivalent to the total emitted power in the far zone. We find that, in contrast to the spin-orbit contributions (reported in a companion paper), the radiation reaction affects the evolution of the spin vectors once spin-spin effects are incorporated.


title: Generalized Extension of Watson's theorem for the series $_{3}F_{2}(1)$

abstract: The $_{3}F_{2}$ hypergeometric function plays a very significant role in the theory of hypergeometric and generalized hypergeometric series. Despite that $_{3}F_{2}$ hypergeometric function has several applications in mathematics, also it has a lot of applications in physics and statistics. The fundamental purpose of this research paper is to find out the explicit expression of the $_{3}F_{2}$ Watson's classical summation theorem of the form: \[ _{3}F_{2}\left[ \begin{array} [c]{ccccc}% a, & b, & c & & \\ & & & ; & 1\\ \frac{1}{2}(a+b+i+1), & 2c+j & & & \end{array} \right] \] with arbitrary $i$ and $j$, where for $i=j=0$, we get the well known Watson's theorem for the series $_{3}F_{2}(1)$.


title: Joint Uplink/Downlink Resource Allocation and Data Offloading in OFDMA-Based Wireless Powered HetNets

abstract: This paper considers joint uplink/downlink of an orthogonal frequency division multiple access (OFDMA)-based heterogeneous network (HetNet) consisting of a single macro base station (MBS), multiple femto base stations (FBSs) and access points (APs) where base stations (BSs) can offload data to APs and each mobile user (MU) is able to harvest the received energy using the simultaneous wireless information and power transfer (SWIPT) technique. We also suppose that the harvested energy of MUs are used for their uplink information transmission. We devise a radio resource allocation (RRA) algorithm to maximize the uplink sum data rate of MUs subject to a minimum required downlink data rate of each MU and maximum allowable transmit power of each BS, AP, and MU. More specifically, both the frequency division duplex (FDD) and time division duplex (TDD) schemes are investigated. The proposed non-convex optimization problems are solved using an iterative algorithm. It is also proved that the proposed algorithm converges to a near-optimal solution. Simulation results illustrate that the TDD scheme improves the performance compared to the FDD scheme. In addition, it is shown that utilizing the data offloading technique improves the uplink sum data rate of MUs compared to the scenario without any AP.


title: Prediction Measures in Nonlinear Beta Regression Models

abstract: Nonlinear models are frequently applied to determine the optimal supply natural gas to a given residential unit based on economical and technical factors, or used to fit biochemical and pharmaceutical assay nonlinear data. In this article we propose PRESS statistics and prediction coefficients for a class of nonlinear beta regression models, namely $P^2$ statistics. We aim at using both prediction coefficients and goodness-of-fit measures as a scheme of model select criteria. In this sense, we introduce for beta regression models under nonlinearity the use of the model selection criteria based on robust pseudo-$R^2$ statistics. Monte Carlo simulation results on the finite sample behavior of both prediction-based model selection criteria $P^2$ and the pseudo-$R^2$ statistics are provided. Three applications for real data are presented. The linear application relates to the distribution of natural gas for home usage in S\~ao Paulo, Brazil. Faced with the economic risk of too overestimate or to underestimate the distribution of gas has been necessary to construct prediction limits and to select the best predicted and fitted model to construct best prediction limits it is the aim of the first application. Additionally, the two nonlinear applications presented also highlight the importance of considering both goodness-of-predictive and goodness-of-fit of the competitive models.


title: Energy Stored by Radiating Systems

abstract: Though commonly used to calculate Q-factor and fractional bandwidth, the energy stored by radiating systems (antennas) is a subtle and challenging concept that has perplexed researchers for over half a century. Here, the obstacles in defining and calculating stored energy in general electromagnetic systems are presented from first principles as well as using demonstrative examples from electrostatics, circuits, and radiating systems. Along the way, the concept of unobservable energy is introduced to formalize such challenges. Existing methods of defining stored energy in radiating systems are then reviewed in a framework based on technical commonalities rather than chronological order. Equivalences between some methods under common assumptions are highlighted, along with the strengths, weaknesses, and unique applications of certain techniques. Numerical examples are provided to compare the relative margin between methods on several radiating structures.


title: On Distinguishability of Black Hole Microstates

abstract: We use the Holevo information to estimate distinguishability of microstates of a black hole in anti-de Sitter space by measurements one can perform on a subregion of a Cauchy surface of the dual conformal field theory. We find that microstates are not distinguishable at all until the subregion reaches a certain size and that perfect distinguishability can be achieved before the subregion covers the entire Cauchy surface. We will compare our results with expectations from the entanglement wedge reconstruction, tensor network models, and the bit threads interpretation of the Ryu-Takayanagi formula.


title: Constraining the clustering transition for colorings of sparse random graphs

abstract: Let $\Omega_q$ denote the set of proper $q$-colorings of the random graph $G_{n,m}, m=dn/2$ and let $H_q$ be the graph with vertex set $\Omega_q$ and an edge $\{\sigma,\tau\}$ where $\sigma,\tau$ are mappings $[n]\to[q]$ iff $h(\sigma,\tau)=1$. Here $h(\sigma,\tau)$ is the Hamming distance $|\{v\in [n]:\sigma(v)\neq\tau(v)\}|$. We show that w.h.p. $H_q$ contains a single giant component containing almost all colorings in $\Omega_q$ if $d$ is sufficiently large and $q\geq \frac{cd}{\log d}$ for a constant $c>3/2$.


title: Inflation with $f(R,\phi)$ in Jordan frame

abstract: We consider an $f(R)$ action that is non-minimally coupled to a massive scalar field. The model closely resembles scalar-tensor theory and by conformal transformation can be transformed to Einstein frame. To avoid the ambiguity of the frame dependence, we obtain an exact analytical solution in Jordan frame and show that the model leads to a period of accelerated expansion with an exit. Further, we compute the scalar and tensor power spectrum for the model and compare them with observations.


title: Spectral analysis of molecular resonances in erbium isotopes: Are they close to semi-Poisson?

abstract: We perform a thorough analysis of the spectral statistics of experimental molecular resonances, of bosonic erbium $^{166}$Er and $^{168}$Er isotopes, produced as a function of magnetic field($B$) by Frisch et al. [Nature 507, (2014) 475], utilizing some recently derived surmises which interpolate between Poisson and GOE and without unfolding. Supplementing this with an analysis using unfolded spectrum, it is shown that the resonances are close to semi-Poisson distribution. There is an earlier claim of missing resonances by Molina et al. [Phys. Rev. E 92, (2015) 042906]. These two interpretations can be tested by more precise measurements in future experiments.


title: Semiclassical theory of spin-orbit torques in disordered multiband electron systems

abstract: We study spin-orbit torques (SOT) in non-degenerate multiband electron systems in the weak disorder limit. In order to have better physical transparency a semiclassical Boltzmann approach equivalent to the Kubo diagrammatic approach in the non-crossing approximation is formulated. This semiclassical framework accounts for the interband- coherence effects induced by both the electric field and static impurity scattering. Using the two-dimensional Rashba ferromagnet as a model system, we show that the antidamping-like SOT arising from disorder-induced interband-coherence effects is very sensitive to the spin structure of disorder and may have the same sign as the intrinsic SOT in the presence of spin-dependent disorder. While the cancellation of this SOT and the intrinsic one occurs only in the case of spin-independent short-range disorder.


title: Viscosity solutions and the minimal surface system

abstract: We give a definition of viscosity solution for the minimal surface system and prove a version of Allard regularity theorem in this setting.


title: Trajectory-free approximation of phase space structures using the trajectory divergence rate

abstract: This paper introduces the trajectory divergence rate, a scalar field which locally gives the instantaneous attraction or repulsion of adjacent trajectories. This scalar field may be used to find highly attracting or repelling invariant manifolds, such as slow manifolds, to rapidly approximating hyperbolic Lagrangian coherent structures, or to provide the local stability of invariant manifolds. This work presents the derivation of the trajectory divergence rate and the related trajectory divergence ratio for 2-dimensional systems, investigates their properties, shows their application to several example systems, and presents their extension to higher dimensions.


title: Variable Screening for High Dimensional Time Series

abstract: Variable selection is a widely studied problem in high dimensional statistics, primarily since estimating the precise relationship between the covariates and the response is of great importance in many scientific disciplines. However, most of theory and methods developed towards this goal for the linear model invoke the assumption of iid sub-Gaussian covariates and errors. This paper analyzes the theoretical properties of Sure Independence Screening (SIS) (Fan and Lv [J. R. Stat. Soc. Ser. B Stat. Methodol. 70 (2008) 849-911]) for high dimensional linear models with dependent and/or heavy tailed covariates and errors. We also introduce a generalized least squares screening (GLSS) procedure which utilizes the serial correlation present in the data. By utilizing this serial correlation when estimating our marginal effects, GLSS is shown to outperform SIS in many cases. For both procedures we prove sure screening properties, which depend on the moment conditions, and the strength of dependence in the error and covariate processes, amongst other factors. Additionally, combining these screening procedures with the adaptive Lasso is analyzed. Dependence is quantified by functional dependence measures (Wu [Proc. Natl. Acad. Sci. USA 102 (2005) 14150-14154]), and the results rely on the use of Nagaev-type and exponential inequalities for dependent random variables. We also conduct simulations to demonstrate the finite sample performance of these procedures, and include a real data application of forecasting the US inflation rate.


title: Tourists' digital footprint in cities: comparing big data sources

abstract: There is little knowledge available on the spatial behaviour of urban tourists, and yet tourists generate an enormous quantity of data (Big Data) when they visit cities. These data sources can be used to track their presence through their activities. The aim of this paper is to analyse the digital footprint of urban tourists through Big Data. Unlike other papers that use a single data source, this article examines three sources of data to reflect different tourism activities in cities: Panoramio (sightseeing), Foursquare (consumption), and Twitter (being connected). Tourist density in the three data sources is compared via maps, correlation analysis (OLS) and spatial self-correlation analysis (Global Moran's I statistic and LISA). Finally the data are integrated using cluster analysis and combining the spatial clusters identified in the LISA analysis in the different data sources. The results show that the data from the three activities are partly spatially redundant and partly complementary, and allow the characterisation of multifunction tourist spaces (with several activities) and spaces specialising in one or various activities (for example, sightseeing and consumption). The case study analysed (Madrid) reveals a significant presence of tourists in the city centre, and increasing specialisation from the centre outwards towards the periphery. The main conclusion of the paper is that it is not sufficient to use one data source to analyse the presence of tourists in cities; several must be used in a complementary manner.


title: Layer-by-layer assembly of patchy particles as a route to non-trivial structures

abstract: We propose a new strategy for robust high-quality self-assembly of non-trivial periodic structures out of patchy particles, and investigate it with Brownian Dynamics (BD) simulations. Its first element is the use of specific patch-patch and shell-shell interactions between the particles, that can be implemented through differential functionalization of patched and shell regions with specific DNA strands. The other key element of our approach is the use of layer-by-layer protocol that allows to avoid a formations of undesired random aggregates. As an example, we design and self-assemble in silico a version of a Double Diamond (DD) lattice in which four particle types are arranged into BCC crystal made of four FCC sub-lattices. The lattice can be further converted to Cubic Diamond (CD) by selective removal of the particles of certain types. Our results demonstrate that by combining the directionality, selectivity of interactions and the layer-by-layer protocol, a high-quality robust self-assembly can be achieved.


title: Brazilian Network of PhDs Working with Probability and Statistics

abstract: Statistical and probabilistic reasoning enlightens our judgments about uncertainty and the chance or beliefs on the occurrence of random events in everyday life. Therefore, there are scientists working with Probability and Statistics in various fields of knowledge, what favors the formation of scientific network collaborations of researchers with different backgrounds. Here, we propose to describe the Brazilian PhDs who work with probability and statistics. In particular, we analyze national and states collaboration networks of such researchers by calculating different metrics. We show that there is a greater concentration of nodes in and around the cites which host Probability and Statistics graduate programs. Moreover, the states that host P & S Doctoral programs are the most central. We also observe a disparity in the size of the states networks. The clustering coefficient of the national network suggests that this network and regional differences especially with respect to states from South-east and North is not cohesive and, probably, it is in a maturing stage


title: Nanophotonic modal dichroism: mode-multiplexed modulators

abstract: As the diffraction limit is approached, device miniaturization to integrate more functionality per area becomes more and more challenging. Here we propose a novel strategy to increase the functionality-per-area by exploiting the modal properties of a waveguide system. With such approach the design of a mode-multiplexed nanophotonic modulator relying on the mode-selective absorption of a patterned Indium-Tin-Oxide is proposed. Full-wave simulations of a device operating at the telecom wavelength of 1550nm show that two modes can be independently modulated, while maintaining performances in line with conventional single-mode ITO modulators reported in the recent literature. The proposed design principles can pave the way to a novel class of mode-multiplexed compact photonic devices able to effectively multiply the functionality-per-area in integrated photonic systems.


title: Reconstruction of the polarization distribution of the Rice-Mele model

abstract: We calculate the gauge invariant cumulants (and moments) associated with the Zak phase in the Rice-Mele model. We reconstruct the underlying probability distribution by maximizing the information entropy and applying the moments as constraints. When the Wannier functions are localized within one unit cell, the probability distribution so obtained corresponds to that of the Wannier function. We show that in the fully dimerized limit the magnitude of the moments are all equal. In this limit, if the on-site interaction is decreased towards zero, the distribution shifts towards the midpoint of the unit cell, but the overall shape of the distribution remains the same. Away from this limit, if alternate hoppings are finite, and the on-site interaction is decreased, the distribution also shifts towards the midpoint of the unit cell, but it does this by changing shape, by becoming asymmetric around the maximum, as well as by shifting. We also follow the probability distribution of the polarization in cycles around the topologically non-trivial point of the model. The distribution moves across to the next unit cell, its shape distorting considerably in the process. If the radius of the cycle is large, the shift of the distribution is accompanied by large variations in the maximum.


title: The pulse of the city through Twitter: relationships between land use and spatiotemporal demographics

abstract: Social network data offer interesting opportunities in urban studies. In this study, we used Twitter data to analyse city dynamics over the course of the day. Users of this social network were grouped according to city zone and time slot in order to analyse the daily dynamics of the city and the relationship between this and land use. First, daytime activity in each zone was compared with activity at night in order to determine which zones showed increased activity in each of the time slots. Then, typical Twitter activity profiles were obtained based on the predominant land use in each zone, indicating how land uses linked to activities were activated during the day, but at different rates depending on the type of land use. Lastly, a multiple regression analysis was performed to determine the influence of the different land uses on each of the major time slots (morning, afternoon, evening and night) through their changing coefficients. Activity tended to decrease throughout the day for most land uses (e.g. offices, education, health and transport), but remained constant in parks and increased in retail and residential zones. Our results show that social network data can be used to improve our understanding of the link between land use and urban dynamics.


title: Large Scale Empirical Risk Minimization via Truncated Adaptive Newton Method

abstract: We consider large scale empirical risk minimization (ERM) problems, where both the problem dimension and variable size is large. In these cases, most second order methods are infeasible due to the high cost in both computing the Hessian over all samples and computing its inverse in high dimensions. In this paper, we propose a novel adaptive sample size second-order method, which reduces the cost of computing the Hessian by solving a sequence of ERM problems corresponding to a subset of samples and lowers the cost of computing the Hessian inverse using a truncated eigenvalue decomposition. We show that while we geometrically increase the size of the training set at each stage, a single iteration of the truncated Newton method is sufficient to solve the new ERM within its statistical accuracy. Moreover, for a large number of samples we are allowed to double the size of the training set at each stage, and the proposed method subsequently reaches the statistical accuracy of the full training set approximately after two effective passes. In addition to this theoretical result, we show empirically on a number of well known data sets that the proposed truncated adaptive sample size algorithm outperforms stochastic alternatives for solving ERM problems.


title: X-rays from Green Pea Analogs

abstract: X-ray observations of two metal-deficient luminous compact galaxies (LCG) (SHOC~486 and SDSS J084220.94+115000.2) with properties similar to the so-called Green Pea galaxies were obtained using the {\emph{Chandra X-ray Observatory}}. Green Pea galaxies are relatively small, compact (a few kpc across) galaxies that get their green color from strong [OIII]$\lambda$5007\AA\ emission, an indicator of intense, recent star formation. These two galaxies were predicted to have the highest observed count rates, using the X-ray luminosity -- star formation rate ($L_X$--SFR) relation for X-ray binaries, from a statistically complete sample drawn from optical criteria. We determine the X-ray luminosity relative to star-formation rate and metallicity for these two galaxies. Neither exhibit any evidence of active galactic nuclei and we suspect the X-ray emission originates from unresolved populations of high mass X-ray binaries. We discuss the $L_X$--SFR--metallicity plane for star-forming galaxies and show that the two LCGs are consistent with the prediction of this relation. This is the first detection of Green Pea analogs in X-rays.


title: Comparison of statistical sampling methods with ScannerBit, the GAMBIT scanning module

abstract: We introduce ScannerBit, the statistics and sampling module of the public, open-source global fitting framework GAMBIT. ScannerBit provides a standardised interface to different sampling algorithms, enabling the use and comparison of multiple computational methods for inferring profile likelihoods, Bayesian posteriors, and other statistical quantities. The current version offers random, grid, raster, nested sampling, differential evolution, Markov Chain Monte Carlo (MCMC) and ensemble Monte Carlo samplers. We also announce the release of a new standalone differential evolution sampler, Diver, and describe its design, usage and interface to ScannerBit. We subject Diver and three other samplers (the nested sampler MultiNest, the MCMC GreAT, and the native ScannerBit implementation of the ensemble Monte Carlo algorithm T-Walk) to a battery of statistical tests. For this we use a realistic physical likelihood function, based on the scalar singlet model of dark matter. We examine the performance of each sampler as a function of its adjustable settings, and the dimensionality of the sampling problem. We evaluate performance on four metrics: optimality of the best fit found, completeness in exploring the best-fit region, number of likelihood evaluations, and total runtime. For Bayesian posterior estimation at high resolution, T-Walk provides the most accurate and timely mapping of the full parameter space. For profile likelihood analysis in less than about ten dimensions, we find that Diver and MultiNest score similarly in terms of best fit and speed, outperforming GreAT and T-Walk; in ten or more dimensions, Diver substantially outperforms the other three samplers on all metrics.

categories: hep-ph astro-ph.CO astro-ph.IM physics.data-an stat.CO


title: Hodge theory in combinatorics

abstract: George Birkhoff proved in 1912 that the number of proper colorings of a finite graph G with n colors is a polynomial in n, called the chromatic polynomial of G. Read conjectured in 1968 that for any graph G, the sequence of absolute values of coefficients of the chromatic polynomial is unimodal: it goes up, hits a peak, and then goes down. Read's conjecture was proved by June Huh in a 2012 paper making heavy use of methods from algebraic geometry. Huh's result was subsequently refined and generalized by Huh and Katz, again using substantial doses of algebraic geometry. Both papers in fact establish log-concavity of the coefficients, which is stronger than unimodality. The breakthroughs of Huh and Huh-Katz left open the more general Rota-Welsh conjecture where graphs are generalized to (not necessarily representable) matroids and the chromatic polynomial of a graph is replaced by the characteristic polynomial of a matroid. The Huh and Huh-Katz techniques are not applicable in this level of generality, since there is no underlying algebraic geometry to which to relate the problem. But in 2015 Adiprasito, Huh, and Katz announced a proof of the Rota-Welsh conjecture based on a novel approach motivated by but not making use of any results from algebraic geometry. The authors first prove that the Rota-Welsh conjecture would follow from combinatorial analogues of the Hard Lefschetz Theorem and Hodge-Riemann relations in algebraic geometry. They then implement an elaborate inductive procedure to prove the combinatorial Hard Lefschetz Theorem and Hodge-Riemann relations using purely combinatorial arguments. We will survey these developments.


title: Compatible extensions and consistent closures: a fuzzy approach

abstract: In this paper $\ast$--compatible extensions of fuzzy relations are studied, generalizing some results obtained by Duggan in case of crisp relations. From this general result are obtained as particular cases fuzzy versions of some important extension theorems for crisp relations (Szpilrajn, Hansson, Suzumura). Two notions of consistent closure of a fuzzy relation are introduced.


title: pix2code: Generating Code from a Graphical User Interface Screenshot

abstract: Transforming a graphical user interface screenshot created by a designer into computer code is a typical task conducted by a developer in order to build customized software, websites, and mobile applications. In this paper, we show that deep learning methods can be leveraged to train a model end-to-end to automatically generate code from a single input image with over 77% of accuracy for three different platforms (i.e. iOS, Android and web-based technologies).


title: Dipole-field-assisted charge extraction in metal-perovskite-metal back-contact solar cells

abstract: Hybrid organic-inorganic halide perovskites are low-cost solution-processable solar cell materials with photovoltaic properties that rival those of crystalline silicon. The perovskite films are typically sandwiched between thin layers of hole and electron transport materials, which efficiently extract photogenerated charges. This affords high-energy conversion efficiencies but results in significant performance and fabrication challenges. Herein we present a simple charge transport layer-free perovskite solar cell (PSC), comprising only a perovskite layer with two interdigitated gold back-contacts. Charge extraction is achieved via self-assembled molecular monolayers (SAMs) and their associated dipole fields at the metal/perovskite interface. Photovoltages of approximately 600 mV generated by SAM-modified PSCs are equivalent to the built-in potential generated by individual dipole layers. Efficient charge extraction results in photocurrents of up to 12.1 mA/cm2 under simulated sunlight, despite a large electrode spacing.


title: On-Site Wireless Power Generation

abstract: Conventional wireless power transfer systems consist of a microwave power generator and a microwave power receiver separated by some distance. To realize efficient power transfer, the system is typically brought to resonance, and the coupled-antenna mode is optimized to reduce radiation into the surrounding space. In this scheme, any modification of the receiver position or of its electromagnetic properties results in the necessity of dynamically tuning the whole system to restore the resonant matching condition. It implies poor robustness to the receiver location and load impedance, as well as additional energy consumption in the control network. In this study, we introduce a new paradigm for wireless power delivery based on which the whole system, including transmitter and receiver and the space in between, forms a unified microwave power generator. In our proposed scenario the load itself becomes part of the generator. Microwave oscillations are created directly at the receiver location, eliminating the need for dynamical tuning of the system within the range of the self-oscillation regime. The proposed concept has relevant connections with the recent interest in parity-time symmetric systems, in which balanced loss and gain distributions enable unusual electromagnetic responses.


title: Horocyclic invariance of Ruelle resonant states for contact Anosov flows in dimension 3

abstract: We show that for smooth contact Anosov flows in dimension 3, the resonant states associated to the first band of Ruelle resonances are distributions that are killed by the unstable derivative.


title: Using Sloppy Models for Constrained Emittance Minimization at the Cornell Electron Storage Ring (CESR)

abstract: In order to minimize the emittance at the Cornell Electron Storage Ring (CESR), we measure and correct the orbit, dispersion, and transverse coupling of the beam. However, this method is limited by finite measurement resolution of the dispersion, and so a new procedure must be used to further reduce the emittance due to dispersion. In order to achieve this, we use a method based upon the theory of sloppy models. We use a model of the accelerator to create the Hessian matrix which encodes the effects of various corrector magnets on the vertical emittance. A singular value decomposition of this matrix yields the magnet combinations which have the greatest effect on the emittance. We can then adjust these magnet "knobs" sequentially in order to decrease the dispersion and the emittance. We present here comparisons of the effectiveness of this procedure in both experiment and simulation using a variety of CESR lattices. We also discuss techniques to minimize changes to parameters we have already corrected.


title: Consistencies and inconsistencies between model selection and link prediction in networks

abstract: A principled approach to understand network structures is to formulate generative models. Given a collection of models, however, an outstanding key task is to determine which one provides a more accurate description of the network at hand, discounting statistical fluctuations. This problem can be approached using two principled criteria that at first may seem equivalent: selecting the most plausible model in terms of its posterior probability; or selecting the model with the highest predictive performance in terms of identifying missing links. Here we show that while these two approaches yield consistent results in most of cases, there are also notable instances where they do not, that is, where the most plausible model is not the most predictive. We show that in the latter case the improvement of predictive performance can in fact lead to overfitting both in artificial and empirical settings. Furthermore, we show that, in general, the predictive performance is higher when we average over collections of models that are individually less plausible, than when we consider only the single most plausible model.


title: High resolution quantum sensing with shaped control pulses

abstract: We investigate the application of amplitude-shaped control pulses for enhancing the time and frequency resolution of multipulse quantum sensing sequences. Using the electronic spin of a single nitrogen vacancy center in diamond and up to 10,000 coherent microwave pulses with a cosine square envelope, we demonstrate 0.6 ps timing resolution for the interpulse delay. This represents a refinement by over 3 orders of magnitude compared to the 2 ns hardware sampling. We apply the method for the detection of external AC magnetic fields and nuclear magnetic resonance signals of carbon-13 spins with high spectral resolution. Our method is simple to implement and especially useful for quantum applications that require fast phase gates, many control pulses, and high fidelity.


title: Radially symmetric connections over round spheres

abstract: We classify the radially symmetric connections in vector bundles over round spheres by proving that they are all parallel.


title: Atomic-scale Structural and Chemical Characterization of Hexagonal Boron Nitride Layers Synthesized at the Wafer-Scale with Monolayer Thickness Control

abstract: Hexagonal boron nitride (h-BN) is a promising two-dimensional insulator with a large band gap and low density of charged impurities that is isostructural and isoelectronic with graphene. Here we report the chemical and atomic-scale structure of CVD-grown wafer-scale (~25 cm2) h-BN sheets ranging in thickness from 1-20 monolayers. Atomic-scale images of h-BN on Au and graphene/Au substrates obtained by scanning tunneling microscopy (STM) reveal high h-BN crystalline quality in monolayer samples. Further characterization of 1-20 monolayer samples indicates uniform thickness for wafer-scale areas; this thickness control is a result of precise control of the precursor flow rate, deposition temperature and pressure. Raman and infrared spectroscopy indicate the presence of B-N bonds and reveal a linear dependence of thickness with growth time. X-ray photoelectron spectroscopy (XPS) shows the film stoichiometry, and the B/N atom ratio in our films is 1 + 0.6% across the range of thicknesses. Electrical current transport in metal/insulator/metal (Au/h-BN/Au) heterostructures indicates that our CVD-grown h-BN films can act as excellent tunnel barriers with a high hard-breakdown field strength. Our results suggest that large-area h-BN films are structurally, chemically and electronically uniform over the wafer scale, opening the door to pervasive application as a dielectric in layered nanoelectronic and nanophotonic heterostructures.


title: The shear and bulk relaxation times from the general correlation functions

abstract: In this paper we present two quantum field theoretical analyses on the shear and bulk relaxation times. First, we discuss how to find Kubo formulas for the shear and the bulk relaxation times. Next, we provide results on the shear viscosity relaxation time obtained within the diagrammatic approach for the massless $\lambda\phi^4$ theory.


title: Universal 3D Wearable Fingerprint Targets: Advancing Fingerprint Reader Evaluations

abstract: We present the design and manufacturing of high fidelity universal 3D fingerprint targets, which can be imaged on a variety of fingerprint sensing technologies, namely capacitive, contact-optical, and contactless-optical. Universal 3D fingerprint targets enable, for the first time, not only a repeatable and controlled evaluation of fingerprint readers, but also the ability to conduct fingerprint reader interoperability studies. Fingerprint reader interoperability refers to how robust fingerprint recognition systems are to variations in the images acquired by different types of fingerprint readers. To build universal 3D fingerprint targets, we adopt a molding and casting framework consisting of (i) digital mapping of fingerprint images to a negative mold, (ii) CAD modeling a scaffolding system to hold the negative mold, (iii) fabricating the mold and scaffolding system with a high resolution 3D printer, (iv) producing or mixing a material with similar electrical, optical, and mechanical properties to that of the human finger, and (v) fabricating a 3D fingerprint target using controlled casting. Our experiments conducted with PIV and Appendix F certified optical (contact and contactless) and capacitive fingerprint readers demonstrate the usefulness of universal 3D fingerprint targets for controlled and repeatable fingerprint reader evaluations and also fingerprint reader interoperability studies.


title: Probing the No-Scale ${\cal F}$-$SU(5)$ One-Parameter Model via Gluino Searches at the LHC2

abstract: In our recent paper entitled "The return of the King: No-Scale ${\cal F}$-$SU(5)$", we showed that the model space supporting the most favorable phenomenology should have been probed in 2016 at the LHC2, with an even further reach into this region of the model in 2017-18. This ideal realm of the one-parameter version of No-Scale ${\cal F}$-$SU(5)$ yields a 1.9-2.3 TeV gluino mass at the very same point where the light Higgs boson mass enters its rather narrow experimentally determined range of $m_h = 125.09 \pm 0.24$ GeV. Given the recent results reported at Moriond 2017 for 36 ${\rm fb}^{-1}$ of luminosity collected in 2016 at the 13 TeV LHC2, we now update the status of the No-Scale ${\cal F}$-$SU(5)$ model space in light of the gluino mass exclusion limits presented by the ATLAS and CMS Collaborations. We illustrate that a resolution could be reached soon as to whether supersymmetry lives in this most critical region of the model space.


title: Hierarchy of instabilities for two counter-streaming magnetized pair beams: influence of field obliquity

abstract: The hierarchy of unstable modes when two counter-streaming pair plasmas interact over a flow-aligned magnetic field has been recently investigated [PoP \textbf{23}, 062122 (2016)]. The analysis is here extended to the case of an arbitrarily tilted magnetic field. The two plasma shells are initially cold and identical. For any angle $\theta \in [0,\pi/2]$ between the field and the initial flow, the hierarchy of unstable modes is numerically determined in terms of the initial Lorentz factor of the shells $\gamma_0$, and the field strength as measured by a parameter denoted $\sigma$. For $\theta=0$, four different kinds of mode are likely to lead the linear phase. The hierarchy simplifies for larger $\theta$'s, partly because the Weibel instability can no longer be cancelled in this regime. For $\theta>0.78$ (44$^\circ$) and in the relativistic regime, the Weibel instability always govern the interaction. In the non-relativistic regime, the hierarchy becomes $\theta$-independent because the interaction turns to be field-independent. As a result, the two-stream instability becomes the dominant one, regardless of the field obliquity.


title: Out-of-equilibrium collective oscillation as phonon condensation in a model protein

abstract: In the first part of the present paper (theoretical), the activation of out-of-equilibrium collective oscillations of a macromolecule is described as a classical phonon condensation phenomenon. If a macromolecule is modeled as an open system, that is, it is subjected to an external energy supply and is in contact with a thermal bath to dissipate the excess energy, the internal nonlinear couplings among the normal modes make the system undergo a non-equilibrium phase transition when the energy input rate exceeds a threshold value. This transition takes place between a state where the energy is incoherently distributed among the normal modes, to a state where the input energy is channeled into the lowest frequency mode entailing a coherent oscillation of the entire molecule. The model put forward in the present work is derived as the classical counterpart of a quantum model proposed long time ago by H. Fr\"ohlich in the attempt to explain the huge speed of enzymatic reactions. In the second part of the present paper (experimental), we show that such a phenomenon is actually possible. Two different and complementary THz near-field spectroscopic techniques, a plasmonic rectenna, and a micro-wire near-field probe, have been used in two different labs to get rid of artefacts. By considering a aqueous solution of a model protein, the BSA (Bovine Serum Albumin), we found that this protein displays a remarkable absorption feature around 0.314 THz, when driven in a stationary out-of-thermal equilibrium state by means of optical pumping. The experimental outcomes are in very good qualitative agreement with the theory developed in the first part, and in excellent quantitative agreement with a theoretical result allowing to identify the observed spectral feature with a collective oscillation of the entire molecule.


title: Completeness of Length-Weighted Sobolev Metrics on the Space of Curves

abstract: In this article we prove completeness results for Sobolev metrics with nonconstant coefficients on the space of immersed curves and on the space of unparametrized curves. We provide necessary as well as sufficient conditions for the coefficients of the Riemannian metric for the metric to be metrically complete and we construct examples of incomplete metrics. This work is an extension of previous work on completeness of Sobolev metrics with constant coefficients.


title: Effects of heterogeneity in site-site couplings for tight-binding models on scale-invariant structures

abstract: We studied the thermodynamic behaviors of non-interacting bosons and fermions trapped by a scale-invariant branching structure of adjustable degree of heterogeneity. The full energy spectrum in tight-binding approximation was analytically solved . We found that the log-periodic oscillation of the specific heat for Fermi gas depended on the heterogeneity of hopping. Also, low dimensional Bose-Einstein condensation occurred only for non-homogeneous setup.


title: Exponential decay of connection probabilities for subcritical Voronoi percolation in $\mathbb{R}^d$

abstract: We prove that for Voronoi percolation on $\mathbb{R}^d$, there exists $p_c\in[0,1]$ such that - for $p<p_c$, there exists $c_p>0$ such that $\mathbb{P}_p[0\text{ connected to distance }n]\leq \exp(-c_p n)$, - there exists $c>0$ such that for $p>p_c$, $\mathbb{P}_p[0\text{ connected to }\infty]\geq c(p-p_c)$. For dimension 2, this result offers a new way of showing that $p_c(2)=1/2$. This paper belongs to a series of papers using the theory of algorithms to prove sharpness of the phase transition.


title: Reissner-Nordstr\"om black holes with non-Abelian hair

abstract: We consider $d\geqslant 4$ Einstein--(extended-)Yang-Mills theory, where the gauge sector is augmented by higher order terms. Linearizing the (extended) Yang-Mills equations on the background of the electric Reissner-Nordstr\"om (RN) black hole, we show the existence of normalizable zero modes, dubbed non-Abelian magnetic stationary clouds. The non-linear realization of these clouds bifurcates the RN family into a branch of static, spherically symmetric, electrically charged and asymptotically flat black holes with non-Abelian hair. Generically, the hairy black holes are thermodynamically preferred over the RN solution, which, in this model, becomes unstable against the formation of non-Abelian hair, for sufficiently large values of the electric charge.


title: Predicting stock market movements using network science: An information theoretic approach

abstract: A stock market is considered as one of the highly complex systems, which consists of many components whose prices move up and down without having a clear pattern. The complex nature of a stock market challenges us on making a reliable prediction of its future movements. In this paper, we aim at building a new method to forecast the future movements of Standard & Poor's 500 Index (S&P 500) by constructing time-series complex networks of S&P 500 underlying companies by connecting them with links whose weights are given by the mutual information of 60-minute price movements of the pairs of the companies with the consecutive 5,340 minutes price records. We showed that the changes in the strength distributions of the networks provide an important information on the network's future movements. We built several metrics using the strength distributions and network measurements such as centrality, and we combined the best two predictors by performing a linear combination. We found that the combined predictor and the changes in S&P 500 show a quadratic relationship, and it allows us to predict the amplitude of the one step future change in S&P 500. The result showed significant fluctuations in S&P 500 Index when the combined predictor was high. In terms of making the actual index predictions, we built ARIMA models. We found that adding the network measurements into the ARIMA models improves the model accuracy. These findings are useful for financial market policy makers as an indicator based on which they can interfere with the markets before the markets make a drastic change, and for quantitative investors to improve their forecasting models.

categories: cs.SI cs.CE cs.IT math.IT physics.soc-ph


title: Quantum field theory of a hyper-complex scalar field on a commutative ring

abstract: Inspired by the structural unification of unitary groups (quantum field theory) with orthogonal groups (relativity) proposed recently through a non-division algebra, we construct a hypercomplex field theory with an internal symmetry that unifies the U(1) compact gauge group with the SO(1,1) noncompact gauge group, using the commutative ring of hypercomplex numbers. From the quantum field theory point of view, the hypercomplex field encodes two charged bosons with opposite charge, and corresponds thus to a neutral compound boson. Furthermore, normal or- dering of operators is not required for controling the vacuum divergences; in an analogy with SUSY, the theory under study contains U(1) boson particles and their hyperbolic SO(1,1) boson partners, whose contributions to the vacuum energy cancel out exactly to a zero value. In fact the present scheme allows us to compare finite measuments of squeezed boson-number statistics obtained with and without normal ordering. Additionally we discuss on the potential applications of the squeezed boson states constructed on the commutative ring, in quantum teleportation and in related areas.


title: On Central-Peripheral Appendage Numbers of Uniform Central Graphs

abstract: In a uniform central graph (UCG) the eccentric verticies of a central vertex is the same for all central verticies. This collection of eccentric verticies is the centered periphery. For a pair of graphs $(C, P)$ the central-peripheral appendage number, $A_{ucg}(C, P)$, is the minimum number verticies needed to be adjoined to the graphs $C$ and $P$ in order to construct a uniform central graph H with center C and centered-periphery P. We compute $A_{ucg}(C, P)$ in terms of the radius and diameter of P and whether or not $C$ is a complete graph. In the process we show $A_{ucg}(C, P)\leq 6$ if $diam(P) > 2$. We also provide structure theorems for UCGs in terms of the centered periphery.


title: Liquid Cloud Storage

abstract: A liquid system provides durable object storage based on spreading redundantly generated data across a network of hundreds to thousands of potentially unreliable storage nodes. A liquid system uses a combination of a large code, lazy repair, and a flow storage organization. We show that a liquid system can be operated to enable flexible and essentially optimal combinations of storage durability, storage overhead, repair bandwidth usage, and access performance.


title: Integrals of motion from quantum toroidal algebras

abstract: We identify the Taylor coefficients of the transfer matrices corresponding to quantum toroidal algebras with the elliptic local and non-local integrals of motion introduced by Kojima, Shiraishi, Watanabe, and one of the authors. That allows us to prove the Litvinov conjectures on the Intermediate Long Wave model. We also discuss the (gl(m),gl(n)) duality of XXZ models in quantum toroidal setting and the implications for the quantum KdV model. In particular, we conjecture that the spectrum of non-local integrals of motion of Bazhanov, Lukyanov, and Zamolodchikov is described by Gaudin Bethe ansatz equations associated to affine sl(2).


title: Capacitated Bounded Cardinality Hub Routing Problem: Model and Solution Algorithm

abstract: In this paper, we address the Bounded Cardinality Hub Location Routing with Route Capacity wherein each hub acts as a transshipment node for one directed route. The number of hubs lies between a minimum and a maximum and the hub-level network is a complete subgraph. The transshipment operations take place at the hub nodes and flow transfer time from a hub-level transporter to a spoke-level vehicle influences spoke- to-hub allocations. We propose a mathematical model and a branch-and-cut algorithm based on Benders decomposition to solve the problem. To accelerate convergence, our solution framework embeds an efficient heuristic producing high-quality solutions in short computation times. In addition, we show how symmetry can be exploited to accelerate and improve the performance of our method.


title: A consistent measure of the merger histories of massive galaxies using close-pair statistics I: Major mergers at $z < 3.5$

abstract: We use a large sample of $\sim 350,000$ galaxies constructed by combining the UKIDSS UDS, VIDEO/CFHT-LS, UltraVISTA/COSMOS and GAMA survey regions to probe the major merging histories of massive galaxies ($>10^{10}\ \mathrm{M}_\odot$) at $0.005 < z < 3.5$. We use a method adapted from that presented in Lopez-Sanjuan et al. (2014) using the full photometric redshift probability distributions, to measure pair $\textit{fractions}$ of flux-limited, stellar mass selected galaxy samples using close-pair statistics. The pair fraction is found to weakly evolve as $\propto (1+z)^{0.8}$ with no dependence on stellar mass. We subsequently derive major merger $\textit{rates}$ for galaxies at $> 10^{10}\ \mathrm{M}_\odot$ and at a constant number density of $n > 10^{-4}$ Mpc$^{-3}$, and find rates a factor of 2-3 smaller than previous works, although this depends strongly on the assumed merger timescale and likelihood of a close-pair merging. Galaxies undergo approximately 0.5 major mergers at $z < 3.5$, accruing an additional 1-4 $\times 10^{10}\ \mathrm{M}_\odot$ in the process. Major merger accretion rate densities of $\sim 2 \times 10^{-4}$ $\mathrm{M}_\odot$ yr$^{-1}$ Mpc$^{-3}$ are found for number density selected samples, indicating that direct progenitors of local massive ($>10^{11}\mathrm{M}_\odot$) galaxies have experienced a steady supply of stellar mass via major mergers throughout their evolution. While pair fractions are found to agree with those predicted by the Henriques et al. (2014) semi-analytic model, the Illustris hydrodynamical simulation fails to quantitatively reproduce derived merger rates. Furthermore, we find major mergers become a comparable source of stellar mass growth compared to star-formation at $z < 1$, but is 10-100 times smaller than the SFR density at higher redshifts.


title: Multivariate generalized Pareto distributions: parametrizations, representations, and properties

abstract: Multivariate generalized Pareto distributions arise as the limit distributions of exceedances over multivariate thresholds of random vectors in the domain of attraction of a max-stable distribution. These distributions can be parametrized and represented in a number of different ways. Moreover, generalized Pareto distributions enjoy a number of interesting stability properties. An overview of the main features of such distributions are given, expressed compactly in several parametrizations, giving the potential user of these distributions a convenient catalogue of ways to handle and work with generalized Pareto distributions.


title: Aqua MODIS Electronic Crosstalk on SMWIR Bands 20 to 26

abstract: Aqua MODIS Moon images obtained with bands 20 to 26 (3.66 - 4.55 and 1.36 - 1.39 $\mu$m) during scheduled lunar events show evidence of electronic crosstalk contamination of the response of detector 1. In this work, we determined the sending bands for each receiving band. We found that the contaminating signal originates, in all cases, from the detector 10 of the corresponding sending band and that the signals registered by the receiving and sending detectors are always read out in immediate sequence. We used the lunar images to derive the crosstalk coefficients, which were then applied in the correction of electronic crosstalk striping artifacts present in L1B images, successfully restoring product quality.


title: Baroclinic Instability and Transitions in a Two-Layer Quasigeostrophic Channel Model

abstract: The main objective of this article is to derive a mathematical theory associated with the nonlinear stability and dynamic transitions of the basic shear flows associated with baroclinic instability, which plays a fundamental role in the dominant mechanism shaping the cyclones and anticyclones that dominate weather in mid-latitudes, as well as the mesoscale eddies that play various roles in oceanic dynamics and the transport of tracers. This article provides a general transition and stability theory for the two-layer quasi-geostrophic model originally derived by Pedlosky \cite{ped1970}. We show that the instability and dynamic transition of the basic shear flow occur only when the Froude number $F>(\gamma^2+1)/2$, where $\gamma$ is the weave number of lowest zonal harmonic. In this case, we derive a precise critical shear velocity $U_c$ and a dynamic transition number $b$, and we show that the dynamic transition and associated flo0w patterns are dictated by the sign of this transition number. In particular, we show that for $b>0$, the system undergoes a continuous transition, leading to spatiotemporal flow oscillatory patterns as the shear velocity $U$ crosses $U_c$; for $b<0$, the system undergoes a jump transition, leading to drastic change and more complex transition patterns in the far field, and to unstable period solutions for $U<U_c$. Furthermore, we show that when the wavenumber of the zonal harmonic $\gamma=1$, the transition number $b$ is always positive and the system always undergoes a continuous dynamic transition leading spatiotemporal oscillations. This suggests that a continuous transition to spatiotemporal patterns is preferable for the shear flow associated with baroclinic instability.


title: Aqueous solution interactions with sex hormone-binding globulin and estradiol: A theoretical investigation

abstract: Sex hormone-binding globulin (SHBG) is a binding protein that regulates availability of steroids hormones in the plasma. Although best known as steroid carrier, studies have associated SHBG in modulating behavioral aspects related to sexual receptivity. Among steroids, estradiol (17\b{eta}-estradiol, oestradiol or E2) is well recognized as the most active endogenous female hormone, exerting important roles in reproductive and nonreproductive functions. Thus, in this study we aimed to employ molecular dynamics (MD) and docking techniques for quantifying the interaction energy between a complex aqueous solution, composed by different salts, SHBG and E2. Due to glucose concentration resembles those observed in diabetic levels, special emphasis was devoted to uncover the main consequences of this carbohydrate on the SHBG and E2 molecules. We also examined possible energetic changes due to solution on the binding energy of SHBG-E2 complex. In this framework, our calculations uncovered a remarkable interaction energy between glucose and SHBG surface. Surprisingly, we also observed solute components movement toward SHBG yielding clusters surrounding the protein. This finding, corroborated by the higher energy and shorter distance found between glucose and SHBG, suggests a scenario in favor of a detainment state. In addition, in spite of protein superficial area increment it does not exerted modification on binding site area nor over binding energy SHBG-E2 complex. Finally, our calculations also highlighted an interaction between E2 and glucose when the hormone was immersed in the solution. In summary, our findings contribute for a better comprehension of both SHBG and E2 interplay with aqueous solution components.


title: Quadratic obstructions to small-time local controllability for scalar-input differential systems

abstract: We consider nonlinear scalar-input differential control systems in the vicinity of an equilibrium. When the linearized system at the equilibrium is controllable, the nonlinear system is smoothly small-time locally controllable, i.e., whatever $m>0$ and $T>0$, the state can reach a whole neighborhood of the equilibrium at time $T$ with controls arbitrary small in $C^m$-norm. When the linearized system is not controllable, we prove that small-time local controllability cannot be recovered from the quadratic expansion and that the following quadratic alternative holds. Either the state is constrained to live within a smooth strict invariant manifold, up to a cubic residual, or the quadratic order adds a signed drift in the evolution with respect to this manifold. In the second case, the quadratic drift holds along an explicit Lie bracket of length $(2k+1)$, it is quantified in terms of an $H^{-k}$-norm of the control, it holds for controls small in $W^{2k,\infty}$-norm. These spaces are optimal for general nonlinear systems and are slightly improved in the particular case of control-affine systems. Unlike other works based on Lie-series formalism, our proof is based on an explicit computation of the quadratic terms by means of appropriate transformations. In particular, it does not require that the vector fields defining the dynamic are smooth. We prove that $C^3$ regularity is sufficient for our alternative to hold. This work underlines the importance of the norm used in the smallness assumption on the control: depending on this choice of functional setting, the same system may or may not be small-time locally controllable, even though the state lives within a finite dimensional space.


title: Probing flavor parameters in the scalar sector and new bounds for the fermion sector

abstract: In this paper we study the Flavor-Changing mediated by Higgs boson within Two-Higgs Doublet Model-III context. We explore the parameter space and, considering recent results, find new limits for the model parameters. We also obtain the total Higgs decay width and the branching ratios for different channels taking limits for the $t\to cV$ and $b\to s\gamma$ processes. Considering different constraints, we estimate the branching ratio for the $h\to t^*_{}c$ in the model as well as the bounds for the $b \to s\gamma$, $h \to \mu\tau$ and $h\to \gamma Z.$ Considering the quark top decays to $W\bar{b},$ we obtained $Br(h \to Wc\bar{b})\sim10^{-3}$.


title: Fair Allocation based on Diminishing Differences

abstract: Ranking alternatives is a natural way for humans to explain their preferences. It is being used in many settings, such as school choice, course allocations and residency matches. In some cases, several `items' are given to each participant. Without having any information on the underlying cardinal utilities, arguing about fairness of allocation requires extending the ordinal item ranking to ordinal bundle ranking. The most commonly used such extension is stochastic dominance (SD), where a bundle X is preferred over a bundle Y if its score is better according to all additive score functions. SD is a very conservative extension, by which few allocations are necessarily fair while many allocations are possibly fair. We propose to make a natural assumption on the underlying cardinal utilities of the players, namely that the difference between two items at the top is larger than the difference between two items at the bottom. This assumption implies a preference extension which we call diminishing differences (DD), where X is preferred over Y if its score is better according to all additive score functions satisfying the DD assumption. We give a full characterization of allocations that are necessarily-proportional or possibly-proportional according to this assumption. Based on this characterization, we present a polynomial-time algorithm for finding a necessarily-DD-proportional allocation if it exists. Using simulations, we show that with high probability, a necessarily-proportional allocation does not exist but a necessarily-DD-proportional allocation exists, and moreover, that allocation is proportional according to the underlying cardinal utilities. We also consider chore allocation under the analogous condition --- increasing-differences.


title: This is not the feedback you have been looking for: nearby optical AGN rarely drive kpc-scale cold-gas outflows

abstract: We study the interstellar Na I $\lambda \lambda 5890, 5895$ (Na D) absorption-line doublet in a nearly-complete sample of $\sim$9900 nearby Seyfert 2 galaxies, in order to quantify the significance of optical AGN activity in driving kpc-scale outflows that can quench star formation. Comparison to a carefully matched sample of $\sim$44,000 control objects indicates that the Seyfert and control population have similar Na D detection rates ($\sim 5-6%$). Only 53 Seyferts (or 0.5% of the population) are found to potentially display galactic-scale winds, compared to 0.8% of the control galaxies. While nearly a third of the Na D outflows observed in our Seyfert 2 galaxies occur around the brightest AGN, both radio and infrared data indicate that star formation could play the dominant role in driving cold-gas outflows in an even higher fraction of the Na D-outflowing Seyfert 2s. Our results indicate that galactic-scale outflows at low redshift are no more frequent in Seyferts than they are in their non-active counterparts, that optical AGN are not significant contributors to the quenching of star formation in the nearby Universe, and that star-formation may actually be the principal driver of outflows even in systems that do host an AGN.


title: Theorem of the heart in negative K-theory for weight structures

abstract: We construct the strong weight complex functor (in the sense of Bondarko) for a stable infinity-category $\underline{C}$ equipped with a bounded weight structure $w$. Along the way we prove that $\underline{C}$ is determined by the infinity-categorical heart of $w$. This allows us to compare the K-theory of $\underline{C}$ and the K-theory of $\underline{Hw}$, the classical heart of $w$. In particular, we prove that $\operatorname{K}_{n}(\underline{C}) \to \operatorname{K}_{n}(\underline{Hw})$ are isomorphisms for $n \le 0$.


title: Living Together: Mind and Machine Intelligence

abstract: In this paper we consider the nature of the machine intelligences we have created in the context of our human intelligence. We suggest that the fundamental difference between human and machine intelligence comes down to \emph{embodiment factors}. We define embodiment factors as the ratio between an entity's ability to communicate information vs compute information. We speculate on the role of embodiment factors in driving our own intelligence and consciousness. We briefly review dual process models of cognition and cast machine intelligence within that framework, characterising it as a dominant System Zero, which can drive behaviour through interfacing with us subconsciously. Driven by concerns about the consequence of such a system we suggest prophylactic courses of action that could be considered. Our main conclusion is that it is \emph{not} sentient intelligence we should fear but \emph{non-sentient} intelligence.


title: Permutation Tests for Infection Graphs

abstract: We formulate and analyze a novel hypothesis testing problem for inferring the edge structure of an infection graph. In our model, a disease spreads over a network via contagion or random infection, where the random variables governing the rates of contracting the disease from neighbors or random infection are independent exponential random variables with unknown rate parameters. A subset of nodes is also censored uniformly at random. Given the statuses of nodes in the network, the goal is to determine the underlying graph. We present a procedure based on permutation testing, and we derive sufficient conditions for the validity of our test in terms of automorphism groups of the graphs corresponding to the null and alternative hypotheses. Further, the test is valid more generally for infection processes satisfying a basic symmetry condition. Our test is easy to compute and does not involve estimating unknown parameters governing the process. We also derive risk bounds for our permutation test in a variety of settings, and motivate our test statistic in terms of approximate equivalence to likelihood ratio testing and maximin tests. We conclude with an application to real data from an HIV infection network.


title: Synaptic Noise Facilitates the Emergence of Self-Organized Criticality in the Caenorhabditis elegans Neuronal Network

abstract: Avalanches with power-law distributed size parameters have been observed in neuronal networks. This observation might be a manifestation of the self-organized criticality (SOC). Yet, the physiological mechanicsm of this behavior is currently unknown. Describing synaptic noise as transmission failures mainly originating from the probabilistic nature of neurotransmitter release, this study investigates the potential of this noise as a mechanism for driving the functional architecture of the neuronal networks towards SOC. To this end, a simple finite state neuron model, with activity dependent and synapse specific failure probabilities, was built based on the known anatomical connectivity data of the nematode Ceanorhabditis elegans. Beginning from random values, it was observed that synaptic noise levels picked out a set of synapses and consequently an active subnetwork which generates power-law distributed neuronal avalanches. The findings of this study brings up the possibility that synaptic failures might be a component of physiological processes underlying SOC in neuronal networks.


title: GP-Unet: Lesion Detection from Weak Labels with a 3D Regression Network

abstract: We propose a novel convolutional neural network for lesion detection from weak labels. Only a single, global label per image - the lesion count - is needed for training. We train a regression network with a fully convolutional architecture combined with a global pooling layer to aggregate the 3D output into a scalar indicating the lesion count. When testing on unseen images, we first run the network to estimate the number of lesions. Then we remove the global pooling layer to compute localization maps of the size of the input image. We evaluate the proposed network on the detection of enlarged perivascular spaces in the basal ganglia in MRI. Our method achieves a sensitivity of 62% with on average 1.5 false positives per image. Compared with four other approaches based on intensity thresholding, saliency and class maps, our method has a 20% higher sensitivity.


title: Concurrent Channel Probing and Data Transmission in Full-duplex MIMO Systems

abstract: An essential step for achieving multiplexing gain in MIMO downlink systems is to collect accurate channel state information (CSI) from the users. Traditionally, CSIs have to be collected before any data can be transmitted. Such a sequential scheme incurs a large feedback overhead, which substantially limits the multiplexing gain especially in a network with a large number of users. In this paper, we propose a novel approach to mitigate the feedback overhead by leveraging the recently developed Full-duplex radios. Our approach is based on the key observation that using Full-duplex radios, when the base-station (BS) is collecting CSI of one user through the uplink channel, it can use the downlink channel to simultaneously transmit data to other (non-interfering) users for which CSIs are already known. By allowing concurrent channel probing and data transmission, our scheme can potentially achieve a higher throughput compared to traditional schemes using Half-duplex radios. The new flexibility introduced by our scheme, however, also leads to fundamental challenges in achieving throughout optimal scheduling. In this paper, we make an initial effort to this important problem by considering a simplified group interference model. We develop a throughput optimal scheduling policy with complexity $O((N/I)^I)$, where $N$ is the number of users and $I$ is the number of user groups. To further reduce the complexity, we propose a greedy policy with complexity $O(N\log N)$ that not only achieves at least 2/3 of the optimal throughput region, but also outperforms any feasible Half-duplex solutions. We derive the throughput gain offered by Full-duplex under different system parameters and show the advantage of our algorithms through numerical studies.


title: Grouped multivariate and functional time series forecasting: An application to annuity pricing

abstract: Age-specific mortality rates are often disaggregated by different attributes, such as sex, state, ethnic group and socioeconomic status. In making social policies and pricing annuity at national and subnational levels, it is important not only to forecast mortality accurately, but also to ensure that forecasts at the subnational level add up to the forecasts at the national level. This motivates recent developments in grouped functional time series methods (Shang and Hyndman, 2017) to reconcile age-specific mortality forecasts. We extend these grouped functional time series forecasting methods to multivariate time series, and apply them to produce point forecasts of mortality rates at older ages, from which fixed-term annuities for different ages and maturities can be priced. Using the regional age-specific mortality rates in Japan obtained from the Japanese Mortality Database, we investigate the one-step-ahead to 15-step-ahead point-forecast accuracy between the independent and grouped forecasting methods. The grouped forecasting methods are shown not only to be useful for reconciling forecasts of age-specific mortality rates at national and subnational levels, but they are also shown to allow improved forecast accuracy. The improved forecast accuracy of mortality rates is of great interest to the insurance and pension industries for estimating annuity prices, in particular at the level of population subgroups, defined by key factors such as sex, region, and socioeconomic grouping.


title: Sagittarius A* High Energy X-ray Flare Properties During NuSTAR Monitoring of the Galactic Center from 2012 to 2015

abstract: Understanding the origin of the flaring activity from the Galactic center supermassive black hole, Sagittarius A*, is a major scientific goal of the NuSTAR Galactic plane survey campaign. We report on the data obtained between July 2012 and April 2015, including 27 observations on Sgr A* with a total exposure of ~ 1 Ms. We found a total of ten X-ray flares detected in the NuSTAR observation window, with luminosities in the range of $L_{3-79~keV}$~$(0.2-4.0) \times 10^{35}~erg~s^{-1}$. With this largest hard X-ray Sgr A* flare dataset to date, we studied the flare spectral properties. Seven flares are detected above 5{\sigma} significance, showing a range of photon indices ({\Gamma} ~ 2.0-2.8) with typical uncertainties of +/-0.5 (90% confidence level). We found no significant spectral hardening for brighter flares as indicated by a smaller sample. The accumulation of all the flare spectra in 1-79 keV can be well fit with an absorbed power-law model with {\Gamma}=2.2+/-0.1, and does not require the existence of a spectral break. The lack of variation in X-ray spectral index with luminosity would point to a single mechanism for the flares and is consistent with the synchrotron scenario. Lastly, we present the quiescent state spectrum of Sgr A*, and derived an upper limit on the quiescent luminosity of Sgr A* above 10 keV to be $L_{Xq, 10-79 keV}$ < $(2.9{\pm}0.2) \times 10^{34}~erg~s^{-1}$.


title: High-resolution Observation of Downflows at One End of a Pre-eruption Filament

abstract: Studying the dynamics of filaments at pre-eruption phase can shed light on the precursor of eruptive events. Such studies in high-resolution (in the order of 0.1") are highly desirable yet very rare so far. In this work, we present a detailed observation of a pre-eruption evolution of a filament obtained by the 1.6 m New Solar Telescope (NST) at Big Bear Solar Observatory (BBSO). One end of the filament is anchored at the sunspot in NOAA active region (AR) 11515, which is well observed by NST H$\alpha$ off-bands four hours before till one hour after the filament eruption. A M1.6 flare is associated with the eruption. We observed persistent downflowing materials along the H$\alpha$ multi-threaded component of the loop towards the AR end during the pre-eruption phase. We traced the trajectories of plasma blobs along the H$\alpha$ threads and obtained the plane-of-sky velocity of 45 km s$^{-1}$ on average. We further estimated the real velocities of the downflows and the altitude of the filament by matching the observed H$\alpha$ threads with magnetic field lines extrapolated from a nonlinear force-free field (NLFFF) model. Observation of chromospheric brightenings (BZs) at the footpoints of the falling plasma blobs is also presented in the paper. The lower limit of the kinetic energy per second of the downflows through the BZs is found to be $\sim$ 10$^{21}$ erg. Larger FOV observations from BBSO full disk H$\alpha$ images show that the AR end of the filament started ascending four hours before the flare. We attribute the observed downflows at the AR end of the filament to the draining effect of the filament rising prior to its eruption. During the slow-rise phase, the downflows continuously drained away $\sim$ 10$^{15}$g mass from the filament over a few hours, which is believed to be essential for the instability at last, and could be an important precursor of eruptive events.


title: Heavy-Heavy-Light-Light correlators in Liouville theory

abstract: We compute four-point functions of two heavy and two "perturbatively heavy" operators in the semiclassical limit of Liouville theory on the sphere. We obtain these "Heavy-Heavy-Light-Light" (HHLL) correlators to leading order in the conformal weights of the light insertions in two ways: (a) via a path integral approach, combining different methods to evaluate correlation functions from complex solutions for the Liouville field, and (b) via the conformal block expansion. This latter approach identifies an integral over the continuum of normalizable states and a sum over an infinite tower of lighter discrete states, whose contribution we extract by analytically continuing standard results to our HHLL setting. The sum over this tower reproduces the sum over those complex saddlepoints of the path integral that contribute to the correlator. Our path integral computations reveal that when the two light operators are inserted at equal time in radial quantization, the leading-order HHLL correlator is independent of their separation, and more generally that at this order there is no short-distance singularity as the two light operators approach each other. The conformal block expansion likewise shows that in the discrete sum short-distance singularities are indeed absent for all intermediate states that contribute. In particular, the Virasoro vacuum block, which would have been singular at short distances, is not exchanged. The separation-independence of equal-time correlators is due to cancelations between the discrete contributions. These features lead to a Lorentzian singularity that, in conformal theories with anti-de Sitter (AdS) duals, would be associated to locality below the AdS scale.


title: The extension of some D(4)-pairs

abstract: In this paper we illustrate the use of the results from [1] proving that $D(4)$-triple $\{a, b, c\}$ with $a < b < a + 57\sqrt{a}$ has a unique extension to a quadruple with a larger element. This furthermore implies that $D(4)$-pair $\{a, b\}$ cannot be extended to a quintuple if $a < b < a + 57\sqrt{a}$.


title: Learning the Morphology of Brain Signals Using Alpha-Stable Convolutional Sparse Coding

abstract: Neural time-series data contain a wide variety of prototypical signal waveforms (atoms) that are of significant importance in clinical and cognitive research. One of the goals for analyzing such data is hence to extract such 'shift-invariant' atoms. Even though some success has been reported with existing algorithms, they are limited in applicability due to their heuristic nature. Moreover, they are often vulnerable to artifacts and impulsive noise, which are typically present in raw neural recordings. In this study, we address these issues and propose a novel probabilistic convolutional sparse coding (CSC) model for learning shift-invariant atoms from raw neural signals containing potentially severe artifacts. In the core of our model, which we call $\alpha$CSC, lies a family of heavy-tailed distributions called $\alpha$-stable distributions. We develop a novel, computationally efficient Monte Carlo expectation-maximization algorithm for inference. The maximization step boils down to a weighted CSC problem, for which we develop a computationally efficient optimization algorithm. Our results show that the proposed algorithm achieves state-of-the-art convergence speeds. Besides, $\alpha$CSC is significantly more robust to artifacts when compared to three competing algorithms: it can extract spike bursts, oscillations, and even reveal more subtle phenomena such as cross-frequency coupling when applied to noisy neural time series.