title: Axion-photon mixing in quantum field theory and vacuum energy

abstract: We analyze axion--photon mixing in the framework of quantum field theory. The condensate structure of the vacuum for mixed fields induces corrections to the oscillation formulae and leads to non-zero energy of the vacuum for the component of the photon mixed with the axion. This energy generates a new effect of the vacuum polarization and it has the state equation of the cosmological constant, $w = -1$. This result holds for any homogeneous and isotropic curved space-time, as well as for diagonal metrics. Numerical estimates of the corrections to the oscillation formulae are presented by considering the intensity of the magnetic field available in the laboratory. Moreover, we estimate the vacuum energy density induced by axion--photon mixing in the Minkowski space-time. A value compatible with that of the energy density of the universe can be obtained for axions with a mass of $(10^{-3}-10^{-2}) eV$ in the presence of the strong magnetic fields that characterize astrophysical objects such as pulsars or neutron stars. In addition, a value of the energy density less than that of the Casimir effect is obtained for magnetic fields used in experiments such as PVLAS. The vacuum polarization induced by this energy could be detected in next experiments and it might provide an indirect proof of the existence of the axion--photon mixing. The quantum field theory effects presented in this work may lead to new methods for studying axion-like particles.


title: Networks with point like nonlinearities

abstract: We study static nonlinear waves in networks described by a nonlinear Schrodinger equation with point-like nonlinearities on metric graphs. Explicit solutions fulfilling vertex boundary conditions are obtained. Spontaneous symmetry breaking caused by bifurcations is found.


title: Efficient n-to-n Collision Detection for Space Debris using 4D AABB Trees (Extended Report)

abstract: Collision detection algorithms are used in aerospace, swarm robotics, automotive, video gaming, dynamics simulation and other domains. As many applications of collision detection run online, timing requirements are imposed on the algorithm runtime: algorithms must, at a minimum, keep up with the passage of time. In practice, this places a limit on the number of objects, n, that can be tracked at the same time. In this paper, we improve the scalability of collision detection, effectively raising the limit n for online object tracking. The key to our approach is the use of a four-dimensional axis-aligned bounding box (AABB) tree, which stores each object's three-dimensional occupancy region in space during a one-dimensional interval of time. This improves efficiency by permitting per-object variable times steps. Further, we describe partitioning strategies that can decompose the 4D AABB tree search into several smaller-dimensional problems that can be solved in parallel. We formalize the collision detection problem and prove our algorithm's correctness. We demonstrate the feasibility of online collision detection for an orbital space debris application, using publicly available data on the full catalog of n=16848 objects provided by www.space-track.org.


title: Physical insights on fractal analysis applied to light curves $\delta$ Scuti stars

abstract: Fractal fingerprints have been found recently in the light curves of several {\delta} Scuti stars observed by CoRoT satellite. This sole fact might pose a problem for the detection of pulsation frequencies using classical pre-whitening techniques, but it is also a potentially rich source for information about physical mechanisms associated to stellar variability. Assuming that a light curve is composed by a superposition of oscillation modes with a fractal background noise, in this work we applied the Coarse Graining Spectral Analysis (CGSA), an FFT based algorithm, which can discriminate in a time series the stochastic fractal power spectra from the harmonic one. We have found that the fractal background component is determining the frequency content extracted using classical pre-whitening techniques in the light curves of {\delta} Scuti stars. This might be crucial to understand the amount of frequencies excited in these kind of pulsating stars. Additionally, CGSA resulted to be relevant in order to extract the oscillation modes, this points to a new criterion to stop the pre-whitening cascade based on the percentage of fractal component in the residuals.


title: Interaction between Molecular Clouds and MeV-TeV Cosmic-ray Protons Escaped from Supernova Remnants

abstract: Recent discovery of the X-ray neutral iron line (Fe I K\alpha at 6.40 keV) around several supernova remnants (SNRs) show that MeV cosmic-ray (CR) protons are distributed around the SNRs and are interacting with neutral gas there. We propose that these MeV CRs are the ones that have been accelerated at the SNRs together with GeV-TeV CRs. In our analytical model, the MeV CRs are still confined in the SNR when the SNR collides with molecular clouds. After the collision, the MeV CRs leak into the clouds and produce the neutral iron line emissions. On the other hand, GeV-TeV CRs had already escaped from the SNRs and emit gamma-rays through interaction with molecular clouds surrounding the SNRs. We apply this model to the SNRs W28 and W44 and show that it can reproduce the observations of the iron line intensities and the gamma-ray spectra. This can be another support of a hadronic scenario for the gamma-ray emissions from these SNRs.


title: SENSEI: Direct-Detection Constraints on Sub-GeV Dark Matter from a Shallow Underground Run Using a Prototype Skipper-CCD

abstract: We present new direct-detection constraints on eV-to-GeV dark matter interacting with electrons using a prototype detector of the Sub-Electron-Noise Skipper-CCD Experimental Instrument. The results are based on data taken in the MINOS cavern at the Fermi National Accelerator Laboratory. We focus on data obtained with two distinct readout strategies. For the first strategy, we read out the Skipper-CCD continuously, accumulating an exposure of 0.177 gram-days. While we observe no events containing three or more electrons, we find a large one- and two-electron background event rate, which we attribute to spurious events induced by the amplifier in the Skipper-CCD readout stage. For the second strategy, we take five sets of data in which we switch off all amplifiers while exposing the Skipper-CCD for 120k seconds, and then read out the data through the best prototype amplifier. We find a one-electron event rate of (3.51 +- 0.10) x 10^(-3) events/pixel/day, which is almost two orders of magnitude lower than the one-electron event rate observed in the continuous-readout data, and a two-electron event rate of (3.18 +0.86 -0.55) x 10^(-5) events/pixel/day. We again observe no events containing three or more electrons, for an exposure of 0.069 gram-days. We use these data to derive world-leading constraints on dark matter-electron scattering for masses between 500 keV to 5 MeV, and on dark-photon dark matter being absorbed by electrons for a range of masses below 12.4 eV.

categories: hep-ex astro-ph.CO hep-ph physics.ins-det


title: A Radio Source Coincident with the Superluminous Supernova PTF10hgi: Evidence for a Central Engine and an Analogue of the Repeating FRB121102?

abstract: We present the detection of an unresolved radio source coincident with the position of the Type I superluminous supernova (SLSN) PTF10hgi ($z=0.098$) about 7.5 years post-explosion, with a flux density of $F_\nu(6\,\,{\rm GHz)}\approx 47.3\ \mu Jy$ and a luminosity of $L_\nu(6\,\,{\rm GHz})\approx 1.1\times 10^{28}$ erg s$^{-1}$ Hz$^{-1}$. This represents the first detection of radio emission coincident with a SLSN on any timescale. We investigate various scenarios for the origin of the radio emission: star formation activity, an active galactic nucleus, and a non-relativistic supernova blastwave. While any of these would be quite novel if confirmed, none appear likely when taken in context of the other properties of the host galaxy, previous radio observations of SLSNe, and the general population of hydrogen-poor SNe. Instead, the radio emission is reminiscent of the quiescent radio source associated with the repeating FRB 121102, which has been argued to be powered by a magnetar born in a SLSN or LGRB explosion several decades ago. We show that the properties of the radio source are consistent with a magnetar wind nebula or an off-axis jet, indicating the presence of a central engine. Our directed search for FRBs from the location of PTF10hgi using 40 min of VLA phased-array data reveals no detections to a limit of $22$ mJy ($10\sigma$; 10 ms duration). We outline several follow-up observations that can conclusively establish the origin of the radio emission.


title: High-pT Signatures in Vector-Leptoquark Models

abstract: We present a detailed analysis of the collider signatures of TeV-scale massive vector bosons motivated by the hints of lepton flavour non-universality observed in $B$-meson decays. We analyse three representations that necessarily appear together in a large class of ultraviolet-complete models: a colour-singlet ($Z'$), a colour-triplet (the $U_1$ leptoquark), and a colour octet ($G'$). Under general assumptions for the interactions of these exotic states with Standard Model fields, including in particular possible right-handed and flavour off-diagonal couplings for the $U_1$, we derive a series of stringent bounds on masses and couplings that constrain a wide range of explicit new-physics models.


title: The Fate of AGB Wind in Massive Galaxies and the ICM

abstract: Asymptotic Giant Branch (AGB) winds from evolved stars not only provide a non-trivial amount of mass and energy return, but also produce dust grains in massive elliptical galaxies. Due to the fast stellar velocity and the high ambient temperature, the wind is thought to form a comet-like tail, similar to Mira in the Local Bubble. Many massive elliptical galaxies and cluster central galaxies host extended dusty cold filaments. The fate of the cold dusty stellar wind and its relation to cold filaments are not well understood. In this work, we carry out both analytical and numerical studies of the interaction between an AGB wind and the surrounding hot gas. We find that the cooling time of the tail is inversely proportional to the ambient pressure. In the absence of cooling, or in low pressure environments (e.g., the outskirts of elliptical galaxies), AGB winds are quickly mixed into the hot gas, and all the AGB winds have similar appearance and head-to-tail ratio. In high pressure environments, such as the Local Bubble and the central regions of massive elliptical galaxies, some of the gas in the mixing layer between the stellar wind and the surrounding hot gas can cool efficiently and cause the tail to become longer. Our simulated tail of Mira itself has similar length and velocity to that observed, and appears similar to the simulated AGB tail in the central regions of massive galaxies. We speculate that instead of thermal instability, the induced condensation at the mixing layer of AGB winds may be the origin of cold filaments in massive galaxies and galaxy clusters. This naturally explains the existence of dust and PAH in the filaments.


title: Source structure and molecular gas properties from high-resolution CO imaging of SPT-selected dusty star-forming galaxies

abstract: We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of high-J CO lines ($J_\mathrm{up}=6$, 7, 8) and associated dust continuum towards five strongly lensed, dusty, star-forming galaxies (DSFGs) at redshift $z = 2.7$-5.7. These galaxies, discovered in the South Pole Telescope survey, are observed at $0.2''$-$0.4''$ resolution with ALMA. Our high-resolution imaging coupled with the lensing magnification provides a measurement of the structure and kinematics of molecular gas in the background galaxies with spatial resolutions down to kiloparsec scales. We derive visibility-based lens models for each galaxy, accurately reproducing observations of four of the galaxies. Of these four targets, three show clear velocity gradients, of which two are likely rotating disks. We find that the reconstructed region of CO emission is less concentrated than the region emitting dust continuum even for the moderate-excitation CO lines, similar to what has been seen in the literature for lower-excitation transitions. We find that the lensing magnification of a given source can vary by 20-50% across the line profile, between the continuum and line, and between different CO transitions. We apply Large Velocity Gradient (LVG) modeling using apparent and intrinsic line ratios between lower-J and high-J CO lines. Ignoring these magnification variations can bias the estimate of physical properties of interstellar medium of the galaxies. The magnitude of the bias varies from galaxy to galaxy and is not necessarily predictable without high resolution observations.


title: Gamma Lines from the Hidden Sector

abstract: We discuss the visibility of gamma lines from dark matter annihilation. We point out a class of theories for dark matter which predict the existence of gamma lines with striking features. In these theories, the final state radiation processes are highly suppressed and one could distinguish easily the gamma lines from the continuum spectrum. We discuss the main experimental bounds and show that one could test the predictions for gamma lines in the near future in the context of simple gauge theories for dark matter.


title: Simplified Models of Flavourful Leptoquarks

abstract: We study the implications of single leptoquark extensions of the Standard Model (SM) under the assumption that their enhanced Yukawa sectors are invariant under global Abelian flavour symmetries already present in SM mass terms. Such symmetries, assumed to be the `residual' subgroups of an ultra-violet flavour theory, have previously been considered in order to predict fermionic mixing angles. Here we focus instead on their effect on the novel flavour structures sourced by the leptoquark representations that address the present ${R}_{K^{(\star)}}$ anomalies in semileptonic rare $B$-decays. Combined with existing flavour data, the residual symmetries prove to be extremely constraining; we find that the (quark-lepton) leptoquark Yukawa couplings fall within $\mathcal{O}(10)$ highly predictive patterns, each with only a single free parameter, when `normal' (SM-like) hierarchies are assumed. In addition, proton decay for the scalar SU(2) triplet representation is naturally avoided in the residual symmetry approach without relying on further model building. Our results indicate that a simultaneous explanation for the ${R}_{K^{(\star)}}$ anomalies and the flavour puzzle may be achieved in a simplified, model-independent formalism.


title: Magic Angle Hierarchy in Twisted Graphene Multilayers

abstract: When two monolayers of graphene are stacked with a small relative twist angle, the resulting band structure exhibits a remarkably flat pair of bands at a sequence of 'magic angles' where correlation effects can induce a host of exotic phases. Here, we study a class of related models of $n$-layered graphene with alternating relative twist angle $\pm \theta$ which exhibit magic angle flat bands coexisting with several Dirac dispersing bands at the Moir\'e K point. Remarkably, we find that the Hamiltonian for the multilayer system can be mapped exactly to a set of decoupled bilayers at {\it different} angles, revealing a remarkable hierarchy mathematically relating all these magic angles to the TBG case. For the trilayer case ($n = 3$), we show that the sequence of magic angle is obtained by multiplying the bilayer magic angles by $\sqrt{2}$, whereas the quadrilayer case ($n = 4$) has two sequences of magic angles obtained by multiplying the bilayer magic angles by the golden ratio $\varphi = (\sqrt{5} + 1)/2 \approx 1.62$ and its inverse. We also show that for larger $n$, we can tune the angle to obtain several narrow (almost flat) bands simultaneously and that for $n \rightarrow \infty$, there is a continuum of magic angles for $\theta \lesssim 2^o$. Furthermore, we show that tuning several perfectly flat bands for a small number of layers is possible if the coupling between different layers is different. The setup proposed here can be readily achieved by repeatedly applying the "tear and stack" method without the need of any extra tuning of the twist angle and has the advantage that the first magic angle is always larger than the bilayer case.


title: Entanglement Spectra of Stabilizer Codes: A Window into Gapped Quantum Phases of Matter

abstract: The entanglement spectrum (ES) provides a barometer of quantum entanglement and encodes physical information beyond that contained in the entanglement entropy. In this paper, we explore the ES of stabilizer codes, which furnish exactly solvable models for a plethora of gapped quantum phases of matter. Studying the ES for stabilizer Hamiltonians in the presence of arbitrary weak local perturbations thus allows us to develop a general framework within which the entanglement features of gapped topological phases can be computed and contrasted. In particular, we study models harboring fracton order, both type-I and type-II, and compare the resulting ES with that of both conventional topological order and of (strong) subsystem symmetry protected topological (SSPT) states. We find that non-local surface stabilizers (NLSS), a set of symmetries of the Hamiltonian which form on the boundary of the entanglement cut, act as purveyors of universal non-local features appearing in the entanglement spectrum. While in conventional topological orders and fracton orders, the NLSS retain a form of topological invariance with respect to the entanglement cut, subsystem symmetric systems---fracton and SSPT phases---additionally show a non-trivial geometric dependence on the entanglement cut, corresponding to the subsystem symmetry. This sheds further light on the interplay between geometric and topological effects in fracton phases of matter and demonstrates that strong SSPT phases harbour a measure of quasi-local entanglement beyond that encountered in conventional SPT phases. We further show that a version of the edge-entanglement correspondence, established earlier for gapped two-dimensional topological phases, also holds for gapped three-dimensional fracton models.


title: Discovery of the First Low-Luminosity Quasar at z > 7

abstract: We report the discovery of a quasar at z = 7.07, which was selected from the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. This quasar, HSC J124353.93+010038.5, has an order of magnitude lower luminosity than do the other known quasars at z > 7. The rest-frame ultraviolet absolute magnitude is M1450 = -24.13 +/- 0.08 mag and the bolometric luminosity is Lbol = (1.4 +/- 0.1) x 10^{46} erg/s. Its spectrum in the optical to near-infrared shows strong emission lines, and shows evidence for a fast gas outflow, as the C IV line is blueshifted and there is indication of broad absorption lines. The Mg II-based black hole mass is Mbh = (3.3 +/- 2.0) x 10^8 Msun, thus indicating a moderate mass accretion rate with an Eddington ratio 0.34 +/- 0.20. It is the first z > 7 quasar with sub-Eddington accretion, besides being the third most distant quasar, known to date. The luminosity and black hole mass are comparable to, or even lower than, those measured for the majority of low-z quasars discovered by the Sloan Digital Sky Survey, and thus this quasar likely represents a z > 7 counterpart to quasars commonly observed in the low-z universe.


title: Initialisation of single spin dressed states using shortcuts to adiabaticity

abstract: We demonstrate the use of shortcuts to adiabaticity protocols for initialisation, readout, and coherent control of dressed states generated by closed-contour, coherent driving of a single spin. Such dressed states have recently been shown to exhibit efficient coherence protection, beyond what their two-level counterparts can offer. Our state transfer protocols yield a transfer fidelity of ~ 99.4(2) % while accelerating the transfer speed by a factor of 2.6 compared to the adiabatic approach. We show bi-directionality of the accelerated state transfer, which we employ for direct dressed state population readout after coherent manipulation in the dressed state manifold. Our results enable direct and efficient access to coherence-protected dressed states of individual spins and thereby offer attractive avenues for applications in quantum information processing or quantum sensing.


title: Pixelated Dark Energy

abstract: We study the phenomenology of a recent string construction with a quantum mechanically stable dark energy. A mild supersymmetry protects the vacuum energy but also allows $O(10 - 100)$ TeV scale superpartner masses. The construction is holographic in the sense that the 4D spacetime is generated from "pixels" originating from five-branes wrapped over metastable five-cycles of the compactification. The cosmological constant scales as $\Lambda \sim 1/N$ in the pixel number. An instability in the construction leads to cosmic expansion. This also causes more five-branes to wind up in the geometry, leading to a slowly decreasing cosmological constant which we interpret as an epoch of inflation followed by (pre-)heating when a rare event occurs in which the number of pixels increases by an order one fraction. The sudden appearance of radiation triggers an exponential increase in the number of pixels. Dark energy has a time varying equation of state with $w_a=-3\Omega_{m,0}(1+w_0)/2$, which is compatible with current bounds, and could be constrained further by future data releases. The pixelated nature of the Universe also implies a large-$l$ cutoff on the angular power spectrum of cosmological observables with $l_{\rm max} \sim O(N)$. We also use this pixel description to study the thermodynamics of de Sitter space, finding rough agreement with effective field theory considerations.


title: The build-up of pseudobulges in a hierarchical universe

abstract: We study the cosmological build-up of pseudobulges using the LGalaxies semi-analytical model for galaxy formation with a new approach for following separately the assembly of classical bulges and pseudobulges. Classical bulges are assumed to be the result of violent processes (i.e., mergers and starbursts), while the formation of pseudobulges is connected to the secular growth of disks. We apply the model to both the Millennium and the Millennium II simulations, in order to study our results across a wide range of stellar masses ($10^{7} - 10^{11.5} M_{\odot}$). We find that $z=0$ pseudobulges mainly reside in galaxies of $ \rm M_{stellar} \sim 10^{10} - 10^{10.5} M_{\odot}$ ($\rm M_{halo} \sim 10^{11.5}-10^{12} M_{\odot}$) and we recover structural properties of these objects (e.g., sizes and bulge-to-total ratios) that are in good agreement with observational results. Tracing their formation history, we find that pseudobulges assembled in galaxies with a very quiet merger history, as opposed to the host galaxies of classical bulges. Regarding the bulge structure, we find that $\sim$ 30\% of the galaxies with a predominant pseudobulge feature a composite structure, hosting both a pseudo and a classical bulge component. The classical component typically constitutes $\sim$10\% of the total bulge galaxy mass. When looking at the properties of the host galaxies, we find that $z = 0$ pseudobulges are hosted by main sequence galaxies, characterized by a stellar population which is generally younger compared to the one of the hosts of classical bulges.


title: A ~60-day super-orbital period originating from the ultraluminous X-ray pulsar in M82

abstract: Ultraluminous X-ray (ULX) pulsars are a new class of object powered by apparent super-critical accretion onto magnetized neutron stars. Three sources in this class identified so far; M82 X-2, NGC 5907 ULX-1 and NGC 7793 P13, have been found to have two properties in common; $\sim1$-s spin periods, and for NGC 5907 ULX-1 and NGC 7793 P13 periodic X-ray flux modulations on timescales of $\sim60-80$ days. M82 X-2 resides in a crowded field that includes the ULX M82 X-1 separated from X-2 by 5", as well as other bright point sources. A 60-day modulation has been observed from the region but the origin has been difficult to identify; both M82 X-1 and X-2 have been suggested as the source. In this paper we present the analysis of a systematic monitoring campaign by Chandra, the only X-ray telescope capable of resolving the crowded field. From a simple Lomb-Scargle periodogram analysis and a more sophisticated Gaussian Process analysis we find that only X-2 exhibits a periodic signal around 60 days supporting previous claims that it is the origin. We also construct a phase-averaged flux profile of the modulations from higher cadence Swift/XRT data and find that the flux variations in the Chandra data are fully consistent with the flux profile. Since the orbit of the neutron star and its companion is known to be 2.5 days, the $\sim60$-day period must be super-orbital in origin. The flux of the modulations varies by a factor of $\sim$100 from minimum to maximum, with no evidence for spectral variations, making the origin difficult to explain.


title: Constructing phase space distributions with internal symmetries

abstract: We discuss an ab initio world-line approach to constructing phase space distributions in systems with internal symmetries. Starting from the Schwinger-Keldysh real time path integral in quantum field theory, we derive the most general extension of the Wigner phase space distribution to include color and spin degrees of freedom in terms of dynamical Grassmann variables. The corresponding Liouville distribution for colored particles, which obey Wong's equation, has only singlet and octet components, while higher moments are fully constrained by the Grassmann algebra. The extension of phase space dynamics to spin is represented by a generalization of the Pauli-Lubanski vector; its time evolution via the Bargmann-Michel-Telegdi equation also follows from the phase space trajectories of the underlying Grassmann coordinates. Our results for the Liouville phase space distribution in systems with both spin and color are of interest in fields as diverse as chiral fluids, finite temperature field theory and polarized parton distribution functions. We also comment on the role of the chiral anomaly in the phase space dynamics of spinning particles.


title: Variations on S-fold CFTs

abstract: A local SL(2,Z) transformation on the Type IIB brane configuration gives rise to an interesting class of superconformal field theories, known as the S-fold CFTs. Previously it has been proposed that the corresponding quiver theory has a link involving the T(U(N)) theory. In this paper, we generalise the preceding result by studying quivers that contain a T(G) link, where G is self-dual under S-duality. In particular, the cases of G = SO(2N), USp'(2N) and G_2 are examined in detail. We propose the theories that arise from an appropriate insertion of an S-fold into a brane system, in the presence of an orientifold threeplane or an orientifold fiveplane. By analysing the moduli spaces, we test such a proposal against its S-dual configuration using mirror symmetry. The case of G_2 corresponds to a novel class of quivers, whose brane construction is not available. We present several mirror pairs, containing G_2 gauge groups, that have not been discussed before in the literature.


title: Search for Semi-Annihilating Dark Matter with Fermi-LAT, H.E.S.S., Planck, and the Cherenkov Telescope Array

abstract: Dark matter annihilations have been fiercely restricted by the non-observation of excess events in indirect detection probes. Typically the interactions that dictate annihilation are also present in the dark matter-nucleon scattering cross section, also severely constrained by direct detection experiments. Semi-annihilations arise as a possible way to ameliorate the tension and even change the standard dark matter relic density calculation. In this work, we derive indirect detection bounds for several semi-annihilation channels including gauge bosons, the Higgs, leptophilic and leptophobic scalars. Our analysis is based on the gamma-ray observations in the direction of Dwarf Spheroidal Galaxies (Fermi-LAT) and the Galactic Center (H.E.S.S.), and Planck measurements of Cosmic Background Radiation. In addition, we derive the prospects for the Cherenkov Telescope Array (CTA) sensitivity to all these semi-annihilation modes.


title: Generating a second-order topological insulator with multiple corner states by periodic driving

abstract: We study the effects of periodic driving on a variant of the Bernevig-Hughes-Zhang (BHZ) model defined on a square lattice. In the absence of driving, the model has both topological and nontopological phases depending on the different parameter values. We also study the anisotropic BHZ model and show that, unlike the isotropic model, it has a nontopological phase which has states localized on only two of the four edges of a finite-sized square. When an appropriate term is added, the edge states get gapped and gapless states appear at the four corners of a square; we have shown that these corner states can be labeled by the eigenvalues of a certain operator. When the system is driven periodically by a sequence of two pulses, multiple corner states may appear depending on the driving frequency and other parameters. We discuss to what extent the system can be characterized by topological invariants such as the Chern number and a diagonal winding number. We have shown that the locations of the jumps in these invariants can be understood in terms of the Floquet operator at both the time-reversal invariant momenta and other momenta which have no special symmetries.


title: Impact of Training Dataset Size on Neural Answer Selection Models

abstract: It is held as a truism that deep neural networks require large datasets to train effective models. However, large datasets, especially with high-quality labels, can be expensive to obtain. This study sets out to investigate (i) how large a dataset must be to train well-performing models, and (ii) what impact can be shown from fractional changes to the dataset size. A practical method to investigate these questions is to train a collection of deep neural answer selection models using fractional subsets of varying sizes of an initial dataset. We observe that dataset size has a conspicuous lack of effect on the training of some of these models, bringing the underlying algorithms into question.


title: $\mathcal{N}=2^*$ gauge theory, free fermions on the torus and Painlev\'e VI

abstract: In this paper we study the extension of Painlev\'e/gauge theory correspondence to circular quivers by focusing on the special case of $SU(2)$ $\mathcal{N}=2^*$ theory. We show that the Nekrasov-Okounkov partition function of this gauge theory provides an explicit combinatorial expression and a Fredholm determinant formula for the tau-function describing isomonodromic deformations of $SL_2$ flat connections on the one-punctured torus. This is achieved by reformulating the Riemann-Hilbert problem associated to the latter in terms of chiral conformal blocks of a free-fermionic algebra. This viewpoint provides the exact solution of the renormalization group flow of the $SU(2)$ $\mathcal{N}=2^*$ theory on self-dual $\Omega$-background and, in the Seiberg-Witten limit, an elegant relation between the IR and UV gauge couplings.


title: Flavor Changing Heavy Higgs Interactions with Leptons at Hadron Colliders

abstract: In a general two Higgs doublet model, we study flavor changing neutral Higgs (FCNH) decays into leptons at hadron colliders, $pp \to \phi^0 \to \tau^\mp\mu^\pm +X$, where $\phi^0$ could be a CP-even scalar ($h^0$, $H^0$) or a CP-odd pseudoscalar ($A^0$). The light Higgs boson $h^0$ is found to resemble closely the Standard Model Higgs boson at the Large Hadron Collider. In the alignment limit of $\cos(\beta-\alpha) \cong 0$ for $h^0$--$H^0$ mixing, FCNH couplings of $h^0$ are naturally suppressed, but such couplings of the heavier $H^0, A^0$ are sustained by $\sin(\beta-\alpha) \simeq 1$. We evaluate physics backgrounds from dominant processes with realistic acceptance cuts and tagging efficiencies. We find promising results for $\sqrt{s} = 14$ TeV, which we extend further to $\sqrt{s} = 27$ TeV and 100 TeV future pp colliders.


title: Near-Unitary Spin Squeezing in $^{171}$Yb

abstract: Spin squeezing can improve atomic precision measurements beyond the standard quantum limit (SQL), and unitary spin squeezing is essential for improving atomic clocks. We report substantial and nearly unitary spin squeezing in $^{171}$Yb, an optical lattice clock atom. The collective nuclear spin of $\sim 10^3$ atoms is squeezed by cavity feedback, using light detuned from the system's resonances to attain unitarity. The observed precision gain over the SQL is limited by state readout to 6.5(4) dB, while the generated states offer a gain of 12.9(6) dB, limited by the curvature of the Bloch sphere. Using a squeezed state within 30% of unitarity, we demonstrate an interferometer that improves the averaging time over the SQL by a factor of 3.7(2). In the future, the squeezing can be simply transferred onto the optical clock transition of $^{171}$Yb.


title: Discretizing Continuous Action Space for On-Policy Optimization

abstract: In this work, we show that discretizing action space for continuous control is a simple yet powerful technique for on-policy optimization. The explosion in the number of discrete actions can be efficiently addressed by a policy with factorized distribution across action dimensions. We show that the discrete policy achieves significant performance gains with state-of-the-art on-policy optimization algorithms (PPO, TRPO, ACKTR) especially on high-dimensional tasks with complex dynamics. Additionally, we show that an ordinal parameterization of the discrete distribution can introduce the inductive bias that encodes the natural ordering between discrete actions. This ordinal architecture further significantly improves the performance of PPO/TRPO.


title: Repairing without Retraining: Avoiding Disparate Impact with Counterfactual Distributions

abstract: When the performance of a machine learning model varies over groups defined by sensitive attributes (e.g., gender or ethnicity), the performance disparity can be expressed in terms of the probability distributions of the input and output variables over each group. In this paper, we exploit this fact to reduce the disparate impact of a fixed classification model over a population of interest. Given a black-box classifier, we aim to eliminate the performance gap by perturbing the distribution of input variables for the disadvantaged group. We refer to the perturbed distribution as a counterfactual distribution, and characterize its properties for common fairness criteria. We introduce a descent algorithm to learn a counterfactual distribution from data. We then discuss how the estimated distribution can be used to build a data preprocessor that can reduce disparate impact without training a new model. We validate our approach through experiments on real-world datasets, showing that it can repair different forms of disparity without a significant drop in accuracy.


title: Sub-ballistic growth of R\'enyi entropies due to diffusion

abstract: We investigate the dynamics of quantum entanglement after a global quench and uncover a qualitative difference between the behavior of the von Neumann entropy and higher R\'enyi entropies. We argue that the latter generically grow \emph{sub-ballistically}, as $\propto\sqrt{t}$, in systems with diffusive transport. We provide strong evidence for this in both a U$(1)$ symmetric random circuit model and in a paradigmatic non-integrable spin chain, where energy is the sole conserved quantity. We interpret our results as a consequence of local quantum fluctuations in conserved densities, whose behavior is controlled by diffusion, and use the random circuit model to derive an effective description. We also discuss the late-time behavior of the second R\'enyi entropy and show that it exhibits hydrodynamic tails with \emph{three distinct power laws} occurring for different classes of initial states.


title: Time-Space tradeoff in deep learning models for crop classification on satellite multi-spectral image time series

abstract: In this article, we investigate several structured deep learning models for crop type classification on multi-spectral time series. In particular, our aim is to assess the respective importance of spatial and temporal structures in such data. With this objective, we consider several designs of convolutional, recurrent, and hybrid neural networks, and assess their performance on a large dataset of freely available Sentinel-2 imagery. We find that the best-performing approaches are hybrid configurations for which most of the parameters (up to 90%) are allocated to modeling the temporal structure of the data. Our results thus constitute a set of guidelines for the design of bespoke deep learning models for crop type classification.


title: A Froggatt-Nielsen flavor model for neutrino physics

abstract: Superheavy neutrinos can, via the seesaw model, provide a mechanism for lepton number violation. If they are combined with flavor violation as characterized by the Froggatt-Nielsen mechanism, then the phenomenology for the neutrinos in oscillation experiments, neutrinoless double beta decay, and other experiments can be described by a relatively few number of parameters. We describe the low-energy neutrino mass matrix and show that the results are consistent with currently available data.


title: A/B Testing in Dense Large-Scale Networks: Design and Inference

abstract: Design of experiments and estimation of treatment effects in large-scale networks, in the presence of strong interference, is a challenging and important problem. Most existing methods' performance deteriorates as the density of the network increases. In this paper, we present a novel strategy for accurately estimating the causal effects of a class of treatments in a dense large-scale network. First, we design an approximate randomized controlled experiment by solving an optimization problem to allocate treatments in the presence of competition among neighboring nodes. Then we apply an importance sampling adjustment to correct for any leftover bias (from the approximation) in estimating average treatment effects. We provide theoretical guarantees, verify robustness in a simulation study, and validate the scalability and usefulness of our procedure in a real-world experiment on a large social network.


title: Initial state description of azimuthally collimated long range correlations in ultrarelativistic light-heavy ion collisions

abstract: It was argued in arXiv:1805.09342 and arXiv:1807.00825 that the systematics of the azimuthal anisotropy coefficients $v_{2,3}$ measured in ultrarelativistic light-heavy ion collisions at RHIC and the LHC can be described in an initial state dilute-dense Color Glass Condensate (CGC) framework. We elaborate here on the discussion in these papers and provide further novel results that strengthen their conclusions. The underlying mathematical framework and numerical techniques employed are very similar to those in the CGC based IP-Glasma model used previously as initial conditions for heavy-ion collisions. The uncertainties in theory/data comparisons for small systems are discussed, with unknowns that are specific to the model distinguished from those that are generic to all models. We present analytical arguments that demonstrate that quantum interference effects such as Bose enhancement and Hanbury-Brown-Twiss correlations of gluons, as well as coherent multiple scattering of gluons in the projectile off color domains in the target, are enhanced in rare events. The quantum origins of the large anisotropies in small systems are corroborated by numerical results for deuteron-gold collisions that show that large anisotropies in rare configurations can occur when the nucleons in the projectile overlap significantly. This is at variance with the classical intuition of hydrodynamical models. We also comment on the consequences of ignoring the many-body color charge correlations of gluons in models that only consider geometrical fluctuations in the energy density.


title: Limiting Probability Measures

abstract: The coordinates along any fixed direction(s), of points on the sphere $S^{n-1}(\sqrt{n})$, roughly follow a standard Gaussian distribution as $n$ approaches infinity. We revisit this classical result from a nonstandard analysis perspective, providing a new proof by working with hyperfinite dimensional spheres. We also set up a nonstandard theory for the asymptotic behavior of integrals over varying domains in general. We obtain a new proof of the Riemann--Lebesgue lemma as a by-product of this theory. We finally show that for any function $f \co \mathbb{R}^k \to \mathbb{R}$ with finite Gaussian moment of an order larger than one, its expectation is given by a Loeb integral integral over a hyperfinite dimensional sphere. Some useful inequalities between high-dimensional spherical means of $f$ and its Gaussian mean are obtained in order to complete the above proof. A review of the requisite nonstandard analysis is provided.


title: Near-field Ultra-wideband mmWave Channel Characterization Using Successive Cancellation Beamspace UCA Algorithm

abstract: Of the wide palette of 5G features, ultra-wide bandwidth and large-scale antenna configuration are regarded as the essential enabling technology components at millimeter wave (mmWave) communication. Accurate knowledge of delay and angle information of multipath components is essential for many applications in mmWave systems. There is a strong need for a low computation-cost channel estimation algorithm for such systems, where typically adopted far-field and narrowband assumptions might be violated. In this work, a generic yet novel beamspace uniform circular array (UCA) beamforming algorithm with successive cancellation scheme is proposed to jointly detect the impinging angle and delay of the multipath components. The proposed algorithm is computationally cheap and it works for ultra-wideband (UWB) systems in the near-field conditions. Both numerical simulations and experimental validation results are provided to demonstrate the effectiveness and robustness of the proposed algorithm, compared to the state-of-art works.


title: Quantum state engineering by non-deterministic noiseless linear amplification

abstract: We address quantum state engineering of single- and two-mode states by means of non-deterministic noiseless linear amplifiers (NLAs) acting on Gaussian states. In particular, we show that NLAs provide an effective scheme to generate highly non-Gaussian and non-classical states. Additionally, we show that the amplification of a two-mode squeezed vacuum state (twin-beam) may highly increase entanglement.


title: Quantified Uncertainty in Thermodynamic Modeling for Materials Design

abstract: Phase fractions, compositions and energies of the stable phases as a function of macroscopic composition, temperature, and pressure (X-T-P) are the principle correlations needed for the design of new materials and improvement of existing materials. They are the outcomes of thermodynamic modeling based on the CALculation of PHAse Diagrams (CALPHAD) approach. The accuracy of CALPHAD predictions vary widely in X-T-P space due to experimental error, model inadequacy and unequal data coverage. In response, researchers have developed frameworks to quantify the uncertainty of thermodynamic property model parameters and propagate it to phase diagram predictions. In previous studies, uncertainty was represented as intervals on phase boundaries (with respect to composition) or invariant reactions (with respect to temperature) and was unable to represent the uncertainty in eutectoid reactions or in the stability of phase regions. In this work, we propose a suite of tools that leverages samples from the multivariate model parameter distribution to represent uncertainty in forms that surpass previous limitations and are well suited to materials design. These representations include the distribution of phase diagrams and their features, as well as the dependence of phase stability and the distributions of phase fraction, composition activity and Gibbs energy on X-T-P location - irrespective of the total number of components. Most critically, the new methodology allows the material designer to interrogate a certain composition and temperature domain and get in return the probability of different phases to be stable, which can positively impact materials design.


title: Eta-quotients of Prime or Semiprime Level and Elliptic Curves

abstract: From the Modularity Theorem proven by Wiles, Taylor, et al, we know that all elliptic curves are modular. It has been shown by Martin and Ono exactly which are represented by eta-quotients, and some examples of elliptic curves represented by modular forms that are linear combinations of eta-quotients have been given by Pathakjee, RosnBrick, and Yoong. In this paper, we first show that eta-quotients which are modular for any congruence subgroup of level $N$ coprime to $6$ can be viewed as modular for $\Gamma_0(N)$. We then categorize when even weight eta-quotients can exist in $M_k (\Gamma_1(p))$ and $M_k (\Gamma_1(pq))$, for distinct primes $p,q$. We conclude by providing some new examples of elliptic curves whose corresponding modular forms can be written as a linear combination of eta-quotients, and describe an algorithmic method for finding additional examples.


title: Tuning friction at material-nanoparticle-liquid interfaces with an external electric field

abstract: The use of electrophoretic forces to tune friction at material-nanoparticle-liquid interfaces with static or low frequency (0.6-50 mHz) electric fields is reported for the first time. External electric fields were employed to reposition negatively charged TiO2 or positively charged Al2O3 nanoparticles suspended in water in directions perpendicular to a planar platinum surface of a quartz crystal microbalance, which was then used to monitor frictional shear forces at the interface. Active electro-tunable control of friction has been demonstrated for both TiO2 and Al2O3 suspensions. For TiO2 suspensions, significant drops in frictional shear forces, not observed for Al2O3, were likely attributed to the presence of molecularly thin interstitial water layers remaining in regions between the TiO2 particles and the substrate. Timescales associated with motion of nanoparticles in directions perpendicular to the surface were also investigated by varying the frequency of the external electric field, and were determined to be similar to those of glass-like or polymeric materials. Overall, the studies reveal that nanoparticles actively driven by electric fields can act as "cantilever-free" atomic force probes capable of "tapping mode" exploration of interfacial properties and nanoscale interactions in geometries inaccessible to optical and micromechanical probes.


title: Adversarial Examples Are a Natural Consequence of Test Error in Noise

abstract: Over the last few years, the phenomenon of adversarial examples --- maliciously constructed inputs that fool trained machine learning models --- has captured the attention of the research community, especially when the adversary is restricted to small modifications of a correctly handled input. Less surprisingly, image classifiers also lack human-level performance on randomly corrupted images, such as images with additive Gaussian noise. In this paper we provide both empirical and theoretical evidence that these are two manifestations of the same underlying phenomenon, establishing close connections between the adversarial robustness and corruption robustness research programs. This suggests that improving adversarial robustness should go hand in hand with improving performance in the presence of more general and realistic image corruptions. Based on our results we recommend that future adversarial defenses consider evaluating the robustness of their methods to distributional shift with benchmarks such as Imagenet-C.


title: Hyperspherical Prototype Networks

abstract: This paper introduces hyperspherical prototype networks, which unify classification and regression with prototypes on hyperspherical output spaces. For classification, a common approach is to define prototypes as the mean output vector over training examples per class. Here, we propose to use hyperspheres as output spaces, with class prototypes defined a priori with large margin separation. We position prototypes through data-independent optimization, with an extension to incorporate priors from class semantics. By doing so, we do not require any prototype updating, we can handle any training size, and the output dimensionality is no longer constrained to the number of classes. Furthermore, we generalize to regression, by optimizing outputs as an interpolation between two prototypes on the hypersphere. Since both tasks are now defined by the same loss function, they can be jointly trained for multi-task problems. Experimentally, we show the benefit of hyperspherical prototype networks for classification, regression, and their combination over other prototype methods, softmax cross-entropy, and mean squared error approaches.


title: The structure of deterministic mass and surface fractals: theory and methods of analyzing small-angle scattering data

abstract: Small-angle scattering (SAS) of X-rays, neutrons or light from ensembles of randomly oriented and placed deterministic fractal structures are studied theoretically. In the standard analysis, a very few parameters can be determined from SAS data: the fractal dimension, and the lower and upper limits of the fractal range. The self-similarity of deterministic structures allows one to obtain additional characteristics of their spatial structures. The paper considers models which can describe accurately SAS from such structures. The developed models of deterministic fractals offer many advantages in describing fractal systems, including the possibility to extract additional structural information, an analytic description of SAS intensity, and effective computational algorithms. Generalized Cantor fractals and few of its variants are used as basic examples to illustrate the above concepts and to model physical samples with mass, surface, and multi-fractal structures. The differences between the deterministic and random fractal structures in analyzing SAS data are emphasized. Several limitations are identified in order to motivate future investigations of deterministic fractal structures.


title: The FFBS Estimation of High Dimensional Panel Data Factor Stochastic Volatility Models

abstract: In this paper, We propose a new style panel data factor stochastic volatility model with observable factors and unobservable factors based on the multivariate stochastic volatility model, which is mainly composed of three parts, such as the mean equation, volatility equation and factor volatility evolution. The stochastic volatility equation is a 1-step forward prediction process with high dimensional parameters to be estimated. Using the Markov Chain Monte Carlo Simulation (MCMC) method, the Forward Filtering Backward Sampling (FFBS) algorithm of the stochastic volatility equation is mainly used to estimate the new model by Kalman Filter Recursive Algorithm (KFRA). The results of numeric simulation and latent factor estimation show that the algorithm possesses robustness and consistency for parameter estimation. This paper makes a comparative analysis of the observable and unobservable factors of internet finance and traditional financial listed companies in the Chinese stock market using the new model and its estimation method. The results show that the influence of observable factors is similar to the two types of listed companies, but the influence of unobservable factors is obviously different.


title: Reparameterizable Subset Sampling via Continuous Relaxations

abstract: Many machine learning tasks require sampling a subset of items from a collection based on a parameterized distribution. The Gumbel-softmax trick can be used to sample a single item, and allows for low-variance reparameterized gradients with respect to the parameters of the underlying distribution. However, stochastic optimization involving subset sampling is typically not reparameterizable. To overcome this limitation, we define a continuous relaxation of subset sampling that provides reparameterization gradients by generalizing the Gumbel-max trick. We use this approach to sample subsets of features in an instance-wise feature selection task for model interpretability, subsets of neighbors to implement a deep stochastic k-nearest neighbors model, and sub-sequences of neighbors to implement parametric t-SNE by directly comparing the identities of local neighbors. We improve performance in all these tasks by incorporating subset sampling in end-to-end training.


title: Spectral caustic in two-color high harmonic generation: role of Coulomb effects

abstract: The high-order harmonic spectrum generated in two-color intense laser field under certain conditions has a pronounced maximum caused by the so called spectral caustic. Using a numerical solution of the 3D time-dependent Schr\"odinger equation we study the width of the maximum and the degree of the generation enhancement due to the caustic as a function of the fundamental intensity, wavelength, and atomic ionization potential. It is shown that the degree of this enhancement can be well quantitatively characterized by a single parameter: the ratio of the radius of the free electronic oscillation in the laser field to the atomic size. This behavior can be explained by the Coulomb attraction of the photoelectron to the parent ion, whereas the effect of this attraction on the width of the maximum is negligible. Choosing the field parameters providing a pronounced enhancement, we calculate the high harmonic macroscopic response taking into account the transient phase-matching. We show that the caustic feature can be used to provide a quasi-monochromatic XUV source with an almost linear frequency vs. time dependence. Such source can be used, in particular, for the time-resolved single-shot XUV imaging.


title: On modifications of the exponential integral with the Mittag-Leffler function

abstract: In this paper we survey the properties of the Schelkunoff modification of the Exponential integral and we generalize it with the Mittag-Leffler function. So doing we get a new special function (as far as we know) that may be relevant in linear viscoelasticity because of its complete monotonicity properties in the time domain. We also consider the generalized sine and cosine integral functions


title: A Discrete Approximation to Gibbs Free Energy of Chemical Reactions is Needed for Accurately Calculating Entropy Production in Mesoscopic Simulations

abstract: In modeling the interior of cells by simulating a reaction-diffusion master equation over a grid of compartments, one employs the assumption that the copy numbers of various chemical species are small, discrete quantities. We show that in this case, textbook expressions for the change in Gibbs free energy accompanying a chemical reaction or diffusion between adjacent compartments become inaccurate. We derive exact expressions for these free energy changes under the assumption of discrete copy numbers and illustrate how these expressions reduce to the textbook expressions under a series of successive approximations leveraging the relative sizes of the stoichiometric coefficients and the copy numbers of the solutes and solvent. Numerical results are presented to corroborate the claim that if the copy numbers are treated as discrete quantities, then only these more exact expressions lead to correct equilibrium behavior. The newly derived expressions are critical for correctly tracking dissipation and entropy production in mesoscopic simulations based on the reaction-diffusion master equation formalism.


title: Spectral Multi-scale Community Detection in Temporal Networks with an Application

abstract: The analysis of temporal networks has a wide area of applications in a world of technological advances. An important aspect of temporal network analysis is the discovery of community structures. Real data networks are often very large and the communities are observed to have a hierarchical structure referred to as multi-scale communities. Changes in the community structure over time might take place either at one scale or across all scales of the community structure. The multilayer formulation of the modularity maximization (MM) method introduced captures the changing multi-scale community structure of temporal networks. This method introduces a coupling between communities in neighboring time layers by allowing inter-layer connections, while different values of the resolution parameter enable the detection of multi-scale communities. However, the range of this parameter's values must be manually selected. When dealing with real life data, communities at one or more scales can go undiscovered if appropriate parameter ranges are not selected. A novel Temporal Multi-scale Community Detection (TMSCD) method overcomes the obstacles mentioned above. This is achieved by using the spectral properties of the temporal network represented as a multilayer network. In this framework we select automatically the range of relevant scales within which multi-scale community partitions are sought.


title: Coping with Selection Effects: A Primer on Regression with Truncated Data

abstract: The finite sensitivity of instruments or detection methods means that data sets in many areas of astronomy, for example cosmological or exoplanet surveys, are necessarily systematically incomplete. Such data sets, where the population being investigated is of unknown size and only partially represented in the data, are called "truncated" in the statistical literature. Truncation can be accounted for through a relatively straightforward modification to the model being fitted in many circumstances, provided that the model can be extended to describe the population of undetected sources. Here I examine the problem of regression using truncated data in general terms, and use a simple example to show the impact of selecting a subset of potential data on the dependent variable, on the independent variable, and on a second dependent variable that is correlated with the variable of interest. Special circumstances in which selection effects are ignorable are noted. I also comment on computational strategies for performing regression with truncated data, as an extension of methods that have become popular for the non-truncated case, and provide some general recommendations.


title: NADA-FLD: A General Relativistic, Multi-dimensional Neutrino-hydrodynamics Code Employing Flux-limited Diffusion

abstract: We present the new code NADA-FLD to solve multi-dimensional neutrino-hydrodynamics in full general relativity (GR) in spherical polar coordinates. The energy-dependent neutrino transport assumes the flux-limited diffusion (FLD) approximation and evolves the neutrino energy densities measured in the frame comoving with the fluid. Operator splitting is used to avoid multi-dimensional coupling of grid cells in implicit integration steps involving matrix inversions. Terms describing lateral diffusion and advection are integrated explicitly using the Allen-Cheng or the Runge-Kutta-Legendre method, which remain stable even in the optically thin regime. We discuss several toy-model problems in one and two dimensions to test the basic functionality and individual components of the transport scheme. We also perform fully dynamic core-collapse supernova (CCSN) simulations in spherical symmetry. For a Newtonian model we find good agreement with the M1 code ALCAR, and for a GR model we reproduce the main effects of GR in CCSNe already found by previous works.


title: On the Transferability of Spectral Graph Filters

abstract: This paper focuses on spectral filters on graphs, namely filters defined as elementwise multiplication in the frequency domain of a graph. In many graph signal processing settings, it is important to transfer a filter from one graph to another. One example is in graph convolutional neural networks (ConvNets), where the dataset consists of signals defined on many different graphs, and the learned filters should generalize to signals on new graphs, not present in the training set. A necessary condition for transferability (the ability to transfer filters) is stability. Namely, given a graph filter, if we add a small perturbation to the graph, then the filter on the perturbed graph is a small perturbation of the original filter. It is a common misconception that spectral filters are not stable, and this paper aims at debunking this mistake. We introduce a space of filters, called the Cayley smoothness space, that contains the filters of state-of-the-art spectral filtering methods, and whose filters can approximate any generic spectral filter. For filters in this space, the perturbation in the filter is bounded by a constant times the perturbation in the graph, and filters in the Cayley smoothness space are thus termed linearly stable. By combining stability with the known property of equivariance, we prove that graph spectral filters are transferable.


title: Fully strained epitaxial Ti_{1-x}Mg_xN(001) layers

abstract: Ti1-xMgxN(001) layers with 0.00 {\leq} x {\leq} 0.49 are deposited on MgO(001) by reactive magnetron co-sputtering from titanium and magnesium targets in 5 mTorr pure N2 at 600 {\deg}C. X-ray diffraction {\omega}-2{\theta} scans, {\omega}-rocking curves, {\phi}-scans, and high resolution reciprocal space maps show that the Ti1-xMgxN layers are rock-salt structure single crystals with a cube-on-cube epitaxial relationship with the substrates: (001)TiMgN{\parallel} (001)MgO and [100]TiMgN{\parallel} [100]MgO. Layers with thickness d = 35-58 nm are fully strained, with an in-plane lattice parameter of 4.212{\pm}0.001 {\AA} matching that of the MgO substrate, while the out-of-plane lattice parameter increases with x from 4.251 {\AA} for TiN(001) to 4.289 {\AA} for Ti0.51Mg0.49N(001). This yields a relaxed lattice parameter for Ti1-xMgxN(001) of ao = (1-x)aTiN + xaMgN - bx(1-x), where aTiN = 4.239 {\AA}, aMgN = 4.345 {\AA}, and the bowing parameter b = 0.113 {\AA}. In contrast, thicker Ti1-xMgxN(001) layers with d = 110-275 nm are partially (pure TiN) or fully (x = 0.37 and 0.39) relaxed, indicating a critical thickness for relaxation of 50-100 nm. The in-plane x-ray coherence length is large (100-400 nm) for fully strained layers with 0.00 {\leq} x {\leq} 0.45 but drops by an order of magnitude for x = 0.49 as the composition approaches the phase stability limit. It is also an order of magnitude smaller for thicker (d {\leq} 110 nm) layers, which is attributed to strain relaxation through the nucleation and growth of misfit dislocations facilitated by glide of threading dislocations.


title: Comprehensive Evaluation of Deep Learning Architectures for Prediction of DNA/RNA Sequence Binding Specificities

abstract: Motivation: Deep learning architectures have recently demonstrated their power in predicting DNA- and RNA-binding specificities. Existing methods fall into three classes: Some are based on Convolutional Neural Networks (CNNs), others use Recurrent Neural Networks (RNNs), and others rely on hybrid architectures combining CNNs and RNNs. However, based on existing studies it is still unclear which deep learning architecture is achieving the best performance. Thus an in-depth analysis and evaluation of the different methods is needed to fully evaluate their relative. Results: In this study, We present a systematic exploration of various deep learning architectures for predicting DNA- and RNA-binding specificities. For this purpose, we present deepRAM, an end-to-end deep learning tool that provides an implementation of novel and previously proposed architectures; its fully automatic model selection procedure allows us to perform a fair and unbiased comparison of deep learning architectures. We find that an architecture that uses k-mer embedding to represent the sequence, a convolutional layer and a recurrent layer, outperforms all other methods in terms of model accuracy. Our work provides guidelines that will assist the practitioner in choosing the best architecture for the task at hand, and provides some insights on the differences between the models learned by convolutional and recurrent networks. In particular, we find that although recurrent networks improve model accuracy, this comes at the expense of a loss in the interpretability of the features learned by the model. Availability and implementation: The source code for deepRAM is available at https://github.com/MedChaabane/deepRAM


title: Big Data Platform Architecture Under The Background of Financial Technology

abstract: With the rise of the concept of financial technology, financial and technology gradually in-depth integration, scientific and technological means to become financial product innovation, improve financial efficiency and reduce financial transaction costs an important driving force. In this context, the new technology platform is from the business philosophy, business model, technical means, sales, internal management, and other dimensions to re-shape the financial industry. In this paper, the existing big data platform architecture technology innovation, adding space-time data elements, combined with the insurance industry for practical analysis, put forward a meaningful product circle and customer circle.


title: Expected $f$-vector of the Poisson Zero Polytope and Random Convex Hulls in the Half-Sphere

abstract: We prove an explicit combinatorial formula for the expected number of faces of the zero polytope of the homogeneous and isotropic Poisson hyperplane tessellation in $\mathbb R^d$. The expected $f$-vector is expressed through the coefficients of the polynomial $$ (1+ (d-1)^2x^2) (1+(d-3)^2 x^2) (1+(d-5)^2 x^2) \ldots. $$ Also, we compute explicitly the expected $f$-vector and the expected volume of the spherical convex hull of $n$ random points sampled uniformly and independently from the $d$-dimensional half-sphere. In the case when $n=d+2$, we compute the probability that this spherical convex hull is a spherical simplex, thus solving an analogue of the Sylvester four-point problem on the half-sphere.


title: Timing Properties of Shocked Accretion Flows around Neutron Stars in Presence of Cooling

abstract: We carry out the first robust numerical simulation of accretion flows on a weakly magnetized neutron star using Smoothed Particle Hydrodynamics (SPH). We follow the Two-Component Advective Flow (TCAF) paradigm for black holes, and focus only on the advective component for the case of a neutron star. This low viscosity sub-Keplerian flow will create a normal boundary layer (or, NBOL) right on the star surface in addition to the centrifugal pressure supported boundary layer (or, CENBOL) present in a black hole accretion. These density jumps could give rise to standing or oscillating shock fronts. During a hard spectral state, the incoming flow has a negligible viscosity causing more sub-Keplerian component as compared to the Keplerian disc component. We show that our simulation of flows with a cooling and a negligible viscosity produces precisely two shocks and strong supersonic winds from these boundary layers. We find that the specific angular momentum of matter dictates the locations and the nature of oscillations of these shocks. For low angular momentum flows, the radial oscillation appears to be preferred. For flows with higher angular momentum, the vertical oscillation appears to become dominant. In all the cases, asymmetries w.r.t. the Z=0 plane are seen and instabilities set in due to the interaction of inflow and outgoing strong winds. Our results capture both the low and high-frequency quasi-periodic oscillations without invoking magnetic fields or any precession mechanism. Most importantly, these solutions directly corroborate observed features of wind dominated high-mass X-ray binaries, such as Cir X-1.


title: An Advanced Conceptual Diagnostic Healthcare Framework for Diabetes and Cardiovascular Disorders

abstract: The data mining along with emerging computing techniques have astonishingly influenced the healthcare industry. Researchers have used different Data Mining and Internet of Things (IoT) for enrooting a programmed solution for diabetes and heart patients. However, still, more advanced and united solution is needed that can offer a therapeutic opinion to individual diabetic and cardio patients. Therefore, here, a smart data mining and IoT (SMDIoT) based advanced healthcare system for proficient diabetes and cardiovascular diseases have been proposed. The hybridization of data mining and IoT with other emerging computing techniques is supposed to give an effective and economical solution to diabetes and cardio patients. SMDIoT hybridized the ideas of data mining, Internet of Things, chatbots, contextual entity search (CES), bio-sensors, semantic analysis and granular computing (GC). The bio-sensors of the proposed system assist in getting the current and precise status of the concerned patients so that in case of an emergency, the needful medical assistance can be provided. The novelty lies in the hybrid framework and the adequate support of chatbots, granular computing, context entity search and semantic analysis. The practical implementation of this system is very challenging and costly. However, it appears to be more operative and economical solution for diabetes and cardio patients.


title: Dependence of Convective Boundary Mixing on Boundary Properties and Turbulence Strength

abstract: Convective boundary mixing is one of the major uncertainties in stellar evolution. In order to study its dependence on boundary properties and turbulence strength in a controlled way, we computed a series of 3D hydrodynamical simulations of stellar convection during carbon burning with a varying boosting factor of the driving luminosity. Our 3D implicit large eddy simulations were computed with the PROMPI code. We performed a mean field analysis of the simulations within the Reynolds-averaged Navier-Stokes framework. Both the vertical RMS velocity within the convective region and the bulk Richardson number of the boundaries are found to scale with the driving luminosity as expected from theory. The positions of the convective boundaries were estimated through the composition profiles across them, and the strength of convective boundary mixing was determined by analysing the boundaries within the framework of the entrainment law. We find that the entrainment is approximately inversely proportional to the bulk Richardson number. Although the entrainment law does not encompass all the processes occurring at boundaries, our results support the use of the entrainment law to describe convective boundary mixing in 1D models, at least for the advanced phases. The next steps and challenges ahead are also discussed.


title: Analytical solution of the incompressible Navier-Stokes equations

abstract: The Navier-Stokes equations, which govern fluid motions, are not resolved yet. This investigation relates to the application of the power series method to the incompressible Navier-Stokes equations. This method involves replacing variables by power series in the equations. The resulting relations allow the computation of the series coefficients from the initial and boundary conditions. All the steps of applying the method to the Navier-Stokes equations are detailed and the exact analytical solution is then established. The last part discusses the domain of convergence of the solution.


title: Voltage Profile Improvement of Distribution Grid by Using a New Control Approach on Injected Reactive Power of Plug-in Electric Vehicle Parking Lots to Grid

abstract: Establishing electrical parking lots will become more important by increasing the use of electric vehicles (EVs). These parking lots not only can be seen as high electrical power consumption loads which can cause the voltage drop at feeders but also they can be used as electrical power plants that can help the main power grid during the load peak hour or any congestion. Therefore, finding an optimum place for the building of these lots, in which the deviation of the voltage at the feeder, power loss in the grid, and cost be minimized, is essential. In this paper, a new method of vehicle-to-grid reactive power support (V2GQ) has been used to add the model of plug-in electric vehicle (PEV) as one part of the objective function to find the optimum place of parking lots. The non-sorting genetic algorithm II (NSGA2) is used here as an optimization algorithm to find the optimum voltage profile based on the location of parking lots. For validation of the purposed method, a 33-bus standard distribution network has been studied.


title: Investigating Limit Order Book Characteristics for Short Term Price Prediction: a Machine Learning Approach

abstract: With the proliferation of algorithmic high-frequency trading in financial markets, the Limit Order Book has generated increased research interest. Research is still at an early stage and there is much we do not understand about the dynamics of Limit Order Books. In this paper, we employ a machine learning approach to investigate Limit Order Book features and their potential to predict short term price movements. This is an initial broad-based investigation that results in some novel observations about LOB dynamics and identifies several promising directions for further research. Furthermore, we obtain prediction results that are significantly superior to a baseline predictor.


title: Generalized Coordinated Multipoint (GCoMP)-Enabled NOMA: Outage, Capacity, and Power Allocation

abstract: A novel generalized coordinated multi-point transmission (GCoMP)-enabled non-orthogonal multiple access (NOMA) scheme is proposed. In particular, the traditional joint transmission CoMP scheme is generalized to be applied for all user-equipments (UEs), i.e. both cell-centre and cell-edge users within the coverage area of cellular base stations (BSs). Furthermore, every BS applies NOMA for all UEs associated to it using the same frequency sub-band (i.e. all UEs associated to a BS forms a single NOMA cluster). To evaluate the proposed scheme, we derive a closed-form expression for the probability of outage for a UE with different orders of BS cooperation. Important insights on the proposed system are extracted by deriving an approximate (asymptotic) expressions for the probability of outage and outage capacity. Furthermore, an optimal transmission power allocation scheme that jointly allocates transmission power fractions from all cooperating BSs to all connected UEs is developed and investigated for the proposed system. Findings show that NOMA with a large number of UEs is feasible when the GCoMP technique is used over all UEs within the network coverage area. Also, the performance degradation caused by a large NOMA cluster size is significantly mitigated by increasing the number of cooperating BSs. In addition, for given feasible system parameters and a given NOMA cluster, the lower the available power budget, the higher is the number of BSs that apply NOMA for their cluster members and the lower the number of BSs that use water-filling for power allocation.


title: Detection of non-thermal hard X-ray emission from the "Fermi bubble" in an external galaxy

abstract: We report new Chandra hard X-ray ($>2\rm~keV$) and JVLA C-band observations of the nuclear superbubble of NGC 3079, an analog of the "Fermi bubble" in our Milky Way. We detect extended hard X-ray emission on the SW side of the galactic nucleus with coherent multi-wavelength features in radio, H$\alpha$, and soft X-ray. The hard X-ray feature has a cone shape with possibly a weak cap, forming a bubble-like structure with a diameter of $\sim1.1\rm~kpc$. A similar extended feature, however, is not detected on the NE side, which is brighter in all other wavelengths such as radio, H$\alpha$, and soft X-ray. Scattered photons from the nuclear region or other nearby point-like X-ray bright sources, inverse Compton emission from cosmic ray electrons via interaction with the cosmic microwave background, or any individually faint stellar X-ray source populations, cannot explain the extended hard X-ray emission on the SW side and the strongly NE/SW asymmetry. A synchrotron emission model, plus a thermal component accounting for the excess at $\sim1\rm~keV$, can well characterize the broadband radio/hard X-ray spectra. The broadband synchrotron spectra do not show any significant cutoff, and even possibly slightly flatten at higher energy. This rules out a loss-limited scenario in the acceleration of the cosmic ray electrons in or around this superbubble. As the first detection of kpc-scale extended hard X-ray emission associated with a galactic nuclear superbubble, the spatial and spectral properties of the multi-wavelength emissions indicate that the cosmic ray leptons responsible for the broad-band synchrotron emission from the SW bubble must be accelerated in situ, instead of transported from the nuclear region of the galaxy.


title: Finitely many physical measures for sectional-hyperbolic attracting sets and statistical stability

abstract: We show that a sectional-hyperbolic attracting set for a H\"older-$C^1$ vector field admits finitely many physical/SRB measures whose ergodic basins cover Lebesgue almost all points of the basin of topological attraction. In addition, these physical measures depend continuously on the flow in the $C^1$ topology, that is, sectional-hyperbolic attracting sets are statistically stable. To prove these results we show that each central-unstable disk in a neighborhood of this class of attracting sets is eventually expanded to contain a ball whose inner radius is uniformly bounded away from zero.


title: Holographic insulator/superconductor phase transition in higher dimensional Gauss-Bonnet gravity

abstract: In this work, we have analytically investigated the insulator/superconductor phase transition in the presence of $d$-dimensional Gauss-Bonnet AdS soliton background. Using the Sturm-Liouville eigenvalue method, we have calculated the value of the critical chemical potential $\mu_c$ in any arbitrary dimension $d\geq 5$. We have then studied the condensation operator values and charge density in terms of the chemical potential and discussed the $d=5, 6, 7$ cases using our general results in $d$ dimensions. Our analytical results agree very well with the numerically findings in the literature.


title: Impact of Data Pruning on Machine Learning Algorithm Performance

abstract: Dataset pruning is the process of removing sub-optimal tuples from a dataset to improve the learning of a machine learning model. In this paper, we compared the performance of different algorithms, first on an unpruned dataset and then on an iteratively pruned dataset. The goal was to understand whether an algorithm (say A) on an unpruned dataset performs better than another algorithm (say B), will algorithm B perform better on the pruned data or vice-versa. The dataset chosen for our analysis is a subset of the largest movie ratings database publicly available on the internet, IMDb [1]. The learning objective of the model was to predict the categorical rating of a movie among 5 bins: poor, average, good, very good, excellent. The results indicated that an algorithm that performed better on an unpruned dataset also performed better on a pruned dataset.


title: One-dimensional backreacting holographic p-wave superconductors

abstract: We analytically as well as numerically study the properties of one-dimensional holographic p- wave superconductors in the presence of backreaction. We employ the Sturm-Liouville eigenvalue problem for the analytical calculation and the shooting method for the numerical investigations. We apply the AdS3/CFT2 correspondence and determine the relation between the critical temperature Tc and the chemical potential \mu for different values of mass m of charged spin-1 field and backreacting parameters. We observe that the data of both analytical and numerical studies are in good agreement. We find out that increasing the backreaction as well as the mass parameter cause the greater values for Tc/ \mu. Therefore, it makes the condensation harder to form. In addition, the analytical and numerical approaches show that the value of the critical exponent \beta is 1/2 which is the same as in the mean field theory. Moreover, both methods confirm the exhibition of a second order phase transition.


title: Abstract I/O Specification

abstract: We recently proposed an approach for the specification and modular formal verification of the interactive (I/O) behavior of programs, based on an embedding of Petri nets into separation logic. While this approach is scalable and modular in terms of the I/O APIs available to a program, enables composing low-level I/O actions into high-level ones, and enables a convenient verification experience, it does not support high-level I/O actions that involve memory manipulation as well as low-level I/O (such as buffered I/O), or that are in fact "virtual I/O" actions that are implemented purely through memory manipulation. Furthermore, it does not allow rewriting an I/O specification into an equivalent one. In this paper, we propose a refined approach that does have these properties. The essential insight is to fix the set of places of the Petri net to be the set of separation logic assertions, thus making available the full power of separation logic for abstractly stating an arbitrary operation's specification in Petri net form, for composing operations into an I/O specification, and for equivalence reasoning on I/O specifications. Our refinement resolves the issue of the justification of the choice of Petri nets over other formalisms such as general state transition systems, in that it "refines them away" into the more essential constructs of separating conjunction and abstract nested triples. To enable a convenient treatment of input operations, we propose the use of prophecy variables to eliminate their non-determinism. We illustrate the approach through a number of example programs, including one where subroutines specified and verified using I/O specifications run as threads communicating through shared memory. The theory and examples of the paper have been machine-checked using the Iris library for program verification in the Coq proof assistant.


title: Renormalization of determinant lines in Quantum Field Theory

abstract: On a compact manifold $M$, we consider the affine space $A$ of non self-adjoint perturbations of some invertible elliptic operator acting on sections of some Hermitian bundle, by some differential operator of lower order. We construct and classify all complex analytic functions on the Fr\'echet space $A$ vanishing exactly over non invertible elements, having minimal order and which are obtained by local renormalizations, a concept coming from quantum field theory, called renormalized determinants. The additive group of local polynomial functionals of finite degrees acts freely and transitively on the space of renormalized determinants. We provide different representations of the renormalized determinants in terms of spectral zeta determinants, Gaussian Free Fields, infinite product and renormalized Feynman amplitudes in perturbation theory in position space \`a la Epstein-Glaser. Specializing to the case of Dirac operators coupled to vector potentials and reformulating our results in terms of determinant line bundles, we prove our renormalized determinants define some complex analytic trivializations of some holomorphic line bundle over $A$ relating our results to a conjectural picture from some unpublished notes by Quillen [52] from April 1989.


title: A High-Dimensional Particle Filter Algorithm

abstract: Online data assimilation in time series models over a large spatial extent is an important problem in both geosciences and robotics. Such models are intrinsically high-dimensional, rendering traditional particle filter algorithms ineffective. Though methods that begin to address this problem exist, they either rely on additional assumptions or lead to error that is spatially inhomogeneous. I present a novel particle-based algorithm for online approximation of the filtering problem on such models, using the fact that each locus affects only nearby loci at the next time step. The algorithm is based on a Metropolis-Hastings-like MCMC for creating hybrid particles at each step. I show simulation results that suggest the error of this algorithm is uniform in both space and time, with a lower bias, though higher variance, as compared to a previously-proposed algorithm.


title: Quantum Brownian oscillator for the stock market

abstract: We pursue the quantum-mechanical challenge to the efficient market hypothesis for the stock market by employing the quantum Brownian motion model. We utilize the quantum Caldeira-Leggett master equation as a possible phenomenological model for the stock-market-prices fluctuations while introducing the external harmonic field for the Brownian particle. Two quantum regimes are of particular interest: the exact regime as well as the approximate regime of the pure decoherence ("recoilless") limit of the Caldeira-Leggett equation. By calculating the standard deviation and the kurtosis for the particle's position observable, we can detect deviations of the quantum-mechanical behavior from the classical counterpart, which bases the efficient market hypothesis. By varying the damping factor, temperature as well as the oscillator's frequency, we are able to provide interpretation of different economic scenarios and possible situations that are not normally recognized by the efficient market hypothesis. Hence we recognize the quantum Brownian oscillator as a possibly useful model for the realistic behavior of stock prices.


title: Modeling and inference of spatio-temporal protein dynamics across brain networks

abstract: Models of misfolded proteins (MP) aim at discovering the bio-mechanical propagation properties of neurological diseases (ND) by identifying plausible associated dynamical systems. Solving these systems along the full disease trajectory is usually challenging, due to the lack of a well defined time axis for the pathology. This issue is addressed by disease progression models (DPM) where long-term progression trajectories are estimated via time reparametrization of individual observations. However, due to their loose assumptions on the dynamics, DPM do not provide insights on the bio-mechanical properties of MP propagation. Here we propose a unified model of spatio-temporal protein dynamics based on the joint estimation of long-term MP dynamics and time reparameterization of individuals observations. The model is expressed within a Gaussian Process (GP) regression setting, where constraints on the MP dynamics are imposed through non--linear dynamical systems. We use stochastic variational inference on both GP and dynamical system parameters for scalable inference and uncertainty quantification of the trajectories. Experiments on simulated data show that our model accurately recovers prescribed rates along graph dynamics and precisely reconstructs the underlying progression. When applied to brain imaging data our model allows the bio-mechanical interpretation of amyloid deposition in Alzheimer's disease, leading to plausible simulations of MP propagation, and achieving accurate predictions of individual MP deposition in unseen data.


title: Optimization of Patterned Surfaces for Improved Superhydrophobicity Through Cost-Effective Large-Scale Computations

abstract: The growing need for creating surfaces with specific wetting properties, such as superhyrdophobic behavior, asks for novel methods for their efficient design. In this work, a fast computational method for the evaluation of patterned superhyrdophobic surfaces is introduced. The hydrophobicity of a surface is quantified in energy terms through an objective function. The increased computational cost led to the parallelization of the method with the Message Passing Interface (MPI) communication protocol that enables calculations on distributed memory systems allowing for parametric investigations at acceptable time frames. The method is demonstrated for a surface consisting of an array of pillars with inverted conical (frustum) geometry. The parallel speedup achieved allows for low cost parametric investigations on the effect of the fine features (curvature and slopes) of the pillars on the superhydophobicity of the surface and consequently for the optimization of superhyrdophobic surfaces.


title: Low-Temperature Quantum Fokker-Planck and Smoluchowski Equations and Their Extension to Multistate Systems

abstract: Simulating electron-nucleus coupled dynamics poses a non-trivial challenge and an important problem in the investigation of ultrafast processes involving coupled electronic and vibrational dynamics. Because irreversibility of the system dynamics results from thermal activation and dissipation caused by the environment, in dynamical studies, it is necessary to include heat bath degrees of freedom in the total system. When the system dynamics involves high-energy electronic transitions, the environment is regarded to be in a low-temperature regime and we must treat it quantum mechanically. In this paper, we present rigorous and versatile approaches for investigating the dynamics of open systems with coupled electronic and vibrational degrees of freedom within a fully quantum mechanical framework. These approaches are based on a quantum Fokker-Planck equation and a quantum Smoluchowski equation employing a heat bath with an Ohmic spectral density, with non-Markovian low-temperature correction terms, and extensions of these equations to the case of multi-state systems. The accuracy of these equations was numerically examined for a single-state Brownian system, while their applicability was examined for multi-state double-well systems by comparing their results with those of the fewest-switch surface hopping and Ehrenfest methods with a classical Markovian Langevin force. Comparison of the transient absorption spectra obtained using these methods clearly reveals the importance of the quantum low-temperature correction terms. These equations allow us to treat non-adiabatic dynamics in an efficient way, while maintaining numerical accuracy. The C++ source codes that we developed, which allow for the treatment of the phase and coordinate space dynamics with any single-state or multi-state potential forms, are provided as Supporting Information.


title: Latent Normalizing Flows for Discrete Sequences

abstract: Normalizing flows are a powerful class of generative models for continuous random variables, showing both strong model flexibility and the potential for non-autoregressive generation. These benefits are also desired when modeling discrete random variables such as text, but directly applying normalizing flows to discrete sequences poses significant additional challenges. We propose a VAE-based generative model which jointly learns a normalizing flow-based distribution in the latent space and a stochastic mapping to an observed discrete space. In this setting, we find that it is crucial for the flow-based distribution to be highly multimodal. To capture this property, we propose several normalizing flow architectures to maximize model flexibility. Experiments consider common discrete sequence tasks of character-level language modeling and polyphonic music generation. Our results indicate that an autoregressive flow-based model can match the performance of a comparable autoregressive baseline, and a non-autoregressive flow-based model can improve generation speed with a penalty to performance.


title: On eigenvalues of a high-dimensional spatial-sign covariance matrix

abstract: This paper investigates limiting properties of eigenvalues of multivariate sample spatial-sign covariance matrices when both the number of variables and the sample size grow to infinity. The underlying p-variate populations are general enough to include the popular independent components model and the family of elliptical distributions. A first result of the paper establishes that the distribution of the eigenvalues converges to a deterministic limit that belongs to the family of generalized Marcenko-Pastur distributions. Furthermore, a new central limit theorem is established for a class of linear spectral statistics. We develop two applications of these results to robust statistics for a high-dimensional shape matrix. First, two statistics are proposed for testing the sphericity. Next, a spectrum-corrected estimator using the sample spatial-sign covariance matrix is proposed. Simulation experiments show that in high dimension, the sample spatial-sign covariance matrix provides a valid and robust tool for mitigating influence of outliers.


title: Personalized Treatment Selection using Causal Heterogeneity

abstract: Randomized experimentation (also known as A/B testing or bucket testing) is widely used in the internet industry to measure the metric impact obtained by different treatment variants. A/B tests identify the treatment variant showing the best performance, which then becomes the chosen or selected treatment for the entire population. However, the effect of a given treatment can differ across experimental units and a personalized approach for treatment selection can greatly improve upon the usual global selection strategy. In this work, we develop a framework for personalization through (i) estimation of heterogeneous treatment effect at either a cohort or member-level, followed by (ii) selection of optimal treatment variants for cohorts (or members) obtained through (deterministic or stochastic) constrained optimization. We perform a two-fold evaluation of our proposed methods. First, a simulation analysis is conducted to study the effect of personalized treatment selection under carefully controlled settings. This simulation illustrates the differences between the proposed methods and the suitability of each with increasing uncertainty. We also demonstrate the effectiveness of the method through a real-life example related to serving notifications at Linkedin. The solution significantly outperformed both heuristic solutions and the global treatment selection baseline leading to a sizable win on top-line metrics like member visits.


title: A Critical Review on the Assumptions of SETI

abstract: The Search for Extraterrestrial Intelligence (SETI) makes certain assumptions which guide all current search programs. To illustrate some, this includes (1) that interstellar flight is not possible (2) that the motivations of interstellar cultures are based largely on anthropomorphic understandings of homo sapiens (3) that the Fermi Paradox and the Drake equation are the starting point (axioms) of all reasoning (4) that definitions of 'life' are based largely on our understanding of homeostasis (5) that radio waves are the most likely method of interstellar communications (6) that unknown single event source signatures detected in space are not amenable to scrutiny due to the demands of the scientific method to be reproducible (7) that such anomalous signatures are either astronomical or communications based in type with no consideration for emissions from advanced industrialisation or propulsion and power technology. These assumptions, and others, have guided the SETI community towards a constrained level of thinking that is equivalent to philosophical dogma. In this paper, we unpack these assumptions, and others, and argue that the potential for life and intelligent life in the Cosmos may be much greater than the SETI community currently appears to conclude. It is also argued that more progress in our underanding of our place in the Cosmos can be made, if the separate disciplines of astronomy, interstellar spacecraft design, SETI, biology and philosophy can work together in a complementary way. Presented at the 47th IAA Symposium on the Search for Extraterrestrial Intelligence. SETI and Society.


title: Stochastic Estimated Risk for Storage Capacity

abstract: Managing data storage growth is of crucial importance to businesses. Poor practices can lead to large data and financial losses. Access to storage information along with timely action, or capacity forecasting, are essential to avoid these losses. In addition, ensuring high accuracy of capacity forecast estimates along with ease of interpretability plays an important role for any customer facing tool. In this paper, we introduce Stochastic Estimated Risk (SER), a tool developed at Nutanix that has been in production. SER shifts the focus from forecasting a single estimate for date of attaining full capacity to predicting the risk associated with running out of storage capacity. Using a Brownian motion with drift model, SER estimates the probability that a system will run out of capacity within a specific time frame. Our results showed that a probabilistic approach is more accurate and credible, for systems with non-linear patterns, compared to a regression or ensemble forecasting models.