title: An Empirical Study on the Design and Evolution of NoSQL Database Schemas

abstract: We study how software engineers design and evolve their domain model when building applications against NoSQL data stores. Specifically, we target Java projects that use object-NoSQL mappers to interface with schema-free NoSQL data stores. Given the source code of ten real-world database applications, we extract the implicit NoSQL database schema. We capture the sizes of the schemas, and investigate whether the schema is denormalized, as is recommended practice in data modeling for NoSQL data stores. Further, we analyze the entire project history, and with it, the evolution history of the NoSQL database schema. In doing so, we conduct the so far largest empirical study on NoSQL schema design and evolution.


title: Thermo-electrical modelling of the ATLAS ITk Strip Detector

abstract: In this paper we discuss the use of linked thermal and electrical network models to predict the behaviour of a complex silicon detector system. We use the silicon strip detector for the ATLAS Phase-II upgrade to demonstrate the application of such a model and its performance. With this example, a thermo-electrical model is used to test design choices, validate specifications, predict key operational parameters such as cooling system requirements, and optimize operational aspects like the temperature profile over the lifetime of the experiment. The model can reveal insights into the interplay of conditions and components in the silicon module, and it is a valuable tool for estimating the headroom to thermal runaway, all with very moderate computational effort.


title: Measuring Node Contribution to Community Structure with Modularity Vitality

abstract: Community-aware centrality is an emerging research area in network science concerned with the importance of nodes in relation to community structure. Measures are a function of a network's structure and a given partition. Previous approaches extend classical centrality measures to account for community structure with little connection to community detection theory. In contrast, we propose cluster-quality vitality measures, i.e., modularity vitality, a community-aware measure which is well-grounded in both centrality and community detection theory. Modularity vitality quantifies positive and negative contributions to community structure, which indicate a node's role as a community bridge or hub. We derive a computationally efficient method of calculating modularity vitality for all nodes in O(M + NC) time, where C is the number of communities. We systematically fragment networks by removing central nodes, and find that modularity vitality consistently outperforms existing community-aware centrality measures. Modularity vitality is over 8 times more effective than the next-best method on a million-node infrastructure network. This result does not generalize to social media communication networks, which exhibit extreme robustness to all community-aware centrality attacks. This robustness suggests that user-based interventions to mitigate misinformation diffusion will be ineffective. Finally, we demonstrate that modularity vitality provides a new approach to community-deception.


title: Local Well-posedness for the Magnetohydrodynamics in the different two liquids case

abstract: We prove the local well-posedness for a two phase problem of magnetohydrodynamics with a sharp interface. The solution is obtained in the maximal regularity space: $H^1_p((0, T), L_q) \cap L_p((0, T), H^2_q)$ with $1 < p, q < \infty$ and 2/p + N/q < 1$, where $N$ is the space dimension


title: Generalized Rational Variable Projection With Application in ECG Compression

abstract: In this paper we develop an adaptive transform-domain technique based on rational function systems. It is of general importance in several areas of signal theory, including filter design, transfer function approximation, system identification, control theory etc. The construction of the proposed method is discussed in the framework of a general mathematical model called variable projection. First we generalize this method by adding dimension type free parameters. Then we deal with the optimization problem of the free parameters. To this order, based on the well-known particle swarm optimization (PSO) algorithm, we develop the multi-dimensional hyperbolic PSO algorithm. It is designed especially for the rational transforms in question. As a result, the system along with its dimension is dynamically optimized during the process. The main motivation was to increase the adaptivity while keeping the computational complexity manageable. We note that the proposed method is of general nature. As a case study the problem of electrocardiogram (ECG) signal compression is discussed. By means of comparison tests performed on the PhysioNet MIT-BIH Arrhythmia database we demonstrate that our method outperforms other transformation techniques.


title: Deterministic Intra-Vehicle Communications: Timing and Synchronization

abstract: As we power through to the future, in-vehicle communications reliance on speed is becoming a challenging predicament. This is mainly due to the ever-increasing number of electronic control units (ECUs), which will continue to drain network capacity, hence further increasing bandwidth demand. For a wired network, a tradeoff between bandwidth requirement, reliability, and cost-effectiveness has been our main motivation in developing a high-speed network architecture that is based on the integration of two time-triggered protocols namely; Time Triggered Ethernet (TT-E) and Time Triggered Controller Area Network (TT-CAN). Therefore, as a visible example of an Internet of Vehicles technology, we present a time triggered communication-based network architecture. The new architecture can provide scalable integration of advanced functionalities, while maintaining safety and high reliability. To comply with the bandwidth requirement, we consider high-speed TT-Ethernet as the main bus (i.e., backbone network) where sub-networks can use more cost-effective and lower bandwidth TT-CAN to communicate with other entities in the network via a gateway. The main challenge in the proposed network architecture has been to resolve interoperability between two entirely different time-triggered protocols, especially in terms of timing and synchronization. In this paper, we first explore the main key drivers of the proposed architecture, which are bandwidth, reliability, and timeliness. We then demonstrate the effectiveness of our gateway design in providing full interoperability between the two time-triggered protocols.


title: Automatically matching topographical measurements of cartridge cases using a record linkage framework

abstract: Firing a gun leaves marks on cartridge cases which purportedly uniquely identify the gun. Firearms examiners typically use a visual examination to evaluate if two cartridge cases were fired from the same gun, and this is a subjective process that has come under scrutiny. Matching can be done in a more reproducible manner using automated algorithms. In this paper, we develop methodology to compare topographical measurements of cartridge cases. We demonstrate the use of a record linkage framework in this context. We compare performance using topographical measurements to older reflectance microscopy images, investigating the extent to which the former produce more accurate comparisons. Using a diverse collection of images of over 1,100 cartridge cases, we find that overall performance is generally improved using topographical data. Some subsets of the data achieve almost perfect predictive performance in terms of precision and recall, while some produce extremely poor performance. Further work needs to be done to assess if examiners face similar difficulties on certain gun and ammunition combinations. For automatic methods, a fuller investigation into their fairness and robustness is necessary before they can be deployed in practice.


title: Spin-polarized two-dimensional electron/hole gases on LiCoO$_2$ layers

abstract: First-principles calculations show the formation of a 2D spin polarized electron (hole) gas on the Li (CoO$_2$) terminated surfaces of finite slabs down to a monolayer of LiCoO$_2$ in remarkable contrast with the bulk band structure stabilized by Li donating its electron to the CoO$_2$ layer forming a Co-$d-t_{2g}^6$ insulator. By mapping the first-principles computational results to a minimal tight-binding models corresponding to a non-chiral 3D generalization of the quadripartite Su-Schriefer-Heeger (SSH4) model, we show that these surface states have topological origin.

categories: cond-mat.str-el cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.other


title: Toward a Schurification of Parking Function Formulas via bijections with Young Tableaux

abstract: This paper contains a partial answer to the open problem 3.11 of \cite{[H2008]}. That is to find an explicit bijection on Schr\"oder paths that inverts the statistics area and bounce. This paper started as an attempt to write the sum over $m$-Schr\"oder paths with a fix number of diagonal steps into Schur functions in the variables $q$ and $t$. Some results have been generalized to parking functions, and some bijections were made with standard Young tableaux giving way to partial combinatorial formulas in the basis $s_\mu(q,t)s_\lambda(X)$ for $\nabla(e_n)$ (respectively, $\nabla^m(e_n)$), when $\mu$ and $\lambda$ are hooks (respectively, $\mu$ is of length one). We also give an explicit algorithm that gives all the Schr\"oder paths related to a Schur function $s_\mu(q,t)$ when $\mu$ is of length one. In a sense, it is a partial decomposition of Schr\"oder paths into crystals.


title: Bio-Inspired Modality Fusion for Active Speaker Detection

abstract: Human beings have developed fantastic abilities to integrate information from various sensory sources exploring their inherent complementarity. Perceptual capabilities are therefore heightened, enabling, for instance, the well-known "cocktail party" and McGurk effects, i.e., speech disambiguation from a panoply of sound signals. This fusion ability is also key in refining the perception of sound source location, as in distinguishing whose voice is being heard in a group conversation. Furthermore, neuroscience has successfully identified the superior colliculus region in the brain as the one responsible for this modality fusion, with a handful of biological models having been proposed to approach its underlying neurophysiological process. Deriving inspiration from one of these models, this paper presents a methodology for effectively fusing correlated auditory and visual information for active speaker detection. Such an ability can have a wide range of applications, from teleconferencing systems to social robotics. The detection approach initially routes auditory and visual information through two specialized neural network structures. The resulting embeddings are fused via a novel layer based on the superior colliculus, whose topological structure emulates spatial neuron cross-mapping of unimodal perceptual fields. The validation process employed two publicly available datasets, with achieved results confirming and greatly surpassing initial expectations.


title: Quasilinear parabolic equations with first order terms and $L^1$-data in moving domains

abstract: The global existence of weak solutions to a class of quasilinear parabolic equations with nonlinearities depending on first order terms and integrable data in a moving domain is investigated. The class includes the $p$-Laplace equation as a special case. Weak solutions are shown to be global by obtaining appropriate estimates on the gradient as well as a suitable version of Aubin-Lions lemma in moving domains.


title: Gravitational waves from Population III binary black holes formed by dynamical capture

abstract: We use cosmological hydrodynamic simulations to study the gravitational wave (GW) signals from high-redshift binary black holes (BBHs) formed by dynamical capture (ex-situ formation channel). We in particular focus on BHs originating from the first generation of massive, metal-poor, so-called Population~III (Pop~III) stars. An alternative (in-situ) formation pathway arises in Pop~III binary stars, whose GW signature has been intensively studied. In our optimistic model, we predict a local GW event rate density for ex-situ BBHs (formed at z>4) of ~0.04 /(yr Gpc^3). This is comparable to or even higher than the conservative predictions of the rate density for in-situ BBHs ~0.01-0.1 /(yr Gpc^3), indicating that the ex-situ formation channel may be as important as the in-situ one for producing GW events. We also evaluate the detectability of our simulated GW events for selected planned GW instruments, such as the Einstein Telescope (ET). For instance, we find the all-sky detection rate with signal-to-noise ratios above 10 to be <100 per year for the xylophone configuration of ET. However, our results are highly sensitive to the sub-grid models for BBH identification and evolution, such that the GW event efficiency (rate) is reduced by a factor of 4 (20) in the pessimistic case. The ex-situ channel of Pop III BBHs deserves further investigation with better modeling of the environments around Pop III-seeded BHs.


title: A review on rigorous derivation of reduced models for fluid-structure interaction systems

abstract: In this paper we review and systematize the mathematical theory on justification of sixth-order thin-film equations as reduced models for various fluid-structure interaction systems in which fluids are lubricating underneath elastic structures. Justification is based on careful examination of energy estimates, weak convergence results of solutions of the original fluid-structure interaction systems to the solution of the sixth-order thin-film equation, and quantitative error estimates which provide even strong convergence results.


title: Oscillation in the temperature profile of the large-scale circulation of turbulent convection induced by a cubic container

abstract: We present observations of oscillations in the shape of the temperature profile of the large-scale circulation (LSC) of turbulent Rayleigh-B{\'e}nard convection. Temperature measurements are broken down into Fourier moments as a function of $\theta-\theta_0$, where $\theta$ is the azimuthal angle in a horizontal plane at mid-height, and $\theta_0$ is the LSC orientation. The oscillation structure is dominated by a 3rd order sine moment and 3rd order cosine moment in a cubic cell. In contrast, these moments are not found to oscillate in a cylindrical cell. This geometry-dependent behavior can be explained by a model that assumes that the heat transported by the LSC is conducted from the thermal boundary layers, and is proportional to pathlength of the LSC along boundary layers at the top and bottom plates. In a non-circular cross-section cell, oscillations of the LSC orientation $\theta_0$ result in an oscillation in the container shape in the reference frame of the LSC, resulting in an oscillation in the pathlength of the LSC at a given $\theta-\theta_0$. In a square-cross-section cell, this model predicts the dominant 3rd order sine moment and 3rd order cosine moment with magnitudes within 50\% of measured values, when using the amplitude of the oscillation of $\theta_0$ as input. A cylindrical cell is special in that the pathlength is independent of $\theta_0$, and so these oscillating moments are not induced. In a cylindrical cell, the model reproduces the sinusoidal mean temperature profile with a sloshing oscillation dominated by the 2nd order sine moment, consistent with previous observations in that geometry.


title: A Time Domain Approach for the Exponential Stability of a Nondissipative Linearized Compressible Flow-Structure PDE System

abstract: We consider a linearized compressible flow structure interaction (FSI) PDE model with a view of analyzing the stability properties of both the compressible flow and plate solution components. In our earlier work, we gave an answer in the affirmative to question of uniform stability for finite energy solutions of said compressible flow-structure system, by means of a ``frequency domain'' approach. However, the frequency domain method of proof in that work is not ``robust'' (insofar as we can see), when one wishes to study longtime behavior of solutions of compressible flow-structure PDE models which track the appearance of the ambient state onto the boundary interface (i.e., $\kappa =1$ in flow PDE component (2)). Nor is a frequency domain approach in this earlier work availing when one wishes to consider the dynamics, in long time, of solutions to physically relevant nonlinear versions of the compressible flow-structure PDE system under present consideration (e.g., the Navier-Stokes nonlinearity in the PDE flow component, or a nonlinearity of Berger/Von Karman type in the plate equation). Accordingly, in the present work, we operate in the time domain by way of obtaining the necessary energy estimates which culminate in an alternative proof for the uniform stability of finite energy compressible flow-structure solutions. This novel time domain proof will be used in our forthcoming paper which addresses the existence of compact global attractors for the corresponding nonlinear coupled system in which the material derivative -- which incorporates the aforesaid ambient state -- of the interaction surface will be taken into account.


title: Finite-horizon Linear Quadratic Control for Networked Control Systems with non-distributed plants

abstract: An optimal control law for networked control systems with a discrete-time linear time-invariant (LTI) system as plant and networks between sensor and controller as well as between controller and actuator is proposed. This controller is designed by solving an optimization problem that is a generalization of the optimization problem used to obtain the Linear Quadratic Regulator (LQR) for deterministic discrete-time LTI systems. The networks are represented by random delays and drop outs of transmitted data packets.


title: Inexpensive surface electromyography sleeve with consistent electrode placement enables dexterous and stable prosthetic control through deep learning

abstract: The dexterity of conventional myoelectric prostheses is limited in part by the small datasets used to train the control algorithms. Variations in surface electrode positioning make it difficult to collect consistent data and to estimate motor intent reliably over time. To address these challenges, we developed an inexpensive, easy-to-don sleeve that can record robust and repeatable surface electromyography from 32 embedded monopolar electrodes. Embedded grommets are used to consistently align the sleeve with natural skin markings (e.g., moles, freckles, scars). The sleeve can be manufactured in a few hours for less than $60. Data from seven intact participants show the sleeve provides a signal-to-noise ratio of 14, a don-time under 11 seconds, and sub-centimeter precision for electrode placement. Furthermore, in a case study with one intact participant, we use the sleeve to demonstrate that neural networks can provide simultaneous and proportional control of six degrees of freedom, even 263 days after initial algorithm training. We also highlight that consistent recordings, accumulated over time to establish a large dataset, significantly improve dexterity. These results suggest that deep learning with a 74-layer neural network can substantially improve the dexterity and stability of myoelectric prosthetic control, and that deep-learning techniques can be readily instantiated and further validated through inexpensive sleeves/sockets with consistent recording locations.


title: Enriched functor categories for functor calculus

abstract: In this paper we present background results in enriched category theory and enriched model category theory necessary for developing model categories of enriched functors suitable for doing functor calculus.


title: Directed searches for continuous gravitational waves from twelve supernova remnants in data from Advanced LIGO's second observing run

abstract: We describe directed searches for continuous gravitational waves from twelve well localized non-pulsing candidate neutron stars in young supernova remnants using data from Advanced LIGO's second observing run. We assumed that each neutron star is isolated and searched a band of frequencies from 15 to 150\,Hz, consistent with frequencies expected from known young pulsars. After coherently integrating spans of data ranging from 12.0 to 55.9 days using the F-statistic and applying data-based vetoes, we found no evidence of astrophysical signals. We set upper limits on intrinsic gravitational wave amplitude in some cases stronger than 10^{-25}, generally about a factor of two better than upper limits on the same objects from Advanced LIGO's first observing run.


title: Preliminary Report: Cerebral blood flow mediates the relationship between progesterone and perceived stress symptoms among female club athletes after mild traumatic brain injury

abstract: Female athletes are severely understudied in the field of concussion research, despite higher prevalence for injuries and tendency to have longer recovery time. Hormonal fluctuations due to normal menstrual cycle (MC) or hormonal contraceptive (HC) use have been shown to impact both post-injury symptoms and neuroimaging measures, but have not been accounted for in concussion studies. In this preliminary study, we compared arterial spin labeling measured cerebral blood flow (CBF) between concussed female club athletes 3-10 days post injury (mTBI) and demographic, HC/MC matched controls (CON). We test whether CBF mediates the relationship between progesterone levels in blood and post-injury symptoms, which may be evidence for progesterone`s role in neuroprotection. We found a significant three-way relationship between progesterone, CBF and perceived stress score (PSS) in the left middle temporal gyrus. Higher progesterone was associated with lower (more normative) PSS, as well as higher (more normative) CBF. CBF mediates 100% of the relationship between progesterone and PSS (Sobel`s p-value=0.017). These findings suggest progesterone may have a neuroprotective role after concussion and highlight the importance of controlling for the effects of sex hormones in future concussion studies.


title: A note on the Erd\H{o}s-Hajnal hypergraph Ramsey problem

abstract: We show that there is an absolute constant $c>0$ such that the following holds. For every $n > 1$, there is a 5-uniform hypergraph on at least $2^{2^{cn^{1/4}}}$ vertices with independence number at most $n$, where every set of 6 vertices induces at most 3 edges. The double exponential growth rate for the number of vertices is sharp. By applying a stepping-up lemma established by the first two authors, analogous sharp results are proved for $k$-uniform hypergraphs. This answers the penultimate open case of a conjecture in Ramsey theory posed by Erd\H{o}s and Hajnal in 1972.


title: Learned Threshold Pruning

abstract: This paper presents a novel differentiable method for unstructured weight pruning of deep neural networks. Our learned-threshold pruning (LTP) method learns per-layer thresholds via gradient descent, unlike conventional methods where they are set as input. Making thresholds trainable also makes LTP computationally efficient, hence scalable to deeper networks. For example, it takes $30$ epochs for LTP to prune ResNet50 on ImageNet by a factor of $9.1$. This is in contrast to other methods that search for per-layer thresholds via a computationally intensive iterative pruning and fine-tuning process. Additionally, with a novel differentiable $L_0$ regularization, LTP is able to operate effectively on architectures with batch-normalization. This is important since $L_1$ and $L_2$ penalties lose their regularizing effect in networks with batch-normalization. Finally, LTP generates a trail of progressively sparser networks from which the desired pruned network can be picked based on sparsity and performance requirements. These features allow LTP to achieve competitive compression rates on ImageNet networks such as AlexNet ($26.4\times$ compression with $79.1\%$ Top-5 accuracy) and ResNet50 ($9.1\times$ compression with $92.0\%$ Top-5 accuracy). We also show that LTP effectively prunes modern \textit{compact} architectures, such as EfficientNet, MobileNetV2 and MixNet.


title: Predicting accreted satellite galaxy masses and accretion redshifts based on globular cluster orbits in the E-MOSAICS simulations

abstract: The ages and metallicities of globular clusters (GCs) are known to be powerful tracers of the properties of their progenitor galaxies, enabling their use in determining the merger histories of galaxies. However, while useful in separating GCs into individual accretion events, the orbits of GC groups themselves have received less attention as probes of their progenitor galaxy properties. In this work, we use simulations of galaxies and their GC systems from the E-MOSAICS project to explore how the present-day orbital properties of GCs are related to the properties of their progenitor galaxies. We find that the orbits of GCs deposited by accretion events are sensitive to the mass and merger redshift of the satellite galaxy. Earlier mergers and larger galaxy masses deposit GCs at smaller median apocentres and lower total orbital energy. The orbital properties of accreted groups of GCs can therefore be used to infer the properties of their progenitor galaxy, though there exists a degeneracy between galaxy mass and accretion time. Combining GC orbits with other tracers (GC ages, metallicities) will help to break the galaxy mass/accretion time degeneracy, enabling stronger constraints on the properties of their progenitor galaxy. In situ GCs generally orbit at lower energies (small apocentres) than accreted GCs, however they exhibit a large tail to high energies and even retrograde orbits (relative to the present-day disc), showing significant overlap with accreted GCs. Applying the results to Milky Way GCs groups suggests a merger redshift $z \sim 1.5$ for the Gaia Sausage/Enceladus and $z>2$ for the `low-energy'/Kraken group, adding further evidence that the Milky Way had two significant mergers in its past.


title: New Rank Records For Elliptic Curves Having Rational Torsion

abstract: We present rank-record breaking elliptic curves having torsion subgroups Z/2Z, Z/3Z, Z/4Z, Z/6Z, and Z/7Z.


title: Breakdown points of penalized and hybrid M-estimators of covariance

abstract: We introduce a class of hybrid M-estimators of multivariate scatter which, analogous to the popular spatial sign covariance matrix (SSCM), possess high breakdown points. We also show that the SSCM can be viewed as an extreme member of this class. Unlike the SSCM, but like the regular M-estimators of scatter, this new class of estimators takes into account the shape of the contours of the data cloud for downweighting observations.


title: Dust depleted inner disks in a large sample of transition disks through long-baseline ALMA observations

abstract: Transition disks with large inner dust cavities are thought to host massive companions. However, the disk structure inside the companion orbit and how material flows toward an actively accreting star remain unclear. We present a high resolution continuum study of inner disks in the cavities of 38 transition disks. Measurements of the dust mass from archival Atacama Large Millimeter/Submillimeter Array observations are combined with stellar properties and spectral energy distributions to assemble a detailed picture of the inner disk. An inner dust disk is detected in 18 of 38 disks in our sample. Of the 14 resolved disks, 9 are significantly misaligned with the outer disk. The near-infrared excess is uncorrelated with the mm dust mass of the inner disk. The size-luminosity correlation known for protoplanetary disks is recovered for the inner disks as well, consistent with radial drift. The inner disks are depleted in dust relative to the outer disk and their dust mass is uncorrelated with the accretion rates. This is interpreted as the result of radial drift and trapping by planets in a low $\alpha$ ($\sim 10^{-3}$) disk, or a failure of the $\alpha$-disk model to describe angular momentum transport and accretion. The only disk in our sample with confirmed planets in the gap, PDS 70, has an inner disk with a significantly larger radius and lower inferred gas-to-dust ratio than other disks in the sample. We hypothesize that these inner disk properties and the detection of planets are due to the gap having only been opened recently by young, actively accreting planets.


title: Transferring Dense Pose to Proximal Animal Classes

abstract: Recent contributions have demonstrated that it is possible to recognize the pose of humans densely and accurately given a large dataset of poses annotated in detail. In principle, the same approach could be extended to any animal class, but the effort required for collecting new annotations for each case makes this strategy impractical, despite important applications in natural conservation, science and business. We show that, at least for proximal animal classes such as chimpanzees, it is possible to transfer the knowledge existing in dense pose recognition for humans, as well as in more general object detectors and segmenters, to the problem of dense pose recognition in other classes. We do this by (1) establishing a DensePose model for the new animal which is also geometrically aligned to humans (2) introducing a multi-head R-CNN architecture that facilitates transfer of multiple recognition tasks between classes, (3) finding which combination of known classes can be transferred most effectively to the new animal and (4) using self-calibrated uncertainty heads to generate pseudo-labels graded by quality for training a model for this class. We also introduce two benchmark datasets labelled in the manner of DensePose for the class chimpanzee and use them to evaluate our approach, showing excellent transfer learning performance.


title: High Rate Communication over One-Bit Quantized Channels via Deep Learning and LDPC Codes

abstract: This paper proposes a method for designing error correction codes by combining a known coding scheme with an autoencoder. Specifically, we integrate an LDPC code with a trained autoencoder to develop an error correction code for intractable nonlinear channels. The LDPC encoder shrinks the input space of the autoencoder, which enables the autoencoder to learn more easily. The proposed error correction code shows promising results for one-bit quantization, a challenging case of a nonlinear channel. Specifically, our design gives a waterfall slope bit error rate even with high order modulation formats such as 16-QAM and 64-QAM despite one-bit quantization. This gain is theoretically grounded by proving that the trained autoencoder provides approximately Gaussian distributed data to the LDPC decoder even though the received signal has non-Gaussian statistics due to the one-bit quantization.


title: A coarse-grid projection method for accelerating incompressible MHD flow simulations

abstract: Coarse grid projection (CGP) is a multiresolution technique for accelerating numerical calculations associated with a set of nonlinear evolutionary equations along with the stiff Poisson equations. In this article we use CGP for the first time to speed up incompressible magnetohydrodynamics (MHD) flow simulations. Accordingly, we solve the nonlinear advection-diffusion equation on a fine mesh, while we execute the electric potential Poisson equation on the corresponding coarsened mesh. Mapping operators connect two grids together. A pressure correction scheme is used to enforce the incompressibility constrain. The study of incompressible flow past a circular cylinder in the presence of Lorentz force is selected as a benchmark problem with a fixed Reynolds number but various Stuart numbers. We consider two different situations. First, we only apply CGP to the electric potential Poisson equation. Second, we apply CGP to the pressure Poisson equation as well. The maximum speedup factors achieved here are approximately 3 and 23 respectively for the first and second situations. For the both situations we examine the accuracy of velocity and vorticity fields as well as the lift and drag coefficients. In general, the results obtained by CGP are in an excellent to reasonable range of accuracy and are significantly consistently more accurate than when we use coarse grids for the discretization of both the advection-diffusion and electric potential Poisson equations.


title: Nonparametric Estimation in the Dynamic Bradley-Terry Model

abstract: We propose a time-varying generalization of the Bradley-Terry model that allows for nonparametric modeling of dynamic global rankings of distinct teams. We develop a novel estimator that relies on kernel smoothing to pre-process the pairwise comparisons over time and is applicable in sparse settings where the Bradley-Terry may not be fit. We obtain necessary and sufficient conditions for the existence and uniqueness of our estimator. We also derive time-varying oracle bounds for both the estimation error and the excess risk in the model-agnostic setting where the Bradley-Terry model is not necessarily the true data generating process. We thoroughly test the practical effectiveness of our model using both simulated and real world data and suggest an efficient data-driven approach for bandwidth tuning.


title: Convexity in Multivalued Harmonic Functions

abstract: We investigate variants of a Three Circles type Theorem in the context of \mathcal{Q}-valued functions. We prove some convexity inequalities related to the L^{2} growth function in the \mathcal{Q}-valued settings. Optimality of these inequalities and comparsion to the case of real valued harmonic functions is also discussed.


title: On the CLT of additive functionals of Markov chains

abstract: In this paper we study the additive functionals of Markov chains via conditioning with respect to both past and future of the chain. We shall point out new sufficient projective conditions, which assure that the variance of partial sums of n consecutive random variables of a stationary Markov chain is linear in n. The paper also addresses the central limit theorem problem and is listing several open questions.


title: Constrained Generative Adversarial Network Ensembles for Sharable Synthetic Data Generation

abstract: The sharing of medical imaging datasets between institutions, and even inside the same institution, is limited by various regulations/legal barriers. Although these limitations are necessities for protecting patient privacy and setting strict boundaries for data ownership, medical research projects that require large datasets suffer considerably as a result. Machine learning has been revolutionized with the emerging deep neural network approaches over recent years, making the data-related limitations even a larger problem as these novel techniques commonly require immense imaging datasets. This paper introduces constrained Generative Adversarial Network ensembles (cGANe) to address this problem by altering the representation of the imaging data, whereas containing the significant information, enabling the reproduction of similar research results elsewhere with the sharable data. Accordingly, a framework representing the generation of a cGANe is described, and the approach is validated for the generation of synthetic 3D brain metastatic region data from T1-weighted contrast-enhanced MRI studies. For 90% brain metastases (BM) detection sensitivity, our previously reported detection algorithm produced on average 9.12 false-positive BM detections per patient after training with the original data, whereas producing 9.53 false-positives after training with the cGANe generated synthetic data. Although the applicability of the introduced approach needs further validation studies with a range of medical imaging data types, the results suggest that the BM-detection algorithm can achieve comparable performance by using cGANe generated synthetic data. Hence, the generalization of the proposed approach for various modalities may occur in the near future.


title: CCD \emph{VI} time-series of the extremely metal poor globular cluster M92: revisiting its variable star population

abstract: We present an analysis of \emph{VI} CCD time-series photometry of the Oo II type globular cluster M92. The variable star population of the cluster is studied with the aim of revising their classifications, identifications, frequency spectra and to select indicators of the parental cluster metallicity and distance. The Fourier decomposition of RR Lyrae light curves lead to the estimation of mean [Fe/H]$_{\rm spec}=-2.20\pm 0.18$, and distance of $8.3\pm0.2$ kpc. Four new variables are reported; one RRd (V40), a multimode SX Phe (V41), a SR (V42) and one RRc (F1) that is most likely not a cluster member. The AC nature of V7 is confirmed. The double mode nature of the RRc star V11 is not confirmed and its amplitude modulations are most likely due to the Blazhko effect. Two modes are found in the known RRc variable V13. It is argued that the variable V30, previously classified as RRab is, in fact, an BL Her-type star not belonging to the cluster. Making use of the $Gaia$-DR2 proper motions, we identified 5012 stars in the field of the cluster, that are very likely cluster members, and for which we possess photometry, enabling the production of a refined Colour-Magnitude Diagram. This also allowed us to identify a few variable stars that do not belong to the cluster. The RR Lyrae pulsation modes on the HB are cleanly separated by the first overtone red edge, a common feature in all Oo II type clusters.


title: Data based two-body current contribution to neutrino-nucleus cross section

abstract: A phenomenological model of two-body current (2p2h) contribution to neutrino cross section is introduced. Predictions of the Valencia model for 2p2h are modified using recent CC0pi measurements from T2K and MINERvA experiments. Our results suggest a significant increase of the 2p2h cross section at neutrino energies bigger than 1 GeV and also a redistribution of 2p2h events as function of energy and momentum transfer. This may have a big impact on neutrino energy reconstruction in neutrino oscillation parameters.


title: Charge Density Modulation and Defect Ordering in Na$_x$MnBi$_y$ magnetic semimetal

abstract: The I-Mn-V antiferromagnet, NaMnBi, develops a very large positive magnetoresistance (MR) up to 10,000% at 2 K and 9 T in crystals showing a semiconductor-to-metal transition (SMT). In the absence of an SMT, a modest (20%) MR is achieved. Here, we show that upon cooling below the magnetic transition, a spatial modulation appears giving rise to new Bragg peaks due to charge and defect ordering in a checkerboard pattern, with two kinds of modulation vectors, $q_1$=($\frac23$, 0, 1) and $q_2$=($\frac23, \frac13, \frac12$). This constitutes a superlattice transition ($T_s$) that lowers the symmetry from the high temperature centrosymmetric P4/nmm to the non-centrosymmetric P$\overline4$m2. In crystals with a large MR, a close to room temperature $T_s$ is observed with $q_1$ appearing first, followed by $q_2$. In crystals with low MR however, $T_s$ is much lower and only $q_1$ is observed. The charge modulation and spin fluctuations may both contribute to the enhancement of MR.


title: Stability of Kerr black holes in generalized hybrid metric-Palatini gravity

abstract: It is shown that the Kerr solution exists in the generalized hybrid metric-Palatini gravity theory and that for certain choices of the function $f(R,\mathcal R)$ that characterizes the theory, the Kerr solution can be stable against perturbations on the scalar degree of freedom of the theory. We start by verifying which are the most general conditions on the function $f(R,\mathcal R)$ that allow for the general relativistic Kerr solution to also be a solution of this theory. We perform a scalar perturbation in the trace of the metric tensor, which in turn imposes a perturbation in both the Ricci and Palatini scalar curvatures. To first order in the perturbation, the equations of motion, namely the field equations and the equation that relates the Ricci and the Palatini curvature scalars, can be rewritten in terms of a fourth-order wave equation for the perturbation $\delta R$ which can be factorized into two second-order massive wave equations for the same variable. The usual ansatz and separation methods are applied and stability bounds on the effective mass of the Ricci scalar perturbation are obtained. These stability regimes are studied case by case and specific forms of the function $f(R,\mathcal R)$ that allow for a stable Kerr solution to exist within the perturbation regime studied are obtained.


title: Model-independent determination of the relative strong-phase difference between $D^0$ and $\bar{D}^0\rightarrow K^0_{S,L}\pi^+\pi^-$ and its impact on the measurement of the CKM angle $\gamma/\phi_3$

abstract: Crucial inputs for a variety of $CP$-violation studies can be determined through the analysis of pairs of quantum-entangled neutral $D$ mesons, which are produced in the decay of the $\psi(3770)$ resonance. The relative strong-phase parameters between $D^0$ and $\bar{D}^0$ in the decays $D^0\rightarrow K^0_{S,L}\pi^+\pi^-$ are studied using 2.93~${\rm fb}^{-1}$ of $e^+e^-$ annihilation data delivered by the BEPCII collider and collected by the BESIII detector at a center-of-mass energy of 3.773 GeV. Results are presented in regions of the phase space of the decay. These are the most precise measurements to date of the strong-phase parameters in $D \to K_{S,L}^0\pi^+\pi^-$ decays. Using these parameters, the associated uncertainty on the Cabibbo-Kobayashi-Maskawa angle $\gamma/\phi_3$ is expected to be between $0.7^\circ$ and $1.2^\circ$, for an analysis using the decay $B^{\pm}\rightarrow DK^{\pm}$, $D\rightarrow K^0_S\pi^+\pi^-$, where $D$ represents a superposition of $D^0$ and $\bar{D^0}$ states. This is a factor of three smaller than that achievable with previous measurements. Furthermore, these results provide valuable input for charm-mixing studies, other measurements of $CP$ violation, and the measurement of strong-phase parameters for other $D$-decay modes.


title: Comment on "Testing claims of the GW170817 binary neutron star inspiral affecting $\beta$-decay rates''

abstract: We show that recently reported radioactive decay data support an earlier observation of an unexpected signal associated with the GW170817 binary neutron star inspiral.


title: A conformal infinity approach to asymptotically $\text{AdS}_2\times S^{n-1}$ spacetimes

abstract: It is well known that the spacetime $\text{AdS}_2\times S^2$ arises as the `near horizon' geometry of the extremal Reisser-Nordstrom solution, and for that reason it has been studied in connection with the AdS/CFT correspondence. Motivated by a conjectural viewpoint of Juan Maldacena, the authors in [4] studied the rigidity of asymptotically $\text{AdS}_2\times S^2$ spacetimes satisfying the null energy condition. In this paper, we take an entirely different and more general approach to the asymptotics based on the notion of conformal infinity. This involves a natural modification of the usual notion of timelike conformal infinity for asymptotically anti-de Sitter spacetimes. As a consequence we are able to obtain a variety of new results, including similar results to those in [4] (but now allowing both higher dimensions and more than two ends) and a version of topological censorship.


title: Improved Algorithm for Min-Cuts in Distributed Networks

abstract: In this thesis, we present fast deterministic algorithm to find small cuts in distributed networks. Finding small min-cuts for a network is essential for ensuring the quality of service and reliability. Throughout this thesis, we use the CONGEST model which is a typical message passing model used to design and analyze algorithms in distributed networks. We survey various algorithmic techniques in the CONGEST model and give an overview of the recent results to find cuts. We also describe elegant graph theoretic ideas like cut spaces and cycle spaces that provide useful intuition upon which our work is built. Our contribution is a novel fast algorith to find small cuts. Our algorithm relies on a new characterization of trees and cuts introduced in this thesis. Our algorithm is built upon several new algorithmic ideas that, when coupled with our characterization of trees and cuts, help us to find the required min-cuts. Our novel techniques include a tree restricted semigroup function (TRSF), a novel sketching technique, and a layered algorithm. TRSF is defined with respect to a spanning tree and is based on a commutative semigroup. This simple yet powerful technique helps us to deterministically find min-cuts of size one (bridges) and min-cuts of size two optimally. Our sketching technique samples a small but relevant vertex set which is enough to find small min-cuts in certain cases. Our layered algorithm finds min-cuts in smaller sub-graphs pivoted by nodes at different levels in a spanning tree and uses them to make the decision about the min-cuts in the complete graph. This is interesting because it enables us to show that even for a global property like finding min-cuts, local information can be exploited in a coordinated manner.


title: Entanglement and boundary entropy in quantum spin chains with arbitrary direction of the boundary magnetic fields

abstract: We calculate the entanglement and the universal boundary entropy (BE) in the critical quantum spin chains, such as the transverse field Ising chain and the XXZ chain, with arbitrary direction of the boundary magnetic field (ADBMF). We determine the boundary universality class that an ADBMF induces. In particular, we show that the induced boundary conformal field theory (BCFT) depends on the point on the Bloch sphere where the boundary magnetic field directs. We show that the classification of the directions boils down to the simple fact that the boundary field breaks the bulk symmetry or does not. We present a procedure to estimate the universal BE, based on the finite-size corrections of the entanglement entropy, that apply to the ADBMF. To calculate the universal BE in the XXZ chain, we use the density matrix renormalization group (DMRG). The transverse field XY chain with ADBMF after Jordan-Wigner (JW) transformation is not a quadratic free fermion Hamiltonian. We map this model to a quadratic free fermion chain by introducing two extra ancillary spins coupled to the main chain at the boundaries, which makes the problem {\it{integrable}}. The eigenstates of the transverse field XY chain can be obtained by proper projection in the enlarged chain. Using this mapping, we are able to calculate the entanglement entropy of the transverse field XY chain using the usual correlation matrix technique up to relatively large sizes.


title: A Recursive Formula for Osculating Curves

abstract: Let $X$ be a smooth complex projective variety. Using a construction devised to Gathmann, we present a recursive formula for some of the Gromov-Witten invariants of $X$. We prove that, when $X$ is homogeneous, this formula gives the number of osculating rational curves at a general point of a general hypersurface of $X$. This generalizes the classical well known pairs of inflexion (asymptotic) lines for surfaces in $\mathbb{P}^{3}$ of Salmon, as well as Darboux's $27$ osculating conics.


title: Jointly Learning to Recommend and Advertise

abstract: Online recommendation and advertising are two major income channels for online recommendation platforms (e.g. e-commerce and news feed site). However, most platforms optimize recommending and advertising strategies by different teams separately via different techniques, which may lead to suboptimal overall performances. To this end, in this paper, we propose a novel two-level reinforcement learning framework to jointly optimize the recommending and advertising strategies, where the first level generates a list of recommendations to optimize user experience in the long run; then the second level inserts ads into the recommendation list that can balance the immediate advertising revenue from advertisers and the negative influence of ads on long-term user experience. To be specific, the first level tackles high combinatorial action space problem that selects a subset items from the large item space; while the second level determines three internally related tasks, i.e., (i) whether to insert an ad, and if yes, (ii) the optimal ad and (iii) the optimal location to insert. The experimental results based on real-world data demonstrate the effectiveness of the proposed framework. We have released the implementation code to ease reproductivity.


title: Landscape of coexisting excitonic states in the insulating single-layer cuprates and nickelates

abstract: We present an $\textit{ab initio}$ study of the excitonic states of a prototypical high-temperature superconductor La$_2$CuO$_4$ and compare them to the isostructural single-layer nickelate La$_2$NiO$_4$. Key difference in the low-energy electronic structure leads to very different excitonic behavior. Excitons in La$_2$CuO$_4$ are delocalized and can freely move in the CuO$_2$ plane without disturbing the antiferromagnetic order. In contrast, in La$_2$NiO$_4$ we find the low-lying excitonic states to be extremely localized, producing a nearly flat dispersion. The theoretically obtained excitonic dispersion and behavior are in excellent agreement with RIXS observations. To classify the excitons we project the electron-hole coupling onto each atomic site including the full manifold of atomic orbitals. We find the excitons to be composed of a linear combination of exciton classes, including Mott-Hubbard, $d-d$, and charge-transfer. The implication of these excitations to the high-$T_c$ pairing mechanism is also discussed.


title: $CP$ violation from $\tau$, $t$ and $b$ dimension-6 Yukawa couplings -- interplay of baryogenesis, EDM and Higgs physics

abstract: We explore the implications of the Standard Model effective field theory (SMEFT) with dimension-six terms involving the Higgs boson and third-generation fermion fields on the rate of Higgs boson production and decay into fermions, on the electric dipole moments (EDMs) of the electron, and on the baryon asymmetry of the Universe. We study the consequences of allowing these additional terms for each flavor separately and for combinations of two flavors. We find that a complex $\tau$ Yukawa coupling can account for the observed baryon asymmetry $Y_B^{\rm obs}$ within current LHC and EDM bounds. A complex $b$ ($t$) Yukawa coupling can account for $4\%$ ($2\%$) of $Y_B^{\rm obs}$, whereas a combination of the two can reach $12\%$. Combining $\tau$ with either $t$ or $b$ enlarges the viable parameter space owing to cancellations in the EDM and in either Higgs production times decay or the total Higgs width, respectively. Interestingly, in such a scenario there exists a region in parameter space where the SMEFT contributions to the electron EDM cancel and collider signal strengths are precisely SM-like, while producing sufficient baryon asymmetry. Measuring $CP$ violation in Higgs decays to $\tau$ leptons is the smoking gun for this scenario.


title: Landscape-induced spatial oscillations in population dynamics

abstract: We study the effect that disturbances in the ecological landscape exert on the spatial distribution of a population that evolves according to the nonlocal FKPP equation. Using both numerical and analytical techniques, we explore the three types of stationary profiles that can develop near abrupt spatial variations in the environmental conditions vital for population growth: sustained oscillations, decaying oscillations and exponential relaxation towards a flat profile. We relate the wavelength and decay length of the wrinkles to the shape of the interaction kernel. Since these spatial structures encode details about the interaction among individuals, we discuss how heterogeneities reveal information that would be hidden in a flat landscape.


title: Reply to Robinson and Michaud, arXiv:2002.08893

abstract: We respond to Robinson and Michaud's (RM) comment (arXiv:2002.08893) on our recent preprint arXiv:2002.06937, in which we discuss recent excesses in low-threshold dark matter searches, and offer a potential unifying dark matter interpretation. We thank RM for their feedback, which highlights the critical need for future measurements to directly calibrate plasmon charge yields for low $\sim$ 10 eV energy depositions. RM objected to our assertion that plasmons generated at energy scales below 100~eV may have a large branching fraction into phonons. As we argue below, the points raised by RM do not invalidate our primary conclusions, as they pertain to a much different energy scale than we discuss in our paper.

categories: hep-ph astro-ph.CO cond-mat.mtrl-sci hep-ex physics.ins-det


title: Loxodromic elements in big mapping class groups via the Hooper-Thurston-Veech construction

abstract: Let $S$ be an infinite-type surface and $p\in S$. We show that the Thurston-Veech construction for pseudo-Anosov elements, adapted for infinite-type surfaces, produces infinitely many loxodromic elements for the action of $Mod(S;p)$ on the loop graph $L(S;p)$ that do not leave any finite-type subsurface $S'\subset S$ invariant. Moreover, in the language of Bavard-Walker, Thurston-Veech's construction produces loxodromic elements of any weight. As a consequence of Bavard and Walker's work, any subgroup of $Mod(S;p)$ containing two "Thurston-Veech loxodromics" of different weight has an infinite-dimensional space of non-trivial quasimorphisms.


title: VAT: Asymptotic Cost Analysis for Multi-Level Key-Value Stores

abstract: Over the past years, there has been an increasing number of key-value (KV) store designs, each optimizing for a different set of requirements. Furthermore, with the advancements of storage technology the design space of KV stores has become even more complex. More recent KV-store designs target fast storage devices, such as SSDs and NVM. Most of these designs aim to reduce amplification during data reorganization by taking advantage of device characteristics. However, until today most analysis of KV-store designs is experimental and limited to specific design points. This makes it difficult to compare tradeoffs across different designs, find optimal configurations and guide future KV-store design. In this paper, we introduce the Variable Amplification- Throughput analysis (VAT) to calculate insert-path amplification and its impact on multi-level KV-store performance.We use VAT to express the behavior of several existing design points and to explore tradeoffs that are not possible or easy to measure experimentally. VAT indicates that by inserting randomness in the insert-path, KV stores can reduce amplification by more than 10x for fast storage devices. Techniques, such as key-value separation and tiering compaction, reduce amplification by 10x and 5x, respectively. Additionally, VAT predicts that the advancements in device technology towards NVM, reduces the benefits from both using key-value separation and tiering.


title: AraBERT: Transformer-based Model for Arabic Language Understanding

abstract: The Arabic language is a morphologically rich language with relatively few resources and a less explored syntax compared to English. Given these limitations, Arabic Natural Language Processing (NLP) tasks like Sentiment Analysis (SA), Named Entity Recognition (NER), and Question Answering (QA), have proven to be very challenging to tackle. Recently, with the surge of transformers based models, language-specific BERT based models have proven to be very efficient at language understanding, provided they are pre-trained on a very large corpus. Such models were able to set new standards and achieve state-of-the-art results for most NLP tasks. In this paper, we pre-trained BERT specifically for the Arabic language in the pursuit of achieving the same success that BERT did for the English language. The performance of AraBERT is compared to multilingual BERT from Google and other state-of-the-art approaches. The results showed that the newly developed AraBERT achieved state-of-the-art performance on most tested Arabic NLP tasks. The pretrained araBERT models are publicly available on https://github.com/aub-mind/arabert hoping to encourage research and applications for Arabic NLP.


title: Self-supervised Representation Learning for Ultrasound Video

abstract: Recent advances in deep learning have achieved promising performance for medical image analysis, while in most cases ground-truth annotations from human experts are necessary to train the deep model. In practice, such annotations are expensive to collect and can be scarce for medical imaging applications. Therefore, there is significant interest in learning representations from unlabelled raw data. In this paper, we propose a self-supervised learning approach to learn meaningful and transferable representations from medical imaging video without any type of human annotation. We assume that in order to learn such a representation, the model should identify anatomical structures from the unlabelled data. Therefore we force the model to address anatomy-aware tasks with free supervision from the data itself. Specifically, the model is designed to correct the order of a reshuffled video clip and at the same time predict the geometric transformation applied to the video clip. Experiments on fetal ultrasound video show that the proposed approach can effectively learn meaningful and strong representations, which transfer well to downstream tasks like standard plane detection and saliency prediction.


title: Quantum Kinetic Equation for spin polarization of massive quarks from pQCD

abstract: We review the recent progress in formulating a quantum kinetic theory for the polarization of spin-$\frac{1}{2}$ massive quarks from the leading-log order of perturbative QCD.


title: NewsStand CoronaViz: A Map Query Interface for Spatio-Temporal and Spatio-Textual Monitoring of Disease Spread

abstract: With the rapid continuing spread of COVID-19, it is clearly important to be able to track the progress of the virus over time in order to be better prepared to anticipate its emergence and spread in new regions as well as declines in its presence in regions thereby leading to or justifying "reopening" decisions. There are many applications and web sites that monitor officially released numbers of cases which are likely to be the most accurate methods for tracking the progress of the virus; however, they will not necessarily paint a complete picture. To begin filling any gaps in official reports, we have developed the NewsStand CoronaViz web application (https://coronaviz.umiacs.io) that can run on desktops and mobile devices that allows users to explore the geographic spread in discussions about the virus through analysis of keyword prevalence in geotagged news articles and tweets in relation to the real spread of the virus as measured by confirmed case numbers reported by the appropriate authorities. NewsStand CoronaViz users have access to dynamic variants of the disease-related variables corresponding to the numbers of confirmed cases, active cases, deaths, and recoveries (where they are provided) via a map query interface. It has the ability to step forward and backward in time using both a variety of temporal window sizes (day, week, month, or combinations thereof) in addition to user-defined varying spatial window sizes specified by direct manipulation actions (e.g., pan, zoom, and hover) as well as textually (e.g., by the name of the containing country, state or province, or county as well as textually-specified spatially-adjacent combinations thereof), and finally by the amount of spatio-temporally-varying news and tweet volume involving COVID-19.


title: Deep Learning in Mining Biological Data

abstract: Recent technological advancements in data acquisition tools allowed life scientists to acquire multimodal data from different biological application domains. Broadly categorized in three types (i.e., sequences, images, and signals), these data are huge in amount and complex in nature. Mining such an enormous amount of data for pattern recognition is a big challenge and requires sophisticated data-intensive machine learning techniques. Artificial neural network-based learning systems are well known for their pattern recognition capabilities and lately their deep architectures - known as deep learning (DL) - have been successfully applied to solve many complex pattern recognition problems. Highlighting the role of DL in recognizing patterns in biological data, this article provides - applications of DL to biological sequences, images, and signals data; overview of open access sources of these data; description of open source DL tools applicable on these data; and comparison of these tools from qualitative and quantitative perspectives. At the end, it outlines some open research challenges in mining biological data and puts forward a number of possible future perspectives.


title: Density Based Outlier Scoring on Kepler Data

abstract: In the present era of large scale surveys, big data presents new challenges to the discovery process for anomalous data. Such data can be indicative of systematic errors, extreme (or rare) forms of known phenomena, or most interestingly, truly novel phenomena which exhibit as-of-yet unobserved behaviors. In this work we present an outlier scoring methodology to identify and characterize the most promising unusual sources to facilitate discoveries of such anomalous data. We have developed a data mining method based on k-Nearest Neighbor distance in feature space to efficiently identify the most anomalous lightcurves. We test variations of this method including using principal components of the feature space, removing select features, the effect of the choice of k, and scoring to subset samples. We evaluate the peformance of our scoring on known object classes and find that our scoring consistently scores rare (<1000) object classes higher than common classes. We have applied scoring to all long cadence lightcurves of quarters 1 to 17 of Kepler's prime mission and present outlier scores for all 2.8 million lightcurves for the roughly 200k objects.


title: Technology dictates algorithms: Recent developments in read alignment

abstract: Massively parallel sequencing techniques have revolutionized biological and medical sciences by providing unprecedented insight into the genomes of humans, animals, and microbes. Modern sequencing platforms generate enormous amounts of genomic data in the form of nucleotide sequences or reads. Aligning reads onto reference genomes enables the identification of individual-specific genetic variants and is an essential step of the majority of genomic analysis pipelines. Aligned reads are essential for answering important biological questions, such as detecting mutations driving various human diseases and complex traits as well as identifying species present in metagenomic samples. The read alignment problem is extremely challenging due to the large size of analyzed datasets and numerous technological limitations of sequencing platforms, and researchers have developed novel bioinformatics algorithms to tackle these difficulties. Importantly, computational algorithms have evolved and diversified in accordance with technological advances, leading to todays diverse array of bioinformatics tools. Our review provides a survey of algorithmic foundations and methodologies across 107 alignment methods published between 1988 and 2020, for both short and long reads. We provide rigorous experimental evaluation of 11 read aligners to demonstrate the effect of these underlying algorithms on speed and efficiency of read aligners. We separately discuss how longer read lengths produce unique advantages and limitations to read alignment techniques. We also discuss how general alignment algorithms have been tailored to the specific needs of various domains in biology, including whole transcriptome, adaptive immune repertoire, and human microbiome studies.


title: Phonons become chiral in the pseudogap phase of cuprates

abstract: The nature of the pseudogap phase of cuprates remains a major puzzle. One of its new signatures is a large negative thermal Hall conductivity $\kappa_{\rm xy}$, which appears for dopings $p$ below the pseudogap critical doping $p^*$, but whose origin is as yet unknown. Because this large $\kappa_{\rm xy}$ is observed even in the undoped Mott insulator La$_2$CuO$_4$, it cannot come from charge carriers, these being localized at $p = 0$. Here we show that the thermal Hall conductivity of La$_2$CuO$_4$ is roughly isotropic, being nearly the same for heat transport parallel and normal to the CuO$_2$ planes, i.e. $\kappa_{\rm zy}(T) \approx \kappa_{\rm xy} (T)$. This shows that the Hall response must come from phonons, these being the only heat carriers able to move as easily normal and parallel to the planes . At $p > p^*$, in both La$_{\rm 1.6-x}$Nd$_{\rm 0.4}$Sr$_x$CuO$_4$ and La$_{\rm 1.8-x}$Eu$_{\rm 0.2}$Sr$_x$CuO$_4$ with $p = 0.24$, we observe no c-axis Hall signal, i.e. $\kappa_{\rm zy}(T) = 0$, showing that phonons have zero Hall response outside the pseudogap phase. The phonon Hall response appears immediately below $p^* = 0.23$, as confirmed by the large $\kappa_{\rm zy}(T)$ signal we find in La$_{1.6-x}$Nd$_{\rm 0.4}$Sr$_x$CuO$_4$ with $p = 0.21$. The microscopic mechanism by which phonons become chiral in cuprates remains to be identified. This mechanism must be intrinsic - from a coupling of phonons to their electronic environment - rather than extrinsic, from structural defects or impurities, as these are the same on both sides of $p^*$. This intrinsic phonon Hall effect provides a new window on quantum materials and it may explain the thermal Hall signal observed in other topologically nontrivial insulators.


title: Gluon Shadowing and Nuclear Entanglement

abstract: Relying on previous results that link entanglement entropy and parton distribution functions in deep inelastic scattering and focusing on the small Bjorken scaling region we present here indications that gluon shadowing might indeed be explained as due to a depletion of the entanglement entropy between observed and unobserved degrees of freedom per nucleon within a nucleus with respect to the free nucleon case. We apply to gluon shadowing the general Page approach to the calculation of the entanglement entropy in bipartite systems, giving physical motivations of the results.


title: Structure-Adaptive Sequential Testing for Online False Discovery Rate Control

abstract: Consider the online testing of a stream of hypotheses where a real--time decision must be made before the next data point arrives. The error rate is required to be controlled at {all} decision points. Conventional \emph{simultaneous testing rules} are no longer applicable due to the more stringent error constraints and absence of future data. Moreover, the online decision--making process may come to a halt when the total error budget, or alpha--wealth, is exhausted. This work develops a new class of structure--adaptive sequential testing (SAST) rules for online false discover rate (FDR) control. A key element in our proposal is a new alpha--investment algorithm that precisely characterizes the gains and losses in sequential decision making. SAST captures time varying structures of the data stream, learns the optimal threshold adaptively in an ongoing manner and optimizes the alpha-wealth allocation across different time periods. We present theory and numerical results to show that the proposed method is valid for online FDR control and achieves substantial power gain over existing online testing rules.


title: Application of dilaton chiral perturbation theory to $N_f=8$, ${\rm SU}(3)$ spectral data

abstract: We extend dilaton chiral perturbation theory (dChPT) to include the taste splittings in the Nambu--Goldstone sector observed in lattice simulations of near-conformal theories with staggered fermions. We then apply dChPT to a recent simulation by the LSD collaboration of the SU(3) gauge theory with 8 fermions in the fundamental representation, which is believed to exhibit near-conformal behavior in the infrared, and in which a light singlet scalar state, nearly degenerate with the pions, has been found. We find that the mesonic sector of this theory can be successfully described by dChPT, including, in particular, the mesonic taste splittings found in the simulation. We confirm that current simulations of this theory are in the "large-mass" regime.


title: Imaging with protons at MedAustron

abstract: Ion beam therapy has become a frequently applied form of cancer therapy over the last years. The advantage of ion beam therapy over conventional radiotherapy using photons is the strongly localized dose deposition, leading to a reduction of dose applied to surrounding healthy tissue. Currently, treatment planning for proton therapy is based on X-ray computed tomography, which entails certain sources of inaccuracy in alculation of the stopping power (SP). A more precise method to acquire the SP is to directly use high energy protons (or other ions such as carbon) and perform proton computed tomography (pCT). With this method, the ions are tracked prior to entering and after leaving the patient and finally their residual energy is measured at the very end. Therefore, an ion imaging demonstrator, comprising a tracking telescope made from double-sided silicon strip detectors and a range telescope as a residual energy detector, was set up. First measurements with this setup were performed at beam tests at MedAustron, a center for ion therapy and research in \mbox{Wiener Neustadt}, \mbox{Austria}. The facility provides three rooms for cancer treatment with proton beams as well as one which is dedicated to non-clinical research. This contribution describes the principle of ion imaging with proton beams in general as well as the design of the experimental setup. Moreover, first results from simulations and recent beam tests as well as ideas for future developments will be presented.


title: BigSurvSGD: Big Survival Data Analysis via Stochastic Gradient Descent

abstract: In many biomedical applications, outcome is measured as a ``time-to-event'' (eg. disease progression or death). To assess the connection between features of a patient and this outcome, it is common to assume a proportional hazards model, and fit a proportional hazards regression (or Cox regression). To fit this model, a log-concave objective function known as the ``partial likelihood'' is maximized. For moderate-sized datasets, an efficient Newton-Raphson algorithm that leverages the structure of the objective can be employed. However, in large datasets this approach has two issues: 1) The computational tricks that leverage structure can also lead to computational instability; 2) The objective does not naturally decouple: Thus, if the dataset does not fit in memory, the model can be very computationally expensive to fit. This additionally means that the objective is not directly amenable to stochastic gradient-based optimization methods. To overcome these issues, we propose a simple, new framing of proportional hazards regression: This results in an objective function that is amenable to stochastic gradient descent. We show that this simple modification allows us to efficiently fit survival models with very large datasets. This also facilitates training complex, eg. neural-network-based, models with survival data.


title: Simultaneous confidence bands for nonparametric regression with partially missing covariates

abstract: In this paper, we consider a weighted local linear estimator based on the inverse selection probability for nonparametric regression with missing covariates at random. The asymptotic distribution of the maximal deviation between the estimator and the true regression function is derived and an asymptotically accurate simultaneous confidence band is constructed. The estimator for the regression function is shown to be oracally efficient in the sense that it is uniformly indistinguishable from that when the selection probabilities are known. Finite sample performance is examined via simulation studies which support our asymptotic theory. The proposed method is demonstrated via an analysis of a data set from the Canada 2010/2011 Youth Student Survey.


title: Assuring the Integrity of Videos from Wireless-based IoT Devices using Blockchain

abstract: Blockchain technology has drawn attention fromvarious communities. The underlying consensus mechanism inBlockchain enables a myriad of applications for the integrityassurance of stored data. In this paper, we utilize Blockchaintechnology to verify the authenticity of a video captured by astreaming IoT device for forensic investigation purposes. Theproposed approach computes the hash of video frames beforethey leave the IoT device and are transferred to a remote basestation. To guarantee the transmission, we ensure that this hashis sent through a TCP-based connection. The hash is then storedon multiple nodes on a permissioned blockchain platform. Incase the video is modified, the discrepancy will be detected byinvestigating the previously stored hash on the blockchain andcomparing it with the hash of the existing frame in question.In this work, we present the prototype as proof-of-concept withexperiment results. The system has been tested on a RaspberryPi with different quality of videos to evaluate performance. Theresults show that the concept can be implemented with moderatevideo resolutions.


title: Communication-Optimal Tilings for Projective Nested Loops with Arbitrary Bounds

abstract: Reducing communication - either between levels of a memory hierarchy or between processors over a network - is a key component of performance optimization (in both time and energy) for many problems, including dense linear algebra, particle interactions, and machine learning. For these problems, which can be represented as nested-loop computations, previous tiling based approaches have been used to find both lower bounds on the communication required to execute them and optimal rearrangements, or blockings, to attain such lower bounds. However, such general approaches have typically assumed the problem sizes are large, an assumption that is often not met in practice. For instance, the classical $(\text{# arithmetic operations})/(\text{cache size})^{1/2}$ lower bound for matrix multiplication is not tight for matrix-vector multiplications, which must read in at least $O(\text{# arithmetic operations})$ words of memory; similar issues occur for almost all convolutions in machine learning applications, which use extremely small filter sizes (and therefore, loop bounds). In this paper, we provide an efficient way to both find and obtain, via an appropriate, efficiently constructible blocking, communication lower bounds and matching tilings which attain these lower bounds for nested loop programs with arbitrary loop bounds that operate on multidimensional arrays in the projective case, where the array indices are subsets of the loop indices. Our approach works on all such problems, regardless of dimensionality, size, memory access patterns, or number of arrays, and directly applies to (among other examples) matrix multiplication and similar dense linear algebra operations, tensor contractions, n-body pairwise interactions, pointwise convolutions, and fully connected layers.


title: Improving Certified Robustness via Statistical Learning with Logical Reasoning

abstract: Intensive algorithmic efforts have been made to enable the rapid improvements of certificated robustness for complex ML models recently. However, current robustness certification methods are only able to certify under a limited perturbation radius. Given that existing pure data-driven statistical approaches have reached a bottleneck, in this paper, we propose to integrate statistical ML models with knowledge (expressed as logical rules) as a reasoning component using Markov logic networks (MLN, so as to further improve the overall certified robustness. This opens new research questions about certifying the robustness of such a paradigm, especially the reasoning component (e.g., MLN). As the first step towards understanding these questions, we first prove that the computational complexity of certifying the robustness of MLN is #P-hard. Guided by this hardness result, we then derive the first certified robustness bound for MLN by carefully analyzing different model regimes. Finally, we conduct extensive experiments on five datasets including both high-dimensional images and natural language texts, and we show that the certified robustness with knowledge-based logical reasoning indeed significantly outperforms that of the state-of-the-arts.


title: Large-scale multiconfiguration Dirac-Hartree-Fock calculations for astrophysics: n=4 levels in P-like ions from Mn~XI to Ni~XIV

abstract: Using the multiconfiguration Dirac-Hartree-Fock and the relativistic configuration interaction methods, a consistent set of transition energies and radiative transition data for the lowest 546 (623, 701, 745) states of the $3p^4 3d$, $3s 3p^2 3d^2$, $3s 3p^3 4p$, $3s 3p^4$, $3s^2 3d^3$, $3s^2 3p^2 3d$, $3s^2 3p^2 4d$, $3s^2 3p^2 4s$, $3p^3 3d^2$, $3p^5$, $3s 3p 3d^3$, $3s 3p^3 3d$, $3s 3p^3 4s$, $3s^2 3p 3d^2$, %$3s^2 3p^2 4f$, $3s^2 3p^2 4p$, $3s^2 3p^3$ configurations in Mn~XI (Fe~XII, Co~XIII, Ni~XIV) is provided. The comparison between calculated excitation energies for the $n=4$ states and available experimental values for Fe XII indicate that the calculations are highly accurate, with uncertainties of only a few hundred cm$^{-1}$. Lines from these states are prominent in the soft X-rays. With the present calculations, several recent new identifications are confirmed. Other identifications involving $3p^2 4d$ levels in Fe~XII that were found questionable are discussed and a few new assignments are recommended. As some $n=4$ states of the other ions also show large discrepancies between experimental and calculated energies, we reassess their identification. The present study provides highly accurate atomic data for the $n=4$ states of P-like ions of astrophysical interest, for which experimental data are scarce.


title: A Time-dependent SIR model for COVID-19 with Undetectable Infected Persons

abstract: In this paper, we conduct mathematical and numerical analyses to address the following crucial questions for COVID-19: (Q1) Is it possible to contain COVID-19? (Q2) When will be the peak and the end of the epidemic? (Q3) How do the asymptomatic infections affect the spread of disease? (Q4) What is the ratio of the population that needs to be infected to achieve herd immunity? (Q5) How effective are the social distancing approaches? (Q6) What is the ratio of the population infected in the long run? For (Q1) and (Q2), we propose a time-dependent susceptible-infected-recovered (SIR) model that tracks 2 time series: (i) the transmission rate at time t and (ii) the recovering rate at time t. Such an approach is more adaptive than traditional static SIR models and more robust than direct estimation methods. Using the data provided by China, we show that the one-day prediction errors for the numbers of confirmed cases are almost in 3%, and the total number of confirmed cases is precisely predicted. Also, the turning point, defined as the day that the transmission rate is less than the recovering rate can be accurately predicted. After that day, the basic reproduction number $R_0$ is less than 1. For (Q3), we extend our SIR model by considering 2 types of infected persons: detectable and undetectable infected persons. Whether there is an outbreak in such a model is characterized by the spectral radius of a 2 by 2 matrix that is closely related to $R_0$. For (Q4), we show that herd immunity can be achieved after at least 1-1/$R_0$ fraction of individuals being infected. For (Q5) and (Q6), we analyze the independent cascade (IC) model for disease propagation in a configuration random graph. By relating the propagation probabilities in the IC model to the transmission rates and recovering rates in the SIR model, we show 2 approaches of social distancing that can lead to a reduction of $R_0$.


title: Assessing Energy Storage Requirements Based on Accepted Risks

abstract: This paper presents a framework for deriving the storage capacity that an electricity system requires in order to satisfy a chosen risk appetite. The framework takes as inputs user-defined event categories, parameterised by peak power-not-served, acceptable number of events per year and permitted probability of exceeding these constraints, and returns as an output the total capacity of storage that is needed. For increased model accuracy, our methodology incorporates multiple nodes with limited transfer capacities, and we provide a foresight-free dispatch policy for application to this setting. Finally, we demonstrate the chance-constrained capacity determination via application to a model of the British network.


title: Quantum Random Number Generator based on Violations of the Free CHSH-3 Inequality

abstract: We describe a protocol for generating random numbers based on the existence of quantum violations of a free Clauser-Horne-Shimony-Holt inequality, namely CHSH-3. Our method uses semidefinite programming relaxations to compute such violations. In a standard setting the CHSH-3 inequality involves two separated qutrits and compatible measurement, that is, commuting with each other, yielding the known quantum bound of $1+\sqrt{11/3} \approx 2.9149$. In our framework, $d$-dimensional quantum systems (qudits) where $d$ is not fixed a priori, and measurement operators possibly not compatible, are allowed. This loss of constraints yields a higher value for the maximum expectation of the CHSH-3 inequality. Based on such upper bound on the violation of CHSH-3, we develop a random number generator of type prepare-and-measure, but with one part. Our protocol generates a maximal entropy and its security is based, through self testing arguments, on the attainability of the maximal violation of the free CHSH-3 for quantum systems.


title: Label-free colorectal cancer screening using deep learning and spatial light interference microscopy (SLIM)

abstract: Current pathology workflow involves staining of thin tissue slices, which otherwise would be transparent, followed by manual investigation under the microscope by a trained pathologist. While the hematoxylin and eosin (H&E) stain is well-established and a cost-effective method for visualizing histology slides, its color variability across preparations and subjectivity across clinicians remain unaddressed challenges. To mitigate these challenges, recently we have demonstrated that spatial light interference microscopy (SLIM) can provide a path to intrinsic, objective markers, that are independent of preparation and human bias. Additionally, the sensitivity of SLIM to collagen fibers yields information relevant to patient outcome, which is not available in H&E. Here, we show that deep learning and SLIM can form a powerful combination for screening applications: training on 1,660 SLIM images of colon glands and validating on 144 glands, we obtained a benign vs. cancer classification accuracy of 99%. We envision that the SLIM whole slide scanner presented here paired with artificial intelligence algorithms may prove valuable as a pre-screening method, economizing the clinician's time and effort.


title: Scaling up Hybrid Probabilistic Inference with Logical and Arithmetic Constraints via Message Passing

abstract: Weighted model integration (WMI) is a very appealing framework for probabilistic inference: it allows to express the complex dependencies of real-world problems where variables are both continuous and discrete, via the language of Satisfiability Modulo Theories (SMT), as well as to compute probabilistic queries with complex logical and arithmetic constraints. Yet, existing WMI solvers are not ready to scale to these problems. They either ignore the intrinsic dependency structure of the problem at all, or they are limited to too restrictive structures. To narrow this gap, we derive a factorized formalism of WMI enabling us to devise a scalable WMI solver based on message passing, MP-WMI. Namely, MP-WMI is the first WMI solver which allows to: 1) perform exact inference on the full class of tree-structured WMI problems; 2) compute all marginal densities in linear time; 3) amortize inference inter query. Experimental results show that our solver dramatically outperforms the existing WMI solvers on a large set of benchmarks.


title: Time of arrival imaging: The proof of concept for a novel medical imaging modality

abstract: It has been shown that with the use of ultra-wideband (UWB) electromagnetic signal and time of arrival (ToA) principle, it is possible to locate medical implants given the permittivity distribution of the body. We propose a new imaging modality using the reverse process to acquire permittivity distributions as a surrogate of human anatomy. In the proposed systems, the locations of the signal source, receiver, and signal shapes are assumed to be known exactly. The measured data is recorded as the time it takes for the signal to travel from the signal source to the signal receiver. The finite-difference-time-domain (FDTD) method is used for the modeling of signal propagation within the phantom, which is used for both simulation and image reconstruction. Image reconstruction is achieved using linear regression on the training pairs, which includes randomly generated images and its corresponding arrival times generated using the FDTD approach. The linear weights of the training images are generated to minimize the difference between the arrival time of the reconstruction image and the measured arrival time. A simulation study using UWB signal with the central frequency of 300 MHz and the Shepp-Logan phantom was carried out. Ten-picosecond timing resolution is used for the simulation and image reconstruction. The quantitative difference between the arrival times of the phantom and the reconstructed image reduced with an increased iteration number. The quantitative error of the reconstructed image reached below 10% after 900 iterations, and 8.4% after 1200 iterations. With additional post-smoothing to suppress the introduced noise pattern through reconstruction, 6.5% error was achieved. In this paper, an approach that utilizes the ToA principle to achieve transmission imaging with radio waves is proposed and validated using a simulation study.


title: Algorithms for subelliptic multipliers in $\mathbb{C}^2$

abstract: We give examples of pseudoconvex domains of finite type in $\mathbb{C}^2$ where the Kohn algorithm for subelliptic estimates fails to yield an effective lower bound for the order of subellipticity in terms of the type. We show how to modify the algorithm to obtain an effective procedure to prove subellipticity on domains of finite type in $\mathbb{C}^2$ with real analytic boundary satisfying a condition slightly stronger than pseudoconvexity. We close with a generalization to higher dimensions.


title: Determination of Latent Dimensionality in International Trade Flow

abstract: Currently, high-dimensional data is ubiquitous in data science, which necessitates the development of techniques to decompose and interpret such multidimensional (aka tensor) datasets. Finding a low dimensional representation of the data, that is, its inherent structure, is one of the approaches that can serve to understand the dynamics of low dimensional latent features hidden in the data. Nonnegative RESCAL is one such technique, particularly well suited to analyze self-relational data, such as dynamic networks found in international trade flows. Nonnegative RESCAL computes a low dimensional tensor representation by finding the latent space containing multiple modalities. Estimating the dimensionality of this latent space is crucial for extracting meaningful latent features. Here, to determine the dimensionality of the latent space with nonnegative RESCAL, we propose a latent dimension determination method which is based on clustering of the solutions of multiple realizations of nonnegative RESCAL decompositions. We demonstrate the performance of our model selection method on synthetic data and then we apply our method to decompose a network of international trade flows data from International Monetary Fund and validate the resulting features against empirical facts from economic literature.


title: Transformers for Limit Order Books

abstract: We introduce a new deep learning architecture for predicting price movements from limit order books. This architecture uses a causal convolutional network for feature extraction in combination with masked self-attention to update features based on relevant contextual information. This architecture is shown to significantly outperform existing architectures such as those using convolutional networks (CNN) and Long-Short Term Memory (LSTM) establishing a new state-of-the-art benchmark for the FI-2010 dataset.


title: Generalized correlation based Imaging for satellites

abstract: We consider imaging of fast moving small objects in space, such as low earth orbit satellites. The imaging system consists of ground based, asynchronous sources of radiation and several passive receivers above the dense atmosphere. We use the cross correlation of the received signals to reduce distortions from ambient medium fluctuations. Imaging with correlations also has the advantage of not requiring any knowledge about the probing pulse and depends weakly on the emitter positions. We account for the target's orbital velocity by introducing the necessary Doppler compensation. We show that over limited imaging regions, a constant Doppler factor can be used, resulting in an efficient data structure for the correlations of the recorded signals. We then investigate and analyze different imaging methods using the cross-correlation data structure. Specifically, we show that using a generalized two point migration of the cross correlation data, the top eigenvector of the migrated data matrix provides superior image resolution compared to the usual single-point migration scheme. We carry out a theoretical analysis that illustrates the role of the two point migration methods as well as that of the inverse aperture in improving resolution. Extensive numerical simulations support the theoretical results and assess the scope of the imaging methodology.


title: Doomsday: A Response to Simpson's Second Question

abstract: The Doomsday Argument (DA) has sparked a variety of opinions. Here I address a key question posed by F. Simpson (2016) that confronts the views of DA proponents and those who, like me, oppose the DA. I agree that typical locations within a complete spatial distribution are calculable using ordinary frequentist probability. But I argue that the temporal probability distribution is unknown: we have records of our past yet are ignorant of our future. It is this asymmetry that upsets the idea of Copernicanism in time. Although frequentist methods do not apply to this asymmetric situation, Bayesian methods do apply. They show that the various Quick Doom and Distant Doom scenarios are equally likely. I conclude that the DA has no predictive power whatsoever.


title: Next-Generation Topology of D-Wave Quantum Processors

abstract: This paper presents an overview of the topology of D-Wave's next-generation quantum processors. It provides examples of minor embeddings and discusses performance of embedding algorithms for the new topology compared to the existing Chimera topology. It also presents some initial performance results for simple, standard Ising model classes of problems.


title: Image Hashing by Minimizing Discrete Component-wise Wasserstein Distance

abstract: Image hashing is one of the fundamental problems that demand both efficient and effective solutions for various practical scenarios. Adversarial autoencoders are shown to be able to implicitly learn a robust, locality-preserving hash function that generates balanced and high-quality hash codes. However, the existing adversarial hashing methods are inefficient to be employed for large-scale image retrieval applications. Specifically, they require an exponential number of samples to be able to generate optimal hash codes and a significantly high computational cost to train. In this paper, we show that the high sample-complexity requirement often results in sub-optimal retrieval performance of the adversarial hashing methods. To address this challenge, we propose a new adversarial-autoencoder hashing approach that has a much lower sample requirement and computational cost. Specifically, by exploiting the desired properties of the hash function in the low-dimensional, discrete space, our method efficiently estimates a better variant of Wasserstein distance by averaging a set of easy-to-compute one-dimensional Wasserstein distances. The resulting hashing approach has an order-of-magnitude better sample complexity, thus better generalization property, compared to the other adversarial hashing methods. In addition, the computational cost is significantly reduced using our approach. We conduct experiments on several real-world datasets and show that the proposed method outperforms the competing hashing methods, achieving up to 10% improvement over the current state-of-the-art image hashing methods. The code accompanying this paper is available on Github (https://github.com/khoadoan/adversarial-hashing).


title: The MAVERIC survey: A hidden pulsar and a black hole candidate in ATCA radio imaging of the globular cluster NGC 6397

abstract: Using a 16.2 hr radio observation by the Australia Telescope Compact Array (ATCA) and archival Chandra data, we found $>5\sigma$ radio counterparts to 4 known and 3 new X-ray sources within the half-light radius ($r_\mathrm{h}$) of the Galactic globular cluster NGC 6397. The previously suggested millisecond pulsar (MSP) candidate, U18, is a steep-spectrum ($S_\nu \propto \nu^\alpha$; $\alpha=-2.0^{+0.4}_{-0.5}$) radio source with a 5.5 GHz flux density of $54.7\pm 4.3~\mathrm{\mu Jy}$. We argue that U18 is most likely a "hidden" MSP that is continuously hidden by plasma shocked at the collision between the winds from the pulsar and companion star. The nondetection of radio pulsations so far is probably the result of enhanced scattering in this shocked wind. On the other hand, we observed 5.5 GHz flux of the known MSP PSR J1740-5340 (U12) to decrease by a factor of $>2.8$ during epochs of 1.4 GHz eclipse, indicating that the radio flux is absorbed in its shocked wind. If U18 is indeed a pulsar whose pulsations are scattered, we note the contrast with U12's flux decrease in eclipse, which argues for two different eclipse mechanisms at the same radio frequency. In addition to U12 and U18, we also found radio associations for 5 other Chandra X-ray sources, four of which are likely background galaxies. The last, U97, which shows strong H$\alpha$ variability, is mysterious; it may be either a quiescent black hole low-mass X-ray binary, or something more unusual.


title: Ultrastrong plasmon-phonon coupling via epsilon-near-zero nanocavities

abstract: Vibrational ultrastrong coupling (USC), where the light-matter coupling strength is comparable to the vibrational frequency of molecules, presents new opportunities to probe the interactions of molecules with zero-point fluctuations, harness cavity-enhanced chemical reactions, and develop novel devices in the mid-infrared regime. Here we use epsilon-near-zero nanocavities filled with a model polar medium (SiO$_2$) to demonstrate USC between phonons and gap plasmons. We present classical and quantum mechanical models to quantitatively describe the observed plasmon-phonon USC phenomena and demonstrate a splitting of up to 50% of the resonant frequency. Our wafer-scale nanocavity platform will enable a broad range of vibrational transitions to be harnessed for USC applications.


title: Hodge Representations

abstract: Hodge representations were introduced by Green-Griffiths-Kerr to classify the Hodge groups of polarized Hodge structures, and the corresponding Mumford-Tate subdomains of a period domain. The purpose of this article is to provide an exposition of how, given a fixed period domain $\mathcal{D}$, to enumerate the Hodge representations corresponding to Mumford-Tate subdomains $D \subset \mathcal{D}$. After reviewing the well-known classical cases that $\mathcal{D}$ is Hermitian symmetric (weight $n=1$, and weight $n=2$ with $p_g = h^{2,0}=1$), we illustrate this in the case that $\mathcal{D}$ is the period domain parameterizing polarized Hodge structures of (effective) weight two Hodge structures with first Hodge number $p_g = h^{2,0} = 2$. We also classify the Hodge representations of Calabi-Yau type, and enumerate the horizontal representations of CY 3-fold type. (The "horizontal" representations those with the property that corresponding domain $D \subset \mathcal{D}$ satisfies the infinitesimal period relation, a.k.a. Griffiths' transversality, and is therefore Hermitian.)