title: Breast Tumor Segmentation and Shape Classification in Mammograms using Generative Adversarial and Convolutional Neural Network

abstract: Mammogram inspection in search of breast tumors is a tough assignment that radiologists must carry out frequently. Therefore, image analysis methods are needed for the detection and delineation of breast masses, which portray crucial morphological information that will support reliable diagnosis. In this paper, we proposed a conditional Generative Adversarial Network (cGAN) devised to segment a breast mass within a region of interest (ROI) in a mammogram. The generative network learns to recognize the breast mass area and to create the binary mask that outlines the breast mass. In turn, the adversarial network learns to distinguish between real (ground truth) and synthetic segmentations, thus enforcing the generative network to create binary masks as realistic as possible. The cGAN works well even when the number of training samples are limited. Therefore, the proposed method outperforms several state-of-the-art approaches. This hypothesis is corroborated by diverse experiments performed on two datasets, the public INbreast and a private in-house dataset. The proposed segmentation model provides a high Dice coefficient and Intersection over Union (IoU) of 94% and 87%, respectively. In addition, a shape descriptor based on a Convolutional Neural Network (CNN) is proposed to classify the generated masks into four mass shapes: irregular, lobular, oval and round. The proposed shape descriptor was trained on Digital Database for Screening Mammography (DDSM) yielding an overall accuracy of 80%, which outperforms the current state-of-the-art.


title: Generalised Markov numbers

abstract: In this paper we introduce generalised Markov numbers and extend the classical Markov theory for the discrete Markov spectrum to the case of generalised Markov numbers. In particular we show recursive properties for these numbers and find corresponding values in the Markov spectrum. Further we construct a counterexample to the generalised Markov uniqueness conjecture. The proposed generalisation is based on geometry of numbers. It substantively uses lattice trigonometry and geometric theory of continued numbers.


title: EPIC 203868608: A low-mass quadruple star system in the Upper Scorpius OB association

abstract: Young multiple star systems provide excellent testing grounds for theories of star formation and evolution. EPIC 203868608 was previously studied~\citep{David2016} as a triple star system in the Upper Scorpius OB association, but the follow-up Keck NIRC2/HIRES/NIRSPAO observations reported here reveal its quadruple nature. We find that the system consists of a double-lined spectroscopic binary (SB2) Aab (M5+M5) and an eclipsing binary (EB) Bab with a total mass that is lower than that of the SB2. Furthermore, we measure the obliquity of the EB using the Doppler tomography technique during the primary eclipse. EPIC 203868608 Bab is likely on an inclined orbit with a projected obliquity of $-57^{+40}_{-36}$ degrees. The inclined orbit is used to constrain the tidal quality factor for low-mass stars and the evolution of the quadruple system. The analytic framework to infer obliquity that has been developed in this paper can be applied to other EB systems as well as transiting planets.


title: Continuity of attractors for $\mathcal{C}^1$ perturbations of a smooth domain

abstract: We consider a family of semilinear parabolic problems with nonlinear boundary conditions \[ \left\{ \begin{aligned} u_t(x,t) &=\Delta u(x,t) -au(x,t) + f(u(x,t)),\ x \in \Omega_\epsilon \mbox{ and } t>0\,,\\ \displaystyle\frac{\partial u}{\partial N}(x,t) &=g(u(x,t)),\ x \in \partial\Omega_\epsilon \mbox{ and } t>0\,, \end{aligned} \right. \] where $\Omega_0 \subset \mathbb{R}^n$ is a smooth (at least $\mathcal{C}^2$) domain , $\Omega_{\epsilon} = h_{\epsilon}(\Omega_0)$ and $h_{\epsilon}$ is a family of diffeomorphisms converging to the identity in the $\mathcal{C}^1$-norm. Assuming suitable regularity and dissipative conditions for the nonlinearites, we show that the problem is well posed for $\epsilon>0$ sufficiently small in a suitable scale of fractional spaces, the associated semigroup has a global attractor $\mathcal{A}_{\epsilon}$ and the family $\{\mathcal{A}_{\epsilon}\}$ is continuous at $\epsilon = 0$.


title: Galaxy detection and identification using deep learning and data augmentation

abstract: We present a method for automatic detection and classification of galaxies which includes a novel data-augmentation procedure to make trained models more robust against the data taken from different instruments and contrast-stretching functions. This method is shown as part of AstroCV, a growing open source computer vision repository for processing and analyzing big astronomical datasets, including high performance Python and C++ algorithms used in the areas of image processing and computer vision. The underlying models were trained using convolutional neural networks and deep learning techniques, which provide better results than methods based on manual feature engineering and SVMs in most of the cases where training datasets are large. The detection and classification methods were trained end-to-end using public datasets such as the Sloan Digital Sky Survey (SDSS), the Galaxy Zoo, and private datasets such as the Next Generation Virgo (NGVS) and Fornax (NGFS) surveys. Training results are strongly bound to the conversion method from raw FITS data for each band into a 3-channel color image. Therefore, we propose data augmentation for the training using 5 conversion methods. This greatly improves the overall galaxy detection and classification for images produced from different instruments, bands and data reduction procedures. The detection and classification methods were trained using the deep learning framework DARKNET and the real-time object detection system YOLO. These methods are implemented in C language and CUDA platform, and makes intensive use of graphical processing units (GPU). Using a single high-end Nvidia GPU card, it can process a SDSS image in 50 milliseconds and a DECam image in less than 3 seconds.


title: Quantum Atmospherics for Materials Diagnosis

abstract: Symmetry breaking states of matter can transmit symmetry breaking to nearby atoms or molecular complexes, perturbing their spectra. We calculate one such effect, involving the `axion electrodynamics' relevant to topological insulators, quantitatively, and identify a signature for T violating superconductivity. We provide an operator framework whereby effects of this kind can be analyzed systematically.


title: A comparison of Alpha Particle and Proton Beam Differential flow in Collisionally Young Solar Wind

abstract: In fast wind or when the local Coulomb collision frequency is low, observations show that solar wind minor ions and ion sub-populations flow with different bulk velocities. Measurements indicate that the drift speed of both alpha particles and proton beams with respect to the bulk or core protons rarely exceeds the local Alfv\'en speed, suggesting that a magnetic instability or other wave-particle process limits their maximum drift. We compare simultaneous alpha particle, proton beam, and proton core observations from instruments on the Wind spacecraft spanning over 20 years. In nearly collisionless solar wind, we find that the normalized alpha particle drift speed is slower than the normalized proton beam speed; no correlation between fluctuations in both species' drifts about their means; and a strong anti-correlation between collisional age and alpha-proton differential flow, but no such correlation with proton beam-core differential flow. Controlling for the collisional dependence, both species' normalized drifts exhibit similar statistical distributions. In the asymptotic, zero Coulomb collision limit, the youngest measured differential flows most strongly correlate with an approximation of the Alfv\'en speed that includes proton pressure anisotropy. In this limit and with this most precise representation, alpha particles drift at 67% and proton beam drift is approximately 105% of the local Alfv\'en speed. We posit that one of two physical explanations is possible. Either (1) an Alfv\'enic process preferentially accelerates or sustains proton beams and not alphas or (2) alpha particles are more susceptible to either an instability or Coulomb drag than proton beams.


title: Accelerated Reinforcement Learning for Sentence Generation by Vocabulary Prediction

abstract: A major obstacle in reinforcement learning-based sentence generation is the large action space whose size is equal to the vocabulary size of the target-side language. To improve the efficiency of reinforcement learning, we present a novel approach for reducing the action space based on dynamic vocabulary prediction. Our method first predicts a fixed-size small vocabulary for each input to generate its target sentence. The input-specific vocabularies are then used at supervised and reinforcement learning steps, and also at test time. In our experiments on six machine translation and two image captioning datasets, our method achieves faster reinforcement learning ($\sim$2.7x faster) with less GPU memory ($\sim$2.3x less) than the full-vocabulary counterpart. The reinforcement learning with our method consistently leads to significant improvement of BLEU scores, and the scores are equal to or better than those of baselines using the full vocabularies, with faster decoding time ($\sim$3x faster) on CPUs.


title: Schur complements of selfadjoint Krein space operators

abstract: Given a bounded selfadjoint operator W on a Krein space H and a closed subspace S of H, the Schur complement of W to S is defined under the hypothesis of weak complementability. A variational characterization of the Schur complement is given and the set of selfadjoint operators W admitting a Schur complement with these variational properties is shown to coincide with the set of S-weakly complementable selfadjoint operators.


title: TVQA: Localized, Compositional Video Question Answering

abstract: Recent years have witnessed an increasing interest in image-based question-answering (QA) tasks. However, due to data limitations, there has been much less work on video-based QA. In this paper, we present TVQA, a large-scale video QA dataset based on 6 popular TV shows. TVQA consists of 152,545 QA pairs from 21,793 clips, spanning over 460 hours of video. Questions are designed to be compositional in nature, requiring systems to jointly localize relevant moments within a clip, comprehend subtitle-based dialogue, and recognize relevant visual concepts. We provide analyses of this new dataset as well as several baselines and a multi-stream end-to-end trainable neural network framework for the TVQA task. The dataset is publicly available at http://tvqa.cs.unc.edu.


title: HASP: A High-Performance Adaptive Mobile Security Enhancement Against Malicious Speech Recognition

abstract: Nowadays, machine learning based Automatic Speech Recognition (ASR) technique has widely spread in smartphones, home devices, and public facilities. As convenient as this technology can be, a considerable security issue also raises -- the users' speech content might be exposed to malicious ASR monitoring and cause severe privacy leakage. In this work, we propose HASP -- a high-performance security enhancement approach to solve this security issue on mobile devices. Leveraging ASR systems' vulnerability to the adversarial examples, HASP is designed to cast human imperceptible adversarial noises to real-time speech and effectively perturb malicious ASR monitoring by increasing the Word Error Rate (WER). To enhance the practical performance on mobile devices, HASP is also optimized for effective adaptation to the human speech characteristics, environmental noises, and mobile computation scenarios. The experiments show that HASP can achieve optimal real-time security enhancement: it can lead an average WER of 84.55% for perturbing the malicious ASR monitoring, and the data processing speed is 15x to 40x faster compared to the state-of-the-art methods. Moreover, HASP can effectively perturb various ASR systems, demonstrating a strong transferability.

categories: cs.CR cs.LG cs.SD eess.AS eess.SP stat.ML


title: Biplanar Foldings

abstract: We are introducing a general framework for the construction of polyhedra and simplicial comlexes that are {\em bifoldable}, i.e. foldable into two two different planes. This vastly generalizes Origami folds known as the Miura pattern, the Eggbox pattern. After describing the framework and its basic features, we give several new examples of infinite triply periodic and fractal bifoldable polyhedra.


title: Data Lakes, Clouds and Commons: A Review of Platforms for Analyzing and Sharing Genomic Data

abstract: Data commons collate data with cloud computing infrastructure and commonly used software services, tools and applications to create biomedical resources for the large-scale management, analysis, harmonization, and sharing of biomedical data. Over the past few years, data commons have been used to analyze, harmonize and share large scale genomics datasets. Data ecosystems can be built by interoperating multiple data commons. It can be quite labor intensive to curate, import and analyze the data in a data commons. Data lakes provide an alternative to data commons and simply provide access to data, with the data curation and analysis deferred until later and delegated to those that access the data. We review software platforms for managing, analyzing and sharing genomic data, with an emphasis on data commons, but also covering data ecosystems and data lakes.


title: The nature of the variable millimetre-selected AGN in the Brightest Cluster Galaxy of Abell 851

abstract: We present the detection of a bright 3mm continuum source in the Brightest Cluster Galaxy (BCG) in Abell 851 (z=0.411) with the NOrthern Extended Millimeter Array (NOEMA). When this detection is compared to other multi-frequency observations across 21cm-100um, including new Arcminute Microkelvin Imager 15GHz observations, we find evidence for a relatively flat, variable core source associated with the BCG. The radio power and amplitude of variability observed in this galaxy is consistent with the cores in lower redshift BCGs in X-ray-selected clusters, and the flat mm-cm spectrum is suggestive of the BCG being a low luminosity AGN archetype. The discovery of this system could provide a basis for a long-term study of the role of low luminosity radio mode 'regulatory' feedback in massive clusters.


title: Pack and Detect: Fast Object Detection in Videos Using Region-of-Interest Packing

abstract: Object detection in videos is an important task in computer vision for various applications such as object tracking, video summarization and video search. Although great progress has been made in improving the accuracy of object detection in recent years due to the rise of deep neural networks, the state-of-the-art algorithms are highly computationally intensive. In order to address this challenge, we make two important observations in the context of videos: (i) Objects often occupy only a small fraction of the area in each video frame, and (ii) There is a high likelihood of strong temporal correlation between consecutive frames. Based on these observations, we propose Pack and Detect (PaD), an approach to reduce the computational requirements of object detection in videos. In PaD, only selected video frames called anchor frames are processed at full size. In the frames that lie between anchor frames (inter-anchor frames), regions of interest (ROIs) are identified based on the detections in the previous frame. We propose an algorithm to pack the ROIs of each inter-anchor frame together into a reduced-size frame. The computational requirements of the detector are reduced due to the lower size of the input. In order to maintain the accuracy of object detection, the proposed algorithm expands the ROIs greedily to provide additional background around each object to the detector. PaD can use any underlying neural network architecture to process the full-size and reduced-size frames. Experiments using the ImageNet video object detection dataset indicate that PaD can potentially reduce the number of FLOPS required for a frame by $4\times$. This leads to an overall increase in throughput of $1.25\times$ on a 2.1 GHz Intel Xeon server with a NVIDIA Titan X GPU at the cost of $1.1\%$ drop in accuracy.


title: Design and Realization of Intersection Traffic Control Simulation System Using Connected Vehicle Technology

abstract: The connected vehicle technology is a remarkable trend in the field of the intelligent transportation system. Since the actual deployment of the connected vehicle system is still lacking hitherto, simulation is widely adopted as the major method of verification in related researches. Although traditional commercial traffic simulation systems perform well in macroscopic simulations, they seem redundant in the circumstance of small-size traffic. In addition, those systems are unable to simulate the communication between vehicles and road-side infrastructures. This study designs a platform for the simulation of intersection traffic control using connected vehicle technology. By providing an abundance of customizable parameters, the platform can simulate at various levels of speed to meet different requirements, serving as a basis for testing further research.


title: HyperML: A Boosting Metric Learning Approach in Hyperbolic Space for Recommender Systems

abstract: This paper investigates the notion of learning user and item representations in non-Euclidean space. Specifically, we study the connection between metric learning in hyperbolic space and collaborative filtering by exploring Mobius gyrovector spaces where the formalism of the spaces could be utilized to generalize the most common Euclidean vector operations. Overall, this work aims to bridge the gap between Euclidean and hyperbolic geometry in recommender systems through metric learning approach. We propose HyperML (Hyperbolic Metric Learning), a conceptually simple but highly effective model for boosting the performance. Via a series of extensive experiments, we show that our proposed HyperML not only outperforms their Euclidean counterparts, but also achieves state-of-the-art performance on multiple benchmark datasets, demonstrating the effectiveness of personalized recommendation in hyperbolic geometry.


title: The Cluster-EAGLE project: a comparison of dynamical mass estimators using simulated clusters

abstract: Forthcoming large-scale spectroscopic surveys will soon provide data on thousands of galaxy clusters. It is important that the systematics of the various mass estimation techniques are well understood and calibrated. We compare three different dynamical mass estimators using the C-EAGLE galaxy clusters, a set of high resolution simulations with resolved galaxies a median total mass, $M_{200c} = 10^{14.7} \, \mathrm{M_\odot}$. We quantify the bias and scatter of the Jeans, virial, and caustic mass estimators using all galaxies with a stellar mass $M_*> 10^9 \, \mathrm{M_\odot}$, both in the ideal 3D case and in the more realistic projected case. On average we find our mass estimates are unbiased, though relative to the true mass within $r_{200c}$ the scatter is large with a range of $0.09$ - $0.15$ dex. We see a slight increase in the scatter when projecting the clusters. Selecting galaxies using the same criteria, we find no significant difference in the mass bias or scatter when comparing results from hydrodynamical and dark matter only simulations. However, selecting galaxies by stellar mass reduces the bias compared to selecting by total mass. Comparing X-ray derived hydrostatic and dynamical masses, the former are ${\sim} 30$ per cent lower. We find a slight dependence between substructure, measured using two different metrics, and mass bias. In conclusion, we find that dynamical mass estimators, when averaged together, are unbiased with a scatter of $0.11 \pm 0.02$ dex when including interloper galaxies and with no prior knowledge of $r_{200c}$.


title: H2CO ortho-to-para ratio in the protoplanetary disk HD 163296

abstract: Ortho-to-para (o/p) ratios of species like water, ammonia and formaldehyde (H2CO) are believed to encode information about the formation history of the molecule. Measurements of o/p ratios in protoplanetary disks could thus be used to constrain their physical and chemical histories. We present the first measurement of the H2CO o/p ratio in a protoplanetary disk, using three ortho and two para lines observed with the Sub-millimeter Array (SMA) combined with one highly resolved measurement of a single H2CO line with ALMA toward the disk around Herbig Ae star HD 163296. We find a disk-averaged H2CO o/p ratio of 1.8-2.8 (depending on the assumed disk structure), corresponding to a spin temperature of 11-22 K. We also derive a rotational temperature of 24 K from the flux ratio of the three ortho lines. The observed spatial distribution, as seen by ALMA, as well as the rotational temperature and the o/p ratio, at the large scales the SMA is most sensitive to, are consistent with a low-temperature formation pathway, most likely grain surface chemistry, of H2CO in this disk.


title: A Limitation of V-Matrix based Methods

abstract: To estimate the conditional probability functions based on the direct problem setting, V-matrix based method was proposed. We construct V-matrix based constrained quadratic programming problems for which the inequality constraints are inconsistent. In particular, we would like to present that the constrained quadratic optimization problem for conditional probability estimation using V-matrix method may not have a consistent solution always.


title: ACCDB: A Collection of Chemistry DataBases for Broad Computational Purposes

abstract: The importance of databases of reliable and accurate data in chemistry has substantially increased in the past two decades. Their main usage is to parametrize electronic structure theory methods, and to assess their capabilities and accuracy for a broad set of chemical problems. The collection we present here -ACCDB- includes data from 16 different research groups, for a total of 44,931 unique reference data points, all at a level of theory significantly higher than DFT, and covering most of the periodic table. It is composed of five databases taken from literature (GMTKN, MGCDB84, Minnesota2015B, DP284, and W4-17), two newly developed reaction energy databases (W4-17-RE, and MN-RE), and a new collection of databases containing transition metals. A set of expandable software tools for its manipulation is also presented here for the first time, as well as a case study where ACCDB is used for benchmarking commercial CPUs for chemistry calculations. To be published in J. Comput. Chem.


title: Statistical and mathematical modeling of spatiotemporal dynamics of stem cells

abstract: Statistical and mathematical modeling are crucial to describe, interpret, compare and predict the behavior of complex biological systems including the organization of hematopoietic stem and progenitor cells in the bone marrow environment. The current prominence of high-resolution and live-cell imaging data provides an unprecedented opportunity to study the spatiotemporal dynamics of these cells within their stem cell niche and learn more about aberrant, but also unperturbed, normal hematopoiesis. However, this requires careful quantitative statistical analysis of the spatial and temporal behavior of cells and the interaction with their microenvironment. Moreover, such quantification is a prerequisite for the construction of hypothesis-driven mathematical models that can provide mechanistic explanations by generating spatiotemporal dynamics that can be directly compared to experimental observations. Here, we provide a brief overview of statistical methods in analyzing spatial distribution of cells, cell motility, cell shapes and cellular genealogies. We also describe cell- based modeling formalisms that allow researchers to simulate emergent behavior in a multicellular system based on a set of hypothesized mechanisms. Together, these methods provide a quantitative workflow for the analytic and synthetic study of the spatiotemporal behavior of hematopoietic stem and progenitor cells.


title: Search for an intrinsic metallicity spread in old globular clusters of the Large Magellanic Cloud

abstract: We report for the first time on the magnitude of the intrinsic [Fe/H] spread among ten old globular clusters (GCs) of the Large Magellanic Cloud (LMC). Such spreads are merely observed in approximately five per cent of the Milky Way GCs and recently gained more attention in theoretical models of GC evolution. We derived metallicities with a typical precision of 0.05 dex < sigma[Fe/H] < 0.20 dex for an average of 14 red giant branch stars per GC from Str\"omgren photometry. The respective, metallicity-sensitive indices have been calibrated to precise and accurate high-dispersion spectroscopy. For all clusters we found null [Fe/H] spreads with a typical uncertainty of 0.04 dex, with the possible exception of NGC1786 that shows an intrinsic dispersion of 0.07+-0.04 dex. The mean, observed standard deviation of the derived metallicities for nearly 40 per cent of our GC sample amounted to smaller than 0.05 dex. At present, we cannot exclude that the remaining GCs also have intrinsic Fe-abundance variations in excess of 0.05 dex, but in order to significantly detect those, the measurement errors on individual [Fe/H]-values would need to be lowered to the 0.03--0.07 dex level. These findings suggest, along with those from ages and light-element abundances, that the LMC GCs studied here are alike to the majority of Galactic GCs.


title: Real-time Optimal Resource Allocation for Embedded UAV Communication Systems

abstract: We consider device-to-device (D2D) wireless information and power transfer systems using an unmanned aerial vehicle (UAV) as a relay-assisted node. As the energy capacity and flight time of UAVs is limited, a significant issue in deploying UAV is to manage energy consumption in real-time application, which is proportional to the UAV transmit power. To tackle this important issue, we develop a real-time resource allocation algorithm for maximizing the energy efficiency by jointly optimizing the energy-harvesting time and power control for the considered (D2D) communication embedded with UAV. We demonstrate the effectiveness of the proposed algorithms as running time for solving them can be conducted in milliseconds.


title: Optimal On-Off Control for a Class of Discrete Event Systems with Real-Time Constraints

abstract: This paper studies an optimal ON-OFF control problem for a class of discrete event systems with real-time constraints. Our goal is to minimize the overall costs, including the operating cost and the wake-up cost, while still guaranteeing the deadline of each individual task. In particular, we consider the homogeneous case in which it takes the same amount of time to serve each task and each task needs to be served by d seconds upon arrival. The problem involves two subproblems: (i) finding the best time to wake up the system and (ii) finding the best time to let the system go to sleep. We study the two subproblems in both off-line and on-line settings. In the off-line case that all task information is known a priori, we combine sample path analysis and dynamic programming to come up with the optimal solution. In the on-line scenario where future task information is completely unknown, we show that the optimal time to wake up the system can be obtained without relying on future task arrivals. We also perform competitive analysis for on-line control and derive the competitive ratios for both deterministic and random controllers.


title: Coverage-Based Designs Improve Sample Mining and Hyper-Parameter Optimization

abstract: Sampling one or more effective solutions from large search spaces is a recurring idea in machine learning, and sequential optimization has become a popular solution. Typical examples include data summarization, sample mining for predictive modeling and hyper-parameter optimization. Existing solutions attempt to adaptively trade-off between global exploration and local exploitation, wherein the initial exploratory sample is critical to their success. While discrepancy-based samples have become the de facto approach for exploration, results from computer graphics suggest that coverage-based designs, e.g. Poisson disk sampling, can be a superior alternative. In order to successfully adopt coverage-based sample designs to ML applications, which were originally developed for 2-d image analysis, we propose fundamental advances by constructing a parameterized family of designs with provably improved coverage characteristics, and by developing algorithms for effective sample synthesis. Using experiments in sample mining and hyper-parameter optimization for supervised learning, we show that our approach consistently outperforms existing exploratory sampling methods in both blind exploration, and sequential search with Bayesian optimization.


title: The Cygnus Loop's Distance, Properties, and Environment Driven Morphology

abstract: The Cygnus Loop is among the brightest and best studied evolved Galactic supernova remnants. However, its distance has remained uncertain thus undermining quantitative understanding about many of its fundamental properties. Here we present moderate-dispersion spectra of stars with projected locations toward the remnant. Spectra of three stars revealed Na I 5890,5896 A and Ca II 3934 A absorption features associated with the remnant's expanding shell, with velocities ranging from -160 to +240 km/s. Combining Gaia DR2 parallax measurements for these stars with other recent observations, we find the distance to the Cygnus Loop's centre is 735 +/-25 pc, only a bit less than the 770 pc value proposed by Minkowski some 60 years ago. Using this new distance, we discuss the remnant's physical properties including size, SN explosion energy, and shock velocities. We also present multi-wavelength emission maps which reveal that, instead of being located in a progenitor wind-driven cavity as has long been assumed, the Cygnus Loop lies in an extended, low density region. Rather than wind-driven cavity walls, these images reveal in unprecedented clarity the sizes and locations of local interstellar clouds with which the remnant is interacting, giving rise to its large-scale morphology.


title: Cartier modules and cyclotomic spectra

abstract: We construct and study a t-structure on p-typical cyclotomic spectra and explain how to recover crystalline cohomology of smooth schemes over perfect fields using this t-structure. Our main tool is a new approach to p-typical cyclotomic spectra via objects we call p-typical topological Cartier modules. Using these, we prove that the heart of the cyclotomic t-structure is the full subcategory of derived V-complete objects in the abelian category of p-typical Cartier modules.


title: Bridging machine learning and cryptography in defence against adversarial attacks

abstract: In the last decade, deep learning algorithms have become very popular thanks to the achieved performance in many machine learning and computer vision tasks. However, most of the deep learning architectures are vulnerable to so called adversarial examples. This questions the security of deep neural networks (DNN) for many security- and trust-sensitive domains. The majority of the proposed existing adversarial attacks are based on the differentiability of the DNN cost function.Defence strategies are mostly based on machine learning and signal processing principles that either try to detect-reject or filter out the adversarial perturbations and completely neglect the classical cryptographic component in the defence. In this work, we propose a new defence mechanism based on the second Kerckhoffs's cryptographic principle which states that the defence and classification algorithm are supposed to be known, but not the key. To be compliant with the assumption that the attacker does not have access to the secret key, we will primarily focus on a gray-box scenario and do not address a white-box one. More particularly, we assume that the attacker does not have direct access to the secret block, but (a) he completely knows the system architecture, (b) he has access to the data used for training and testing and (c) he can observe the output of the classifier for each given input. We show empirically that our system is efficient against most famous state-of-the-art attacks in black-box and gray-box scenarios.


title: 3D classical and quantum stable structures of dissipative systems

abstract: We study the properties of classical and quantum stable structures in a 3D parameter space corresponding to the dissipative kicked top. This is a model system in quantum and classical chaos that gives a starting point for many body examples. We are able to identify the influence of these structures in the spectra and eigenstates of the corresponding (super)operators. This provides with a complementary view with respect to the typical 2D parameter space systems found in the literature. Many properties of the eigenstates, like its localization behaviour can be generalized to this higher dimensional parameter space and spherical phase space topology. Moreover we find a 3D phenomenon -- generalizable to more dimensions -- that we call the coalescence-separation of (q)ISSs, whose main consequence is a marked enhancement of quantum localization. This could be of relevance for systems which have attracted a lot of attention very recently.


title: Approximating Bimatrix Nash Equilibrium Via Trilinear Minimax

abstract: The Bimatrix Nash Equilibrium (NE) for $m \times n$ real matrices $R$ and $C$, denoted as the {\it Row} and {\it Column} players, is characterized as follows: Let $\Delta =S_m \times S_n$, where $S_k$ denotes the unit simplex in $\mathbb{R}^k$. For a given point $p=(x,y) \in \Delta$, define $R[p]=x^TRy$ and $C[p]=x^TCy$. Consequently, there exists a subset $\Delta_* \subset \Delta$ such that for any $p_*=(x_*,y_*) \in \Delta_*$, $\max_{p \in \Delta, y=y_*}R[p]=R[p_*]$ and $\max_{p \in \Delta, x=x_* } C[p]=C[p_*]$. The computational complexity of bimatrix NE falls within the class of {\it PPAD-complete}. Although the von Neumann Minimax Theorem is a special case of bimatrix NE, we introduce a novel extension termed {\it Trilinear Minimax Relaxation} (TMR) with the following implications: Let $\lambda^*=\min_{\alpha \in S_{2}} \max_{p \in \Delta} (\alpha_1 R[p]+ \alpha_2C[p])$ and $\lambda_*=\max_{p \in \Delta} \min_{\alpha \in S_{2}} (\alpha_1 R[p]+ \alpha_2C[p])$. $\lambda^* \geq \lambda_*$. $\lambda^*$ is computable as a linear programming in $O(mn)$ time, ensuring $\max_{p_* \in \Delta_*}\min \{R[p_*], C[p_*]\} \leq \lambda^*$, meaning that in any Nash Equilibrium it is not possible to have both players' payoffs to exceed $\lambda^*$. $\lambda^*=\lambda_*$ if and only if there exists $p^* \in \Delta$ such that $\lambda^*= \min\{R[p^*], C[p^*]\}$. Such a $p^*$ serves as an approximate Nash Equilibrium. We analyze the cases where such $p^*$ exists and is computable. Even when $\lambda^* > \lambda_*$, we derive approximate Nash Equilibria. In summary, the aforementioned properties of TMR and its efficient computational aspects underscore its significance and relevance for Nash Equilibrium, irrespective of the computational complexity associated with bimatrix Nash Equilibrium. Finally, we extend TMR to scenarios involving three or more players.


title: On one generalization of skew tent maps

abstract: We generalize in this work the properties of the conjugacy of skew tent maps. It is known that the conjugacy $h$ from a skew tent map $g_1$ to $g_2$ is differentiable at a point $x^*$ if and only if there exists left and right limits $\lim\limits_{n\rightarrow \infty}h_n'(x^*-)$ and $\lim\limits_{n\rightarrow \infty}h_n'(x^*+)$, where $h_n$ is a piecewise linear function, which coincides with $h$ at $g_1^{-n}(0)$, and all whose kinks belong to $g_1^{-n}(0)$. The attempts to generalize this result to some reacher class of unimodal maps is natural. For this reason we introduce the class of piecewise linear maps, all whose kinks are in the complete pre-image of $0$ and study the relation of the differential properties of their conjugacy with ones of the mentioned approximation~$(h_n)_{n\geq 1}$.


title: Q-zero range has random walking shocks

abstract: ... but no other surprises in the zero range world. We check all nearest neighbour 1-dimensional asymmetric zero range processes for random walking product shock measures as demonstrated already for a few cases in the literature. We find the totally asymmetric version of the celebrated q-zero range process as the only new example besides an already known model of doubly infinite occupation numbers and exponentially increasing jump rates. We also examine the interaction of shocks, which appears somewhat more involved for q-zero range than in the already known cases.


title: Extending Mandelbox Fractals with Shape Inversions

abstract: The Mandelbox is a recently discovered class of escape-time fractals which use a conditional combination of reflection, spherical inversion, scaling, and translation to transform points under iteration. In this paper we introduce a new extension to Mandelbox fractals which replaces spherical inversion with a more generalized shape inversion. We then explore how this technique can be used to generate new fractals in 2D, 3D, and 4D.


title: Three-Stage Speaker Verification Architecture in Emotional Talking Environments

abstract: Speaker verification performance in neutral talking environment is usually high, while it is sharply decreased in emotional talking environments. This performance degradation in emotional environments is due to the problem of mismatch between training in neutral environment while testing in emotional environments. In this work, a three-stage speaker verification architecture has been proposed to enhance speaker verification performance in emotional environments. This architecture is comprised of three cascaded stages: gender identification stage followed by an emotion identification stage followed by a speaker verification stage. The proposed framework has been evaluated on two distinct and independent emotional speech datasets: in-house dataset and Emotional Prosody Speech and Transcripts dataset. Our results show that speaker verification based on both gender information and emotion information is superior to each of speaker verification based on gender information only, emotion information only, and neither gender information nor emotion information. The attained average speaker verification performance based on the proposed framework is very alike to that attained in subjective assessment by human listeners.


title: Unexpected sawtooth artifact in beat-to-beat pulse transit time measured from patient monitor data

abstract: Object: It is increasingly popular to collect as much data as possible in the hospital setting from clinical monitors for research purposes. However, in this setup the data calibration issue is often not discussed and, rather, implicitly assumed, while the clinical monitors might not be designed for the data analysis purpose. We hypothesize that this calibration issue for a secondary analysis may become an important source of artifacts in patient monitor data. We test an off-the-shelf integrated photoplethysmography (PPG) and electrocardiogram (ECG) monitoring device for its ability to yield a reliable pulse transit time (PTT) signal. Approach: This is a retrospective clinical study using two databases: one containing 35 subjects who underwent laparoscopic cholecystectomy, another containing 22 subjects who underwent spontaneous breathing test in the intensive care unit. All data sets include recordings of PPG and ECG using a commonly deployed patient monitor. We calculated the PTT signal offline. Main Results: We report a novel constant oscillatory pattern in the PTT signal and identify this pattern as a sawtooth artifact. We apply an approach based on the de-shape method to visualize, quantify and validate this sawtooth artifact. Significance: The PPG and ECG signals not designed for the PTT evaluation may contain unwanted artifacts. The PTT signal should be calibrated before analysis to avoid erroneous interpretation of its physiological meaning.

categories: q-bio.QM cs.LG eess.SP physics.data-an stat.AP


title: Tunable Quantum Hall Edge Conduction in Bilayer Graphene through Spin-Orbit Interaction

abstract: Bilayer graphene, in the presence of a one-sided spin-orbit interaction (SOI) induced by a suitably chosen substrate, is predicted to exhibit unconventional Quantum Hall states. The new states arise due to strong SOI-induced splittings of the eight zeroth Landau levels, which are strongly layer-polarized, residing fully or partially on one of the two graphene layers. In particular, an Ising SOI in the meV scale is sufficient to invert the Landau level order between the $n=0$ and $n=1$ orbital levels under moderately weak magnetic fields $B \lesssim 10$\~T. Furthermore, when the Ising field opposes the $B$ field, the order of the spin-polarized levels can also be inverted. We show that, under these conditions, three different compensated electron-hole phases, with equal concentrations of electrons and holes, can occur at $\nu = 0$ filling. The three phases have distinct edge conductivity values. One of the phases is especially interesting, since its edge conduction can be turned on and off by switching the sign of the interlayer bias.


title: Langevin equations in the small-mass limit: Higher-order approximations

abstract: We study the small-mass (overdamped) limit of Langevin equations for a particle in a potential and/or magnetic field with matrix-valued and state-dependent drift and diffusion. We utilize a bootstrapping argument to derive a hierarchy of approximate equations for the position degrees of freedom that are able to achieve accuracy of order $m^{\ell/2}$ over compact time intervals for any $\ell\in\mathbb{Z}^+$. This generalizes prior derivations of the homogenized equation for the position degrees of freedom in the $m\to 0$ limit, which result in order $m^{1/2}$ approximations. Our results cover bounded forces, for which we prove convergence in $L^p$ norms, and unbounded forces, in which case we prove convergence in probability.


title: Logarithmic W-algebras and Argyres-Douglas theories at higher rank

abstract: Families of vertex algebras associated to nilpotent elements of simply-laced Lie algebras are constructed. These algebras are close cousins of logarithmic W-algebras of Feigin and Tipunin and they are also obtained as modifications of semiclassical limits of vertex algebras appearing in the context of $S$-duality for four-dimensional gauge theories. In the case of type $A$ and principal nilpotent element the character agrees precisely with the Schur-Index formula for corresponding Argyres-Douglas theories with irregular singularities. For other nilpotent elements they are identified with Schur-indices of type IV Argyres-Douglas theories. Further, there is a conformal embedding pattern of these vertex operator algebras that nicely matches the RG-flow of Argyres-Douglas theories as discussed by Buican and Nishinaka.


title: Neural Comic Style Transfer: Case Study

abstract: The work by Gatys et al. [1] recently showed a neural style algorithm that can produce an image in the style of another image. Some further works introduced various improvements regarding generalization, quality and efficiency, but each of them was mostly focused on styles such as paintings, abstract images or photo-realistic style. In this paper, we present a comparison of how state-of-the-art style transfer methods cope with transferring various comic styles on different images. We select different combinations of Adaptive Instance Normalization [11] and Universal Style Transfer [16] models and confront them to find their advantages and disadvantages in terms of qualitative and quantitative analysis. Finally, we present the results of a survey conducted on over 100 people that aims at validating the evaluation results in a real-life application of comic style transfer.


title: Nonlinear Mixed Modes in Red Giants

abstract: Turbulent motions in the convective envelope of red giants excite a rich spectrum of solar-like oscillation modes. Observations by CoRoT and Kepler have shown that the mode amplitudes increase dramatically as the stars ascend the red giant branch, i.e., as the frequency of maximum power, $\nu_\mathrm{max}$, decreases. Most studies nonetheless assume that the modes are well described by the linearized fluid equations. We investigate to what extent the linear approximation is justified as a function of stellar mass $M$ and $\nu_\mathrm{max}$, focusing on dipole mixed modes with frequency near $\nu_\mathrm{max}$. A useful measure of a mode's nonlinearity is the product of its radial wavenumber and its radial displacement, $k_r \xi_r$ (i.e., its shear). We show that $k_r \xi_r \propto \nu_\mathrm{max}^{-9/2}$, implying that the nonlinearity of mixed modes increases significantly as a star evolves. The modes are weakly nonlinear ($k_r \xi_r > 10^{-3}$) for $\nu_\mathrm{max} \lesssim 150 \, \mu\mathrm{Hz}$ and strongly nonlinear ($k_r \xi_r > 1$) for $\nu_\mathrm{max} \lesssim 30 \, \mu\mathrm{Hz}$, with only a mild dependence on $M$ over the range we consider ($1.0 - 2.0 M_\odot$). A weakly nonlinear mixed mode can excite secondary waves in the stellar core through the parametric instability, resulting in enhanced, but partial, damping of the mode. By contrast, a strongly nonlinear mode breaks as it propagates through the core and is fully damped there. Evaluating the impact of nonlinear effects on observables such as mode amplitudes and linewidths requires large mode network simulations. We plan to carry out such calculations in the future and investigate whether nonlinear damping can explain why some red giants exhibit dipole modes with unusually small amplitudes, known as depressed modes.


title: Attention-based Audio-Visual Fusion for Robust Automatic Speech Recognition

abstract: Automatic speech recognition can potentially benefit from the lip motion patterns, complementing acoustic speech to improve the overall recognition performance, particularly in noise. In this paper we propose an audio-visual fusion strategy that goes beyond simple feature concatenation and learns to automatically align the two modalities, leading to enhanced representations which increase the recognition accuracy in both clean and noisy conditions. We test our strategy on the TCD-TIMIT and LRS2 datasets, designed for large vocabulary continuous speech recognition, applying three types of noise at different power ratios. We also exploit state of the art Sequence-to-Sequence architectures, showing that our method can be easily integrated. Results show relative improvements from 7% up to 30% on TCD-TIMIT over the acoustic modality alone, depending on the acoustic noise level. We anticipate that the fusion strategy can easily generalise to many other multimodal tasks which involve correlated modalities. Code available online on GitHub: https://github.com/georgesterpu/Sigmedia-AVSR


title: On the forced Euler and Navier-Stokes equations: Linear damping and modified scattering

abstract: We study the asymptotic behavior of the forced linear Euler and nonlinear Navier-Stokes equations close to Couette flow in a periodic channel. As our main result we show that for smooth time-periodic forcing linear inviscid damping persists, i.e. the velocity field (weakly) asymptotically converges. However, stability and scattering to the transport problem fail in $H^{s}, s>-1$. We further show that this behavior is consistent with the nonlinear Euler equations and that a similar result also holds for the nonlinear Navier-Stokes equations. Hence, these results provide an indication that nonlinear inviscid damping may still hold in Sobolev regularity in the above sense despite the Gevrey regularity instability results of [Deng-Masmoudi 2018].


title: Measuring the De Sitter precession with a new Earth's satellite to the $\mathbf{\simeq 10^{-5}}$ level: a proposal

abstract: The inclination $I$ of an Earth's satellite in polar orbit undergoes a secular De Sitter precession of $-7.6$ milliarcseconds per year for a suitable choice of the initial value of its non-circulating node $\Omega$. The competing long-periodic harmonic rates of change of $I$ due to the even and odd zonal harmonics of the geopotential vanish for either a circular or polar orbit, while no secular rates occur at all. This may open up, in principle, the possibility of measuring the geodesic precession in the weak-field limit with an accurately tracked satellite by improving the current bound of $9\times 10^{-4}$ from Lunar Laser Ranging, which, on the other hand, may be even rather optimistic, by one order of magnitude, or, perhaps, even better. The most insidious competing effects are due to the solid and ocean components of the $K_1$ tide since their perturbations have nominal huge amplitudes and the same temporal pattern of the De Sitter signature. They vanish for polar orbits. Departures of $\simeq 10^{-5}-10^{-3}~\textrm{deg}$ from the ideal polar geometry allow to keep the $K_1$ tidal perturbations to a sufficiently small level. Most of the other gravitational and non-gravitational perturbations vanish for the proposed orbital configuration, while the non-vanishing ones either have different temporal signatures with respect to the De Sitter effect or can be modeled with sufficient accuracy. In order to meet the proposed goal, the measurement accuracy of $I$ should be better than $\simeq 35~\textrm{microarcseconds}=0.034~\textrm{milliarcseconds}$ over, say, 5 yr.


title: Quantum algorithms and lower bounds for convex optimization

abstract: While recent work suggests that quantum computers can speed up the solution of semidefinite programs, little is known about the quantum complexity of more general convex optimization. We present a quantum algorithm that can optimize a convex function over an $n$-dimensional convex body using $\tilde{O}(n)$ queries to oracles that evaluate the objective function and determine membership in the convex body. This represents a quadratic improvement over the best-known classical algorithm. We also study limitations on the power of quantum computers for general convex optimization, showing that it requires $\tilde{\Omega}(\sqrt n)$ evaluation queries and $\Omega(\sqrt{n})$ membership queries.


title: Decomposition of the Feynman kernel for a particle in a box

abstract: We study the decomposition of the Feynman kernel for a particle in a box with $1/\sin^{2}\theta$ potential to find that the wellknown phase factor $-1$, which is correct for the case of the free particle, for reflection at boundaries should be generalized depending on the parameter of the potential.


title: Deep Joint Source-Channel Coding for Wireless Image Transmission

abstract: We propose a joint source and channel coding (JSCC) technique for wireless image transmission that does not rely on explicit codes for either compression or error correction; instead, it directly maps the image pixel values to the complex-valued channel input symbols. We parameterize the encoder and decoder functions by two convolutional neural networks (CNNs), which are trained jointly, and can be considered as an autoencoder with a non-trainable layer in the middle that represents the noisy communication channel. Our results show that the proposed deep JSCC scheme outperforms digital transmission concatenating JPEG or JPEG2000 compression with a capacity achieving channel code at low signal-to-noise ratio (SNR) and channel bandwidth values in the presence of additive white Gaussian noise (AWGN). More strikingly, deep JSCC does not suffer from the ``cliff effect'', and it provides a graceful performance degradation as the channel SNR varies with respect to the SNR value assumed during training. In the case of a slow Rayleigh fading channel, deep JSCC learns noise resilient coded representations and significantly outperforms separation-based digital communication at all SNR and channel bandwidth values.


title: Every Elementary Higher Topos has a Natural Number Object

abstract: We prove that every elementary $(\infty,1)$-topos has a natural number object. We achieve this by defining the loop space of the circle and showing that we can construct a natural number object out of it. Part of the proof involves showing that various definitions of natural number objects (Lawvere, Freyd and Peano) agree with each other in an elementary $(\infty,1)$-topos. As part of this effort we also study the internal object of contractibility in $(\infty,1)$-categories, which is of independent interest. Finally, we discuss various applications of natural number objects. In particular, we use it to define internal sequential colimits in an elementary $(\infty,1)$-topos.


title: Local and non-local multiplicative Poisson vertex algebras and differential-difference equations

abstract: We develop the notions of multiplicative Lie conformal and Poisson vertex algebras, local and non-local, and their connections to the theory of integrable differential-difference Hamiltonian equations. We establish relations of these notions to $q$-deformed $W$-algebras and lattice Poisson algebras. We introduce the notion of Adler type pseudodifference operators and apply them to integrability of differential-difference Hamiltonian equations.


title: Yoneda Lemma for Elementary Higher Toposes

abstract: We prove the Yoneda lemma inside an elementary higher topos, generalizing the Yonda lemma for spaces.


title: Combining Many-objective Radiomics and 3-dimensional Convolutional Neural Network through Evidential Reasoning to Predict Lymph Node Metastasis in Head and Neck Cancer

abstract: Lymph node metastasis (LNM) is a significant prognostic factor in patients with head and neck cancer, and the ability to predict it accurately is essential for treatment optimization. PET and CT imaging are routinely used for LNM identification. However, uncertainties of LNM always exist especially for small size or reactive nodes. Radiomics and deep learning are the two preferred imaging-based strategies for node malignancy prediction. Radiomics models are built based on handcrafted features, and deep learning can learn the features automatically. We proposed a hybrid predictive model that combines many-objective radiomics (MO-radiomics) and 3-dimensional convolutional neural network (3D-CNN) through evidential reasoning (ER) approach. To build a more reliable model, we proposed a new many-objective radiomics model. Meanwhile, we designed a 3D-CNN that fully utilizes spatial contextual information. Finally, the outputs were fused through the ER approach. To study the predictability of the two modalities, three models were built for PET, CT, and PET&CT. The results showed that the model performed best when the two modalities were combined. Moreover, we showed that the quantitative results obtained from the hybrid model were better than those obtained from MO-radiomics and 3D-CNN.


title: Recovering a Single Community with Side Information

abstract: We study the effect of the quality and quantity of side information on the recovery of a hidden community of size $K=o(n)$ in a graph of size $n$. Side information for each node in the graph is modeled by a random vector with the following features: either the dimension of the vector is allowed to vary with $n$, while log-likelihood ratio (LLR) of each component with respect to the node label is fixed, or the LLR is allowed to vary and the vector dimension is fixed. These two models represent the variation in quality and quantity of side information. Under maximum likelihood detection, we calculate tight necessary and sufficient conditions for exact recovery of the labels. We demonstrate how side information needs to evolve with $n$ in terms of either its quantity, or quality, to improve the exact recovery threshold. A similar set of results are obtained for weak recovery. Under belief propagation, tight necessary and sufficient conditions for weak recovery are calculated when the LLRs are constant, and sufficient conditions when the LLRs vary with $n$. Moreover, we design and analyze a local voting procedure using side information that can achieve exact recovery when applied after belief propagation. The results for belief propagation are validated via simulations on finite synthetic data-sets, showing that the asymptotic results of this paper can also shed light on the performance at finite $n$.


title: Join-the-Shortest Queue Diffusion Limit in Halfin-Whitt Regime: Sensitivity on the Heavy-traffic Parameter

abstract: Consider a system of $N$ parallel single-server queues with unit-exponential service time distribution and a single dispatcher where tasks arrive as a Poisson process of rate $\lambda(N)$. When a task arrives, the dispatcher assigns it to one of the servers according to the Join-the-Shortest Queue (JSQ) policy. Eschenfeldt and Gamarnik (Math. Oper. Res., 43(3):867-886, 2018) identified a novel limiting diffusion process that arises as the weak-limit of the appropriately scaled occupancy measure of the system under the JSQ policy in the Halfin-Whitt regime, where $(N - \lambda(N)) / \sqrt{N} \to \beta > 0$ as $N \to \infty$. The analysis of this diffusion goes beyond the state of the art techniques, and even proving its ergodicity is non-trivial, and was left as an open question. Recently, exploiting a generator expansion framework via the Stein's method, Braverman (arXiv:1801.05121, 2018) established its exponential ergodicity, and adapting a regenerative approach, Banerjee and Mukherjee (Ann. Appl. Probab., 29(2):1262-1309, 2018) analyzed the tail properties of the stationary distribution and path fluctuations of the diffusion. However, the analysis of the bulk behavior of the stationary distribution, viz., the moments, remained intractable until this work. In this paper, we perform a thorough analysis of the bulk behavior of the stationary distribution of the diffusion process, and discover that it exhibits different qualitative behavior, depending on the value of the heavy-traffic parameter $\beta$. Moreover, we obtain precise asymptotic laws of the centered and scaled steady state distribution, as $\beta$ tends to 0 and $\infty$. Of particular interest, we also establish a certain intermittency phenomena in the $\beta\to \infty$ regime and a surprising distributional convergence result in the $\beta\to 0$ regime.


title: Predicting Smoking Events with a Time-Varying Semi-Parametric Hawkes Process Model

abstract: Health risks from cigarette smoking -- the leading cause of preventable death in the United States -- can be substantially reduced by quitting. Although most smokers are motivated to quit, the majority of quit attempts fail. A number of studies have explored the role of self-reported symptoms, physiologic measurements, and environmental context on smoking risk, but less work has focused on the temporal dynamics of smoking events, including daily patterns and related nicotine effects. In this work, we examine these dynamics and improve risk prediction by modeling smoking as a self-triggering process, in which previous smoking events modify current risk. Specifically, we fit smoking events self-reported by 42 smokers to a time-varying semi-parametric Hawkes process (TV-SPHP) developed for this purpose. Results show that the TV-SPHP achieves superior prediction performance compared to related and existing models, with the incorporation of time-varying predictors having greatest benefit over longer prediction windows. Moreover, the impact function illustrates previously unknown temporal dynamics of smoking, with possible connections to nicotine metabolism to be explored in future work through a randomized study design. By more effectively predicting smoking events and exploring a self-triggering component of smoking risk, this work supports development of novel or improved cessation interventions that aim to reduce death from smoking.


title: The Carnegie-Chicago Hubble Program. V. The Distances to NGC 1448 and NGC 1316 via the Tip of the Red Giant Branch

abstract: The Carnegie-Chicago Hubble Program (CCHP) is re-calibrating the extragalactic SN Ia distance scale using exclusively Population II stars. This effort focuses on the Tip of the Red Giant Branch (TRGB) method, whose systematics are entirely independent of the Population I Cepheid-based determinations that have long served as calibrators for the SN Ia distance scale. We present deep Hubble Space Telescope imaging of the low surface-density and low line-of-sight reddening halos of two galaxies, NGC 1448 and NGC 1316, each of which have been hosts to recent SN Ia events. Provisionally anchoring the TRGB zero-point to the geometric distance to the Large Magellanic Cloud derived from detached eclipsing binaries, we measure extinction-corrected distance moduli of 31.23 +/-0.04 (stat) +/- 0.06 (sys) mag for NGC 1448 and 31.37 +/- 0.04 (stat) and +/- 0.06 (sys) mag for NGC 1316, respectively, giving metric distances of 17.7 +/- 0.3 (stat) +/- 0.5 (sys) Mpc, and 18.8 +/- 0.3 (stat) +/- 0.5 (sys) Mpc. We find agreement between our result and the available Cepheid distance for NGC 1448; for NGC 1316, where there are relatively few published distances based on direct measurements, we find that our result is consistent with the published SN Ia distances whose absolute scales are set from other locally-determined methods such as Cepheids. For NGC 1448 and NGC 1316, our distances are some of the most precise (and systematically accurate) measurements with errors at 1.7 (2.8) % and 1.6 (2.7) % levels, respectively.


title: On the wellposedness of some McKean models with moderated or singular diffusion coefficient

abstract: We investigate the well-posedness problem related to two models of nonlinear McKean Stochastic Differential Equations with some local interaction in the diffusion term. First, we revisit the case of the McKean-Vlasov dynamics with moderate interaction, previously studied by Meleard and Jourdain in [16], under slightly weaker assumptions, by showing the existence and uniqueness of a weak solution using a Sobolev regularity framework instead of a Holder one. Second, we study the construction of a Lagrangian Stochastic model endowed with a conditional McKean diffusion term in the velocity dynamics and a nondegenerate diffusion term in the position dynamics.


title: A plethora of new R Coronae Borealis stars discovered from a dedicated spectroscopic follow-up survey

abstract: It is increasingly suspected that the rare R Coronae Borealis (RCB) stars - hydrogen-deficient and carbon-rich supergiant stars - are the products of mergers of CO/He white-dwarf binary systems in the intermediate mass regime ($0.6<M_{Tot}<1.2 M_{\odot}$). Only 77 RCB stars are currently known in our Galaxy while up to 1000 were expected. It is necessary to find more of these peculiar and diverse stars to understand their origin and evolutionary path. We are undertaking such a dedicated search. We plan to follow up spectroscopically 2356 targets of interest that were carefully selected using the all sky 2MASS and WISE surveys. We have observed nearly 500 of these targets using optical low-resolution spectrographs. These spectra were compared to synthetic spectra from a new grid of MARCs hydrogen-deficient atmospheric models. Classical RCB stars' photospheric temperatures range mostly from 4000 K to about 8500 K, and therefore their spectra look very different showing the presence of carbon molecules C$_2$ and CN up to $\sim$6800 K and solely atomic absorption lines above that. We have put in place a series of criteria to distinguish RCB stars from other AGB carbon-rich stars. We found 45 new RCB stars, including 30 Cold ($4000<T_{eff}<6800$ K), 14 Warm ($6800<T_{eff}<8500$ K) and one hot RCB ($T_{eff}>15000$ K). Forty of these belong to the Milky Way and five are located in the Magellanic Clouds. We also confirmed that the long lasting candidate KDM 5651 is indeed a new Magellanic RCB star, increasing the total number of Magellanic Cloud RCB stars to 30. We have increased by $\sim$50\% the total number of RCB stars known, now reaching 147. We also include a list of 14 strong RCB candidates, most certainly observed during a dust obscuration phase. From the detection efficiency and success rate so far, we estimate that there should be no more than 500 RCB stars/HdC stars in the Milky Way.


title: Probing color reconnection with underlying event observables at the LHC energies

abstract: In this work we study the underlying event (UE) activity as a function of the highest jet transverse momentum ($p_{\rm T}^{\rm jet}$) in terms of the number and summed $p_{\rm T}$ densities of charged particles in the azimuthal region transverse to the $p_{\rm T}^{\rm jet}$ direction. The UE observables normalised to the INEL multiplicity obey an approximate Koba-Nielsen-Olesen (KNO) scaling. Based on PYTHIA~8.2 simulations of $pp$ collisions at LHC energies, we show that the remaining imperfection of the KNO scaling is due to the increasing importance of multiple partonic interactions (MPI) at higher $\sqrt{s}$. Motivated by this, we studied the UE activity considering charged particles within different $p_{\rm T}$ intervals in $pp$ collisions at $\sqrt{s}=7$\,TeV. We demonstrate that the saturations of both number and summed $p_{\rm T}$ densities, which are commonly claimed, are only observed for low-\pt charged particles ($0.5<p_{\rm T}<2$\,GeV/$c$). Moreover, for the $p_{\rm T}$-integrated case ($p_{\rm T}>0.5$\,GeV/$c$) the summed $p_{\rm T}$ density is not sensitive to the variation of color reconnection (CR), however at low-p_{\rm T} it is reduced with increasing CR, whereas an opposite behaviour is found at intermediate-$p_{\rm T}$ ($2<p_{\rm T}<10$\,GeV/$c$). Finally, we show that CR produces flow-like behaviour only in the UE region and the effects are reduced with increasing $p_{\rm T}^{\rm jet}$ due to the hardering of UE. The outcomes encourage the measurement of inclusive and identified charged particle $p_{\rm T}$ spectra (over a wide range of $p_{\rm T}$) associated to UE aimed at better understanding the similarities between $pp$ and heavy-ion data discovered at the LHC.


title: Threshold dynamics for corotational wave maps

abstract: We study the dynamics of corotational wave maps from $\mathbb R^{1+2} \rightarrow \mathbb S^2$ at threshold energy. It is known that topologically trivial wave maps with energy $< 8\pi$ are global and scatter to a constant map. In this work, we prove that a corotational wave map with energy equal to $8\pi$ is globally defined and scatters in one time direction, and in the other time direction, either the map is globally defined and scatters, or the map breaks down in finite time and converges to a superposition of two harmonic maps. The latter behavior stands in stark contrast to higher equivariant wave maps with threshold energy which have been proven to be globally defined for all time. Using techniques developed in this paper, we also construct a corotational wave map with energy $= 8\pi$ which blows up in finite time. The blow-up solution we construct provides the first example of a minimal topologically trivial non-dispersing solution to the full wave map evolution.


title: To the rescue of Copenhagen interpretation

abstract: A recent paper "Single-world interpretations of quantum theory cannot be self-consistent" [arXiv:1604.07422] by D. Frauchiger and R. Renner has attracted a considerable interest of a broader physics audience and shortly elicited a number of replies. In spite of the objections that ensued, we find that significant part of the controversial initial claim has not been refuted - on the contrary, arguments presented both in the initial paper and in the replies seem to leave no longer room for the basic Copenhagen interpretation of quantum mechanics. We revisit the controversy trying to pinpoint the exact phenomenon that lies in its core, pointing out to some aspects of the problem that seem to have been mostly overlooked. Taking all conclusions into account, we here propose novel ways to approach the ostensible paradox, in which the Copenhagen interpretation is naturally preserved. Whereas arguments presented in Frauchiger and Renner paper and the replies that followed revolve over whether Wigner's friend would perceive only one or both of the possible outcomes, we offer the third, pretty obvious solution - namely, that Wigner's friend, in the well known setup, would subjectively perceive {\it none} of the potentially conflicting outcomes. As a side note, we also argue that, contrary to the often stated beliefs, hypothetical real experimental realization of Wigner's friend thought experiment actually would not provide us with any new insight into the quantum mechanics.


title: Glimpses of Space-Time Beyond the Singularities Using Supercomputers

abstract: A fundamental problem of Einstein's theory of classical general relativity is the existence of singularities such as the big bang. All known laws of physics end at these boundaries of classical space-time. Thanks to recent developments in quantum gravity, supercomputers are now playing an important role in understanding the resolution of big bang and black hole singularities. Using supercomputers, explorations of the very genesis of space and time from quantum geometry are revealing a novel picture of what lies beyond classical singularities and the new physics of the birth of our universe.


title: Pathwise solutions for fully nonlinear first- and second-order partial differential equations with multiplicative rough time dependence

abstract: The notes are an overview of part of the theory of pathwise weak solutions to two classes of scalar fully nonlinear first- and second-order degenerate parabolic partial differential equations with multiplicative rough time dependence, a special case being Brownian. These are Hamilton-Jacobi-Isaacs-Bellman and quasilinear divergence form equations including multi-dimensional scalar conservation laws. If the time dependence is `regular', the weak solutions are respectively the viscosity and entropy/kinetic solutions. The results presented here are about the wellposedness of the solutions. Some concrete applications are also discussed. The material for the first class of problems are part of the ongoing development of the theory in collaboration with P.-L. Lions.


title: Geometry of Deep Learning for Magnetic Resonance Fingerprinting

abstract: Current popular methods for Magnetic Resonance Fingerprint (MRF) recovery are bottlenecked by the heavy storage and computation requirements of a dictionary-matching (DM) step due to the growing size and complexity of the fingerprint dictionaries in multi-parametric quantitative MRI applications. In this paper we study a deep learning approach to address these shortcomings. Coupled with a dimensionality reduction first layer, the proposed MRF-Net is able to reconstruct quantitative maps by saving more than 60 times in memory and computations required for a DM baseline. Fine-grid manifold enumeration i.e. the MRF dictionary is only used for training the network and not during image reconstruction. We show that the MRF-Net provides a piece-wise affine approximation to the Bloch response manifold projection and that rather than memorizing the dictionary, the network efficiently clusters this manifold and learns a set of hierarchical matched-filters for affine regression of the NMR characteristics in each segment.


title: Discrete channel surfaces

abstract: We present a definition of discrete channel surfaces in Lie sphere geometry, which reflects several properties for smooth channel surfaces. Various sets of data, defined at vertices, on edges or on faces, are associated with a discrete channel surface that may be used to reconstruct the underlying particular discrete Legendre map. As an application we investigate isothermic discrete channel surfaces and prove a discrete version of Vessiot's Theorem.


title: Misreading EPR: Variations on an Incorrect Theme

abstract: Notwithstanding its great influence in modern physics, the EPR thought-experiment has been explained incorrectly a surprising number of times.


title: Survey on UAV Cellular Communications: Practical Aspects, Standardization Advancements, Regulation, and Security Challenges

abstract: The rapid growth of consumer Unmanned Aerial Vehicles (UAVs) is creating promising new business opportunities for cellular operators. On the one hand, UAVs can be connected to cellular networks as new types of user equipment, therefore generating significant revenues for the operators that can guarantee their stringent service requirements. On the other hand, UAVs offer the unprecedented opportunity to realize UAV-mounted flying base stations that can dynamically reposition themselves to boost coverage, spectral efficiency, and user quality of experience. Indeed, the standardization bodies are currently exploring possibilities for serving commercial UAVs with cellular networks. Industries are beginning to trial early prototypes of flying base stations or user equipments, while academia is in full swing researching mathematical and algorithmic solutions to address interesting new problems arising from flying nodes in cellular networks. In this article, we provide a comprehensive survey of all of these developments promoting smooth integration of UAVs into cellular networks. Specifically, we survey (i) the types of consumer UAVs currently available off-the-shelf, (ii) the interference issues and potential solutions addressed by standardization bodies for serving aerial users with the existing terrestrial base stations, (iii) the challenges and opportunities for assisting cellular communications with UAV-based flying relays and base stations, (iv) the ongoing prototyping and test bed activities, (v) the new regulations being developed to manage the commercial use of UAVs, and (vi) the cyber-physical security of UAV-assisted cellular communications.


title: Convergence and machine learning predictions of Monkhorst-Pack k-points and plane-wave cut-off in high-throughput DFT calculations

abstract: In this work, we developed an automatic convergence procedure for k-points and plane wave cut-off in density functional (DFT) calculations and applied it to more than 30000 materials. The computational framework for automatic convergence can take a user-defined input as a convergence criterion. For k-points, we converged energy per cell (EPC) to 0.001 eV/cell tolerance and compared the results with those obtained using an energy per atom (EPA) convergence criteria of 0.001 eV/atom. From the analysis of our results, we could relate k-point density and plane wave cut-off to material parameters such as density, the slope of bands, number of band-crossings, the maximum plane-wave cut-off used in pseudopotential generation, crystal systems, and the number of unique species in materials. We also identified some material species that would require more careful convergence than others. Moreover, we statistically investigated the dependence of k-points and cutoff on exchange-correlation functionals. We utilized all this data to train machine learning models to predict the k-point line density and plane-wave cut-off for generalized materials. This would provide users with a good starting point towards converged DFT calculations. The code used, and the converged data are available on the following websites: https://jarvis.nist.gov/, and https://github.com/usnistgov/jarvis .


title: Identifying Direct Collapse Black Hole Seeds through their Small Host Galaxies

abstract: Observations of high-redshift quasars indicate that super massive black holes (SMBHs) with masses greater than ${\sim}10^9~M_\odot$ were assembled within the first billion years after the Big Bang. It is unclear how such massive black holes formed so early. One possible explanation is that these SMBHs were seeded by `heavy' direct collapse black holes (DCBHs) with masses of $M_{\rm BH}\approx10^5~M_\odot$, but observations have not yet confirmed or refuted this scenario. In this Letter, we utilize a cosmological N-body simulation to demonstrate that before they grow roughly an order of magnitude in mass, DCBHs will have black hole mass to halo mass ratios much higher than expected for black hole remnants of Population III (Pop III) stars which have grown to the same mass (${\sim}10^6~M_\odot$). We also show that when $T_{\rm vir}\approx 10^4~{\rm K}$ halos (the potential sites of DCBH formation) merge with much larger nearby halos ($M_{\rm h} > 10^{10}~M_\odot$), they almost always orbit their larger host halos with a separation of a few kpc, which is sufficient to be spatially resolved with future X-ray and infrared telescopes. Thus, we propose that a future X-ray mission such as Lynx combined with infrared observations will be able to distinguish high-redshift DCBHs from smaller black hole seeds due to the unusually high black hole mass to stellar mass ratios of the faintest observed quasars, with inferred BH masses below ${\sim} 10^6 M_\odot$.


title: Lorenz-Mie scattering of focused light via complex focus fields: an analytic treatment

abstract: The Lorenz-Mie scattering of a wide class of focused electromagnetic fields off spherical particles is studied. The focused fields in question are constructed through complex focal displacements, leading to closed-form expressions that can exhibit several interesting physical properties, such as orbital and/or spin angular momentum, spatially-varying polarization, and a controllable degree of focusing. These fields constitute complete bases that can be considered as nonparaxial extensions of the standard Laguerre-Gauss beams and the recently proposed polynomials-of-Gaussians beams. Their analytic form turns out to lead also to closed-form expressions for their multipolar expansion. Such expansion can be used to compute the field scattered by a spherical particle and the resulting forces and torques exerted on it, for any relative position between the field's focus and the particle.


title: Token Curated Registries - A Game Theoretic Approach

abstract: Token curated registries (TCRs) have been proposed recently as an approach to create and maintain high quality lists of resources or recommendations in a decentralized manner. Applications range from maintaining registries of web domains for advertising purposes (e.g., adChain) or restaurants, consumer products, etc. The registry is maintained through a combination of candidate applications requiring a token deposit, challenges based on token staking and token-weighted votes with a redistribution of tokens occurring as a consequence of the vote. We present a simplified mathematical model of a TCR and its challenge and voting process analyze it from a game-theoretic perspective. We derive some insights into conditions with respect to the quality of a candidate under which challenges occur, and under which the outcome is reject or accept. We also show that there are conditions under which the outcome may not be entirely predictable in the sense that everyone voting for accept and everyone voting for reject could both be Nash Equilibria outcomes. For such conditions, we also explore when a particular strategy profile may be payoff dominant. We identify ways in which our modeling can be extended and also some implications of our model with respect to the composition of TCRs.


title: Coronal properties of Swift/BAT-selected Seyfert 1 AGN observed with NuSTAR

abstract: The NuSTAR observatory, with its high sensitivity in hard X-rays, has enabled detailed broadband modeling of the X-ray spectra of Active Galactic Nuclei (AGN), thereby allowing constraints to be placed on the high-energy cutoff of the X-ray coronal continuum. We investigate the spectral properties of a sample of 46 NuSTAR-observed Seyfert 1 AGN selected from the Swift/BAT 70-month hard X-ray survey. Our measurements of the high-energy cutoff of the continuum from modeling the NuSTAR X-ray spectra are used to map out the temperature - compactness ($\theta-l$) plane for AGN coronae. We find that most of the coronae lie clustered near the boundary for runaway pair production, suggesting that annihilation and pair production act to regulate the temperature of the corona. We discuss the implications of coronae whose high-energy cutoff may indicate a low coronal temperature on the heating and thermalization mechanisms in the corona.


title: Accurate Quantum Logic Gates by Spin Echo in Rydberg Atoms

abstract: Scalable quantum computing is based on realizable accurate quantum gates. For neutral atoms, it is an outstanding challenge to design a high-fidelity two-qubit entangling gate without resorting to difficult techniques like shaping laser pulses or cooling atoms to motional ground states. By using spin echo to suppress the blockade error, we propose an easily realizable controlled-phase Rydberg quantum gate of high intrinsic fidelity. In the context of spin echo, we show that the fundamental blockade error of the traditional Rydberg gate, on the order of $\epsilon\sim10^{-3}$, actually results from two `clockwise' rotations of Rabi frequencies $\bar\Omega_\pm=V\pm \sqrt{V^2+\Omega^2}$. In our `echo' sequence, such an error can be suppressed to the order of $\epsilon^2$ by adding two `anticlockwise' rotations with frequencies $-\bar\Omega_\pm$. With the blockade error effectively removed, the error caused by Rydberg state decay becomes the final fundamental limit to the gate accuracy, which in principle, can be reduced beyond the level of $10^{-5}$. Furthermore, due to the small population $\epsilon$ involved in the `echo' process, the spin-echo gate is robust against the variation of Rydberg blockade caused by the drift of the qubits, so that it can still be much more accurate than that of a traditional Rydberg gate even for qubits cooled only to the sub-mK regime.


title: Multi-finger binary search trees

abstract: We study multi-finger binary search trees (BSTs), a far-reaching extension of the classical BST model, with connections to the well-studied $k$-server problem. Finger search is a popular technique for speeding up BST operations when a query sequence has locality of reference. BSTs with multiple fingers can exploit more general regularities in the input. In this paper we consider the cost of serving a sequence of queries in an optimal (offline) BST with $k$ fingers, a powerful benchmark against which other algorithms can be measured. We show that the $k$-finger optimum can be matched by a standard dynamic BST (having a single root-finger) with an $O(\log{k})$ factor overhead. This result is tight for all $k$, improving the $O(k)$ factor implicit in earlier work. Furthermore, we describe new online BSTs that match this bound up to a $(\log{k})^{O(1)}$ factor. Previously only the "one-finger" special case was known to hold for an online BST (Iacono, Langerman, 2016; Cole et al., 2000). Splay trees, assuming their conjectured optimality (Sleator and Tarjan, 1983), would have to match our bounds for all $k$. Our online algorithms are randomized and combine techniques developed for the $k$-server problem with a multiplicative-weights scheme for learning tree metrics. To our knowledge, this is the first time when tools developed for the $k$-server problem are used in BSTs. As an application of our $k$-finger results, we show that BSTs can efficiently serve queries that are close to some recently accessed item. This is a (restricted) form of the unified property (Iacono, 2001) that was previously not known to hold for any BST algorithm, online or offline.


title: Upgrade plans and ageing studies for the CMS muon system in preparation of HL-LHC

abstract: The CMS muon system operates gas-based detectors. Upgrades of the detectors and trigger components are needed to cope with increasingly challenging conditions of the HL-LHC. New irradiation tests are performed to ensure that the muon detectors will survive the harsher conditions and operate reliably. The new CERN GIF++ (Gamma Irradiation Facility) allowed to perform aging tests. We present results in terms of system performance under large backgrounds and after accelerated aging tests. New detectors will be added to improve the performance in the critical forward region: gas electron multiplier (GEM) detectors will be installed in LS2 in the region 1.6 $<$ $|\eta|$ $<$ 2.4, aiming at suppressing background triggers while maintaining high trigger efficiency. Further enhancements are foreseen with a second GEM station and with two stations of new generation RPCs, having low resistivity electrodes. These detectors will combine tracking and triggering capabilities and can stand particle rates up to few kHz/cm$^2$. In addition, to take advantage of the pixel tracking coverage extension, a new GEM detector covering up to $|\eta|$ = 2.8. is foreseen behind the new forward calorimeter.


title: Turbulence modulation in channel flow of finite-size spheroidal particles

abstract: We study the effect of particle shape on the turbulence in suspensions of spheroidal particles at volume fraction $\phi = 10\%$ and show how the near-wall particle dynamics deeply changes with the particle aspect ratio and how this affects the global suspension behavior. The turbulence reduces with the aspect ratio of oblate particles, leading to drag reduction with respect to the single phase flow for particles with aspect ratio $\mathcal{AR}\leq1/3$, when the significant reduction in Reynolds shear stress is more than the compensation by the additional stresses, induced by the solid phase. Oblate particles are found to avoid the region close to the wall, travelling parallel to it with small angular velocities, while preferentially sampling high-speed fluid in the wall region. Prolate particles, also tend to orient parallel to the wall and avoid its vicinity. Their reluctancy to rotate around spanwise axis reduce the wall-normal velocity fluctuation of the flow and therefore the turbulence Reynolds stress similar to oblates; however, they undergo rotations in wall-parallel planes which increases the additional solid stresses due to their relatively larger angular velocities compared to the oblates. These larger additional stresses compensates for the reduction in turbulence activity and leads to a wall-drag similar to that of single-phase flows. Spheres on the other hand, form a layer close to the wall with large angular velocities in spanwise direction, which increases the turbulence activity in addition to exerting the largest solid stresses on the suspension, in comparison to the other studied shapes. Spherical particles therefore increase the wall-drag with respect to the single-phase flow.


title: Factorization of a class of composed polynomials

abstract: In this paper, we provide the degree distribution of irreducible factors of the composed polynomial $f(L(x))$ over $\mathbb F_q$, where $f(x)\in \mathbb F_q[x]$ is irreducible and $L(x)\in \mathbb F_q[x]$ is a linearized polynomial. We further provide some applications of our main result, including lower bounds for the number of irreducible factors of $f(L(x))$, constructions of high degree irreducible polynomials and the explicit factorization of $f(x^q-x)$ under certain conditions on $f(x)$.


title: Zero-diffusion Limit for Aggregation Equations over Bounded Domains

abstract: We establish the zero-diffusion limit for both continuous and discrete aggregation models over convex and bounded domains. Compared with a similar zero-diffusion limit derived in [44], our approach is different and relies on a coupling method connecting PDEs with their underlying SDEs. Moreover, our result relaxes the regularity assumptions on the interaction and external potentials and improves the convergence rate (in terms of the diffusion coefficient). The particular rate we derive is shown to be consistent with numerical computations.


title: Widespread orthorhombic fluctuations in the (Sr,Na)Fe$_2$As$_2$ family of superconductors

abstract: We report comprehensive pair distribution function measurements of the hole-doped iron-based superconductor system Sr$_{1-x}$Na$_{x}$Fe$_2$As$_2$. Structural refinements performed as a function of temperature and length scale reveal orthorhombic distortions of the instantaneous local structure across a large region of the phase diagram possessing average tetragonal symmetry, indicative of fluctuating nematicity. These nematic fluctuations are present up to high doping levels ($x \gtrsim 0.48$, near optimal superconductivity) and high temperatures (above room temperature for $x = 0$, decreasing to 150 K for $x = 0.48$), with a typical length scale of 1-3 nm. This work highlights the ubiquity of nematic fluctuations in a representative iron-based superconductor and provides important details about the evolution of these fluctuations across the phase diagram.


title: Sample Efficient Stochastic Gradient Iterative Hard Thresholding Method for Stochastic Sparse Linear Regression with Limited Attribute Observation

abstract: We develop new stochastic gradient methods for efficiently solving sparse linear regression in a partial attribute observation setting, where learners are only allowed to observe a fixed number of actively chosen attributes per example at training and prediction times. It is shown that the methods achieve essentially a sample complexity of $O(1/\varepsilon)$ to attain an error of $\varepsilon$ under a variant of restricted eigenvalue condition, and the rate has better dependency on the problem dimension than existing methods. Particularly, if the smallest magnitude of the non-zero components of the optimal solution is not too small, the rate of our proposed {\it Hybrid} algorithm can be boosted to near the minimax optimal sample complexity of {\it full information} algorithms. The core ideas are (i) efficient construction of an unbiased gradient estimator by the iterative usage of the hard thresholding operator for configuring an exploration algorithm; and (ii) an adaptive combination of the exploration and an exploitation algorithms for quickly identifying the support of the optimum and efficiently searching the optimal parameter in its support. Experimental results are presented to validate our theoretical findings and the superiority of our proposed methods.


title: Lepton non-universality in $B$ decays and fermion mass structure

abstract: We consider the possibility that the neutral-current $B$ anomalies are due to radiative corrections generated by Yukawa interactions of quarks and leptons with new vector-like quark and lepton electroweak doublets and new Standard Model singlet scalars. We show that the restricted interactions needed can result from an underlying Abelian family symmetry and that the same symmetry can give rise to an acceptable pattern of quark and charged lepton masses and mixings, providing a bridge between the non-universality observed in the B-sector and that of the fermion mass matrices. We construct two simple models, one with a single singlet scalar in which the flavour changing comes from quark and lepton mixing and one with an additional scalar in which the flavour changing can come from both fermion and scalar mixing. We show that for the case the new quarks are much heavier than the new leptons and scalars the $B$ anomalies can be due to box diagrams with couplings in the perturbative regime consistent with the bounds coming from $B_s- \bar B_s$, $K- \bar K$ and $D- \bar D$ mixing as well as other lepton family number violating processes. The new states can be dark matter candidates and, in the two scalar model with a light scalar of O(60) GeV and vector-like lepton of O(100) GeV, there can be a simultaneous explanation of the B-anomalies, the muon anomalous magnetic moment and the dark matter abundance.


title: The Maximum Size of $(k,l)$-Sum-Free Sets in Cyclic Groups

abstract: A subset $A$ of a finite abelian group $G$ is called $(k,l)$-sum-free if the sum of $k$ (not-necessarily-distinct) elements of $A$ never equals the sum of $l$ (not-necessarily-distinct) elements of $A$. We find an explicit formula for the maximum size of a $(k,l)$-sum-free subset in $G$ for all $k$ and $l$ in the case when $G$ is cyclic by proving that it suffices to consider $(k,l)$-sum-free intervals in subgroups of $G$. This simplifies and extends earlier results by Hamidoune and Plagne and by Bajnok.


title: Analyzing collisions in classical mechanics using mass-momentum diagrams

abstract: We show the value of mass-momentum diagrams for analyzing collision problems in classical mechanics in one dimension. Collisions are characterized by the coefficient of restitution and the momentum of the interacting particles both before and after the collision. All those quantities are presented in the mass-momentum diagrams without the need to do any calculations. We also show that the same diagrams can be used to investigate collisions with respect to the center-of-mass system. In this case, also, we do not need to do any calculations to obtain the momentum. Since we give an alternative way of looking at collisions in classical mechanics, this article is aimed at undergraduate-level students.


title: Special Configurations in Anchored Rectangle Packings

abstract: Given a finite set S in $[0,1]^2$ including the origin, an anchored rectangle packing is a set of non-overlapping rectangles in the unit square where each rectangle has a point of S as its left-bottom corner and contains no point of S in its interior. Allen Freedman conjectured in the 1960's one can always find an anchored rectangle packing with total area at least $1/2$. We verify the conjecture for point configurations whose relative positions belong to certain classes of permutations.


title: Energy-preserving integration of non-canonical Hamiltonian systems by continuous-stage methods

abstract: As is well known, energy is generally deemed as one of the most important physical invariants in many conservative problems and hence it is of remarkable interest to consider numerical methods which are able to preserve it. In this paper, we are concerned with the energy-preserving integration of non-canonical Hamiltonian systems by continuous-stage methods. Algebraic conditions in terms of the Butcher coefficients for ensuring the energy preservation, symmetry and quadratic-Casimir preservation respectively are presented. With the presented condition and in use of orthogonal expansion techniques, the construction of energy-preserving integrators is examined. A new class of energy-preserving integrators which is symmetric and of order $2m$ is constructed. Some numerical results are reported to verify our theoretical analysis and show the effectiveness of our new methods.


title: An Analysis of Hierarchical Text Classification Using Word Embeddings

abstract: Efficient distributed numerical word representation models (word embeddings) combined with modern machine learning algorithms have recently yielded considerable improvement on automatic document classification tasks. However, the effectiveness of such techniques has not been assessed for the hierarchical text classification (HTC) yet. This study investigates the application of those models and algorithms on this specific problem by means of experimentation and analysis. We trained classification models with prominent machine learning algorithm implementations---fastText, XGBoost, SVM, and Keras' CNN---and noticeable word embeddings generation methods---GloVe, word2vec, and fastText---with publicly available data and evaluated them with measures specifically appropriate for the hierarchical context. FastText achieved an ${}_{LCA}F_1$ of 0.893 on a single-labeled version of the RCV1 dataset. An analysis indicates that using word embeddings and its flavors is a very promising approach for HTC.


title: Multi-view Factorization AutoEncoder with Network Constraints for Multi-omic Integrative Analysis

abstract: Multi-omic data provides multiple views of the same patients. Integrative analysis of multi-omic data is crucial to elucidate the molecular underpinning of disease etiology. However, multi-omic data has the "big p, small N" problem (the number of features is large, but the number of samples is small), it is challenging to train a complicated machine learning model from the multi-omic data alone and make it generalize well. Here we propose a framework termed Multi-view Factorization AutoEncoder with network constraints to integrate multi-omic data with domain knowledge (biological interactions networks). Our framework employs deep representation learning to learn feature embeddings and patient embeddings simultaneously, enabling us to integrate feature interaction network and patient view similarity network constraints into the training objective. The whole framework is end-to-end differentiable. We applied our approach to the TCGA Pan-cancer dataset and achieved satisfactory results to predict disease progression-free interval (PFI) and patient overall survival (OS) events. Code will be made publicly available.


title: Tight MIP formulations for bounded length cyclic sequences

abstract: We study cyclic binary strings with bounds on the lengths of the intervals of consecutive ones and zeros. This is motivated by scheduling problems where such binary strings can be used to represent the state (on/off) of a machine. In this context the bounds correspond to minimum and maximum lengths of on- or off-intervals, and cyclic strings can be used to model periodic schedules. Extending results for non-cyclic strings is not straight forward. We present a non-trivial tight compact extended network flow formulation, as well as valid inequalities in the space of the state and start-up variables some of which are shown to be facet-defining. Applying a result from disjunctive programming, we also convert the extended network flow formulation into an extended formulation over the space of the state and start-up variables.