title: Output-sensitive algorithms for sumset and sparse polynomial multiplication

abstract: We present randomized algorithms to compute the sumset (Minkowski sum) of two integer sets, and to multiply two univariate integer polynomials given by sparse representations. Our algorithm for sumset has cost softly linear in the combined size of the inputs and output. This is used as part of our sparse multiplication algorithm, whose cost is softly linear in the combined size of the inputs, output, and the sumset of the supports of the inputs. As a subroutine, we present a new method for computing the coefficients of a sparse polynomial, given a set containing its support. Our multiplication algorithm extends to multivariate Laurent polynomials over finite fields and rational numbers. Our techniques are based on sparse interpolation algorithms and results from analytic number theory.


title: Multi Stage based Time Series Analysis of User Activity on Touch Sensitive Surfaces in Highly Noise Susceptible Environments

abstract: This article proposes a multistage framework for time series analysis of user activity on touch sensitive surfaces in noisy environments. Here multiple methods are put together in multi stage framework; including moving average, moving median, linear regression, kernel density estimation, partial differential equations and Kalman filter. The proposed three stage filter consisting of partial differential equation based denoising, Kalman filter and moving average method provides ~25% better noise reduction than other methods according to Mean Squared Error (MSE) criterion in highly noise susceptible environments. Apart from synthetic data, we also obtained real world data like hand writing, finger/stylus drags etc. on touch screens in the presence of high noise such as unauthorized charger noise or display noise and validated our algorithms. Furthermore, the proposed algorithm performs qualitatively better than the existing solutions for touch panels of the high end hand held devices available in the consumer electronics market qualitatively.


title: Global root bracketing method with adaptive mesh refinement

abstract: An efficient method for finding all real roots of a univariate function in a given bounded domain is formulated. The proposed method uses adaptive mesh refinement to locate bracketing intervals based on bisection criterion for root finding. Each bracketing interval encloses one root. An adaptive form of error is introduced to enclose roots in a desired tolerance based on how much close the roots are. Detecting roots with even multiplicity, which is regarded out of the realm of bracketing methods, becomes possible with the method proposed in this paper. Also, strategies for finding odd-multiple roots with the least number of function evaluations is proposed. Adaptive mesh refinement lead to considerable reduction in function evaluations in comparison to previous global root bracketing methods. The reliability of the proposed method is being illustrated by several examples.


title: Kiloparsec-scale jets in three radio-loud narrow-line Seyfert 1 galaxies

abstract: We have discovered kiloparsec-scale extended radio emission in three narrow-line Seyfert 1 galaxies (NLS1s) in sub-arcsecond resolution 9 GHz images from the Karl G. Jansky Very Large Array (VLA). We find all sources show two-sided, mildly core-dominated jet structures with diffuse lobes dominated by termination hotspots. These span 20-70 kpc with morphologies reminiscent of FR II radio galaxies, while the extended radio luminosities are intermediate between FR I and FR II sources. In two cases the structure is linear, while a $45^{\circ}$ bend is apparent in the third. Very Long Baseline Array images at 7.6 GHz reveal parsec-scale jet structures, in two cases with extended structure aligned with the inner regions of the kiloparsec-scale jets. Based on this alignment, the ratio of the radio core luminosity to the optical luminosity, the jet/counter-jet intensity and extension length ratios, and moderate core brightness temperatures ($\lesssim10^{10}$ K), we conclude these jets are mildly relativistic ($\beta\lesssim0.3$, $\delta\sim1$-$1.5$) and aligned at moderately small angles to the line of sight (10-15$^{\circ}$). The derived kinematic ages of $\sim10^6$-$10^7$ y are much younger than radio galaxies but comparable to other NLS1s. Our results increase the number of radio-loud NLS1s with known kiloparsec-scale extensions from seven to ten and suggest that such extended emission may be common, at least among the brightest of these sources.


title: The Spatial Structure of Young Stellar Clusters. II. Total Young Stellar Populations

abstract: We investigate the intrinsic stellar populations (estimated total numbers of OB and pre-main-sequence stars down to 0.1 Mo) that are present in 17 massive star-forming regions (MSFRs) surveyed by the MYStIX project. The study is based on the catalog of >31,000 MYStIX Probable Complex Members with both disk-bearing and disk-free populations, compensating for extinction, nebulosity, and crowding effects. Correction for observational sensitivities is made using the X-ray Luminosity Function (XLF) and the near-infrared Initial Mass Function (IMF)--a correction that is often not made by infrared surveys of young stars. The resulting maps of the projected structure of the young stellar populations, in units of intrinsic stellar surface density, allow direct comparison between different regions. Several regions have multiple dense clumps, similar in size and density to the Orion Nebula Cluster. The highest projected density of ~34,000 stars/pc^2 is found in the core of the RCW38 cluster. Histograms of surface density show different ranges of values in different regions, supporting the conclusion of Bressert et al. (2010, B10) that no universal surface-density threshold can distinguish between clustered and distributed star-formation. However, a large component of the young stellar population of MSFRs resides in dense environments of 200-10,000 stars/pc^2 (including within the nearby Orion molecular clouds), and we find that there is no evidence for the B10 conclusion that such dense regions form an extreme "tail" of the distribution. Tables of intrinsic populations for these regions are used in our companion study of young cluster properties and evolution.


title: The Origin of the Extragalactic Gamma-Ray Background and Implications for Dark-Matter Annihilation

abstract: The origin of the extragalactic $\gamma$-ray background (EGB) has been debated for some time. { The EGB comprises the $\gamma$-ray emission from resolved and unresolved extragalactic sources, such as blazars, star-forming galaxies and radio galaxies, as well as radiation from truly diffuse processes.} This letter focuses on the blazar source class, the most numerous detected population, and presents an updated luminosity function and spectral energy distribution model consistent with the blazar observations performed by the {\it Fermi} Large Area Telescope (LAT). We show that blazars account for 50$^{+12}_{-11}$\,\% of the EGB photons ($>$0.1\,GeV), and that {\it Fermi}-LAT has already resolved $\sim$70\,\% of this contribution. Blazars, and in particular low-luminosity hard-spectrum nearby sources like BL Lacs, are responsible for most of the EGB emission above 100\,GeV. We find that the extragalactic background light, which attenuates blazars' high-energy emission, is responsible for the high-energy cut-off observed in the EGB spectrum. Finally, we show that blazars, star-forming galaxies and radio galaxies can naturally account for the amplitude and spectral shape of the background in the 0.1--820\,GeV range, leaving only modest room for other contributions. This allows us to set competitive constraints on the dark-matter annihilation cross section.


title: Driving extreme variability: The evolving corona and evidence for jet launching in Markarian 335

abstract: Variations in the X-ray emission from the narrow line Seyfert 1 galaxy, Markarian 335 (Mrk 335), are studied on both long and short timescales through observations made between 2006 and 2013 with XMM-Newton, Suzaku and NuSTAR. Changes in the geometry and energetics of the corona that give rise to this variability are inferred through measurements of the relativistically blurred reflection seen from the accretion disc. On long timescales, we find that during the high flux epochs the corona has expanded, covering the inner regions of the accretion disc out to a radius of 26(-7,+10)rg. The corona contracts to within 12rg and 5rg in the intermediate and low flux epochs, respectively. While the earlier high flux observation made in 2006 is consistent with a corona extending over the inner part of the accretion disc, a later high flux observation that year revealed that the X-ray source had become collimated into a vertically-extended jet-like corona and suggested relativistic motion of material upward. On short timescales, we find that an X-ray flare during a low flux epoch in 2013 corresponded to a reconfiguration from a slightly extended corona to one much more compact, within just 2~3rg of the black hole. There is evidence that during the flare itself, the spectrum softened and the corona became collimated and slightly extended vertically as if a jet-launching event was aborted. Understanding the evolution of the X-ray emitting corona may reveal the underlying mechanism by which the luminous X-ray sources in AGN are powered.


title: Cosmic Web Reconstruction through Density Ridges: Method and Algorithm

abstract: The detection and characterization of filamentary structures in the cosmic web allows cosmologists to constrain parameters that dictates the evolution of the Universe. While many filament estimators have been proposed, they generally lack estimates of uncertainty, reducing their inferential power. In this paper, we demonstrate how one may apply the Subspace Constrained Mean Shift (SCMS) algorithm (Ozertem and Erdogmus (2011); Genovese et al. (2012)) to uncover filamentary structure in galaxy data. The SCMS algorithm is a gradient ascent method that models filaments as density ridges, one-dimensional smooth curves that trace high-density regions within the point cloud. We also demonstrate how augmenting the SCMS algorithm with bootstrap-based methods of uncertainty estimation allows one to place uncertainty bands around putative filaments. We apply the SCMS method to datasets sampled from the P3M N-body simulation, with galaxy number densities consistent with SDSS and WFIRST-AFTA and to LOWZ and CMASS data from the Baryon Oscillation Spectroscopic Survey (BOSS). To further assess the efficacy of SCMS, we compare the relative locations of BOSS filaments with galaxy clusters in the redMaPPer catalog, and find that redMaPPer clusters are significantly closer (with p-values $< 10^{-9}$) to SCMS-detected filaments than to randomly selected galaxies.


abstract: New telescopes like the Square Kilometre Array (SKA) will push into a new sensitivity regime and expose systematics, such as direction-dependent effects, that could previously be ignored. Current methods for handling such systematics rely on alternating best estimates of instrumental calibration and models of the underlying sky, which can lead to inadequate uncertainty estimates and biased results because any correlations between parameters are ignored. These deconvolution algorithms produce a single image that is assumed to be a true representation of the sky, when in fact it is just one realization of an infinite ensemble of images compatible with the noise in the data. In contrast, here we report a Bayesian formalism that simultaneously infers both systematics and science. Our technique, Bayesian Inference for Radio Observations (BIRO), determines all parameters directly from the raw data, bypassing image-making entirely, by sampling from the joint posterior probability distribution. This enables it to derive both correlations and accurate uncertainties, making use of the flexible software MEQTREES to model the sky and telescope simultaneously. We demonstrate BIRO with two simulated sets of Westerbork Synthesis Radio Telescope data sets. In the first, we perform joint estimates of 103 scientific (flux densities of sources) and instrumental (pointing errors, beamwidth and noise) parameters. In the second example, we perform source separation with BIRO. Using the Bayesian evidence, we can accurately select between a single point source, two point sources and an extended Gaussian source, allowing for 'super-resolution' on scales much smaller than the synthesized beam.


title: Assembling Fibonacci Anyons From a $\mathbb{Z}_3$ Parafermion Lattice Model

abstract: Recent concrete proposals suggest it is possible to engineer a two-dimensional bulk phase supporting non-Abelian Fibonacci anyons out of Abelian fractional quantum Hall systems. The low-energy degrees of freedom of such setups can be modeled as $\mathbb{Z}_3$ parafermions "hopping" on a two-dimensional lattice. We use the density matrix renormalization group to study a model of this type interpolating between the decoupled-chain, triangular-lattice, and square-lattice limits. The results show clear evidence of the Fibonacci phase over a wide region of the phase diagram, most notably including the isotropic triangular lattice point. We also study the broader phase diagram of this model and show that elsewhere it supports an Abelian state with semionic excitations.


title: A Dark Year for Tidal Disruption Events

abstract: The disruption of a main-sequence star by a supermassive black hole results in the initial production of an extended debris stream that winds repeatedly around the black hole, producing a complex three-dimensional figure that may self-intersect. Both analytical work and simulations have shown that typical encounters generate streams that are extremely thin. In this paper we show that this implies that even small relativistic precessions attributed to black hole spin can induce deflections that prevent the stream from self-intersecting even after many windings. Additionally, hydrodynamical simulations have demonstrated that energy is deposited very slowly via hydrodynamic processes alone, resulting in the liberation of very little gravitational binding energy in the absence of stream-stream collisions. This naturally leads to a "dark period" in which the flare is not observable for some time, persisting for up to a dozen orbital periods of the most bound material, which translates to years for disruptions around black holes with mass $\sim 10^{7} M_{\odot}$. We find that more-massive black holes tend to have more violent stream self-intersections, resulting in short viscous times that lead to prompt accretion onto the black hole. For these tidal disruption events (TDEs), the accretion rate onto the black hole should still closely follow the original fallback rate after a fixed delay time $t_{\rm delay}$. For lower black hole masses ($M_{\rm h} \lesssim 10^{6}$), we find that flares are typically slowed down by about an order of magnitude, and because the accretion rates for TDEs about higher-mass black hole are already sub-Eddington, this results in the majority of TDEs being sub-Eddington at peak. This also implies that current searches for TDEs are biased towards prompt flares, with slowed flares likely having been unidentified. [abridged]


title: Tasting the SU(5) nature of supersymmetry at the LHC

abstract: We elaborate on a recently found SU(5) relation confined to the up-(s)quark flavour space, that remains immune to large quantum corrections up to the TeV scale. We investigate the possibilities opened by this new window on the GUT scale in order to find TeV-scale SU(5) tests realizable at the LHC. These SU(5) tests appear as relations among observables involving either flavour violation or chirality flip in the up-(s)quark sector. The power of these tests is systematically evaluated using a frequentist, p-value based criterion. SU(5) tests in the Heavy supersymmetry (SUSY), Natural supersymmetry and Top-charm supersymmetry spectra are investigated. The latter scenario features light stops and scharms and is well-motivated from various five-dimensional constructions. A variety of SU(5) tests is obtained, involving techniques of top polarimetry, charm-tagging, or Higgs detection from SUSY cascade decays. We find that O(10) to O(100) events are needed to obtain 50% of relative precision at 3-sigma significance for all of these tests. In addition, we propose a set of precision measurements in ultraperipheral collisions in order to search for the flavour-changing dipole operators of Heavy supersymmetry.


title: Gamma-ray novae as probes of relativistic particle acceleration at non-relativistic shocks

abstract: The Fermi LAT discovery that classical novae produce >100 MeV gamma-rays establishes that shocks and relativistic particle acceleration are key features of these events. These shocks are likely to be radiative due to the high densities of the nova ejecta at early times coincident with the gamma-ray emission. Thermal X-rays radiated behind the shock are absorbed by neutral gas and reprocessed into optical emission, similar to Type IIn (interacting) supernovae. Gamma-rays are produced by collisions between relativistic protons with the nova ejecta (hadronic scenario) or Inverse Compton/bremsstrahlung emission from relativistic electrons (leptonic scenario), where in both scenarios the efficiency for converting relativistic particle energy into LAT gamma-rays is at most a few tens of per cent. The ratio of gamma-ray and optical luminosities, L_gam/L_opt, thus sets a lower limit on the fraction of the shock power used to accelerate relativistic particles, e_nth. The measured values of L_gam/L_opt for two classical novae, V1324 Sco and V339 Del, constrains e_nth > 1e-2 and > 1e-3, respectively. Inverse Compton models for the gamma-ray emission are disfavored given the low electron acceleration efficiency, e_nth ~ 1e-4-1e-3, inferred from observations of Galactic cosmic rays and particle-in-cell (PIC) numerical simulations. A fraction > 100(0.01/e_nth) and > 10(0.01/e_nth) per cent of the optical luminosity is powered by shocks in V1324 Sco and V339 Del, respectively. Such high fractions challenge standard models that instead attribute all nova optical emission to the direct outwards transport of thermal energy released near the white dwarf surface.


title: A photometric selection of White Dwarf candidates in SDSS DR10

abstract: We present a method which uses cuts in colour-colour and reduced proper motion-colour space to select white dwarfs without the recourse to spectroscopy while allowing an adjustable compromise between completeness and efficiency. Rather than just producing a list of white dwarf candidates, our method calculates a probability of being a white dwarf (Pwd) for any object with available multi band photometry and proper motion. We applied this method to all objects in the SDSS DR10 photometric footprint and to a few selected sources in DR7 which did not have reliable photometry in DR9 or DR10. This application results in a sample of 61969 DR10 and 3799 DR7 photometric sources with calculated Pwd from which it is possible to select a sample of ~23000 high-fidelity white dwarf candidates with Teff >~ 7000 K and <= 19. This sample contains over 14000 high confidence white dwarfs candidates which have not yet received spectroscopic follow-up. These numbers show that, to date, the spectroscopic coverage of white dwarfs in the SDSS photometric footprint is, on average, only ~40% complete. While we describe here in detail the application of our selection to the SDSS catalogue, the same method could easily be applied to other multi colour, large area surveys. We also publish a list of 8701 bright (<= 19) WDs with SDSS spectroscopy, of which 1781 are new identifications in DR9/10.


title: Naturalness in the Dark at the LHC

abstract: We revisit the Twin Higgs scenario as a "dark" solution to the little hierarchy problem, identify the structure of a minimal model and its viable parameter space, and analyze its collider implications. In this model, dark naturalness generally leads to Hidden Valley phenomenology. The twin particles, including the top partner, are all Standard-Model-neutral, but naturalness favors the existence of twin strong interactions -- an asymptotically-free force that confines not far above the Standard Model QCD scale -- and a Higgs portal interaction. We show that, taken together, these typically give rise to exotic decays of the Higgs to twin hadrons. Across a substantial portion of the parameter space, certain twin hadrons have visible and often displaced decays, providing a potentially striking LHC signature. We briefly discuss appropriate experimental search strategies.


title: Uncertainties in The Interstellar Extinction Curve and the Cepheid Distance to M101

abstract: I revisit the Cepheid-distance determination to the nearby spiral galaxy M101 (Pinwheel Galaxy) of Shappee & Stanek (2011), in light of several recent investigations questioning the shape of the interstellar extinction curve at $\lambda \approx 8,000$ \AA (i.e. I-band). I find that the relatively steep extinction ratio $A_{I}/E(V-I)=1.1450$ (Fitzpatrick & Massa 2007) is slightly favoured relative to $A_{I}/E(V-I)=1.2899$ (Fitzpatrick 1999) and significantly favoured relative the historically canonical value of $A_{I}/E(V-I)=1.4695$ (Cardelli et al. 1989). The steeper extinction curves, with lower values of $A_{I}/E(V-I)$, yield fits with reduced scatter, metallicity-dependences to the dereddened Cepheid luminosities that are closer to values inferred in the local group, and that are less sensitive to the choice of reddening cut imposed in the sample selection. The increase in distance modulus to M101 when using the preferred extinction curve is ${\Delta}{\mu} \sim 0.06$ mag, resulting in an estimate of the distance modulus to M101 relative to the LMC of $ {\Delta}\mu_{\rm{LMC}} \approx 10.72 \pm 0.03$ (stat). The best-fit metallicity-dependence is $dM_{I}/d\rm{[O/H]} \approx (-0.38 \pm 0.14$ (stat)) mag dex$^{-1}$.


title: Overturning the Case for Gravitational Powering in the Prototypical Cooling Lyman-alpha Nebula

abstract: The Nilsson et al. (2006) Lyman-alpha nebula has often been cited as the most plausible example of a Lyman-alpha nebula powered by gravitational cooling. In this paper, we bring together new data from the Hubble Space Telescope and the Herschel Space Observatory as well as comparisons to recent theoretical simulations in order to revisit the questions of the local environment and most likely power source for the Lyman-alpha nebula. In contrast to previous results, we find that this Lyman-alpha nebula is associated with 6 nearby galaxies and an obscured AGN that is offset by $\sim$4"$\approx$30 kpc from the Lyman-alpha peak. The local region is overdense relative to the field, by a factor of $\sim$10, and at low surface brightness levels the Lyman-alpha emission appears to encircle the position of the obscured AGN, highly suggestive of a physical association. At the same time, we confirm that there is no compact continuum source located within $\sim$2-3"$\approx$15-23 kpc of the Lyman-alpha peak. Since the latest cold accretion simulations predict that the brightest Lyman-alpha emission will be coincident with a central growing galaxy, we conclude that this is actually a strong argument against, rather than for, the idea that the nebula is gravitationally-powered. While we may be seeing gas within cosmic filaments, this gas is primarily being lit up, not by gravitational energy, but due to illumination from a nearby buried AGN.


title: An Emerging Wolf-Rayet Massive Star Cluster in NGC 4449

abstract: We present a panchromatic investigation of the partially-embedded, emerging massive cluster Source 26 (= S26) in NGC 4449 with optical spectra obtained at Apache Point Observatory and archival Hubble, Spitzer, and Herschel Space Telescope images. First identified as a radio continuum source with a thermal component due to ionized material, the massive cluster S26 also exhibits optical Wolf-Rayet (WR) emission lines that reveal a large evolved massive star population. We find that S26 is host to $\sim$240 massive stars, of which $\sim$18 are Wolf-Rayet stars; the relative populations are roughly consistent with other observed massive star forming clusters and galaxies. We construct SEDs over two spatial scales (roughly 100 pc and 300 pc) that clearly exhibit warm dust and polycyclic aromatic hydrocarbon (PAH) emission. The best fit dust and grain models reveal that both the intensity of the exciting radiation and PAH grain destruction increase toward the cluster center. Given that the timescale of evacuation is important for the future dynamical evolution of the cluster, it is important to determine whether O- and WR stars can evacuate the material gradually before supernova do so on a much faster timescale. With a minimum age of $\approx$ 3 Myr, it is clear that S26 has not yet fully evacuated its natal material, which indicates that unevolved O-type stars alone do not provide sufficient feedback to remove the gas and dust. We hypothesize that the feedback of WR stars in this cluster may be necessary for clearing the material from the gravitational potential of the cluster. We find S26 is similar to Emission Line Clusters observed in the Antennae Galaxies and may be considered a younger analog to 30 Doradus in the LMC.


title: How to form planetesimals from mm-sized chondrules and chondrule aggregates

abstract: The size distribution of asteroids and Kuiper belt objects in the solar system is difficult to reconcile with a bottom-up formation scenario due to the observed scarcity of objects smaller than $\sim$100 km in size. Instead, planetesimals appear to form top-down, with large $100-1000$ km bodies forming from the rapid gravitational collapse of dense clumps of small solid particles. In this paper we investigate the conditions under which solid particles can form dense clumps in a protoplanetary disk. We use a hydrodynamic code to model the interaction between solid particles and the gas inside a shearing box inside the disk, considering particle sizes from sub-millimeter-sized chondrules to meter-sized rocks. We find that particles down to millimeter sizes can form dense particle clouds through the run-away convergence of radial drift known as the streaming instability. We make a map of the range of conditions (strength of turbulence, particle mass-loading, disk mass, and distance to the star) which are prone to producing dense particle clumps. Finally, we estimate the distribution of collision speeds between mm-sized particles. We calculate the rate of sticking collisions and obtain a robust upper limit on the particle growth timescale of $\sim$$10^5$ years. This means that mm-sized chondrule aggregates can grow on a timescale much smaller than the disk accretion timescale ($\sim$$10^6 - 10^7$ years). Our results suggest a pathway from the mm-sized grains found in primitive meteorites to fully formed asteroids. We speculate that asteroids may form from a positive feedback loop in which coagualation leads to particle clumping driven by the streaming instability. This clumping, in turn reduces collision speeds and enhances coagulation.} Future simulations should model coagulation and the streaming instability together to explore this feedback loop further.


title: Virasoro Conformal Blocks and Thermality from Classical Background Fields

abstract: We show that in 2d CFTs at large central charge, the coupling of the stress tensor to heavy operators can be re-absorbed by placing the CFT in a non-trivial background metric. This leads to a more precise computation of the Virasoro conformal blocks between heavy and light operators, which are shown to be equivalent to global conformal blocks evaluated in the new background. We also generalize to the case where the operators carry U(1) charges. The refined Virasoro blocks can be used as the seed for a new Virasoro block recursion relation expanded in the heavy-light limit. We comment on the implications of our results for the universality of black hole thermality in $AdS_3$, or equivalently, the eigenstate thermalization hypothesis for $CFT_2$ at large central charge.


title: Composition of the Fermi-LAT isotropic gamma-ray background intensity: Emission from extragalactic point sources and dark matter annihilations

abstract: A new estimation of the isotropic diffuse gamma-ray background (IGRB) observed by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope (Fermi) has been presented for 50 months of data, in the energy range 100 MeV-820 GeV and for different modelings of the Galactic foreground. We attempt here the interpretation of the Fermi-LAT IGRB data in terms of the gamma-ray unresolved emission from different extragalactic populations. We find very good fits to the experimental IGRB, obtained with theoretical predictions for the emission from active galactic nuclei and star-forming galaxies. In addition, we probe a possible emission coming from the annihilation of weakly interacting dark matter (DM) particles in the halo of our Galaxy. We set stringent limits on its annihilation cross section into gamma rays, which are about the thermal relic value for a wide range of DM masses. We also identify regions in the DM mass and annihilation cross section parameter space which can significantly improve the fit to the IGRB data. Our analysis is conducted within the different IGRB data sets obtained from different models for the Galactic emission, which is shown to add a significant ambiguity on the IGRB interpretation.


title: Observations of the Icy Universe

abstract: Freeze-out of the gas phase elements onto cold grains in dense interstellar and circumstellar media builds up ice mantles consisting of molecules that are mostly formed in situ (H2O, NH3, CO2, CO, CH3OH, and more). This review summarizes the detected infrared spectroscopic ice features and compares the abundances across Galactic, extragalactic, and solar system environments. A tremendous amount of information is contained in the ice band profiles. Laboratory experiments play a critical role in the analysis of the observations. Strong evidence is found for distinct ice formation stages, separated by CO freeze out at high densities. The ice bands have proven to be excellent probes of the thermal history of their environment. The evidence for the long-held idea that processing of ices by energetic photons and cosmic rays produces complex molecules is weak. Recent state of the art observations show promise for much progress in this area with planned infrared facilities.


title: Electrically-Charged Lifshitz Spacetimes, and Hyperscaling Violations

abstract: Electrically-charged Lifshitz spacetimes are hard to come by. In this paper, we construct a class of such solutions in five dimensional Einstein gravity coupled to Maxwell and $SU(2)$ Yang-Mills fields. The solutions are electrically-charged under the Maxwell field, whose equation is sourced by the Yang-Mills instanton(-like) configuration living in the hyperbolic four-space of the Lifshitz spacetime. We then introduce a dilaton and construct charged and colored Lifshitz spacetimes with hyperscaling violations. We obtain a class of exact Lifshitz black holes. We also perform similar constructions in four dimensions.


title: Maximum nonlocality and minimum uncertainty using magic states

abstract: We prove that magic states from the Clifford hierarchy give optimal solutions for tasks involving nonlocality and entropic uncertainty with respect to Pauli measurements. For both the nonlocality and uncertainty tasks, stabilizer states are the worst possible pure states so our solutions have an operational interpretation as being highly non-stabilizer. The optimal strategy for a qudit version of the Clauser-Horne-Shimony-Holt (CHSH) game in prime dimensions is achieved by measuring maximally entangled states that are isomorphic to single-qudit magic states. These magic states have an appealingly simple form and our proof shows that they are "balanced" with respect to all but one of the mutually unbiased stabilizer bases. Of all equatorial qudit states, magic states minimize the average entropic uncertainties for collision entropy and also, for small prime dimensions, min-entropy -- a fact that may have implications for cryptography.


title: Exotic topological density waves in cold atomic Rydberg fermions

abstract: Versatile controllability of interactions in ultracold atomic and molecular gases has now reached an unprecedented era where quantum correlations and unconventional many-body phases can be studied with no corresponding analogs in solid state systems. Recent experiments in Rydberg atomic gases have achieved exquisite control over non-local interactions, allowing novel quantum phases unreachable with the usual local interactions in atomic systems. Here, we study Rydberg dressed atomic fermions in a three dimensional optical lattice predicting the existence of hitherto unheard-of exotic mixed topological density wave phases. We show that varying spatial range of the non-local interaction leads to a rich phase diagram containing various bond density waves, with unexpected spontaneous time-reversal symmetry breaking. Quasiparticles in these chiral phases experience emergent gauge fields and form three dimensional quantum Hall and Weyl semimetal states. Remarkably, certain density waves even exhibit mixed topologies beyond the existing topological classification. Experimental signatures of density waves and their topological properties are predicted in time-of-flight measurements.


title: Interstellar Medium Mitigation Techniques in Pulsar Timing Arrays

abstract: Pulsar Timing Arrays use a set of millisecond pulsars in an attempt to directly detect nanohertz gravitational waves. For this purpose, high precision timing of the pulsars is essential and ultimately a precision of the order of ~100 ns is required. Propagation effects in the interstellar medium cause the radio emission from a pulsar to be dispersed and scattered, introducing time variable delays of the pulses on their way to Earth. If these delays are not properly corrected for, they may cause significant errors in the timing analysis of a pulsar. These proceedings will review the effects of the interstellar medium on pulse arrival times and present some of the techniques used to mitigate the associated time delays from the pulsar signal. Correcting for these delays is essential to providing a higher timing precision and hence to increasing the array's sensitivity to gravitational waves.


title: Confronting uncertainties in stellar physics: calibrating convective overshooting with eclipsing binaries

abstract: As part of a larger program aimed at better quantifying the uncertainties in stellar computations, we attempt to calibrate the extent of convective overshooting in low to intermediate mass stars by means of eclipsing binary systems. We model 12 such systems, with component masses between 1.3 and 6.2 solar masses, using the detailed binary stellar evolution code STARS, producing grids of models in both metallicity and overshooting parameter. From these, we determine the best fit parameters for each of our systems. For three systems, none of our models produce a satisfactory fit. For the remaining systems, no single value for the convective overshooting parameter fits all the systems, but most of our systems can be well described with an overshooting parameter between 0.09 and 0.15, corresponding to an extension of the mixed region above the core of about 0.1-0.3 pressure scale heights. Of the nine systems where we are able to obtain a good fit, seven can be reasonably well fit with a single parameter of 0.15. We find no evidence for a trend of the extent of overshooting with either mass or metallicity, though the data set is of limited size. We repeat our calculations with a second evolution code, MESA, and we find general agreement between the two codes. For the extension of the mixed region above the convective core required by the MESA models is about 0.15-0.4 pressure scale heights. For the system EI Cep, we find that MESA gives an overshooting region that is larger than the STARS one by about 0.1 pressure scale heights for the primary, while for the secondary the difference is only 0.05 pressure scale heights.


title: Passive galaxies as tracers of cluster environments at z~2

abstract: Even 10 billion years ago, the cores of the first galaxy clusters are often found to host a characteristic population of massive galaxies with already suppressed star formation. Here we search for distant cluster candidates at z~2 using massive passive galaxies as tracers. With a sample of ~40 spectroscopically confirmed passive galaxies at 1.3<z<2.1, we tune photometric redshifts of several thousands passive sources in the full 2 sq.deg. COSMOS field. This allows us to map their density in redshift slices, probing the large scale structure in the COSMOS field as traced by passive sources. We report here on the three strongest passive galaxy overdensities that we identify in the redshift range 1.5<z<2.5. While the actual nature of these concentrations is still to be confirmed, we discuss their identification procedure, and the arguments supporting them as candidate galaxy clusters (likely mid-10^13 M_sun range). Although this search approach is likely biased towards more evolved structures, it has the potential to select still rare, cluster-like environments close to their epoch of first appearance, enabling new investigations of the evolution of galaxies in the context of structure growth.


title: SU(2)-symmetry in a realistic spin-fermion model for cuprate superconductors

abstract: We consider the Pseudo-Gap (PG) state of high-$T_c$ superconductors in form of a composite order parameter fluctuating between 2$p_F$-charge ordering and superconducting (SC) pairing. In the limit of linear dispersion and at the hotspots, both order parameters are related by a SU(2) symmetry and the eight hotspot model of Efetov $\textit{et al.}$ [Nat. Phys. $\bf{9}$, 442 (2013)] is recovered. In the general case however, curvature terms of the dispersion will break this symmetry and the degeneracy between both states is lifted. Taking the full momentum dependence of the order parameter into account, we measure the strength of this SU(2) symmetry breaking over the full Brillouin zone. For realistic dispersion relations including curvature we find generically that the SU(2) symmetry breaking is small and robust to the fermiology and that the symmetric situation is restored in the large paramagnon mass and coupling limit. Comparing the level splitting for different materials we propose a scenario that could account for the competition between the PG and the SC states in the phase diagram of high-$T_c$ superconductors.


title: Measurement of the ratio of inclusive cross sections $\sigma (p\bar{p} \rightarrow Z+2~b \text{jets}) / \sigma (p\bar{p} \rightarrow Z+ \text{2 jets})$ in $p\bar{p}$ collisions at $\sqrt s=1.96$ TeV

abstract: We measure the ratio of cross sections, $\sigma(p\bar{p}\rightarrow Z+2~b~\text{jets})$/$\sigma(p\bar{p}\rightarrow Z+\text{2 jets})$, for associated production of a $Z$ boson with at least two jets with transverse momentum $p_T^{\rm jet} > 20$ GeV and pseudorapidity $|\eta^{\rm jet}| < 2.5$. This measurement uses data corresponding to an integrated luminosity of 9.7 fb$^{-1}$ collected by the \DO experiment in Run II of Fermilab's Tevatron \ppbar Collider at a center-of-mass energy of 1.96 TeV. The measured integrated ratio of $0.0236\pm0.0032\left(\mbox{stat}\right)\pm0.0035\left(\mbox{syst}\right)$ is in agreement with predictions from next-to-leading-order perturbative QCD and the Monte Carlo event generators {\sc pythia} and {\sc alpgen}.


title: Prospects for a Narrow Line MOT in YO

abstract: In addition to being suitable for laser cooling and trapping in a magneto-optical trap (MOT) using a relatively broad ($\sim$5 MHz) transition, the molecule YO possesses a narrow-line transition. This forbidden transition between the $X^{2}\Sigma$ and $A'^{2}\Delta_{3/2}$ states has linewidth $\sim$2$\pi\times$160 kHz. After cooling in a MOT on the $X^{2}\Sigma$ to $A^{2}\Pi_{1/2}$ (orange) transition, the narrow (red) transition can be used to further cool the sample, requiring only minimal additions to the first stage system. The narrow line cooling stage will bring the temperature from $\sim$1 mK to $\sim$10 $\mu$K, significantly advancing the frontier on direct cooling achievable for molecules.


title: Turing jumps through provability

abstract: Fixing some computably enumerable theory $T$, the Friedman-Goldfarb-Harrington (FGH) theorem says that over elementary arithmetic, each $\Sigma_1$ formula is equivalent to some formula of the form $\Box_T \varphi$ provided that $T$ is consistent. In this paper we give various generalizations of the FGH theorem. In particular, for $n>1$ we relate $\Sigma_{n}$ formulas to provability statements $[n]_T^{\sf True}\varphi$ which are a formalization of "provable in $T$ together with all true $\Sigma_{n+1}$ sentences". As a corollary we conclude that each $[n]_T^{\sf True}$ is $\Sigma_{n+1}$-complete. This observation yields us to consider a recursively defined hierarchy of provability predicates $[n+1]^\Box_T$ which look a lot like $[n+1]_T^{\sf True}$ except that where $[n+1]_T^{\sf True}$ calls upon the oracle of all true $\Sigma_{n+2}$ sentences, the $[n+1]^\Box_T$ recursively calls upon the oracle of all true sentences of the form $\langle n \rangle_T^\Box\phi$. As such we obtain a `syntax-light' characterization of $\Sigma_{n+1}$ definability whence of Turing jumps which is readily extended beyond the finite. Moreover, we observe that the corresponding provability predicates $[n+1]_T^\Box$ are well behaved in that together they provide a sound interpretation of the polymodal provability logic ${\sf GLP}_\omega$.


title: Finitely Balanced Sequences and Plasticity of 1-Dimensional Tilings

abstract: We relate a balancing property of letters for bi-infinite sequences to the invariance of the resulting 1-dimensional tiling dynamics under changes in the lengths of the tiles. If the language of the sequence space is finitely balanced, then all length changes in the corresponding tiling space result in topological conjugacies, up to an overall rescaling.


title: The nonlinear steepest descent method for Riemann-Hilbert problems of low regularity

abstract: We prove a nonlinear steepest descent theorem for Riemann-Hilbert problems with Carleson jump contours and jump matrices of low regularity and slow decay. We illustrate the theorem by deriving the long-time asymptotics for the mKdV equation in the similarity sector for initial data with limited decay and regularity.


title: Systematic Uncertainties in the Spectroscopic Measurements of Neutron-Star Masses and Radii from Thermonuclear X-ray Bursts. III. Absolute Flux Calibration

abstract: Many techniques for measuring neutron star radii rely on absolute flux measurements in the X-rays. As a result, one of the fundamental uncertainties in these spectroscopic measurements arises from the absolute flux calibrations of the detectors being used. Using the stable X-ray burster, GS 1826-238, and its simultaneous observations by Chandra HETG/ACIS-S and RXTE/PCA as well as by XMM-Newton EPIC-pn and RXTE/PCA, we quantify the degree of uncertainty in the flux calibration by assessing the differences between the measured fluxes during bursts. We find that the RXTE/PCA and the Chandra gratings measurements agree with each other within their formal uncertainties, increasing our confidence in these flux measurements. In contrast, XMM-Newton EPIC-pn measures 14.0$\pm$0.3 % less flux than the RXTE/PCA. This is consistent with the previously reported discrepancy with the flux measurements of EPIC-pn, compared to EPIC-MOS1, MOS2 and ACIS-S detectors. We also show that any intrinsic time dependent systematic uncertainty that may exist in the calibration of the satellites has already been implicity taken into account in the neutron star radius measurements.


title: Renormalisation of Newton's constant

abstract: The problem of obtaining a gauge independent beta function for Newton's constant is addressed. By a specific parameterisation of metric fluctuations a gauge independent functional integral is constructed for the semiclassical theory around an arbitrary Einstein space. The effective action then has the property that only physical polarisations of the graviton contribute, while all other modes cancel with the functional measure. We are then able to compute a gauge independent beta function for Newton's constant in $d$-dimensions to one-loop order. No Landau pole is present provided $N_g < 18$, where $N_g = d(d-3)/2$ is the number of polarisations of the graviton. While adding a large number of matter fields can change this picture, the absence of a pole persists for the particle content of the standard model in four spacetime dimensions.


title: Commissioning of the Inner Tracker of the KLOE-2 experiment

abstract: The KLOE-2 experiment is undergoing commissioning at the DA$\Phi$NE $e^{+}e^{-}$ collider of the Frascati National Laboratory of the INFN, after the integration of new detectors in the former KLOE apparatus. The Inner Tracker, a very light detector (material budget $<2\%X_{0}$), is one of the new subdetectors and it is composed of 4 cylindrical triple-GEM with a stereo X-V strips/pads readout. The commissioning phase of the Inner Tracker consists in the characterization of the detector response and in its performance evaluation. The method used to evaluate its detection efficiency is reported, together with some preliminary results.


title: Investigation of Kelvin-like solid foams for potential engineering applications: An attractive set of geometrical and thermo-hydraulic properties

abstract: Open cell foams have diverse industrial applications e.g. heat exchangers, structured reactors, filtration due to their unique properties such as high porosity and high specific surface area. In order to theoretically determine the geometric specific surface area and relationships between geometrical parameters of isotropic open cell foams, a generalized mathematical correlation was developed. For this purpose the tetrakaidecahedron geometry was used and different shapes of strut cross-sections of foam structures were taken explicitly into account. The derived correlation to predict geometrical properties can be easily extended to different strut shapes. 3-D numerical simulations at pore scale were performed to study the pressure drop characteristics and effective thermal conductivity. Fluid flow through open cell foam was performed in three different regimes: Darcy regime, transition regime and inertia regime. Importance of geometrical properties on fluid flow characteristics and their inclusion in the proposed correlations for predicting pressure drop is discussed. Can Ergun parameters have constant numerical values or not is also extensively discussed. Three different correlations were derived to predict the effective thermal conductivity for both, isotropic and anisotropic open cell foams. Geometrical parameters of foam matrix were introduced in the correlations to predict effective thermal conductivity.


title: Mirror Symmetry and the Classification of Orbifold del Pezzo Surfaces

abstract: We state a number of conjectures that together allow one to classify a broad class of del Pezzo surfaces with cyclic quotient singularities using mirror symmetry. We prove our conjectures in the simplest cases. The conjectures relate mutation-equivalence classes of Fano polygons with Q-Gorenstein deformation classes of del Pezzo surfaces.


title: Minimality and mutation-equivalence of polygons

abstract: We introduce a concept of minimality for Fano polygons. We show that, up to mutation, there are only finitely many Fano polygons with given singularity content, and give an algorithm to determine the mutation-equivalence classes of such polygons. This is a key step in a program to classify orbifold del Pezzo surfaces using mirror symmetry. As an application, we classify all Fano polygons such that the corresponding toric surface is qG-deformation-equivalent to either (i) a smooth surface; or (ii) a surface with only singularities of type 1/3(1,1).


title: Higher-degree Dirac Currents of Twistor and Killing Spinors in Supergravity Theories

abstract: We show that higher degree Dirac currents of twistor and Killing spinors correspond to the hidden symmetries of the background spacetime which are generalizations of conformal Killing and Killing vector fields respectively. They are the generalizations of 1-form Dirac currents to higher degrees which are used in constructing the bosonic supercharges in supergravity theories. In the case of Killing spinors, we find that the equations satisfied by the higher degree Dirac currents are related to Maxwell-like and Duffin-Kemmer-Petiau equations. Correspondence between the Dirac currents and harmonic forms for parallel and pure spinor cases is determined. We also analyze the supergravity twistor and Killing spinor cases in 10 and 11-dimensional supergravity theories and find that although different inner product classes induce different involutions on spinors, the higher degree Dirac currents still correspond to the hidden symmetries of the spacetime.


title: The free boundary of variational inequalities with gradient constraints

abstract: In this paper we prove that the free boundary of some variational inequalities with gradient constraints is as regular as the tangent bundle of the boundary of the domain. To this end, we study a generalized notion of ridge of a domain in the plane, which is the set of singularities of the distance function in the p-norm to the boundary of the domain.


title: Semantics for Locking Specifications

abstract: To prevent concurrency errors, programmers need to obey a locking discipline. Annotations that specify that discipline, such as Java's @GuardedBy, are already widely used. Unfortunately, their semantics is expressed informally and is consequently ambiguous. This article highlights such ambiguities and formalizes the semantics of @GuardedBy in two alternative ways, building on an operational semantics for a small concurrent fragment of a Java-like language. It also identifies when such annotations are actual guarantees against data races. Our work aids in understanding the annotations and supports the development of sound formal tools that verify or infer such annotations.


title: The regularity of some vector-valued variational inequalities with gradient constraints

abstract: We prove the optimal regularity for some class of vector-valued variational inequalities with gradient constraints. We also give a new proof for the optimal regularity of some scalar variational inequalities with gradient constraints. In addition, we prove that some class of variational inequalities with gradient constraints are equivalent to an obstacle problem, both in the scalar and vector-valued case.


title: An unconditionally stable discontinuous Galerkin method for the elastic Helmholtz equations with large frequency

abstract: In this paper we propose and analyze an interior penalty discontinuous Galerkin (IP-DG) method using piecewise linear polynomials for the elastic Helmholtz equations with the first order absorbing boundary condition. It is proved that the sesquilinear form for the problem satisfies a generalized weak coercivity property, which immediately infers a stability estimate for the solution of the differential problem in all frequency regimes. It is also proved that the proposed IP-DG method is unconditionally stable with respect to both frequency $\omega$ and mesh size $h$. Sub-optimal order (with respect to $h$) error estimates in the broken $H^1$-norm and in the $L^2$-norm are obtained in all mesh regimes. These estimate improve to optimal order when the mesh size $h$ is restricted to the pre-asymptotic regime (i.e., $\omega^\beta h =O(1)$ for some $1\leq \beta<2$). The novelties of the proposed IP-DG method include: first, the method penalizes not only the jumps of the function values across the element edges but also the jumps of the normal derivatives; second, the penalty parameters are taken as complex numbers with positive imaginary parts. In order to establish the desired unconditional stability estimate for the numerical solution, the main idea is to exploit a (simple) property of linear functions to overcome the main difficulty caused by non-Hermitian nature and strong indefiniteness of the Helmholtz-type problem. The error estimate is then derived using a nonstandard technique adapted from \cite{Feng_Wu09}. Numerical experiments are also presented to validate the theoretical results and to numerically examine the pollution effect (with respect to $\omega$) in the error bounds.


title: Energy dependence of acceptance-corrected dielectron excess mass spectrum at mid-rapidity in Au+Au collisions at $\sqrt{s_{NN}} = 19.6$ and 200 GeV

abstract: The acceptance-corrected dielectron excess mass spectra, where the known hadronic sources have been subtracted from the inclusive dielectron mass spectra, are reported for the first time at mid-rapidity $|y_{ee}|<1$ in minimum-bias Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 and 200 GeV. The excess mass spectra are consistently described by a model calculation with a broadened $\rho$ spectral function for $M_{ee}<1.1$ GeV/$c^{2}$. The integrated dielectron excess yield at $\sqrt{s_{NN}}$ = 19.6 GeV for $0.4<M_{ee}<0.75$ GeV/$c^2$, normalized to the charged particle multiplicity at mid-rapidity, has a value similar to that in In+In collisions at $\sqrt{s_{NN}}$ = 17.3 GeV. For $\sqrt{s_{NN}}$ = 200 GeV, the normalized excess yield in central collisions is higher than that at $\sqrt{s_{NN}}$ = 17.3 GeV and increases from peripheral to central collisions. These measurements indicate that the lifetime of the hot, dense medium created in central Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV is longer than those in peripheral collisions and at lower energies.


title: Hausdorff volume in non equiregular sub-Riemannian manifolds

abstract: In this paper we study the Hausdorff volume in a non equiregular sub-Riemannian manifold and we compare it with a smooth volume. We first give the Lebesgue decomposition of the Hausdorff volume. Then we study the regular part, show that it is not commensurable with the smooth volume, and give conditions under which it is a Radon measure. We finally give a complete characterization of the singular part. We illustrate our results and techniques on numerous examples and cases (e.g. to generic sub-Riemannian structures).


title: NANOGrav Constraints on Gravitational Wave Bursts with Memory

abstract: Among efforts to detect gravitational radiation, pulsar timing arrays are uniquely poised to detect "memory" signatures, permanent perturbations in spacetime from highly energetic astrophysical events such as mergers of supermassive black hole binaries. The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) observes dozens of the most stable millisecond pulsars using the Arecibo and Green Bank radio telescopes in an effort to study, among other things, gravitational wave memory. We herein present the results of a search for gravitational wave bursts with memory (BWMs) using the first five years of NANOGrav observations. We develop original methods for dramatically speeding up searches for BWM signals. In the directions of the sky where our sensitivity to BWMs is best, we would detect mergers of binaries with reduced masses of $10^9$ $M_\odot$ out to distances of 30 Mpc; such massive mergers in the Virgo cluster would be marginally detectable. We find no evidence for BWMs. However, with our non-detection, we set upper limits on the rate at which BWMs of various amplitudes could have occurred during the time spanned by our data--e.g., BWMs with amplitudes greater than $10^{-13}$ must occur at a rate less than 1.5 yr$^{-1}$.


title: Dynamics of airflow in a short inhalation

abstract: During a rapid inhalation, such as a sniff, the flow in the airways accelerates and decays quickly. The consequences for flow development and convective trans- port of an inhaled gas were investigated in a subject geometry extending from the nose to the bronchi. The progress of flow transition and the advance of an inhaled non-absorbed gas were determined using highly resolved simulations of a sniff 0.5 s long, 1 litre per second peak flow, 364 ml inhaled volume. In the nose, the distribution of airflow evolved through three phases: (i) an initial transient of about 50 ms, roughly the filling time for a nasal volume, (ii) quasi-equilibrium over the majority of the inhalation, and (iii) a terminating phase. Flow transition commenced in the supraglottic region within 20ms, resulting in large- amplitude fluctuations persisting throughout the inhalation; in the nose, fluctuations that arose nearer peak flow were of much reduced intensity and diminished in the flow decay phase. Measures of gas concentration showed non-uniform build-up and wash-out of the inhaled gas in the nose. At the carina, the form of the temporal concentration profile reflected both shear dispersion and airway filling defects owing to recirculation regions.


title: Most linear flows on $\mathbb{R}^d$ are Benford

abstract: A necessary and sufficient condition ("exponential nonresonance") is established for every signal obtained from a linear flow on $\mathbb{R}^d$ by means of a linear observable to either vanish identically or else exhibit a strong form of Benford's Law (logarithmic distribution of significant digits). The result extends and unifies all previously known (sufficient) conditions. Exponential nonresonance is shown to be typical for linear flows, both from a topological and a measure-theoretical point of view.


title: Universal distribution of magnetic anisotropy of impurities in ordered and disordered nano-grains

abstract: We examine the distribution of the magnetic anisotropy (MA) experienced by a magnetic impurity embedded in a metallic nano-grain. As an example of a generic magnetic impurity with partially filled $d$-shell, we study the case of $d^{1}$ impurities imbedded into ordered and disordered Au nano-grains, described in terms of a realistic band structure. Confinement of the electrons induces a magnetic anisotropy that is large, and can be characterized by 5 real parameters, coupling to the quadrupolar moments of the spin. In ordered (spherical) nano-grains, these parameters exhibit symmetrical structures and reflect the symmetry of the underlying lattice, while for disordered grains they are randomly distributed and, - for stronger disorder, - their distribution is found to be characterized by random matrix theory. As a result, the probability of having small magnetic anisotropies $K_L$ is suppressed below a characteristic scale $\Delta_E$, which we predict to scale with the number of atoms $N$ as $\Delta_E\sim 1/N^{3/2}$. This gives rise to anomalies in the specific heat and the susceptibility at temperatures $T\sim \Delta_E$ and produces distinct structures in the magnetic excitation spectrum of the clusters, that should be possible to detect experimentally.


title: The VMC Survey - XIV. First results on the look-back time star-formation rate tomography of the Small Magellanic Cloud

abstract: We analyse deep images from the VISTA survey of the Magellanic Clouds in the YJKs filters, covering 14 sqrdeg (10 tiles), split into 120 subregions, and comprising the main body and Wing of the Small Magellanic Cloud (SMC). We apply a colour--magnitude diagram reconstruction method that returns their best-fitting star formation rate SFR(t), age-metallicity relation (AMR), distance and mean reddening, together with 68% confidence intervals. The distance data can be approximated by a plane tilted in the East-West direction with a mean inclination of 39 deg, although deviations of up to 3 kpc suggest a distorted and warped disk. After assigning to every observed star a probability of belonging to a given age-metallicity interval, we build high-resolution population maps. These dramatically reveal the flocculent nature of the young star-forming regions and the nearly smooth features traced by older stellar generations. They document the formation of the SMC Wing at ages <0.2 Gyr and the peak of star formation in the SMC Bar at 40 Myr. We clearly detect periods of enhanced star formation at 1.5 Gyr and 5 Gyr. The former is possibly related to a new feature found in the AMR, which suggests ingestion of metal-poor gas at ages slightly larger than 1 Gyr. The latter constitutes a major period of stellar mass formation. We confirm that the SFR(t) was moderately low at even older ages.


title: Locally smeared operator product expansions in scalar field theory

abstract: We propose a new locally smeared operator product expansion to decompose nonlocal operators in terms of a basis of smeared operators. The smeared operator product expansion formally connects nonperturbative matrix elements determined numerically using lattice field theory to matrix elements of nonlocal operators in the continuum. These nonperturbative matrix elements do not suffer from power-divergent mixing on the lattice, which significantly complicates calculations of quantities such as the moments of parton distribution functions, provided the smearing scale is kept fixed in the continuum limit. The presence of this smearing scale complicates the connection to the Wilson coefficients of the standard operator product expansion and requires the construction of a suitable formalism. We demonstrate the feasibility of our approach with examples in real scalar field theory.


title: Exploiting Big Data in Logistics Risk Assessment via Bayesian Nonparametrics

abstract: In cargo logistics, a key performance measure is transport risk, defined as the deviation of the actual arrival time from the planned arrival time. Neither earliness nor tardiness is desirable for customer and freight forwarders. In this paper, we investigate ways to assess and forecast transport risks using a half-year of air cargo data, provided by a leading forwarder on 1336 routes served by 20 airlines. Interestingly, our preliminary data analysis shows a strong multimodal feature in the transport risks, driven by unobserved events, such as cargo missing flights. To accommodate this feature, we introduce a Bayesian nonparametric model -- the probit stick-breaking process (PSBP) mixture model -- for flexible estimation of the conditional (i.e., state-dependent) density function of transport risk. We demonstrate that using simpler methods, such as OLS linear regression, can lead to misleading inferences. Our model provides a tool for the forwarder to offer customized price and service quotes. It can also generate baseline airline performance to enable fair supplier evaluation. Furthermore, the method allows us to separate recurrent risks from disruption risks. This is important, because hedging strategies for these two kinds of risks are often drastically different.


title: Embedding degenerate graphs of small bandwidth

abstract: We develop a tool for embedding almost spanning degenerate graphs of small bandwidth. As an application, we extend the blow-up lemma to degenerate graphs of small bandwidth, the bandwidth theorem to degenerate graphs, and make progress on a conjecture of Burr and Erd\H{o}s on Ramsey number of degenerate graphs.


title: Nonlocal photon correlations and violation of Bell inequalities for spatially separated classical light fields

abstract: It is theoretically and experimentally shown that photons emitted by statistically independent incoherent classical light sources and measured in the far field in spatially separated modes may display spatial correlations akin to path-entanglement of photons produced by quantum sources. By measuring higher order photon-correlations at different locations, i.e., $m$ photons in one mode and one photon in another mode, we experimentally demonstrate for $m \geq 6$ a violation of Bell-type inequalities for spatial degrees of freedom. The spatial correlations among the photons can be understood from state projection where the detection of the first $m$ photons projects the sources onto a state which emits the subsequent photon in a strongly correlated manner. From this perspective the entanglement and violation of Bell's inequalities appears as a consequence of nonvanishing cross correlations between noncommuting quadrature phase components of the two spatially separated fields after $m$ photons have been recorded. In this way we show that classical systems may produce spatial field correlations violating local realistic theories.


title: Optimizing affinity-based binary hashing using auxiliary coordinates

abstract: In supervised binary hashing, one wants to learn a function that maps a high-dimensional feature vector to a vector of binary codes, for application to fast image retrieval. This typically results in a difficult optimization problem, nonconvex and nonsmooth, because of the discrete variables involved. Much work has simply relaxed the problem during training, solving a continuous optimization, and truncating the codes a posteriori. This gives reasonable results but is quite suboptimal. Recent work has tried to optimize the objective directly over the binary codes and achieved better results, but the hash function was still learned a posteriori, which remains suboptimal. We propose a general framework for learning hash functions using affinity-based loss functions that uses auxiliary coordinates. This closes the loop and optimizes jointly over the hash functions and the binary codes so that they gradually match each other. The resulting algorithm can be seen as a corrected, iterated version of the procedure of optimizing first over the codes and then learning the hash function. Compared to this, our optimization is guaranteed to obtain better hash functions while being not much slower, as demonstrated experimentally in various supervised datasets. In addition, our framework facilitates the design of optimization algorithms for arbitrary types of loss and hash functions.


title: Rapid Sequence Identification of Potential Pathogens Using Techniques from Sparse Linear Algebra

abstract: The decreasing costs and increasing speed and accuracy of DNA sample collection, preparation, and sequencing has rapidly produced an enormous volume of genetic data. However, fast and accurate analysis of the samples remains a bottleneck. Here we present D$^{4}$RAGenS, a genetic sequence identification algorithm that exhibits the Big Data handling and computational power of the Dynamic Distributed Dimensional Data Model (D4M). The method leverages linear algebra and statistical properties to increase computational performance while retaining accuracy by subsampling the data. Two run modes, Fast and Wise, yield speed and precision tradeoffs, with applications in biodefense and medical diagnostics. The D$^{4}$RAGenS analysis algorithm is tested over several datasets, including three utilized for the Defense Threat Reduction Agency (DTRA) metagenomic algorithm contest.


title: A speed and departure time optimization algorithm for the Pollution-Routing Problem

abstract: We propose a new speed and departure time optimization algorithm for the Pollution-Routing Problem (PRP), which runs in quadratic time and returns a certified optimal schedule. This algorithm is embedded into an iterated local search-based metaheuristic to achieve a combined speed, scheduling and routing optimization. The start of the working day is set as a decision variable for individual routes, thus enabling a better assignment of human resources to required demands. Some routes that were evaluated as unprofitable can now appear as viable candidates later in the day, leading to a larger search space and further opportunities of distance optimization via better service consolidation. Extensive computational experiments on available PRP benchmark instances demonstrate the good performance of the algorithms. The flexible departure times from the depot contribute to reduce the operational costs by 8.36% on the considered instances.


title: Incoherent excitation and switching of spin states in exciton-polariton condensates

abstract: We investigate, theoretically and numerically, the spin dynamics of a two-component exciton-polariton condensate created and sustained by non-resonant spin-polarized optical pumping of a semiconductor microcavity. Using the open-dissipative mean-field model, we show that the existence of well defined phase-locked steady states of the condensate may lead to efficient switching and control of spin (polarization) states with a non-resonant excitation. Spatially inhomogeneous pulsed excitations can cause symmetry breaking in the pseudo-spin structure of the condensate and lead to formation of non-trivial spin textures. Our model is universally applicable to two weakly coupled polariton condensates, and therefore can also describe the behaviour of condensate populations and phases in 'double-well' type potentials.


title: On the construction of bivariate linear exponential distribution with FGM family

abstract: In this paper we propose a Farlie-Gumbel-Morgenstern (FGM) family of bivariate linear exponential distributions generated from given marginal's. Therefore, properties of FGM are analogous to properties of bivariate distributions. We study some important statistical properties and results for the new distribution.


title: Forward Compton Scattering with weak neutral current: constraints from sum rules

abstract: We generalize forward real Compton amplitude to the case of the interference of the electromagnetic and weak neutral current, formulate a low-energy theorem, relate the new amplitudes to the interference structure functions and obtain a new set of sum rules. We address a possible new sum rule that relates the product of the axial charge and magnetic moment of the nucleon to the 0th moment of the structure function $g_5(\nu,0)$. We apply the GDH and the finite energy sum rule for constraining the dispersive $\gamma Z$-box correction to the proton's weak charge.


title: The Merger of Small and Large Black Holes

abstract: We present simulations of binary black holes mergers in which, after the common outer horizon has formed, the marginally outer trapped surfaces (MOTSs) corresponding to the individual black holes continue to approach and eventually penetrate each other. This has very interesting consequences according to recent results in the theory of MOTSs. Uniqueness and stability theorems imply that two MOTSs which touch with a common outer normal must be identical. This suggests a possible dramatic consequence of the collision between a small and large black hole. If the penetration were to continue to completion then the two MOTSs would have to coalesce, by some combination of the small one growing and the big one shrinking. Here we explore the relationship between theory and numerical simulations, in which a small black hole has halfway penetrated a large one.


title: Experimental Constraints on {\gamma}-ray Pulsar Gap Models and the Pulsar GeV to Pulsar Wind Nebula TeV Connection

abstract: The pulsar emission mechanism in the gamma-ray energy band is poorly understood. Currently, there are several models under discussion in the pulsar community. These models can be constrained by studying the collective properties of a sample of pulsars, which became possible with the large sample of gamma-ray pulsars discovered by the Fermi Large Area Telescope (Fermi-LAT). In this paper we develop a new experimental multi-wavelength technique to determine the beaming factor $\left( f_\Omega \right)$ dependance on spin-down luminosity of a set of GeV pulsars. This technique requires three input parameters: pulsar spin-down luminosity, pulsar phase-averaged GeV flux and TeV or X-ray flux from the associated Pulsar Wind Nebula (PWN). The analysis presented in this paper uses the PWN TeV flux measurements to study the correlation between $f_\Omega$ and $\dot{E}$. The measured correlation has some features that favor the Outer Gap model over the Polar Cap, Slot Gap and One Pole Caustic models for pulsar emission in the energy range of 0.1 to 100 GeV, but one must keep in mind that these simulated models failed to explain many of the most important pulsar population characteristics. A tight correlation between the pulsar GeV emission and PWN TeV emission was also observed, which suggests the possibility of a linear relationship between the two emission mechanisms. In this paper we also discuss a possible mechanism to explain this correlation.


title: The Asteroseismology of ZZ Ceti star GD1212

abstract: The ZZ Ceti star GD 1212 was detected to have 19 independent modes from the two-wheel-controlled Kepler Spacecraft in 2014. By asymptotic analysis, we identify most of pulsation modes. We find out two set of complete triplets, and four sets of doublet which are interpreted as rotation modes with $l=1$. For the other five modes, the four modes $f_{13}$, $f_{15}$, $f_{16}$ and $f_{4}$ are identified as ones with $l=2$; and the mode $f_{7}$ is identified to be the one with $l=1$. Meanwhile we derive a mean rotation period of $6.65\pm0.21$ h for GD 1212 according to the rotation splitting. Using the method of matching the observed periods to theoretical ones, we obtain the best-fitting model with the four parameters as $M_{\rm{*}}/M_{\rm{\odot}} = 0.775$, $T_{\rm{eff}} = 11400$ K, $\log (M_{\rm{H}}/M_{\rm{*}}) = -5.0$, $\log (M_{\rm{He}}/M_{\rm{*}})=-2.5$ for GD 1212. We find that due to the gradient of C/O abundance in the interior of white dwarf, some modes can not propagate to the stellar interior, which leads to the period spacing of the adjacent modes to become large. This feature is just proven by the observational data from GD 1212. All of these imply that GD 1212 may be evolved from an intermediate mass star.


title: Reconstruction of a fully anisotropic elasticity tensor from knowledge of displacement fields

abstract: We present explicit reconstruction algorithms for fully anisotropic unknown elasticity tensors from knowledge of a finite number of internal displacement fields, with applications to transient elastography. Under certain rank-maximality assumptions satified by the strain fields, explicit algebraic reconstruction formulas are provided. A discussion ensues on how to fulfill these assumptions, describing the range of validity of the approach. We also show how the general method can be applied to more specific cases such as the transversely isotropic one.


title: Towards Quantized Number Theory: Spectral Operators and an Asymmetric Criterion for the Riemann Hypothesis

abstract: This research expository article contains a survey of earlier work (in \S2--\S4) but also contains a main new result (in \S5), which we first describe. Given $c \geq 0$, the spectral operator $\mathfrak{a} = \mathfrak{a}_c$ can be thought of intuitively as the operator which sends the geometry onto the spectrum of a fractal string of dimension not exceeding $c$. Rigorously, it turns out to coincide with a suitable quantization of the Riemann zeta function $\zeta = \zeta(s)$: $\mathfrak{a} = \zeta (\partial)$, where $\partial = \partial_c$ is the infinitesimal shift of the real line acting on the weighted Hilbert space $L^2 (\mathbb{R}, e^{-2ct} dt)$. In this paper, we establish a new asymmetric criterion for the Riemann hypothesis, expressed in terms of the invertibility of the spectral operator for all values of the dimension parameter $c \in (0, 1/2)$ (i.e., for all $c$ in the left half of the critical interval $(0,1)$). This corresponds (conditionally) to a mathematical (and perhaps also, physical) "phase transition" occurring in the midfractal case when $c= 1/2$. Both the universality and the non-universality of $\zeta = \zeta (s)$ in the right (resp., left) critical strip $\{1/2 < \text{Re}(s) < 1 \}$ (resp., $\{0 < \text{Re}(s) < 1/2 \}$) play a key role in this context. These new results are presented in \S5. In \S2, we briefly discuss earlier joint work on the complex dimensions of fractal strings, while in \S3 and \S4, we survey earlier related work of the author with H. Maier and with H. Herichi, respectively, in which were established symmetric criteria for the Riemann hypothesis, expressed respectively in terms of a family of natural inverse spectral problems for fractal strings of Minkowski dimension $D \in (0,1),$ with $D \neq 1/2$, and of the quasi-invertibility of the family of spectral operators $\mathfrak{a}_c$ (with $c \in (0,1), c \neq 1/2$).


title: The extension class and KMS states for Cuntz--Pimsner algebras of some bi-Hilbertian bimodules

abstract: For bi-Hilbertian $A$-bimodules, in the sense of Kajiwara--Pinzari--Watatani, we construct a Kasparov module representing the extension class defining the Cuntz--Pimsner algebra. The construction utilises a singular expectation which is defined using the $C^*$-module version of the Jones index for bi-Hilbertian bimodules. The Jones index data also determines a novel quasi-free dynamics and KMS states on these Cuntz--Pimsner algebras.


title: Particle Trapping and Streaming Instability in Vortices

abstract: We analyse the concentration of solid particles in vortices created and sustained by radial buoyancy in protoplanetary disks, i.e. baroclinic vortex growth. Besides the gas drag acting on particles we also allow for back-reaction from dust onto the gas. This becomes important when the local dust-to-gas ratio approaches unity. In our 2D, local, shearing sheet simulations we see high concentrations of grains inside the vortices for a broad range of Stokes numbers, ${\rm St}$. An initial dust-to-gas ratio of 1:100 can easily be reversed to 100:1 for ${\rm St}=1$. The increased dust-to-gas ratio triggers the streaming instability, thus counter-intuitively limiting the maximal achievable overdensities. We find that particle trapping inside vortices opens the possibility for gravity-assisted planetesimal formation even for small particles ($\rm{St}=0.01$) and low initial dust-to-gas ratios (1:$10^4$)


title: Automated Reidemeister Moves: A Numerical Approach to the Unknotting Problem

abstract: In mathematics, a knot is a single strand of string crossed over itself any number of times, and connected at the ends. The Reidemeister Moves have been proven to be the three core moves necessary to fully untangle a knot. Some knots can be untangled to a loop (the unknot), while others are fundamentally knotted. We define four generalized moves based on the Reidemeister Moves. These moves have the capability of untangling knots which must be made more complicated before they can be simplified. With these moves, we construct a computer program which reads the two-dimensional projection of a knot in its Gauss Code notation and untangles it to the fewest possible number of crossings. Due to the properties of the Gauss Code notation, the program runs efficiently with minimal computation time, compared to currently existing untangling programs. We have tested it on all possible Gauss Codes of up to 50 crossings, and the program successfully determines a simplest possible projection of each knot. The results are consistent with our predictions about how the moves work in untangling a knot.


title: Topologically Protected Dynamics of Spin Textures

abstract: We study current-induced dynamics of spin textures in thin magnetic nanowires. We derive effective equations of motion describing the dynamics of the domain-wall soft modes associated with topological defects. Because the magnetic domain walls are topological objects, these equations are universal and depend only on a few parameters. We obtain spin spiral domain-wall structure in ferromagnetic wires with Dzyaloshinskii-Moriya interaction and critical current dependence on this interaction. We also find the most efficient way to move the domain walls by resonant current pulses and propose a procedure to determine their dynamics by measuring the voltage induced by a moving domain wall. Based on translationally non-invariant nanowires, we show how to make prospective magnetic memory nanodevices much more energy efficient.


title: Delivering SKA Science

abstract: The SKA will be capable of producing a stream of science data products that are Exa-scale in terms of their storage and processing requirements. This Google-scale enterprise is attracting considerable international interest and excitement from within the industrial and academic communities. In this chapter we examine the data flow, storage and processing requirements of a number of key SKA survey science projects to be executed on the baseline SKA1 configuration. Based on a set of conservative assumptions about trends for HPC and storage costs, and the data flow process within the SKA Observatory, it is apparent that survey projects of the scale proposed will potentially drive construction and operations costs beyond the current anticipated SKA1 budget. This implies a sharing of the resources and costs to deliver SKA science between the community and what is contained within the SKA Observatory. A similar situation was apparent to the designers of the LHC more than 10 years ago. We propose that it is time for the SKA project and community to consider the effort and process needed to design and implement a distributed SKA science data system that leans on the lessons of other projects and looks to recent developments in Cloud technologies to ensure an affordable, effective and global achievement of SKA science goals.


title: Spatial Spectrum and Energy Efficiency of Random Cellular Networks

abstract: It is a great challenge to evaluate the network performance of cellular mobile communication systems. In this paper, we propose new spatial spectrum and energy efficiency models for Poisson-Voronoi tessellation (PVT) random cellular networks. To evaluate the user access the network, a Markov chain based wireless channel access model is first proposed for PVT random cellular networks. On that basis, the outage probability and blocking probability of PVT random cellular networks are derived, which can be computed numerically. Furthermore, taking into account the call arrival rate, the path loss exponent and the base station (BS) density in random cellular networks, spatial spectrum and energy efficiency models are proposed and analyzed for PVT random cellular networks. Numerical simulations are conducted to evaluate the network spectrum and energy efficiency in PVT random cellular networks.


title: An analogue of the L\'{e}vy-Hin\v{c}in formula for bi-free infinitely divisible distributions

abstract: In this paper, we derive the bi-free analogue of the L\'{e}vy-Hin\v{c}in formula for compactly supported planar probability measures which are infinitely divisible with respect to the additive bi-free convolution introduced by Voiculescu. We also provide examples of bi-free infinitely divisible distributions with their bi-free L\'{e}vy-Hin\v{c}in representations. Furthermore, we construct the bi-free L\'{e}vy processes and the additive bi-free convolution semigroups generated by compactly supported planar probability measures.


title: Parametric Estimation from Approximate Data: Non-Gaussian Diffusions

abstract: We study the problem of parameters estimation in Indirect Observability contexts, where $X_t \in R^r$ is an unobservable stationary process parametrized by a vector of unknown parameters and all observable data are generated by an approximating process $Y^{\varepsilon}_t$ which is close to $X_t$ in $L^4$ norm. We construct consistent parameter estimators which are smooth functions of the sub-sampled empirical mean and empirical lagged covariance matrices computed from the observable data. We derive explicit optimal sub-sampling schemes specifying the best paired choices of sub-sampling time-step and number of observations. We show that these choices ensure that our parameter estimators reach optimized asymptotic $L^2$-convergence rates, which are constant multiples of the $L^4$ norm $|| Y^{\varepsilon}_t - X_t ||$.


title: Multistatic Cloud Radar Systems: Joint Sensing and Communication Design

abstract: In a multistatic cloud radar system, receive sensors measure signals sent by a transmit element and reflected from a target and possibly clutter, in the presence of interference and noise. The receive sensors communicate over non-ideal backhaul links with a fusion center, or cloud processor, where the presence or absence of the target is determined. The backhaul architecture can be characterized either by an orthogonal-access channel or by a non-orthogonal multiple-access channel. Two backhaul transmission strategies are considered, namely compress-and-forward (CF), which is well suited for the orthogonal-access backhaul, and amplify-and-forward (AF), which leverages the superposition property of the non-orthogonal multiple-access channel. In this paper, the joint optimization of the sensing and backhaul communication functions of the cloud radar system is studied. Specifically, the transmitted waveform is jointly optimized with backhaul quantization in the case of CF backhaul transmission and with the amplifying gains of the sensors for the AF backhaul strategy. In both cases, the information-theoretic criterion of the Bhattacharyya distance is adopted as a metric for the detection performance. Algorithmic solutions based on successive convex approximation are developed under different assumptions on the available channel state information (CSI). Numerical results demonstrate that the proposed schemes outperform conventional solutions that perform separate optimizations of the waveform and backhaul operation, as well as the standard distributed detection approach.


title: Intrinsic Radiation in Lutetium Based PET Detector: Advantages and Disadvantages

abstract: Lutetium (Lu) based scintillators such as LSO and LYSO, are widely used in modern PET detectors due to their high stopping power for 511 keV gamma rays, high light yield and short decay time. However, 2.6% of naturally occurring Lu is 176Lu, a long-lived radioactive element including a beta decay and three major simultaneous gamma decays. This phenomenon introduces random events to PET systems that affects the system performance. On the other hand, the advantages of intrinsic radiation of 176Lu (IRL) continues to be exploited. In this paper, research literatures about IRL in PET detectors are reviewed. Details about the adverse effects of IRL to PET and their solutions, as well as the useful applications are presented and discussed.


title: $B \to \pi \ell \nu$ and $B_s \to K \ell \nu$ form factors and $|V_{ub}|$ from 2+1-flavor lattice QCD with domain-wall light quarks and relativistic heavy quarks

abstract: We calculate the $B \to\pi\ell\nu$ and $B_s \to K \ell\nu$ form factors in dynamical lattice QCD. We use the (2+1)-flavor RBC-UKQCD gauge-field ensembles generated with the domain-wall fermion and Iwasaki gauge actions. For the $b$ quarks we use the anisotropic clover action with a relativistic heavy-quark interpretation. We analyze two lattice spacings $a \approx 0.11, 0.086$ fm and unitary pion masses as light as $M_\pi \approx 290$ MeV. We simultaneously extrapolate our numerical results to the physical light-quark masses and to the continuum and interpolate in the pion/kaon energy using SU(2) "hard-pion" chiral perturbation theory. We provide complete error budgets for the form factors $f_+(q^2)$ and $f_0(q^2)$ at three momenta that span the $q^2$ range accessible in our numerical simulations. We extrapolate these results to $q^2 = 0$ using a model-independent $z$-parametrization and present our final form factors as the $z$-coefficients and the matrix of correlations between them. Our results agree with other lattice determinations using staggered light quarks and provide important independent cross-checks. Both $B \to\pi\ell\nu$ and $B_s \to K \ell\nu$ decays enable a determination of the CKM matrix element $|V_{ub}|$. To illustrate this, we perform a combined $z$-fit of our numerical $B\to\pi\ell\nu$ form-factor data with the experimental branching-fraction measurements leaving the relative normalization as a free parameter; we obtain $|V_{ub}| = 3.61(32) \times 10^{-3}$, where the error includes statistical and systematic uncertainties. This approach can be applied to $B_s\to K \ell\nu$ decay to determine $|V_{ub}|$ once the process has been measured experimentally. Finally, in anticipation of future measurements, we make predictions for $B \to \pi\ell\nu$ and $B_s\to K \ell\nu$ Standard-Model differential branching fractions and forward-backward asymmetries.


title: Ultralong Faraday laser as an optical frequency standard

abstract: In this letter, we introduce the concept and experimentally demonstrate an ultralong Faraday laser as an optical frequency standard in principle. The ultralong Faraday laser is based on the Faraday anomalous dispersion optical filter (FADOF) with ultra-narrow bandwidth and the ultralong fiber extended cavity of $800$ m. The ultra-narrow FADOF is based on atomic transition line of isotope $^{87}$Rb, which has an ultra-narrow bandwidth of $26.0$ MHz and a transmission of $23.6\%$ at $780$ nm. Fibers of length $150$ m and $800$ m are used as ultralong fiber extended cavities, which provide optical feedback and give extremely small FSR of $0.667$ MHz and $0.125$ MHz, respectively. The mechanism of the proposed ultralong Faraday laser is to combine FADOF's ultra-narrow bandwidth and ultralong cavity's small free spectral range to limit the lasing frequency within FADOF bandwidth covered by the semiconductor gain. The active lasing frequency of the ultralong Faraday laser is determined by the center frequency of FADOF transmission, which is corresponding to atomic transition $5^{2}S_{1/2},\ F = 2 \ \rightarrow\ 5^{2}P_{3/2},\ F^{\prime} = 2,\ 3$ of isotope $^{87}$Rb. The Allan deviation of the fractional frequency of the ultralong Faraday laser output signal in 0.1 s -- 1 s measuring time is around $1\times10^{-11}$.


title: Higgs Inflation, Seesaw Physics and Fermion Dark Matter

abstract: We present an inflationary model in which the Standard Model Higgs doublet field with non-minimal coupling to gravity drives inflation, and the effective Higgs potential is stabilized by new physics which includes a dark matter particle and right-handed neutrinos for the seesaw mechanism. All of the new particles are fermions, so that the Higgs doublet is the unique inflaton candidate. With central values for the masses of the top quark and the Higgs boson, the renormalization group improved Higgs potential is employed to yield the scalar spectral index $n_s \simeq 0.968$, the tensor-to-scalar ratio $r \simeq 0.003$, and the running of the spectral index $\alpha=dn_s/d \ln k \simeq -5.2 \times 10^{-4}$ for the number of e-folds $N_0=60$ ($n_s \simeq 0.962$, $r \simeq 0.004$, and $\alpha \simeq -7.5 \times 10^{-4}$ for $N_0=50$). The fairly low value of $r \simeq 0.003$ predicted in this class of models means that the ongoing space and land based experiments are not expected to observe gravity waves generated during inflation. [Dedicated to the memory of Dr. Paul Weber (1947 - 2015). Paul was an exceptional human being and a very special friend who will be sorely missed.]


title: Wireless Information and Power Transfer in Relay Systems with Multiple Antennas and Interference

abstract: In this paper, an energy harvesting dual-hop relaying system without/with the presence of co-channel interference (CCI) is investigated. Specifically, the energy constrained multi-antenna relay node is powered by either the information signal of the source or via the signal receiving from both the source and interferer. In particular, we first study the outage probability and ergodic capacity of an interference free system, and then extend the analysis to an interfering environment. To exploit the benefit of multiple antennas, three different linear processing schemes are investigated, namely, 1) Maximum ratio combining/maximal ratio transmission (MRC/MRT), 2) Zero-forcing/MRT (ZF/MRT) and 3) Minimum mean-square error/MRT (MMSE/MRT). For all schemes, both the systems outage probability and ergodic capacity are studied, and the achievable diversity order is also presented. In addition, the optimal power splitting ratio minimizing the outage probability is characterized. Our results show that the implementation of multiple antennas increases the energy harvesting capability, hence, significantly improves the systems performance. Moreover, it is demonstrated that the CCI could be potentially exploited to substantially boost the performance, while the choice of a linear processing scheme plays a critical role in determining how much gain could be extracted from the CCI.


title: Modern Studies of the Deuteron: from the Lab Frame to the Light Front

abstract: We review the recent progress made in studies of deuteron structure at small internucleon distances. This progress is largely facilitated by the new generation of experiments in deuteron electrodisintegration carried out at unprecedentedly high momentum transfer. The theoretical analysis of these data confirms the onset of the high energy eikonal regime in the scattering process which allows one to separate long range nuclear effects from the effects genuinely related to the short distance structure of the deuteron. Our conclusion is that for the first time the deuteron is probed at relative momenta beyond 300 MeV/c without dominating long range effects. As a result, at these large nucleon momenta the cross section is sensitive to the nuclear dynamics at sub-fermi distances. Due to large internal momenta involved we are dealing with the relativistic bound state that is best described by the light-cone momentum distribution of nucleons in the deuteron. We present the first attempt of extracting the deuteron light-cone momentum distribution function from data and discuss the importance of this quantity for studies of QCD structure of the bound nucleon in deep inelastic scattering off the deuteron. We conclude the review giving an outlook of the next generation of high energy experiments which will extend our reach to much smaller distances in the deuteron.


title: Generalized Yang-Mills Theory and Gravity

abstract: We propose a generalization of Yang-Mills theory for which the symmetry algebra does not have to be factorized as mutually commuting algebras of a finite-dimensional Lie algebra and the algebra of functions on base space. The algebra of diffeomorphism can be constructed as an example, and a class of gravity theories can be interpreted as generalized Yang-Mills theories. These theories in general include a graviton, a dilaton and a rank-2 antisymmetric field, although Einstein gravity is also included as a special case. We present calculations suggesting that the connection in scattering amplitudes between Yang-Mills theory and gravity via BCJ duality can be made more manifest in this formulation.


title: Communication Theoretic Data Analytics

abstract: Widespread use of the Internet and social networks invokes the generation of big data, which is proving to be useful in a number of applications. To deal with explosively growing amounts of data, data analytics has emerged as a critical technology related to computing, signal processing, and information networking. In this paper, a formalism is considered in which data is modeled as a generalized social network and communication theory and information theory are thereby extended to data analytics. First, the creation of an equalizer to optimize information transfer between two data variables is considered, and financial data is used to demonstrate the advantages. Then, an information coupling approach based on information geometry is applied for dimensionality reduction, with a pattern recognition example to illustrate the effectiveness. These initial trials suggest the potential of communication theoretic data analytics for a wide range of applications.


title: Quasi-Periodic Schr\"odinger Cocycles with Positive Lyapunov Exponent are not Open in the Smooth Topology

abstract: One knows that the set of quasi-periodic Schr\"odinger cocycles with positive Lyapunov exponent is open and dense in the analytic topology. In this paper, we construct cocycles with positive Lyapunov exponent which can be approximated by ones with zero Lyapunov exponent in the space of ${\cal C}^ l$ ($1 \le l\le \infty$) smooth quasi-periodic cocycles. As a consequence, the set of quasi-periodic Schr\"odinger cocycles with positive Lyapunov exponent is not ${\cal C}^ l$ open.


title: Combining Alphas via Bounded Regression

abstract: We give an explicit algorithm and source code for combining alpha streams via bounded regression. In practical applications typically there is insufficient history to compute a sample covariance matrix (SCM) for a large number of alphas. To compute alpha allocation weights, one then resorts to (weighted) regression over SCM principal components. Regression often produces alpha weights with insufficient diversification and/or skewed distribution against, e.g., turnover. This can be rectified by imposing bounds on alpha weights within the regression procedure. Bounded regression can also be applied to stock and other asset portfolio construction. We discuss illustrative examples.