title: Tromino: Demand and DRF Aware Multi-Tenant Queue Manager for Apache Mesos Cluster

abstract: Apache Mesos, a two-level resource scheduler, provides resource sharing across multiple users in a multi-tenant cluster environment. Computational resources (i.e., CPU, memory, disk, etc. ) are distributed according to the Dominant Resource Fairness (DRF) policy. Mesos frameworks (users) receive resources based on their current usage and are responsible for scheduling their tasks within the allocation. We have observed that multiple frameworks can cause fairness imbalance in a multiuser environment. For example, a greedy framework consuming more than its fair share of resources can deny resource fairness to others. The user with the least Dominant Share is considered first by the DRF module to get its resource allocation. However, the default DRF implementation, in Apache Mesos' Master allocation module, does not consider the overall resource demands of the tasks in the queue for each user/framework. This lack of awareness can result in users without any pending task receiving more resource offers while users with a queue of pending tasks starve due to their high dominant shares. We have developed a policy-driven queue manager, Tromino, for an Apache Mesos cluster where tasks for individual frameworks can be scheduled based on each framework's overall resource demands and current resource consumption. Dominant Share and demand awareness of Tromino and scheduling based on these attributes can reduce (1) the impact of unfairness due to a framework specific configuration, and (2) unfair waiting time due to higher resource demand in a pending task queue. In the best case, Tromino can significantly reduce the average waiting time of a framework by using the proposed Demand-DRF aware policy.


title: Exploring the Fairness and Resource Distribution in an Apache Mesos Environment

abstract: Apache Mesos, a cluster-wide resource manager, is widely deployed in massive scale at several Clouds and Data Centers. Mesos aims to provide high cluster utilization via fine grained resource co-scheduling and resource fairness among multiple users through Dominant Resource Fairness (DRF) based allocation. DRF takes into account different resource types (CPU, Memory, Disk I/O) requested by each application and determines the share of each cluster resource that could be allocated to the applications. Mesos has adopted a two-level scheduling policy: (1) DRF to allocate resources to competing frameworks and (2) task level scheduling by each framework for the resources allocated during the previous step. We have conducted experiments in a local Mesos cluster when used with frameworks such as Apache Aurora, Marathon, and our own framework Scylla, to study resource fairness and cluster utilization. Experimental results show how informed decision regarding second level scheduling policy of frameworks and attributes like offer holding period, offer refusal cycle and task arrival rate can reduce unfair resource distribution. Bin-Packing scheduling policy on Scylla with Marathon can reduce unfair allocation from 38\% to 3\%. By reducing unused free resources in offers we bring down the unfairness from to 90\% to 28\%. We also show the effect of task arrival rate to reduce the unfairness from 23\% to 7\%.


title: Time-varying Autoregression with Low Rank Tensors

abstract: We present a windowed technique to learn parsimonious time-varying autoregressive models from multivariate timeseries. This unsupervised method uncovers interpretable spatiotemporal structure in data via non-smooth and non-convex optimization. In each time window, we assume the data follow a linear model parameterized by a system matrix, and we model this stack of potentially different system matrices as a low rank tensor. Because of its structure, the model is scalable to high-dimensional data and can easily incorporate priors such as smoothness over time. We find the components of the tensor using alternating minimization and prove that any stationary point of this algorithm is a local minimum. We demonstrate on a synthetic example that our method identifies the true rank of a switching linear system in the presence of noise. We illustrate our model's utility and superior scalability over extant methods when applied to several synthetic and real-world example: two types of time-varying linear systems, worm behavior, sea surface temperature, and monkey brain datasets.


title: R-mode oscillations in accreting white dwarfs in cataclysmic variables

abstract: Some cataclysmic variables show short-period ($\sim$200 to $\sim$2000 s) light variations, which are attributable to nonradial pulsations of the accreting primary white dwarf. We regard these periodic variations as r-mode (global Rossby-wave) oscillations which are presumably excited mechanically and confined into hydrogen-rich layers in the white dwarf. Often observed modulations in amplitudes and frequencies even during an observing season may be interpreted as beatings among densely distributed r-mode frequencies. Making the r-mode frequency distribution of a model to be consistent (by adjusting the rotation frequency) with the observed pulsation frequencies, we obtain the best-fit rotation rate for each case. The rotation periods thus obtained for 17 accreting white dwarfs lie between 4 and 10 minutes. Pulsation frequencies observed sometime after a dwarf-nova outburst tend to be higher than the pre-outburst values, so that the best-fit rotation rate for the post-outburst frequencies is also higher. This indicates a spin-up of the hydrogen-rich layers to occur by an enhanced accretion at the outburst. We also discuss the relation between the rotation frequencies of pulsating accreting white dwarfs and the orbital periods.


title: Group gradings on finite dimensional incidence algebras

abstract: In this work, we classify the group gradings on finite-dimensional incidence algebras over a field, where the field has characteristic zero, or the characteristic is greater than the dimension of the algebra, or the grading group is abelian. Moreover, we investigate the structure of $G$-graded $(D_1,D_2)$-bimodules, where $G$ is an abelian group, and $D_1$ and $D_2$ are the group algebra of finite subgroups of $G$. As a consequence, we can provide a more profound structure result concerning the group gradings on the incidence algebras, and we can classify their isomorphism classes of group gradings.


title: A Causality-Guided Prediction of the TED Talk Ratings from the Speech-Transcripts using Neural Networks

abstract: Automated prediction of public speaking performance enables novel systems for tutoring public speaking skills. We use the largest open repository---TED Talks---to predict the ratings provided by the online viewers. The dataset contains over 2200 talk transcripts and the associated meta information including over 5.5 million ratings from spontaneous visitors to the website. We carefully removed the bias present in the dataset (e.g., the speakers' reputations, popularity gained by publicity, etc.) by modeling the data generating process using a causal diagram. We use a word sequence based recurrent architecture and a dependency tree based recursive architecture as the neural networks for predicting the TED talk ratings. Our neural network models can predict the ratings with an average F-score of 0.77 which largely outperforms the competitive baseline method.


title: Bayesian semiparametric analysis of multivariate continuous responses, with variable selection

abstract: This article presents an approach to Bayesian semiparametric inference for Gaussian multivariate response regression. We are motivated by various small and medium dimensional problems from the physical and social sciences. The statistical challenges revolve around dealing with the unknown mean and variance functions and in particular, the correlation matrix. To tackle these problems, we have developed priors over the smooth functions and a Markov chain Monte Carlo algorithm for inference and model selection. Specifically, Dirichlet process mixtures of Gaussian distributions are used as the basis for a cluster-inducing prior over the elements of the correlation matrix. The smooth, multidimensional means and variances are represented using radial basis function expansions. The complexity of the model, in terms of variable selection and smoothness, is then controlled by spike-slab priors. A simulation study is presented, demonstrating performance as the response dimension increases. Finally, the model is fit to a number of real world datasets. An R package, scripts for replicating synthetic and real data examples, and a detailed description of the MCMC sampler are available in the supplementary materials online.


title: General-purpose quantum circuit simulator with Projected Entangled-Pair States and the quantum supremacy frontier

abstract: Recent advances on quantum computing hardware have pushed quantum computing to the verge of quantum supremacy. Random quantum circuits are outstanding candidates to demonstrate quantum supremacy, which could be implemented on a quantum device that supports nearest-neighbour gate operations on a two-dimensional configuration. Here we show that using the Projected Entangled-Pair States algorithm, a tool to study two-dimensional strongly interacting many-body quantum systems, we can realize an effective general-purpose simulator of quantum algorithms. This technique allows to quantify precisely the memory usage and the time requirements of random quantum circuits, thus showing the frontier of quantum supremacy. With this approach we can compute the full wave-function of the system, from which single amplitudes can be sampled with unit fidelity. Applying this general quantum circuit simulator we measured amplitudes for a $7\times 7$ lattice of qubits with depth $1+40+1$ and double-precision numbers in 31 minutes using less than $93$ TB memory on the Tianhe-2 supercomputer.


title: A SiO J = 5 - 4 Survey Toward Massive Star Formation Regions

abstract: We performed a survey in the SiO $J=5\rightarrow4$ line toward a sample of 199 Galactic massive star-forming regions at different evolutionary stages with the SMT 10 m and CSO 10.4 m telescopes. The sample consists of 44 infrared dark clouds (IRDCs), 86 protostellar candidates, and 69 young \HII\ regions. We detected SiO $J=5\rightarrow4$ line emission in 102 sources, with a detection rate of 57\%, 37\%, and 65\% for IRDCs, protostellar candidates, and young \HII\ regions, respectively. We find both broad line with Full Widths at Zero Power (FWZP) $>$ 20 \kms and narrow line emissons of SiO in objects at various evolutionary stages, likely associated with high-velocity shocks and low-velocity shocks, respectively. The SiO luminosities do not show apparent differences among various evolutionary stages in our sample. We find no correlation between the SiO abundance and the luminosity-to-mass ratio, indicating that the SiO abundance does not vary significantly in regions at different evolutionary stages of star formation.


title: Vortices Instead of Equilibria in MinMax Optimization: Chaos and Butterfly Effects of Online Learning in Zero-Sum Games

abstract: We establish that algorithmic experiments in zero-sum games "fail miserably" to confirm the unique, sharp prediction of maxmin equilibration. Contradicting nearly a century of economic thought that treats zero-sum games nearly axiomatically as the exemplar symbol of economic stability, we prove that no meaningful prediction can be made about the day-to-day behavior of online learning dynamics in zero-sum games. Concretely, Multiplicative Weights Updates (MWU) with constant step-size is Lyapunov chaotic in the dual (payoff) space. Simply put, let's assume that an observer asks the agents playing Matching-Pennies whether they prefer Heads or Tails (and by how much in terms of aggregate payoff so far). The range of possible answers consistent with any arbitrary small set of initial conditions blows up exponentially with time everywhere in the payoff space. This result is robust both algorithmically as well as game theoretically: 1) Algorithmic robustness: Chaos is robust to agents using any of a general sub-family of Follow-the-Regularized-Leader (FTRL) algorithms, the well known regret-minimizing dynamics, even when agents mix-and-match dynamics, use different or slowly decreasing step-sizes. 2) Game theoretic robustness: Chaos is robust to all affine variants of zero-sum games (strictly competitive games), network variants with arbitrary large number of agents and even to competitive settings beyond these. Our result is in stark contrast with the time-average convergence of online learning to (approximate) Nash equilibrium, a result widely reported as "(weak) convergence to equilibrium".


title: On the CR analogue of Frankel conjecture and a smooth representative of the first Kohn-Rossi cohomology group

abstract: In this note, we first give a criterion of pseudo-Einstein contact forms and then affirm the CR analogue of Frankel conjecture in a closed, spherical, strictly pseudoconvex CR manifold of nonnegative pseudohermitian curvature on the space of smooth representatives of the first Kohn-Rossi cohomology group. Moreover, we obtain the CR Frankel conjecture in a closed, spherical, strictly pseudoconvex CR manifold with the vanishing first Kohn-Rossi cohomology group. In particular, this conjecture holds in a spherical boundary of the Stein manifold.


title: The existence of full dimensional invariant tori for 1-dimensional nonlinear wave equation

abstract: In this paper we prove the existence and linear stability of full dimensional tori with subexponential decay for 1-dimensional nonlinear wave equation with external parameters, which relies on the method of KAM theory and the idea proposed by Bourgain \cite{BJFA2005}.


title: Anisotropic spin fluctuations in detwinned FeSe

abstract: Superconductivity in FeSe emerges from a nematic phase that breaks four-fold rotational symmetry in the iron plane. This phase may arise from orbital ordering, spin fluctuations, or hidden magnetic quadrupolar order. Here we use inelastic neutron scattering on a mosaic of single crystals of FeSe detwinned by mounting on a BaFe2As2 substrate to demonstrate that spin excitations are most intense at the antiferromagnetic wave vectors QAF = (1, 0) at low energies E = 6-11 meV in the normal state. This two-fold (C2) anisotropy is reduced at lower energies 3-5 meV, indicating a gapped four-fold (C4) mode. In the superconducting state, however, the strong nematic anisotropy is again reflected in the spin resonance (E = 3.7 meV) at QAF with incommensurate scattering around 5-6 meV. Our results highlight the extreme electronic anisotropy of the nematic phase of FeSe and are consistent with a highly anisotropic superconducting gap driven by spin fluctuations.


title: Homological dimensions of smooth crossed products

abstract: In this paper we provide upper estimates for the global projective dimensions of smooth crossed products $\mathscr{S}(G, A; \alpha)$ for $G = \mathbb{R}$ and $G = \mathbb{T}$ and a self-induced Fr\'echet-Arens-Michael algebra $A$. In order to do this, we provide a powerful generalization of methods which are used in the works of Ogneva and Helemskii.


title: Existence of solutions for critical Choquard problem with singular coefficients

abstract: In this paper, we investigate the following fractional Choquard type equation: \[ (- \Delta)_p^s\, u = \lambda\frac{|u|^{r-2}u}{|x|^\alpha}\,+\gamma \big(\int_\Omega \frac{|u|^q}{|x-y|^\mu}dy\big) |u|^{q-2}u \ \ \text{in } \Omega,\ \ u = 0 \ \text{in } \R^N \setminus \Omega, \] where $\Omega$ is a bounded domain in $\R^N$ with Lipschitz boundary, $p>1$, $0<s<1$, $N>sp$, $0\leq\alpha\leq sp$, $0<\mu<N$,$\lambda, \gamma>0$, $p\leq r\leq p^*_\alpha$, $p\leq 2q\leq 2p_{\mu,s}^*$, $p_\alpha^*=\frac{(N-\alpha)p}{N-sp}$ and $p_{\mu,s}^*=\frac{(N-\frac{\mu}{2})p}{N-sp}$ are the fractional critical Hardy-Sobolev and the critical exponents in the sense of Hardy-Littlewood-Sobolev inequality, respectively. Under some suitable assumptions, positive and sign-changing solutions are obtained.


title: Hadron Spectra Parameters within the Non-Extensive Approach

abstract: We investigate how the non-extensive approach works in high-energy physics. Transverse momentum ($p_T$) spectra of several hadrons are fitted by various non-extensive momentum distributions and by the Boltzmann--Gibbs statistics.~It is shown that some non-extensive distributions can be transferred one into another.~We find explicit hadron mass and center-of-mass energy {scaling both in the temperature and in the non-extensive parameter, $q$,} in proton--proton and heavy-ion collisions. We find that the temperature depends linearly, but the Tsallis $q$ follows a logarithmic dependence on the collision energy in proton--proton collisions. In the nucleus--nucleus collisions, on~the other hand, $T$ and $q$ correlate linearly, as was predicted in our previous work.


title: Mechanical Purcell Filters for Microwave Quantum Machines

abstract: In circuit quantum electrodynamics, measuring the state of a superconducting qubit introduces a loss channel which can enhance spontaneous emission through the Purcell effect, thus decreasing qubit lifetime. This decay can be mitigated by performing the measurement through a Purcell filter which forbids signal propagation at the qubit transition frequency. If the filter is also well-matched at the readout cavity frequency, it will protect the qubit from decoherence channels without sacrificing measurement speed. We propose and analyze design for a mechanical Purcell filter, which we also fabricate and characterize at room temperature. The filter is comprised of an array of nanomechanical resonators in thin-film lithium niobate, connected in a ladder topology, with series and parallel resonances arranged to produce a bandpass response. The modest footprint, steep band edges, and absence of cross-talk in these filters make them a novel and appealing alternative to analogous electromagnetic versions currently used in microwave quantum machines.


title: Relative Severi inequality for fibrations of maximal Albanese dimension over curves

abstract: Let $f: X \to B$ be a relatively minimal fibration of maximal Albanese dimension from a variety $X$ of dimension $n \ge 2$ to a curve $B$ defined over an algebraically closed field of characteristic zero. We prove that $K_{X/B}^n \ge 2n! \chi_f$, which was conjectured by Barja in [2]. Via the strategy outlined in [5], it also leads to a new proof of the Severi inequality for varieties of maximal Albanese dimension. Moreover, when the equality holds and $\chi_f > 0$, we prove that the general fiber $F$ of $f$ has to satisfy the Severi equality that $K_F^{n-1} = 2(n-1)! \chi(F, \omega_F)$. We also prove some sharper results of the same type under extra assumptions.


title: Pulsar tests of the gravitational Lorentz violation

abstract: Pulsars are precision celestial clocks. When being put in a binary, the ticking conveys the secret of underlying spacetime geometrodynamics. We use pulsars to test if the gravitational interaction possesses a tiny deviation from Einstein's General Relativity (GR). In the framework of Standard-Model Extension (SME), we systematically search for Lorentz-violating operators cataloged by (a) the minimal couplings of mass dimension 4, (b) the CPT symmetry of mass dimension 5, and (c) the gravitational weak equivalence principle (GWEP) of mass dimension 8. No deviation from GR was found yet.


title: Checking Robustness Against Snapshot Isolation

abstract: Transactional access to databases is an important abstraction allowing programmers to consider blocks of actions (transactions) as executing in isolation. The strongest consistency model is {\em serializability}, which ensures the atomicity abstraction of transactions executing over a sequentially consistent memory. Since ensuring serializability carries a significant penalty on availability, modern databases provide weaker consistency models, one of the most prominent being \emph{snapshot isolation}. In general, the correctness of a program relying on serializable transactions may be broken when using weaker models. However, certain programs may also be insensitive to consistency relaxations, i.e., all their properties holding under serializability are preserved even when they are executed over a weak consistent database and without additional synchronization. In this paper, we address the issue of verifying if a given program is {\em robust against snapshot isolation}, i.e., all its behaviors are serializable even if it is executed over a database ensuring snapshot isolation. We show that this verification problem is polynomial time reducible to a state reachability problem in transactional programs over a sequentially consistent shared memory. This reduction opens the door to the reuse of the classic verification technology for reasoning about weakly-consistent programs. In particular, we show that it can be used to derive a proof technique based on Lipton's reduction theory that allows to prove programs robust.


title: Generating Logical Forms from Graph Representations of Text and Entities

abstract: Structured information about entities is critical for many semantic parsing tasks. We present an approach that uses a Graph Neural Network (GNN) architecture to incorporate information about relevant entities and their relations during parsing. Combined with a decoder copy mechanism, this approach provides a conceptually simple mechanism to generate logical forms with entities. We demonstrate that this approach is competitive with the state-of-the-art across several tasks without pre-training, and outperforms existing approaches when combined with BERT pre-training.


title: Stopping time signatures for some algorithms in cryptography

abstract: We consider the normalized distribution of the overall running times of some cryptographic algorithms, and what information they reveal about the algorithms. Recent work of Deift, Menon, Olver, Pfrang, and Trogdon has shown that certain numerical algorithms applied to large random matrices exhibit a characteristic distribution of running times, which depends only on the algorithm but are independent of the choice of probability distributions for the matrices. Different algorithms often exhibit different running time distributions, and so the histograms for these running time distributions provide a time-signature for the algorithms, making it possible, in many cases, to distinguish one algorithm from another. In this paper we extend this analysis to cryptographic algorithms, and present examples of such algorithms with time-signatures that are indistinguishable, and others with time-signatures that are clearly distinct.


title: Convolutions on Spherical Images

abstract: Applying convolutional neural networks to spherical images requires particular considerations. We look to the millennia of work on cartographic map projections to provide the tools to define an optimal representation of spherical images for the convolution operation. We propose a representation for deep spherical image inference based on the icosahedral Snyder equal-area (ISEA) projection, a projection onto a geodesic grid, and show that it vastly exceeds the state-of-the-art for convolution on spherical images, improving semantic segmentation results by 12.6%.


title: Inter-plane satellite matching in dense LEO constellations

abstract: Dense constellations of Low Earth Orbit (LEO) small satellites are envisioned to make extensive use of the inter-satellite link (ISL). Within the same orbital plane, the inter-satellite distances are preserved and the links are rather stable. In contrast, the relative motion between planes makes the inter-plane ISL challenging. In a dense set-up, each spacecraft has several satellites in its coverage volume, but the time duration of each of these links is small and the maximum number of active connections is limited by the hardware. We analyze the matching problem of connecting satellites using the inter-plane ISL for unicast transmissions. We present and evaluate the performance of two solutions to the matching problem with any number of orbital planes and up to two transceivers: a heuristic solution with the aim of minimizing the total cost; and a Markovian solution to maintain the on-going connections as long as possible. The Markovian algorithm reduces the time needed to solve the matching up to 1000x and 10x with respect to the optimal solution and to the heuristic solution, respectively, without compromising the total cost. Our model includes power adaptation and optimizes the network energy consumption as the exemplary cost in the evaluations, but any other QoS-oriented KPI can be used instead.


title: Deep learning to generate in silico chemical property libraries and candidate molecules for small molecule identification in complex samples

abstract: Comprehensive and unambiguous identification of small molecules in complex samples will revolutionize our understanding of the role of metabolites in biological systems. Existing and emerging technologies have enabled measurement of chemical properties of molecules in complex mixtures and, in concert, are sensitive enough to resolve even stereoisomers. Despite these experimental advances, small molecule identification is inhibited by (i) chemical reference libraries representing <1% of known molecules, limiting the number of possible identifications, and (ii) the lack of a method to generate candidate matches directly from experimental features (i.e. without a library). To this end, we developed a variational autoencoder (VAE) to learn a continuous numerical, or latent, representation of molecular structure to expand reference libraries for small molecule identification. We extended the VAE to include a chemical property decoder, trained as a multitask network, in order to shape the latent representation such that it assembles according to desired chemical properties. The approach is unique in its application to small molecule identification, with its focus on m/z and CCS, paired with its training paradigm, which involved a cascade of transfer learning iterations. This allows the network to learn as much as possible at each stage, enabling success with progressively smaller datasets without overfitting. Once trained, the network can rapidly predict chemical properties directly from structure, as well as generate candidate structures with desired chemical properties. Additionally, the ability to generate novel molecules along manifolds, defined by chemical property analogues, positions DarkChem as highly useful in a number of application areas, including metabolomics and small molecule identification, drug discovery and design, chemical forensics, and beyond.


title: Position Paper: From Multi-Agent Pathfinding to Pipe Routing

abstract: The 2D Multi-Agent Path Finding (MAPF) problem aims at finding collision-free paths for a number of agents, from a set of start locations to a set of goal positions in a known 2D environment. MAPF has been studied in theoretical computer science, robotics, and artificial intelligence over several decades, due to its importance for robot navigation. It is currently experiencing significant scientific progress due to its relevance in automated warehousing (such as those operated by Amazon) and in other contemporary application areas. In this paper, we demonstrate that many recently developed MAPF algorithms apply more broadly than currently believed in the MAPF research community. In particular, we describe the 3D Pipe Routing (PR) problem, which aims at placing collision-free pipes from given start locations to given goal locations in a known 3D environment. The MAPF and PR problems are similar: a solution to a MAPF instance is a set of blocked cells in x-y-t space, while a solution to the corresponding PR instance is a set of blocked cells in x-y-z space. We show how to use this similarity to apply several recently developed MAPF algorithms to the PR problem, and discuss their performance on abstract PR instances. We also discuss further research necessary to tackle real-world pipe-routing instances of interest to industry today. This opens up a new direction of industrial relevance for the MAPF research community.


title: Dual-branch residual network for lung nodule segmentation

abstract: An accurate segmentation of lung nodules in computed tomography (CT) images is critical to lung cancer analysis and diagnosis. However, due to the variety of lung nodules and the similarity of visual characteristics between nodules and their surroundings, a robust segmentation of nodules becomes a challenging problem. In this study, we propose the Dual-branch Residual Network (DB-ResNet) which is a data-driven model. Our approach integrates two new schemes to improve the generalization capability of the model: 1) the proposed model can simultaneously capture multi-view and multi-scale features of different nodules in CT images; 2) we combine the features of the intensity and the convolution neural networks (CNN). We propose a pooling method, called the central intensity-pooling layer (CIP), to extract the intensity features of the center voxel of the block, and then use the CNN to obtain the convolutional features of the center voxel of the block. In addition, we designed a weighted sampling strategy based on the boundary of nodules for the selection of those voxels using the weighting score, to increase the accuracy of the model. The proposed method has been extensively evaluated on the LIDC dataset containing 986 nodules. Experimental results show that the DB-ResNet achieves superior segmentation performance with an average dice score of 82.74% on the dataset. Moreover, we compared our results with those of four radiologists on the same dataset. The comparison showed that our average dice score was 0.49% higher than that of human experts. This proves that our proposed method is as good as the experienced radiologist.


title: Strategic Bayesian Asset Allocation

abstract: Strategic asset allocation requires an investor to select stocks from a given basket of assets. The perspective of our investor is to maximize risk-adjusted alpha returns relative to a benchmark index. Historical returns are used to provide inputs into an optimization algorithm. Our approach uses Bayesian regularization to not only provide stock selection but also optimal sequential portfolio weights. By incorporating investor preferences with a number of different regularization penalties we extend the approaches of Black (1992) and Puelz (2015). We tailor standard sparse MCMC algorithms to calculate portfolio weights and perform selection. We illustrate our methodology on stock selection from the SP100 stock index and from the top fifty holdings of two hedge funds Renaissance Technologies and Viking Global. Finally, we conclude with directions for future research.


title: Evaluation of Docker Containers for Scientific Workloads in the Cloud

abstract: The HPC community is actively researching and evaluating tools to support execution of scientific applications in cloud-based environments. Among the various technologies, containers have recently gained importance as they have significantly better performance compared to full-scale virtualization, support for microservices and DevOps, and work seamlessly with workflow and orchestration tools. Docker is currently the leader in containerization technology because it offers low overhead, flexibility, portability of applications, and reproducibility. Singularity is another container solution that is of interest as it is designed specifically for scientific applications. It is important to conduct performance and feature analysis of the container technologies to understand their applicability for each application and target execution environment. This paper presents a (1) performance evaluation of Docker and Singularity on bare metal nodes in the Chameleon cloud (2) mechanism by which Docker containers can be mapped with InfiniBand hardware with RDMA communication and (3) analysis of mapping elements of parallel workloads to the containers for optimal resource management with container-ready orchestration tools. Our experiments are targeted toward application developers so that they can make informed decisions on choosing the container technologies and approaches that are suitable for their HPC workloads on cloud infrastructure. Our performance analysis shows that scientific workloads for both Docker and Singularity based containers can achieve near-native performance. Singularity is designed specifically for HPC workloads. However, Docker still has advantages over Singularity for use in clouds as it provides overlay networking and an intuitive way to run MPI applications with one container per rank for fine-grained resources allocation.


title: A novel algorithm for segmentation of leukocytes in peripheral blood

abstract: In the detection of anemia, leukemia and other blood diseases, the number and type of leukocytes are essential evaluation parameters. However, the conventional leukocyte counting method is not only quite time-consuming but also error-prone. Consequently, many automation methods are introduced for the diagnosis of medical images. It remains difficult to accurately extract related features and count the number of cells under the variable conditions such as background, staining method, staining degree, light conditions and so on. Therefore, in order to adapt to various complex situations, we consider RGB color space, HSI color space, and the linear combination of G, H and S components, and propose a fast and accurate algorithm for the segmentation of peripheral blood leukocytes in this paper. First, the nucleus of leukocyte was separated by using the stepwise averaging method. Then based on the interval-valued fuzzy sets, the cytoplasm of leukocyte was segmented by minimizing the fuzzy divergence. Next, post-processing was carried out by using the concave-convex iterative repair algorithm and the decision mechanism of candidate mask sets. Experimental results show that the proposed method outperforms the existing non-fuzzy sets methods. Among the methods based on fuzzy sets, the interval-valued fuzzy sets perform slightly better than interval-valued intuitionistic fuzzy sets and intuitionistic fuzzy sets.


title: Offloading Deadline-Constrained Cellular Traffic

abstract: In this work we study the problem of hard-deadline constrained data offloading in cellular networks. A single-Base-Station (BS) single-frequency-channel downlink system is studied where users request the same packet from the BS at the beginning of each time slot. Packets have a hard deadline of one time slot. The slot is divided into two phases. Out of those users having high channel gain allowing them to decode the packet in the first phase, one is chosen to rebroadcast it to the remaining users in the second phase. This gives the remaining users a second opportunity to potentially decode this packet before the deadline passes. By this, the BS has offloaded the packet to a `local network of users' which eliminates unnecessary BS retransmissions. The problem is modeled as a rate-adaptation and scheduling optimization problem to maximize the duration of this second phase such that each user receives a certain percentage of the packets. We show that the proposed algorithm has a polynomial complexity in the number of users with optimal performance.


title: Personality and Behavior in Role-based Online Games

abstract: Both offline and online human behaviors are affected by personality. Of special interests are online games, where players have to impersonate specific roles and their behaviors are extensively tracked by the game. In this paper, we propose to study the relationship between players' personality and game behavior in League of Legends (LoL), one of the most popular Multiplayer Online Battle Arena (MOBA) games. We use linear mixed effects (LME) models to describe relationships between players' personality traits (measured by the Five Factor Model) and two major aspects of the game: the impersonated roles and in-game actions. On the one hand, we study relationships within the game environment by modeling role attributes from match behaviors and vice versa. On the other hand, we analyze the relationship between a player's five personality traits and their game behavior by showing significant correlations between each personality trait and the set of corresponding behaviors. Our findings suggest that personality and behavior are highly entangled and provide a new perspective to understand how personality can affect behavior in role-based online games.


title: Clustering with Similarity Preserving

abstract: Graph-based clustering has shown promising performance in many tasks. A key step of graph-based approach is the similarity graph construction. In general, learning graph in kernel space can enhance clustering accuracy due to the incorporation of nonlinearity. However, most existing kernel-based graph learning mechanisms is not similarity-preserving, hence leads to sub-optimal performance. To overcome this drawback, we propose a more discriminative graph learning method which can preserve the pairwise similarities between samples in an adaptive manner for the first time. Specifically, we require the learned graph be close to a kernel matrix, which serves as a measure of similarity in raw data. Moreover, the structure is adaptively tuned so that the number of connected components of the graph is exactly equal to the number of clusters. Finally, our method unifies clustering and graph learning which can directly obtain cluster indicators from the graph itself without performing further clustering step. The effectiveness of this approach is examined on both single and multiple kernel learning scenarios in several datasets.


title: Boundary-limited and glassy-like phonon thermal conduction in EtMe$_3$Sb[Pd(dmit)$_2$]$_2$

abstract: In molecular-based quantum-spin-liquid candidate EtMe$_3$Sb[Pd(dmit)$_2$]$_2$ with two-dimensional $S$=1/2 triangular lattice, a finite residual linear term in the thermal conductivity, $\kappa_0/T\equiv\kappa/T (T \rightarrow 0)$, has been observed and attributed to the presence of itinerant gapless excitations. Here we show that the data of $\kappa$ measured in several single crystals are divided into two groups with and without the residual linear term. In the first group with finite $\kappa_0/T$, the phonon thermal conductivity $\kappa_{ph}$ is comparable to that of other organic compounds. In these crystals, the phonon mean free path $\ell_{ph}$ saturates at low temperatures, being limited by sample size. On the other hand, in the second group with zero $\kappa_0/T$, $\kappa_{ph}$ is one order of magnitude smaller than that in the first group, comparable to that of amorphous solids. In contrast to the first group, $\ell_{ph}$ shows a glassy-like non-saturating behavior at low temperatures. These results suggest that the crystals with long $\ell_{ph}$ are required to discuss the magnetic excitations by thermal conductivity measurements.


title: Effect of Boundary Constraints on the Nonlinear Flapping of Filaments Animated by Follower Forces

abstract: Elastically driven filaments subjected to animating compressive follower forces provide a synthetic way to mimic the oscillatory beating of active biological filaments such as eukaryotic cilia. The dynamics of such active filaments can, under favorable conditions, exhibit stable time-periodic responses that result due to the interplay of elastic buckling instabilities, geometric constraints, boundary conditions, and dissipation due to fluid drag. In this paper, we use a continuum elastic rod model to estimate the critical follower force required for the onset of the stable time-periodic flapping oscillations in pre-stressed rods subjected to fluid drag. The pre-stress is generated by imposing either clamped-clamped or clamped-pinned boundary constraints and the results are compared with those of clamped-free case, which is without pre-stress. We find that the critical value increases with the initial slack--that quantifies the pre-stress, and strongly depends on the type of the constraints at the boundaries. The frequency of oscillations far from the onset, however, depends primarily on the magnitude of the follower force, not on the boundary constraints. Interestingly, oscillations for the clamped-pinned case are observed only when the follower forces are directed towards the clamped end. This finding can be exploited to design a mechanical switch to initiate or quench the oscillations by reversing the direction of the follower force or altering the boundary conditions.


title: The mechanical and electrical properties of direct-spun carbon nanotube mat-epoxy composites

abstract: Composites of direct-spun carbon nanotube (CNT) mats and epoxy are manufactured and tested in order to determine their mechanical and electrical properties. The mats are spun directly from a floating catalyst, chemical vapour deposition reactor. The volume fraction of epoxy is varied widely by suitable dilution of the epoxy resin with acetone. Subsequent evaporation of the acetone, followed by a cure cycle, leads to composites of varying volume fraction of CNT, epoxy and air. The modulus, strength, electrical conductivity and piezoresistivity of the composites are measured. The CNT mats and their composites exhibit an elastic-plastic stress-strain response under uniaxial tensile loading, and the degree of anisotropy is assessed by testing specimens in 0{\deg}, 45{\deg} and 90{\deg} directions with respect to the draw direction of mat manufacture. The electrical conductivity scales linearly with CNT volume fraction, irrespective of epoxy volume fraction. In contrast, the modulus and strength depend upon both CNT and epoxy volume fractions in a non-linear manner. The macroscopic moduli of the CNT mat-epoxy composites are far below the Voigt bound based on the modulus of CNT walls and epoxy. A micromechanical model is proposed to relate the macroscopic modulus and yield strength of a CNT mat-epoxy composite to the microstructure.


title: Parallel memory-efficient all-at-once algorithms for the sparse matrix triple products in multigrid methods

abstract: Multilevel/multigrid methods is one of the most popular approaches for solving a large sparse linear system of equations, typically, arising from the discretization of partial differential equations. One critical step in the multilevel/multigrid methods is to form coarse matrices through a sequence of sparse matrix triple products. A commonly used approach for the triple products explicitly involves two steps, and during each step a sparse matrix-matrix multiplication is employed. This approach works well for many applications with a good computational efficiency, but it has a high memory overhead since some auxiliary matrices need to be temporarily stored for accomplishing the calculations. In this work, we propose two new algorithms that construct a coarse matrix with taking one pass through the input matrices without involving any auxiliary matrices for saving memory. The new approaches are referred to as "all-at-once" and "merged all-at-once", and the traditional method is denoted as "two-step". The all-at-once and the merged all-at-once algorithms are implemented based on hash tables in PETSc as part of this work with a careful consideration on the performance in terms of the compute time and the memory usage. We numerically show that the proposed algorithms and their implementations are perfectly scalable in both the compute time and the memory usage with up to 32,768 processor cores for a model problem with 27 billions of unknowns. The scalability is also demonstrated for a realistic neutron transport problem with over 2 billion unknowns on a supercomputer with 10,000 processor cores. Compared with the traditional two-step method, the all-at-once and the merged all-at-once algorithms consume much less memory for both the model problem and the realistic neutron transport problem meanwhile they are able to maintain the computational efficiency.


title: Efficient Estimation For The Cox Proportional Hazards Cure Model

abstract: While analysing time-to-event data, it is possible that a certain fraction of subjects will never experience the event of interest and they are said to be cured. When this feature of survival models is taken into account, the models are commonly referred to as cure models. In the presence of covariates, the conditional survival function of the population can be modelled by using cure model which depends on the probability of being uncured (incidence) and the conditional survival function of the uncured subjects (latency), and a combination of logistic regression and Cox proportional hazards (PH) regression is used to model the incidence and latency respectively. In this paper, we have shown the asymptotic normality of the profile likelihood estimator via asymptotic expansion of the profile likelihood and obtain the explicit form of the variance estimator with an implicit function in the profile likelihood. We have also shown the efficient score function based on projection theory and the profile likelihood score function are equal. Our contribution in this paper is that we have expressed the efficient information matrix as the variance of the profile likelihood score function. A simulation study suggests that the estimated standard errors from bootstrap samples (SMCURE package) and the profile likelihood score function (our approach) are providing similar and comparable results. The numerical result of our proposed method is also shown by using the melanoma data from SMCURE R-package (Cai et al., 2012) and we compare the results with the output obtained from SMCURE package.


title: The weakness of the pigeonhole principle under hyperarithmetical reductions

abstract: The infinite pigeonhole principle for 2-partitions ($\mathsf{RT}^1_2$) asserts the existence, for every set $A$, of an infinite subset of $A$ or of its complement. In this paper, we study the infinite pigeonhole principle from a computability-theoretic viewpoint. We prove in particular that $\mathsf{RT}^1_2$ admits strong cone avoidance for arithmetical and hyperarithmetical reductions. We also prove the existence, for every $\Delta^0_n$ set, of an infinite low${}_n$ subset of it or its complement. This answers a question of Wang. For this, we design a new notion of forcing which generalizes the first and second-jump control of Cholak, Jockusch and Slaman.


title: Interpretation of interconnection of alternative macroscopic approaches associated with the Abraham-Minkowski con-troversy using microscopic description

abstract: Interconnection of alternative equations of conservation of the momentum related with different interpretations of equations of the electromagnetic field is an essential aspect of the Abraham-Minkowski controversy, which is associated in the first place with rivalry of alternative definitions of optical momentum of the electromagnetic wave in a medium. Assuming significance of alternative equations of momentum conservation in accord with their interconnection, potential functional equivalence of alternative descriptions of optical momentum is examined in the context of the radiation pressure. In accordance with conducted examination, the material derivative of the momentum of a closed medium-field system is preserved in transformations linking alternative equations of momentum conservation, whereas the optical momentum is changed in compliance with alternative approaches. It is found that Abraham's, Minkowski's, and the total momenta are inherently linked with specific ponderomotive forces, which are not interchangeable in spite of interconnection of involved forms of field equations and equations of conservation of the momentum. Interconnection of alternative macroscopic approaches is clarified using microscopic interpretation of the energy and the momentum of the electromagnetic wave in a host medium. Contradictory references to microscopic aspects associated with alternative approaches, including assumptions in regard to manifestation of the Lorentz force on the level of macroscopic description and nature of inertial properties of electromagnetic wave in a host medium, are evaluated in this context. It leads to clarification of properties of ponderomotive forces as well as concepts of momentum and pseudomomentum of the electromagnetic wave in a host medium, which are essential for explanation of interconnection of alternative macroscopic approaches.


title: SRT22 does not imply RT22 in omega-models

abstract: We complete a 40-year old program on the computability-theoretic analysis of Ramsey's theorem, starting with Jockusch in 1972, and improving a result of Chong, Slaman and Yang in 2014. Given a set $X$, let $[X]^n$ be the collection of all $n$-element subsets of $X$. Ramsey's theorem for $n$-tuples asserts the existence, for every finite coloring of $[\omega]^n$, of an infinite set $X \subseteq \omega$ such that $[X]^n$ is monochromatic. The meta-mathematical study of Ramsey has a rich history, with several long-standing open problems and seminal theorems, including Seetapun's theorem in 1995 and Liu's theorem in 2012 about Ramsey's theorem for pairs. The remaining question about the study of Ramsey's theorem from a computational viewpoint was the relation between Ramsey's theorem for pairs ($\mathsf{RT}^2_2$) and its restriction to stable colorings ($\mathsf{SRT}^2_2$), that is, colorings admitting a limit behavior. Chong, Slaman and Yang first proved that $\mathsf{SRT}^2_2$ does not formally imply $\mathsf{RT}^2_2$ in a proof-theoretic sense, using non-standard models of reverse mathematics. In this article, we answer the open question whether this non-implication also holds within the framework of computability theory. More precisely, we construct a $\omega$-model of $\mathsf{SRT}^2_2$ which is not a model of $\mathsf{RT}^2_2$. For this, we design a new notion of effective forcing refining Mathias forcing using the notion of largeness classes.


title: Hot Corona Properties of Swift/BAT detected AGN

abstract: Using a sample of 208 broad-line active galactic nuclei (AGNs) from Swift/BAT AGN Spectroscopic Survey in ultra-hard X-ray band ($14-195$ keV), the hot corona properties are investigated, i.e. the fraction of gravitational energy dissipated in the hot corona and the hard X-ray photon index. The bolometric luminosity, \lb, is calculated from host-corrected luminosity at 5100 \AA. Virial supermassive black hole masses (SMBH, \mbh) are calculated from the $\rm H\beta$ line width and the corresponding broad line region size-luminosity empirical relation at 5100 \AA. We find a strong anti-correlation between the fraction of energy released in corona ($F_x \equiv L_{14-195 keV}/L_{\rm Edd}$) and the Eddington ratio ($ \varepsilon \equiv L_{\rm Bol}/L_{\rm Edd}$), $F_x \propto \varepsilon^{-0.60\pm 0.1}$. It is found that this fraction also has a correlation with the SMBH mass, $F_x \propto \varepsilon^{-0.74\pm 0.14} M_{\rm BH}^{-0.30\pm 0.03}$. Assuming that magnetic buoyancy and feild reconnection lead to the formation of a hot corona, our result favours the shear stress tensor being a proportion of the gas pressure. For our entire sample, it is found that the hard X-ray photon index $\Gamma$ has a weak but significant correlation with the Eddington ratio, $ \Gamma=2.17+0.21\log \varepsilon$. However, this correlation is not robust because the relation is not statistically significant for its subsample of 32 RM AGNs with relatively reliable $M_{\rm BH}$ or its subsample of 166 AGNs with single-epoch $M_{\rm BH}$. We do not find a statistically significant relation between the photon index and the Eddington ratio taking into account an additional dependence on $F_x$.


title: Quantum fractionalism: the Born rule as a consequence of the complex Pythagorean theorem

abstract: Everettian Quantum Mechanics, or the Many Worlds Interpretation, lacks an explanation for quantum probabilities. We show that the values given by the Born rule equal projection factors, describing the contraction of Lebesgue measures in orthogonal projections from the complex line of a quantum state to eigenspaces of an observable. Unit total probability corresponds to a complex Pythagorean theorem: the measure of a subset of the complex line is the sum of the measures of its projections on all eigenspaces. To show that projection factors can work as probabilities, we postulate the existence of a continuum infinity of identical quantum universes, all with the same quasi-classical worlds. In a measurement, these factors give the relative amounts of worlds with each result, which we associate to frequentist and Bayesian probabilities. This solves the probability problem of Everett's theory, allowing its preferred basis problem to be solved as well, and may help settle questions about the nature of probability.


title: A statistical mechanical analysis on the bound state solution of an energy-dependent deformed Hulth\'en potential energy

abstract: In this article, we investigate the bound state solution of the Klein Gordon equation under mixed vector and scalar coupling of an energy-dependent deformed Hulth\'en potential in D-dimensions. We obtain a transcendental equation after we impose the boundary conditions. We calculate energy spectra in four different limits and in arbitrary dimension via the Newton-Raphson method. Then, we use a statistical method, namely canonical partition function, and discuss the thermodynamic properties of the system in a comprehensive way. We find out that some of the thermodynamic properties overlap with each other, some of them do not.


title: Fully tunable and switchable coupler for photonic routing in quantum detection and modulation

abstract: Photonic routing is a key building block of many optical applications challenging its development. We report a 2$\times$2 photonic coupler with splitting ratio switchable by a low-voltage electronic signal with 10~GHz bandwidth and tens of nanoseconds latency. The coupler can operate at any splitting ratio ranging from 0:100 to 100:0 with the extinction ratio of 26 dB in optical bandwidth of 1.3 THz. We show sub-nanosecond switching between arbitrary coupling regimes including balanced 50:50 beam splitter, 0:100 switch, and a photonic tap. The core of the device is based on Mach-Zehnder interferometer in a dual-wavelength configuration allowing real-time phase lock with long-term sub-degree stability at single-photon level. Using the reported coupler, we demonstrate for the first time the perfectly balanced time-multiplexed device for photon-number-resolving detectors and also the active preparation of a photonic temporal qudit state up to four time bins. Verified long-term stable operation of the coupler at the single photon level makes it suitable for wide application range in quantum information processing and quantum optics in general.


title: Formation of disc galaxies around z~2

abstract: We present combined evolution of morphological and stellar properties of galaxies on the two sides of z=2 (2.0<z<4.0 and 1.5<z<2.0) in CDFS, with ground-based spectroscopic redshifts. We perform bulge-disc decomposition on their images in J and H filters, from the 3DHST Legacy Survey obtained using HST/WFC3. Combining morphological information with stellar properties, we provide a detailed account of the formation/growth of discs and spheroids around z~2. The fraction of 2-component (bulge+disc) systems increases from ~46% for z>2 to ~70% for z<2, compensating for the fall in population of pure discs and pure spheroids. All quiescent outliers of our full sample on the main-sequence, are 2-component systems, belonging to the lower redshift range (z<2). The doubling of stellar mass of 2-component systems and decrease in their SFR by the same factor, suggests that mechanisms involved in morphological transformations are also responsible for the quenching of their star formation activity. Interestingly, while there is substantial increase in the size (~2.5 times) and mass (~5 times) of pure discs, from z>2 to z<2, pure spheroids maintain roughly the same values. Additionally, while bulge hosting discs witness an expansion in scale length (~1.3 times), their bulge sizes as well as bulge to total light ratio see no evolution, suggesting that z~2 is pre-dominantly a disc formation period.


title: Optimal unidirectional amplification induced by optical gain in optomechanical systems

abstract: We propose a three-mode optomechanical system to realize optical nonreciprocal transmission with unidirectional amplification, where the system consists of two coupled cavities and one mechanical resonator which interacts with only one of the cavities. Additionally, the optical gain is introduced into the optomechanical cavity. It is found that for a strong optical input, the optical transmission coefficient can be greatly amplified in a particular direction and suppressed in the opposite direction. The expressions of the optimal transmission coefficient and the corresponding isolation ratio are given analytically. Our results pave a way to design high-quality nonreciprocal devices based on optomechanical systems.


title: A bound for Hall's criterion for nilpotence in semi-abelian categories

abstract: In this paper, we focus on Hall's criterion for nilpotence in semi-abelian categories, and we improve the bound of Gray's main theorem of [3,Theorem 3.4] (see Main Theorem). And this bound is best possible.


title: Effective Spin-2 Quasi-particles at Linear Dispersive Five-fold Degenerate Points with Tunable Topological Chern Numbers

abstract: In this work, which is based on spin-2 vectors and traceless spin-2 tensors, an effective Hamiltonian is constructed with a linearly dispersive five-fold degenerate point with spin-2 vector-momentum couplings. For the model without spin-2 vector-tensor coupling, the topological Chern numbers of five bands are calculated as 4, 2, 0, -2, -4. After including spin-2 vector-tensor coupling, separate topological Chern numbers are obtained for different couplings. A cubic lattice of atoms with five internal states is designed to realize two five-fold degenerate points. The Chern numbers of the bands can be changed by tuning the coupling coefficient. In this work we propose a theoretical design to obtain spin-2 quasi-particles.

categories: cond-mat.quant-gas cond-mat.mtrl-sci gr-qc hep-th physics.optics


title: Noncommutative Choquet theory

abstract: We introduce a new and extensive theory of noncommutative convexity along with a corresponding theory of noncommutative functions. We establish noncommutative analogues of the fundamental results from classical convexity theory, and apply these ideas to develop a noncommutative Choquet theory that generalizes much of classical Choquet theory. The central objects of interest in noncommutative convexity are noncommutative convex sets. The category of compact noncommutative sets is dual to the category of operator systems, and there is a robust notion of extreme point for a noncommutative convex set that is dual to Arveson's notion of boundary representation for an operator system. We identify the C*-algebra of continuous noncommutative functions on a compact noncommutative convex set as the maximal C*-algebra of the operator system of continuous noncommutative affine functions on the set. In the noncommutative setting, unital completely positive maps on this C*-algebra play the role of representing measures in the classical setting. The continuous convex noncommutative functions determine an order on the set of unital completely positive maps that is analogous to the classical Choquet order on probability measures. We characterize this order in terms of the extensions and dilations of the maps, providing a powerful new perspective on the structure of completely positive maps on operator systems. Finally, we establish a noncommutative generalization of the Choquet-Bishop-de Leeuw theorem asserting that every point in a compact noncommutative convex set has a representing map that is supported on the extreme boundary. In the separable case, we obtain a corresponding integral representation theorem.


title: A new example of a limit variety of monoids

abstract: A variety of universal algebras is called limit if it is non-finitely based but all its proper subvarieties are finitely based. Until recently, only two explicit examples of limit varieties of monoids, constructed by Jackson, were known. Recently Zhang and Luo found the third example of such a variety. In our work, one more example of a limit variety of monoids is given.


title: Bivariate Semialgebraic Splines

abstract: Semialgebraic splines are bivariate splines over meshes whose edges are arcs of algebraic curves. They were first considered by Wang, Chui, and Stiller. We compute the dimension of the space of semialgebraic splines in two extreme cases. If the polynomials defining the edges span a three-dimensional space of polynomials, then we compute the dimensions from the dimensions for a corresponding rectilinear mesh. If the mesh is sufficiently generic, we give a formula for the dimension of the spline space valid in large degree and bound how large the degree must be for the formula to hold. We also study the dimension of the spline space in examples which do not satisfy either extreme. The results are derived using commutative and homological algebra.


title: Observation of gain-pinned dissipative solitons in a microcavity laser

abstract: We demonstrate an experimental approach to create dissipative solitons in a microcavity laser. In particular, we shape the spatial gain profile of a quasi-one-dimensional microcavity laser with a nonresonant, pulsed optical pump to create spatially localised structures, called gain-pinned dissipative solitons that exist due to the balance of gain and nonlinear losses and are confined to a diffraction-limited volume. The ultrafast formation dynamics and decay of the gain-pinned solitons are probed directly, showing that they are created on a picosecond timescale, orders of magnitude faster than laser cavity solitons. All of the experimentally observed features and dynamics are reconstructed by using a standard complex Ginzburg-Landau model.


title: Suitable weak solutions for the co-rotational Beris-Edwards system in dimension three

abstract: In this paper, we establish the global existence of a suitable weak solution to the co-rotational Beris-Edwards $Q$-tensor system modeling the hydrodynamic motion of nematic liquid crystals with either Landau-De Gennes bulk potential in $\mathbb R^3$ or Ball-Majumdar bulk potential in $\mathbb{T}^3$, a system coupling the forced incompressible Navier-Stokes equation with a dissipative, parabolic system of Q-tensor $Q$ in $\mathbb R^3$, which is shown to be smooth away from a closed set $\Sigma$ whose $1$-dimensional parabolic Hausdorff measure is zero.


title: Effects of active fluctuations on energetics of a colloidal particle: superdiffusion, dissipation and entropy production

abstract: We consider a colloidal particle immersed in an active bath and derive a Smoluchowski equation that governs the dynamics of colloidal particle. We address this as active Smoluchowski equation. Our analysis based on this active Smoluchowski equation shows a short time superdiffusive behavior that strongly depends on the activity. Our model also predicts a non-monotonic dependence of mean energy dissipation against time, a signature of activity-induced dynamics. By introducing a frequency-dependent effective temperature, we show that the mean rate of entropy production is time dependent unlike a thermal system. The prime reason for these anomalies is the absence of any fluctuation-dissipation theorem for the active noise. We also comment on how microscopic details of activity can reverse the trends for mean energy dissipation and mean rate of entropy production.


title: Existence of a weak-disorder non-Fermi liquid fixed point in the hydrodynamic regime of two-dimensional nematic quantum criticality

abstract: Role of quenched randomness in metallic quantum criticality is one of the long standing problems in condensed matter physics. An aspect of the fundamental difficulties lies in the fact that such nonmagnetic disorders lead effective interactions between abundant soft modes near the Fermi surface to be drastically enhanced particularly in the diffusive regime, where the perturbative framework does not work. Here, we revisit the problem of dirty quantum critical metals in a different angle, focusing on the hydrodynamic regime instead of the diffusive regime near the non-Fermi liquid quantum critical point. More concretely, we study effects of mutual correlations between quantum critical nematic fluctuations and weak localization corrections, and show the existence of a weak-disorder non-Fermi liquid fixed point, based on the renormalization group (RG) analysis up to the two-loop order. The two-loop order RG analysis suggests that the absence of quantum coherence in two-particle composite excitations weakens the role of weak localization corrections and allows a weakly disordered non-Fermi liquid metallic state in the hydrodynamic regime of the nematic quantum critical point. Although this dirty non-Fermi liquid metallic state may not be the true infrared stable fixed point at zero temperature, expected to be characterized by the diffusive Ohmic regime, we argue that this weak-disorder non-Fermi liquid metallic fixed point would govern the region of an intermediate energy scale, described by effective hydrodynamics of nematic quantum criticality. In this respect we believe that this research would be an important step in understanding the $T-$linear electrical resistivity as a characteristic feature of non-Fermi liquids and the origin of unconventional superconductivity from effective hydrodynamics of quantum criticality.


title: The Geometry of Deep Networks: Power Diagram Subdivision

abstract: We study the geometry of deep (neural) networks (DNs) with piecewise affine and convex nonlinearities. The layers of such DNs have been shown to be {\em max-affine spline operators} (MASOs) that partition their input space and apply a region-dependent affine mapping to their input to produce their output. We demonstrate that each MASO layer's input space partitioning corresponds to a {\em power diagram} (an extension of the classical Voronoi tiling) with a number of regions that grows exponentially with respect to the number of units (neurons). We further show that a composition of MASO layers (e.g., the entire DN) produces a progressively subdivided power diagram and provide its analytical form. The subdivision process constrains the affine maps on the (exponentially many) power diagram regions to greatly reduce their complexity. For classification problems, we obtain a formula for a MASO DN's decision boundary in the input space plus a measure of its curvature that depends on the DN's nonlinearities, weights, and architecture. Numerous numerical experiments support and extend our theoretical results.


title: Predicting and Forecasting the Price of Constituents and Index of Cryptocurrency Using Machine Learning

abstract: At present, cryptocurrencies have become a global phenomenon in financial sectors as it is one of the most traded financial instruments worldwide. Cryptocurrency is not only one of the most complicated and abstruse fields among financial instruments, but it is also deemed as a perplexing problem in finance due to its high volatility. This paper makes an attempt to apply machine learning techniques on the index and constituents of cryptocurrency with a goal to predict and forecast prices thereof. In particular, the purpose of this paper is to predict and forecast the close (closing) price of the cryptocurrency index 30 and nine constituents of cryptocurrencies using machine learning algorithms and models so that, it becomes easier for people to trade these currencies. We have used several machine learning techniques and algorithms and compared the models with each other to get the best output. We believe that our work will help reduce the challenges and difficulties faced by people, who invest in cryptocurrencies. Moreover, the obtained results can play a major role in cryptocurrency portfolio management and in observing the fluctuations in the prices of constituents of cryptocurrency market. We have also compared our approach with similar state of the art works from the literature, where machine learning approaches are considered for predicting and forecasting the prices of these currencies. In the sequel, we have found that our best approach presents better and competitive results than the best works from the literature thereby advancing the state of the art. Using such prediction and forecasting methods, people can easily understand the trend and it would be even easier for them to trade in a difficult and challenging financial instrument like cryptocurrency.


title: A Variational Approach on Level sets and Linear Convergence of Variable Bregman Proximal Gradient Method for Nonconvex Optimization Problems

abstract: We develop a new variational approach on level sets aiming towards convergence rate analysis of a variable Bregman proximal gradient (VBPG) method for a broad class of nonsmooth and nonconvex optimization problems. With this new approach, we are able to extend the concepts of Bregman proximal mapping and their corresponding Bregman proximal envelops, Bregman proximal gap function to nonconvex setting. Properties of these mappings and functions are examined. An aim of this work is to provide a solid foundation on which further design and analysis of VBPG for more general nonconvex optimization problems are possible. Another aim is to provide a unified theory on linear convergence of VBPG with a particular interest towards proximal gradient methods. Centrol to our analysis for achieving the above goals is an error bound in terms of level sets and subdifferentials (level-set subdifferential error bound) along with its links to other level-set error bounds. As a consequence, we have established a number of positive results. These newly established results not only enable us to show that any accumulation of the sequence generated by VBPG is at least a critical point of the limiting subdifferential or even a critical point of the proximal subdifferential with a fixed Bregman function in each iteration, but also provide a fresh perspective that allows us to explore inner-connections among many known sufficient conditions for linear convergence of various first-order methods. Along the way, we are able to derive a number of verifiable conditions for level-set error bounds to hold, obtain linear convergence of VBPG, and derive necessary conditions and sufficient conditions for linear convergence relative to a level set for nonsmooth and nonconvex optimization problems.


title: Inference for Change Points in High Dimensional Data via Self-Normalization

abstract: This article considers change point testing and estimation for a sequence of high-dimensional data. In the case of testing for a mean shift for high-dimensional independent data, we propose a new test which is based on $U$-statistic in Chen and Qin (2010) and utilizes the self-normalization principle [Shao (2010), Shao and Zhang (2010)]. Our test targets dense alternatives in the high-dimensional setting and involves no tuning parameters. To extend to change point testing for high-dimensional time series, we introduce a trimming parameter and formulate a self-normalized test statistic with trimming to accommodate the weak temporal dependence. On the theory front, we derive the limiting distributions of self-normalized test statistics under both the null and alternatives for both independent and dependent high-dimensional data. At the core of our asymptotic theory, we obtain weak convergence of a sequential U-statistic based process for high-dimensional independent data, and weak convergence of sequential trimmed U-statistic based processes for high-dimensional linear processes, both of which are of independent interests. Additionally, we illustrate how our tests can be used in combination with wild binary segmentation to estimate the number and location of multiple change points. Numerical simulations demonstrate the competitiveness of our proposed testing and estimation procedures in comparison with several existing methods in the literature.


title: Investigating two counting methods of the holographic complexity

abstract: We investigated the distinction between two kinds of "Complexity equals Action"(CA) conjecture counting methods which are separately provided by Brown $ et\, al. $ and Lehner $et\, al.$ separately. For the late-time CA complexity growth rate, we show that the difference between two counting methods only comes from the boundary term of the segments on the horizon. However, both counting methods give the identical late-time result. Our proof is general, independent of the underlying theories of higher curvature gravity as well as the explicit stationary spacetime background. To be specific, we calculate the late-time action growth rate in SAdS black hole for F(Ricci) gravity, and show that these two methods actually give the same result. Moreover, by using the Iyer-Wald formalism, we find that the full action rate within the WDW patch can be expressed as some boundary integrations, and the final contribution only comes from the boundary on singularity. Although the definitions of the mass of black hole has been modified in F(Ricci) gravity, its late-time result has the same form with that of SAdS black hole in Einstein gravity.


title: Efficient Profile Maximum Likelihood for Universal Symmetric Property Estimation

abstract: Estimating symmetric properties of a distribution, e.g. support size, coverage, entropy, distance to uniformity, are among the most fundamental problems in algorithmic statistics. While each of these properties have been studied extensively and separate optimal estimators are known for each, in striking recent work, Acharya et al. 2016 showed that there is a single estimator that is competitive for all symmetric properties. This work proved that computing the distribution that approximately maximizes \emph{profile likelihood (PML)}, i.e. the probability of observed frequency of frequencies, and returning the value of the property on this distribution is sample competitive with respect to a broad class of estimators of symmetric properties. Further, they showed that even computing an approximation of the PML suffices to achieve such a universal plug-in estimator. Unfortunately, prior to this work there was no known polynomial time algorithm to compute an approximate PML and it was open to obtain a polynomial time universal plug-in estimator through the use of approximate PML. In this paper we provide a algorithm (in number of samples) that, given $n$ samples from a distribution, computes an approximate PML distribution up to a multiplicative error of $\exp(n^{2/3} \mathrm{poly} \log(n))$ in time nearly linear in $n$. Generalizing work of Acharya et al. 2016 on the utility of approximate PML we show that our algorithm provides a nearly linear time universal plug-in estimator for all symmetric functions up to accuracy $\epsilon = \Omega(n^{-0.166})$. Further, we show how to extend our work to provide efficient polynomial-time algorithms for computing a $d$-dimensional generalization of PML (for constant $d$) that allows for universal plug-in estimation of symmetric relationships between distributions.

categories: cs.DS cs.IT cs.LG math.IT stat.CO stat.ML


title: Far-field subwavelength resolution imaging by spatial spectrum sampling

abstract: Imaging below the diffraction limit is always a public interest because of the restricted resolution of conventional imaging systems. To beat the limit, evanescent harmonics decaying in space must participate in the imaging process. Here, we introduce the method of spatial spectrum sampling, a novel far-field superresolution imaging method for microwave and terahertz regime. Strong dispersion and momentum conservation allow the spoof surface plasmon polaritons (SSP) structure to become a sensitive probe for spatial harmonics. This enables that the spatial information of the targets including both propagating and evanescent components, can be extracted by tuning and recording SSP in the far field. Then, the subwavelength resolution is constructed by the inversed Fourier transform of the sampled spatial spectrum. Using the modified subwavelength metallic grating as the spoof plasmonic structure, a far-field resolution of 0.17 wavelength is numerically and experimentally verified, and two-dimensional imaging ability is also fully discussed. The imaging ability and flexibility can be further optimizing the SSP structures. We are confident that our working mechanism will have great potentials in the superresolution imaging applications in the microwave and terahertz frequency range


title: The P-LOOP Estimator: Covariate Adjustment for Paired Experiments

abstract: In paired experiments, participants are grouped into pairs with similar characteristics, and one observation from each pair is randomly assigned to treatment. Because of both the pairing and the randomization, the treatment and control groups should be well balanced; however, there may still be small chance imbalances. It may be possible to improve the precision of the treatment effect estimate by adjusting for these imbalances. Building on related work for completely randomized experiments, we propose the P-LOOP (paired leave-one-out potential outcomes) estimator for paired experiments. We leave out each pair and then impute its potential outcomes using any prediction algorithm. The imputation method is flexible; for example, we could use lasso or random forests. While similar methods exist for completely randomized experiments, covariate adjustment methods in paired experiments are relatively understudied. A unique trade-off exists for paired experiments, where it can be unclear whether to factor in pair assignments when making adjustments. We address this issue in the P-LOOP estimator by automatically deciding whether to account for the pairing when imputing the potential outcomes. By addressing this trade-off, the method has the potential to improve precision over existing methods.


title: Spatially Constrained Spectral Clustering Algorithms for Region Delineation

abstract: Regionalization is the task of dividing up a landscape into homogeneous patches with similar properties. Although this task has a wide range of applications, it has two notable challenges. First, it is assumed that the resulting regions are both homogeneous and spatially contiguous. Second, it is well-recognized that landscapes are hierarchical such that fine-scale regions are nested wholly within broader-scale regions. To address these two challenges, first, we develop a spatially constrained spectral clustering framework for region delineation that incorporates the tradeoff between region homogeneity and spatial contiguity. The framework uses a flexible, truncated exponential kernel to represent the spatial contiguity constraints, which is integrated with the landscape feature similarity matrix for region delineation. To address the second challenge, we extend the framework to create fine-scale regions that are nested within broader-scaled regions using a greedy, recursive bisection approach. We present a case study of a terrestrial ecology data set in the United States that compares the proposed framework with several baseline methods for regionalization. Experimental results suggest that the proposed framework for regionalization outperforms the baseline methods, especially in terms of balancing region contiguity and homogeneity, as well as creating regions of more similar size, which is often a desired trait of regions.


title: Two Dimensional Representations of the Braid Group B3 and its Burau Representation Squared

abstract: We list the irreducible two dimensional complex representations of the Braid group B3 in elementary way. Then, we make a decomposition of the square of its irreducible Burau representation.


title: Towards Safety-Aware Computing System Design in Autonomous Vehicles

abstract: Recently, autonomous driving development ignited competition among car makers and technical corporations. Low-level automation cars are already commercially available. But high automated vehicles where the vehicle drives by itself without human monitoring is still at infancy. Such autonomous vehicles (AVs) rely on the computing system in the car to to interpret the environment and make driving decisions. Therefore, computing system design is essential particularly in enhancing the attainment of driving safety. However, to our knowledge, no clear guideline exists so far regarding safety-aware AV computing system and architecture design. To understand the safety requirement of AV computing system, we performed a field study by running industrial Level-4 autonomous driving fleets in various locations, road conditions, and traffic patterns. The field study indicates that traditional computing system performance metrics, such as tail latency, average latency, maximum latency, and timeout, cannot fully satisfy the safety requirement for AV computing system design. To address this issue, we propose a `safety score' as a primary metric for measuring the level of safety in AV computing system design. Furthermore, we propose a perception latency model, which helps architects estimate the safety score of given architecture and system design without physically testing them in an AV. We demonstrate the use of our safety score and latency model, by developing and evaluating a safety-aware AV computing system computation hardware resource management scheme.


title: A Seq-to-Seq Transformer Premised Temporal Convolutional Network for Chinese Word Segmentation

abstract: The prevalent approaches of Chinese word segmentation task almost rely on the Bi-LSTM neural network. However, the methods based the Bi-LSTM have some inherent drawbacks: hard to parallel computing, little efficient in applying the Dropout method to inhibit the Overfitting and little efficient in capturing the character information at the more distant site of a long sentence for the word segmentation task. In this work, we propose a sequence-to-sequence transformer model for Chinese word segmentation, which is premised a type of convolutional neural network named temporal convolutional network. The model uses the temporal convolutional network to construct an encoder, and uses one layer of fully-connected neural network to build a decoder, and applies the Dropout method to inhibit the Overfitting, and captures the character information at the distant site of a sentence by adding the layers of the encoder, and binds Conditional Random Fields model to train parameters, and uses the Viterbi algorithm to infer the final result of the Chinese word segmentation. The experiments on traditional Chinese corpora and simplified Chinese corpora show that the performance of Chinese word segmentation of the model is equivalent to the performance of the methods based the Bi-LSTM, and the model has a tremendous growth in parallel computing than the models based the Bi-LSTM.


title: Infrafiltration Theorem and Some Inductive Sequence of Models of Generalized Second-Order Dedekind Theory of Real Numbers With Exponentially Increasing Powers

abstract: The paper is devoted to construction of some closed inductive sequence of models of the generalized second-order Dedekind theory of real numbers with exponentially increasing powers. These models are not isomorphic whereas all models of the standard second-order Dedekind theory are. The main idea in passing to generalized models is to consider instead of superstructures with the single common set-theoretical equality and the single common set-theoretical belonging superstructures with several generalized equalities and several generalized belongings for first and second orders. The basic tools for the presented construction are the infraproduct of collection of mathematical systems different from the factorized Lo\'s ultraproduct and the corresponding generalized infrafiltration theorem. As its auxiliary corollary we obtain the generalized compactness theorem for the generalized second-order language.


title: On characterizing the critical graphs for matching Ramsey numbers

abstract: Given simple graphs $H_{1},H_{2},\ldots,H_{c}$, the Ramsey number $r(H_{1},H_{2},\ldots,H_{c})$ is the smallest positive integer $n$ such that every edge-colored $K_{n}$ with $c$ colors contains a subgraph in color $i$ isomorphic to $H_{i}$ for some $i\in\{1,2,\ldots,c\}$. The critical graphs for $r(H_1,H_2,\ldots,H_c)$ are edge-colored complete graphs on $r(H_1,H_2,\ldots,H_c)-1$ vertices with $c$ colors which contain no subgraphs in color $i$ isomorphic to $H_{i}$ for any $i\in \{1,2,\ldots,c\}$. For $n_1\geq n_2\geq \ldots\geq n_c\geq 1$, Cockayne and Lorimer (The Ramsey number for stripes, {\it J.\ Austral.\ Math.\ Soc.} \textbf{19} (1975), 252--256.) showed that $r(n_{1}K_{2},n_{2}K_{2},\ldots,n_{c}K_{2})=n_{1}+1+ \sum\limits_{i=1}^c(n_{i}-1)$, in which $n_{i}K_{2}$ is a matching of size $n_{i}$. Using the Gallai-Edmonds Theorem, we characterized all the critical graphs for $r(n_{1}K_{2},n_{2}K_{2},\ldots,n_{c}K_{2})$, implying a new proof for this Ramsey number.


title: Four-term progression free sets with three-term progressions in all large subsets

abstract: This paper is mainly concerned with sets which do not contain four-term arithmetic progressions, but are still very rich in three term arithmetic progressions, in the sense that all sufficiently large subsets contain at least one such progression. We prove that there exists a positive constant $c$ and a set $A \subset \mathbb F_q^n$ which does not contain a four-term arithmetic progression, with the property that for every subset $A' \subset A$ with $|A'| \geq |A|^{1-c}$, $A'$ contains a nontrivial three term arithmetic progression. We derive this from a more general quantitative Roth-type theorem in random subsets of $\mathbb{F}_{q}^{n}$, which improves a result of Kohayakawa-Luczak-R\"odl/Tao-Vu. We also discuss a similar phenomenon over the integers, where we show that for all $\epsilon >0$, and all sufficiently large $N \in \mathbb N$, there exists a four-term progression-free set $A$ of size $N$ with the property that for every subset $A' \subset A$ with $|A'| \gg \frac{1}{(\log N)^{1-\epsilon}} \cdot N$ contains a nontrivial three term arithmetic progression. Finally, we include another application of our methods, showing that for sets in $\mathbb{F}_{q}^{n}$ or $\mathbb{Z}$ the property of "having nontrivial three-term progressions in all large subsets" is almost entirely uncorrelated with the property of "having large additive energy".


title: Acceleration of solar wind particles by traveling interplanetary shocks

abstract: The acceleration of thermal solar wind protons at spherical interplanetary shocks driven by coronal mass ejections is investigated. The solar wind velocity distribution is represented using $\kappa$-functions, which are transformed in response to simulated shock transitions in the fixed-frame flow speed, plasma number density, and temperature. These heated solar wind distributions are specified as source spectra at the shock from which particles with sufficient energy can be injected into the diffusive shock acceleration process. It is shown that for shock-accelerated spectra to display the classically expected power-law indices associated with the compression ratio, diffusion length scales must exceed the width of the compression region. The maximum attainable energies of shock-accelerated spectra are found to be limited by the transit times of interplanetary shocks, while spectra may be accelerated to higher energies in the presence of higher levels of magnetic turbulence or at faster-moving shocks. Indeed, simulations suggest fast-moving shocks are more likely to produce very high-energy particles, while strong shocks, associated with harder shock-accelerated spectra, are linked to higher intensities of energetic particles. The prior heating of the solar wind distribution is found to complement shock acceleration in reproducing the intensities of typical energetic storm particle events, especially where injection energies are high. Moreover, simulations of $\sim$0.2 to 1 MeV proton intensities are presented that naturally reproduce the observed flat energy spectra prior to shock passages. Energetic particles accelerated from the solar wind, aided by its prior heating, are shown to contribute substantially to intensities during energetic storm particle events.


title: Non-Autoregressive Neural Text-to-Speech

abstract: In this work, we propose ParaNet, a non-autoregressive seq2seq model that converts text to spectrogram. It is fully convolutional and brings 46.7 times speed-up over the lightweight Deep Voice 3 at synthesis, while obtaining reasonably good speech quality. ParaNet also produces stable alignment between text and speech on the challenging test sentences by iteratively improving the attention in a layer-by-layer manner. Furthermore, we build the parallel text-to-speech system and test various parallel neural vocoders, which can synthesize speech from text through a single feed-forward pass. We also explore a novel VAE-based approach to train the inverse autoregressive flow (IAF) based parallel vocoder from scratch, which avoids the need for distillation from a separately trained WaveNet as previous work.


title: The Mirkovic-Vilonen basis and Duistermaat-Heckman measures

abstract: Using the geometric Satake correspondence, the Mirkovic-Vilonen cycles in the affine Grasssmannian give bases for representations of a semisimple group G . We prove that these bases are "perfect", i.e. compatible with the action of the Chevelley generators of the positive half of the Lie algebra g. We compute this action in terms of intersection multiplicities in the affine Grassmannian. We prove that these bases stitch together to a basis for the algebra C[N] of regular functions on the unipotent subgroup. We compute the multiplication in this MV basis using intersection multiplicities in the Beilinson-Drinfeld Grassmannian, thus proving a conjecture of Anderson. In the third part of the paper, we define a map from C[N] to a convolution algebra of measures on the dual of the Cartan subalgebra of g. We characterize this map using the universal centralizer space of G. We prove that the measure associated to an MV basis element equals the Duistermaat-Heckman measure of the corresponding MV cycle. This leads to a proof of a conjecture of Muthiah. Finally, we use the map to measures to compare the MV basis and Lusztig's dual semicanonical basis. We formulate conjectures relating the algebraic invariants of preprojective algebra modules (which underlie the dual semicanonical basis) and geometric invariants of MV cycles. In the appendix, we use these ideas to prove that the MV basis and the dual semicanonical basis do not coincide in SL_6.


title: Random products of matrices: a dynamical point of view

abstract: We study random products of matrices in SL_2(C) from the point of view of holomorphic dynamics. For non-elementary measures with finite first moment we obtain the exponential convergence towards the stationary measure in Sobolev norm. As a consequence we obtain the exponentially fast equidistribution of forward images of points towards the stationary measure. We also give a new proof of the Central Limit Theorem for the norm cocycle under a second moment condition, originally due to Benoist-Quint, and obtain some general regularity results for stationary measures.


title: A polynomial approach to the Collatz conjecture

abstract: The Collatz conjecture is explored using polynomials based on a binary numeral system. It is shown that the degree of the polynomials, on average, decreases after a finite number of steps of the Collatz operation, which provides a weak proof of the conjecture by using induction with respect to the degree of the polynomials.


title: Spontaneous formation of polar superfluid droplets in a p-wave interacting Bose gas

abstract: We study the quantum fluctuations in the condensates of a mixture of bosonic atoms and molecules with interspecies p-wave interaction. Our analysis shows that the quantum phase of coexisting atomic and molecular condensates is unstable at the mean-field level. Unlike the mixture of s-wave interaction, the Lee-Huang-Yang correction of p-wave interaction is unexpectedly found here to exhibit an opposite sign with respect to its mean-field term above a critical particle density. This quantum correction to the mean-field energy provides a remarkable mechanism to self-stabilize the phase. The order parameter of this superfluid phase carries opposite finite momenta for the two atomic species while the molecular component is a polar condensate. Such a correlated order spontaneously breaks a rich set of global U(1) gauge, atomic spin, spatial rotation and translation, and time-reversal symmetries. For potential experimental observation, the phenomenon of anisotropic polar superfluid droplets is predicted to occur, when the particle number is kept finite.


title: Robustness Against Outliers For Deep Neural Networks By Gradient Conjugate Priors

abstract: We analyze a new robust method for the reconstruction of probability distributions of observed data in the presence of output outliers. It is based on a so-called gradient conjugate prior (GCP) network which outputs the parameters of a prior. By rigorously studying the dynamics of the GCP learning process, we derive an explicit formula for correcting the obtained variance of the marginal distribution and removing the bias caused by outliers in the training set. Assuming a Gaussian (input-dependent) ground truth distribution contaminated with a proportion $\varepsilon$ of outliers, we show that the fitted mean is in a $c e^{-1/\varepsilon}$-neighborhood of the ground truth mean and the corrected variance is in a $b\varepsilon$-neighborhood of the ground truth variance, whereas the uncorrected variance of the marginal distribution can even be infinite. We explicitly find $b$ as a function of the output of the GCP network, without a priori knowledge of the outliers (possibly input-dependent) distribution. Experiments with synthetic and real-world data sets indicate that the GCP network fitted with a standard optimizer outperforms other robust methods for regression.


title: Mixing and time-dependent CP violation in beauty at LHCb

abstract: Recent measurements of the time-dependent CP violation are presented. The decays of $B_{s}^{0}$ mesons to $J/\psi\,K^+ K^-$ and $J/\psi\,\pi^+ \pi^-$ final states are used to measure CP-violating parameters with proton-proton collision data, corresponding to an integrated luminosity of 1.9 fb$^{-1}$, collected by the LHCb detector at a centre-of-mass energy of 13 TeV in 2015 and 2016.


title: Percolation on a maximally disassortative network

abstract: We propose a maximally disassortative (MD) network model which realizes a maximally negative degree-degree correlation, and study its percolation transition to discuss the effect of a strong degree-degree correlation on the percolation critical behaviors. Using the generating function method for bipartite networks, we analytically derive the percolation threshold and the order parameter critical exponent, $\beta$. For the MD scale-free networks, whose degree distribution is $P(k) \sim k^{-\gamma}$, we show that the exponent, $\beta$, for the MD networks and corresponding uncorrelated networks are same for $\gamma>3$ but are different for $2<\gamma<3$. A strong degree-degree correlation significantly affects the percolation critical behavior in heavy-tailed scale-free networks. Our analytical results for the critical exponents are numerically confirmed by a finite-size scaling argument.


title: Nonexistence of solutions for Dirichlet problems with supercritical growth in tubular domains

abstract: We deal with Dirichlet problems of the form $$ \Delta u+f(u)=0 \mbox{ in }\Omega,\qquad u=0\ \mbox{ on }\partial \Omega $$ where $\Omega$ is a bounded domain of $\mathbb{R}^n$, $n\ge 3$, and $f$ has supercritical growth from the viewpoint of Sobolev embedding. In particular, we consider the case where $\Omega$ is a tubular domain $T_\varepsilon(\Gamma_k)$ with thickness $\varepsilon>0$ and centre $\Gamma_k$, a $k$-dimensional, smooth, compact submanifold of $\mathbb{R}^n$. Our main result concerns the case where $k=1$ and $\Gamma_k$ is contractible in itself. In this case we prove that the problem does not have nontrivial solutions for $\varepsilon>0$ small enough. When $k\ge 2$ or $\Gamma_k$ is noncontractible in itself we obtain weaker nonexistence results. Some examples show that all these results are sharp for what concerns the assumptions on $k$ and $f$.


title: A parallel & automatically tuned algorithm for multispectral image deconvolution

abstract: In the era of big data, radio astronomical image reconstruction algorithms are challenged to estimate clean images given limited computing resources and time. This article is driven by the need for large scale image reconstruction for the future Square Kilometre Array (SKA), which will become in the next decades the largest low and intermediate frequency radio telescope in the world. This work proposes a scalable wideband deconvolution algorithm called MUFFIN, which stands for "MUlti Frequency image reconstruction For radio INterferometry". MUFFIN estimates the sky images in various frequency bands given the corresponding dirty images and point spread functions. The reconstruction is achieved by minimizing a data fidelity term and joint spatial and spectral sparse analysis regularization terms. It is consequently non-parametric w.r.t. the spectral behaviour of radio sources. MUFFIN algorithm is endowed with a parallel implementation and an automatic tuning of the regularization parameters, making it scalable and well suited for big data applications such as SKA. Comparisons between MUFFIN and the state-of-the-art wideband reconstruction algorithm are provided.