title: Graph Memory Networks for Molecular Activity Prediction

abstract: Molecular activity prediction is critical in drug design. Machine learning techniques such as kernel methods and random forests have been successful for this task. These models require fixed-size feature vectors as input while the molecules are variable in size and structure. As a result, fixed-size fingerprint representation is poor in handling substructures for large molecules. In addition, molecular activity tests, or a so-called BioAssays, are relatively small in the number of tested molecules due to its complexity. Here we approach the problem through deep neural networks as they are flexible in modeling structured data such as grids, sequences and graphs. We train multiple BioAssays using a multi-task learning framework, which combines information from multiple sources to improve the performance of prediction, especially on small datasets. We propose Graph Memory Network (GraphMem), a memory-augmented neural network to model the graph structure in molecules. GraphMem consists of a recurrent controller coupled with an external memory whose cells dynamically interact and change through a multi-hop reasoning process. Applied to the molecules, the dynamic interactions enable an iterative refinement of the representation of molecular graphs with multiple bond types. GraphMem is capable of jointly training on multiple datasets by using a specific-task query fed to the controller as an input. We demonstrate the effectiveness of the proposed model for separately and jointly training on more than 100K measurements, spanning across 9 BioAssay activity tests.


title: High rate tests of the photon detection system for the LHCb RICH Upgrade

abstract: The photon detection system for the LHCb RICH Upgrade consists of an array of multianode photomultiplier tubes (MaPMTs) read out by custom-built modular electronics. The behaviour of the whole chain was studied at CERN using a pulsed laser. Threshold scans were performed in order to study the MaPMT pulse-height spectra at high event rates and different photon intensities. The results show a reduction in photon detection efficiency at 900 V bias voltage, marked by a 20 % decrease in the single-photon peak height, when increasing the event rate from 100 kHz to 20 MHz. This reduction was not observed at 1000 V bias voltage.


title: Disks Around Merging Binary Black Holes: From GW150914 to Supermassive Black Holes

abstract: We perform magnetohydrodynamic simulations in full general relativity of disk accretion onto nonspinning black hole binaries with mass ratio 36:29. We survey different disk models which differ in their scale height, total size and magnetic field to quantify the robustness of previous simulations on the initial disk model. Scaling our simulations to LIGO GW150914 we find that such systems could explain possible gravitational wave and electromagnetic counterparts such as the Fermi GBM hard X-ray signal reported 0.4s after GW150915 ended. Scaling our simulations to supermassive binary black holes, we find that observable flow properties such as accretion rate periodicities, the emergence of jets throughout inspiral, merger and post-merger, disk temperatures, thermal frequencies, and the time-delay between merger and the boost in jet outflows that we reported in earlier studies display only modest dependence on the initial disk model we consider here.

categories: astro-ph.HE astro-ph.GA astro-ph.SR gr-qc


title: Ising versus $SU(2)_2$ string-net ladder

abstract: We consider the string-net model obtained from $SU(2)_2$ fusion rules. These fusion rules are shared by two different sets of anyon theories. In this work, we study the competition between the two corresponding non-Abelian quantum phases in the ladder geometry. A detailed symmetry analysis shows that the nontrivial low-energy sector corresponds to the transverse-field cluster model that displays a critical point described by the $so(2)_1$ conformal field theory. Other sectors are obtained by freezing spins in this model.


title: The Antesonic Condition for the Explosion of Core-Collapse Supernovae I: Spherically Symmetric Polytropic Models: Stability & Wind Emergence

abstract: Shock revival in core-collapse supernovae (CCSNe) may be due to the neutrino mechanism. While it is known that in a neutrino-powered CCSN, explosion begins when the neutrino luminosity of the proto-neutron star exceeds a critical value, the physics of this condition in time-dependent, multidimensional simulations are not fully understood. \citet{Pejcha2012} found that an `antesonic condition' exists for time-steady spherically symmetric models, potentially giving a physical explanation for the critical curve observed in simulations. In this paper, we extend that analysis to time-dependent, spherically symmetric polytropic models. We verify the critical antesonic condition in our simulations, showing that models exceeding it drive transonic winds whereas models below it exhibit steady accretion. In addition, we find that (1) high spatial resolution is needed for accurate determination of the antesonic ratio and shock radius at the critical curve, and that low resolution simulations systematically underpredict these quantities, making explosion more difficult at lower resolution; (2) there is an important physical connection between the critical mass accretion rate at explosion and the mass loss rate of the post-explosion wind: the two are directly proportional at criticality, implying that, at criticality, the wind kinetic power is tied directly to the accretion power; (3) the value of the post-shock adiabatic index $\Gamma$ has a large effect on the length and time scales of the post-bounce evolution of the explosion larger values of $\Gamma$ result in a longer transition from the accretion to wind phases.


title: Simplicity from Recoil: The Three-Loop Soft Function and Factorization for the Energy-Energy Correlation

abstract: We derive an operator based factorization theorem for the energy-energy correlation (EEC) observable in the back-to-back region, allowing the cross section to be written as a convolution of hard, jet and soft functions. We prove the equivalence of the soft functions for the EEC and color singlet transverse-momentum resummation to all-loop order, and give their analytic result to three-loops. Large logarithms appearing in the perturbative expansion of the EEC can be resummed to all orders using renormalization group evolution. We give analytic results for all required anomalous dimensions to three-loop order, providing the first example of a transverse-momentum (recoil) sensitive $e^+e^-$ event shape whose anomalous dimensions are known at this order. The EEC can now be computed to next-to-next-to-next-to-leading logarithm matched to next-to-next-to-leading order, making it a prime candidate for precision QCD studies and extractions of the strong coupling constant. We anticipate that our factorization theorem will also be crucial for understanding non-perturbative power corrections for the EEC, and their relationship to those appearing in other observables.


title: Classical Discrete Time Crystals

abstract: The spontaneous breaking of time-translation symmetry in periodically driven quantum systems leads to a new phase of matter: discrete time crystals (DTC). This phase exhibits collective subharmonic oscillations that depend upon an interplay of non-equilibrium driving, many-body interactions, and the breakdown of ergodicity. However, subharmonic responses are also a well-known feature of classical dynamical systems ranging from predator-prey models to Faraday waves and AC-driven charge density waves. This raises the question of whether these classical phenomena display the same rigidity characteristic of a quantum DTC. In this work, we explore this question in the context of periodically driven Hamiltonian dynamics coupled to a finite-temperature bath, which provides both friction and, crucially, noise. Focusing on one-dimensional chains, where in equilibrium any transition would be forbidden at finite temperature, we provide evidence that the combination of noise and interactions drives a sharp, first-order dynamical phase transition between a discrete time-translation invariant phase and an activated classical discrete time crystal (CDTC) in which time-translation symmetry is broken out to exponentially-long time scales. Power-law correlations are present along a first-order line which terminates at a critical point. We analyze the transition by mapping it to the locked-to-sliding transition of a DC-driven charge density wave. Our work points to a classical limit for quantum time crystals, and raises several intriguing questions concerning the non-equilibrium universality class of the CDTC critical point.

categories: cond-mat.stat-mech cond-mat.dis-nn cond-mat.mes-hall cond-mat.str-el quant-ph


title: LOFAR/H-ATLAS: The low-frequency radio luminosity - star-formation rate relation

abstract: Radio emission is a key indicator of star-formation activity in galaxies, but the radio luminosity-star formation relation has to date been studied almost exclusively at frequencies of 1.4 GHz or above. At lower radio frequencies the effects of thermal radio emission are greatly reduced, and so we would expect the radio emission observed to be completely dominated by synchrotron radiation from supernova-generated cosmic rays. As part of the LOFAR Surveys Key Science project, the Herschel-ATLAS NGP field has been surveyed with LOFAR at an effective frequency of 150 MHz. We select a sample from the MPA-JHU catalogue of SDSS galaxies in this area: the combination of Herschel, optical and mid-infrared data enable us to derive star-formation rates (SFRs) for our sources using spectral energy distribution fitting, allowing a detailed study of the low-frequency radio luminosity--star-formation relation in the nearby Universe. For those objects selected as star-forming galaxies (SFGs) using optical emission line diagnostics, we find a tight relationship between the 150 MHz radio luminosity ($L_{150}$) and SFR. Interestingly, we find that a single power-law relationship between $L_{150}$ and SFR is not a good description of all SFGs: a broken power law model provides a better fit. This may indicate an additional mechanism for the generation of radio-emitting cosmic rays. Also, at given SFR, the radio luminosity depends on the stellar mass of the galaxy. Objects which were not classified as SFGs have higher 150-MHz radio luminosity than would be expected given their SFR, implying an important role for low-level active galactic nucleus activity.


title: Effective theory approach to the Schrodinger-Poisson problem in semiconductor Majorana devices

abstract: We propose a method for solving the Schrodinger-Poisson problem that can be efficiently implemented in realistic 3D tight-binding models of semiconductor-based Majorana devices. The method is based on two key ideas: i) For a given geometry, the Poisson problem is only solved once (for each local orbital) and the results are stored as an interaction tensor; using this Green's function scheme, the Poisson component of the iteration procedure is reduced to a few simple summations. ii) The 3D problem is mapped into an effective multi-orbital 1D problem with molecular orbitals calculated self-consistently as the transverse modes of an infinite wire with the same electrostatic potential as the local electrostatic potential of the finite 3D device. To demonstrate the capabilities of our approach, we calculate the response of the system to an external magnetic field, the dependence of the effective chemical potential on the work function difference, and the dependence of the effective semiconductor-superconductor coupling on the applied gate potential. We find that different low-energy bands are characterized by similar effective couplings, which results in induced gap features characterized by a single energy scale. We also find that electrostatic effects are responsible for a partial suppression of the Majorana energy splitting oscillations. Finally, we show that a position-dependent work function difference can produce a non-homogeneous effective potential that is not affected by the screening due to the superconductor and is only partially suppressed by the charge inside the wire. In turn, this potential can induce trivial low-energy states that mimic the phenomenology of Majorana zero modes. Thus, any position-dependent work function difference (even at the 1% level) along the nanowire must be avoided through carefully engineered semiconductor-superconductor interfaces.


title: Variations of the stellar initial mass function in semi-analytical models II: the impact of Cosmic Ray regulation

abstract: Recent studies proposed that cosmic rays (CR) are a key ingredient in setting the conditions for star formation, thanks to their ability to alter the thermal and chemical state of dense gas in the UV-shielded cores of molecular clouds. In this paper, we explore their role as regulators of the stellar initial mass function (IMF) variations, using the semi-analytic model for GAlaxy Evolution and Assembly (GAEA). The new model confirms our previous results obtained using the integrated galaxy-wide IMF (IGIMF) theory: both variable IMF models reproduce the observed increase of $\alpha$-enhancement as a function of stellar mass and the measured $z=0$ excess of dynamical mass-to-light ratios with respect to photometric estimates assuming a universal IMF. We focus here on the mismatch between the photometrically-derived ($M^{\rm app}_{\star}$) and intrinsic ($M_{\star}$) stellar masses, by analysing in detail the evolution of model galaxies with different values of $M_{\star}/M^{\rm app}_{\star}$. We find that galaxies with small deviations (i.e. formally consistent with a universal IMF hypothesis) are characterized by more extended star formation histories and live in less massive haloes with respect to the bulk of the galaxy population. While the IGIMF theory does not change significantly the mean evolution of model galaxies with respect to the reference model, a CR-regulated IMF implies shorter star formation histories and higher peaks of star formation for objects more massive than $10^{10.5} M_\odot$. However, we also show that it is difficult to unveil this behaviour from observations, as the key physical quantities are typically derived assuming a universal IMF.


title: Approaches to linear local gauge-invariant observables in inflationary cosmologies

abstract: We review and relate two recent complementary constructions of linear local gauge-invariant observables for cosmological perturbations in generic spatially flat single-field inflationary cosmologies. After briefly discussing their physical significance, we give explicit, covariant and mutually invertible transformations between the two sets of observables, thus resolving any doubts about their equivalence. In this way, we get a geometric interpretation and show the completeness of both sets of observables, while previously each of these properties was available only for one of them.


title: The Next Generation Fornax Survey (NGFS): II. The Central Dwarf Galaxy Population

abstract: We present a photometric study of the dwarf galaxy population in the core region ($< r_{\rm vir}/4$) of the Fornax galaxy cluster based on deep $u'g'i'$ photometry from the Next Generation Fornax Cluster Survey. All imaging data were obtained with the Dark Energy Camera mounted on the 4-meter Blanco telescope at the Cerro-Tololo Interamerican Observatory. We identify 258 dwarf galaxy candidates with luminosities $-17 < M_{g'} < -8$ mag, corresponding to typical stellar masses of $9.5\gtrsim \log{\cal M}_{\star}/M_\odot \gtrsim 5.5$, reaching $\sim\!3$ mag deeper in point-source luminosity and $\sim\!4$ mag deeper in surface-brightness sensitivity compared to the classic Fornax Cluster Catalog. Morphological analysis shows that surface-brightness profiles are well represented by single-component S\'ersic models with average S\'ersic indices of $\langle n\rangle_{u',g',i'}=(0.78-0.83) \pm 0.02$, and average effective radii of $\langle r_e\rangle_{u',g',i'}\!=(0.67-0.70) \pm 0.02$ kpc. Color-magnitude relations indicate a flattening of the galaxy red sequence at faint galaxy luminosities, similar to the one recently discovered in the Virgo cluster. A comparison with population synthesis models and the galaxy mass-metallicity relation reveals that the average faint dwarf galaxy is likely older than ~5 Gyr. We study galaxy scaling relations between stellar mass, effective radius, and stellar mass surface density over a stellar mass range covering six orders of magnitude. We find that over the sampled stellar mass range several distinct mechanisms of galaxy mass assembly can be identified: i) dwarf galaxies assemble mass inside the half-mass radius up to $\log{\cal M}_{\star}$ ~8.0, ii) isometric mass assembly in the range $8.0 < \log{\cal M}_{\star}/M_\odot < 10.5$, and iii) massive galaxies assemble stellar mass predominantly in their halos at $\log{\cal M}_{\star}$ ~10.5 and above.


title: The Astropy Project: Building an inclusive, open-science project and status of the v2.0 core package

abstract: The Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project.


title: UV-Luminous, Star-Forming Hosts of z~2 Reddened Quasars in the Dark Energy Survey

abstract: We present the first rest-frame UV population study of 17 heavily reddened, high-luminosity (E(B-V)$_{\rm{QSO}}\gtrsim$ 0.5; L$_{\rm{bol}}>$ 10$^{46}$ergs$^{-1}$) broad-line quasars at $1.5 < z < 2.7$. We combine the first year of deep, optical, ground-based observations from the Dark Energy Survey (DES) with the near infrared VISTA Hemisphere Survey (VHS) and UKIDSS Large Area Survey (ULAS) data, from which the reddened quasars were initially identified. We demonstrate that the significant dust reddening towards the quasar in our sample allows host galaxy emission to be detected at the rest-frame UV wavelengths probed by the DES photometry. By exploiting this reddening effect, we disentangle the quasar emission from that of the host galaxy via spectral energy distribution (SED) fitting. We find evidence for a relatively unobscured, star-forming host galaxy in at least ten quasars, with a further three quasars exhibiting emission consistent with either star formation or scattered light. From the rest-frame UV emission, we derive instantaneous, dust-corrected star formation rates (SFRs) in the range 25 < SFR$_{\rm{UV}}$ < 365 M$_{\odot}$yr$^{-1}$, with an average SFR$_{\rm{UV}}$ = 130 $\pm$ 95 M$_{\odot}$yr$^{-1}$. We find a broad correlation between SFR$_{\rm{UV}}$ and the bolometric quasar luminosity. Overall, our results show evidence for coeval star formation and black hole accretion occurring in luminous, reddened quasars at the peak epoch of galaxy formation.


title: Imaging the water snowline in a protostellar envelope with H$^{13}$CO$^+$

abstract: Snowlines are key ingredients for planet formation. Providing observational constraints on the locations of the major snowlines is therefore crucial for fully connecting planet compositions to their formation mechanism. Unfortunately, the most important snowline, that of water, is very difficult to observe directly in protoplanetary disks due to its close proximity to the central star. Based on chemical considerations, HCO$^+$ is predicted to be a good chemical tracer of the water snowline, because it is particularly abundant in dense clouds when water is frozen out. This work maps the optically thin isotopologue H$^{13}$CO$^+$ ($J=3-2$) toward the envelope of the low-mass protostar NGC1333-IRAS2A (observed with NOEMA at ~0.9" resolution), where the snowline is at larger distance from the star than in disks. The H$^{13}$CO$^+$ emission peaks ~2" northeast of the continuum peak, whereas the previously observed H$_2^{18}$O shows compact emission on source. Quantitative modeling shows that a decrease in H$^{13}$CO$^+$ abundance by at least a factor of six is needed in the inner ~360 AU to reproduce the observed emission profile. Chemical modeling predicts indeed a steep increase in HCO$^+$ just outside the water snowline; the 50% decrease in gaseous H$_2$O at the snowline is not enough to allow HCO$^+$ to be abundant. This places the water snowline at 225 AU, further away from the star than expected based on the 1D envelope temperature structure for NGC1333-IRAS2A. In contrast, DCO$^+$ observations show that the CO snowline is at the expected location, making an outburst scenario unlikely. The spatial anticorrelation of the H$^{13}$CO$^+$ and H$_2^{18}$O emission provide a proof of concept that H$^{13}$CO$^+$ can be used as a tracer of the water snowline.


title: Scale-invariance and scale-breaking in parity-invariant three-dimensional QCD

abstract: We present a numerical study of three-dimensional two-color QCD with $N=0, 2, 4, 8$ and 12 flavors of massless two-component fermions using parity-preserving improved Wilson-Dirac fermions. A finite volume analysis provides strong evidence for the presence of $Sp(N)$ symmetry-breaking bilinear condensate when $N \le 2$ and its absence for $N \ge 8$. A weaker evidence for the bilinear condensate is shown for $N=4$. We estimate the critical number of flavors below which scale-invariance is broken by the bilinear condensate to be between $N=4$ and 6.


title: Vaex: Big Data exploration in the era of Gaia

abstract: We present a new Python library called vaex, to handle extremely large tabular datasets, such as astronomical catalogues like the Gaia catalogue, N-body simulations or any other regular datasets which can be structured in rows and columns. Fast computations of statistics on regular N-dimensional grids allows analysis and visualization in the order of a billion rows per second. We use streaming algorithms, memory mapped files and a zero memory copy policy to allow exploration of datasets larger than memory, e.g. out-of-core algorithms. Vaex allows arbitrary (mathematical) transformations using normal Python expressions and (a subset of) numpy functions which are lazily evaluated and computed when needed in small chunks, which avoids wasting of RAM. Boolean expressions (which are also lazily evaluated) can be used to explore subsets of the data, which we call selections. Vaex uses a similar DataFrame API as Pandas, a very popular library, which helps migration from Pandas. Visualization is one of the key points of vaex, and is done using binned statistics in 1d (e.g. histogram), in 2d (e.g. 2d histograms with colormapping) and 3d (using volume rendering). Vaex is split in in several packages: vaex-core for the computational part, vaex-viz for visualization mostly based on matplotlib, vaex-jupyter for visualization in the Jupyter notebook/lab based in IPyWidgets, vaex-server for the (optional) client-server communication, vaex-ui for the Qt based interface, vaex-hdf5 for hdf5 based memory mapped storage, vaex-astro for astronomy related selections, transformations and memory mapped (column based) fits storage. Vaex is open source and available under MIT license on github, documentation and other information can be found on the main website: https://vaex.io, https://docs.vaex.io or https://github.com/maartenbreddels/vaex


title: Larkin-Ovchinikov superfluidity in time-reversal symmetric bilayer Fermi gases

abstract: Larkin-Ovchinnikov (LO) state which combines the superfluidity and spatial periodicity of pairing order parameter and exhibits the supersolid properties has been attracting intense attention in both condensed matter physics and ultracold atoms. Conventionally, realization of LO state from an intrinsic s-wave interacting system necessitates to break the time-reversal (TR) and sometimes spatial-inversion (SI) symmetries. Here we report a novel prediction that the LO state can be realized in a TR and SI symmetric system representing a bilayer Fermi gas subjected to a laserassisted interlayer tunneling. We show that the intralayer s-wave atomic interaction acts effectively like a p-wave interaction in the pseudospin space. This provides distinctive pairing effects in the present system with pseudspin spin-orbit coupling, and leads to a spontaneous density-modulation of the pairing order predicted in a very broad parameter regime. Unlike the conventional schemes, our results do not rely on the spin imbalance or external Zeeman fields, showing a highly feasible way to observe the long-sought-after LO superfluid phase using the laser-assisted bilayer Fermi gases.


title: Exotic Sterile Neutrinos and Pseudo-Goldstone Phenomenology

abstract: We study the phenomenology of a light (GeV scale) sterile neutrino sector and the pseudo-Goldstone boson (not the majoron) associated with a global symmetry in this sector that is broken at a high scale. Such scenarios can be motivated from considerations of singlet fermions from a hidden sector coupling to active neutrinos via heavy right-handed seesaw neutrinos, effectively giving rise to a secondary, low-energy seesaw framework. Such scenarios allow for rich phenomenology with observable implications for cosmology, dark matter, and direct searches, involving novel sterile neutrino dark matter production mechanisms from the pseudo-Goldstone-mediated scattering or decay, modifications of BBN bounds on sterile neutrinos, suppression of canonical sterile neutrino decay channels at direct search experiments, late injection of an additional population of neutrinos in the Universe after neutrino decoupling, and measurable dark radiation at BBN or CMB decoupling.


title: An ALMA [CII] survey of 27 quasars at z>5.94

abstract: We present a survey of the [CII] 158 $\mu$m line and underlying far-infrared (FIR) dust continuum emission in a sample of 27 z>6 quasars using the Atacama Large Millimeter Array (ALMA) at ~1" resolution. The [CII] line was significantly detected (at >5-sigma) in 23 sources (85%). We find typical line luminosities of $L_{\rm [CII]}=10^{9-10}$ L$_\odot$, and an average line width of ~385 km/s. The [CII]-to-far-infrared luminosity ratio ([CII]/FIR) in our sources span one order of magnitude, highlighting a variety of conditions in the star-forming medium. Four quasar host galaxies are clearly resolved in their [CII] emission on a few kpc scales. Basic estimates of the dynamical masses of the host galaxies give masses between $2\times10^{10}$ and $2\times10^{11}$ M$_\odot$, i.e., more than an order of magnitude below what is expected from local scaling relations, given the available limits on the masses of the central black holes ($>3\times10^8$ M$_\odot$, assuming Eddington-limited accretion). In stacked ALMA [CII] spectra of individual sources in our sample, we find no evidence of a deviation from a single Gaussian profile. The quasar luminosity does not strongly correlate with either the [CII] luminosity or equivalent width. This survey (with typical on-source integration times of 8 min) showcases the unparalleled sensitivity of ALMA at millimeter wavelengths, and offers a unique reference sample for the study of the first massive galaxies in the universe.


title: Boundary Optimizing Network (BON)

abstract: Despite all the success that deep neural networks have seen in classifying certain datasets, the challenge of finding optimal solutions that generalize still remains. In this paper, we propose the Boundary Optimizing Network (BON), a new approach to generalization for deep neural networks when used for supervised learning. Given a classification network, we propose to use a collaborative generative network that produces new synthetic data points in the form of perturbations of original data points. In this way, we create a data support around each original data point which prevents decision boundaries from passing too close to the original data points, i.e. prevents overfitting. We show that BON improves convergence on CIFAR-10 using the state-of-the-art Densenet. We do however observe that the generative network suffers from catastrophic forgetting during training, and we therefore propose to use a variation of Memory Aware Synapses to optimize the generative network (called BON++). On the Iris dataset, we visualize the effect of BON++ when the generator does not suffer from catastrophic forgetting and conclude that the approach has the potential to create better boundaries in a higher dimensional space.


title: Evidence for magnetic-field-induced decoupling of superconducting bilayers in La$_{2-x}$Ca$_{1+x}$Cu$_{2}$O$_{6}$

abstract: We report a study of magnetic susceptibility and electrical resistivity as a function of temperature and magnetic field in superconducting crystals of La$_{2-x}$Ca$_{1+x}$Cu$_{2}$O$_{6}$ with $x=0.10$ and 0.15 and transition temperature $T_{c}^{\rm m} = 54$ K (determined from the susceptibility). When an external magnetic field is applied perpendicular to the CuO$_2$ bilayers, the resistive superconducting transition measured with currents flowing perpendicular to the bilayers is substantially lower than that found with currents flowing parallel to the bilayers. Intriguingly, this anisotropic behavior is quite similar to that observed for the magnetic irreversibility points with the field applied either perpendicular or parallel to the bilayers. We discuss the results in the context of other studies that have found evidence for the decoupling of superconducting layers induced by a perpendicular magnetic field.


title: On monomial ideals and their socles

abstract: For a finite subset $M\subset [x_1,\ldots,x_d]$ of monomials, we describe how to constructively obtain a monomial ideal $I\subseteq R = K[x_1,\ldots,x_d]$ such that the set of monomials in $\text{Soc}(I)\setminus I$ is precisely $M$, or such that $\overline{M}\subseteq R/I$ is a $K$-basis for the the socle of $R/I$. For a given $M$ we obtain a natural class of monomials $I$ with this property. This is done by using solely the lattice structure of the monoid $[x_1,\ldots,x_d]$. We then present some duality results by using anti-isomorphisms between upsets and downsets of $(\mathbb Z^d,\preceq)$. Finally, we define and analyze zero-dimensional monomial ideals of $R$ of type $k$, where type $1$ are exactly the Artinian Gorenstein ideals, and describe the structure of such ideals that correspond to order-generic antichains in $\mathbb Z^d$.


title: Personalizing deep learning models for automatic sleep staging

abstract: Despite continued advancement in machine learning algorithms and increasing availability of large data sets, there is still no universally acceptable solution for automatic sleep staging of human sleep recordings. One reason is that a skilled neurophysiologist scoring brain recordings of a sleeping person implicitly adapts his/her staging to the individual characteristics present in the brain recordings. Trying to incorporate this adaptation step in an automatic scoring algorithm, we introduce in this paper a method for personalizing a general sleep scoring model. Starting from a general convolutional neural network architecture, we allow the model to learn individual characteristics of the first night of sleep in order to quantify sleep stages of the second night. While the original neural network allows to sleep stage on a public database with a state of the art accuracy, personalizing the model further increases performance (on the order of two percentage points on average, but more for difficult subjects). This improvement is particularly present in subjects where the original algorithm did not perform well (typically subjects with accuracy less than $80\%$). Looking deeper, we find that optimal classification can be achieved when broad knowledge of sleep staging in general (at least 20 separate nights) is combined with subject-specific knowledge. We hypothesize that this method will be very valuable for improving scoring of lower quality sleep recordings, such as those from wearable devices.


title: Exploiting random lead times for significant inventory cost savings

abstract: We study the classical single-item inventory system in which unsatisfied demands are backlogged. Replenishment lead times are random, independent identically distributed, causing orders to cross in time. We develop a new inventory policy to exploit implications of lead time randomness and order crossover, and evaluate its performance by asymptotic analysis and simulations. Our policy does not follow the basic principle of Constant Base Stock (CBS) policy, or more generally, (s,S) and (r,Q) policies, which is to keep the inventory position within a fixed range. Instead, it uses the current inventory level (= inventory-on-hand minus backlog) to set a dynamic target for inventory in-transit, and place orders to follow this target. Our policy includes CBS policy as a special case, under a particular choice of a policy parameter. We show that our policy can significantly reduce the average inventory cost compared with CBS policy. Specifically, we prove that if the lead time is exponentially distributed, then under our policy, with properly chosen policy parameters, the expected (absolute) inventory level scales as $o(\sqrt{r})$, as the demand rate $r\to\infty$. In comparison, it is known to scale as $\Theta(\sqrt{r})$ under CBS policy. In particular, this means that, as $r\to\infty$, the average inventory cost under our policy vanishes in comparison with that under CBS policy. Furthermore, our simulations show that the advantage of our policy remains to be substantial under non-exponential lead time distributions, and may even be greater than under exponential distribution. We also use simulations to compare GBS to an optimal policy for some cases where computing the optimal cost is tractable. The results show that our policy removes a majority of excess costs of CBS policy over the minimum cost, leading to much smaller optimality gaps.


title: Millimeter mapping at z~1: dust-obscured bulge building and disk growth

abstract: A randomly chosen star in today's Universe is most likely to live in a galaxy with a stellar mass between that of the Milky Way and Andromeda. Yet it remains uncertain how the structural evolution of these bulge-disk systems proceeded. Most of the unobscured star formation we observe building Andromdeda progenitors at 0.7<z<1.5 occurs in disks, but >90% of their star formation is reprocessed by dust and remains unaccounted for. Here we map 500micron dust continuum emission in an Andromeda progenitor at z=1.25 to probe where it is growing through dust-obscured star formation. Combining resolved dust measurements from the NOEMA interferometer with Hubble Space Telescope Halpha maps and multicolor imaging (including new UV data from the HDUV survey), we find a bulge growing by dust-obscured star formation: while the unobscured star formation is centrally suppressed, the dust continuum is centrally concentrated, filling in the ring-like structures evident in the Halpha and UV emission. Reflecting this, the dust emission is more compact than the optical/UV tracers of star formation with r_e(dust)=3.4kpc, r_e(Halpha)/r_e(dust)=1.4, and r_e(UV)/r_e(dust)=1.8. Crucially, however, the bulge and disk of this galaxy are building simultaneously; although the dust emission is more compact than the rest-optical emission (r_e(optical)/r_e(dust)=1.4), it is somewhat less compact than the stellar mass (r_e(M_*)/r_e(dust)=0.9). Taking the 500micron emission as a tracer of star formation, the expected structural evolution of this galaxy can be accounted for by star formation: it will grow in size by Delta(r_e)/Delta(M_*)~0.3 and central surface density by Delta(Sigma_cen)/Delta(M_*)~0.9. Finally, our observations are consistent with a picture in which merging and disk instabilities drive gas to the center of galaxies, boosting global star formation rates above the main sequence and building bulges.


title: Current Status of Warm Inflation

abstract: Warm inflation is an inflationary scenario in which a thermal bath coexists with the inflaton during inflation. This is unlike standard cold inflation in which the Universe is effectively devoid of particles during inflation. The thermal bath in warm inflation is maintained by the dissipation of the inflaton's energy through its couplings to other fields. Many models of warm inflation have been proposed and their predictions have been compared with cosmological data. Certain models of inflation that are disallowed in the context of cold inflation by the data are allowed in the warm inflationary scenario, and vice versa.


title: Reionization and Cosmic Dawn: theory and simulations

abstract: We highlight recent progress in the sophistication and diversification of cosmic dawn and reionization simulations. The application of these modeling tools to current observations has allowed us narrow down the timing of reionization, which we now know to within dz ~ 1 for the bulk of reionization. The strongest constraints come from the optical depth to the CMB measured with the {\it Planck} satellite and the first detection of ongoing reionization from the spectra of the z=7.1 QSOs ULASJ1120+0641. However, we still know virtually nothing about the astrophysical sources during the first billion years. The revolution in our understanding will be led by upcoming interferometric observations of the cosmic 21-cm signal. The properties of the sources and sinks of UV and X-ray photons are encoded in the 3D patterns of the signal. The development of Bayesian parameter recovery techniques, which tap into the wealth of the 21-cm signal, will soon usher in an era of precision astrophysical cosmology.


title: Higher order recurrences and row sequences of Hermite-Pad\'e approximation

abstract: We obtain extensions of the Poincar\'e and Perron theorems for higher order recurrence relations and apply them to obtain an inverse type theorem for row sequences of (type II) Hermite-Pad\'e approximation of a vector of formal power series.


title: Towards General Distributed Resource Selection

abstract: The advantages of distributing workloads and utilizing multiple distributed resources are now well established. The type and degree of heterogeneity of distributed resources is increasing, and thus determining how to distribute the workloads becomes increasingly difficult, in particular with respect to the selection of suitable resources. We formulate and investigate the resource selection problem in a way that it is agnostic of specific task and resource properties, and which is generalizable to range of metrics. Specifically, we developed a model to describe the requirements of tasks and to estimate the cost of running that task on an arbitrary resource using baseline measurements from a reference machine. We integrated our cost model with the Condor matchmaking algorithm to enable resource selection. Experimental validation of our model shows that it provides execution time estimates with 157-171% error on XSEDE resources and 18-31% on OSG resources. We use the task execution cost model to select resources for a bag-of-tasks of up to 1024 GROMACS MD simulations across the target resources. Experiments show that using the model's estimates reduces the workload's time-to-completion up to ~85% when compared to the random distribution of workload across the same resources.


title: Doubly charmed baryon production in heavy ion collisions

abstract: We give an estimate of $\Xi_{cc}^{++}$ production rate and transverse momentum spectra in relativistic heavy ion collisions. We use Boltzmann transport equations to describe the dynamical evolution of charm quarks and diquarks inside quark-gluon plasma. In-medium formation and dissociation rates of charm diquarks are calculated from potential non-relativistic QCD for the diquark sector. We solve the transport equations by Monte Carlo simulations. For $2.76$ TeV Pb-Pb collisions with $0-10\%$ centrality, the number of $\Xi_{cc}^{++}$ produced in the transverse momentum range $0-5$ GeV and rapidity from $-1$ to $1$ is roughly $0.02$ per collision. We repeat the calculation with a melting temperature $250$ MeV above which no diquarks can be formed. The number of $\Xi_{cc}^{++}$ produced in the same kinematic region is about $0.0125$ per collision. We discuss how to study diquarks at finite temperature on a lattice and construct the anti-triplet free energy in a gauge invariant but path dependent way. We also comment on extensions of the calculation to other doubly heavy baryons and doubly heavy tetraquarks and the feasibility of experimental measurements.


title: The effect of RNA stiffness on the self-assembly of virus particles

abstract: Under many in vitro conditions, some small viruses spontaneously encapsidate a single stranded (ss) RNA into a protein shell called the capsid. While viral RNAs are found to be compact and highly branched because of long distance base-pairing between nucleotides, recent experiments reveal that in a head-to-head competition between a ssRNA with no secondary or higher order structure and a viral RNA, the capsid proteins preferentially encapsulate the linear polymer! In this paper, we study the impact of genome stiffness on the encapsidation free energy of the complex of RNA and capsid proteins. We show that an increase in effective chain stiffness because of base-pairing could be the reason why under certain conditions linear chains have an advantage over branched chains when it comes to encapsidation efficiency. While branching makes the genome more compact, RNA base-pairing increases the effective Kuhn length of the RNA molecule, which could result in an increase of the free energy of RNA confinement, that is, the work required to encapsidate RNA, and thus less efficient packaging.


title: Quiddity sequences for $\mathrm{SL}_3$-frieze patterns

abstract: The notion of a $(k,n)$-frieze pattern was introduced by the author as a generalisation of the classical frieze patterns. In this article we describe connections between classes of $(3,n)$-frieze patterns and classes of $\mathrm{SL}_3$-frieze patterns. We introduce the idea of a superimposed triangulation and clarify how superimposed triangulations may be used to understand quiddity sequences for $\mathrm{SL}_3$-frieze patterns.


title: Piezoelectric Polar Nano Regions and Relaxation-Coupled Resonances in Relaxor Ferroelectrics

abstract: \indent It is a generally accepted fact that the unique dielectric properties of relaxor ferroelectrics are related to the formation of polar nanoregion (PNRs). Less well recognized is the corollary that, because they are polar and therefore lack inversion symmetry, PNRs are also piezoelectric at the nanoscale and can therefore behave as nanoresonators. Using the particular relaxor ferroelectric K$_{\tt1-x}$Li$_{\tt x}$TaO$_{\tt 3}$ (KLT), we show that, when electrically excited into oscillation, these piezoelectric nanoresonators can drive macroscopic electro-mechanical resonances. Unexpectedly however, pairs of coupled resonances are observed, with one of the two exhibiting a characteristic Fano-like lineshape. The complex resonance spectra can be described equally well by two alternative but complementary models both involving two resonances coupled through a relaxation: a purely classical one based on two coupled damped harmonic oscillators and a semi-classical based on two discrete excitations coupled to each other through a continuum. Together, they provide complementary perspectives on the underlying physics of the system. Both reproduce the rapid evolution of the resonance spectrum across three wide temperature ranges, including a phase transition range. In the high temperature range, the coupling between modes is due to the collective $\pi$ relaxation of the lithium ions within PNRs and in the phase transition range to "heterophase relaxation" of the surrounding lattice between its high temperature cubic and low temperature tetragonal phases. The coupling is suppressed in the intermediate range of the collective $\pi/2$ relaxation of the lithium ions. Incidentally, the measured dielectric spectra are shown to bear a surprising but justifiable resemblance to the optical spectra of certain atomic vapors that exhibit electromagnetically induced transparency.


title: The clustering of the SDSS-IV extended Baryon Oscillation Spectroscopic Survey DR14 quasar sample: anisotropic clustering analysis in configuration-space

abstract: We explore the cosmological implications of anisotropic clustering measurements of the quasar sample from Data Release 14 of the Sloan Digital Sky Survey IV Extended Baryon Oscillation Spectroscopic Survey (eBOSS) in configuration space. The $\sim 147,000$ quasar sample observed by eBOSS offers a direct tracer of the density field and bridges the gap of previous BAO measurements between redshift $0.8<z<2.2$. By analysing the two-point correlation function characterized by clustering wedges $\xi_{\rm w_i}(s)$ and multipoles $\xi_{\ell}(s)$, we measure the angular diameter distance, Hubble parameter and cosmic structure growth rate. We define a systematic error budget for our measurements based on the analysis of $N$-body simulations and mock catalogues. Based on the DR14 large scale structure quasar sample at the effective redshift $z_{\rm eff}=1.52$, we find the growth rate of cosmic structure $f\sigma_8(z_{\rm eff})=0.396\pm 0.079$, and the geometric parameters $D_{\rm V}(z)/r_{\rm d}=26.47\pm 1.23$, and $F_{\rm AP}(z)=2.53\pm 0.22$, where the uncertainties include both statistical and systematic errors. These values are in excellent agreement with the best-fitting standard ${\rm \Lambda CDM}$ model to the latest cosmic microwave background data from Planck.


title: What constitutes elemental shape information for biological vision?

abstract: We do not yet understand how the vertebrate visual system provides for recognition of ob- jects. Countless experiments have been performed to examine the contribution of cues such as color, texture, and shadowing, but the most important cues are the con- tours of the outer boundary. Most objects that we can name can be identified as a silhouette, or equally well as a line drawing of the boundary. This has long been ap- preciated, so it is somewhat surprising that after more than a century of experimental research, we have not yet established how our visual system encodes this shape in- formation.


title: 1I/2017 'Oumuamua-like Interstellar Asteroids as Possible Messengers from the Dead Stars

abstract: Discovery of the first interstellar asteroid (ISA) - 1I/2017 'Oumuamua - raised a number of questions regarding its origin. Many of them relate to its lack of cometary activity, suggesting refractory composition of 'Oumuamua. Here we explore the possibility that 'Oumuamua-like ISAs are produced in tidal disruption events (TDEs) of refractory planetoids (asteroids, dwarf planets, etc.) by the white dwarfs (WDs). This idea is supported by existing spectroscopic observations of metal-polluted WDs, hinting at predominantly volatile-poor composition of accreted material. We show that such TDEs sourced by realistic planetary systems (including a population of >1000 km planetoids and massive perturbers - Neptune-to-Saturn mass planets) can eject to interstellar space up to 30% of planetary mass involved in them. Collisional fragmentation, caused by convergent vertical motion of the disrupted planetoid's debris inside the Roche sphere of the WD, channels most of the original mass into 0.1-1 km fragments, similar to 'Oumuamua. Such size spectrum of ISAs (very different from the top-heavy distributions expected in other scenarios) implies that planetary TDEs can account for a significant fraction (up to ~30% under optimistic assumptions) of the ISAs. This figure is based on existing observations of WD metal pollution and accounts for observational biases by using realistic models of circum-WD planetary systems. ISAs should exhibit kinematic characteristics similar to old, dynamically hot Galactic populations; we interpret 'Oumuamua's slow Galactic motion as a statistical fluctuation. ISA ejection in individual planetary TDEs is highly anisotropic, resulting in large fluctuations of their space density. We also show that other ISA production mechanisms involving stellar remnants - direct ejection by massive planets around WDs and SN explosions - have difficulty explaining 'Oumuamua-like ISAs.


title: Superconductor Electronics Fabrication Process with MoN$_x$ Kinetic Inductors and Self-Shunted Josephson Junctions

abstract: Recent progress in superconductor electronics fabrication has enabled single-flux-quantum (SFQ) digital circuits with close to one million Josephson junctions (JJs) on 1-cm$^2$ chips. Increasing the integration scale further is challenging because of the large area of SFQ logic cells, mainly determined by the area of resistively shunted Nb/AlO$_x$-Al/Nb JJs and geometrical inductors utilizing multiple layers of Nb. To overcome these challenges, we are developing a fabrication process with self-shunted high-J$_c$ JJs and compact thin-film MoN$_x$ kinetic inductors instead of geometrical inductors. We present fabrication details and properties of MoN$_x$ films with a wide range of T$_c$, including residual stress, electrical resistivity, critical current, and magnetic field penetration depth {\lambda}$_0$. As kinetic inductors, we implemented Mo$_2$N films with T$_c$ about 8 K, {\lambda}$_0$ about 0.51 {\mu}m, and inductance adjustable in the range from 2 to 8 pH/sq. We also present data on fabrication and electrical characterization of Nb-based self-shunted JJs with AlO$_x$ tunnel barriers and J$_c$ = 0.6 mA/{\mu}m$^2$, and with 10-nm thick Si$_{1-x}$Nb$_x$ barriers, with x from 0.03 to 0.15, fabricated on 200-mm wafers by co-sputtering. We demonstrate that the electron transport mechanism in Si$_{1-x}$Nb$_x$ barriers at x < 0.08 is inelastic resonant tunneling via chains of multiple localized states. At larger x, their Josephson characteristics are strongly dependent on x and residual stress in Nb electrodes, and in general are inferior to AlO$_x$ tunnel barriers.


title: Spitzer Matching survey of the UltraVISTA ultra-deep Stripes (SMUVS): Full-mission IRAC Mosaics and Catalogs

abstract: This paper describes new deep 3.6 and 4.5 micron imaging of three UltraVISTA near-infrared survey stripes within the COSMOS field. The observations were carried out with Spitzer's Infrared Array Camera (IRAC) for the Spitzer Matching Survey of the Ultra-VISTA Deep Stripes (SMUVS). In this work we present our data reduction techniques, and document the resulting mosaics, coverage maps, and catalogs in both IRAC passbands for the three easternmost UltraVISTA survey stripes, covering a combined area of about 0.66 square degrees, of which 0.45 square degrees have at least 20 hr integration time. SMUVS reaches point-source sensitivities of about 25.0 AB mag at both 3.6 and 4.5 microns with a significance of 4-sigma accounting for both survey sensitivity and source confusion. To this limit the SMUVS catalogs contain a total of about 350,000 sources, each of which is detected significantly in at least one IRAC band. Because of its uniform and high sensitivity, relatively large area coverage, and the wide array of ancillary data available in COSMOS, the SMUVS survey will be useful for a large number of cosmological investigations. We will make all images and catalogues described herein publicly available via the Spitzer Science Center.


title: Thermodynamics of a class of regular black holes with a generalized uncertainty principle

abstract: In this article, we present a study on thermodynamics of a class of regular black holes. Such a class includes Bardeen and Hayward regular black holes. We obtained thermodynamic quantities like the Hawking temperature, entropy, and heat capacity for the entire class. As part of an effort to indicate some physical observable to distinguish regular black holes from singular black holes, we suggest that regular black holes are colder than singular black holes. Besides, contrary to the Schwarzschild black hole, that class of regular black holes may be thermodynamically stable. From a generalized uncertainty principle, we also obtained the quantum-corrected thermodynamics for the studied class. Such quantum corrections provide a logarithmic term for the quantum-corrected entropy.


title: Tensor network ranks

abstract: In problems involving approximation, completion, denoising, dimension reduction, estimation, interpolation, modeling, order reduction, regression, etc, we argue that the near-universal practice of assuming that a function, matrix, or tensor (which we will see are all the same object in this context) has \emph{low rank} may be ill-justified. There are many natural instances where the object in question has high rank with respect to the classical notions of rank: matrix rank, tensor rank, multilinear rank --- the latter two being the most straightforward generalizations of the former. To remedy this, we show that one may vastly expand these classical notions of ranks: Given any undirected graph $G$, there is a notion of $G$-rank associated with $G$, which provides us with as many different kinds of ranks as there are undirected graphs. In particular, the popular tensor network states in physics (e.g., \textsc{mps}, \textsc{ttns}, \textsc{peps}) may be regarded as functions of a specific $G$-rank for various choices of $G$. Among other things, we will see that a function, matrix, or tensor may have very high matrix, tensor, or multilinear rank and yet very low $G$-rank for some $G$. In fact the difference is in the orders of magnitudes and the gaps between $G$-ranks and these classical ranks are arbitrarily large for some important objects in computer science, mathematics, and physics. Furthermore, we show that there is a $G$ such that almost every tensor has $G$-rank exponentially lower than its rank or the dimension of its ambient space.


title: Magnetoresistance when Spin Effects on Conduction are Weak

abstract: This paper considers certain materials, including topological insulators, where spin rotation symmetry is broken much more strongly than time reversal symmetry. When these materials are in the diffusive regime, i.e. when they have disorder that is strong enough to cause an electron to scatter many times while crossing a sample, electrons and holes move in pairs that have zero spin and are insensitive to spin physics. Working within this spinless scenario, we show that Fourier transforming the magnetoconductance with respect to external magnetic field obtains a curve describing the area distribution of loops traced by electrons and holes within the sample. We present loop area distributions of Landau levels, weak (anti)localization, conduction governed by Levy flights, and linear-in-field resistance. Of these four the last two are new results. Comparing these distributions, we argue that the linear-in-field resistance seen in some topological insulators is caused by the same diffusive scattering that causes weak antilocalization. The difference is that linear-in-field resistance materials retain a level of quantum coherence that is usually seen only on the surface of 2-D wires or in ring geometries. In an appendix we include some speculative material about linear-in-temperature resistance.


title: On Division Polynomial PIT and Supersingularity

abstract: For an elliptic curve $E$ over a finite field $\F_q$, where $q$ is a prime power, we propose new algorithms for testing the supersingularity of $E$. Our algorithms are based on the Polynomial Identity Testing (PIT) problem for the $p$-th division polynomial of $E$. In particular, an efficient algorithm using points of high order on $E$ is given.


title: Symbolic relative entropy in quantifying nonlinear dynamics of equalities-involved heartbeats

abstract: Symbolic relative entropy, an efficient nonlinear complexity parameter measuring probabilistic divergences of symbolic sequences, is proposed in our nonlinear dynamics analysis of heart rates considering equal states. Equalities are not rare in discrete heartbeats because of the limits of resolution of signals collection, and more importantly equal states contain underlying important cardiac regulation information which is neglected by some chaotic deterministic parameters and temporal asymmetric measurements. The relative entropy of symbolization associated with equal states has satisfied nonlinear dynamics complexity detections in heartbeats and shows advantages to some nonlinear dynamics parameters without considering equalities. Researches on cardiac activities suggest the highest probabilistic divergence of the healthy young heart rates and highlight the facts that heart diseases and aging reduce the nonlinear dynamical complexity of heart rates.


title: Increased sensitivity of spin noise spectroscopy using homodyne detection in $n$-doped GaAs

abstract: We implement the homodyne detection scheme for an increase of the polarimetric sensitivity in spin noise spectroscopy. Controlling the laser intensity of the local oscillator, which is guided around the sample and does not perturb the measured spin system, we are able to improve the signal-to-noise ratio. The opportunity of additional amplification of the measured signal strength allows us to reduce the probe laser intensity incident onto the sample and therefore to approach the non-perturbative regime. The efficiency of this scheme with signal enhancement by more than a factor of 3 at low probe powers is demonstrated on bulk $n$-doped GaAs where the reduced electron-spin relaxation rate is shown experimentally. Additionally, the control of the optical phase provides us with the possibility to switch between the measurement of Faraday rotation and ellipticity without changes in the optical setup.


title: Probing the single-particle character of rotational states in $^{19}$F using a short-lived isomeric beam

abstract: A beam containing a substantial component of both the $J^{\pi}=5^+$, $T_{1/2}=162$ ns isomeric state of $^{18}$F and its $1^+$, 109.77-min ground state has been utilized to study members of the ground-state rotational band in $^{19}$F through the neutron transfer reaction $(d$,$p)$ in inverse kinematics. The resulting spectroscopic strengths confirm the single-particle nature of the 13/2$^+$ band-terminating state. The agreement between shell-model calculations, using an interaction constructed within the $sd$ shell, and our experimental results reinforces the idea of a single-particle/collective duality in the descriptions of the structure of atomic nuclei.


title: Evorus: A Crowd-powered Conversational Assistant Built to Automate Itself Over Time

abstract: Crowd-powered conversational assistants have been shown to be more robust than automated systems, but do so at the cost of higher response latency and monetary costs. A promising direction is to combine the two approaches for high quality, low latency, and low cost solutions. In this paper, we introduce Evorus, a crowd-powered conversational assistant built to automate itself over time by (i) allowing new chatbots to be easily integrated to automate more scenarios, (ii) reusing prior crowd answers, and (iii) learning to automatically approve response candidates. Our 5-month-long deployment with 80 participants and 281 conversations shows that Evorus can automate itself without compromising conversation quality. Crowd-AI architectures have long been proposed as a way to reduce cost and latency for crowd-powered systems; Evorus demonstrates how automation can be introduced successfully in a deployed system. Its architecture allows future researchers to make further innovation on the underlying automated components in the context of a deployed open domain dialog system.


title: The optical afterglow of the short gamma-ray burst associated with GW170817

abstract: The binary neutron star merger GW170817 was the first multi-messenger event observed in both gravitational and electromagnetic waves. The electromagnetic signal began approximately 2 seconds post-merger with a weak, short burst of gamma-rays, which was followed over the next hours and days by the ultraviolet, optical and near-infrared emission from a radioactively- powered kilonova. Later, non-thermal rising X-ray and radio emission was observed. The low luminosity of the gamma-rays and the rising non-thermal flux from the source at late times could indicate that we are outside the opening angle of the beamed relativistic jet. Alternatively, the emission could be arising from a cocoon of material formed from the interaction between a jet and the merger ejecta. Here we present late-time optical detections and deep near-infrared limits on the emission from GW170817 at 110 days post-merger. Our new observations are at odds with expectations of late-time emission from kilonova models, being too bright and blue. Instead, the emission arises from the interaction between the relativistic ejecta of GW170817 and the interstellar medium. We show that this emission matches the expectations of a Gaussian structured relativistic jet, which would have launched a high luminosity short GRB to an aligned observer. However, other jet structure or cocoon models can also match current data - the future evolution of the afterglow will directly distinguish the origin of the emission.


title: Electroweak phase transitions in multi-Higgs models: the case of Trinification-inspired THDSM

abstract: The rich vacuum structure of multi-Higgs extensions of the Standard Model (SM) may have interesting cosmological implications for the electroweak phase transition (EWPT). As an important example of such class of models, we consider a particularly simple low-energy SM-like limit of a recently proposed Grand-Unified Trinification model with the scalar sector composed of two Higgs doublets and a complex singlet and with a global $\mathrm{U}(1)$ family symmetry. The fermion sector of this model is extended with a family of vector-like quarks which enhances CP violation. With the current study, we aim at exploring the generic vacuum structure and uncovering the features of the EWPT in this model relevant for cosmology. We show the existence of different phase transition patterns providing strong departure from thermal equilibrium. Most of these observations are not specific to the considered model and may generically be expected in other multi-Higgs extensions of the SM.


title: Absorber Model: the Halo-like model for the Lyman-alpha forest

abstract: We present a semi-analytic model for the Lyman-$\alpha$ forest that is inspired by the Halo Model. This model is built on the absorption line decomposition of the forest. Flux correlations are decomposed into those within each absorption line (the 1-absorber term) and those between separate lines (the 2-absorber term), treating the lines as biased tracers of the underlying matter fluctuations. While the nonlinear exponential mapping between optical depth and flux requires an infinite series of moments to calculate any statistic, we show that this series can be resumed (truncating at the desired order in the linear matter overdensity). We focus on the $z=2-3$ line-of-sight power spectrum. Our model finds that 1-absorber term dominates the power on all scales, with most of its contribution coming from HI columns of $10^{14}-10^{15}\;\mathrm{cm^{-2}}$, while the smaller 2-absorber contribution comes from lower columns that trace overdensities of a few. The prominence of the 1-absorber correlations indicates that the line-of-sight power spectrum is shaped principally by the lines' number densities and their absorption profiles, with correlations between lines contributing to a lesser extent. We present intuitive formulae for the effective optical depth as well as the large-scale limits of 1-absorber and 2-absorber terms, which simplify to integrals over the HI column density distribution with different equivalent-width weightings. With minimalist models for the bias of absorption systems and their peculiar velocity broadening, our model predicts values for the density bias and velocity gradient bias that are consistent with those found in simulations.


title: Violable Contracts and Governance for Blockchain Applications

abstract: We examine blockchain technologies, especially smart contracts, as a platform for decentralized applications. By providing a basis for consensus, blockchain promises to upend business models that presuppose a central authority. However, blockchain suffers from major shortcomings arising from an over-regimented way of organizing computation that limits its prospects. We propose a sociotechnical, yet computational, perspective that avoids those shortcomings. A centerpiece of our vision is the notion of a declarative, violable contract in contradistinction to smart contracts. This new way of thinking enables flexible governance, by formalizing organizational structures; verification of correctness without obstructing autonomy; and a meaningful basis for trust.


title: The Colorado Ultraviolet Transit Experiment (CUTE): A dedicated cubesat mission to study exoplanetary mass loss and magnetic fields

abstract: The Colorado Ultraviolet Transit Experiment (CUTE) is a near-UV (2550 - 3300 Angstrom) 6U cubesat mission designed to monitor transiting hot Jupiters to quantify their atmospheric mass loss and magnetic fields. CUTE will probe both atomic (Mg and Fe) and molecular (OH) lines for evidence of enhanced transit absorption, and to search for evidence of early ingress due to bow shocks ahead of the planet's orbital motion. As a dedicated mission, CUTE will observe more than 100 spectroscopic transits of hot Jupiters over a nominal seven month mission. This represents the equivalent of more than 700 orbits of the only other instrument capable of these measurements, the Hubble Space Telescope. CUTE efficiently utilizes the available cubesat volume by means of an innovative optical design to achieve a projected effective area of 28 sq. cm, low instrumental background, and a spectral resolving power of 3000 over the primary science bandpass. These performance characteristics enable CUTE to discern transit depths between 0.1 - 1% in individual spectral absorption lines. We present the CUTE optical and mechanical design, a summary of the science motivation and expected results, and an overview of the projected fabrication, calibration and launch timeline.


title: Improvements in Gravitational-Wave Sky Localization with Expanded Networks of Interferometers

abstract: A milestone of multi-messenger astronomy has been achieved with the detection of gravitational waves from a binary neutron star merger accompanied by observations of several associated electromagnetic counterparts. Joint observations can reveal details of the engines that drive the electromagnetic and gravitational-wave emission. However, locating and identify an electromagnetic counterparts to a gravitational-wave event is heavily reliant on localization of the source through gravitational-wave information. We explore the sky localization of a simulated set of neutron star mergers as the worldwide network of gravitational-wave detectors evolves through the next decade, performing the first such study for neutron star -- black hole binary sources. Currently, three detectors are observing with additional detectors in Japan and India expected to become operational in the coming years. With three detectors, we recover a median neutron star -- black hole binary sky localization of 60 deg$^2$ at the 90\% credible level. As all five detectors become operational, sources can be localized to a median of 11 deg$^2$ on the sky.


title: Wulff shapes and a characterization of simplices via a Bezout type inequality

abstract: Inspired by a fundamental theorem of Bernstein, Kushnirenko, and Khovanskii we study the following Bezout type inequality for mixed volumes $$ V(L_1,\dots,L_{n})V_n(K)\leq V(L_1,K[{n-1}])V(L_2,\dots, L_{n},K). $$ We show that the above inequality characterizes simplices, i.e. if $K$ is a convex body satisfying the inequality for all convex bodies $L_1, \dots, L_n \subset {\mathbb R}^n$, then $K$ must be an $n$-dimensional simplex. The main idea of the proof is to study perturbations given by Wulff shapes. In particular, we prove a new theorem on differentiability of the support function of the Wulff shape, which is of independent interest. In addition, we study the Bezout inequality for mixed volumes introduced in arXiv:1507.00765 . We introduce the class of weakly decomposable convex bodies which is strictly larger than the set of all polytopes that are non-simplices. We show that the Bezout inequality in arXiv:1507.00765 characterizes weakly indecomposable convex bodies.


title: Strong- and weak-universal critical behaviour of a mixed-spin Ising model with triplet interactions on the Union Jack (centered square) lattice

abstract: The mixed spin-1/2 and spin-S Ising model on the Union Jack (centered square) lattice with four different three-spin (triplet) interactions and the uniaxial single-ion anisotropy is exactly solved by establishing a rigorous mapping equivalence with the corresponding zero-field (symmetric) eight-vertex model on a dual square lattice. A rigorous proof of the aforementioned exact mapping equivalence is provided by two independent approaches exploiting either a graph-theoretical or spin representation of the zero-field eight-vertex model. An influence of the interaction anisotropy as well as the uniaxial single-ion anisotropy on phase transitions and critical phenomena is examined in particular. It is shown that the considered model exhibits a strong-universal critical behaviour with constant critical exponents when considering the isotropic model with four equal triplet interactions or the anisotropic model with one triplet interaction differing from the other three. The anisotropic models with two different triplet interactions, which are pairwise equal to each other, contrarily exhibit a weak-universal critical behaviour with critical exponents continuously varying with a relative strength of the triplet interactions as well as the uniaxial single-ion anisotropy. It is evidenced that the variations of critical exponents of the mixed-spin Ising models with the integer-valued spins S differ basically from their counterparts with the half-odd-integer spins S.


title: Geometrical picture of photocounting measurements

abstract: We revisit the representation of generalized quantum observables by establishing a geometric picture in terms of their positive operator-valued measures (POVMs). This leads to a clear geometric interpretation of Born's rule by introducing the concept of contravariant operator-valued measures. Our approach is applied to the theory of array detectors, which is a challenging task as the finite dimensionality of the POVM substantially restricts the available information about quantum states. Our geometric technique allows for a direct estimation of expectation values of different observables, which are typically not accessible with such detection schemes. In addition, we also demonstrate the applicability of our method to quantum-state reconstruction with unbalanced homodyne detection.


title: A Unified Enumeration of 1-dimension Garden Algebras and Valise Adinkras

abstract: In the study of supersymmetry in one dimension, various works enumerate sets of generators of garden algebras $GR(d,N)$ (and equivalently, valise Adinkras) for special cases $N = d = 4$ and $N = d = 8$, using group-theoretic methods and computer computation. We complement this work by enumerating the objects for arbitrary $N$ and $d$ via a formula in a streamlined manner.


title: Graph-Based Radio Resource Management for Vehicular Networks

abstract: This paper investigates the resource allocation problem in device-to-device (D2D)-based vehicular communications, based on slow fading statistics of channel state information (CSI), to alleviate signaling overhead for reporting rapidly varying accurate CSI of mobile links. We consider the case when each vehicle-to-infrastructure (V2I) link shares spectrum with multiple vehicle-to-vehicle (V2V) links. Leveraging the slow fading statistical CSI of mobile links, we maximize the sum V2I capacity while guaranteeing the reliability of all V2V links. We propose a graph-based algorithm that uses graph partitioning tools to divide highly interfering V2V links into different clusters before formulating the spectrum sharing problem as a weighted 3-dimensional matching problem, which is then solved through adapting a high-performance approximation algorithm.


title: Wigner function of noninteracting trapped fermions

abstract: We study analytically the Wigner function $W_N({\bf x},{\bf p})$ of $N$ noninteracting fermions trapped in a smooth confining potential $V({\bf x})$ in $d$ dimensions. At zero temperature, $W_N({\bf x},{\bf p})$ is constant over a finite support in the phase space $({\bf x},{\bf p})$ and vanishes outside. Near the edge of this support, we find a universal scaling behavior of $W_N({\bf x},{\bf p})$ for large $N$. The associated scaling function is independent of the precise shape of the potential as well as the spatial dimension $d$. We further generalize our results to finite temperature $T>0$. We show that there exists a low temperature regime $T \sim e_N/b$ where $e_N$ is an energy scale that depends on $N$ and the confining potential $V({\bf x})$, where the Wigner function at the edge again takes a universal scaling form with a $b$-dependent scaling function. This temperature dependent scaling function is also independent of the potential as well as the dimension $d$. Our results generalize to any $d\geq 1$ and $T \geq 0$ the $d=1$ and $T=0$ results obtained by Bettelheim and Wiegman [Phys. Rev. B ${\bf 84}$, 085102 (2011)].

categories: cond-mat.stat-mech cond-mat.quant-gas math-ph math.MP


title: Diversification, economies of scope, and exports growth of Chinese firms

abstract: In the 1990s, China started a process of structural reforms and of trade liberalization, which was followed by the accession to the World Trade Organization (WTO) in 2001. In this paper, we analyze trade patterns of Chinese firms for the period 2000-2006, characterized by a notable increase in exports volumes. Theoretically, in a more open economy, firms are expected to move from the production of a set of less-competitive products towards more internationally competitive ones, which implies specialization. We study several stylized facts on the distribution of Chinese firms trade and growth rates, and we analyze whether firms have diversified or specialized their trade patterns between 2000 and 2006. We show that Chinese export patterns are very heterogeneous, that the volatility of growth rates depends on the level of exports, and that volatility is stronger after trade liberalization. Both, diversification in products and destinations have a positive impact on trade growth, but diversification of destinations has a stronger effect. We conclude that the success of Chinese exports is not only due to an increase in the intensive margin, related to the existence of economies of scale, but also due to an increase in the extensive margin, related to the existence of economies of scope.


title: Topology and homotopy of lattice isomorphic arrangements

abstract: We prove the existence of lattice isomorphic line arrangements having $\pi_1$-equivalent or homotopy-equivalent complements and non homeomorphic embeddings in the complex projective plane. We also provide two explicit examples, one is formed by real-complexified arrangements while the second is not.


title: On the definition of Suzuki groups over rings

abstract: The definition of Suzuki groups over rings is given by means of an explicit description as a difference-algebraic group. For a (not necessarily perfect) field with more than two elements this construction produces a simple group.


title: Generative Sensing: Transforming Unreliable Sensor Data for Reliable Recognition

abstract: This paper introduces a deep learning enabled generative sensing framework which integrates low-end sensors with computational intelligence to attain a high recognition accuracy on par with that attained with high-end sensors. The proposed generative sensing framework aims at transforming low-end, low-quality sensor data into higher quality sensor data in terms of achieved classification accuracy. The low-end data can be transformed into higher quality data of the same modality or into data of another modality. Different from existing methods for image generation, the proposed framework is based on discriminative models and targets to maximize the recognition accuracy rather than a similarity measure. This is achieved through the introduction of selective feature regeneration in a deep neural network (DNN). The proposed generative sensing will essentially transform low-quality sensor data into high-quality information for robust perception. Results are presented to illustrate the performance of the proposed framework.


title: P-MOD: Secure Privilege-Based Multilevel Organizational Data-Sharing in Cloud Computing

abstract: Cloud computing has changed the way enterprises store, access and share data. Data is constantly being uploaded to the cloud and shared within an organization built on a hierarchy of many different individuals that are given certain data access privileges. With more data storage needs turning over to the cloud, finding a secure and efficient data access structure has become a major research issue. With different access privileges, individuals with more privileges (at higher levels of the hierarchy) are granted access to more sensitive data than those with fewer privileges (at lower levels of the hierarchy). In this paper, a Privilege-based Multilevel Organizational Data-sharing scheme~(P-MOD) is proposed that incorporates a privilege-based access structure into an attribute-based encryption mechanism to handle these concerns. Each level of the privilege-based access structure is affiliated with an access policy that is uniquely defined by specific attributes. Data is then encrypted under each access policy at every level to grant access to specific data users based on their data access privileges. An individual ranked at a certain level can decrypt the ciphertext (at that specific level) if and only if that individual owns a correct set of attributes that can satisfy the access policy of that level. The user may also decrypt the ciphertexts at the lower levels with respect to the user's level. Security analysis shows that P-MOD is secure against adaptively chosen plaintext attack assuming the DBDH assumption holds.The comprehensive performance analysis demonstrates that P-MOD is more efficient in computational complexity and storage space than the existing schemes in secure data sharing within an organization.


title: Towards Multi-Object Detection and Tracking in Urban Scenario under Uncertainties

abstract: Urban-oriented autonomous vehicles require a reliable perception technology to tackle the high amount of uncertainties. The recently introduced compact 3D LIDAR sensor offers a surround spatial information that can be exploited to enhance the vehicle perception. We present a real-time integrated framework of multi-target object detection and tracking using 3D LIDAR geared toward urban use. Our approach combines sensor occlusion-aware detection method with computationally efficient heuristics rule-based filtering and adaptive probabilistic tracking to handle uncertainties arising from sensing limitation of 3D LIDAR and complexity of the target object movement. The evaluation results using real-world pre-recorded 3D LIDAR data and comparison with state-of-the-art works shows that our framework is capable of achieving promising tracking performance in the urban situation.


title: Term Relevance Feedback for Contextual Named Entity Retrieval

abstract: We address the role of a user in Contextual Named Entity Retrieval (CNER), showing (1) that user identification of important context-bearing terms is superior to automated approaches, and (2) that further gains are possible if the user indicates the relative importance of those terms. CNER is similar in spirit to List Question answering and Entity disambiguation. However, the main focus of CNER is to obtain user feedback for constructing a profile for a class of entities on the fly and use that to retrieve entities from free text. Given a sentence, and an entity selected from that sentence, CNER aims to retrieve sentences that have entities similar to query entity. This paper explores obtaining term relevance feedback and importance weighting from humans in order to improve a CNER system. We report our findings based on the efforts of IR researchers as well as crowdsourced workers.


title: Joint Resource Allocation and Antenna Selection In the Uplink of OFDMA Networks

abstract: The problem of joint sub-channel, power control and antenna selection in the uplink of OFDMA networks is investigated. The corresponding optimization problem is a non convex mixed integer non-linear problem(MINLP). To tackle the problem, an iterative two-step optimization problem is devised. First,~Having fixed the selected antenna as well as the associated power of each user,~a sub-channel assignment problem using the so-called Hungarian method is proposed. Second, the joint power and antenna selection problem is extended into a more elegant form through adding a penalty term to the objective function forcing the relaxed boolean variables associated with antenna selection problem to take binary values. To this end, the problem is restated as a difference of two concave function (DC programming) and is effectively solved in an iterative manner. The foregoing procedure is iteratively repeated until approaching a fixed network throughput. Simulation results demonstrate the superiority of the proposed method w.r.t. the exhaustive search with fixed power control mechanism.


title: The clustering of the SDSS-IV extended Baryon Oscillation Spectroscopic Survey DR14 quasar sample: structure growth rate measurement from the anisotropic quasar power spectrum in the redshift range $0.8<z<2.2$

abstract: We analyse the clustering of the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey Data Release 14 quasar sample (DR14Q). We measure the redshift space distortions using the power spectrum monopole, quadrupole and hexadecapole inferred from 148,659 quasars between redshifts 0.8 and 2.2 covering a total sky footprint of 2112.9 deg$^2$. We constrain the logarithmic growth of structure times the amplitude of dark matter density fluctuations, $f\sigma_8$, and the Alcock-Paczynski dilation scales which allow constraints to be placed on the angular diameter distance $D_A(z)$ and the Hubble $H(z)$ parameter. At the effective redshift of $z_{\rm eff}=1.52$, $f\sigma_8(z_{\rm eff})=0.420\pm0.076$, $H(z_{\rm eff})=[162\pm 12]\, (r_s^{\rm fid}/r_s)\,{\rm km\, s}^{-1}{\rm Mpc}^{-1}$, and $D_A(z_{\rm eff})=[1.85\pm 0.11]\times10^3\,(r_s/r_s^{\rm fid})\,{\rm Mpc}$, where $r_s$ is the comoving sound horizon at the baryon drag epoch and the superscript `fid' stands for its fiducial value. The errors take into account the full error budget, including systematics and statistical contributions. These results are in full agreement with the current $\Lambda$-Cold Dark Matter ($\Lambda$CDM) cosmological model inferred from Planck measurements. Finally, we compare our measurements with other eBOSS companion papers and find excellent agreement, demonstrating the consistency and complementarity of the different methods used for analysing the data.


title: DCASE 2017 Task 1: Acoustic Scene Classification Using Shift-Invariant Kernels and Random Features

abstract: Acoustic scene recordings are represented by different types of handcrafted or Neural Network-derived features. These features, typically of thousands of dimensions, are classified in state of the art approaches using kernel machines, such as the Support Vector Machines (SVM). However, the complexity of training these methods increases with the dimensionality of these input features and the size of the dataset. A solution is to map the input features to a randomized lower-dimensional feature space. The resulting random features can approximate non-linear kernels with faster linear kernel computation. In this work, we computed a set of 6,553 input features and used them to compute random features to approximate three types of kernels, Gaussian, Laplacian and Cauchy. We compared their performance using an SVM in the context of the DCASE Task 1 - Acoustic Scene Classification. Experiments show that both, input and random features outperformed the DCASE baseline by an absolute 4%. Moreover, the random features reduced the dimensionality of the input by more than three times with minimal loss of performance and by more than six times and still outperformed the baseline. Hence, random features could be employed by state of the art approaches to compute low-storage features and perform faster kernel computations.


title: A Trip to the Moon: Personalized Animated Movies for Self-reflection

abstract: Self-tracking physiological and psychological data poses the challenge of presentation and interpretation. Insightful narratives for self-tracking data can motivate the user towards constructive self-reflection. One powerful form of narrative that engages audience across various culture and age groups is animated movies. We collected a week of self-reported mood and behavior data from each user and created in Unity a personalized animation based on their data. We evaluated the impact of their video in a randomized control trial with a non-personalized animated video as control. We found that personalized videos tend to be more emotionally engaging, encouraging greater and lengthier writing that indicated self-reflection about moods and behaviors, compared to non-personalized control videos.


title: L-spaces, left-orderability and two-bridge knots

abstract: We show that the 3-fold cyclic branched cover of any genus 2 two-bridge knot $K_{[-2q,2s,-2t,2l]}$ is an L-space and its fundamental group is not left-orderable. Therefore the family of 3-fold cyclic branched cover of any genus 2 two-bridge knot $K_{[-2q,2s,-2t,2l]}$ verifies the $L$-space conjecture. We also show that if $K_{[2k,-2l]}$ is a 2-bridge knot with $k\geq 2$, $l>0$, then the fundamental group of the 5-fold cyclic branched cover of $K_{[2k,-2l]}$ is not left-orderable, which will complete the proof that the fundamental group of the 5-fold cyclic branched cover of any genus one two-bridge knot is not left-orderable.


title: Stable Marriage with Multi-Modal Preferences

abstract: We introduce a generalized version of the famous Stable Marriage problem, now based on multi-modal preference lists. The central twist herein is to allow each agent to rank its potentially matching counterparts based on more than one "evaluation mode" (e.g., more than one criterion); thus, each agent is equipped with multiple preference lists, each ranking the counterparts in a possibly different way. We introduce and study three natural concepts of stability, investigate their mutual relations and focus on computational complexity aspects with respect to computing stable matchings in these new scenarios. Mostly encountering computational hardness (NP-hardness), we can also spot few islands of tractability and make a surprising connection to the \textsc{Graph Isomorphism} problem.


title: Microwave-optical two-photon excitation of Rydberg states

abstract: We report efficient microwave-optical two photon excitation of Rb Rydberg atoms in a magneto optical trap. This approach allows the excitation of normally inaccessible states and provides a path toward excitation of high angular momentum states. The efficiency stems from the elimination of the Doppler width, the use of a narrow band pulsed laser, and the enormous electric dipole matrix element connecting the intermediate and final states of the transition. The excitation is efficient in spite of the low optical and microwave powers, of order 1 kW and 1 mW, respectively. This is an application of the large dipole coupling strengths between Rydberg states to achieve two photon excitation of Rydberg atoms.


title: Traveling salesman problem across dense cities

abstract: Consider~\(n\) nodes~\(\{X_i\}_{1 \leq i \leq n}\) distributed independently across~\(N\) cities contained with the unit square~\(S\) according to a distribution~\(f.\) Each city is modelled as an~\(r_n \times r_n\) square contained within~\(S\) and let~\(TSPC_n\) denote the length of the minimum length cycle containing all the~\(n\) nodes, corresponding to the traveling salesman problem (TSP). We obtain variance estimates for~\(TSPC_n\) and prove that if the cities are well-connected and densely populated in a certain sense, then~\(TSPC_n\) appropriately centred and scaled converges to zero in probability. We also obtain large deviation type estimates for~\(TSPC_n.\) Using the proof techniques, we alternately obtain corresponding results for the length~\(TSP_n\) of the minimum length cycle in the unconstrained case, when the nodes are independently distributed throughout the unit square~\(S.\)


title: Exact moments of the Sachdev-Ye-Kitaev model up to order $1/N^2$

abstract: We analytically evaluate the moments of the spectral density of the $q$-body Sachdev-Ye-Kitaev (SYK) model, and obtain order $1/N^2$ corrections for all moments, where $N$ is the total number of Majorana fermions. To order $1/N$, moments are given by those of the weight function of the Q-Hermite polynomials. Representing Wick contractions by rooted chord diagrams, we show that the $1/N^2$ correction for each chord diagram is proportional to the number of triangular loops of the corresponding intersection graph, with an extra grading factor when $q$ is odd. Therefore the problem of finding $1/N^2$ corrections is mapped to a triangle counting problem. Since the total number of triangles is a purely graph-theoretic property, we can compute them for the $q=1$ and $q=2$ SYK models, where the exact moments can be obtained analytically using other methods, and therefore we have solved the moment problem for any $q$ to $1/N^2$ accuracy. The moments are then used to obtain the spectral density of the SYK model to order $1/N^2$. We also obtain an exact analytical result for all contraction diagrams contributing to the moments, which can be evaluated up to eighth order. This shows that the Q-Hermite approximation is accurate even for small values of $N$.


title: Minimum spanning trees across dense cities

abstract: Consider~\(n\) nodes distributed independently across~\(N\) cities contained with the unit square~\(S\) according to a distribution~\(f.\) Each city is modelled as an~\(r_n \times r_n\) square contained within~\(S\) and~\(MSTC_n\) denotes the length of the minimum spanning tree containing all the~\(n\) nodes. We use approximation methods to obtain variance estimates for~\(MSTC_n\) and prove that if the cities are well-connected and densely populated in a certain sense, then~\(MSTC_n\) appropriately centred and scaled converges to zero in probability. Using the proof techniques, we alternately derive corresponding results for the length~\(MST_n\) of the minimum spanning tree for the usual case when the nodes are independently distributed throughout the unit square~\(S.\) In particular, we obtain that the variance of~\(MST_n\) grows at most as a power of the logarithm of~\(n\) and use a subsequence argument to get almost sure convergence of~\(MST_n\) appropriately centred and scaled.


title: Seismology with optical links: enabling a global network for submarine earthquake monitoring

abstract: Earthquake monitoring across the globe is currently achieved with networks of seismic stations. The data from these networks have been instrumental in advancing our understanding of the Earth's interior structure and dynamic behaviour. However, almost all seismic stations are located on land and earthquakes of magnitude smaller than 4 at the bottom of the oceans remain largely undetected. Here we show that ordinary telecommunication optical fibre links can detect seismic events when combined with state-of-the-art frequency metrology techniques. We have detected earthquakes over terrestrial and submarine optical fibre links with length ranging from 75 to 535 km and a geographical distance from the earthquake's epicentre ranging from 25 to 18,500 km. In contrast to existing commercial reflectometry-based acoustic sensing methods used widely in the oil and gas industry, which are limited to only a few tens of kilometres, the technique presented here can be extended over thousands of kilometres, paving the way for detection of remote underwater earthquakes. By using the proposed technique on the existing extensive submarine optical fibre infrastructure, which already criss-crosses the seas and oceans, a global seismic network for real-time detection of underwater earthquakes could be implemented. The ability to detect off-shore earthquakes closer to the source could also enable a cost-effective solution for early detection of tsunamis.


title: Measuring the binary black hole mass spectrum with an astrophysically motivated parameterization

abstract: Gravitational-wave detections have revealed a previously unknown population of stellar mass black holes with masses above $20\, M_{\odot}$. These observations provide a new way to test models of stellar evolution for massive stars. By considering the astrophysical processes likely to determine the shape of the binary black hole mass spectrum, we construct a parameterized model to capture key spectral features that relate gravitational-wave data to theoretical stellar astrophysics. In particular, we model the signature of pulsational pair-instability supernovae, which are expected to cause all stars with initial mass $100\, M_{\odot}\lesssim M \lesssim 150\, M_{\odot}$ to form $\sim 40\, M_{\odot}$ black holes. This would cause a cut-off in the black hole mass spectrum along with an excess of black holes near $40\, M_{\odot}$. We carry out a simulated data study to illustrate some of the stellar physics that can be inferred using gravitational-wave measurements of binary black holes and demonstrate several such inferences that might be made in the near future. First, we measure the minimum and maximum stellar black hole mass. Second, we infer the presence of a peak due to pair-instability supernovae. Third, we measure the black hole mass ratio distribution. Finally, we show how inadequate models of the black hole mass spectrum lead to biased estimates of the merger rate and the amplitude of the stochastic gravitational-wave background.


title: Mass-structure of weighted real trees

abstract: Rooted, weighted continuum random trees are used to describe limits of sequences of random discrete trees. Formally, they are random quadruples $(\mathcal{T},d,r,p)$, where $(\mathcal{T},d)$ is a tree-like metric space, $r\in\mathcal{T}$ is a distinguished root, and $p$ is a probability measure on this space. The underlying branching structure is carried implicitly in the metric $d$. We explore various ways of describing the interaction between branching structure and mass in $(\mathcal{T},d,r,p)$ in a way that depends on $d$ only by way of this branching structure. We introduce a notion of mass-structure equivalence and show that two rooted, weighted $\mathbb{R}$-trees are equivalent in this sense if and only if the discrete hierarchies derived by i.i.d. sampling from their weights, in a manner analogous to Kingman's paintbox, have the same distribution. We introduce a family of trees, called "interval partition trees" that serve as representatives of mass-structure equivalence classes, and which naturally represent the laws of the aforementioned hierarchies.


title: Novel Impossibility Results for Group-Testing

abstract: In this work we prove non-trivial impossibility results for perhaps the simplest non-linear estimation problem, that of {\it Group Testing} (GT), via the recently developed Madiman-Tetali inequalities. Group Testing concerns itself with identifying a hidden set of $d$ defective items from a set of $n$ items via $t$ {disjunctive/pooled} measurements ("group tests"). We consider the linear sparsity regime, i.e. $d = \delta n$ for any constant $\delta >0$, a hitherto little-explored (though natural) regime. In a standard information-theoretic setting, where the tests are required to be non-adaptive and a small probability of reconstruction error is allowed, our lower bounds on $t$ are the {\it first} that improve over the classical counting lower bound, $t/n \geq H(\delta)$, where $H(\cdot)$ is the binary entropy function. As corollaries of our result, we show that (i) for $\delta \gtrsim 0.347$, individual testing is essentially optimal, i.e., $t \geq n(1-o(1))$; and (ii) there is an {adaptivity gap}, since for $\delta \in (0.3471,0.3819)$ known {adaptive} GT algorithms require fewer than $n$ tests to reconstruct ${\cal D}$, whereas our bounds imply that the best nonadaptive algorithm must essentially be individual testing of each element. Perhaps most importantly, our work provides a framework for combining combinatorial and information-theoretic methods for deriving non-trivial lower bounds for a variety of non-linear estimation problems.


title: Revealed Price Preference: Theory and Empirical Analysis

abstract: To determine the welfare implications of price changes in demand data, we introduce a revealed preference relation over prices. We show that the absence of cycles in this relation characterizes a consumer who trades off the utility of consumption against the disutility of expenditure. Our model can be applied whenever a consumer's demand over a strict subset of all available goods is being analyzed; it can also be extended to settings with discrete goods and nonlinear prices. To illustrate its use, we apply our model to a single-agent data set and to a data set with repeated cross-sections. We develop a novel test of linear hypotheses on partially identified parameters to estimate the proportion of the population who are revealed better off due to a price change in the latter application. This new technique can be used for nonparametric counterfactual analysis more broadly.