title: On minima of sum of theta functions and Mueller-Ho Conjecture

abstract: Let $z=x+iy \in \mathbb{H}:=\{z= x+ i y\in\mathbb{C}: y>0\}$ and $ \theta (s;z)=\sum_{(m,n)\in\mathbb{Z}^2 } e^{-s \frac{\pi }{y }|mz+n|^2}$ be the theta function associated with the lattice $\Lambda ={\mathbb Z}\oplus z{\mathbb Z}$. In this paper we consider the following pair of minimization problems $$ \min_{ \mathbb{H} } \theta (2;\frac{z+1}{2})+\rho\theta (1;z),\;\;\rho\in[0,\infty),$$ $$ \min_{ \mathbb{H} } \theta (1; \frac{z+1}{2})+\rho\theta (2; z),\;\;\rho\in[0,\infty),$$ where the parameter $\rho\in[0,\infty)$ represents the competition of two intertwining lattices. We find that as $\rho$ varies the optimal lattices admit a novel pattern: they move from rectangular (the ratio of long and short side changes from $\sqrt3$ to 1), square, rhombus (the angle changes from $\pi/2$ to $\pi/3$) to hexagonal; furthermore, there exists a closed interval of $\rho$ such that the optimal lattices is always square lattice. This is in sharp contrast to optimal lattice shapes for single theta function ($\rho=\infty$ case), for which the hexagonal lattice prevails. As a consequence, we give a partial answer to optimal lattice arrangements of vortices in competing systems of Bose-Einstein condensates as conjectured (and numerically and experimentally verified) by Mueller-Ho \cite{Mue2002}.


title: Unmanned Aerial Systems for Wildland and Forest Fires

abstract: Wildfires represent an important natural risk causing economic losses, human death and important environmental damage. In recent years, we witness an increase in fire intensity and frequency. Research has been conducted towards the development of dedicated solutions for wildland and forest fire assistance and fighting. Systems were proposed for the remote detection and tracking of fires. These systems have shown improvements in the area of efficient data collection and fire characterization within small scale environments. However, wildfires cover large areas making some of the proposed ground-based systems unsuitable for optimal coverage. To tackle this limitation, Unmanned Aerial Systems (UAS) were proposed. UAS have proven to be useful due to their maneuverability, allowing for the implementation of remote sensing, allocation strategies and task planning. They can provide a low-cost alternative for the prevention, detection and real-time support of firefighting. In this paper we review previous work related to the use of UAS in wildfires. Onboard sensor instruments, fire perception algorithms and coordination strategies are considered. In addition, we present some of the recent frameworks proposing the use of both aerial vehicles and Unmanned Ground Vehicles (UV) for a more efficient wildland firefighting strategy at a larger scale.


title: Discrete Approximations and Optimal Control of Nonsmooth Perturbed Sweeping Processes

abstract: The main goal of this paper is developing the method of discrete approximations to derive necessary optimality conditions for a class of constrained sweeping processes with nonsmooth perturbations. Optimal control problems for sweeping processes have been recently recognized among the most interesting and challenging problems in modern control theory for discontinuous differential inclusions with irregular dynamics and implicit state constrained, while deriving necessary optimality conditions for their local minimizers have been significantly based on the smoothness of controlled dynamic perturbations. To overcome these difficulties, we use the method of discrete approximations and employ advanced tools of second-order variational analysis. This approach allows us to obtain new necessary optimality conditions for nonsmooth and nonconvex discrete-time problems of the sweeping type. Then we employ the obtained conditions and the strong convergence of discrete approximations to establish novel results for original nonsmooth sweeping control problems that include extended Euler-Lagrange and maximization conditions for local minimizers. Finally, we present applications of the obtained results to solving a controlled mobile robot model with a nonsmooth sweeping dynamics that is of some practical interest.


title: Enhancement of acoustic spin pumping by acoustic distributed Bragg reflector cavity

abstract: Surface acoustic waves (SAWs) in the GHz frequency range can inject spin currents dynamically into adjacent nonmagnetic layers via spin pumping effect associated with ferromagnetic resonance. Here, we demonstrate an enhancement of acoustic ferromagnetic resonance and spin current generation by a pair of SAW reflector gratings, which form an acoustic analogue of the distributed Bragg reflector cavity. In the experiment, we confirmed 2.04 $\pm$ 0.02 times larger SAW power absorption in a device with cavity than in case of no acoustic cavity. We confirmed up to 2.96 $\pm$ 0.02 times larger spin current generation by measuring electric voltages generated by the inverse Edelstein effect (IEE) at the interface between Cu and Bi$_2$O$_3$. The results suggest that acoustic cavities would be useful to enhance the conversion efficiency in SAW driven coupled magnon-phonon dynamics.


title: Synonymy = Translational Equivalence

abstract: Synonymy and translational equivalence are the relations of sameness of meaning within and across languages. As the principal relations in wordnets and multi-wordnets, they are vital to computational lexical semantics, yet the field suffers from the absence of a common formal framework to define their properties and relationship. This paper proposes a unifying treatment of these two relations, which is validated by experiments on existing resources. In our view, synonymy and translational equivalence are simply different types of semantic identity. The theory establishes a solid foundation for critically re-evaluating prior work in cross-lingual semantics, and facilitating the creation, verification, and amelioration of lexical resources.


title: An efficient algorithm for weakly compressible flows in spherical geometries

abstract: This study proposes an algorithm for modeling compressible flows in spherical shells in nearly incompressible and weakly compressible regimes based on an implicit direction splitting approach. The method retains theoretically expected convergence rates and remains stable for extremely small values of the characteristic Mach number. The staggered spatial discretization on the MAC stencil, commonly used in numerical methods for incompressible Navier-Stokes equations, was found to be convenient for the discretization of the compressible Navier-Stokes equations written in the non-conservative form in terms of the primitive variables. This approach helped to avoid the high-frequency oscillations without any artificial stabilization terms. Nonlinear Picard iterations with the splitting error reduction were also implemented to allow one to obtain a solution of the fully nonlinear system of equations. These results, alongside excellent parallel performance, prove the viability of the direction splitting approach in large-scale high-resolution high-performance simulations of atmospheric and oceanic flows.

categories: math.NA cs.NA physics.ao-ph physics.comp-ph


title: A Swarm of Simple Robots Constructing Planar Shapes

abstract: We present a new version of our previously proposed algorithm enabling a swarm of robots to construct a desired shape from objects in the plane. We also describe a hardware realization for this system which makes use of simple and readily sourced components. We refer to the task as planar construction which is the gathering of ambient objects into some desired shape. As an example application, a swarm of robots could use this algorithm to not only gather waste material into a pile, but shape that pile into a line for easy collection. The shape is specified by an image known as the scalar field. The scalar field serves an analogous role to the template pheromones that guide the construction of complex natural structures such as termite mounds. In addition to describing the algorithm and hardware platform, we develop some performance insights using a custom simulation environment and present experimental results on physical robots.


title: LNMap: Departures from Isomorphic Assumption in Bilingual Lexicon Induction Through Non-Linear Mapping in Latent Space

abstract: Most of the successful and predominant methods for bilingual lexicon induction (BLI) are mapping-based, where a linear mapping function is learned with the assumption that the word embedding spaces of different languages exhibit similar geometric structures (i.e., approximately isomorphic). However, several recent studies have criticized this simplified assumption showing that it does not hold in general even for closely related languages. In this work, we propose a novel semi-supervised method to learn cross-lingual word embeddings for BLI. Our model is independent of the isomorphic assumption and uses nonlinear mapping in the latent space of two independently trained auto-encoders. Through extensive experiments on fifteen (15) different language pairs (in both directions) comprising resource-rich and low-resource languages from two different datasets, we demonstrate that our method outperforms existing models by a good margin. Ablation studies show the importance of different model components and the necessity of non-linear mapping.


title: Customer Data Sharing Platform: A Blockchain-Based Shopping Cart

abstract: We propose a new free eCommerce platform with blockchains that allows customers to connect to the seller directly, share personal data without losing control and ownership of it and apply it to the domain of shopping cart. Our new platform provides a solution to four important problems: private payment, ensuring privacy and user control, and incentives for sharing. It allows the trade to be open, transparent with immutable transactions that can be used for settling any disputes. The paper presents a case study of applying the framework for a shopping cart as one of the enterprise nodes of MultiChain which provides trading in ethers controlled by smart contracts and also collects user profile data and allows them to receive rewards for sharing their data with other business enterprises. It tracks who shared what, with whom, when, by what means and for what purposes in a verifiable fashion. The user data from the repository is converted into an open data format and shared via stream in the blockchain so that other nodes can efficiently process and use the data. The smart contract verifies and executes the agreed terms of use of the data and transfers digital tokens as a reward to the customer. The smart contract imposes double deposit collateral to ensure that all participants act honestly.


title: Hierarchy-Based Algorithms for Minimizing Makespan under Precedence and Communication Constraints

abstract: We consider the classic problem of scheduling jobs with precedence constraints on a set of identical machines to minimize the makespan objective function. Understanding the exact approximability of the problem when the number of machines is a constant is a well-known question in scheduling theory. Indeed, an outstanding open problem from the classic book of Garey and Johnson asks whether this problem is NP-hard even in the case of 3 machines and unit-length jobs. In a recent breakthrough, Levey and Rothvoss gave a $(1+\epsilon)$-approximation algorithm, which runs in nearly quasi-polynomial time, for the case when job have unit lengths. However, a substantially more difficult case where jobs have arbitrary processing lengths has remained open. We make progress on this more general problem. We show that there exists a $(1+\epsilon)$-approximation algorithm (with similar running time as that of Levey and Rothvoss) for the non-migratory setting: when every job has to be scheduled entirely on a single machine, but within a machine the job need not be scheduled during consecutive time steps. Further, we also show that our algorithmic framework generalizes to another classic scenario where, along with the precedence constraints, the jobs also have communication delay constraints. Both of these fundamental problems are highly relevant to the practice of datacenter scheduling.


title: Permutation tests under a rotating sampling plan with clustered data

abstract: Consider a population consisting of clusters of sampling units, evolving temporally, spatially, or according to other dynamics. We wish to monitor the evolution of its means, medians, or other parameters. For administrative convenience and informativeness, clustered data are often collected via a rotating plan. Under rotating plans, the observations in the same clusters are correlated, and observations on the same unit collected on different occasions are also correlated. Ignoring this correlation structure may lead to invalid inference procedures. Accommodating cluster structure in parametric models is difficult or will have a high level of misspecification risk. In this paper, we explore exchangeability in clustered data collected via a rotating sampling plan to develop a permutation scheme for testing various hypotheses of interest. We also introduce a semiparametric density ratio model to facilitate the multiple population structure in rotating sampling plans. The combination ensures the validity of the inference methods while extracting maximum information from the sampling plan. A simulation study indicates that the proposed tests firmly control the type I error whether or not the data are clustered. The use of the density ratio model improves the power of the tests.


title: Vassiliev Invariants for Flows Via Chern-Simons Perturbation Theory

abstract: The perturbative expansion of Chern-Simons gauge theory leads to invariants of knots and links, the finite type invariants or Vassiliev invariants. It has been proven that at any order in perturbation theory the resulting expression is an invariant of that order. Bott-Taubes integrals on configuration spaces are introduced in the present context to write Feynman diagrams at a given order in perturbation theory in a geometrical and topological setting. One of the consequences of the configuration space formalism is that the resulting amplitudes are given in cohomological terms. This cohomological structure can be used to translate Bott-Taubes integrals into Chern-Simons perturbative amplitudes and vice versa. In this article this program is performed up to third order in the coupling constant. This expands some work previously worked out by Thurston. Finally we take advantage of these results to incorporate in the formalism a smooth and divergenceless vector field on the $3$-manifold. The Bott-Taubes integrals obtained are used for constructing higher-order asymptotic Vassiliev invariants extending the work of Komendarczyk and Voli\'c.


title: Exponential decay for damped Klein-Gordon equations on asymptotically cylindrical and conic manifolds

abstract: We study the decay of the global energy for the damped Klein-Gordon equation on non-compact manifolds with finitely many cylindrical and subconic ends up to bounded perturbation. We prove that under the Geometric Control Condition, the decay is exponential, and that under the weaker Network Control Condition, the decay is logarithmic, by developing the global Carleman estimate with multiple weights.


title: Radiation Campaign of HPK Prototype LGAD sensors for the High-Granularity Timing Detector (HGTD)

abstract: We report on the results of a radiation campaign with neutrons and protons of Low Gain Avalanche Detectors (LGAD) produced by Hamamatsu (HPK) as prototypes for the High-Granularity Timing Detector (HGTD) in ATLAS. Sensors with an active thickness of 50~$\mu$m were irradiated in steps of roughly 2$\times$ up to a fluence of $3\times10^{15}~\mathrm{n_{eq}cm^{-2}}$. As a function of the fluence, the collected charge and time resolution of the irradiated sensors will be reported for operation at $-30^{\circ}$.


title: Organic Narrative Charts

abstract: Storyline visualizations display the interactions of groups and entities and their development over time. Existing approaches have successfully adopted the general layout from hand-drawn illustrations to automatically create similar depictions. Ward Shelley is the author of several diagrammatic paintings that show the timeline of art-related subjects, such as Downtown Body, a history of art scenes. His drawings include many stylistic elements that are not covered by existing storyline visualizations, like links between entities, splits and merges of streams, and tags or labels to describe the individual elements. We present a visualization method that provides a visual mapping for the complex relationships in the data, creates a layout for their display, and adopts a similar styling of elements to imitate the artistic appeal of such illustrations. We compare our results to the original drawings and provide an open-source authoring tool prototype.


title: Empower Entity Set Expansion via Language Model Probing

abstract: Entity set expansion, aiming at expanding a small seed entity set with new entities belonging to the same semantic class, is a critical task that benefits many downstream NLP and IR applications, such as question answering, query understanding, and taxonomy construction. Existing set expansion methods bootstrap the seed entity set by adaptively selecting context features and extracting new entities. A key challenge for entity set expansion is to avoid selecting ambiguous context features which will shift the class semantics and lead to accumulative errors in later iterations. In this study, we propose a novel iterative set expansion framework that leverages automatically generated class names to address the semantic drift issue. In each iteration, we select one positive and several negative class names by probing a pre-trained language model, and further score each candidate entity based on selected class names. Experiments on two datasets show that our framework generates high-quality class names and outperforms previous state-of-the-art methods significantly.


title: Positron scattering and transport in liquid helium

abstract: In previous papers we have proposed a method for the ab initio calculation of fully differential cross-sections for electron scattering in liquids and applied it to liquid argon, xenon and krypton. In this paper, we extend the procedure to the consideration of positron scattering in liquid helium, which is complicated by the annihilation process as well as the fact that the electron definition for the region "owned" by a target atom used previously does not have a positron analogue. We explore several physically motivated definitions to obtain effective positron scattering in the dense fluid. We find that our calculations of a pure helium system cannot precisely match experimental measurements, however by including a small admixture (<0.1%) of an impurity, we can obtain reasonable agreement in the dense gas phase. In contrast, our calculations do not match well to the liquid phase measurements. This provides motivation to explore further multiple scattering effects in the theory.


title: Numerical study of space charge electric field inside Resistive Plate Chamber

abstract: Resistive plate chamber (RPC) is one of the state-of-the-art particle detection technology for the High Energy Physics (HEP) experiments. The basic operating mechanism of an RPC involves ionization of gas due to the passage of charged particles followed by electron transport , avalanche, and subsequent electromagnetic induction on readout strips due to the movement of the electrons and ions. Especially during streamer mode of operation, the electric field applied to the RPC can get significantly modified due to the presence of large number of electrons and ions. In this study, we have worked on dominant issues related to the estimation of electric field due to the space charge arising out of the presence of electrons, ions within an RPC. For this purpose we have considered two approaches: representation of the space charge cloud as (a) a collection of ring charges, and (b) as a collection of line charges. The results from these different methods have been compared with results available in the literature.


title: A GPS spoofing detection and classification correlator-based technique using the LASSO

abstract: This work proposes a global navigation satellite system (GNSS) spoofing detection and classification technique for single antenna receivers. We formulate an optimization problem at the baseband correlator domain by using the Least Absolute Shrinkage and Selection Operator (LASSO). We model correlator tap outputs of the received signal to form a dictionary of triangle-shaped functions and leverage sparse signal processing to choose a decomposition of shifted matching triangles from said dictionary. The optimal solution of this minimization problem discriminates the presence of a potential spoofing attack peak by observing a decomposition of two different code-phase values (authentic and spoofed) in a sparse vector output. We use a threshold to mitigate false alarms. Furthermore, we present a variation of the minimization problem by enhancing the dictionary to a higher-resolution of shifted triangles. The proposed technique can be implemented as an advanced fine-acquisition monitoring tool to aid in the tracking loops for spoofing mitigation. In our experiments, we are able to distinguish authentic and spoofer peaks from synthetic data simulations and from a real dataset, namely, the Texas Spoofing Test Battery (TEXBAT). The proposed method achieves 0.3% detection error rate (DER) for a spoofer attack in nominal signal-to-noise ratio (SNR) conditions for an authentic-over-spoofer power of 3 dB.


title: Unifying the dynamical effects of quantum and classical noises

abstract: We develop a new master equation as a unified description of the effects of both quantum noise (system-bath interaction) and classical noise on a system's dynamics, using a two-dimensional series expansion method. When quantum and classical noises are both present, their combined effect on a system's dynamics is not necessarily a simple sum of the two individual effects. Thus previous master equations for open systems and those for classical noise, even when jointly used, may not capture the full physics. Our formalism can determine whether there is interference between quantum and classical noises and will be able to capture and describe such interference if there is any (in a perturbative manner). We find that, interestingly, second-order interference between quantum and classical noises vanishes identically. This work thus also serves to justify simple additive treatments of quantum and classical noises, especially in the weak coupling regime. For a Zeeman-splitted atom in a stochastic magnetic field interacting with an optical cavity, we use the formalism developed herein to find the overall decoherence rate between the atom's energy levels.


title: Normal coordinates in a system of coupled oscillators and influence of the masses of the springs

abstract: Experimental analysis of the motion in a system of two coupled oscillators with arbitrary initial conditions was performed and the normal coordinates were obtained directly. The system consisted of two gliders moving on an air track, joined together by a spring and joined by two other springs to the fixed ends. From the positions of the center of mass and the relative distance, acquired by analysis of the digital video of the experiment, normal coordinates were obtained, and by a non linear fit the normal frequencies were also obtained. It is shown that although the masses of the springs are relatively small compared to that of the gliders, it is necessary to take them into consideration to improve the agreement with the experimental results. This experimental-theoretical proposal is targeted to an undergraduate laboratory.


title: An Auto-Encoder Strategy for Adaptive Image Segmentation

abstract: Deep neural networks are powerful tools for biomedical image segmentation. These models are often trained with heavy supervision, relying on pairs of images and corresponding voxel-level labels. However, obtaining segmentations of anatomical regions on a large number of cases can be prohibitively expensive. Thus there is a strong need for deep learning-based segmentation tools that do not require heavy supervision and can continuously adapt. In this paper, we propose a novel perspective of segmentation as a discrete representation learning problem, and present a variational autoencoder segmentation strategy that is flexible and adaptive. Our method, called Segmentation Auto-Encoder (SAE), leverages all available unlabeled scans and merely requires a segmentation prior, which can be a single unpaired segmentation image. In experiments, we apply SAE to brain MRI scans. Our results show that SAE can produce good quality segmentations, particularly when the prior is good. We demonstrate that a Markov Random Field prior can yield significantly better results than a spatially independent prior. Our code is freely available at https://github.com/evanmy/sae.


title: A fast methodology for large-scale focusing inversion of gravity and magnetic data using the structured model matrix and the $2D$ fast Fourier transform

abstract: Focusing inversion of potential field data for the recovery of sparse subsurface structures from surface measurement data on a uniform grid is discussed. For the uniform grid the model sensitivity matrices exhibit block Toeplitz Toeplitz block structure, by blocks for each depth layer of the subsurface. Then, through embedding in circulant matrices, all forward operations with the sensitivity matrix, or its transpose, are realized using the fast two dimensional Fourier transform. Simulations demonstrate that this fast inversion algorithm can be implemented on standard desktop computers with sufficient memory for storage of volumes up to size $n \approx 1M$. The linear systems of equations arising in the focusing inversion algorithm are solved using either Golub Kahan bidiagonalization or randomized singular value decomposition algorithms in which all matrix operations with the sensitivity matrix are implemented using the fast Fourier transform. These two algorithms are contrasted for efficiency for large-scale problems with respect to the sizes of the projected subspaces adopted for the solutions of the linear systems. The presented results confirm earlier studies that the randomized algorithms are to be preferred for the inversion of gravity data, and that it is sufficient to use projected spaces of size approximately $m/8$, for data sets of size $m$. In contrast, the Golub Kahan bidiagonalization leads to more efficient implementations for the inversion of magnetic data sets, and it is again sufficient to use projected spaces of size approximately $m/8$. Moreover, it is sufficient to use projected spaces of size $m/20$ when $m$ is large, $m \approx 50000$, to reconstruct volumes with $n \approx 1M$. Simulations support the presented conclusions and are verified on the inversion of a practical magnetic data set that is obtained over the Wuskwatim Lake region in Manitoba, Canada.


title: End-to-end NILM System Using High Frequency Data and Neural Networks

abstract: Improving energy efficiency is a necessity in the fight against climate change. Non Intrusive Load Monitoring (NILM) systems give important information about the household consumption that can be used by the electric utility or the end users. In this work the implementation of an end-to-end NILM system is presented, which comprises a custom high frequency meter and neural-network based algorithms. The present article presents a novel way to include high frequency information as input of neural network models by means of multivariate time series that include carefully selected features. Furthermore, it provides a detailed assessment of the generalization error and shows that this class of models generalize well to new instances of seen-in-training appliances. An evaluation database formed of measurements in two Uruguayan homes is collected and discussion on general unsupervised approaches is provided.


title: Laurent skew orthogonal polynomials and related symplectic matrices

abstract: Particular class of skew orthogonal polynomials are introduced and investigated, which possess Laurent symmetry. They are also shown to appear as eigenfunctions of symplectic generalized eigenvalue problems. The modification of these polynomials gives some symplectic eigenvalue problem and the corresponding matrix is shown to be equivalent to butterfly matrix, which is a canonical form of symplectic matrices.


title: Synergistic CPU-FPGA Acceleration of Sparse Linear Algebra

abstract: This paper describes REAP, a software-hardware approach that enables high performance sparse linear algebra computations on a cooperative CPU-FPGA platform. REAP carefully separates the task of organizing the matrix elements from the computation phase. It uses the CPU to provide a first-pass re-organization of the matrix elements, allowing the FPGA to focus on the computation. We introduce a new intermediate representation that allows the CPU to communicate the sparse data and the scheduling decisions to the FPGA. The computation is optimized on the FPGA for effective resource utilization with pipelining. REAP improves the performance of Sparse General Matrix Multiplication (SpGEMM) and Sparse Cholesky Factorization by 3.2X and 1.85X compared to widely used sparse libraries for them on the CPU, respectively.


title: Interactive Rainbow Score: A Visual-centered Multimodal Flute Tutoring System

abstract: Learning to play an instrument is intrinsically multimodal, and we have seen a trend of applying visual and haptic feedback in music games and computer-aided music tutoring systems. However, most current systems are still designed to master individual pieces of music; it is unclear how well the learned skills can be generalized to new pieces. We aim to explore this question. In this study, we contribute Interactive Rainbow Score, an interactive visual system to boost the learning of sight-playing, the general musical skill to read music and map the visual representations to performance motions. The key design of Interactive Rainbow Score is to associate pitches (and the corresponding motions) with colored notation and further strengthen such association via real-time interactions. Quantitative results show that the interactive feature on average increases the learning efficiency by 31.1%. Further analysis indicates that it is critical to apply the interaction in the early period of learning.


title: Improving Vertical Positioning Accuracy with the Weighted Multinomial Logistic Regression Classifier

abstract: In this paper, a method of improving vertical positioning accuracy with the Global Positioning System (GPS) information and barometric pressure values is proposed. Firstly, we clear null values for the raw data collected in various environments, and use the 3$\sigma$-rule to identify outliers. Secondly, the Weighted Multinomial Logistic Regression (WMLR) classifier is trained to obtain the predicted altitude of outliers. Finally, in order to verify its effect, we compare the MLR method, the WMLR method, and the Support Vector Machine (SVM) method for the cleaned dataset which is regarded as the test baseline. The numerical results show that the vertical positioning accuracy is improved from 5.9 meters (the MLR method), 5.4 meters (the SVM method) to 5 meters (the WMLR method) for 67% test points.


title: Measurements of Absolute Branching Fractions of Fourteen Exclusive Hadronic $D$ Decays to $\eta$

abstract: Using $2.93\,\rm fb^{-1}$ of $e^+e^-$ collision data taken at a center-of-mass energy of 3.773\,GeV with the BESIII detector, we report the first measurements of the absolute branching fractions of fourteen hadronic $D^{0(+)}$ decays to exclusive final states with an $\eta$, e.g., $D^0\to K^-\pi^+\eta$, $K^0_S\pi^0\eta$, $K^+K^-\eta$, $K^0_SK^0_S\eta$, $K^-\pi^+\pi^0\eta$, $K^0_S\pi^+\pi^-\eta$, $K^0_S\pi^0\pi^0\eta$, and $\pi^+\pi^-\pi^0\eta$; $D^+\to K^0_S\pi^+\eta$, $K^0_SK^+\eta$, $K^-\pi^+\pi^+\eta$, $K^0_S\pi^+\pi^0\eta$, $\pi^+\pi^+\pi^-\eta$, and $\pi^+\pi^0\pi^0\eta$. Among these decays, the $D^0\to K^-\pi^+\eta$ and $D^+\to K^0_S\pi^+\eta$ decays have the largest branching fractions, which are $\mathcal{B} (D^0\to K^-\pi^+\eta )=(1.853\pm0.025_{\rm stat}\pm0.031_{\rm syst})\%$ and $\mathcal{B}(D^+\to K^0_S\pi^+\eta)=(1.309\pm0.037_{\rm stat}\pm0.031_{\rm syst})\%$, respectively. We also determine the $CP$ asymmetries for the six decays with highest event yields. No evidence of $CP$ violation is found.


title: Dynamic evolution of major element chemistry in protoplanetary disks and its implications for chondrite formation

abstract: Chondrites are the likely building blocks of Earth, and identifying the group of chondrite that best represents Earth is a key to resolving the state of the early Earth. The origin of chondrites, however, remains controversial partly because of their puzzling major element compositions, some exhibiting depletion in Al, Ca, and Mg. Based on a new thermochemical evolution model of protoplanetary disks, we show that planetesimals with depletion patterns similar to ordinary and enstatite chondrites can originate at 1-2 AU just outside where enstatite evaporates. Around the "evaporation front" of enstatite, the large inward flow of refractory minerals, including forsterite, takes place with a high pebble concentration, and the loss of those minerals result in depletion in Al, Ca, and Mg. When evaporated solid grains re-condense onto pebbles, the concentration of pebbles is further enhanced, potentially triggering the streaming instability. Planetesimals with the composition of ordinary and enstatite chondrites can thus be naturally created in the terrestrial region. The preferential loss of forsterite also creates an enhancement of Mg/Si and heavy Si isotopes just inside the potential source region for ordinary and enstatite chondrites. Earth, which shows both features, may originate just inside where ordinary and enstatite chondrites were born.


title: Neural Additive Models: Interpretable Machine Learning with Neural Nets

abstract: Deep neural networks (DNNs) are powerful black-box predictors that have achieved impressive performance on a wide variety of tasks. However, their accuracy comes at the cost of intelligibility: it is usually unclear how they make their decisions. This hinders their applicability to high stakes decision-making domains such as healthcare. We propose Neural Additive Models (NAMs) which combine some of the expressivity of DNNs with the inherent intelligibility of generalized additive models. NAMs learn a linear combination of neural networks that each attend to a single input feature. These networks are trained jointly and can learn arbitrarily complex relationships between their input feature and the output. Our experiments on regression and classification datasets show that NAMs are more accurate than widely used intelligible models such as logistic regression and shallow decision trees. They perform similarly to existing state-of-the-art generalized additive models in accuracy, but are more flexible because they are based on neural nets instead of boosted trees. To demonstrate this, we show how NAMs can be used for multitask learning on synthetic data and on the COMPAS recidivism data due to their composability, and demonstrate that the differentiability of NAMs allows them to train more complex interpretable models for COVID-19.


title: Machine Learning the Cosmic Curvature in a Model-independent Way

abstract: In this work, we achieve the determination of the cosmic curvature $\Omega_K$ in a cosmological model-independent way, by using the Hubble parameter measurements $H(z)$ and type Ia supernovae (SNe Ia). In our analysis, two nonlinear interpolating tools are used to reconstruct the Hubble parameter, one is the Artificial Neural Network (ANN) method, and the other is the Gaussian process (GP) method. We find that $\Omega_K$ based on the GP method can be greatly influenced by the prior of $H_0$, while the ANN method can overcome this. Therefore, the ANN method may have more advantages than GP in the measurement of the cosmic curvature. Based on the ANN method, we find a spatially open universe is preferred by the current $H(z)$ and SNe Ia data, and the difference between our result and the value inferred from Planck CMB is $1.6\sigma$. In order to test the reliability of the ANN method, and the potentiality of the future gravitational waves (GW) standard sirens in the measurement of the cosmic curvature, we constrain $\Omega_K$ using the simulated Hubble parameter and GW standard sirens in a model-independent way. We find that the ANN method is reliable and unbiased, and the error of $\Omega_K$ is $\sim0.186$ when 100 GW events with electromagnetic counterparts are detected, which is $\sim56\%$ smaller than that constrained from the Pantheon SNe Ia. Therefore, the data-driven method based on ANN has potential in the measurement of the cosmic curvature.


title: A generalized information criterion for high-dimensional PCA rank selection

abstract: Principal component analysis (PCA) is the most commonly used statistical procedure for dimension reduction. An important issue for applying PCA is to determine the rank, which is the number of dominant eigenvalues of the covariance matrix. The Akaike information criterion (AIC) and Bayesian information criterion (BIC) are among the most widely used rank selection methods. Both use the number of free parameters for assessing model complexity. In this work, we adopt the generalized information criterion (GIC) to propose a new method for PCA rank selection under the high-dimensional framework. The GIC model complexity takes into account the sizes of covariance eigenvalues and can be better adaptive to practical applications. Asymptotic properties of GIC are derived and the selection consistency is established under the generalized spiked covariance model.


title: HST survey of the Orion Nebula Cluster in the H$_2$O 1.4 $\mu$m absorption band: I. A census of substellar and planetary mass objects

abstract: In order to obtain a complete census of the stellar and sub-stellar population, down to a few M$_{Jup}$ in the $\sim1$ Myr old Orion Nebula Cluster, we used the infrared channel of the Wide Field Camera 3 of the Hubble Space Telescope with the F139M and F130N filters. These bandpasses correspond to the $1.4 \mu$m H$_2$O absorption feature and an adjacent line-free continuum region. Out of $4,504$ detected sources, $3,352$ (about $75\%$) appear fainter than m$_{130}=14$ (Vega mag) in the F130N filter, a brightness corresponding to the hydrogen-burning limit mass (M$\simeq 0.072 M_\odot$) at $\sim 1$ Myr. Of these, however, only $742$ sources have a negative F130M-139N color index, indicative of the presence of H$_2$O vapor in absorption, and can therefore be classified as bona-fide M and L dwarfs, with effective temperatures T$\lesssim 2850$ K at an assumed $1$ Myr cluster age. On our color-magnitude diagram, this population of sources with H$_2$O absorption appears clearly distinct from the larger background population of highly reddened stars and galaxies with positive F130M-F139N color index, and can be traced down to the sensitivity limit of our survey, m$_{130}\simeq 21.5$, corresponding to a $1$ Myr old $\simeq 3 $M$_{Jup}$, planetary mass object under about 2 magnitudes of visual extinction. Theoretical models of the BT-Settl family predicting substellar isochrones of $1, 2$ and $3$ Myr (down to $\sim 1 $M$_{Jup}$) fail to reproduce the observed H$_2$O color index at M$\lesssim 20 $M$_{Jup}$. We perform a Bayesian analysis to determine extinction, mass and effective temperature of each sub-stellar member of our sample, together with its membership probability.


title: On $q$-Isomonodromic Deformations and $q$-Nekrasov Functions

abstract: We construct a fundamental system of a $q$-difference Lax pair of rank $N$ in terms of 5d Nekrasov functions with $q=t$. Our fundamental system degenerates by the limit $q\to 1$ to a fundamental system of a differential Lax pair, which yields the Fuji-Suzuki-Tsuda system. We introduce tau functions of our system as Fourier transforms of 5d Nekrasov functions. Using asymptotic expansions of the fundamental system at $0$ and $\infty$, we obtain several determinantal identities of the tau functions.


title: Effects of Dark Energy anisotropic stress on the matter power spectrum

abstract: We study the effects of dark energy (DE) anisotropic stress on features of the matter power spectrum (PS). We employ the Parametrized Post-Friedmannian (PPF) formalism to emulate an effective DE, and model its anisotropic stress properties through a two-parameter equation that governs its overall amplitude ($g_0$) and transition scale ($c_g$). For the background cosmology, we consider different equations of state to model DE including a constant $w_0$ parameter, and models that provide thawing (CPL) and freezing (nCPL) behaviors. We first constrain these parameters by using the Pantheon, BAO, $H_0$ and CMB Planck data. Then, we analyze the role played by these parameters in the linear PS. In order for the anisotropic stress not to provoke deviations larger than $10\%$ and $5\%$ with respect to the $\Lambda$CDM PS at $k \sim 0.01 \,h/\text{Mpc}$, the parameters have to be in the range $-0.30< g_0 < 0.32$, $0 \leq c_g^2 < 0.01$ and $-0.15 < g_0 < 0.16$, $0 \leq c_g^2 < 0.01$, respectively. Additionally, we compute the leading nonlinear corrections to the PS using standard perturbation theory in real and redshift space, showing that the differences with respect to the $\Lambda$CDM are enhanced, especially for the quadrupole and hexadecapole RSD multipoles.


title: EmbraceNet for Activity: A Deep Multimodal Fusion Architecture for Activity Recognition

abstract: Human activity recognition using multiple sensors is a challenging but promising task in recent decades. In this paper, we propose a deep multimodal fusion model for activity recognition based on the recently proposed feature fusion architecture named EmbraceNet. Our model processes each sensor data independently, combines the features with the EmbraceNet architecture, and post-processes the fused feature to predict the activity. In addition, we propose additional processes to boost the performance of our model. We submit the results obtained from our proposed model to the SHL recognition challenge with the team name "Yonsei-MCML."


title: Technological improvement rate estimates for all technologies: Use of patent data and an extended domain description

abstract: In this work, we attempt to provide a comprehensive granular account of the pace of technological change. More specifically, we survey estimated yearly performance improvement rates for nearly all definable technologies for the first time. We do this by creating a correspondence of all patents within the US patent system to a set of technology domains. A technology domain is a body of patented inventions achieving the same technological function using the same knowledge and scientific principles. We obtain a set of 1757 domains using an extension of the previously defined classification overlap method (COM). These domains contain 97.14% of all patents within the entire US patent system. From the identified patent sets, we calculated the average centrality of the patents in each domain to estimate their improvement rates, following a methodology tested in prior work. The estimated improvement rates vary from a low of 1.9% per year for the Mechanical Skin treatment - Hair Removal and wrinkles domain to a high of 228.8% per year for the Network management - client-server applications domain. We developed a one-line descriptor identifying the technological function achieved and the underlying knowledge base for the largest 50, fastest 20 as well as slowest 20 of these domains, which cover more than forty percent of the patent system. In general, the rates of improvement were not a strong function of the patent set size and the fastest improving domains are predominantly software-based. We make available an online system that allows for automated searching for domains and improvement rates corresponding to any technology of interest to researchers, strategists and policy formulators.

categories: econ.GN cs.CY physics.soc-ph q-fin.EC q-fin.GN q-fin.PM


title: HST survey of the Orion Nebula Cluster in the H$_2$O 1.4 $\mu$m absorption band: II. The substellar IMF down to planetary masses

abstract: We exploit the ability of the Hubble Space Telescope to probe near infrared water absorption present in the atmosphere of low-mass stars, brown dwarfs and planetary mass objects to create a very pure sample of Orion Nebula Cluster (ONC) members, not affected by contamination from background stars and galaxies which lack water absorption. Thanks to these data we infer the Initial Mass Function (IMF) of the ONC in the $0.005 - 1.4$M$_{\odot}$ regime, i.e. down to few Jupiter masses. The young age of the ONC, $\sim1$ Myr, provides a snapshot of the outcome of star formation for the present-day conditions (metallicity, temperature, pressure) of typical Milky Way disk molecular clouds. We demonstrate that the IMF of the ONC is well described by either a log-normal function or a broken power-law, with parameter values qualitatively in agreement with the canonical Chabrier or Kroupa forms for the Milky Way disk IMF. This continuity in the mass distribution provides clues to the fact that the same physical processes may be regulating formation of stars, brown dwarfs, and planetary mass objects. Both the canonical IMF forms under-predict the observed number of very low mass members (below $0.1$ M$_\odot$), a regime where our data allows more precise constraints. Nevertheless, we do not observe a rise or secondary peak in the brown dwarfs or planetary mass regimes. Our study thus contradicts findings based on broad-band near infrared ground-based photometry, which predict an extremely high number of free-floating planets, but likely suffer from unaccounted background contamination.


title: Provably-secure symmetric private information retrieval with quantum cryptography

abstract: Private information retrieval (PIR) is a database query protocol that provides user privacy, in that the user can learn a particular entry of the database of his interest but his query would be hidden from the data centre. Symmetric private information retrieval (SPIR) takes PIR further by additionally offering database privacy, where the user cannot learn any additional entries of the database. Unconditionally secure SPIR solutions with multiple databases are known classically, but are unrealistic because they require long shared secret keys between the parties for secure communication and shared randomness in the protocol. Here, we propose using quantum key distribution (QKD) instead for a practical implementation, which can realise both the secure communication and shared randomness requirements. We prove that QKD maintains the security of the SPIR protocol and that it is also secure against any external eavesdropper. We also show how such a classical-quantum system could be implemented practically, using the example of a two-database SPIR protocol with keys generated by measurement device-independent QKD. Through key rate calculations, we show that such an implementation is feasible at the metropolitan level with current QKD technology.


title: Revisiting Pre-Trained Models for Chinese Natural Language Processing

abstract: Bidirectional Encoder Representations from Transformers (BERT) has shown marvelous improvements across various NLP tasks, and consecutive variants have been proposed to further improve the performance of the pre-trained language models. In this paper, we target on revisiting Chinese pre-trained language models to examine their effectiveness in a non-English language and release the Chinese pre-trained language model series to the community. We also propose a simple but effective model called MacBERT, which improves upon RoBERTa in several ways, especially the masking strategy that adopts MLM as correction (Mac). We carried out extensive experiments on eight Chinese NLP tasks to revisit the existing pre-trained language models as well as the proposed MacBERT. Experimental results show that MacBERT could achieve state-of-the-art performances on many NLP tasks, and we also ablate details with several findings that may help future research. Resources available: https://github.com/ymcui/MacBERT


title: Overcoming optical contrast limit of mesoscale light focusing by means of retrograde-reflection photonic nanojet

abstract: The physical origin of subwavelength photonic nanojet in retrograde-reflection mode (retro-PNJ) is theoretically considered. This specific type of photonic nanojet emerges upon sequential double focusing of a plane optical wave by a transparent dielectric microparticle located near a flat reflecting mirror. For the first time to the best of our knowledge, we report a unique property of retro-PNJ for increasing its focal length and intensity using microparticles with the optical contrast exceeding limiting value (n > 2) that fundamentally distinguishes retro-PNJ from regular PNJ behavior. This may drastically increase the trapping potential of PNJ-based optical tweezers.


title: HST survey of the Orion Nebula Cluster in the H$_2$O 1.4 $\mu$m absorption band: III. The population of sub-stellar binary companions

abstract: We present new results concerning the sub-stellar binary population in the Orion Nebula Cluster (ONC). Using the Karhunen-Lo\`{e}ve Image Projection (KLIP) algorithm, we have reprocessed images taken with the IR channel of the Wide Field Camera 3 mounted on the Hubble Space Telescope to unveil faint close companions in the wings of the stellar PSFs. Starting with a sample of $1392$ bona-fide not saturated cluster members, we detect $39$ close-pairs cluster candidates with separation $0.16''-0.77''$. The primary masses span a range M$_p$ $\sim 0.015-1.27$ M$_{\odot}$ whereas for the companions we derive M$_c$ $\sim 0.004-0.54$ M$_{\odot}$. Of these $39$ binary systems, $18$ were already known while the remaining $21$ are new detections. Correcting for completeness and combining our catalog with previously detected ONC binaries, we obtain an overall binary fraction of $11.5\% \pm 0.9\%$. Compared to other star forming regions, our multiplicity function is $\sim 2$ smaller than e.g. Taurus, while compared to the binaries in the field we obtain comparable values. We analyze the mass function of the binaries, finding differences between the mass distribution of binaries and single stars and between primary and companion mass distributions. The mass ratio shows a bottom-heavy distribution with median value of $M_c/M_p \sim 0.25$. Overall our results suggest that ONC binaries may represent a template for the typical population of field binaries, supporting the hypothesis that the ONC may be regarded as a most typical star forming region in the Milky Way.


title: Hybrid Adaptive Evolutionary Algorithm for Multi-objective Optimization

abstract: The major difficulty in Multi-objective Optimization Evolutionary Algorithms (MOEAs) is how to find an appropriate solution that is able to converge towards the true Pareto Front with high diversity. Most existing methodologies, which have demonstrated their niche on various practical problems involving two and three objectives, face significant challenges in the dependency of the selection of the EA parameters. Moreover, the process of setting such parameters is considered time-consuming, and several research works have tried to deal with this problem. This paper proposed a new Multi-objective Algorithm as an extension of the Hybrid Adaptive Evolutionary algorithm (HAEA) called MoHAEA. MoHAEA allows dynamic adaptation of the application of operator probabilities (rates) to evolve with the solution of the multi-objective problems combining the dominance- and decomposition-based approaches. MoHAEA is compared with four states of the art MOEAs, namely MOEA/D, pa$\lambda$-MOEA/D, MOEA/D-AWA, and NSGA-II on ten widely used multi-objective test problems. Experimental results indicate that MoHAEA outperforms the benchmark algorithms in terms of how it is able to find a well-covered and well-distributed set of points on the Pareto Front.


title: Intelligent networking with Mobile Edge Computing: Vision and Challenges for Dynamic Network Scheduling

abstract: Mobile edge computing (MEC) has been considered as a promising technique for internet of things (IoT). By deploying edge servers at the proximity of devices, it is expected to provide services and process data at a relatively low delay by intelligent networking. However, the vast edge servers may face great challenges in terms of cooperation and resource allocation. Furthermore, intelligent networking requires online implementation in distributed mode. In such kinds of systems, the network scheduling can not follow any previously known rule due to complicated application environment. Then statistical learning rises up as a promising technique for network scheduling, where edges dynamically learn environmental elements with cooperations. It is expected such learning based methods may relieve deficiency of model limitations, which enhance their practical use in dynamic network scheduling. In this paper, we investigate the vision and challenges of the intelligent IoT networking with mobile edge computing. From the systematic viewpoint, some major research opportunities are enumerated with respect to statistical learning.


title: Dynamic Anomaly Detection with High-fidelity Simulators: A Convex Optimization Approach

abstract: The main objective of this article is to develop scalable dynamic anomaly detectors when high-fidelity simulators of power systems are at our disposal. On the one hand, mathematical models of these high-fidelity simulators are typically "intractable" to apply existing model-based approaches. On the other hand, pure data-driven methods developed primarily in the machine learning literature neglect our knowledge about the underlying dynamics of the systems. In this study, we combine tools from these two mainstream approaches to develop a diagnosis filter that utilizes the knowledge of both the dynamical system as well as the simulation data of the high-fidelity simulators. The proposed diagnosis filter aims to achieve two desired features: (i) performance robustness with respect to model mismatch; (ii) high scalability. To this end, we propose a tractable (convex) optimization-based reformulation in which decisions are the filter parameters, the model-based information introduces feasible sets, and the data from the simulator forms the objective function to-be-minimized regarding the effect of model mismatch on the filter performance. To validate the theoretical results, we implement the developed diagnosis filter in DIgSILENT PowerFactory to detect false data injection attacks on the Automatic Generation Control measurements in the three-area IEEE 39-bus system.


title: Detection of multipolar orders in the spin-orbit-entangled 5d Mott insulator Ba2MgReO6

abstract: In electronic solids with strong spin-orbit interactions (SOIs), the spin and orbital degrees of freedom of an electron are quantum mechanically entangled, which may result in an exotic multipolar order instead of a conventional dipolar order such as a magnetic order. Such a higher-degree order is called "hidden order" because of difficulties in experimental detection. Moreover, the number of candidate compounds is limited, especially rare in d electron systems, in which an interplay between SOIs and Coulomb interactions is expected to cause rich physics. Here, we employ state-of-the-art synchrotron X-ray diffraction techniques on a high-quality single crystal to probe subtle symmetry breaking induced by a multipolar order. We unequivocally demonstrate that the double-perovskite Ba2MgReO6 exhibits successive transitions to quadrupolar and then dipolar orders upon cooling, which is consistent with a theory considering SOIs. Our findings are a significant step towards understanding the intriguing physics of multipoles realized by spin-orbit-entangled 5d electrons.


title: Hierarchies of holonomy groupoids for foliated bundles

abstract: We give a new construction of the holonomy groupoid of a regular foliation in terms of a partial connection on a diffeological principal bundle of germs of transverse parametrisations. We extend these ideas to construct a novel holonomy groupoid for any foliated bundle, which we prove sits at the top of a hierarchy of diffeological jet holonomy groupoids associated to the foliated bundle.


title: Task-Feature Collaborative Learning with Application to Personalized Attribute Prediction

abstract: As an effective learning paradigm against insufficient training samples, Multi-Task Learning (MTL) encourages knowledge sharing across multiple related tasks so as to improve the overall performance. In MTL, a major challenge springs from the phenomenon that sharing the knowledge with dissimilar and hard tasks, known as negative transfer, often results in a worsened performance. Though a substantial amount of studies have been carried out against the negative transfer, most of the existing methods only model the transfer relationship as task correlations, with the transfer across features and tasks left unconsidered. Different from the existing methods, our goal is to alleviate negative transfer collaboratively across features and tasks. To this end, we propose a novel multi-task learning method called Task-Feature Collaborative Learning (TFCL). Specifically, we first propose a base model with a heterogeneous block-diagonal structure regularizer to leverage the collaborative grouping of features and tasks and suppressing inter-group knowledge sharing. We then propose an optimization method for the model. Extensive theoretical analysis shows that our proposed method has the following benefits: (a) it enjoys the global convergence property and (b) it provides a block-diagonal structure recovery guarantee. As a practical extension, we extend the base model by allowing overlapping features and differentiating the hard tasks. We further apply it to the personalized attribute prediction problem with fine-grained modeling of user behaviors. Finally, experimental results on both simulated dataset and real-world datasets demonstrate the effectiveness of our proposed method


title: Span-based Localizing Network for Natural Language Video Localization

abstract: Given an untrimmed video and a text query, natural language video localization (NLVL) is to locate a matching span from the video that semantically corresponds to the query. Existing solutions formulate NLVL either as a ranking task and apply multimodal matching architecture, or as a regression task to directly regress the target video span. In this work, we address NLVL task with a span-based QA approach by treating the input video as text passage. We propose a video span localizing network (VSLNet), on top of the standard span-based QA framework, to address NLVL. The proposed VSLNet tackles the differences between NLVL and span-based QA through a simple yet effective query-guided highlighting (QGH) strategy. The QGH guides VSLNet to search for matching video span within a highlighted region. Through extensive experiments on three benchmark datasets, we show that the proposed VSLNet outperforms the state-of-the-art methods; and adopting span-based QA framework is a promising direction to solve NLVL.


title: CoronaVis: A Real-time COVID-19 Tweets Data Analyzer and Data Repository

abstract: Due to the nature of the data and public interaction, twitter is becoming more and more useful to understand and model various events. The goal of CoronaVis is to use tweets as the information shared by the people to visualize topic modeling, study subjectivity, and to model the human emotions during the COVID-19 pandemic. The main objective is to explore the psychology and behavior of the societies at large which can assist in managing the economic and social crisis during the ongoing pandemic as well as the after-effects of it. The novel coronavirus (COVID-19) pandemic forced people to stay at home to reduce the spread of the virus by maintaining social distancing. However, social media is keeping people connected both locally and globally. People are sharing information (e.g. personal opinions, some facts, news, status, etc.) on social media platforms which can be helpful to understand the various public behavior such as emotions, sentiments, and mobility during the ongoing pandemic. In this work, we develop a live application to observe the tweets on COVID-19 generated from the USA. In this paper, we have generated various data analytics over a period of time to study the changes in topics, subjectivity, and human emotions. We also share a cleaned and processed dataset named CoronaVis Twitter dataset (focused on the United States) available to the research community at https://github.com/mykabir/COVID19. This will enable the community to find more useful insights and create different applications and models to fight with COVID-19 pandemic and future pandemics as well.


title: On the complexity of Winner Verification and Candidate Winner for Multiwinner Voting Rules

abstract: The Chamberlin-Courant and Monroe rules are fundamental and well-studied rules in the literature of multi-winner elections. The problem of determining if there exists a committee of size k that has a Chamberlin-Courant (respectively, Monroe) score of at most r is known to be NP-complete. We consider the following natural problems in this setting: a) given a committee S of size k as input, is it an optimal k-sized committee, and b) given a candidate c and a committee size k, does there exist an optimal k-sized committee that contains c? In this work, we resolve the complexity of both problems for the Chamberlin-Courant and Monroe voting rules in the settings of rankings as well as approval ballots. We show that verifying if a given committee is optimal is coNP-complete whilst the latter problem is complete for $\Theta_{2}^{P}$. We also demonstrate efficient algorithms for the second problem when the input consists of single-peaked rankings. Our contribution fills an essential gap in the literature for these important multi-winner rules.


title: Pomeron-LQCD model of $J/\Psi$ photo-production on the nucleon

abstract: Based on the vector meson dominance assumption, a Hamiltonian model has been developed to investigate $J/\Psi$ photo-production reaction on the nucleon by using the $J/\Psi$-nucleon potential extracted from a lattice QCD calculation of Phys. Rev. D{\bf 82}, 091501 (2010). It is found that the predicted total cross sections are comparable to the recent data of $J/\Psi$ photo-production reaction from Jefferson Laboratory. The model is then extended to include the two-gluon exchange amplitude modeled by Donnachie and Lanshoff within Regge Phenomenology. The resulting Pomeron-LQCD model can then explain the data up to invariant mass $W=$ 300 GeV. Future improvements needed to reduce the uncertainties of the predictions are discussed. The need of an accurate extraction of $J/\Psi$-N potential at short distances from LQCD is illustrated.


title: An Averaging Processes on Hypergraphs

abstract: Consider the following iterated process on a hypergraph $H$. Each vertex $v$ has an initial vertex weight. At each step, we uniformly at random select an edge $F$ in $H$, and for each vertex $v$ in $F$ we replace the weight of $v$ by the average value of the vertex weights over all vertices in $F$. This is a generalization of an interactive process on graphs, first proposed by Aldous and Lanoue. In this paper, we use the eigenvalues of a Laplacian for hypergraphs to bound the rate of convergence for the iterated averaging process.


title: Multi-layer local optima networks for the analysis of advanced local search-based algorithms

abstract: A Local Optima Network (LON) is a graph model that compresses the fitness landscape of a particular combinatorial optimization problem based on a specific neighborhood operator and a local search algorithm. Determining which and how landscape features affect the effectiveness of search algorithms is relevant for both predicting their performance and improving the design process. This paper proposes the concept of multi-layer LONs as well as a methodology to explore these models aiming at extracting metrics for fitness landscape analysis. Constructing such models, extracting and analyzing their metrics are the preliminary steps into the direction of extending the study on single neighborhood operator heuristics to more sophisticated ones that use multiple operators. Therefore, in the present paper we investigate a twolayer LON obtained from instances of a combinatorial problem using bitflip and swap operators. First, we enumerate instances of NK-landscape model and use the hill climbing heuristic to build the corresponding LONs. Then, using LON metrics, we analyze how efficiently the search might be when combining both strategies. The experiments show promising results and demonstrate the ability of multi-layer LONs to provide useful information that could be used for in metaheuristics based on multiple operators such as Variable Neighborhood Search.


title: Revisiting Round-Trip Translation for Quality Estimation

abstract: Quality estimation (QE) is the task of automatically evaluating the quality of translations without human-translated references. Calculating BLEU between the input sentence and round-trip translation (RTT) was once considered as a metric for QE, however, it was found to be a poor predictor of translation quality. Recently, various pre-trained language models have made breakthroughs in NLP tasks by providing semantically meaningful word and sentence embeddings. In this paper, we employ semantic embeddings to RTT-based QE. Our method achieves the highest correlations with human judgments, compared to previous WMT 2019 quality estimation metric task submissions. While backward translation models can be a drawback when using RTT, we observe that with semantic-level metrics, RTT-based QE is robust to the choice of the backward translation system. Additionally, the proposed method shows consistent performance for both SMT and NMT forward translation systems, implying the method does not penalize a certain type of model.


title: Families of sequences with good family complexity and cross-correlation measure

abstract: In this paper we study pseudorandomness of a family of sequences in terms of two measures, the family complexity ($f$-complexity) and the cross-correlation measure of order $\ell$. We consider sequences not only on binary alphabet but also on $k$-symbols ($k$-ary) alphabet. We first generalize some known methods on construction of the family of binary pseudorandom sequences. We prove a bound on the $f$-complexity of a large family of binary sequences of Legendre-symbols of certain irreducible polynomials. We show that this family as well as its dual family have both a large family complexity and a small cross-correlation measure up to a rather large order. Next, we present another family of binary sequences having high $f$-complexity and low cross-correlation measure. Then we extend the results to the family of sequences on $k$-symbols alphabet.


title: Evaluating Transformer-Based Multilingual Text Classification

abstract: As NLP tools become ubiquitous in today's technological landscape, they are increasingly applied to languages with a variety of typological structures. However, NLP research does not focus primarily on typological differences in its analysis of state-of-the-art language models. As a result, NLP tools perform unequally across languages with different syntactic and morphological structures. Through a detailed discussion of word order typology, morphological typology, and comparative linguistics, we identify which variables most affect language modeling efficacy; in addition, we calculate word order and morphological similarity indices to aid our empirical study. We then use this background to support our analysis of an experiment we conduct using multi-class text classification on eight languages and eight models.


title: DS-FACTO: Doubly Separable Factorization Machines

abstract: Factorization Machines (FM) are powerful class of models that incorporate higher-order interaction among features to add more expressive power to linear models. They have been used successfully in several real-world tasks such as click-prediction, ranking and recommender systems. Despite using a low-rank representation for the pairwise features, the memory overheads of using factorization machines on large-scale real-world datasets can be prohibitively high. For instance on the criteo tera dataset, assuming a modest $128$ dimensional latent representation and $10^{9}$ features, the memory requirement for the model is in the order of $1$ TB. In addition, the data itself occupies $2.1$ TB. Traditional algorithms for FM which work on a single-machine are not equipped to handle this scale and therefore, using a distributed algorithm to parallelize the computation across a cluster is inevitable. In this work, we propose a hybrid-parallel stochastic optimization algorithm DS-FACTO, which partitions both the data as well as parameters of the factorization machine simultaneously. Our solution is fully de-centralized and does not require the use of any parameter servers. We present empirical results to analyze the convergence behavior, predictive power and scalability of DS-FACTO.


title: Highly tunable polarization-engineered two-dimensional electron gas in ${\epsilon}$-AlGaO3 / ${\epsilon}$-Ga2O3 heterostructures

abstract: We report on the modeling of polarization-induced two-dimensional electron gas (2DEG) formation at ${\epsilon}$-AlGaO3 / ${\epsilon}$-Ga2O3 heterointerface and the effect of spontaneous polarization (Psp) reversal on 2DEG density in ${\epsilon}$-Ga2O3 /${\epsilon}$-AlGaO3 / ${\epsilon}$-Ga2O3 double heterostructures. Density-functional theory (DFT) is utilized to calculate the material properties of ${\epsilon}$-Ga2O3 and ${\epsilon}$-AlGaO3 alloys. Using Schrodinger-Poisson solver along with DFT calculated parameters, the 2DEG density is calculated as a function of barrier type and thickness. By optimizing the layer thicknesses of ${\epsilon}$-Ga2O3/${\epsilon}$-AlGaO3/${\epsilon}$-Ga2O3 heterostructures, charge contrast ratios exceeding 1600 are obtained. This computational study indicates the high potential for ${\epsilon}$-Ga2O3-based heterostructure devices for non-volatile memories and neuromorphic applications.


title: On the design of non-autonomous fixed-time controllers with a predefined upper bound of the settling time

abstract: This paper aims to introduce a design methodology to stabilize a chain of integrators in a fixed-time with predefined Upper Bound for the Settling-Time (UBST). This approach is based on time-varying gains (time-base generator) that become singular as the time approaches the predefined convergence time. We present the conditions such that every nonzero trajectory convergence exactly at the predefined time with feedback laws that are simpler than in previous approaches. Additionally, we present results such that the origin is reached before the singularity occurs, making this approach realizable in practice. A significant contribution, since current results, based on time-varying gains, required the time-varying gain to tend to infinity as the time approaches the prescribed/predefined-time. Moreover, our approach guarantees fixed-time convergence with predefined UBST even in the presence of bounded disturbances (with known bound) and provides a methodology to re-design an autonomous fixed-time controller where the estimate of the UBST for the closed-loop system is very conservative into a less conservative one with predefined \textit{UBST} that can be set arbitrarily tight. We present numerical examples to illustrate the effectiveness of the approach together with comparisons with existing autonomous fixed-time controllers.


title: Wavelength-scale Optical Parametric Oscillators

abstract: Despite recent progress in nonlinear optics in wavelength-scale resonators, there are still open questions on the possibility of parametric oscillation in such resonators. We present a general approach to predict the behavior and estimate the oscillation threshold of multi-mode subwavelength and wavelength-scale optical parametric oscillators (OPOs). As an example, we propose an OPO based on Mie-type multipolar resonances, and we demonstrate that due to the low-Q nature of multipolar modes in wavelength-scale resonators, there is a nonlinear interaction between these modes. As a result, the OPO threshold, compared to the single-mode case, can be reduced by a factor that is significantly larger than the number of interacting modes. The multi-mode interaction can also lead to a phase transition manifested through a sudden change in the parametric gain as well as the oscillation threshold, which can be utilized for enhanced sensing. We establish an explicit connection between the second-harmonic generation efficiency and the OPO threshold. This allows us to estimate the OPO threshold based on measured or simulated second-harmonic generation in different classes of resonators, such as bound states in the continuum and inversely designed resonators. Our approach for analyzing and modeling miniaturized OPOs can open unprecedented opportunities for classical and quantum nonlinear photonics.


title: Rheology of a concentrated suspension of spherical squirmers: monolayer in simple shear flow

abstract: A concentrated, vertical monolayer of identical spherical squirmers, which may be bottom-heavy, and which are subjected to a linear shear flow, is modelled computationally by two different methods: Stokesian dynamics, and a lubrication-theory-based method. Inertia is negligible. The aim is to compute the effective shear viscosity and, where possible, the normal stress differences as functions of the areal fraction of spheres $\phi$, the squirming parameter $\beta$ (proportional to the ratio of a squirmer's active stresslet to its swimming speed), the ratio $Sq$ of swimming speed to a typical speed of the shear flow, the bottom-heaviness parameter $G_{bh}$, the angle $\alpha$ that the shear flow makes with the horizontal, and two parameters that define the repulsive force that is required computationally to prevent the squirmers from overlapping when their distance apart is less than a critical value $\epsilon a$, where $\epsilon$ is very small and $a$ is the sphere radius. The Stokesian dynamics method allows the rheological quantities to be computed for values of $\phi$ up to $0.75$; the lubrication-theory method can be used for $\phi> 0.5$. A major finding of this work is that, despite very different assumptions, the two methods of computation give overlapping results for viscosity as a function of $\phi$ in the range $0.5 < \phi < 0.75$. This suggests that lubrication theory, based on near-field interactions alone, contains most of the relevant physics, and that taking account of interactions with more distant particles than the nearest is not essential to describe the dominant physics.


title: Linguistic Resources for Bhojpuri, Magahi and Maithili: Statistics about them, their Similarity Estimates, and Baselines for Three Applications

abstract: Corpus preparation for low-resource languages and for development of human language technology to analyze or computationally process them is a laborious task, primarily due to the unavailability of expert linguists who are native speakers of these languages and also due to the time and resources required. Bhojpuri, Magahi, and Maithili, languages of the Purvanchal region of India (in the north-eastern parts), are low-resource languages belonging to the Indo-Aryan (or Indic) family. They are closely related to Hindi, which is a relatively high-resource language, which is why we compare with Hindi. We collected corpora for these three languages from various sources and cleaned them to the extent possible, without changing the data in them. The text belongs to different domains and genres. We calculated some basic statistical measures for these corpora at character, word, syllable, and morpheme levels. These corpora were also annotated with parts-of-speech (POS) and chunk tags. The basic statistical measures were both absolute and relative and were exptected to indicate of linguistic properties such as morphological, lexical, phonological, and syntactic complexities (or richness). The results were compared with a standard Hindi corpus. For most of the measures, we tried to the corpus size the same across the languages to avoid the effect of corpus size, but in some cases it turned out that using the full corpus was better, even if sizes were very different. Although the results are not very clear, we try to draw some conclusions about the languages and the corpora. For POS tagging and chunking, the BIS tagset was used to manually annotate the data. The POS tagged data sizes are 16067, 14669 and 12310 sentences, respectively, for Bhojpuri, Magahi and Maithili. The sizes for chunking are 9695 and 1954 sentences for Bhojpuri and Maithili, respectively.


title: Probing the initial conditions of high-mass star formation -- IV. Gas dynamics and NH$_2$D chemistry in high-mass precluster and protocluster clumps

abstract: The initial stage of star formation is a complex area study because of its high density and low temperature. Under such conditions, many molecules become depleted from the gas phase by freezing out onto dust grains. However, the deuterated species could remain gaseous and are thus ideal tracers. We investigate the gas dynamics and NH$_2$D chemistry in eight massive pre/protocluster clumps. We present NH$_2$D 1$_{11}$-1$_{01}$ (at 85.926 GHz), NH$_3$ (1, 1) and (2, 2) observations in the eight clumps using the PdBI and the VLA, respectively. We find that the distribution between deuterium fractionation and kinetic temperature shows a number density peak at around $T_{\rm kin}=16.1$ K, and the NH$_2$D cores are mainly located at a temperature range of 13.0 to 22.0 K. We detect seven instances of extremely high deuterium fractionation of $1.0 \leqslant D_{\rm frac} \leqslant 1.41$. We find that the NH$_2$D emission does not appear to coincide exactly with either dust continuum or NH$_3$ peak positions, but often surrounds the star-formation active regions. This suggests that the NH$_{2}$D has been destroyed by the central young stellar object (YSO) due to its heating. The detected NH$_2$D lines are very narrow with a median width of $\rm 0.98\pm0.02 km/s$. The extracted NH$_2$D cores are gravitationally bound ($\alpha_{\rm vir} < 1$), are likely prestellar or starless, and can potentially form intermediate-mass or high-mass stars. Using NH$_3$ (1, 1) as a dynamical tracer, we find very complicated dynamical movement, which can be explained by a combined process with outflow, rotation, convergent flow, collision, large velocity gradient, and rotating toroids. High deuterium fractionation strongly depends on the temperature condition. NH$_2$D is a poor evolutionary indicator of high-mass star formation in evolved stages, but a useful tracer in the starless and prestellar cores.


title: BURT: BERT-inspired Universal Representation from Twin Structure

abstract: Pre-trained contextualized language models such as BERT have shown great effectiveness in a wide range of downstream Natural Language Processing (NLP) tasks. However, the effective representations offered by the models target at each token inside a sequence rather than each sequence and the fine-tuning step involves the input of both sequences at one time, leading to unsatisfying representations of various sequences with different granularities. Especially, as sentence-level representations taken as the full training context in these models, there comes inferior performance on lower-level linguistic units (phrases and words). In this work, we present BURT (BERT inspired Universal Representation from Twin Structure) that is capable of generating universal, fixed-size representations for input sequences of any granularity, i.e., words, phrases, and sentences, using a large scale of natural language inference and paraphrase data with multiple training objectives. Our proposed BURT adopts the Siamese network, learning sentence-level representations from natural language inference dataset and word/phrase-level representations from paraphrasing dataset, respectively. We evaluate BURT across different granularities of text similarity tasks, including STS tasks, SemEval2013 Task 5(a) and some commonly used word similarity tasks, where BURT substantially outperforms other representation models on sentence-level datasets and achieves significant improvements in word/phrase-level representation.


title: Maximum entanglement of formation for a two-mode Gaussian state over passive operations

abstract: We quantify the maximum amount of entanglement of formation (EoF) that can be achieved by continuous-variable states under passive operations, which we refer to as EoF-potential. Focusing, in particular, on two-mode Gaussian states we derive analytical expressions for the EoF-potential for specific classes of states. For more general states, we demonstrate that this quantity can be upper-bounded by the minimum amount of squeezing needed to synthesize the Gaussian modes, a quantity called squeezing of formation. Our work, thus, provides a new link between non-classicality of quantum states and the non-classicality of correlations.


title: Prediction Model Based on Integrated Political Economy System: The Case of US Presidential Election

abstract: This paper studies an integrated system of political and economic systems from a systematic perspective to explore the complex interaction between them, and specially analyzes the case of the US presidential election forecasting. Based on the signed association networks of industrial structure constructed by economic data, our framework simulates the diffusion and evolution of opinions during the election through a kinetic model called the Potts Model. Remarkably, we propose a simple and efficient prediction model for the US presidential election, and meanwhile inspire a new way to model the economic structure. Findings also highlight the close relationship between economic structure and political attitude. Furthermore, the case analysis in terms of network and economy demonstrates the specific features and the interaction between political tendency and industrial structure in a particular period, which is consistent with theories in politics and economics.


title: Structures and Properties of $\beta$-Titanium Doping Trace Transition Metal Elements: a Density Functional Theory Study

abstract: We systematically calculate the structure, formation enthalpy, formation free energy, elastic constants and electronic structure of Ti$_{0.98}$X$_{0.02}$ system by density functional theory (DFT) simulations to explore the effect of transition metal X (X=Ag, Cd, Co, Cr, Cu, Fe, Mn, Mo, Nb, Ni, Pd, Rh, Ru, Tc, and Zn) on the stability mechanism of $\beta$-titanium. Based on our calculations, the results of formation enthalpy and free energy show that adding trace X is beneficial to the thermodynamic stability of $\beta$-titanium. This behavior is well explained by the density of state (DOS). However, the tetragonal shear moduli of Ti$_{0.98}$X$_{0.02}$ systems are negative, indicating that $\beta$-titanium doping with a low concentration of X is still elastically unstable at 0 K. Therefore, we theoretically explain that $\beta$-titanium doping with trace transition metal X is unstable in the ground state.


title: Near-infrared observations of active asteroid (3200) Phaethon reveal no evidence for hydration

abstract: Asteroid (3200) Phaethon is an active near-Earth asteroid and the parent body of the Geminid Meteor Shower. Because of its small perihelion distance, Phaethon's surface reaches temperatures sufficient to destabilize hydrated materials. We conducted rotationally resolved spectroscopic observations of this asteroid, mostly covering the northern hemisphere and the equatorial region, beyond 2.5-micron to search for evidence of hydration on its surface. Here we show that the observed part of Phaethon does not exhibit the 3-micron hydrated mineral absorption (within 2-sigma). These observations suggest that Phaethon's modern activity is not due to volatile sublimation or devolatilization of phyllosilicates on its surface. It is possible that the observed part of Phaethon was originally hydrated and has since lost volatiles from its surface via dehydration, supporting its connection to the Pallas family, or it was formed from anhydrous material.


title: Data Augmentation for Spoken Language Understanding via Pretrained Language Models

abstract: The training of spoken language understanding (SLU) models often faces the problem of data scarcity. In this paper, we put forward a data augmentation method using pretrained language models to boost the variability and accuracy of generated utterances. Furthermore, we investigate and propose solutions to two previously overlooked semi-supervised learning scenarios of data scarcity in SLU: i) Rich-in-Ontology: ontology information with numerous valid dialogue acts is given; ii) Rich-in-Utterance: a large number of unlabelled utterances are available. Empirical results show that our method can produce synthetic training data that boosts the performance of language understanding models in various scenarios.


title: Asymptotic Properties of High-Dimensional Random Forests

abstract: As a flexible nonparametric learning tool, the random forests algorithm has been widely applied to various real applications with appealing empirical performance, even in the presence of high-dimensional feature space. Unveiling the underlying mechanisms has led to some important recent theoretical results on the consistency of the random forests algorithm and its variants. However, to our knowledge, almost all existing works concerning random forests consistency in high dimensional setting were established for various modified random forests models where the splitting rules are independent of the response; a few exceptions assume simple data generating models with binary features. In light of this, in this paper we derive the consistency rates for the random forests algorithm associated with the sample CART splitting criterion, which is the one used in the original version of the algorithm, in a general high-dimensional nonparametric regression setting through a bias-variance decomposition analysis. Our new theoretical results show that random forests can indeed adapt to high dimensionality and allow for discontinuous regression function. Our bias analysis characterizes explicitly how the random forests bias depends on the sample size, tree height, and column subsampling parameter. Some limitations on our current results are also discussed.


title: Rethink the Connections among Generalization, Memorization and the Spectral Bias of DNNs

abstract: Over-parameterized deep neural networks (DNNs) with sufficient capacity to memorize random noise can achieve excellent generalization performance, challenging the bias-variance trade-off in classical learning theory. Recent studies claimed that DNNs first learn simple patterns and then memorize noise; some other works showed a phenomenon that DNNs have a spectral bias to learn target functions from low to high frequencies during training. However, we show that the monotonicity of the learning bias does not always hold: under the experimental setup of deep double descent, the high-frequency components of DNNs diminish in the late stage of training, leading to the second descent of the test error. Besides, we find that the spectrum of DNNs can be applied to indicating the second descent of the test error, even though it is calculated from the training set only.


title: Impressive optoelectronic and thermoelectric properties of two-dimensional XI$_2$ (X=Sn, Si): a first principle study

abstract: Two-dimensional (2D) metal halides have received more attention because of their electronic and optoelectronic properties. Recently, researchers are interested to investigate the thermoelectric properties of metal halide monolayers because of their ultralow lattice conductivity, high Seebeck coefficient and figure of merit. Here, we have investigated thermoelectric and optoelectronic properties of XI$_2$ (X=Sn and Si) monolayers with the help of density functional theory and Boltzmann transport equation. The structural parameters have been optimized with relaxation of atomic positions. Excellent thermoelectric and optical properties have been obtained for both SnI$_2$ and SiI$_2$ monolayers. For SnI$_2$ an indirect bandgap of 2.06 eV was observed and the absorption peak was found at 4.68 eV. For this the highest ZT value of 0.84 for p-type doping at 600K has been calculated. Similarly, for SiI$_2$ a comparatively low indirect bandgap of 1.63 eV was observed, and the absorption peak was obtained at 4.86 eV. The calculated ZT product for SiI$_2$ was 0.87 at 600K. Both the crystals having high absorbance and ZT value suggest that they can be promising candidates for optoelectronic and thermoelectric devices.


title: Zero-shot topic generation

abstract: We present an approach to generating topics using a model trained only for document title generation, with zero examples of topics given during training. We leverage features that capture the relevance of a candidate span in a document for the generation of a title for that document. The output is a weighted collection of the phrases that are most relevant for describing the document and distinguishing it within a corpus, without requiring access to the rest of the corpus. We conducted a double-blind trial in which human annotators scored the quality of our machine-generated topics along with original human-written topics associated with news articles from The Guardian and The Huffington Post. The results show that our zero-shot model generates topic labels for news documents that are on average equal to or higher quality than those written by humans, as judged by humans.


title: Continuous time digital search tree and a border aggregation model

abstract: We consider the continuous-time version of the random digital search tree, and construct a coupling with a border aggregation model as studied in Thacker and Volkov (2018), showing a relation between the height of the tree and the time required for aggregation. This relation carries over to the corresponding discrete-time models. As a consequence we find a very precise asymptotic result for the time to aggregation, using recent results by Drmota et al.\ (2020) for the digital search tree.


title: Optical parametric oscillation in silicon carbide nanophotonics

abstract: Silicon carbide (SiC) is rapidly emerging as a leading platform for the implementation of nonlinear and quantum photonics. Here, we find that commercial SiC, which hosts a variety of spin qubits, possesses low optical absorption that can enable SiC integrated photonics with quality factors exceeding $10^7$. We fabricate microring resonators with quality factors as high as 1.1 million, and observe low-threshold (8.5 $\pm$ 0.5 mW) optical parametric oscillation as well as optical microcombs spanning 200 nm. Our demonstration is an essential milestone in the development of photonic devices that harness the unique optical properties of SiC, paving the way toward the monolithic integration of nonlinear photonics with spin-based quantum technologies.


title: Video Contents Understanding using Deep Neural Networks

abstract: We propose a novel application of Transfer Learning to classify video-frame sequences over multiple classes. This is a pre-weighted model that does not require to train a fresh CNN. This representation is achieved with the advent of "deep neural network" (DNN), which is being studied these days by many researchers. We utilize the classical approaches for video classification task using object detection techniques for comparison, such as "Google Video Intelligence API" and this study will run experiments as to how those architectures would perform in foggy or rainy weather conditions. Experimental evaluation on video collections shows that the new proposed classifier achieves superior performance over existing solutions.


title: $n$-th parafermion $\mathcal{W}_N$ characters from $U(N)$ instanton counting on ${\mathbb {C}}^2/{\mathbb {Z}}_n$

abstract: We propose, following the AGT correspondence, how the $\mathcal{W}^{\, para}_{N, n}$ ($n$-th parafermion $\mathcal{W}_N$) minimal model characters are obtained from the $U(N)$ instanton counting on ${\mathbb {C}}^2/{\mathbb {Z}}_n$ with $\Omega$-deformation by imposing specific conditions which remove the minimal model null states.


title: Preconditioned Legendre spectral Galerkin methods for the non-separable elliptic equation

abstract: The Legendre spectral Galerkin method of self-adjoint second order elliptic equations usually results in a linear system with a dense and ill-conditioned coefficient matrix. In this paper, the linear system is solved by a preconditioned conjugate gradient (PCG) method where the preconditioner $M$ is constructed by approximating the variable coefficients with a ($T$+1)-term Legendre series in each direction to a desired accuracy. A feature of the proposed PCG method is that the iteration step increases slightly with the size of the resulting matrix when reaching a certain approximation accuracy. The efficiency of the method lies in that the system with the preconditioner $M$ is approximately solved by a one-step iterative method based on the ILU(0) factorization. The ILU(0) factorization of $M\in \mathbb{R}^{(N-1)^d\times(N-1)^d}$ can be computed using $\mathcal{O}(T^{2d} N^d)$ operations, and the number of nonzeros in the factorization factors is of $\mathcal{O}(T^{d} N^d)$, $d=1,2,3$. To further speed up the PCG method, an algorithm is developed for fast matrix-vector multiplications by the resulting matrix of Legendre-Galerkin spectral discretization, without the need to explicitly form it. The complexity of the fast matrix-vector multiplications is of $\mathcal{O}(N^d (\log N)^2)$. As a result, the PCG method has a $\mathcal{O}(N^d (\log N)^2)$ total complexity for a $d$ dimensional domain with $(N-1)^d$ unknows, $d=1,2,3$. Numerical examples are given to demonstrate the efficiency of proposed preconditioners and the algorithm for fast matrix-vector multiplications.


title: Energy Balancing of Covariate Distributions

abstract: Bias in causal comparisons has a direct correspondence with distributional imbalance of covariates between treatment groups. Weighting strategies such as inverse propensity score weighting attempt to mitigate bias by either modeling the treatment assignment mechanism or balancing specified covariate moments. This paper introduces a new weighting method, called energy balancing, which instead aims to balance weighted covariate distributions. By directly targeting distributional imbalance, the proposed weighting strategy can be flexibly utilized in a wide variety of causal analyses, including the estimation of average treatment effects and individualized treatment rules. Our energy balancing weights (EBW) approach has several advantages over existing weighting techniques. First, it offers a model-free and robust approach for obtaining covariate balance that does not require tuning parameters, obviating the need for modeling decisions of secondary nature to the scientific question at hand. Second, since this approach is based on a genuine measure of distributional balance, it provides a means for assessing the balance induced by a given set of weights for a given dataset. Finally, the proposed method is computationally efficient and has desirable theoretical guarantees under mild conditions. We demonstrate the effectiveness of this EBW approach in a suite of simulation experiments, and in studies on the safety of right heart catheterization and the effect of indwelling arterial catheters.


title: Optimal Study Design for Reducing Variances of Coefficient Estimators in Change-Point Models

abstract: In longitudinal studies, we observe measurements of the same variables at different time points to track the changes in their pattern over time. In such studies, scheduling of the data collection waves (i.e. time of participants' visits) is often pre-determined to accommodate ease of project management and compliance. Hence, it is common to schedule those visits at equally spaced time intervals. However, recent publications based on simulated experiments indicate that the power of studies and the precision of model parameter estimators is related to the participants' visiting schemes. In this paper, we consider the longitudinal studies that investigate the changing pattern of a disease outcome, (e.g. the accelerated cognitive decline of senior adults). Such studies are often analyzed by the broken-stick model, consisting of two segments of linear models connected at an unknown change-point. We formulate this design problem into a high-dimensional optimization problem and derive its analytical solution. Based on this solution, we propose an optimal design of the visiting scheme that maximizes the power (i.e. reduce the variance of estimators) to identify the onset of accelerated decline. Using both simulation studies and evidence from real data, we demonstrate our optimal design outperforms the standard equally-spaced design. Applying our novel design to plan the longitudinal studies, researchers can improve the power of detecting pattern change without collecting extra data.


title: Sequential Bayesian Parameter Estimation of Stochastic Dynamic Load Models

abstract: In this paper we focus on the parameter estimation of dynamic load models with stochastic terms, in particular, load models where protection settings are uncertain, such as in aggregated air conditioning units. We show how the uncertainty in the aggregated protection characteristics can be formulated as a stochastic differential equation with process noise. We cast the parameter inversion within a Bayesian parameter estimation framework, and we present methods to include process noise. We demonstrate the benefits of considering stochasticity in the parameter estimation and the risks of ignoring it.