title: DistDGL: Distributed Graph Neural Network Training for Billion-Scale Graphs

abstract: Graph neural networks (GNN) have shown great success in learning from graph-structured data. They are widely used in various applications, such as recommendation, fraud detection, and search. In these domains, the graphs are typically large, containing hundreds of millions of nodes and several billions of edges. To tackle this challenge, we develop DistDGL, a system for training GNNs in a mini-batch fashion on a cluster of machines. DistDGL is based on the Deep Graph Library (DGL), a popular GNN development framework. DistDGL distributes the graph and its associated data (initial features and embeddings) across the machines and uses this distribution to derive a computational decomposition by following an owner-compute rule. DistDGL follows a synchronous training approach and allows ego-networks forming the mini-batches to include non-local nodes. To minimize the overheads associated with distributed computations, DistDGL uses a high-quality and light-weight min-cut graph partitioning algorithm along with multiple balancing constraints. This allows it to reduce communication overheads and statically balance the computations. It further reduces the communication by replicating halo nodes and by using sparse embedding updates. The combination of these design choices allows DistDGL to train high-quality models while achieving high parallel efficiency and memory scalability. We demonstrate our optimizations on both inductive and transductive GNN models. Our results show that DistDGL achieves linear speedup without compromising model accuracy and requires only 13 seconds to complete a training epoch for a graph with 100 million nodes and 3 billion edges on a cluster with 16 machines. DistDGL is now publicly available as part of DGL:https://github.com/dmlc/dgl/tree/master/python/dgl/distributed.


title: We Can Detect Your Bias: Predicting the Political Ideology of News Articles

abstract: We explore the task of predicting the leading political ideology or bias of news articles. First, we collect and release a large dataset of 34,737 articles that were manually annotated for political ideology -left, center, or right-, which is well-balanced across both topics and media. We further use a challenging experimental setup where the test examples come from media that were not seen during training, which prevents the model from learning to detect the source of the target news article instead of predicting its political ideology. From a modeling perspective, we propose an adversarial media adaptation, as well as a specially adapted triplet loss. We further add background information about the source, and we show that it is quite helpful for improving article-level prediction. Our experimental results show very sizable improvements over using state-of-the-art pre-trained Transformers in this challenging setup.


title: Collision free motion planning on a wedge of circles

abstract: We exhibit an algorithm with continuous instructions for two robots moving without collisions on a track shaped as a wedge of three circles. We show that the topological complexity of the configuration space associated with this problem is 3. The topological complexity is a homotopy invariant that can be thought of as the minimum number of continuous instructions required to describe the movement of the robots between any initial configuration to any final one without collisions. The algorithm presented is optimal in the sense that it requires exactly 3 continuous instructions.


title: Self-attention aggregation network for video face representation and recognition

abstract: Models based on self-attention mechanisms have been successful in analyzing temporal data and have been widely used in the natural language domain. We propose a new model architecture for video face representation and recognition based on a self-attention mechanism. Our approach could be used for video with single and multiple identities. To the best of our knowledge, no one has explored the aggregation approaches that consider the video with multiple identities. The proposed approach utilizes existing models to get the face representation for each video frame, e.g., ArcFace and MobileFaceNet, and the aggregation module produces the aggregated face representation vector for video by taking into consideration the order of frames and their quality scores. We demonstrate empirical results on a public dataset for video face recognition called IJB-C to indicate that the self-attention aggregation network (SAAN) outperforms naive average pooling. Moreover, we introduce a new multi-identity video dataset based on the publicly available UMDFaces dataset and collected GIFs from Giphy. We show that SAAN is capable of producing a compact face representation for both single and multiple identities in a video. The dataset and source code will be publicly available.


title: On Determining and Qualifying the Number of Superstates in Aggregation of Markov Chains

abstract: Many studies involving large Markov chains require determining a smaller representative (aggregated) chains. Each {\em superstate} in the representative chain represents a {\em group of related} states in the original Markov chain. Typically, the choice of number of superstates in the aggregated chain is ambiguous, and based on the limited prior know-how. In this paper we present a structured methodology of determining the best candidate for the number of superstates. We achieve this by comparing aggregated chains of different sizes. To facilitate this comparison we develop and quantify a notion of {\em marginal return}. Our notion captures the decrease in the {\em heterogeneity} within the group of the {\em related} states (i.e., states represented by the same superstate) upon a unit increase in the number of superstates in the aggregated chain. We use Maximum Entropy Principle to justify the notion of marginal return, as well as our quantification of heterogeneity. Through simulations on synthetic Markov chains, where the number of superstates are known apriori, we show that the aggregated chain with the largest marginal return identifies this number. In case of Markov chains that model real-life scenarios we show that the aggregated model with the largest marginal return identifies an inherent structure unique to the scenario being modelled; thus, substantiating on the efficacy of our proposed methodology.


title: Using Information to Amplify Competition

abstract: I characterize the consumer-optimal market segmentation in competitive markets where multiple firms selling differentiated products to consumers with unit demand. This segmentation is public---in that each firm observes the same market segments---and takes a simple form: in each market segment, there is a dominant firm favored by all consumers in that segment. By segmenting the market, all but the dominant firm maximally compete to poach the consumer's business, setting price to equal marginal cost. Information, thus, is being used to amplify competition. This segmentation simultaneously generates an efficient allocation and delivers to each firm its minimax profit.


title: Non-Stationary Stochastic Global Optimization Algorithms

abstract: Gomez proposes a formal and systematic approach for characterizing stochastic global optimization algorithms. Using it, Gomez formalizes algorithms with a fixed next-population stochastic method, i.e., algorithms defined as stationary Markov processes. These are the cases of standard versions of hill-climbing, parallel hill-climbing, generational genetic, steady-state genetic, and differential evolution algorithms. This paper continues such a systematic formal approach. First, we generalize the sufficient conditions convergence lemma from stationary to non-stationary Markov processes. Second, we develop Markov kernels for some selection schemes. Finally, we formalize both simulated-annealing and evolutionary-strategies using the systematic formal approach.


title: ASSURE: RTL Locking Against an Untrusted Foundry

abstract: Semiconductor design companies are integrating proprietary intellectual property (IP) blocks to build custom integrated circuits (IC) and fabricate them in a third-party foundry. Unauthorized IC copies cost these companies billions of dollars annually. While several methods have been proposed for hardware IP obfuscation, they operate on the gate-level netlist, i.e., after the synthesis tools embed the semantic information into the netlist. We propose ASSURE to protect hardware IP modules operating on the register-transfer level (RTL) description. The RTL approach has three advantages: (i) it allows designers to obfuscate IP cores generated with many different methods (e.g., hardware generators, high-level synthesis tools, and pre-existing IPs). (ii) it obfuscates the semantics of an IC before logic synthesis; (iii) it does not require modifications to EDA flows. We perform a cost and security assessment of ASSURE.


title: Do Language Embeddings Capture Scales?

abstract: Pretrained Language Models (LMs) have been shown to possess significant linguistic, common sense, and factual knowledge. One form of knowledge that has not been studied yet in this context is information about the scalar magnitudes of objects. We show that pretrained language models capture a significant amount of this information but are short of the capability required for general common-sense reasoning. We identify contextual information in pre-training and numeracy as two key factors affecting their performance and show that a simple method of canonicalizing numbers can have a significant effect on the results.


title: Explicit universal minimal constants for polynomial growth of groups

abstract: Shalom and Tao showed that a polynomial upper bound on the size of a single, large enough ball in a Cayley graph implies that the underlying group has a nilpotent subgroup with index and degree of polynomial growth both bounded effectively. The third and fourth authors proved the optimal bound on the degree of polynomial growth of this subgroup, at the expense of making some other parts of the result ineffective. In the present paper we prove the optimal bound on the degree of polynomial growth without making any losses elsewhere. As a consequence, we show that there exist explicit positive numbers $\varepsilon_d$ such that in any group with growth at least a polynomial of degree $d$, the growth is at least $\varepsilon_dn^d$. We indicate some applications in probability; in particular, we show that the gap at $1$ for the critical probability for Bernoulli site percolation on a Cayley graph, recently proven to exist by Panagiotis and Severo, is at least $\exp\bigl\{-\exp\bigl\{17 \exp\{100 \cdot 8^{100}\}\bigr\}\bigr\}$.


title: Optical analogue of Dresselhaus spin-orbit interaction in photonic graphene

abstract: The concept of gauge fields plays a significant role in many areas of physics from particle physics and cosmology to condensed matter systems, where gauge potentials are a natural consequence of electromagnetic fields acting on charged particles and are of central importance in topological states of matter. Here, we report on the experimental realization of a synthetic non-Abelian gauge field for photons in a honeycomb microcavity lattice. We show that the effective magnetic field associated with TE-TM splitting has the symmetry of Dresselhaus spin-orbit interaction around Dirac points in the dispersion, and can be regarded as an SU(2) gauge field. The symmetry of the field is revealed in the optical spin Hall effect (OSHE), where under resonant excitation of the Dirac points precession of the photon pseudospin around the field direction leads to the formation of two spin domains. Furthermore, we observe that the Dresselhaus field changes its sign in the same Dirac valley on switching from s to p bands in good agreement with the tight binding modelling. Our work demonstrating a non-Abelian gauge field for light on the microscale paves the way towards manipulation of photons via spin on a chip.


title: Automatic Particle Trajectory Classification in Plasma Simulations

abstract: Numerical simulations of plasma flows are crucial for advancing our understanding of microscopic processes that drive the global plasma dynamics in fusion devices, space, and astrophysical systems. Identifying and classifying particle trajectories allows us to determine specific on-going acceleration mechanisms, shedding light on essential plasma processes. Our overall goal is to provide a general workflow for exploring particle trajectory space and automatically classifying particle trajectories from plasma simulations in an unsupervised manner. We combine pre-processing techniques, such as Fast Fourier Transform (FFT), with Machine Learning methods, such as Principal Component Analysis (PCA), k-means clustering algorithms, and silhouette analysis. We demonstrate our workflow by classifying electron trajectories during magnetic reconnection problem. Our method successfully recovers existing results from previous literature without a priori knowledge of the underlying system. Our workflow can be applied to analyzing particle trajectories in different phenomena, from magnetic reconnection, shocks to magnetospheric flows. The workflow has no dependence on any physics model and can identify particle trajectories and acceleration mechanisms that were not detected before.


title: ComStreamClust: a communicative multi-agent approach to text clustering in streaming data

abstract: Topic detection is the task of determining and tracking hot topics in social media. Twitter is arguably the most popular platform for people to share their ideas with others about different issues. One such prevalent issue is the COVID-19 pandemic. Detecting and tracking topics on these kinds of issues would help governments and healthcare companies deal with this phenomenon. In this paper, we propose a novel, multi-agent, communicative clustering approach, so-called ComStreamClust for clustering sub-topics inside a broader topic, e.g., COVID-19. The proposed approach is parallelizable, and can simultaneously handle several data-point. The LaBSE sentence embedding is used to measure the semantic similarity between two tweets. ComStreamClust has been evaluated on two datasets: the COVID-19 and the FA CUP. The results obtained from ComStreamClust approve the effectiveness of the proposed approach when compared to existing methods.


title: Google Landmark Recognition 2020 Competition Third Place Solution

abstract: We present our third place solution to the Google Landmark Recognition 2020 competition. It is an ensemble of global features only Sub-center ArcFace models. We introduce dynamic margins for ArcFace loss, a family of tune-able margin functions of class size, designed to deal with the extreme imbalance in GLDv2 dataset. Progressive finetuning and careful postprocessing are also key to the solution. Our two submissions scored 0.6344 and 0.6289 on private leaderboard, both ranking third place out of 736 teams.


title: Identifying Melanoma Images using EfficientNet Ensemble: Winning Solution to the SIIM-ISIC Melanoma Classification Challenge

abstract: We present our winning solution to the SIIM-ISIC Melanoma Classification Challenge. It is an ensemble of convolutions neural network (CNN) models with different backbones and input sizes, most of which are image-only models while a few of them used image-level and patient-level metadata. The keys to our winning are: (1) stable validation scheme (2) good choice of model target (3) carefully tuned pipeline and (4) ensembling with very diverse models. The winning submission scored 0.9600 AUC on cross validation and 0.9490 AUC on private leaderboard.


title: MoCo-CXR: MoCo Pretraining Improves Representation and Transferability of Chest X-ray Models

abstract: Contrastive learning is a form of self-supervision that can leverage unlabeled data to produce pretrained models. While contrastive learning has demonstrated promising results on natural image classification tasks, its application to medical imaging tasks like chest X-ray interpretation has been limited. In this work, we propose MoCo-CXR, which is an adaptation of the contrastive learning method Momentum Contrast (MoCo), to produce models with better representations and initializations for the detection of pathologies in chest X-rays. In detecting pleural effusion, we find that linear models trained on MoCo-CXR-pretrained representations outperform those without MoCo-CXR-pretrained representations, indicating that MoCo-CXR-pretrained representations are of higher-quality. End-to-end fine-tuning experiments reveal that a model initialized via MoCo-CXR-pretraining outperforms its non-MoCo-CXR-pretrained counterpart. We find that MoCo-CXR-pretraining provides the most benefit with limited labeled training data. Finally, we demonstrate similar results on a target Tuberculosis dataset unseen during pretraining, indicating that MoCo-CXR-pretraining endows models with representations and transferability that can be applied across chest X-ray datasets and tasks.


title: Local Connectivity in Centroid Clustering

abstract: Clustering is a fundamental task in unsupervised learning, one that targets to group a dataset into clusters of similar objects. There has been recent interest in embedding normative considerations around fairness within clustering formulations. In this paper, we propose 'local connectivity' as a crucial factor in assessing membership desert in centroid clustering. We use local connectivity to refer to the support offered by the local neighborhood of an object towards supporting its membership to the cluster in question. We motivate the need to consider local connectivity of objects in cluster assignment, and provide ways to quantify local connectivity in a given clustering. We then exploit concepts from density-based clustering and devise LOFKM, a clustering method that seeks to deepen local connectivity in clustering outputs, while staying within the framework of centroid clustering. Through an empirical evaluation over real-world datasets, we illustrate that LOFKM achieves notable improvements in local connectivity at reasonable costs to clustering quality, illustrating the effectiveness of the method.


title: Black hole or Gravastar? The GW190521 case

abstract: The existence of cosmological compact objects with very strong gravity is a prediction of General Relativity and an exact solution of the Einstein equations. These objects are called black holes and recently we had the first observations of them. However, the theory of black hole formation has some disadvantages. In order to avoid these, some scientists suggest the existence of gravastars (gravitation vacuum stars), an alternative stellar model which seems to solve the problems of the black hole theory. In this work we compare black holes and gravastars using a wide range of the literature and we emphasize the properties of gravastars, which are consistent with the current cosmological observations. Also, we propose gravastars as the solution of the "pair-instability" effect and a possible explanation for the observed masses of the compact objects, before the collapse, from the gravitational signal GW190521, since in the formation of a gravastar there aren't mass restrictions.


title: Medical Image Harmonization Using Deep Learning Based Canonical Mapping: Toward Robust and Generalizable Learning in Imaging

abstract: Conventional and deep learning-based methods have shown great potential in the medical imaging domain, as means for deriving diagnostic, prognostic, and predictive biomarkers, and by contributing to precision medicine. However, these methods have yet to see widespread clinical adoption, in part due to limited generalization performance across various imaging devices, acquisition protocols, and patient populations. In this work, we propose a new paradigm in which data from a diverse range of acquisition conditions are "harmonized" to a common reference domain, where accurate model learning and prediction can take place. By learning an unsupervised image to image canonical mapping from diverse datasets to a reference domain using generative deep learning models, we aim to reduce confounding data variation while preserving semantic information, thereby rendering the learning task easier in the reference domain. We test this approach on two example problems, namely MRI-based brain age prediction and classification of schizophrenia, leveraging pooled cohorts of neuroimaging MRI data spanning 9 sites and 9701 subjects. Our results indicate a substantial improvement in these tasks in out-of-sample data, even when training is restricted to a single site.


title: Frequency Multiplexed Optical Entangled Source based on the Pockels Effect

abstract: In the recent years important experimental advances in resonant electro-optic modulators as high efficient sources for coherent frequency combs and as devices for quantum information transfer have been realized, where strong optical and microwave mode coupling were achieved. These features suggest electro-optic based devices as candidates for entangled optical frequency comb sources. In the present I study the generation of entangled optical frequency combs in mm-sized resonant electro-optic modulators. These devices profit from the experimentally proven advantages such as nearly constant optical free spectral ranges over several gigahertz, high optical and microwave quality factors. The generation of frequency multiplexed quantum channels with spectral bandwidth in the range of a MHz for conservative parameter values paves the way towards novel uses in long distant hybrid quantum networks, quantum key distribution, enhanced optical metrology and quantum computing.


title: A Knowledge-Driven Approach to Classifying Object and Attribute Coreferences in Opinion Mining

abstract: Classifying and resolving coreferences of objects (e.g., product names) and attributes (e.g., product aspects) in opinionated reviews is crucial for improving the opinion mining performance. However, the task is challenging as one often needs to consider domain-specific knowledge (e.g., iPad is a tablet and has aspect resolution) to identify coreferences in opinionated reviews. Also, compiling a handcrafted and curated domain-specific knowledge base for each domain is very time consuming and arduous. This paper proposes an approach to automatically mine and leverage domain-specific knowledge for classifying objects and attribute coreferences. The approach extracts domain-specific knowledge from unlabeled review data and trains a knowledgeaware neural coreference classification model to leverage (useful) domain knowledge together with general commonsense knowledge for the task. Experimental evaluation on realworld datasets involving five domains (product types) shows the effectiveness of the approach.


title: Learning Which Features Matter: RoBERTa Acquires a Preference for Linguistic Generalizations (Eventually)

abstract: One reason pretraining on self-supervised linguistic tasks is effective is that it teaches models features that are helpful for language understanding. However, we want pretrained models to learn not only to represent linguistic features, but also to use those features preferentially during fine-turning. With this goal in mind, we introduce a new English-language diagnostic set called MSGS (the Mixed Signals Generalization Set), which consists of 20 ambiguous binary classification tasks that we use to test whether a pretrained model prefers linguistic or surface generalizations during fine-tuning. We pretrain RoBERTa models from scratch on quantities of data ranging from 1M to 1B words and compare their performance on MSGS to the publicly available RoBERTa-base. We find that models can learn to represent linguistic features with little pretraining data, but require far more data to learn to prefer linguistic generalizations over surface ones. Eventually, with about 30B words of pretraining data, RoBERTa-base does demonstrate a linguistic bias with some regularity. We conclude that while self-supervised pretraining is an effective way to learn helpful inductive biases, there is likely room to improve the rate at which models learn which features matter.


title: Inverse Dynamics vs. Forward Dynamics in Direct Transcription Formulations for Trajectory Optimization

abstract: Benchmarks of state-of-the-art rigid-body dynamics libraries report better performance solving the inverse dynamics problem than the forward alternative. Those benchmarks encouraged us to question whether that computational advantage would translate to direct transcription, where calculating rigid-body dynamics and their derivatives accounts for a significant share of computation time. In this work, we implement an optimization framework where both approaches for enforcing the system dynamics are available. We evaluate the performance of each approach for systems of varying complexity, for domains with rigid contacts. Our tests reveal that formulations using inverse dynamics converge faster, require less iterations, and are more robust to coarse problem discretization. These results indicate that inverse dynamics should be preferred to enforce the nonlinear system dynamics in simultaneous methods, such as direct transcription.


title: A range characterization of the single-quadrant ADRT

abstract: This work characterizes the range of the single-quadrant approximate discrete Radon transform (ADRT) of square images. The characterization follows from a set of linear constraints on the codomain. We show that for data satisfying these constraints, the exact and fast inversion formula [Rim, Appl. Math. Lett. 102 106159, 2020] yields a square image in a stable manner. The range characterization is obtained by first showing that the ADRT is a bijection between images supported on infinite half-strips, then identifying the linear subspaces that stay finitely supported under the inversion formula.


title: Local Models For Rapoport-Zink Spaces For Local P-Shtukas

abstract: This article provides a ``local'' complementary to the previous results concerning the local models for the moduli stacks of ``global'' $G$-shtukas. Here we study the geometry of Rapoport-Zink spaces for local $P$-shtukas by constructing local models for them. We further discuss certain applications, including some results related to the theory of formal nearby cycles associated to these spaces and the semi-simple trace of Frobenius on the corresponding sheaves.


title: Superconducting junction with tri-component pairing gap functions

abstract: We study a superconducting hetro-junction with one side characterized by the unconventional chiral $p$-wave gap function $p_x\pm ip_y$ and the other side the conventional $s$-wave one. Though a relative phase of $\pm \frac{\pi}{2}$ between any two components of gap functions is favored in the junction region, mutual phase differences cannot achieve $\pm \frac{\pi}{2}$ simultaneously, which results in frustration. Based on a Ginzburg-Landau free energy analysis, the frustrated pattern is determined to be $s+ i\eta_1 (e^{ i\eta_2 \varphi/2}p_x +\eta_3 e^{- i\eta_2 \varphi/2}p_y)$ with $\eta_j=\pm 1$ ($j=1,2,3$), where $\varphi$ is the phase difference between the $p_x$- and $p_y$-wave gap functions. Furthermore, we find that the junction exhibits an anisotropic magnetoelectric effect, manifesting itself as an anisotropic spin magnetization along the edge of the junction.


title: Finite Rigid Sets in Flip Graphs

abstract: We show that for most pairs of surfaces, there exists a finite subgraph of the flip graph of the first surface so that any injective homomorphism of this finite subgraph into the flip graph of the second surface can be extended uniquely to an injective homomorphism between the two flip graphs. Combined with a result of Aramayona-Koberda-Parlier, this implies that any such injective homomorphism of this finite set is induced by an embedding of the surfaces. We also include images of several flip graphs in an appendix.


title: Global solvability and propagation of regularity of sums of squares on compact manifolds

abstract: We investigate global solvability, in the framework of smooth functions and Schwartz distributions, of certain sums of squares of vector fields defined on a product of compact Riemannian manifolds $T \times G$, where $G$ is further assumed to be a Lie group. As in a recent article due to the authors, our analysis is carried out in terms of a system of left-invariant vector fields on $G$ naturally associated with the operator under study, a simpler object which nevertheless conveys enough information about the original operator so as to fully encode its solvability. As a welcome side effect of the tools developed for our main purpose, we easily prove a general result on propagation of regularity for such operators.


title: ArXiving Before Submission Helps Everyone

abstract: We claim, and present evidence, that allowing arXiv publication before a conference or journal submission benefits researchers, especially early career, as well as the whole scientific community. Specifically, arXiving helps professional identity building, protects against independent re-discovery, idea theft and gate-keeping; it facilitates open research result distribution and reduces inequality. The advantages dwarf the drawbacks -- mainly the relative increase in acceptance rate of papers of well-known authors -- which studies show to be marginal. Analyzing the pros and cons of arXiving papers, we conclude that requiring preprints be anonymous is nearly as detrimental as not allowing them. We see no reasons why anyone but the authors should decide whether to arXiv or not.


title: Canonical connections on sub-Riemannian manifolds with constant symbol

abstract: We give formulas for a canonical choice of grading and compatible affine connection on a sub-Riemannian manifold. This grading and affine connection is available on any sub-Riemannian manifold with constant symbol and is based on Morimoto's normalization condition for a Cartan connection in such a setting. We completely compute these structures for contact manifolds of constant symbol, including the cases where the connections of Tanaka-Webster-Tanno are not defined. We also give an original intrinsic grading on sub-Riemannian (2,3,5)-manifolds. We use this grading to give a flatness theorem for such manifolds and explain how the grading and connection from Morimoto's theory compare.


title: On the Colloidal Phase of the Homogeneous Electron Fluid

abstract: We provide semi-rigorous arguments that the Homogeneous Electron Fluid (HEF) has a colloidal phase separating the Wigner Crystal from the high density fluid phase. Near the crossover between crystal and fluid ground state energies, the argument is quite general and valid for practically any quantum transition between a crystal and a more amorphous phase. In this regime, the colloid is a gel and its "Goldstone" modes are flows of irregular fluid droplets separated by crystalline walls. A metal insulator transition occurs when a single bubble of fluid spans the entire system. Beyond this transition the colloid is a sol and its properties depend on the existence of meta-stable finite crystallites with negative surface tension. If these exist, the sol phase has lower energy than the homogeneous fluid. In the two dimensional HEF, Kivelson and Spivak argued that such negative surface tension objects always exist, at least in the form of stripes. We argue that the existence of stripes also implies the existence of finite negative tension elliptical crystallites and the detailed competition between the finite crystallite phase and striped phases is difficult to calculate. We also provide weaker arguments that finite crystallites exist in three dimensions. This provides evidence for our claim that the gapless excitations at non-zero wavenumber, observed in the numerical calculations of\cite{haule}\cite{gapless} are quantum remnants of these crystallites (Bosonic quasi-particles) in the limit of vanishing surface tension. Finally, we suggest a Landau mean field theory for the second order quantum phase transition between the fluid and sol phases.


title: Existence of solutions for deterministic bilevel games under a general Bayesian approach

abstract: In 1996, Mallozzi and Morgan [33] proposed a new model for Stackelberg games which we refer here to as the Bayesian approach. The leader has only partial information about how followers select their reaction among possibly multiple optimal ones. This partial information is modeled as a decision-dependent distribution, the so-called belief of the leader. In this work, we formalize the setting of this approach for bilevel games admitting multiple leaders and we provide new results of existence of solutions. We pay particular attention to the fundamental case of linear bilevel problems, which has not been studied before, and which main difficulty is given by possible variations in the dimension of the reaction set of the follower. Our main technique to address this difficulty is based on a stronger notion of continuity for set-valued maps that we call rectangular continuity, and which is verified by the solution set of parametric linear problems. Finally, we provide some numerical experiments to address linear bilevel problems under the Bayesian approach.


title: Quantitative Argument Summarization and Beyond: Cross-Domain Key Point Analysis

abstract: When summarizing a collection of views, arguments or opinions on some topic, it is often desirable not only to extract the most salient points, but also to quantify their prevalence. Work on multi-document summarization has traditionally focused on creating textual summaries, which lack this quantitative aspect. Recent work has proposed to summarize arguments by mapping them to a small set of expert-generated key points, where the salience of each key point corresponds to the number of its matching arguments. The current work advances key point analysis in two important respects: first, we develop a method for automatic extraction of key points, which enables fully automatic analysis, and is shown to achieve performance comparable to a human expert. Second, we demonstrate that the applicability of key point analysis goes well beyond argumentation data. Using models trained on publicly available argumentation datasets, we achieve promising results in two additional domains: municipal surveys and user reviews. An additional contribution is an in-depth evaluation of argument-to-key point matching models, where we substantially outperform previous results.


title: Game-theoric approach to decision-making problem for blockchain mining

abstract: It is an important decision-making problem for a miner in the blockchain networks if he/she participates in the mining so that he/she earns a reward by creating a new block earlier than other miners. We formulate this decision-making problem as a noncooperative game, because the probability of creating a block depends not only on one's own available computational resources, but also those of other miners. Through theoretical and numerical analyses, we show a hysteresis phenomenon of Nash equilibria depending on the reward and a jump phenomenon of miner decisions by a slight change in reward. We also show that the reward for which miners decide not to participate in the mining becomes smaller as the number of miners increases.


title: Early Abandoning PrunedDTW and its application to similarity search

abstract: The Dynamic Time Warping ("DTW") distance is widely used in time series analysis, be it for classification, clustering or similarity search. However, its quadratic time complexity prevents it from scaling. Strategies, based on early abandoning DTW or skipping its computation altogether thanks to lower bounds, have been developed for certain use cases such as nearest neighbour search. But vectorization and approximation aside, no advance was made on DTW itself until recently with the introduction of PrunedDTW. This algorithm, able to prune unpromising alignments, was later fitted with early abandoning. We present a new version of PrunedDTW, "EAPrunedDTW", designed with early abandon in mind from the start, and able to early abandon faster than before. We show that EAPrunedDTW significantly improves the computation time of similarity search in the UCR Suite, and renders lower bounds dispensable.


title: Electrical properties of thermal oxide scales on pure iron in liquid lead-bismuth eutectic

abstract: The impedance behavior of pre-oxidized iron in liquid lead-bismuth eutectic (LBE) at 200 oC is studied using electrochemical impedance spectroscopy. The structures and resistance of oxide grown on iron oxidized in air at different temperatures and durations are compared. The results show that the resistance of the oxide film increases with increasing oxidizing temperature, due to the formation of a thicker scale and fewer defects. At the same temperature (600 oC), increasing the oxidation time can also reduce the defect concentration in the oxide film and improve the impedance of the oxide scale in LBE.


title: Distributionally Robust Local Non-parametric Conditional Estimation

abstract: Conditional estimation given specific covariate values (i.e., local conditional estimation or functional estimation) is ubiquitously useful with applications in engineering, social and natural sciences. Existing data-driven non-parametric estimators mostly focus on structured homogeneous data (e.g., weakly independent and stationary data), thus they are sensitive to adversarial noise and may perform poorly under a low sample size. To alleviate these issues, we propose a new distributionally robust estimator that generates non-parametric local estimates by minimizing the worst-case conditional expected loss over all adversarial distributions in a Wasserstein ambiguity set. We show that despite being generally intractable, the local estimator can be efficiently found via convex optimization under broadly applicable settings, and it is robust to the corruption and heterogeneity of the data. Experiments with synthetic and MNIST datasets show the competitive performance of this new class of estimators.


title: The generating pairs of the 2-transitive groups

abstract: Given a finite group $G$. The generating pair $(H,a)$ of $G$, that is, $H<G$ and $a\in G$ such that $\langle a,H\rangle=G$. In this paper, we introduce the definition of FF-subgroup to characterize the generating pairs of the symmetric groups, alternating groups and projective groups $PSL(2,q)$. This gives a partial answer to an open problem of J. Andr\'e and J. Ara\'ujo and P. J. Cameron.


title: Causal learning with sufficient statistics: an information bottleneck approach

abstract: The inference of causal relationships using observational data from partially observed multivariate systems with hidden variables is a fundamental question in many scientific domains. Methods extracting causal information from conditional independencies between variables of a system are common tools for this purpose, but are limited in the lack of independencies. To surmount this limitation, we capitalize on the fact that the laws governing the generative mechanisms of a system often result in substructures embodied in the generative functional equation of a variable, which act as sufficient statistics for the influence that other variables have on it. These functional sufficient statistics constitute intermediate hidden variables providing new conditional independencies to be tested. We propose to use the Information Bottleneck method, a technique commonly applied for dimensionality reduction, to find underlying sufficient sets of statistics. Using these statistics we formulate new additional rules of causal orientation that provide causal information not obtainable from standard structure learning algorithms, which exploit only conditional independencies between observable variables. We validate the use of sufficient statistics for structure learning both with simulated systems built to contain specific sufficient statistics and with benchmark data from regulatory rules previously and independently proposed to model biological signal transduction networks.


title: Ambiguous Persuasion: An Ex-Ante Formulation

abstract: Consider a persuasion game where both the sender and receiver are ambiguity averse with maxmin expected utility (MEU) preferences and the sender can choose to design an ambiguous information structure. This paper studies the game with an ex-ante formulation: The sender first commits to a (possibly ambiguous) information structure and then the receiver best responds by choosing an ex-ante message-contingent action plan. Under this formulation, I show it is never strictly beneficial for the sender to use an ambiguous information structure as opposed to a standard (unambiguous) information structure. This result is shown to be robust to the receiver having non-MEU Uncertainty Averse preferences but not to the sender having non-MEU preferences.


title: Koopman Operator, Geometry, and Learning

abstract: We provide a framework for learning of dynamical systems rooted in the concept of representations and Koopman operators. The interplay between the two leads to the full description of systems that can be represented linearly in a finite dimension, based on the properties of the Koopman operator spectrum. The geometry of state space is connected to the notion of representation, both in the linear case - where it is related to joint level sets of eigenfunctions - and in the nonlinear representation case. As shown here, even nonlinear finite-dimensional representations can be learned using the Koopman operator framework, leading to a new class of representation eigenproblems. The connection to learning using neural networks is given. An extension of the Koopman operator theory to "static" maps between different spaces is provided. The effect of the Koopman operator spectrum on Mori-Zwanzig type representations is discussed.


title: Rectification of interleavings and a persistent Whitehead theorem

abstract: The homotopy interleaving distance, a distance between persistent spaces, was introduced by Blumberg and Lesnick and shown to be universal, in the sense that it is the largest homotopy-invariant distance for which sublevel-set filtrations of close-by real-valued functions are close-by. There are other ways of constructing homotopy-invariant distances, but not much is known about the relationships between these choices. We show that other natural distances differ from the homotopy interleaving distance in at most a multiplicative constant, and prove versions of the persistent Whitehead theorem, a conjecture of Blumberg and Lesnick that relates morphisms that induce interleavings in persistent homotopy groups to stronger homotopy-invariant notions of interleaving.


title: MAF: Multimodal Alignment Framework for Weakly-Supervised Phrase Grounding

abstract: Phrase localization is a task that studies the mapping from textual phrases to regions of an image. Given difficulties in annotating phrase-to-object datasets at scale, we develop a Multimodal Alignment Framework (MAF) to leverage more widely-available caption-image datasets, which can then be used as a form of weak supervision. We first present algorithms to model phrase-object relevance by leveraging fine-grained visual representations and visually-aware language representations. By adopting a contrastive objective, our method uses information in caption-image pairs to boost the performance in weakly-supervised scenarios. Experiments conducted on the widely-adopted Flickr30k dataset show a significant improvement over existing weakly-supervised methods. With the help of the visually-aware language representations, we can also improve the previous best unsupervised result by 5.56%. We conduct ablation studies to show that both our novel model and our weakly-supervised strategies significantly contribute to our strong results.


title: Efficient Wasserstein Natural Gradients for Reinforcement Learning

abstract: A novel optimization approach is proposed for application to policy gradient methods and evolution strategies for reinforcement learning (RL). The procedure uses a computationally efficient Wasserstein natural gradient (WNG) descent that takes advantage of the geometry induced by a Wasserstein penalty to speed optimization. This method follows the recent theme in RL of including a divergence penalty in the objective to establish a trust region. Experiments on challenging tasks demonstrate improvements in both computational cost and performance over advanced baselines.


title: Torsion Subgroups of Groups with Quadratic Dehn Function

abstract: We construct the first examples of finitely presented groups with quadratic Dehn function containing a finitely generated infinite torsion subgroup. These examples are "optimal" in the sense that the Dehn function of any such finitely presented group must be at least quadratic. Moreover, we show that for any $n\geq2^{48}$ such that $n$ is either odd or divisible by $2^9$, any infinite free Burnside group with exponent $n$ is a quasi-isometrically embedded subgroup of a finitely presented group with quadratic Dehn function satisfying the Congruence Extension Property.


title: Miniscope3D: optimized single-shot miniature 3D fluorescence microscopy

abstract: Miniature fluorescence microscopes are a standard tool in systems biology. However, widefield miniature microscopes capture only 2D information, and modifications that enable 3D capabilities increase the size and weight and have poor resolution outside a narrow depth range. Here, we achieve the 3D capability by replacing the tube lens of a conventional 2D Miniscope with an optimized multifocal phase mask at the objective's aperture stop. Placing the phase mask at the aperture stop significantly reduces the size of the device, and varying the focal lengths enables a uniform resolution across a wide depth range. The phase mask encodes the 3D fluorescence intensity into a single 2D measurement, and the 3D volume is recovered by solving a sparsity-constrained inverse problem. We provide methods for designing and fabricating the phase mask and an efficient forward model that accounts for the field-varying aberrations in miniature objectives. We demonstrate a prototype that is 17 mm tall and weighs 2.5 grams, achieving 2.76 $\mu$m lateral, and 15 $\mu$m axial resolution across most of the 900x700x390 $\mu m^3$ volume at 40 volumes per second. The performance is validated experimentally on resolution targets, dynamic biological samples, and mouse brain tissue. Compared with existing miniature single-shot volume-capture implementations, our system is smaller and lighter and achieves a more than 2x better lateral and axial resolution throughout a 10x larger usable depth range. Our microscope design provides single-shot 3D imaging for applications where a compact platform matters, such as volumetric neural imaging in freely moving animals and 3D motion studies of dynamic samples in incubators and lab-on-a-chip devices.


title: Landau Quantization and Highly Mobile Fermions in an Insulator

abstract: In strongly correlated materials, quasiparticle excitations can carry fractional quantum numbers. An intriguing possibility is the formation of fractionalized, charge-neutral fermions, e.g., spinons and fermionic excitons, that result in neutral Fermi surfaces and Landau quantization in an insulator. While previous experiments in quantum spin liquids, topological Kondo insulators, and quantum Hall systems have hinted at charge-neutral Fermi surfaces, evidence for their existence remains far from conclusive. Here we report experimental observation of Landau quantization in a two dimensional (2D) insulator, i.e., monolayer tungsten ditelluride (WTe$_{2}$), a large gap topological insulator. Using a detection scheme that avoids edge contributions, we uncover strikingly large quantum oscillations in the monolayer insulator's magnetoresistance, with an onset field as small as ~ 0.5 tesla. Despite the huge resistance, the oscillation profile, which exhibits many periods, mimics the Shubnikov-de Haas oscillations in metals. Remarkably, at ultralow temperatures the observed oscillations evolve into discrete peaks near 1.6 tesla, above which the Landau quantized regime is fully developed. Such a low onset field of quantization is comparable to high-mobility conventional two-dimensional electron gases. Our experiments call for further investigation of the highly unusual ground state of the WTe$_{2}$ monolayer. This includes the influence of device components and the possible existence of mobile fermions and charge-neutral Fermi surfaces inside its insulating gap.


title: A BERT-based Distractor Generation Scheme with Multi-tasking and Negative Answer Training Strategies

abstract: In this paper, we investigate the following two limitations for the existing distractor generation (DG) methods. First, the quality of the existing DG methods are still far from practical use. There is still room for DG quality improvement. Second, the existing DG designs are mainly for single distractor generation. However, for practical MCQ preparation, multiple distractors are desired. Aiming at these goals, in this paper, we present a new distractor generation scheme with multi-tasking and negative answer training strategies for effectively generating \textit{multiple} distractors. The experimental results show that (1) our model advances the state-of-the-art result from 28.65 to 39.81 (BLEU 1 score) and (2) the generated multiple distractors are diverse and show strong distracting power for multiple choice question.


title: Non-Einstein relative Yamabe metrics

abstract: In this paper, we give a sufficient condition for a positive constant scalar curvature metric on a manifold with boundary to be a relative Yamabe metric, which is a natural relative version of the classical Yamabe metric. We also give examples of non-Einstein relative Yamabe metrics with positive scalar curvature.


title: Multi-Objective Bayesian Optimisation and Joint Inversion for Active Sensor Fusion

abstract: A critical decision process in data acquisition for mineral and energy resource exploration is how to efficiently combine a variety of sensor types and to minimize total cost. We propose a probabilistic framework for multi-objective optimisation and inverse problems given an expensive cost function for allocating new measurements. This new method is devised to jointly solve multi-linear forward models of 2D-sensor data and 3D-geophysical properties using sparse Gaussian Process kernels while taking into account the cross-variances of different parameters. Multiple optimisation strategies are tested and evaluated on a set of synthetic and real geophysical data. We demonstrate the advantages on a specific example of a joint inverse problem, recommending where to place new drill-core measurements given 2D gravity and magnetic sensor data, the same approach can be applied to a variety of remote sensing problems with linear forward models - ranging from constraints limiting surface access for data acquisition to adaptive multi-sensor positioning.


title: Meta-Active Learning for Node Response Prediction in Graphs

abstract: Meta-learning is an important approach to improve machine learning performance with a limited number of observations for target tasks. However, when observations are unbalancedly obtained, it is difficult to improve the performance even with meta-learning methods. In this paper, we propose an active learning method for meta-learning on node response prediction tasks in attributed graphs, where nodes to observe are selected to improve performance with as few observed nodes as possible. With the proposed method, we use models based on graph convolutional neural networks for both predicting node responses and selecting nodes, by which we can predict responses and select nodes even for graphs with unseen response variables. The response prediction model is trained by minimizing the expected test error. The node selection model is trained by maximizing the expected error reduction with reinforcement learning. We demonstrate the effectiveness of the proposed method with 11 types of road congestion prediction tasks.


title: AI Song Contest: Human-AI Co-Creation in Songwriting

abstract: Machine learning is challenging the way we make music. Although research in deep generative models has dramatically improved the capability and fluency of music models, recent work has shown that it can be challenging for humans to partner with this new class of algorithms. In this paper, we present findings on what 13 musician/developer teams, a total of 61 users, needed when co-creating a song with AI, the challenges they faced, and how they leveraged and repurposed existing characteristics of AI to overcome some of these challenges. Many teams adopted modular approaches, such as independently running multiple smaller models that align with the musical building blocks of a song, before re-combining their results. As ML models are not easily steerable, teams also generated massive numbers of samples and curated them post-hoc, or used a range of strategies to direct the generation, or algorithmically ranked the samples. Ultimately, teams not only had to manage the "flare and focus" aspects of the creative process, but also juggle them with a parallel process of exploring and curating multiple ML models and outputs. These findings reflect a need to design machine learning-powered music interfaces that are more decomposable, steerable, interpretable, and adaptive, which in return will enable artists to more effectively explore how AI can extend their personal expression.


title: A Discussion on Stabilization of Frequency Control for Power Systems

abstract: How to practically maintain the frequency stability of large-scale power systems by a decentralized way is a simple but non-trivial question. In other words, is it possible to design any local controller without understanding the other controlled areas and with less understanding of network structure? With respect to the special properties of physical interaction between the local areas, this paper suggests two existing theories for tackling this issue. Firstly, passivity theory is shown to be a candidate for frequency control problem using swing equation. Based on the passivity of swing dynamics, it is possible to guarantee the system stability by designing for each local area a passive controller. We further extend the passivity approach to the hierarchically decentralized control system with unknown communication delay. Secondly, we discuss the application of generalized frequency variable (GFV) to the frequency control problem using area-control-error. Each local controller is designed such that each local subsystem follows a nominal model set. Utilizing GFV theory, we present a triad of conditions that sufficiently guarantee the system stability. The conditions can be tested conveniently by a limited set of inequalities established from the GFV and the eigenvalues of the physical interaction matrix. The effectiveness, limitation, and challenge of two theories are discussed by design examples with numerical simulations.


title: Evidence for a Monolayer Excitonic Insulator

abstract: The interplay between topology and correlations can generate a variety of quantum phases, many of which remain to be explored. Recent advances have identified monolayer WTe2 as a promising material for doing so in a highly tunable fashion. The ground state of this two-dimensional (2D) crystal can be electrostatically tuned from a quantum spin Hall insulator (QSHI) to a superconductor. However, much remains unknown about the gap-opening mechanism of the insulating state. Here we report evidence that the QSHI is also an excitonic insulator (EI), arising from the spontaneous formation of electron-hole bound states (excitons). We reveal the presence of an intrinsic insulating state at the charge neutrality point (CNP) in clean samples and confirm the correlated nature of this charge-neutral insulator by tunneling spectroscopy. We provide evidence against alternative scenarios of a band insulator or a localized insulator and support the existence of an EI phase in the clean limit. These observations lay the foundation for understanding a new class of correlated insulators with nontrivial topology and identify monolayer WTe2 as a promising candidate for exploring quantum phases of ground-state excitons.


title: A Progressive Conditional Generative Adversarial Network for Generating Dense and Colored 3D Point Clouds

abstract: In this paper, we introduce a novel conditional generative adversarial network that creates dense 3D point clouds, with color, for assorted classes of objects in an unsupervised manner. To overcome the difficulty of capturing intricate details at high resolutions, we propose a point transformer that progressively grows the network through the use of graph convolutions. The network is composed of a leaf output layer and an initial set of branches. Every training iteration evolves a point vector into a point cloud of increasing resolution. After a fixed number of iterations, the number of branches is increased by replicating the last branch. Experimental results show that our network is capable of learning and mimicking a 3D data distribution, and produces colored point clouds with fine details at multiple resolutions.


title: Magnetically regulated disk-formation in the inner 100 au region of the Class 0 young stellar object OMC-3/MMS~6 resolved by JVLA and ALMA

abstract: We have carried out polarization calibration for archival JVLA ($\sim$9 mm) full polarization observations towards the Class 0 young stellar object (YSO) OMC-3/MMS 6 (also known as HOPS-87), and then compared with the archival ALMA 1.2 mm observations. We found that the innermost $\sim$100 au region of OMC-3/MMS 6 is likely very optically thick (e.g., $\tau\gg$1) at $\sim$1 mm wavelength such that the dominant polarization mechanism is dichroic extinction. It is marginally optically thin (e.g., $\tau\lesssim$1) at $\sim$9 mm wavelength such that the JVLA observations can directly probe the linearly polarized emission from non-spherical dust. Assuming that the projected long axis of dust grains is aligned perpendicular to magnetic field (B-field) lines, we propose that the overall B-field topology resembles an hourglass shape, while this "hourglass" appears $\sim$40$^{\circ}$ inclined with respect to the previously reported outflow axis. The geometry of this system is consistent with a magnetically regulated dense (pseudo-)disk. Based on the observed 29.45 GHz flux density and assuming a dust absorption opacity $\kappa^{abs}_{29.45\,GHz}=$0.0096 cm$^{2} $g$^{-1}$, the derived overall dust mass within a $\sim$43 au radius is $\sim$14000 $M_{\oplus}$. From this case study, it appears to us that some previous 9 mm surveys towards Class 0/I YSOs might have systematically underestimated dust masses by one order of magnitude, owing to that they assumed the too high dust absorption opacity ($\sim$0.1 cm$^{2}$ g$^{-1}$) for $\sim$9 mm wavelengths but without self-consistently considering the dust scattering opacity.


title: Experimental demonstration of novel optical Brownian ratchet by controllable phase profile of light

abstract: Brownian ratchet has emerged as a promising tool for understanding motion mechanism of molecules and proteins, and dynamically manipulating particles in non-equilibrium thermodynamics state. Here, we propose and experimentally demonstrate a new type of optical Brownian ratchets, which is generated by using controllable phase profiles in holographic optical trapping system. The potential energy profiles are dynamically switched in on-off mode, controlled by Labview software. Experimental results reveal that not only high speed but also large step distance can be easily achieved in Brownian ratchet, in which the average velocity of forward motion is about 28 {\mu}m/s and step distance is about 42 {\mu}m. More importantly, its motion path, step distance and velocity can be easily changed by manipulating its holograms, which exhibits high flexibility, easy control and excellent capability. This technique provides new opportunities for exploring non-equilibrium dynamics at the nanoscale level, and developing novel nanoparticle manipulation.


title: Exploiting Local Optimality in Metaheuristic Search

abstract: A variety of strategies have been proposed for overcoming local optimality in metaheuristic search. This paper examines characteristics of moves that can be exploited to make good decisions about steps that lead away from a local optimum and then lead toward a new local optimum. We introduce strategies to identify and take advantage of useful features of solution history with an adaptive memory metaheuristic, to provide rules for selecting moves that offer promise for discovering improved local optima. Our approach uses a new type of adaptive memory based on a construction called exponential extrapolation. The memory operates by means of threshold inequalities that ensure selected moves will not lead to a specified number of most recently encountered local optima. Associated thresholds are embodied in choice rule strategies that further exploit the exponential extrapolation concept. Together these produce a threshold based Alternating Ascent (AA) algorithm that opens a variety of research possibilities for exploration.


title: Reconstruction of Quantitative Susceptibility Maps from Phase of Susceptibility Weighted Imaging with Cross-Connected $\Psi$-Net

abstract: Quantitative Susceptibility Mapping (QSM) is a new phase-based technique for quantifying magnetic susceptibility. The existing QSM reconstruction methods generally require complicated pre-processing on high-quality phase data. In this work, we propose to explore a new value of the high-pass filtered phase data generated in susceptibility weighted imaging (SWI), and develop an end-to-end Cross-connected $\Psi$-Net (C$\Psi$-Net) to reconstruct QSM directly from these phase data in SWI without additional pre-processing. C$\Psi$-Net adds an intermediate branch in the classical U-Net to form a $\Psi$-like structure. The specially designed dilated interaction block is embedded in each level of this branch to enlarge the receptive fields for capturing more susceptibility information from a wider spatial range of phase images. Moreover, the crossed connections are utilized between branches to implement a multi-resolution feature fusion scheme, which helps C$\Psi$-Net capture rich contextual information for accurate reconstruction. The experimental results on a human dataset show that C$\Psi$-Net achieves superior performance in our task over other QSM reconstruction algorithms.


title: New P$c$N and AP$c$N functions over finite fields

abstract: Functions with low $c$-differential uniformity were proposed in $2020$ and attracted lots of attention, especially the P$c$N and AP$c$N functions, due to their applications in cryptography. The objective of this paper is to study P$c$N and AP$c$N functions. As a consequence, we propose a class of P$c$N functions and four classes of AP$c$N functions by using the cyclotomic technique and the switch method. In addition, four classes of P$c$N or AP$c$N functions are presented by virtue of (generalized) AGW criterion.


title: RNN Training along Locally Optimal Trajectories via Frank-Wolfe Algorithm

abstract: We propose a novel and efficient training method for RNNs by iteratively seeking a local minima on the loss surface within a small region, and leverage this directional vector for the update, in an outer-loop. We propose to utilize the Frank-Wolfe (FW) algorithm in this context. Although, FW implicitly involves normalized gradients, which can lead to a slow convergence rate, we develop a novel RNN training method that, surprisingly, even with the additional cost, the overall training cost is empirically observed to be lower than back-propagation. Our method leads to a new Frank-Wolfe method, that is in essence an SGD algorithm with a restart scheme. We prove that under certain conditions our algorithm has a sublinear convergence rate of $O(1/\epsilon)$ for $\epsilon$ error. We then conduct empirical experiments on several benchmark datasets including those that exhibit long-term dependencies, and show significant performance improvement. We also experiment with deep RNN architectures and show efficient training performance. Finally, we demonstrate that our training method is robust to noisy data.


title: Tight Bounds for a Class of Data-Driven Distributionally Robust Risk Measures

abstract: This paper expands the notion of robust moment problems to incorporate distributional ambiguity using Wasserstein distance as the ambiguity measure. The classical Chebyshev-Cantelli (zeroth partial moment) inequalities, Scarf and Lo (first partial moment) bounds, and semideviation (second partial moment) in one dimension are investigated. The infinite dimensional primal problems are formulated and the simpler finite dimensional dual problems are derived. A principal motivating question is how does data-driven distributional ambiguity affect the moment bounds. Towards answering this question, some theory is developed and computational experiments are conducted for specific problem instances in inventory control and portfolio management. Finally some open questions and suggestions for future research are discussed.


title: Beating the House: Fast Simulation of Dissipative Quantum Systems with Ensemble Rank Truncation

abstract: We introduce a new technique for the simulation of dissipative quantum systems. This method is composed of an approximate decomposition of the Lindblad equation into a Kraus map, from which one can define an ensemble of wavefunctions. Using principal component analysis, this ensemble can be truncated to a manageable size without sacrificing numerical accuracy. We term this method \emph{Ensemble Rank Truncation} (ERT), and find that in the regime of weak coupling, this method is able to outperform existing wavefunction Monte-Carlo methods by an order of magnitude in both accuracy and speed. We also explore the possibility of combining ERT with approximate techniques for simulating large systems (such as Matrix Product States (MPS)), and show that in many cases this approach will be more efficient than directly expressing the density matrix in its MPS form. We expect the ERT technique to be of practical interest when simulating dissipative systems for quantum information, metrology and thermodynamics.


title: Efficient evaluation of two-center Gaussian integrals in periodic systems

abstract: By using Poisson's summation formula, we calculate periodic integrals over Gaussian basis functions by partitioning the lattice summations between the real and reciprocal space, where both sums converge exponentially fast with a large exponent. We demonstrate that the summation can be performed efficiently to calculate 2-center Gaussian integrals over various kernels including overlap, kinetic, and Coulomb. The summation in real space is performed using an efficient flavor of the McMurchie-Davidson Recurrence Relation (MDRR). The expressions for performing summation in the reciprocal space are also derived and implemented. The algorithm for reciprocal space summation allows us to reuse several terms and leads to significant improvement in efficiency when highly contracted basis functions with large exponents are used. We find that the resulting algorithm is only between a factor of 5 to 15 slower than that for molecular integrals, indicating the very small number of terms needed in both the real and reciprocal space summations. An outline of the algorithm for calculating 3-center Coulomb integrals is also provided.


title: Complete solutions of Toda equations and cyclic Higgs bundles over non-compact surfaces

abstract: On a Riemann surface with a holomorphic $r$-differential, one can naturally define a Toda equation and a cyclic Higgs bundle with a grading. A solution of the Toda equation is equivalent to a harmonic metric of the Higgs bundle for which the grading is orthogonal. Here we focus on a general non-compact Riemann surface with an $r$-differential which is not necessarily meromorphic at infinity. We introduce the notion of complete solution of the Toda equation, and we prove the existence and uniqueness of a complete solution by using techniques for both Toda equations and harmonic bundles. Moreover, we show some quantitative estimates of the complete solution.


title: Quantum Monte Carlo study of topological phases on a spin analogue of Benalcazar-Bernevig-Hughes model

abstract: We study the higher-order topological spin phases based on a spin analogue of Benalcazar-Bernevig-Hughes model in two dimensions using large-scale quantum Monte Carlo simulations. A continuous N\'eel-valence bond solid quantum phase transition is revealed by tuning the ratio between dimerized spin couplings, namely, the weak and strong exchange couplings. Through the finite-size scaling analysis, we identify the phase critical points, and consequently, map out the full phase diagrams in related parameter spaces. Particularly, we find that the valence bond solid phase can be a higher-order topological spin phase, which has a gap for spin excitations in the bulk while demonstrates characteristic gapless spin modes at corners of open lattices. We further discuss the connection between the higher-order topological spin phases and the electronic correlated higher-order phases, and find both of them possess gapless spin corner modes that are protected by higher-order topology. Our result exemplifies higher-order physics in the correlated spin systems and will contribute to further understandings of the many-body higher-order topological phenomena.


title: Distribution of Si II $\lambda$6355 Velocities of Type Ia Supernovae and Implications for Asymmetric Explosions

abstract: The ejecta velocity is a very important parameter in studying the structure and properties of Type Ia supernovae (SNe Ia). It is also a candidate key parameter in improving the utility of SNe Ia for cosmological distance determinations. Here we study the velocity distribution of a sample of 311 SNe Ia from the kaepora database. The velocities are derived from the Si II $\lambda$6355 absorption line in optical spectra measured at (or extrapolated to) the time of peak brightness. We statistically show that the observed velocity has a bimodal Gaussian distribution consisting of two groups of SNe Ia: Group I with a lower but narrower scatter ($\mu_1 = 11000 \text{km s}^{-1}$, $\sigma_1 = 700 \text{km s}^{-1}$), and Group II with a higher but broader scatter ($\mu_2 = 12300 \text{km s}^{-1}$, $\sigma_2 = 1800 \text{km s}^{-1}$). The population ratio of Group I to Group II is 201:110 (65%:35%). There is substantial degeneracy between the two groups, but for SNe Ia with velocity $v > 12000 \text{km s}^{-1}$, the distribution is dominated by Group II. The true origin of the two components is unknown, though there could be that naturally there exist two intrinsic velocity distributions as observed. However, we try to use asymmetric geometric models through statistical simulations to reproduce the observed distribution assuming all SNe Ia share the same intrinsic distribution. In the two cases we consider, 35\% of SNe Ia are considered to be asymmetric in Case 1, and all SNe Ia are asymmetric in Case 2. Simulations for both cases can reproduce the observed velocity distribution but require a significantly large portion ($>35\%$) of SNe Ia to be asymmetric. In addition, the Case 1 result is consistent with recent polarization observations that SNe Ia with higher Si II $\lambda$6355 velocity tend to be more polarized.


title: Locality Preserving Dense Graph Convolutional Networks with Graph Context-Aware Node Representations

abstract: Graph convolutional networks (GCNs) have been widely used for representation learning on graph data, which can capture structural patterns on a graph via specifically designed convolution and readout operations. In many graph classification applications, GCN-based approaches have outperformed traditional methods. However, most of the existing GCNs are inefficient to preserve local information of graphs -- a limitation that is especially problematic for graph classification. In this work, we propose a locality-preserving dense GCN with graph context-aware node representations. Specifically, our proposed model incorporates a local node feature reconstruction module to preserve initial node features into node representations, which is realized via a simple but effective encoder-decoder mechanism. To capture local structural patterns in neighbourhoods representing different ranges of locality, dense connectivity is introduced to connect each convolutional layer and its corresponding readout with all previous convolutional layers. To enhance node representativeness, the output of each convolutional layer is concatenated with the output of the previous layer's readout to form a global context-aware node representation. In addition, a self-attention module is introduced to aggregate layer-wise representations to form the final representation. Experiments on benchmark datasets demonstrate the superiority of the proposed model over state-of-the-art methods in terms of classification accuracy.


title: Dynamic Multiple Nanoparticle Trapping using Metamaterial Plasmonic Tweezers

abstract: Optical manipulation has attracted remarkable interest owing to its versatile and non-invasive nature. However, conventional optical trapping remains inefficient for the nanoscopic world. The emergence of plasmonics in recent years has brought a revolutionary change in overcoming limitations due to diffraction and the requirements for high trapping laser powers. Among the near-field optical trapping cavity-based systems, Fano resonant optical tweezers have a robust trapping capability. In this work, we experimentally demonstrate sequential trapping of 20 nm particles through the use of metamaterial plasmonic optical tweezers. We investigate the multiple trapping via trap stiffness measurements for various trapping positions at low and high incident laser intensities. Our configuration could be used as a light-driven nanoscale sorting device under low laser excitation. Our results provide an alternative approach to trap multiple nanoparticles at distinct hotspots, enabling new ways to control mass transport on the nanoscale.


title: VMSMO: Learning to Generate Multimodal Summary for Video-based News Articles

abstract: A popular multimedia news format nowadays is providing users with a lively video and a corresponding news article, which is employed by influential news media including CNN, BBC, and social media including Twitter and Weibo. In such a case, automatically choosing a proper cover frame of the video and generating an appropriate textual summary of the article can help editors save time, and readers make the decision more effectively. Hence, in this paper, we propose the task of Video-based Multimodal Summarization with Multimodal Output (VMSMO) to tackle such a problem. The main challenge in this task is to jointly model the temporal dependency of video with semantic meaning of article. To this end, we propose a Dual-Interaction-based Multimodal Summarizer (DIMS), consisting of a dual interaction module and multimodal generator. In the dual interaction module, we propose a conditional self-attention mechanism that captures local semantic information within video and a global-attention mechanism that handles the semantic relationship between news text and video from a high level. Extensive experiments conducted on a large-scale real-world VMSMO dataset show that DIMS achieves the state-of-the-art performance in terms of both automatic metrics and human evaluations.


title: The Cayley graphs associated with some quasi-perfect Lee codes are Ramanujan graphs

abstract: Let $\Z_n[i]$ be the ring of Gaussian integers modulo a positive integer $n$. Very recently, Camarero and Mart\'{i}nez [IEEE Trans. Inform. Theory, {\bf 62} (2016), 1183--1192], showed that for every prime number $p>5$ such that $p\equiv \pm 5 \pmod{12}$, the Cayley graph $\mathcal{G}_p=\textnormal{Cay}(\Z_p[i], S_2)$, where $S_2$ is the set of units of $\Z_p[i]$, induces a 2-quasi-perfect Lee code over $\Z_p^m$, where $m=2\lfloor \frac{p}{4}\rfloor$. They also conjectured that $\mathcal{G}_p$ is a Ramanujan graph for every prime $p$ such that $p\equiv 3 \pmod{4}$. In this paper, we solve this conjecture. Our main tools are Deligne's bound from 1977 for estimating a particular kind of trigonometric sum and a result of Lov\'{a}sz from 1975 (or of Babai from 1979) which gives the eigenvalues of Cayley graphs of finite Abelian groups. Our proof techniques may motivate more work in the interactions between spectral graph theory, character theory, and coding theory, and may provide new ideas towards the famous Golomb--Welch conjecture on the existence of perfect Lee codes.


title: Gauge-induced Floquet topological states in photonic waveguides

abstract: Tremendous efforts have been devoted to the search for exotic topological states, which usually exist at an interface between lattices with differing topological invariants according to the bulk-edge correspondence. Here, we show a new finding of topological states localized at the interface between two gauge-shifted Floquet photonic lattices, despite the same topological order across the entire structure. The quasienergy band structures reveal that these new interface modes belong to the Floquet {\pi} modes, which are further found to enable a robust one-way propagation thanks to the flexible control of the Floquet gauge. The intriguing propagations of these {\pi}-interface modes are experimentally verified in a silicon waveguides platform at near-infrared wavelengths, which show both broad working bandwidth and high tolerance to the structural fluctuations. Our approach provides a new route for light manipulations with robust behaviours in Floquet engineering and beyond.


title: Band Structure for Thermal Conduction in Multilayered Systems. Thermal Crystals

abstract: In this paper we solve the Cattaneo-Vernotte Equation for a periodic heterostructure made of alternate layers of different materials. The solutions describe thermal waves traveling in a periodic system, and it allows us to introduce the concept of thermal crystals. We show that the dispersion relation shows the characteristics of a band-structure, however the corresponding Bloch wave vector is always complex corresponding to pseudo-bands, unlike what happens in photonic or acoustic crystals. In this context, we also discuss the use of the Floquet-Bloch theorem for thermal waves. The case of finite layered structures is also analyzed showing the possibility of changing the temperature and heat flux by introducing defects opening the possibility of thermal management through the pseudo-band structure.


title: Topological Hall effect and the magnetic states of Nowotney chimney ladder compound Cr$_{11}$Ge$_{19}$}

abstract: We have investigated the magnetic and charge transport properties of single crystals of Nowotney Chimney Ladder compound Cr$_{11}$Ge$_{19}$ and mapped out a comprehensive phase diagram reflecting the complicated interplay between the Dzyaloshinskii-Moriya (DM) interaction, the dipolar interaction, and the magnetic anisotropy. We have identified a set of interesting magnetic phases and attributed a finite topological Hall effect to the recently discovered bi-skyrmion phase. These data also suggest the existence of an anti-skyrmion state at finite fields for temperatures just below the magnetic ordering temperature, $T_c$, as indicated by a distinct change in sign of the topological Hall effect. Above $T_c$, we discovered a region of enhanced magnetic response corresponding to a disordered phase likely existing near the ferromagnetic critical point under small magnetic fields. Strong spin chirality fluctuations are demonstrated by the large value of the topological Hall resistivity persisting up to 1 T which is most likely due to the existence of the DM interaction. We argue that changes to the topological Hall effect correspond to different topological spin textures that are controlled by magnetic dipolar and DM interactions that vary in importance with temperature.


title: Deep Learning Prediction of Adverse Drug Reactions Using Open TG-GATEs and FAERS Databases

abstract: With the advancements in Artificial intelligence (AI) and the accumulation of healthrelated big data, it has become increasingly feasible and commonplace to leverage machine learning technologies to analyze clinical and omics metadata to assess the possibility of adverse drug reactions or events (ADRs) in the course of drug discovery. Here, we have described a novel approach that combined drug-induced gene expression profile from Open TG-GATEs (Toxicogenomics Project-Genomics Assisted Toxicity Evaluation Systems) and ADR occurrence information from FAERS (FDA [Food and Drug Administration] Adverse Events Reporting System) database to predict the likelihood of ADRs. We generated a total of 14 models using Deep Neural Networks (DNN) to predict different ADRs; in the validation tests, our models achieved a mean accuracy of 85.71%, indicating that our approach successfully and consistently predicted ADRs for a wide range of drugs. As an example, we have described the ADR model in the context of Duodenal ulcer. We believe that our models will help predict the likelihood of ADRs while testing novel pharmaceutical compounds, and will be useful for researchers in drug discovery.


title: On the origin of X-ray oxygen emission-lines in obscured AGN

abstract: We present the Catalog of High REsolution Spectra of Obscured Sources (CHRESOS) from the XMM-Newton Science Archive. It comprises the emission-line luminosities of H- and He-like transitions from C to Si, and the Fe 3C and Fe 3G L-shell ones. Here, we concentrate on the soft X-ray OVII(f) and OVIII Ly_alpha emission lines to shed light onto the physical processes with which their formation can be related to: active galactic nucleus vs. star forming regions. We compare their luminosity with that of two other important oxygen key lines [OIII]5007A, in the optical, and [OIV]25.89mic, in the IR. We also test OVII(f) and OVIIILy_alpha luminosities against that of continuum bands in the IR and hard X-rays, which point to different ionization processes. We probe into those processes by analyzing photoionization and colisional ionization model predictions upon our lines. We show that both scenarios can explain the formation and observed intensities of OVII(f) and OVIII Ly_alpha. By analyzing the relationships between OVII(f) and OVIII Ly_alpha, and all other observables: [OIII]5007A, [OIV]25.89mic emission lines, and MIR-12mic, FIR-60mic, FIR-100mic, 2-10 keV and 14-195 keV continuum bands, we conclude that the AGN radiation field is mainly responsible of the soft X-ray oxygen excitation.


title: Dark supernova remnant

abstract: An almost perfect round hole of CO-line emission with a diameter of 3.7 pc was found in a molecular cloud (MC) centered on G35.75-0.25 ($l = 35^\circ.75, b = -0^\circ.25$) at radial velocity of 28 km s$^{-1}$. The hole is quiet in radio continuum emission, unlike the usual supernova remnants (SNR), and the molecular edge is only weakly visible in 8 and 24 $\mu$m dust emissions. The hole may be either a fully evolved molecular bubble around a young stellar object (YSO), or a relic of a radio-quiet SNR that has already stopped expansion after rapid evolution in the dense MC as a buried SNR. Because G35.75 exhibits quite different properties from YSO-driven bubbles of the same size, we prefer the latter interpretation. Existence of such a "dark" SNR would affect the estimation of the supernova rate, and therefore the star formation history in the Galaxy.


title: Heat kernel upper bounds for symmetric Markov semigroups

abstract: It is well known that Nash-type inequalities for symmetric Dirichlet forms are equivalent to on-diagonal heat kernel upper bounds for the associated symmetric Markov semigroups. In this paper, we show that both imply (and hence are equivalent to) off-diagonal heat kernel upper bounds under some mild assumptions. Our approach is based on a new generalized Davies' method. Our results extend that of \cite{CKS} for Nash-type inequalities with power order considerably and also extend that of \cite{Gri} for second order differential operators on a complete non-compact manifold.


title: Unweighted linear congruences with distinct coordinates and the Varshamov--Tenengolts codes

abstract: In this paper, we first give explicit formulas for the number of solutions of unweighted linear congruences with distinct coordinates. Our main tools are properties of Ramanujan sums and of the discrete Fourier transform of arithmetic functions. Then, as an application, we derive an explicit formula for the number of codewords in the Varshamov--Tenengolts code $VT_b(n)$ with Hamming weight $k$, that is, with exactly $k$ $1$'s. The Varshamov--Tenengolts codes are an important class of codes that are capable of correcting asymmetric errors on a $Z$-channel. As another application, we derive Ginzburg's formula for the number of codewords in $VT_b(n)$, that is, $|VT_b(n)|$. We even go further and discuss connections to several other combinatorial problems, some of which have appeared in seemingly unrelated contexts. This provides a general framework and gives new insight into all these problems which might lead to further work.


title: Rhythmic Control of Automated Traffic -- Part II: Grid Network Rhythm and Online Routing

abstract: Connected and automated vehicle (CAV) technology is providing urban transportation managers tremendous opportunities for better operation of urban mobility systems. However, there are significant challenges in real-time implementation, as the computational time of the corresponding operations optimization model increases exponentially with increasing vehicle numbers. Following the companion paper (Chen et al., 2020) which proposes a novel automated traffic control scheme for isolated intersections, this study proposes a network-level real-time traffic control framework for CAVs on grid networks. The proposed framework integrates a rhythmic control (RC) method with an online routing algorithm to realize collisionfree control of all CAVs on a network and achieve superior performance in average vehicle delay, network traffic throughput, and computational scalability. Specifically, we construct a preset network rhythm that all CAVs can follow to move on the network and avoid collisions at all intersections. Based on the network rhythm, we then formulate online routing for the CAVs as a mixed integer linear program, which optimizes the entry times of CAVs at all entrances of the network and their time-space routings in real time. We provide a sufficient condition that the linear programming relaxation of the online routing model yields an optimal integer solution. Extensive numerical tests are conducted to show the performance of the proposed operations management framework under various scenarios. It is illustrated that the framework is capable of achieving negligible delays and increased network throughput. Furthermore, the computational time results are also promising. The CPU time for solving a collision-free control optimization problem with 2,000 vehicles is only 0.3 s on an ordinary personal computer.


title: Semiclassical asymptotics for a class of singular Schr\"odinger operators

abstract: Let $\Omega \subset \mathbb{R}^d$ be bounded with $C^1$ boundary. In this paper we consider Schr\"odinger operators $-\Delta+ W$ on $\Omega$ with $W(x)\approx\mathrm{dist}(x, \partial\Omega)^{-2}$ as $\mathrm{dist}(x, \partial\Omega)\to 0$. Under weak assumptions on $W$ we derive a two-term asymptotic formula for the sum of the eigenvalues of such operators.


title: Achilles Heels for AGI/ASI via Decision Theoretic Adversaries

abstract: As progress in AI continues to advance, it is important to know how advanced systems will make choices and in what ways they may fail. Machines can already outsmart humans in some domains, and understanding how to safely build ones which may have capabilities at or above the human level is of particular concern. One might suspect that artificially generally intelligent (AGI) and artificially superintelligent (ASI) will be systems that humans cannot reliably outsmart. As a challenge to this assumption, this paper presents the Achilles Heel hypothesis which states that even a potentially superintelligent system may nonetheless have stable decision-theoretic delusions which cause them to make irrational decisions in adversarial settings. In a survey of key dilemmas and paradoxes from the decision theory literature, a number of these potential Achilles Heels are discussed in context of this hypothesis. Several novel contributions are made toward understanding the ways in which these weaknesses might be implanted into a system.


title: Gradient-based Analysis of NLP Models is Manipulable

abstract: Gradient-based analysis methods, such as saliency map visualizations and adversarial input perturbations, have found widespread use in interpreting neural NLP models due to their simplicity, flexibility, and most importantly, their faithfulness. In this paper, however, we demonstrate that the gradients of a model are easily manipulable, and thus bring into question the reliability of gradient-based analyses. In particular, we merge the layers of a target model with a Facade that overwhelms the gradients without affecting the predictions. This Facade can be trained to have gradients that are misleading and irrelevant to the task, such as focusing only on the stop words in the input. On a variety of NLP tasks (text classification, NLI, and QA), we show that our method can manipulate numerous gradient-based analysis techniques: saliency maps, input reduction, and adversarial perturbations all identify unimportant or targeted tokens as being highly important. The code and a tutorial of this paper is available at http://ucinlp.github.io/facade.


title: MMH* with arbitrary modulus is always almost-universal

abstract: Universal hash functions, discovered by Carter and Wegman in 1979, are of great importance in computer science with many applications. MMH$^*$ is a well-known $\triangle$-universal hash function family, based on the evaluation of a dot product modulo a prime. In this paper, we introduce a generalization of MMH$^*$, that we call GMMH$^*$, using the same construction as MMH$^*$ but with an arbitrary integer modulus $n>1$, and show that GMMH$^*$ is $\frac{1}{p}$-almost-$\triangle$-universal, where $p$ is the smallest prime divisor of $n$. This bound is tight.


title: Factorizable Graph Convolutional Networks

abstract: Graphs have been widely adopted to denote structural connections between entities. The relations are in many cases heterogeneous, but entangled together and denoted merely as a single edge between a pair of nodes. For example, in a social network graph, users in different latent relationships like friends and colleagues, are usually connected via a bare edge that conceals such intrinsic connections. In this paper, we introduce a novel graph convolutional network (GCN), termed as factorizable graph convolutional network(FactorGCN), that explicitly disentangles such intertwined relations encoded in a graph. FactorGCN takes a simple graph as input, and disentangles it into several factorized graphs, each of which represents a latent and disentangled relation among nodes. The features of the nodes are then aggregated separately in each factorized latent space to produce disentangled features, which further leads to better performances for downstream tasks. We evaluate the proposed FactorGCN both qualitatively and quantitatively on the synthetic and real-world datasets, and demonstrate that it yields truly encouraging results in terms of both disentangling and feature aggregation. Code is publicly available at https://github.com/ihollywhy/FactorGCN.PyTorch.


title: Global HI asymmetries in IllustrisTNG: a diversity of physical processes disturb the cold gas in galaxies

abstract: Observations of the cold neutral atomic hydrogen (HI) in and around disc galaxies have revealed that spatial and kinematic asymmetries are commonplace, and are reflected in the global HI spectra. We use the TNG100 box from the IllustrisTNG suite of cosmological simulations to study the conditions under which these asymmetries may arise in current theoretical galaxy formation models. We find that more than 50% of the sample has at least a 10% difference in integrated flux between the high- and low-velocity half of the spectrum, thus the typical TNG100 galaxy has an HI profile that is not fully symmetric. We find that satellite galaxies are a more asymmetric population than centrals, consistent with observational results. Using halo mass as a proxy for environment, this trend appears to be driven by the satellite population within the virial radius of haloes more massive than $10^{13} M_{\odot}$, typical of medium/large groups. We show that, while the excess of HI asymmetry in group satellites is likely driven by ram pressure, the bulk of the asymmetric HI profiles observed in TNG100 are driven by physical processes able to affect both the central and satellite populations. Our results highlight how asymmetries are not driven solely by environment, and multiple physical processes can produce the same asymmetric shape in global HI spectra.


title: Ab initio determination on the thermal evolution of the Earth's core

abstract: Earth's magnetic field is generated by the liquid outer core and sensitively depends on the thermal conductivity of the core. The dominant component of the Earth's core is Fe and Ni. However, current estimates on FeNi mixtures have not been previously tested at high pressures. In this paper, ab initio simulations were first applied to calculations of the thermal and electrical conductivities of FeNi mixtures at Earth's outer core conditions. Compared with the results for pure Fe, the addition of Ni decreases the thermal conductivity (12.30 W/m/K on average) along the adiabatic curve in the outer core. Based on the restriction of the entropy production rate or Joule losses, the existence of Ni prolongs the age of the inner core. The age of the inner core is 0.66 Ga with pure Fe and 0.67 Ga with an FeNi mixture when heat flow at the core-mantle boundary is 12 TW. In contrast, we observe that Ni decreases the thickness of thermal stratification in the outer core by analyzing the effective temperature gradient. After inner core solidification, the thickness of thermal stratification is 417.02 km with pure Fe and 320.12 km with an FeNi mixture when the cooling rate at the core-mantle boundary is 126 K/Ga.


title: Tunable dynamical magnetoelectric effect in antiferromagnetic topological insulator MnBi$_2$Te$_4$ films

abstract: More than forty years ago, axion was postulated as an elementary particle with a low mass and weak interaction in particle physics to solve the strong $\mathcal{CP}$ (charge conjugation and parity) puzzle. Axions are also considered as a possible component of dark matter of the universe. However, the existence of axions in nature has not been confirmed. Interestingly, axions arise as pseudoscalar fields derived from the Chern-Simons theory in condensed matter physics. In antiferromagnetic insulators, the axion field can become dynamical induced by spin-wave excitations and exhibits rich exotic phenomena, such as, the chiral magnetic effect, axionic polariton and so on. However, the study of the dynamical axion field is rare due to the lack of real materials. Recently, MnBi$_2$Te$_4$ was discovered to be an antiferromagnetic topological insulator with a quantized axion field protected by the inversion symmetry $\mathcal{P}$ and the magnetic-crystalline symmetry $\mathcal{S}$. Here, we studied MnBi$_2$Te$_4$ films in which both the $\mathcal{P}$ and $\mathcal{S}$ symmetries are spontaneously broken and found that the dynamical axion field and largely tunable dynamical magnetoelectric effects can be realized through tuning the thickness of MnBi$_2$Te$_4$ films, the temperature and the element substitution. Our results open a broad avenue to study axion dynamics in antiferromagnetic topological insulator MnBi$_2$Te$_4$ and related materials, and also is hopeful to promote the research of dark matter.


title: FPRAS via MCMC where it mixes torpidly (and very little effort)

abstract: Is Fully Polynomial-time Randomized Approximation Scheme (FPRAS) for a problem via an MCMC algorithm possible when it is known that rapid mixing provably fails? We introduce several weight-preserving maps for the eight-vertex model on planar and on bipartite graphs, respectively. Some are one-to-one, while others are holographic which map superpositions of exponentially many states from one setting to another, in a quantum-like many-to-many fashion. In fact we introduce a set of such mappings that forms a group in each case. Using some holographic maps and their compositions we obtain FPRAS for the eight-vertex model at parameter settings where it is known that rapid mixing provably fails due to an intrinsic barrier. This FPRAS is indeed the same MCMC algorithm, except its state space corresponds to superpositions of the given states, where rapid mixing holds. FPRAS is also given for torus graphs for parameter settings where natural Markov chains are known to mix torpidly. Our results show that the eight-vertex model is the first problem with the provable property that while NP-hard to approximate on general graphs (even #P-hard for planar graphs in exact complexity), it possesses FPRAS on both bipartite graphs and planar graphs in substantial regions of its parameter space.


title: Throughput Analysis of Small Cell Networks under D-TDD and FFR

abstract: Dynamic time-division duplex (D-TDD) has emerged as an effective solution to accommodate the unaligned downlink and uplink traffic in small cell networks. However, the flexibility of traffic configuration also introduces additional inter-cell interference. In this letter, we study the effectiveness of applying fractional frequency reuse (FFR) as an interference coordination technique for D-TDD small cell networks. We derive the analytical expressions of downlink and uplink mean packet throughput (MPT), then study a network parameter optimization problem to maximize MPT while guaranteeing each user's throughput. Numerical results corroborate the benefits of the proposed FFR-based D-TDD in terms of improving throughput.