title: Minimally Intersective Polynomials with Arbitrarily Many Quadratic Factors

abstract: Given a natural number $n \geq 4$ we show that there exists infinitely many polynomials $f_{n}(x):= \prod_{i=1}^{n} (x^{2} - a_{i})$ such that (i) $f_{n}(x)$ has a root modulo every positive integer, (ii) $f_{n}(x)$ has no rational roots, and (iii) every proper divisor of $f_{n}(x)$ fails to have root modulo some positive integer. We exhibit a process to explicitly construct such $f_{n}$ and this process demonstrates that the set of natural numbers $a_{n}$, such that the polynomial $f_{n}(x):= \prod_{i=1}^{n} (x^{2} - a_{i})$ satisfies the properties (i), (ii) and (iii), is of positive asymptotic density in $\mathbb{N}$.


title: Entity-level Factual Consistency of Abstractive Text Summarization

abstract: A key challenge for abstractive summarization is ensuring factual consistency of the generated summary with respect to the original document. For example, state-of-the-art models trained on existing datasets exhibit entity hallucination, generating names of entities that are not present in the source document. We propose a set of new metrics to quantify the entity-level factual consistency of generated summaries and we show that the entity hallucination problem can be alleviated by simply filtering the training data. In addition, we propose a summary-worthy entity classification task to the training process as well as a joint entity and summary generation approach, which yield further improvements in entity level metrics.


title: An abstract critical point theorem with applications to elliptic problems with combined nonlinearities

abstract: We prove an abstract critical point theorem based on a cohomological index theory that produces pairs of nontrivial critical points with nontrivial higher critical groups. This theorem yields pairs of nontrivial solutions that are neither local minimizers nor of mountain pass type for problems with combined nonlinearities. Applications are given to subcritical and critical $p$-Laplacian problems, Kirchhoff type nonlocal problems, and critical fractional $p$-Laplacian problems.


title: Efficient Carpooling and Toll Pricing for Autonomous Transportation

abstract: In this paper, we address the existence and computation of competitive equilibrium in the transportation market for autonomous carpooling first proposed by [Ostrovsky and Schwarz, 2019]. At equilibrium, the market organizes carpooled trips over a transportation network in a socially optimal manner and sets the corresponding payments for individual riders and toll prices on edges. The market outcome ensures individual rationality, stability of carpooled trips, budget balance, and market clearing properties under heterogeneous rider preferences. We show that the question of market's existence can be resolved by proving the existence of an integer optimal solution of a linear programming problem. We characterize conditions on the network topology and riders' disutility for carpooling under which a market equilibrium can be computed in polynomial time. This characterization relies on ideas from the theory of combinatorial auctions and minimum cost network flow problem. Finally, we characterize a market equilibrium that achieves strategyproofness and maximizes welfare of individual riders.


title: Densely Nested Top-Down Flows for Salient Object Detection

abstract: With the goal of identifying pixel-wise salient object regions from each input image, salient object detection (SOD) has been receiving great attention in recent years. One kind of mainstream SOD methods is formed by a bottom-up feature encoding procedure and a top-down information decoding procedure. While numerous approaches have explored the bottom-up feature extraction for this task, the design on top-down flows still remains under-studied. To this end, this paper revisits the role of top-down modeling in salient object detection and designs a novel densely nested top-down flows (DNTDF)-based framework. In every stage of DNTDF, features from higher levels are read in via the progressive compression shortcut paths (PCSP). The notable characteristics of our proposed method are as follows. 1) The propagation of high-level features which usually have relatively strong semantic information is enhanced in the decoding procedure; 2) With the help of PCSP, the gradient vanishing issues caused by non-linear operations in top-down information flows can be alleviated; 3) Thanks to the full exploration of high-level features, the decoding process of our method is relatively memory efficient compared against those of existing methods. Integrating DNTDF with EfficientNet, we construct a highly light-weighted SOD model, with very low computational complexity. To demonstrate the effectiveness of the proposed model, comprehensive experiments are conducted on six widely-used benchmark datasets. The comparisons to the most state-of-the-art methods as well as the carefully-designed baseline models verify our insights on the top-down flow modeling for SOD. The code of this paper is available at https://github.com/new-stone-object/DNTD.


title: On the Mathematics of Swarming: Emergent Behavior in Alignment Dynamics

abstract: We overview recent developments in the study of alignment hydrodynamics, driven by a general class of symmetric communication kernels. A main question of interest is to characterize the emergent behavior of such systems, which we quantify in terms of the spectral gap of a weighted Laplacian associated with the alignment operator. Our spectral analysis of energy fluctuation covers both long-range and short-range kernels and does not require thermal equilibrium (no closure for the pressure). In particular, in the prototypical case of metric-based short-range kernels, the spectral gap admits a lower-bound expressed in terms of the discrete Fourier coefficients of the radial kernel, which enables us to quantify an emerging flocking behavior for non-vacuous solutions. These large-time behavior results apply as long as the solutions remain smooth. It is known that global smooth solutions exist in one and two spatial dimensions, subject to sub-critical initial data. We settle the question for arbitrary dimension, obtaining non-trivial initial threshold conditions which guarantee existence of multiD global smooth solutions.


title: Tunable magnetism in ferroelectric {\alpha}-In2Se3 by hole-doping

abstract: Two-dimensional (2D) multiferroics attract intensive investigations because of underlying science and their potential applications. Although many 2D systems have been observed/predicted to be ferroelectric or ferromagnetic, 2D materials with both ferroic properties are still scarce. By using first-principles calculations, we predict that hole-doping can induce robust ferromagnetism in 2D ferroelectric {\alpha}-In2Se3 due to its unique flat band structure, and the Curie temperature (TC) can be much higher than room temperature. Moreover, the doping concentration, strain, and number of layers can effectively modulate the magnetic moment and the TC of the material. Interestingly, strong magnetoelectric coupling is found at the surface of hole doped multilayer {\alpha}-In2Se3, which allows non-volatile electric control of magnetization. Our work provides a feasible approach for designing/searching 2D multiferroics with great potential in future device applications, such as memory devices and sensors.


title: From Extreme Multi-label to Multi-class: A Hierarchical Approach for Automated ICD-10 Coding Using Phrase-level Attention

abstract: Clinical coding is the task of assigning a set of alphanumeric codes, referred to as ICD (International Classification of Diseases), to a medical event based on the context captured in a clinical narrative. The latest version of ICD, ICD-10, includes more than 70,000 codes. As this is a labor-intensive and error-prone task, automatic ICD coding of medical reports using machine learning has gained significant interest in the last decade. Existing literature has modeled this problem as a multi-label task. Nevertheless, such multi-label approach is challenging due to the extremely large label set size. Furthermore, the interpretability of the predictions is essential for the endusers (e.g., healthcare providers and insurance companies). In this paper, we propose a novel approach for automatic ICD coding by reformulating the extreme multi-label problem into a simpler multi-class problem using a hierarchical solution. We made this approach viable through extensive data collection to acquire phrase-level human coder annotations to supervise our models on learning the specific relations between the input text and predicted ICD codes. Our approach employs two independently trained networks, the sentence tagger and the ICD classifier, stacked hierarchically to predict a codeset for a medical report. The sentence tagger identifies focus sentences containing a medical event or concept relevant to an ICD coding. Using a supervised attention mechanism, the ICD classifier then assigns each focus sentence with an ICD code. The proposed approach outperforms strong baselines by large margins of 23% in subset accuracy, 18% in micro-F1, and 15% in instance based F-1. With our proposed approach, interpretability is achieved not through implicitly learned attention scores but by attributing each prediction to a particular sentence and words selected by human coders.


title: Multi-Agent Reinforcement Learning of 3D Furniture Layout Simulation in Indoor Graphics Scenes

abstract: In the industrial interior design process, professional designers plan the furniture layout to achieve a satisfactory 3D design for selling. In this paper, we explore the interior graphics scenes design task as a Markov decision process (MDP) in 3D simulation, which is solved by multi-agent reinforcement learning. The goal is to produce furniture layout in the 3D simulation of the indoor graphics scenes. In particular, we firstly transform the 3D interior graphic scenes into two 2D simulated scenes. We then design the simulated environment and apply two reinforcement learning agents to learn the optimal 3D layout for the MDP formulation in a cooperative way. We conduct our experiments on a large-scale real-world interior layout dataset that contains industrial designs from professional designers. Our numerical results demonstrate that the proposed model yields higher-quality layouts as compared with the state-of-art model. The developed simulator and codes are available at \url{https://github.com/CODE-SUBMIT/simulator2}.


title: Time-lag between disk and corona radiation leads to hysteresis effect observed in black-hole X-ray binary MAXI J1348-630

abstract: Accretion is an essential physical process in black-hole X-ray binaries (BHXRBs) and active galactic nuclei. The properties of accretion flows and their radiation were originally considered to be uniquely determined by the mass accretion rate of the disk; however, the ``hysteresis effect'' observed during outbursts of nearly all BHXRBs seriously challenges this paradigm. The hysteresis effect is referred to that the hard-to-soft state transition in the fast-rise stage occurs at much higher luminosity than the soft-to-hard state transition in the slow-decay stage. That is, the same source can show different spectral/temporal properties at the same luminosity. Phenomenologically, this effect is also represented as the so-called ``q''-shaped hardness-intensity diagram, which has been proposed as a unified scene for BHXRBs. However, there is still a lack of quantitative theoretical interpretation and observational understanding on the ``q''-diagram. Here, we present a detailed time-lag analysis of a recently found BHXRB, MAXI J1348-630, intensively monitored by Insight-HXMT over a broad energy band (1-150 keV). We find the first observational evidence that the observed time-lag between radiations of the accretion disk and the corona leads naturally to the hysteresis effect and the ``q''-diagram. Moreover, complemented by the quasi-simultaneous Swift data, we achieve a panorama of the accretion flow: the hard X-ray outburst from the corona heats and subsequently induces the optical brightening in the outer disk with nearly no lag; thereafter, the enhanced accretion in the outer disk propagates inward, generating the delayed soft X-ray outburst at the viscous timescale of ~ 8-12 days.


title: Understanding algorithmic collusion with experience replay

abstract: In an infinitely repeated pricing game, pricing algorithms based on artificial intelligence (Q-learning) may consistently learn to charge supra-competitive prices even without communication. Although concerns on algorithmic collusion have arisen, little is known on underlying factors. In this work, we experimentally analyze the dynamics of algorithms with three variants of experience replay. Algorithmic collusion still has roots in human preferences. Randomizing experience yields prices close to the static Bertrand equilibrium and higher prices are easily restored by favoring the latest experience. Moreover, relative performance concerns also stabilize the collusion. Finally, we investigate the scenarios with heterogeneous agents and test robustness on various factors.


title: Learning Fair Representations for Recommendation: A Graph-based Perspective

abstract: As a key application of artificial intelligence, recommender systems are among the most pervasive computer aided systems to help users find potential items of interests. Recently, researchers paid considerable attention to fairness issues for artificial intelligence applications. Most of these approaches assumed independence of instances, and designed sophisticated models to eliminate the sensitive information to facilitate fairness. However, recommender systems differ greatly from these approaches as users and items naturally form a user-item bipartite graph, and are collaboratively correlated in the graph structure. In this paper, we propose a novel graph based technique for ensuring fairness of any recommendation models. Here, the fairness requirements refer to not exposing sensitive feature set in the user modeling process. Specifically, given the original embeddings from any recommendation models, we learn a composition of filters that transform each user's and each item's original embeddings into a filtered embedding space based on the sensitive feature set. For each user, this transformation is achieved under the adversarial learning of a user-centric graph, in order to obfuscate each sensitive feature between both the filtered user embedding and the sub graph structures of this user. Finally, extensive experimental results clearly show the effectiveness of our proposed model for fair recommendation. We publish the source code at https://github.com/newlei/FairGo.


title: Closing the loop: nonlinear Taylor vortex flow through the lens of resolvent analysis

abstract: We present an optimization-based method to efficiently calculate accurate nonlinear models of Taylor vortex flow. We use the resolvent formulation of McKeon & Sharma (2010) to model these Taylor vortex solutions by treating the nonlinearity not as an inherent part of the governing equations but rather as a triadic constraint which must be satisfied by the model solution. We exploit the low rank linear dynamics of the system to calculate an efficient basis for our solution, the coefficients of which are then calculated through an optimization problem where the cost function to be minimized is the triadic consistency of the solution with itself as well as with the input mean flow. Our approach constitutes, what is to the best of our knowledge, the first fully nonlinear and self-sustaining, resolvent-based model described in the literature. We compare our results to direct numerical simulation of Taylor Couette flow at up to five times the critical Reynolds number, and show that our model accurately captures the structure of the flow. Additionally, we find that as the Reynolds number increases the flow undergoes a fundamental transition from a classical weakly nonlinear regime, where the forcing cascade is strictly down scale, to a fully nonlinear regime characterized by the emergence of an inverse (up scale) forcing cascade. Triadic contributions from the inverse and traditional cascade destructively interfere implying that the accurate modeling of a certain Fourier mode requires knowledge of its immediate harmonic and sub-harmonic. We show analytically that this finding is a direct consequence of the structure of the quadratic nonlinearity of the governing equations formulated in Fourier space. Finally, we show that using our model solution as an initial condition to a higher Reynolds number DNS significantly reduces the time to convergence.


title: Improved Point Transformation Methods For Self-Supervised Depth Prediction

abstract: Given stereo or egomotion image pairs, a popular and successful method for unsupervised learning of monocular depth estimation is to measure the quality of image reconstructions resulting from the learned depth predictions. Continued research has improved the overall approach in recent years, yet the common framework still suffers from several important limitations, particularly when dealing with points occluded after transformation to a novel viewpoint. While prior work has addressed this problem heuristically, this paper introduces a z-buffering algorithm that correctly and efficiently handles occluded points. Because our algorithm is implemented with operators typical of machine learning libraries, it can be incorporated into any existing unsupervised depth learning framework with automatic support for differentiation. Additionally, because points having negative depth after transformation often signify erroneously shallow depth predictions, we introduce a loss function to penalize this undesirable behavior explicitly. Experimental results on the KITTI data set show that the z-buffer and negative depth loss both improve the performance of a state of the art depth-prediction network.


title: An Expansion Formula for Decorated Super-Teichm\"uller Spaces

abstract: Motivated by the definition of super-Teichm\"uller spaces, and Penner-Zeitlin's recent extension of this definition to decorated super-Teichm\"uller space, as examples of super Riemann surfaces, we use the super Ptolemy relations to obtain formulas for super $\lambda$-lengths associated to arcs in a bordered surface. In the special case of a disk, we are able to give combinatorial expansion formulas for the super $\lambda$-lengths associated to diagonals of a polygon in the spirit of Ralf Schiffler's $T$-path formulas for type $A$ cluster algebras. We further connect our formulas to the super-friezes of Morier-Genoud, Ovsienko, and Tabachnikov, and obtain partial progress towards defining super cluster algebras of type $A_n$. In particular, following Penner-Zeitlin, we are able to get formulas (up to signs) for the $\mu$-invariants associated to triangles in a triangulated polygon, and explain how these provide a step towards understanding odd variables of a super cluster algebra.


title: Stochastic Spatio-Temporal Optimization for Control and Co-Design of Systems in Robotics and Applied Physics

abstract: Correlated with the trend of increasing degrees of freedom in robotic systems is a similar trend of rising interest in Spatio-Temporal systems described by Partial Differential Equations (PDEs) among the robotics and control communities. These systems often exhibit dramatic under-actuation, high dimensionality, bifurcations, and multimodal instabilities. Their control represents many of the current-day challenges facing the robotics and automation communities. Not only are these systems challenging to control, but the design of their actuation is an NP-hard problem on its own. Recent methods either discretize the space before optimization, or apply tools from linear systems theory under restrictive linearity assumptions in order to arrive at a control solution. This manuscript provides a novel sampling-based stochastic optimization framework based entirely in Hilbert spaces suitable for the general class of \textit{semi-linear} SPDEs which describes many systems in robotics and applied physics. This framework is utilized for simultaneous policy optimization and actuator co-design optimization. The resulting algorithm is based on variational optimization, and performs joint episodic optimization of the feedback control law and the actuation design over episodes. We study first and second order systems, and in doing so, extend several results to the case of second order SPDEs. Finally, we demonstrate the efficacy of the proposed approach with several simulated experiments on a variety of SPDEs in robotics and applied physics including an infinite degree-of-freedom soft robotic manipulator.


title: The warm-up phase in massive star-forming cores around RCW 120

abstract: We study molecular emission in a massive condensation at the border of the HII region RCW 120, paying particular attention to the Core 1 and Core 2 objects, the most massive fragments of the condensation found previously by ALMA. The latter fragment was previously suggested to host a high-mass analogue of Class 0 young stellar object. We present spectra of molecular emission in the 1 mm range made with the APEX telescope. We detect CH$_3$OH and C$^{34}$S lines in Core 1 and Core 2. The CH$_3$CN series and the SO$_2$ lines are only found in Core 2. We estimate gas physical parameters using methanol lines and obtain gas temperature less than 100 K in both regions. Molecular hydrogen number density in Core 2 is in the range of $10^5-10^7$ cm$^{-3}$ and is more uncertain in Core 1. However, the detection of the CH$_3$CN lines corresponding to highly excited transitions ($E_{\rm u}> 400$~K) in Core~2 indicates that the region contains hot gas, while the abundances of CH$_3$OH, CS, SO$_2$ and CH$_3$CN are quite low for a hot core stage. We propose that Core 2 is in the warm-up phase prior to the establishing of the hot gas chemistry. We suggest that Core 2 is in the beginning of the hot core stage. There are no detected CH$_3$CN lines in Core 1, therefore, it might be on an even less evolved evolutionary stage.


title: Characterising the Non-Equilibrium Dynamics of a Neural Cell

abstract: We examine the dynamical evolution of the state of a neurone, with particular care to the non-equilibrium nature of the forces influencing its movement in state space. We combine non-equilibrium statistical mechanics and dynamical systems theory to characterise the nature of the neural resting state, and its relationship to firing. The stereotypical shape of the action potential arises from this model, as well as bursting dynamics, and the non-equilibrium phase transition from resting to spiking. Geometric properties of the system are discussed, such as the birth and shape of the neural limit cycle, which provide a complementary understanding of these dynamics. This provides a multiscale model of the neural cell, from molecules to spikes, and explains various phenomena in a unified manner. Some more general notions for damped oscillators, birth-death processes, and stationary non-equilibrium systems are included.

categories: q-bio.NC cond-mat.stat-mech math.DS nlin.AO


title: Effect of grid sensitivity on the performance of wall adapting SGS models for LES of swirling and separating reattaching flows

abstract: The present study assesses the performance of the Wall Adapting SGS models along with the Dynamic Smagorinsky model for flows involving separation, reattachment and swirl. Due to the simple geometry and wide application in a variety of engineering systems, the Backward-Facing Step (BFS) geometry and Confined Swirling Flow (CSF) geometry are invoked in the present case. The calculation of the SGS stresses employs three models, namely, the Dynamic Smagorinsky model, the Wall Adapting Local Eddy viscosity (WALE) model and the Dynamic WALE model. For studying the effect of the grid sensitivity, the simulations are performed over two sets of grids with different resolutions based on the non-dimensional wall distance parameter ( y+ ). Grids corresponding to y+ = 70 and y+ = 20 are employed for the subsonic flow over the BFS while grids corresponding to y+ = 40 and y+ = 20 are employed for supersonic flow over the BFS and for confined swirling flow geometry. The validation against the experimental results includes the mean flow fields and the turbulent stresses obtained for each case. The results reveal that for the fine grid (y+ = 20), the near wall eddy viscosity profile for the WALE model is better than both the Dynamic WALE and the Dynamic Smagorinsky model. The difference between the predictions of the coarse and fine grids for Dynamic Smagorinsky and the WALE model is high whereas, the Dynamic WALE model is almost insensitive to the grid resolutions considered for the present case. The mean velocity and pressure values as well as the turbulent quantities predicted by the Dynamic WALE model are closest to the experimental values for all the cases.


title: Some continuous and discontinuous Galerkin methods and structure preservation for incompressible flows

abstract: In this paper, we present consistent and inconsistent discontinuous Galerkin methods for incompressible Euler and Navier-Stokes equations with the kinematic pressure, Bernoulli function and EMAC function. Semi- and fully discrete energy stability of the proposed dG methods are proved in a unified fashion. Conservation of total energy, linear and angular momentum is discussed with both central and upwind fluxes. Numerical experiments are presented to demonstrate our findings and compare our schemes with conventional schemes in the literature in both unsteady and steady problems. Numerical results show that global conservation of the physical quantities may not be enough to demonstrate the performance of the schemes, and our schemes are competitive and able to capture essential physical features in several benchmark problems.


title: Classically Verifiable NIZK for QMA with Preprocessing

abstract: We propose three constructions of classically verifiable non-interactive zero-knowledge proofs and arguments (CV-NIZK) for QMA in various preprocessing models. - We construct a CV-NIZK for QMA in the quantum secret parameter model where a trusted setup sends a quantum proving key to the prover and a classical verification key to the verifier. It is information theoretically sound and zero-knowledge. - Assuming the quantum hardness of the learning with errors problem, we construct a CV-NIZK for QMA in a model where a trusted party generates a CRS and the verifier sends an instance-independent quantum message to the prover as preprocessing. This model is the same as one considered in the recent work by Coladangelo, Vidick, and Zhang (CRYPTO '20). Our construction has the so-called dual-mode property, which means that there are two computationally indistinguishable modes of generating CRS, and we have information theoretical soundness in one mode and information theoretical zero-knowledge property in the other. This answers an open problem left by Coladangelo et al, which is to achieve either of soundness or zero-knowledge information theoretically. To the best of our knowledge, ours is the first dual-mode NIZK for QMA in any kind of model. - We construct a CV-NIZK for QMA with quantum preprocessing in the quantum random oracle model. This quantum preprocessing is the one where the verifier sends a random Pauli-basis states to the prover. Our construction uses the Fiat-Shamir transformation. The quantum preprocessing can be replaced with the setup that distributes Bell pairs among the prover and the verifier, and therefore we solve the open problem by Broadbent and Grilo (FOCS '20) about the possibility of NIZK for QMA in the shared Bell pair model via the Fiat-Shamir transformation.


title: An Enhanced Adversarial Network with Combined Latent Features for Spatio-Temporal Facial Affect Estimation in the Wild

abstract: Affective Computing has recently attracted the attention of the research community, due to its numerous applications in diverse areas. In this context, the emergence of video-based data allows to enrich the widely used spatial features with the inclusion of temporal information. However, such spatio-temporal modelling often results in very high-dimensional feature spaces and large volumes of data, making training difficult and time consuming. This paper addresses these shortcomings by proposing a novel model that efficiently extracts both spatial and temporal features of the data by means of its enhanced temporal modelling based on latent features. Our proposed model consists of three major networks, coined Generator, Discriminator, and Combiner, which are trained in an adversarial setting combined with curriculum learning to enable our adaptive attention modules. In our experiments, we show the effectiveness of our approach by reporting our competitive results on both the AFEW-VA and SEWA datasets, suggesting that temporal modelling improves the affect estimates both in qualitative and quantitative terms. Furthermore, we find that the inclusion of attention mechanisms leads to the highest accuracy improvements, as its weights seem to correlate well with the appearance of facial movements, both in terms of temporal localisation and intensity. Finally, we observe the sequence length of around 160\,ms to be the optimum one for temporal modelling, which is consistent with other relevant findings utilising similar lengths.


title: Investigation of flow structures in a turbulent separating flow using hybrid RANS-LES model

abstract: Hybrid Reynolds-averaged Navier Stokes large eddy simulation (RANS LES) methods have become popular for simulation of massively separated flows at high Reynolds numbers due to their reduced computational cost and good accuracy. The current study aims to examine the performance of LES and hybrid RANS-LES model for a given grid resolution. For better assessment and contrast of model performance, both mean and instantaneous flow fields have been investigated. For studying instantaneous flow, proper orthogonal decomposition has been used. Current analysis shows that hybrid RANS-LES is capable of achieving similar accuracy in prediction of both mean and instantaneous flow fields at a very coarse grid as compared to LES. Focusing mostly on the practical applications of computation, most of the attention has been given to the prediction of one-point flow statistics and little consideration has been put to two-point statistics. Here, two-point statistics has been considered using POD to investigate unsteady turbulent flow.


title: Layer-Dependent Electronic and Optical Properties of 2D Black Phosphorus: Fundamentals and Engineering

abstract: In 2D materials, the quantum confinement and van der Waals-type interlayer interactions largely govern the fundamental electronic and optical properties, and the dielectric screening plays a dominant role in the excitonic properties. This suggests strongly layer-dependent properties, and a central topic is to characterize and control the interlayer interactions in 2D materials and heterostructures. Black phosphorus is an emerging 2D semiconductor with unusually strong interlayer interactions and widely tunable direct bandgaps from the monolayer to the bulk, offering us an ideal platform to probe the layer-dependent properties and the crossover from 2D to 3D (i.e., the scaling effects). In this review, we present a comprehensive and thorough summary of the fundamental physical properties of black phosphorus, with a special focus on the layer-dependence characters, including the electronic band structures, optical absorption and photoluminescence, and excitonic properties, as well as the band structure engineering by means of electrical gating, strain, and electrochemical intercalation. Finally, we give an outlook for the future research.


title: On the Optimal Duration of Spectrum Leases in Exclusive License Markets with Stochastic Demand

abstract: This paper addresses the following question which is of interest in designing efficient exclusive-use spectrum licenses sold through spectrum auctions. Given a system model in which customer demand, revenue, and bids of wireless operators are characterized by stochastic processes and an operator is interested in joining the market only if its expected revenue is above a threshold and the lease duration is below a threshold, what is the optimal lease duration which maximizes the net customer demand served by the wireless operators? Increasing or decreasing lease duration has many competing effects; while shorter lease duration may increase the efficiency of spectrum allocation, longer lease duration may increase market competition by incentivizing more operators to enter the market. We formulate this problem as a two-stage Stackelberg game consisting of the regulator and the wireless operators and design efficient algorithms to find the Stackelberg equilibrium of the entire game. These algorithms can also be used to find the Stackelberg equilibrium under some generalizations of our model. Using these algorithms, we obtain important numerical results and insights that characterize how the optimal lease duration varies with respect to market parameters in order to maximize the spectrum utilization. A few of our numerical results are non-intuitive as they suggest that increasing market competition may not necessarily improve spectrum utilization. To the best of our knowledge, this paper presents the first mathematical approach to optimize the lease duration of spectrum licenses.


title: Categorising Software Contexts: Research-in-Progress

abstract: A growing number of researchers suggest that software process must be tailored to a project's context to achieve maximal performance. Researchers have studied 'context' in an ad-hoc way, with focus on those contextual factors that appear to be of significance. The result is that we have no useful basis upon which to contrast and compare studies. We are currently researching a theoretical basis for software context for the purpose of tailoring and note that a deeper consideration of the meaning of the term 'context' is required before we can proceed. In this paper, we examine the term and present a model based on insights gained from our initial categorisation of contextual factors from the literature. We test our understanding by analysing a further six documents. Our contribution thus far is a model that we believe will support a theoretical operationalisation of software context for the purpose of process tailoring.


title: Simulations of Dynamical Gas-Dust Circumstellar Disks: Going Beyond the Epstein Regime

abstract: In circumstellar disks, the size of dust particles varies from submicron to several centimeters, while planetesimals have sizes of hundreds of kilometers. Therefore, various regimes for the aerodynamic drag between solid bodies and gas can be realized in these disks, depending on the grain sizes and velocities: Epstein, Stokes, and Newton, as well as transitional regimes between them. For small bodies moving in the Epstein regime, the time required to establish the constant relative velocity between the gas and bodies can be much less than the dynamical time scale for the problem - the time for the rotation of the disk about the central body. In addition, the dust may be concentrated in individual regions of the disk, making it necessary to take into account the transfer of momentum between the dust and gas. It is shown that, for a system of equations for gas and monodisperse dust, a semi-implicit first-order approximation scheme in time in which the interphase interaction is calculated implicitly, while other forces, such as the pressure gradient and gravity are calculated explicitly, is suitable for stiff problems with intense interphase interactions and for computations of the drag in non-linear regimes. The piece-wise drag coefficient widely used in astrophysical simulations has a discontinuity at some values of the Mach and Knudsen numbers that are realized in a circumstellar disk. A continuous drag coefficient is presented, which corresponds to experimental dependences obtained for various drag regimes.


title: Can Massive MIMO Support URLLC?

abstract: We investigate the feasibility of using Massive MIMO to support URLLC in both coherence interval based and 3GPP compliant pilot settings. We consider grant-free uplink transmission with MMSE receiver and adopt 3GPP channel models. In the coherence interval based pilot setting, by extensive system level simulations, we find that using a Massive MIMO base station with 128 antennas and MMSE receiver, URLLC requirements can be achieved in Urban Macro (UMa) Non-Line of Sight (NLoS) with orthogonal pilots and Neyman-Pearson detector. However, in the 3GPP compliant pilot setting, even by using the covariance matrix of Physical Resource Block (PRB) subcarriers for active UE detection and channel estimation as well as open-loop power control, we find that URLLC requirements are still challenging to achieve due to the insufficient pilot length and pilot symbol location regulations in a PRB.


title: Solving the linear transport equation by a deep neural network approach

abstract: In this paper, we study the linear transport model by adopting the deep learning method, in particular the deep neural network (DNN) approach. While the interest of using DNN to study partial differential equations is arising, here we adapt it to study kinetic models, in particular the linear transport model. Moreover, theoretical analysis on the convergence of the neural network and its approximated solution towards the analytic solution is shown. We demonstrate the accuracy and effectiveness of the proposed DNN method in the numerical experiments.


title: Determination of responses of liquid xenon to low energy electron and nuclear recoils using the PandaX-II detector

abstract: We report a systematic determination of the responses of PandaX-II, a dual phase xenon time projection chamber detector, to low energy recoils. The electron recoil (ER) and nuclear recoil (NR) responses are calibrated, respectively, with injected tritiated methane or $^{220}$Rn source, and with $^{241}$Am-Be neutron source, within an energy range from $1-25$ keV (ER) and $4-80$ keV (NR), under the two drift fields of 400 and 317 V/cm. An empirical model is used to fit the light yield and charge yield for both types of recoils. The best fit models can well describe the calibration data. The systematic uncertainties of the fitted models are obtained via statistical comparison against the data.


title: Robust and Differentially Private Mean Estimation

abstract: In statistical learning and analysis from shared data, which is increasingly widely adopted in platforms such as federated learning and meta-learning, there are two major concerns: privacy and robustness. Each participating individual should be able to contribute without the fear of leaking one's sensitive information. At the same time, the system should be robust in the presence of malicious participants inserting corrupted data. Recent algorithmic advances in learning from shared data focus on either one of these threats, leaving the system vulnerable to the other. We bridge this gap for the canonical problem of estimating the mean from i.i.d. samples. We introduce PRIME, which is the first efficient algorithm that achieves both privacy and robustness for a wide range of distributions. We further complement this result with a novel exponential time algorithm that improves the sample complexity of PRIME, achieving a near-optimal guarantee and matching a known lower bound for (non-robust) private mean estimation. This proves that there is no extra statistical cost to simultaneously guaranteeing privacy and robustness.


title: Gravitational perturbations in the Newman-Penrose formalism: Applications to wormholes

abstract: In this work we study the problem of linear stability of gravitational perturbations in stationary and spherically symmetric wormholes. For this purpose, we employ the Newman-Penrose formalism which is well-suited for treating gravitational radiation in General Relativity, as well as the geometrical aspect of this theory. With this method we obtain a "master equation" that describes the behavior of gravitational perturbations that are of odd-parity in the Regge-Wheeler gauge. This equation is later applied to a specific class of Morris-Thorne wormholes and also to the metric of an asymptotically flat scalar field wormhole. The analysis of the equations that these space-times yield reveals that there are no unstable vibrational modes generated by the type of perturbations here studied.


title: On the Sample Complexity of Stability Constrained Imitation Learning

abstract: We study the following question in the context of imitation learning for continuous control: how are the underlying stability properties of an expert policy reflected in the sample-complexity of an imitation learning task? We provide the first results showing that a surprisingly granular connection can be made between the underlying expert system's incremental gain stability, a novel measure of robust convergence between pairs of system trajectories, and the dependency on the task horizon $T$ of the resulting generalization bounds. In particular, we propose and analyze incremental gain stability constrained versions of behavior cloning and a DAgger-like algorithm, and show that the resulting sample-complexity bounds naturally reflect the underlying stability properties of the expert system. As a special case, we delineate a class of systems for which the number of trajectories needed to achieve $\varepsilon$-suboptimality is sublinear in the task horizon $T$, and do so without requiring (strong) convexity of the loss function in the policy parameters. Finally, we conduct numerical experiments demonstrating the validity of our insights on both a simple nonlinear system for which the underlying stability properties can be easily tuned, and on a high-dimensional quadrupedal robotic simulation.


title: Optimal Spectrum Partitioning and Licensing in Tiered Access under Stochastic Market Models

abstract: We consider the problem of partitioning a spectrum band into M channels of equal bandwidth, and then further assigning these M channels into P licensed channels and M-P unlicensed channels. Licensed channels can be accessed both for licensed and opportunistic use following a tiered structure which has a higher priority for licensed use. Unlicensed channels can be accessed only for opportunistic use. We address the following question in this paper. Given a market setup, what values of M and P maximize the net spectrum utilization of the spectrum band? While this problem is of fundamental nature, it is highly relevant practically, e.g., in the context of partitioning the recently proposed Citizens Broadband Radio Service band. If M is too high or too low, it may decrease spectrum utilization due to limited channel capacity or due to wastage of channel capacity, respectively. If P is too high (low), it will not incentivize the wireless operators who are primarily interested in unlicensed channels (licensed channels) to join the market. These tradeoffs are captured in our optimization problem which manifests itself as a two-stage Stackelberg game. We design an algorithm to solve the Stackelberg game and hence find the optimal M and P. The algorithm design also involves an efficient Monte Carlo integrator to evaluate the expected value of the involved random variables like spectrum utilization and operators' revenue. We also benchmark our algorithms using numerical simulations.


title: Light Curve Calculations for Triple Microlensing Systems

abstract: We present a method to compute the magnification of a finite source star lensed by a triple lens system based on the image boundary (contour integration) method. We describe a new procedure to obtain continuous image boundaries from solutions of the tenth-order polynomial obtained from the lens equation. Contour integration is then applied to calculate the image areas within the image boundaries, which yields the magnification of a source with uniform brightness. We extend the magnification calculation to limb-darkened stars approximated with a linear profile. In principle, this method works for all multiple lens systems, not just triple lenses. We also include an adaptive sampling and interpolation method for calculating densely covered light curves. The C++ source code and a corresponding Python interface are publicly available.


title: Symmetry enriched phases of quantum circuits

abstract: Quantum circuits consisting of random unitary gates and subject to local measurements have been shown to undergo a phase transition, tuned by the rate of measurement, from a state with volume-law entanglement to an area-law state. From a broader perspective, these circuits generate a novel ensemble of quantum many-body states at their output. In this paper, we characterize this ensemble and classify the phases that can be established as steady states. Symmetry plays a nonstandard role in that the physical symmetry imposed on the circuit elements does not on its own dictate the possible phases. Instead, it is extended by dynamical symmetries associated with this ensemble to form an enlarged symmetry. Thus, we predict phases that have no equilibrium counterpart and could not have been supported by the physical circuit symmetry alone. We give the following examples. First, we classify the phases of a circuit operating on qubit chains with $\mathbb{Z}_2$ symmetry. One striking prediction, corroborated with numerical simulation, is the existence of distinct volume-law phases in one dimension, which nonetheless support true long-range order. We furthermore argue that owing to the enlarged symmetry, this system can in principle support a topological area-law phase, protected by the combination of the circuit symmetry and a dynamical permutation symmetry. Second, we consider a Gaussian fermionic circuit that only conserves fermion parity. Here the enlarged symmetry gives rise to a $U(1)$ critical phase at moderate measurement rates and a Kosterlitz-Thouless transition to area-law phases. We comment on the interpretation of the different phases in terms of the capacity to encode quantum information. We discuss close analogies to the theory of spin glasses pioneered by Edwards and Anderson as well as crucial differences that stem from the quantum nature of the circuit ensemble.

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


title: On p-parts of character degrees

abstract: In this paper, we get the sharp bound for $|G/O_p(G)|_p$ under the assumption that either $p^2 \nmid \chi(1)$ for all $\chi \in {\rm Irr}(G)$ or $p^2 \nmid \phi(1)$ for all $\phi \in {\rm IBr}_p(G)$. This would settle two conjectures raised by Lewis, Navarro, Tiep, and Tong-Viet.


title: Facing to Latency of Hyperledger Fabric for Blockchain-enabled IoT: Modeling and Analysis

abstract: Hyperledger Fabric (HLF), one of the most popular private blockchains, has recently received attention for blockchain-enabled Internet of Things (IoT). However, for IoT applications to handle time-sensitive data, the processing latency in HLF has emerged as a new challenge. In this article, therefore, we establish a practical HLF latency model for HLF-enabled IoT. We first discuss the structure and transaction flow of HLF-enabled IoT. After implementing real HLF, we capture the latencies that each transaction experiences and show that the total latency of HLF can be modeled as a Gamma distribution, which is validated by conducting a goodness-of-fit test (i.e., the Kolmogorov-Smirnov (KS) test). We also provide the parameter values of the modeled latency distribution for various HLF environments. Furthermore, we explore the impacts of three important HLF parameters including transaction generation rate, block size, and block-generation timeout on the HLF latency. As a result, this article provides design insights on minimizing the latency for HLF-enabled IoT.


title: Characterization of probability measures based on Q-independent generalized random fields

abstract: Prakasa Rao (Studia Sci. Math. Hungar. 11 (1976) 277-282) studied a characterization of probability distributions for linear functions of independent generalized random fields. These results are extended to Q-independent generalized random fields. It is known that independence implies Q-independence but the converse is not true.


title: Gaussian Kernelized Self-Attention for Long Sequence Data and Its Application to CTC-based Speech Recognition

abstract: Self-attention (SA) based models have recently achieved significant performance improvements in hybrid and end-to-end automatic speech recognition (ASR) systems owing to their flexible context modeling capability. However, it is also known that the accuracy degrades when applying SA to long sequence data. This is mainly due to the length mismatch between the inference and training data because the training data are usually divided into short segments for efficient training. To mitigate this mismatch, we propose a new architecture, which is a variant of the Gaussian kernel, which itself is a shift-invariant kernel. First, we mathematically demonstrate that self-attention with shared weight parameters for queries and keys is equivalent to a normalized kernel function. By replacing this kernel function with the proposed Gaussian kernel, the architecture becomes completely shift-invariant with the relative position information embedded using a frame indexing technique. The proposed Gaussian kernelized SA was applied to connectionist temporal classification (CTC) based ASR. An experimental evaluation with the Corpus of Spontaneous Japanese (CSJ) and TEDLIUM 3 benchmarks shows that the proposed SA achieves a significant improvement in accuracy (e.g., from 24.0% WER to 6.0% in CSJ) in long sequence data without any windowing techniques.


title: Geometric Structure of Affine Deligne-Lusztig Varieties for $GL_3$

abstract: In this paper we study the geometric structure of affine Deligne-Lusztig varieties for $GL_3$ and $b$ basic. We completely determine the irreducible components of the affine Deligne-Lusztig variety. In particular, we classify the cases where all of the irreducible components are classical Deligne-Lusztig varieties times finite-dimensional affine spaces. If this is the case, then the irreducible components are pairwise disjoint.


title: Subgrid multiscale stabilized finite element analysis of non-Newtonian Power-law model fully coupled with Advection-Diffusion-Reaction equations

abstract: This article presents stability and convergence analyses of subgrid multiscale stabilized finite element formulation of non-Newtonian power-law fluid flow model strongly coupled with variable coefficients Advection-Diffusion-Reaction ($VADR$) equation. Considering the highly non-linear viscosity coefficient as solute concentration dependent makes the coupling two way. The stabilized formulation of the transient coupled system is developed based upon time dependent subscales, which ensures inherent consistency of the method. The proposed algebraic expressions of the stabilization parameters appropriately shape up the apriori and aposteriori error estimates. Both the shear thinning and shear thickening properties, indicated by different power-law indices are properly highlighted in theoretical derivations as well as in numerical validations. The numerical experiments carried out for different combinations of small and large Reynolds numbers and power-law indices establish far better performance of time dependent $ASGS$ method in approximating the solution of this coupled system for all the cases over the other well known stabilized finite element methods.


title: Data Poisoning Attacks and Defenses to Crowdsourcing Systems

abstract: A key challenge of big data analytics is how to collect a large volume of (labeled) data. Crowdsourcing aims to address this challenge via aggregating and estimating high-quality data (e.g., sentiment label for text) from pervasive clients/users. Existing studies on crowdsourcing focus on designing new methods to improve the aggregated data quality from unreliable/noisy clients. However, the security aspects of such crowdsourcing systems remain under-explored to date. We aim to bridge this gap in this work. Specifically, we show that crowdsourcing is vulnerable to data poisoning attacks, in which malicious clients provide carefully crafted data to corrupt the aggregated data. We formulate our proposed data poisoning attacks as an optimization problem that maximizes the error of the aggregated data. Our evaluation results on one synthetic and two real-world benchmark datasets demonstrate that the proposed attacks can substantially increase the estimation errors of the aggregated data. We also propose two defenses to reduce the impact of malicious clients. Our empirical results show that the proposed defenses can substantially reduce the estimation errors of the data poisoning attacks.


title: PIC Simulation of Double Plasma Resonance and Zebra Pattern of Solar Radio Bursts

abstract: Latest study reports that plasma emission can be generated by energetic electrons of DGH distribution via the electron cyclotron maser instability (ECMI) in plasmas characterized by a large ratio of plasma oscillation frequency to electron gyro-frequency ($\omega_{pe}/\Omega_{ce}$). In this study, on the basis of the ECMI-plasma emission mechanism, we examine the double plasma resonance (DPR) effect and the corresponding plasma emission at both harmonic (H) and fundamental (F) bands using PIC simulations with various $\omega_{pe}/\Omega_{ce}$. This allows us to directly simulate the feature of zebra pattern (ZP) observed in solar radio bursts for the first time. We find that (1) the simulations reproduce the DPR effect nicely for the upper hybrid (UH) and Z modes, as seen from their variation of intensity and linear growth rate with $\omega_{pe}/\Omega_{ce}$, (2) the intensity of the H emission is stronger than that of the F emission by $\sim$ 2 orders of magnitude and vary periodically with increasing $\omega_{pe}/\Omega_{ce}$, while the F emission is too weak to be significant, therefore we suggest that it is the H emission accounting for solar ZPs, (3) the peak-valley contrast of the total intensity of H is $\sim 4$, and the peak lies around integer values of $\omega_{pe}/\Omega_{ce}$ (= 10 and 11) for the present parameter setup. We also evaluate the effect of energy of energetic electrons on the characteristics of ECMI-excited waves and plasma radiation. The study provides novel insight on the physical origin of ZPs of solar radio bursts.


title: Deep Neural Networks based Invisible Steganography for Audio-into-Image Algorithm

abstract: In the last few years, steganography has attracted increasing attention from a large number of researchers since its applications are expanding further than just the field of information security. The most traditional method is based on digital signal processing, such as least significant bit encoding. Recently, there have been some new approaches employing deep learning to address the problem of steganography. However, most of the existing approaches are designed for image-in-image steganography. In this paper, the use of deep learning techniques to hide secret audio into the digital images is proposed. We employ a joint deep neural network architecture consisting of two sub-models: the first network hides the secret audio into an image, and the second one is responsible for decoding the image to obtain the original audio. Extensive experiments are conducted with a set of 24K images and the VIVOS Corpus audio dataset. Through experimental results, it can be seen that our method is more effective than traditional approaches. The integrity of both image and audio is well preserved, while the maximum length of the hidden audio is significantly improved.


title: Critical Effect of Bottom Electrode on Ferroelectricity of Epitaxial Hf0.5Zr0.5O2 Thin Films

abstract: Epitaxial orthorhombic Hf0.5Zr0.5O2 (HZO) films on La0.67Sr0.33MnO3 (LSMO) electrodes show robust ferroelectricity, with high polarization, endurance and retention. However, no similar results have been achieved using other perovskite electrodes so far. Here, LSMO and other perovskite electrodes are compared. A small amount of orthorhombic phase and low polarization is found in HZO films grown on La-doped BaSnO3 and Nb-doped SrTiO3, while null amounts of orthorhombic phase and polarization are detected in films on LaNiO3 and SrRuO3. The critical effect of the electrode on the stabilized phases is not consequence of differences in the electrode lattice parameter. The interface is critical, and engineering the HZO bottom interface on just a few monolayers of LSMO permits the stabilization of the orthorhombic phase. Furthermore, while the specific divalent ion (Sr or Ca) in the manganite is not relevant, reducing the La content causes a severe reduction of the amount of orthorhombic phase and the ferroelectric polarization in the HZO film.


title: Connection Problem for an Extension of $q$-Hypergeometric Systems

abstract: We give an example of solutions of the connection problem associated with a certain system of linear $q$-difference equations recently introduced by Park. The result contains a connection formulas of the $q$-Lauricella hypergeometric function $\varphi_{D}$ and those of the $q$-generalized hypergeometric function ${}_{N+1}\varphi_{N}$ as special cases.


title: Geometric and algebraic approaches to quantum theory

abstract: We show how to formulate physical theory taking as a starting point the set of states (geometric approach). We discuss the relation of this formulation to the conventional approach to classical and quantum mechanics and the theory of complex systems. The equations of motion and the formulas for probabilities of physical quantities are analyzed. A heuristic proof of decoherence in our setting is used to justify the formulas for probabilities. We show that any physical theory theory can be obtained from classical theory if we restrict the set of observables. This remark can be used to construct models with any prescribed group of symmetries; one can hope that this construction leads to new interesting models that cannot be build in the conventional framework. The geometric approach can be used to formulate quantum theory in terms of Jordan algebras, generalizing the algebraic approach to quantum theory. The scattering theory can be formulated in geometric approach.


title: QSO photometric redshifts using machine learning and neural networks

abstract: The scientific value of the next generation of large continuum surveys would be greatly increased if the redshifts of the newly detected sources could be rapidly and reliably estimated. Given the observational expense of obtaining spectroscopic redshifts for the large number of new detections expected, there has been substantial recent work on using machine learning techniques to obtain photometric redshifts. Here we compare the accuracy of the predicted photometric redshifts obtained from Deep Learning(DL) with the k-Nearest Neighbour (kNN) and the Decision Tree Regression (DTR) algorithms. We find using a combination of near-infrared, visible and ultraviolet magnitudes, trained upon a sample of SDSS QSOs, that the kNN and DL algorithms produce the best self-validation result with a standard deviation of {\sigma} = 0.24. Testing on various sub-samples, we find that the DL algorithm generally has lower values of {\sigma}, in addition to exhibiting a better performance in other measures. Our DL method, which uses an easy to implement off-the-shelf algorithm with no filtering nor removal of outliers, performs similarly to other, more complex, algorithms, resulting in an accuracy of {\Delta}z < 0.1$ up to z ~ 2.5. Applying the DL algorithm trained on our 70,000 strong sample to other independent (radio-selected) datasets, we find {\sigma} < 0.36 over a wide range of radio flux densities. This indicates much potential in using this method to determine photometric redshifts of quasars detected with the Square Kilometre Array.


title: Voltage control of skyrmions: creation, annihilation and zero-magnetic field stablization

abstract: Voltage manipulation of skyrmions is a promising path towards low-energy spintronic devices. Here, voltage effects on skyrmions in a GdOx/Gd/Co/Pt heterostructure are observed experimentally. The results show that the skyrmion density can be both enhanced and depleted by the application of an electric field, along with the ability, at certain magnetic fields to completely switch the skyrmion state on and off. Further, a zero magnetic field skyrmion state can be stablized under a negative bias voltage using a defined voltage and magnetic field sequence. The voltage effects measured here occur on a few-second timescale, suggesting an origin in voltage-controlled magnetic anisotropy rather than ionic effects. By investigating the skyrmion nucleation rate as a function of temperature, we extract the energy barrier to skyrmion nucleation in our sample. Further, micromagnetic simulations are used to explore the effect of changing the anisotropy and Dzyaloshinskii-Moriya interaction on skyrmion density. Our work demonstrates the control of skyrmions by voltages, showing functionalities desirable for commercial devices.


title: Optical design of a multi-object fiber-fed spectrograph system for Southern Spectroscopic Survey Telescope

abstract: Southern Spectroscopic Survey Telescope (SSST) is a wide-field spectroscopic survey telescope that China plans to build in Chile in the next few years. As an instrument for astronomical spectroscopic survey, the multi-object and fiber-fed spectrograph (MOFFS) is one of the most important scientific instruments for SSST. In this paper, we present a recommended optical design for the MOFFS system based on the Volume Phase Holographic Gratings (VPHG). The whole design philosophy and procedure, including the analytic method to determine the initial structure, optimization procedures of the VPHG and the camera groups, are demonstrated in detail. The numerical results of the final obtained spectrograph show a superior imaging quality and a relatively high transmittance for the whole working waveband and the field of view. The design method proposed in this paper can provide a reference for the design of MOFFS accommodated in spectroscopic survey telescopes.


title: A maximum entropy model of bounded rational decision-making with prior beliefs and market feedback

abstract: Bounded rationality is an important consideration stemming from the fact that agents often have limits on their processing abilities, making the assumption of perfect rationality inapplicable to many real tasks. We propose an information-theoretic approach to the inference of agent decisions under Smithian competition. The model explicitly captures the boundedness of agents (limited in their information-processing capacity) as the cost of information acquisition for expanding their prior beliefs. The expansion is measured as the Kullblack-Leibler divergence between posterior decisions and prior beliefs. When information acquisition is free, the homo economicus agent is recovered, while in cases when information acquisition becomes costly, agents instead revert to their prior beliefs. The maximum entropy principle is used to infer least-biased decisions based upon the notion of Smithian competition formalised within the Quantal Response Statistical Equilibrium framework. The incorporation of prior beliefs into such a framework allowed us to systematically explore the effects of prior beliefs on decision-making in the presence of market feedback, as well as importantly adding a temporal interpretation to the framework. We verified the proposed model using Australian housing market data, showing how the incorporation of prior knowledge alters the resulting agent decisions. Specifically, it allowed for the separation of past beliefs and utility maximisation behaviour of the agent as well as the analysis into the evolution of agent beliefs.

categories: cs.IT cs.AI econ.TH math.IT physics.soc-ph q-fin.CP


title: Information Carrier and Resource Optimization of Counterfactual Quantum Communication

abstract: Counterfactual quantum communication is unique in its own way that allows remote parties to transfer information without sending any message carrier in the channel. Although no message carrier travels in the channel at the time of successful information transmission, it is impossible to transmit information faster than the speed of light, thus without an information carrier. In this paper, we address an important question What carries the information in counterfactual quantum communication? and optimize the resource efficiency of the counterfactual quantum communication in terms of the number of channels used, time consumed to transmit 1-bit classical information between two remote parties, and the number of qubits required to accomplish the counterfactual quantum communication.


title: Testing Lotka's Law and Pattern of Author Productivity in the Scholarly Publications of Artificial Intelligence

abstract: Artificial intelligence has changed our day to day life in multitude ways. AI technology is rearing itself as a driving force to be reckoned with in the largest industries in the world. AI has already engulfed our educational system, our businesses and our financial establishments. The future is definite that machines with artificial intelligence will soon be captivating over trained manual work that now is mostly cared by humans. Machines can carry out human-like tasks by new inputs as artificial intelligence makes it possible for machines to learn from experience. AI data from web of science database from 2008 to 2017 have been mapped to depict the average growth rate, relative growth rate, contribution made by authors in the view of research productivity, authorship pattern and collaboration of AI literature. The Lotka's law on authorship productivity of AI literature has been tested to confirm the applicability of the law to the present data set. A K-S test was applied to measure the degree of agreement between the distribution of the observed set of data against the inverse general power relationship and the theoretical value of {\alpha} =2. It is found that the inverse square law of Lotka follow as such.


title: A LAMOST BHB Catalog and Kinematics Therein I: Catalog and Halo Properties

abstract: We collect a sample of stars observed both in LAMOST and Gaia which have colors implying a temperature hotter than 7000 K. We train a machine learning algorithm on LAMOST spectroscopic data which has been tagged with stellar classifications and metallicities, and use this machine to construct a catalog of Blue Horizontal Branch stars (BHBs) with metallicity information. Another machine is trained using Gaia parallaxes to predict absolute magnitudes for these stars. The final catalog of 13,693 BHBs is thought to be about 86\% pure, with $\sigma_{[Fe/H]}\sim$0.35 dex and $\sigma_{G}\sim$0.31 mag. These values are confirmed via comparison to globular clusters, although a covariance error seems to affect our magnitude and abundance estimates. We analyze a subset of this catalog in the Galactic Halo. We find that BHB populations in the outer halo appear redder, which could imply a younger population, and that the metallicity gradient is relatively flat around [Fe/H] = -1.9 dex over our sample footprint. We find that our metal rich BHB stars are on more radial velocity dispersion dominated orbits ($\beta \sim 0.70$) at all radii than our metal poor BHB stars ($\beta \sim 0.62$).


title: Formalization of Bohr's contextuality within theory of open quantum systems

abstract: In quantum physics, the notion of contextuality has a variety of interpretations which are typically associated with the names of their inventors, say Bohr, Bell, Kochen and Specker, and recently Dzhafarov. In fact, Bohr was the first who pointed to contextuality of quantum measurements as a part of formulation of his principle of complementarity. (Instead of "contextuality", he considered dependence on "experimental conditions.") Unfortunately, the contextuality counterpart of the complementarity principle was overshadowed by the issue of incompatibility of observables. And the interest for contextuality of quantum measurements rose again only in connection with the Bell inequality. The original Bohr's contextuality, as contextuality of each quantum measurement, was practically forgotten. It was highlighted in the works of the author of this paper, with applications both to physics and cognition. In this note, the theory of open quantum systems is applied to formalization of Bohr's contextuality within the the scheme of indirect measurements. This scheme is widely used in quantum information theory and it leads to the theory of quantum instruments (Davis-Lewis-Ozawa). In this scheme, Bohr's viewpoint on contextuality of quantum measurements is represented in the formal mathematical framework.


title: Link Prediction Approach to Recommender Systems

abstract: The problem of recommender system is very popular with myriad available solutions. A novel approach that uses the link prediction problem in social networks has been proposed in the literature that model the typical user-item information as a bipartite network in which link prediction would actually mean recommending an item to a user. The standard recommender system methods suffer from the problems of sparsity and scalability. Since link prediction measures involve computations pertaining to small neighborhoods in the network, this approach would lead to a scalable solution to recommendation. One of the issues in this conversion is that link prediction problem is modelled as a binary classification task whereas the problem of recommender systems is solved as a regression task in which the rating of the link is to be predicted. We overcome this issue by predicting top k links as recommendations with high ratings without predicting the actual rating. Our work extends similar approaches in the literature by focusing on exploiting the probabilistic measures for link prediction. Moreover, in the proposed approach, prediction measures that utilize temporal information available on the links prove to be more effective in improving the accuracy of prediction. This approach is evaluated on the benchmark 'Movielens' dataset. We show that the usage of temporal probabilistic measures helps in improving the quality of recommendations. Temporal random-walk based measure T_Flow improves recommendation accuracy by 4% and Temporal cooccurrence probability measure improves prediction accuracy by 10% over item-based collaborative filtering method in terms of AUROC score.


title: Hierarchical Attention Fusion for Geo-Localization

abstract: Geo-localization is a critical task in computer vision. In this work, we cast the geo-localization as a 2D image retrieval task. Current state-of-the-art methods for 2D geo-localization are not robust to locate a scene with drastic scale variations because they only exploit features from one semantic level for image representations. To address this limitation, we introduce a hierarchical attention fusion network using multi-scale features for geo-localization. We extract the hierarchical feature maps from a convolutional neural network (CNN) and organically fuse the extracted features for image representations. Our training is self-supervised using adaptive weights to control the attention of feature emphasis from each hierarchical level. Evaluation results on the image retrieval and the large-scale geo-localization benchmarks indicate that our method outperforms the existing state-of-the-art methods. Code is available here: \url{https://github.com/YanLiqi/HAF}.


title: What drives adsorption of ions on surface of nanodiamonds in aqueous solutions?

abstract: It is not yet clear what drives the adsorption of ions on detonation nanodiamonds (DNDs), which plays a critical role on the loading (unloading) of chemotherapeutic drugs on (from) the surface of DNDs in their targeted therapy applications. Furthermore, effects of adsorbed ions on the hydration layers of water around DNDs with different surface chemistries have not been studied yet. Through a series of Molecular Dynamics simulations, we found out that the law of matching water affinity generally explains well the adsorption patterns of ions onto the surface functional groups of DNDs. Depending on whether the water affinity of the ion matches with that of the surface functional group or not, the former predominantly forms either Contact Ion-Pair (CIP) or Solvent-shared Ion-Pair (SIP) with the latter. In this regard, Na$^{+}$ and Mg$^{2+}$ have the highest tendencies to form, respectively, CIP and SIP associations with ${-}$COO$^{-}$ functional groups. In the extreme case of 84 ${-}$COO$^{-}$ groups on DND${-}$COOH, however, we observed few Mg$^{2+}$${-}$COO$^{-}$ CIP associations, for which we have proposed a hypothesis based on the entropy gains. Furthermore, Mg$^{2+}$ and to a lesser extent Ca$^{2+}$ in cooperation with ${-}$COO$^{-}$ functional groups on the surface of charged DND${-}$COOH lead to relatively high residence times of water in the first hydration layer of DND. This study also provides a firsthand molecular level insight about the preferential orientation of water in the vicinity of positively charged DND${-}$H, on which prior experimental studies have not yet reached a consensus.


title: Edge Sparse Basis Network: A Deep Learning Framework for EEG Source Localization

abstract: EEG source localization is an important technical issue in EEG analysis. Despite many numerical methods existed for EEG source localization, they all rely on strong priors and the deep sources are intractable. Here we propose a deep learning framework using spatial basis function decomposition for EEG source localization. This framework combines the edge sparsity prior and Gaussian source basis, called Edge Sparse Basis Network (ESBN). The performance of ESBN is validated by both synthetic data and real EEG data during motor tasks. The results suggest that the supervised ESBN outperforms the traditional numerical methods in synthetic data and the unsupervised fine-tuning provides more focal and accurate localizations in real data. Our proposed deep learning framework can be extended to account for other source priors, and the real-time property of ESBN can facilitate the applications of EEG in brain-computer interfaces and clinics.


title: Coherently controlled quantum features in a coupled interferometric scheme

abstract: Over the last several decades, entangled photon pairs generated by spontaneous parametric down conversion processes in both second-order and third-order nonlinear optical materials have been intensively studied for various quantum features such as Bell inequality violation and anticorrelation. In an interferometric scheme, anticorrelation results from photon bunching based on randomness when entangled photon pairs coincidently impinge on a beam splitter. Compared with post-measurement-based probabilistic confirmation, a coherence version has been recently proposed using the wave nature of photons. Here, the origin of quantum features in a coupled interferometric scheme is investigated using pure coherence optics. In addition, a deterministic method of entangled photon-pair generation is proposed for on-demand coherence control of quantum processing.


title: The Diverse Morphology, Stellar Population, and Black Hole Scaling Relations of the Host Galaxies of Nearby Quasars

abstract: We present rest-frame $B$ and $I$ imaging of 35 low-redshift ($z < 0.5$) Palomar-Green quasars using the Hubble Space Telescope Wide Field Camera 3. We perform multi-component two-dimensional image decomposition to separate the host galaxy from its bright active nucleus, characterize its morphology, and measure its photometric properties. Special care is devoted to quantifying the structural parameters of the galaxy bulge, determine its $B-I$ color, and estimate its stellar mass. Roughly half of the sample, comprising the less luminous ($L_{5100} \lesssim 10^{45}\,\mathrm{erg\,s^{-1}}$) but most high Eddington ratio quasars, reside in disk galaxies that are often barred and possess pseudo bulges. The large stellar masses, large effective radii, and faint surface brightnesses suggest that the host galaxies of the most luminous quasars are mostly ellipticals. Major mergers constitute only a minority ($\lesssim 20\%$) of our sample. Our quasar sample roughly obeys the scaling relations between black hole mass and host galaxy (bulge, core, total) stellar mass. Hosts with black holes more massive than $\sim 10^8\,M_\odot$ behave similarly to classical bulges and early-type galaxies, while those with less massive black holes, particular the narrow-line Seyfert 1s, are consistent with pseudo bulges in late-type galaxies. The host galaxy bulges, irrespective of whether they are classical or pseudo, follow the relatively tight inverse relation between effective radius and mean effective surface brightness of inactive classical bulges and ellipticals. We argue that pseudo bulges experience recent or ongoing nuclear star formation.


title: Non-approximate Inference for Collective Graphical Models on Path Graphs via Discrete Difference of Convex Algorithm

abstract: The importance of aggregated count data, which is calculated from the data of multiple individuals, continues to increase. Collective Graphical Model (CGM) is a probabilistic approach to the analysis of aggregated data. One of the most important operations in CGM is maximum a posteriori (MAP) inference of unobserved variables under given observations. Because the MAP inference problem for general CGMs has been shown to be NP-hard, an approach that solves an approximate problem has been proposed. However, this approach has two major drawbacks. First, the quality of the solution deteriorates when the values in the count tables are small, because the approximation becomes inaccurate. Second, since continuous relaxation is applied, the integrality constraints of the output are violated. To resolve these problems, this paper proposes a new method for MAP inference for CGMs on path graphs. First we show that the MAP inference problem can be formulated as a (non-linear) minimum cost flow problem. Then, we apply Difference of Convex Algorithm (DCA), which is a general methodology to minimize a function represented as the sum of a convex function and a concave function. In our algorithm, important subroutines in DCA can be efficiently calculated by minimum convex cost flow algorithms. Experiments show that the proposed method outputs higher quality solutions than the conventional approach.


title: Modeling of Magneto-Conductivity of Bismuth Selenide -- A Topological Insulator

abstract: We report the magneto-conductivity analysis of Bi2Se3 single crystal at different temperatures in a magnetic field range of 14Tesla. The single crystals are grown by the self-flux method and characterized through X-ray diffraction, Scanning Electron Microscopy, and Raman Spectroscopy. The single crystals show magnetoresistance (MR) of around 380 percent at a magnetic field of 14T and a temperature of 5K. The Hikami Larkin Nagaoka (HLN) equation has been used to fit the magneto-conductivity (MC) data. However, the HLN fitted curve deviates at higher magnetic fields above 1 Tesla, suggesting that the role of surface driven conductivity suppresses with an increasing magnetic field. This article proposes a speculative model comprising of surface-driven HLN and added quantum diffusive and bulk carriers driven classical terms. The model successfully explains the MC of the Bi2Se3 single crystal at various temperatures (5 to 200K) and applied magnetic fields (up to 14Tesla).


title: SeaPearl: A Constraint Programming Solver guided by Reinforcement Learning

abstract: The design of efficient and generic algorithms for solving combinatorial optimization problems has been an active field of research for many years. Standard exact solving approaches are based on a clever and complete enumeration of the solution set. A critical and non-trivial design choice with such methods is the branching strategy, directing how the search is performed. The last decade has shown an increasing interest in the design of machine learning-based heuristics to solve combinatorial optimization problems. The goal is to leverage knowledge from historical data to solve similar new instances of a problem. Used alone, such heuristics are only able to provide approximate solutions efficiently, but cannot prove optimality nor bounds on their solution. Recent works have shown that reinforcement learning can be successfully used for driving the search phase of constraint programming (CP) solvers. However, it has also been shown that this hybridization is challenging to build, as standard CP frameworks do not natively include machine learning mechanisms, leading to some sources of inefficiencies. This paper presents the proof of concept for SeaPearl, a new CP solver implemented in Julia, that supports machine learning routines in order to learn branching decisions using reinforcement learning. Support for modeling the learning component is also provided. We illustrate the modeling and solution performance of this new solver on two problems. Although not yet competitive with industrial solvers, SeaPearl aims to provide a flexible and open-source framework in order to facilitate future research in the hybridization of constraint programming and machine learning.


title: Maximum weighted induced forests and trees: New formulations and a computational comparative review

abstract: Given a graph $G=(V,E)$ with a weight $w_v$ associated with each vertex $v\in V$, the maximum weighted induced forest problem (MWIF) consists of encountering a maximum weighted subset $V'\subseteq V$ of the vertices such that $V'$ induces a forest. This NP-hard problem is known to be equivalent to the minimum weighted feedback vertex set problem, which has applicability in a variety of domains. The closely related maximum weighted induced tree problem (MWIT), on the other hand, requires that the subset $V'\subseteq V$ induces a tree. We propose two new integer programming formulations with an exponential number of constraints and branch-and-cut procedures for MWIF. Computational experiments using benchmark instances are performed comparing several formulations, including the newly proposed approaches and those available in the literature, when solved by a standard commercial mixed integer programming solver. More specifically, five formulations are compared, two compact (i.e., with a polynomial number of variables and constraints) ones and the three others with an exponential number of constraints. The experiments show that a new formulation for the problem based on directed cutset inequalities for eliminating cycles (DCUT) offers stronger linear relaxation bounds earlier in the search process. The results also indicate that the other new formulation, denoted tree with cycle elimination (TCYC), outperforms those available in the literature when it comes to the average times for proving optimality for the small instances, especially the more challenging ones. Additionally, this formulation can achieve much lower average times for solving the larger random instances that can be optimally solved. Furthermore, we show how the formulations for MWIF can be easily extended for MWIT. Such extension allowed us to compare the optimal solution values of the two problems.


title: Strain localization and yielding dynamics in disordered collagen networks

abstract: Collagen is the most abundant extracellular-matrix protein in mammals and the main structural and load-bearing element of connective tissues. Collagen networks show remarkable strain-stiffening properties which tune the mechanical functions of tissues and regulate cell behaviours. Linear and non-linear mechanical properties of in-vitro disordered collagen networks have been widely studied using rheology for a range of self-assembly conditions in recent years. However, a one to one correlation between the onset of macroscopic network failure and local deformations inside the sample is yet to be established in these systems. Here, using shear rheology and in-situ high-resolution boundary imaging, we study the yielding dynamics of in-vitro reconstituted networks of uncrosslinked type-I collagen. We find that in the non-linear regime, the differential shear modulus ($K$) of the network initially increases with applied strain and then begins to drop as the network starts to yield beyond a critical strain (yield strain). Measurement of the local velocity profile using colloidal tracer particles reveals that beyond the peak of $K$, strong strain-localization and slippage between the network and the rheometer plate sets in that eventually leads to a detachment. We generalize this observation for a range of collagen concentrations, applied strain ramp rates, as well as different network architectures obtained by varying the polymerization temperature. Furthermore, by fitting the stress vs strain data with a continuum affine network model, we map out a state diagram showing the dependence of yield-stain and -stress on the reduced persistence length and mesh size of the network. Our findings can have broad implications in tissue engineering, particularly, in designing highly resilient biological scaffolds.


title: DeepMuD: Multi-user Detection for Uplink Grant-Free NOMA IoT Networks via Deep Learning

abstract: In this letter, we propose a deep learning-aided multi-user detection (DeepMuD) in uplink non-orthogonal multiple access (NOMA) to empower the massive machine-type communication where an offline-trained Long Short-Term Memory (LSTM)-based network is used for multi-user detection. In the proposed DeepMuD, a perfect channel state information (CSI) is also not required since it is able to perform a joint channel estimation and multi-user detection with the pilot responses, where the pilot-to-frame ratio is very low. The proposed DeepMuD improves the error performance of the uplink NOMA significantly and outperforms the conventional detectors (even with perfect CSI). Moreover, this gain becomes superb with the increase in the number of Internet of Things (IoT) devices. Furthermore, the proposed DeepMuD has a flexible detection and regardless of the number of IoT devices, the multi-user detection can be performed. Thus, an arbitrary number of IoT devices can be served without a signaling overhead, which enables the grant-free communication.


title: Network Size Estimation in Small-World Networks under Byzantine Faults

abstract: We study the fundamental problem of counting the number of nodes in a sparse network (of unknown size) under the presence of a large number of Byzantine nodes. We assume the full information model where the Byzantine nodes have complete knowledge about the entire state of the network at every round (including random choices made by all the nodes), have unbounded computational power, and can deviate arbitrarily from the protocol. Our main contribution is a randomized distributed algorithm that estimates the size of a network under the presence of a large number of Byzantine nodes. In particular, our algorithm estimates the size of a sparse, "small-world", expander network with up to $O(n^{1 - \delta})$ Byzantine nodes, where $n$ is the (unknown) network size and $\delta$ can be be any arbitrarily small (but fixed) positive constant. Our algorithm outputs a (fixed) constant factor estimate of $\log(n)$ with high probability; the correct estimate of the network size will be known to a large fraction ($(1 - \epsilon)$-fraction, for any fixed positive constant $\epsilon$) of the honest nodes. Our algorithm is fully distributed, lightweight, and simple to implement, runs in $O(\log^3{n})$ rounds, and requires nodes to send and receive messages of only small-sized messages per round; any node's local computation cost per round is also small.


title: Learning Continuous Exponential Families Beyond Gaussian

abstract: We address the problem of learning of continuous exponential family distributions with unbounded support. While a lot of progress has been made on learning of Gaussian graphical models, we still lack scalable algorithms for reconstructing general continuous exponential families modeling higher-order moments of the data beyond the mean and the covariance. Here, we introduce a computationally efficient method for learning continuous graphical models based on the Interaction Screening approach. Through a series of numerical experiments, we show that our estimator maintains similar requirements in terms of accuracy and sample complexity scalings compared to alternative approaches such as maximization of conditional likelihood, while considerably improving upon the algorithm's run-time.


title: Minimizing false negative rate in melanoma detection and providing insight into the causes of classification

abstract: Our goal is to bridge human and machine intelligence in melanoma detection. We develop a classification system exploiting a combination of visual pre-processing, deep learning, and ensembling for providing explanations to experts and to minimize false negative rate while maintaining high accuracy in melanoma detection. Source images are first automatically segmented using a U-net CNN. The result of the segmentation is then used to extract image sub-areas and specific parameters relevant in human evaluation, namely center, border, and asymmetry measures. These data are then processed by tailored neural networks which include structure searching algorithms. Partial results are then ensembled by a committee machine. Our evaluation on the largest skin lesion dataset which is publicly available today, ISIC-2019, shows improvement in all evaluated metrics over a baseline using the original images only. We also showed that indicative scores computed by the feature classifiers can provide useful insight into the various features on which the decision can be based.


title: Unsupervised Clustering of Time Series Signals using Neuromorphic Energy-Efficient Temporal Neural Networks

abstract: Unsupervised time series clustering is a challenging problem with diverse industrial applications such as anomaly detection, bio-wearables, etc. These applications typically involve small, low-power devices on the edge that collect and process real-time sensory signals. State-of-the-art time-series clustering methods perform some form of loss minimization that is extremely computationally intensive from the perspective of edge devices. In this work, we propose a neuromorphic approach to unsupervised time series clustering based on Temporal Neural Networks that is capable of ultra low-power, continuous online learning. We demonstrate its clustering performance on a subset of UCR Time Series Archive datasets. Our results show that the proposed approach either outperforms or performs similarly to most of the existing algorithms while being far more amenable for efficient hardware implementation. Our hardware assessment analysis shows that in 7 nm CMOS the proposed architecture, on average, consumes only about 0.005 mm^2 die area and 22 uW power and can process each signal with about 5 ns latency.


title: The classical $\beta$-ensembles with $\beta$ proportional to $1/N$: from loop equations to Dyson's disordered chain

abstract: In the classical $\beta$-ensembles of random matrix theory, setting $\beta = 2 \alpha/N$ and taking the $N \to \infty$ limit gives a statistical state depending on $\alpha$. Using the loop equations for the classical $\beta$-ensembles, we study the corresponding eigenvalue density, its moments, covariances of monomial linear statistics, and the moments of the leading $1/N$ correction to the density. From earlier literature the limiting eigenvalue density is known to be related to classical functions. Our study gives a unifying mechanism underlying this fact, identifying in particular the Gauss hypergeometric differential equation determining the Stieltjes transform of the limiting density in the Jacobi case. Our characterisation of the moments and covariances of monomial linear statistics is through recurrence relations. Also, we extend recent work which begins with the $\beta$-ensembles in the high temperature limit and constructs a family of tridiagonal matrices referred to as $\alpha$-ensembles, obtaining a random anti-symmetric tridiagonal matrix with i.i.d.~gamma distributed random variables. From this we are able to supplement analytic results obtained by Dyson in the study of the so-called type I disordered chain.

categories: math.PR cond-mat.stat-mech math-ph math.MP math.SP


title: Low Resource Audio-to-Lyrics Alignment From Polyphonic Music Recordings

abstract: Lyrics alignment in long music recordings can be memory exhaustive when performed in a single pass. In this study, we present a novel method that performs audio-to-lyrics alignment with a low memory consumption footprint regardless of the duration of the music recording. The proposed system first spots the anchoring words within the audio signal. With respect to these anchors, the recording is then segmented and a second-pass alignment is performed to obtain the word timings. We show that our audio-to-lyrics alignment system performs competitively with the state-of-the-art, while requiring much less computational resources. In addition, we utilise our lyrics alignment system to segment the music recordings into sentence-level chunks. Notably on the segmented recordings, we report the lyrics transcription scores on a number of benchmark test sets. Finally, our experiments highlight the importance of the source separation step for good performance on the transcription and alignment tasks. For reproducibility, we publicly share our code with the research community.


title: A new take on the inflationary quintessence

abstract: The quintessence field coupled to the cosmic neutrino background (CNB) has been widely discussed as an alternative mechanism to address the coincidence problem. As it is well known, it is possible to extend such models to obtain quintessential inflation, that is, to incorporate inflationary stage as well. Taking an alternative route, one can start from the well established inflationary models and obtain successful quintessence models at the expense of coupling with the CNB. To Follow this route, we use a slightly reformulated model addressed in PRD95, 123521 (2017). This particular model assumes $\mathcal{Z}_2$ symmetry for both scalar field potential and coupling term, which then breaks down in course of the cosmological evolution. For our discussion, however, the $\mathcal{Z}_2$ symmetry of the potential is not mandatory the model to work. The conventional mechanism of particle production by the oscillating inflaton field (and their subsequent thermalization) remains operative. It is plain to see that the proposed construction can be easily applied for many successful models of inflation to incorporate dark energy at the expense of coupling with the CNB. We address the issue of neutrino nuggets from the quantum field theory point of view. Namely, these nuggets are considered as bound states caused basically by the Yukawa force, which arises in the framework of linear perturbation theory due to exchange of virtual quanta of quintessence field between the neutrinos.


title: Towards a mathematical theory of trajectory inference

abstract: We devise a theoretical framework and a numerical method to infer trajectories of a stochastic process from samples of its temporal marginals. This problem arises in the analysis of single cell RNA-sequencing data, which provide high dimensional measurements of cell states but cannot track the trajectories of the cells over time. We prove that for a class of stochastic processes it is possible to recover the ground truth trajectories from limited samples of the temporal marginals at each time-point, and provide an efficient algorithm to do so in practice. The method we develop, Global Waddington-OT (gWOT), boils down to a smooth convex optimization problem posed globally over all time-points involving entropy-regularized optimal transport. We demonstrate that this problem can be solved efficiently in practice and yields good reconstructions, as we show on several synthetic and real datasets.

categories: stat.ML cs.LG math.OC math.PR math.ST stat.TH


title: Clustering by quantum annealing on three-level quantum elements qutrits

abstract: Clustering is grouping of data by the proximity of some properties. We report on the possibility of increasing the efficiency of clustering of points in a plane using artificial quantum neural networks after the replacement of the two-level neurons called qubits represented by the spins S = 1/2 by the three-level neurons called qutrits represented by the spins S = 1. The problem has been solved by the slow adiabatic change of the Hamiltonian in time. The methods for controlling a qutrit system using projection operators have been developed and the numerical simulation has been performed. The Hamiltonians for two well-known cluster-ing methods, one-hot encoding and k-means ++, have been built. The first method has been used to partition a set of six points into three or two clusters and the second method, to partition a set of nine points into three clusters and seven points into four clusters. The simulation has shown that the clustering problem can be ef-fectively solved on qutrits represented by the spins S = 1. The advantages of clustering on qutrits over that on qubits have been demonstrated. In particular, the number of qutrits required to represent data points is smaller than the number of qubits by a factor of log2 N / log3 N . Since, for qutrits, it is easier to partition the data points into three clusters rather than two ones, the approximate hierarchical procedure of data partition-ing into a larger number of clusters is accelerated. At the exact data partition into more than three clusters, it has been proposed to number the clusters by the numbers of states of the corresponding multi-spin subsys-tems, instead of using the numbers of individual spins. This reduces even more the number of qutrits (N log3 K instead of NK ) required to implement the algorithm.