title: HTCC: Haskell to Handel-C Compiler

abstract: Functional programming languages, such as Haskell, enable simple, concise, and correct-by-construction hardware development. HTCC compiles a subset of Haskell to Handel-C language with hardware output. Moreover, HTCC generates VHDL, Verilog, EDIF, and SystemC programs. The design of HTCC compiler includes lexical, syntax and semantic analyzers. HTCC automates a transformational derivation methodology to rapidly produce hardware that maps onto Field Programmable Gate Arrays (FPGAs) . HTCC is generated using ANTLR compiler-compiler tool and supports an effective integrated development environment. This paper presents the design rationale and the implementation of HTCC. Several sample generations of first-class and higher-order functions are presented. In-addition, a compilation case-study is presented for the XTEA cipher. The investigation comprises a thorough evaluation and performance analysis. The targeted FPGAs include Cyclone II, Stratix IV, and Virtex-6 from Altera and Xilinx.


title: Pakistan's Internet Voting Experiment

abstract: Pakistan recently conducted small-scale trials of a remote Internet voting system for overseas citizens. In this contribution, we report on the experience: we document the unique combination of sociopolitical, legal, and institutional factors motivating this exercise. We describe the system and it's reported vulnerabilities, and we also highlight new issues pertaining to materiality. If this system is deployed in the next general elections, as seems likely, this development would constitute the largest enfranchised diaspora in the world. Our goal in this paper, therefore, is to provide comprehensive insight into Pakistan's experiment with Internet voting, emphasize outstanding challenges, and identify directions for future research.


title: Flatter is better: Percentile Transformations for Recommender Systems

abstract: It is well known that explicit user ratings in recommender systems are biased towards high ratings, and that users differ significantly in their usage of the rating scale. Implementers usually compensate for these issues through rating normalization or the inclusion of a user bias term in factorization models. However, these methods adjust only for the central tendency of users' distributions. In this work, we demonstrate that lack of \textit{flatness} in rating distributions is negatively correlated with recommendation performance. We propose a rating transformation model that compensates for skew in the rating distribution as well as its central tendency by converting ratings into percentile values as a pre-processing step before recommendation generation. This transformation flattens the rating distribution, better compensates for differences in rating distributions, and improves recommendation performance. We also show a smoothed version of this transformation designed to yield more intuitive results for users with very narrow rating distributions. A comprehensive set of experiments show improved ranking performance for these percentile transformations with state-of-the-art recommendation algorithms in four real-world data sets.


title: Delta -- new logic programming language and Delta-methodology for p-computable programs on Turing Complete Languages

abstract: In paper describes the new logic programming language Delta, which have a many good properties. Delta-programs is p-computable, verifiable and can translation on other languages. Also we describe the Delta-methodology for constructing p-computable programs in high-level languages such as PHP, Java, JavaScript, C++, Pascal, Delphi, Python, Solidity and other. We would like to especially note the use of the Delta methodology for creating Smart Contracts and for Internet of things. We change the concept of the formula and define D-formulas(or Delta programs) are special list-formulas. Then we define the execution of a program how is the process of checking truth D-formula on a dynamic model. Main idea our paper consider program how list-formula from another formulas on dynamic models. And we created by iterations new Delta-programs use simple base formulas for this. Also we entered a dynamic models how models where we save final values of variables when check formula on this model.


title: A Novel Approach for Detection and Ranking of Trendy and Emerging Cyber Threat Events in Twitter Streams

abstract: We present a new machine learning and text information extraction approach to detection of cyber threat events in Twitter that are novel (previously non-extant) and developing (marked by significance with respect to similarity with a previously detected event). While some existing approaches to event detection measure novelty and trendiness, typically as independent criteria and occasionally as a holistic measure, this work focuses on detecting both novel and developing events using an unsupervised machine learning approach. Furthermore, our proposed approach enables the ranking of cyber threat events based on an importance score by extracting the tweet terms that are characterized as named entities, keywords, or both. We also impute influence to users in order to assign a weighted score to noun phrases in proportion to user influence and the corresponding event scores for named entities and keywords. To evaluate the performance of our proposed approach, we measure the efficiency and detection error rate for events over a specified time interval, relative to human annotator ground truth.


title: Hierarchical Sequence to Sequence Voice Conversion with Limited Data

abstract: We present a voice conversion solution using recurrent sequence to sequence modeling for DNNs. Our solution takes advantage of recent advances in attention based modeling in the fields of Neural Machine Translation (NMT), Text-to-Speech (TTS) and Automatic Speech Recognition (ASR). The problem consists of converting between voices in a parallel setting when {\it $<$source,target$>$} audio pairs are available. Our seq2seq architecture makes use of a hierarchical encoder to summarize input audio frames. On the decoder side, we use an attention based architecture used in recent TTS works. Since there is a dearth of large multispeaker voice conversion databases needed for training DNNs, we resort to training the network with a large single speaker dataset as an autoencoder. This is then adapted for the smaller multispeaker voice conversion datasets available for voice conversion. In contrast with other voice conversion works that use $F_0$, duration and linguistic features, our system uses mel spectrograms as the audio representation. Output mel frames are converted back to audio using a wavenet vocoder.


title: Topology and geometry of molecular conformational spaces and energy landscapes

abstract: Understanding the geometry and topology of configuration or conformational spaces of molecules has relevant applications in chemistry and biology such as the proteins folding problem, drug design and the structure activity relationship problem. Despite their relevance, configuration spaces of molecules are only partially understood. In this paper we discuss both theoretical and computational approaches to the configuration spaces of molecules and their associated energy landscapes. Our mathematical approach shows that when symmetries of the molecules are taken into account, configuration spaces of molecules give rise to certain principal bundles and orbifolds. We also make use of a variety of geometric and topological tools for data analysis to study the topology and geometry of these spaces.


title: Classification Schemas for Artificial Intelligence Failures

abstract: In this paper we examine historical failures of artificial intelligence (AI) and propose a classification scheme for categorizing future failures. By doing so we hope that (a) the responses to future failures can be improved through applying a systematic classification that can be used to simplify the choice of response and (b) future failures can be reduced through augmenting development lifecycles with targeted risk assessments.


title: Modern CNNs for IoT Based Farms

abstract: Recent introduction of ICT in agriculture has brought a number of changes in the way farming is done. This means use of Internet of Things(IoT), Cloud Computing(CC), Big Data (BD) and automation to gain better control over the process of farming. As the use of these technologies in farms has grown exponentially with massive data production, there is need to develop and use state-of-the-art tools in order to gain more insight from the data within reasonable time. In this paper, we present an initial understanding of Convolutional Neural Network (CNN), the recent architectures of state-of-the-art CNN and their underlying complexities. Then we propose a classification taxonomy tailored for agricultural application of CNN. Finally, we present a comprehensive review of research dedicated to applications of state-of-the-art CNNs in agricultural production systems. Our contribution is in two-fold. First, for end users of agricultural deep learning tools, our benchmarking finding can serve as a guide to selecting appropriate architecture to use. Second, for agricultural software developers of deep learning tools, our in-depth analysis explains the state-of-the-art CNN complexities and points out possible future directions to further optimize the running performance.


title: Architectures of Exoplanetary Systems. I: A Clustered Forward Model for Exoplanetary Systems around Kepler's FGK Stars

abstract: Observations of exoplanetary systems provide clues about the intrinsic distribution of planetary systems, their architectures, and how they formed. We develop a forward modelling framework for generating populations of planetary systems and "observed" catalogues by simulating the Kepler detection pipeline (SysSim). We compare our simulated catalogues to the Kepler DR25 catalogue of planet candidates, updated to include revised stellar radii from Gaia DR2. We constrain our model based on the observed 1D marginal distributions of orbital periods, period ratios, transit depths, transit depth ratios, transit durations, transit duration ratios, and transit multiplicities. Models assuming planets with independent periods and sizes do not adequately account for the properties of the multiplanet systems. Instead, a clustered point process model for exoplanet periods and sizes provides a significantly better description of the Kepler population, particularly the observed multiplicity and period ratio distributions. We find that $0.56^{+0.18}_{-0.15}$ of FGK stars have at least one planet larger than $0.5 R_\oplus$ between 3 and 300 d. Most of these planetary systems ($\sim 98\%$) consist of one or two clusters with a median of three planets per cluster. We find that the Kepler dichotomy is evidence for a population of highly inclined planetary systems and is unlikely to be solely due to a population of intrinsically single planet systems. We provide a large ensemble of simulated physical and observed catalogues of planetary systems from our models, as well as publicly available code for generating similar catalogues given user-defined parameters.


title: Image-Driven Biophysical Tumor Growth Model Calibration

abstract: We present a novel formulation for the calibration of a biophysical tumor growth model from a single-time snapshot, MRI scan of a glioblastoma patient. Tumor growth models are typically nonlinear parabolic partial differential equations (PDEs). Thus, we have to generate a second snapshot to be able to extract significant information from a single patient snapshot. We create this two-snapshot scenario as follows. We use an atlas (an average of several scans of healthy individuals) as a substitute for an earlier, pretumor, MRI scan of the patient. Then, using the patient scan and the atlas, we combine image-registration algorithms and parameter estimation algorithms to achieve a better estimate of the healthy patient scan and the tumor growth parameters that are consistent with the data. Our scheme is based on our recent work (Scheufele et al, "Biophysically constrained diffeomorphic image registration, Tumor growth, Atlas registration, Adjoint-based methods, Parallel algorithms", CMAME, 2018), but apply a different and novel scheme where the tumor growth simulation in contrast to the previous work is executed in the patient brain domain and not in the atlas domain yielding more meaningful patient-specific results. As a basis, we use a PDE-constrained optimization framework. We derive a modified Picard-iteration-type solution strategy in which we alternate between registration and tumor parameter estimation in a new way. In addition, we consider an $\ell_1$ sparsity constraint on the initial condition for the tumor and integrate it with the new joint inversion scheme. We solve the subproblems with a reduced-space, inexact Gauss-Newton-Krylov/quasi-Newton methods. We present results using real brain data with synthetic tumor data that show that the new scheme reconstructs the tumor parameters in a more accurate and reliable way compared to our earlier scheme.


title: Model-free Control of Chaos with Continuous Deep Q-learning

abstract: The OGY method is one of control methods for a chaotic system. In the method, we have to calculate a stabilizing periodic orbit embedded in its chaotic attractor. Thus, we cannot use this method in the case where a precise mathematical model of the chaotic system cannot be identified. In this case, the delayed feedback control proposed by Pyragas is useful. However, even in the delayed feedback control, we need the mathematical model to determine a feedback gain that stabilizes the periodic orbit. To overcome this problem, we propose a model-free reinforcement learning algorithm to the design of a controller for the chaotic system. In recent years, model-free reinforcement learning algorithms with deep neural networks have been paid much attention to. Those algorithms make it possible to control complex systems. However, it is known that model-free reinforcement learning algorithms are not efficient because learners must explore their control policies over the entire state space. Moreover, model-free reinforcement learning algorithms with deep neural networks have the disadvantage in taking much time to learn their control optimal policies. Thus, we propose a data-based control policy consisting of two steps, where we determine a region including the stabilizing periodic orbit first, and make the controller learn an optimal control policy for its stabilization. In the proposed method, the controller efficiently explores its control policy only in the region.


title: CADS: Core-Aware Dynamic Scheduler for Multicore Memory Controllers

abstract: Memory controller scheduling is crucial in multicore processors, where DRAM bandwidth is shared. Since increased number of requests from multiple cores of processors becomes a source of bottleneck, scheduling the requests efficiently is necessary to utilize all the computing power these processors offer. However, current multicore processors are using traditional memory controllers, which are designed for single-core processors. They are unable to adapt to changing characteristics of memory workloads that run simultaneously on multiple cores. Existing schedulers may disrupt locality and bank parallelism among data requests coming from different cores. Hence, novel memory controllers that consider and adapt to the memory access characteristics, and share memory resources efficiently and fairly are necessary. We introduce Core-Aware Dynamic Scheduler (CADS) for multicore memory controller. CADS uses Reinforcement Learning (RL) to alter its scheduling strategy dynamically at runtime. Our scheduler utilizes locality among data requests from multiple cores and exploits parallelism in accessing multiple banks of DRAM. CADS is also able to share the DRAM while guaranteeing fairness to all cores accessing memory. Using CADS policy, we achieve 20% better cycles per instruction (CPI) in running memory intensive and compute intensive PARSEC parallel benchmarks simultaneously, and 16% better CPI with SPEC 2006 benchmarks.


title: The Intrinsic Temperature and Radiative-Convective Boundary Depth in the Atmospheres of Hot Jupiters

abstract: In giant planet atmosphere modelling, the intrinsic temperature $T_\mathrm{int}$ and radiative-convective boundary (RCB) are important lower boundary conditions. Often in one-dimensional radiative-convective models and in three-dimensional general circulation models it is assumed that $T_\mathrm{int}$ is similar to that of Jupiter itself, around 100 K, which yields a RCB around 1 kbar for hot Jupiters. In this work, we show that the inflated radii, and hence high specific entropy interiors, of hot Jupiters suggest much higher $T_\mathrm{int}$ values. Assuming the effect is primarily due to current heating (rather than delayed cooling), we derive an equilibrium relation between $T_\mathrm{eq}$ and $T_\mathrm{int}$, showing that the latter can take values as high as 700 K. In response, the RCB moves upward in the atmosphere. Using one-dimensional radiative-convective atmosphere models, we find RCBs of only a few bars, rather than the kilobar typically supposed. This much shallower RCB has important implications for the atmospheric structure, vertical and horizontal circulation, interpretations of phase curves, and the effect of deep cold traps on cloud formation.


title: A Survey on the Use of Preferences for Virtual Machine Placement in Cloud Data Centers

abstract: With the rapid development of virtualization techniques, cloud data centers allow for cost effective, flexible, and customizable deployments of applications on virtualized infrastructure. Virtual machine (VM) placement aims to assign each virtual machine to a server in the cloud environment. VM Placement is of paramount importance to the design of cloud data centers. Typically, VM placement involves complex relations and multiple design factors as well as local policies that govern the assignment decisions. It also involves different constituents including cloud administrators and customers that might have disparate preferences while opting for a placement solution. Thus, it is often valuable to not only return an optimized solution to the VM placement problem but also a solution that reflects the given preferences of the constituents. In this paper, we provide a detailed review on the role of preferences in the recent literature on VM placement. We further discuss key challenges and identify possible research opportunities to better incorporate preferences within the context of VM placement.


title: Symplectic homology of convex domains and Clarke's duality

abstract: We prove that the Floer complex that is associated with a convex Hamiltonian function on $\mathbb{R}^{2n}$ is isomorphic to the Morse complex of Clarke's dual action functional that is associated with the Fenchel-dual Hamiltonian. This isomorphism preserves the action filtrations. As a corollary, we obtain that the symplectic capacity from the symplectic homology of a convex domain with smooth boundary coincides with the minimal action of closed characteristics on its boundary.


title: Persistent atomic frequency comb based on Zeeman sub-levels of an erbium-doped crystal waveguide

abstract: Long-lived sub-levels of the electronic ground-state manifold of rare-earth ions in crystals can be used as atomic population reservoirs for photon echo-based quantum memories. We measure the dynamics of the Zeeman sub-levels of erbium ions that are doped into a lithium niobate waveguide, finding population lifetimes at cryogenic temperatures as long as seconds. Then, using these levels, we prepare and characterize atomic frequency combs, which can serve as a memory for quantum light at 1532 nm wavelength. The results allow predicting a 0.1\% memory efficiency, mainly limited by unwanted background absorption that we conjecture to be caused by the coupling between two-level systems (TLS) and erbium spins. Hence, while it should be possible to create an AFC-based quantum memory in Er$^{3+}$:Ti$^{3+}$:LiNbO$_3$, improved crystal growth together with optimized AFC preparation will be required to make it suitable for applications in quantum communication.


title: Heavy metals in intermediate He-rich hot subdwarfs: The chemical composition of HZ44 and HD127493

abstract: Hot subluminous stars can be spectroscopically classified as subdwarf B (sdB) and O (sdO) stars. While the latter are predominantly hydrogen deficient, the former are mostly helium deficient. The atmospheres of most sdOs are almost devoid of hydrogen, whereas a small group of hot subdwarf stars of mixed H/He composition exists, showing extreme metal abundance anomalies. Whether such intermediate helium-rich (iHe) subdwarf stars provide an evolutionary link between the dominant classes is an open question. The presence of strong Ge, Sn, and Pb lines in the UV spectrum of HZ$\,$44 suggests a strong enrichment of heavy elements in this iHe-sdO star and calls for a detailed quantitative spectral analysis focusing on trans-iron elements. Non-LTE model atmospheres calculated with TLUSTY are combined with high-quality optical, UV and FUV spectra of HZ$\,$44 and its hotter sibling HD$\,$127493 to determine their atmospheric parameters and metal abundance patterns. By collecting atomic data from literature we succeeded to determine abundances of 29 metals in HZ$\,$44, including the trans-iron elements Ga, Ge, As, Se, Zr, Sn, and Pb and provide upper limits for 10 other metals. This makes it the best described hot subdwarf in terms of chemical composition. For HD$\,$127493 the abundance of 15 metals, including Ga, Ge, and Pb and upper limits for another 16 metals were derived. Heavy elements turn out to be overabundant by one to four orders of magnitude with respect to the Sun. Zr and Pb are among the most enriched elements. The C, N, and O abundance for both stars can be explained by nucleosynthesis of hydrogen burning in the CNO cycle along with their helium enrichment. On the other hand, the heavy-element anomalies are unlikely to be caused by nucleosynthesis. Instead diffusion processes are evoked with radiative levitation overcoming gravitational settlement of the heavy elements.


title: Characteristic classes of bundles of K3 manifolds and the Nielsen realization problem

abstract: Let $K$ be the K3 manifold. In this note, we discuss two methods to prove that certain generalized Miller--Morita--Mumford classes for smooth bundles with fiber $K$ are non-zero. As a consequence, we fill a gap in a paper of the first author, and prove that the homomorphism $Diff(K)\to \pi_0 Diff(K)$ does not split. One of the two methods of proof uses a result of Franke on the stable cohomology of arithmetic groups that strengthens work of Borel, and may be of independent interest.


title: Patient-specific Conditional Joint Models of Shape, Image Features and Clinical Indicators

abstract: We propose and demonstrate a joint model of anatomical shapes, image features and clinical indicators for statistical shape modeling and medical image analysis. The key idea is to employ a copula model to separate the joint dependency structure from the marginal distributions of variables of interest. This separation provides flexibility on the assumptions made during the modeling process. The proposed method can handle binary, discrete, ordinal and continuous variables. We demonstrate a simple and efficient way to include binary, discrete and ordinal variables into the modeling. We build Bayesian conditional models based on observed partial clinical indicators, features or shape based on Gaussian processes capturing the dependency structure. We apply the proposed method on a stroke dataset to jointly model the shape of the lateral ventricles, the spatial distribution of the white matter hyperintensity associated with periventricular white matter disease, and clinical indicators. The proposed method yields interpretable joint models for data exploration and patient-specific statistical shape models for medical image analysis.


title: Neutron Scattering Study of Magnetic Anisotropy in a Tetragonal Antiferromagnet Bi$_2$CuO$_4$

abstract: We present a comprehensive study of magnon excitations in the tetragonal easy-plane anti-ferromagnet Bi$_2$CuO$_4$ using inelastic neutron scattering and spin wave analyses. The nature of low energy magnons, and hence the anisotropy in this material, has been controversial. We show unambiguously that the low energy magnon spectrum consists of a gapped and a gapless mode, which we attribute to out-of-plane and in-plane spin fluctuations, respectively. We modelled the observed magnon spectrum using linear spin wave analysis of a minimal anisotropic spin model motivated by the lattice symmetry. By studying the magnetic field dependence of the (1, 0, 0) Bragg peak intensity and the in-plane magnon intensity, we observed a spin-flop transition in the $ab$ plane at $\sim0.4$~T which directly indicates the existence of a small in-plane anisotropy that is classically forbidden. It is only by taking into account magnon zero-point fluctuations beyond the linear spin wave approximation, we could explain this in-plane anisotropy and its magnitude, the latter of which is deduced from critical field of the spin-flop transition. The microscopic origins of the observed anisotropic interactions are also discussed. We found that our data is inconsistent with a large Dzyaloshinskii-Moriya interaction, which suggests a potential departure of Bi$_2$CuO$_4$ from the conventional theories of magnetic anisotropy for other cuprates.


title: Epitaxial stabilization of thin films of the frustrated Ge-based spinels

abstract: Frustrated magnets can host numerous exotic many-body quantum and topological phenomena. GeNi$_2$O$_4$ is a three dimensional $S=1$ frustrated magnet with an unusual two-stage transition to the two-dimensional antiferromagnetic ground state, while GeCu$_2$O$_4$ is a high-pressure phase with a strongly tetragonally elongated spinel structure and magnetic lattice formed by $S=1/2$ CuO$_2$ linear chains with frustrated exchange interactions and exotic magnetic behavior. Here we report on the first thin-film epitaxial stabilization of these two compounds. Developed growth mode, surface morphology, crystal structure and copper valence state were characterized by in-situ reflection high-energy electron diffraction, atomic force microscopy, X-ray reflectivity, X-ray diffraction, X-ray photoelectron spectroscopy and resonant X-ray absorption spectroscopy. Our results pave an alternative route to the comprehensive investigation of the puzzling magnetic properties of these compounds and exploration of novel emergent features driven by strain.


title: Product Aesthetic Design: A Machine Learning Augmentation

abstract: Aesthetics are critically important to market acceptance. In the automotive industry, an improved aesthetic design can boost sales by 30% or more. Firms invest heavily in designing and testing aesthetics. A single automotive "theme clinic" can cost over $100,000, and hundreds are conducted annually. We propose a model to augment the commonly-used aesthetic design process by predicting aesthetic scores and automatically generating innovative and appealing product designs. The model combines a probabilistic variational autoencoder (VAE) with adversarial components from generative adversarial networks (GAN) and a supervised learning component. We train and evaluate the model with data from an automotive partner-images of 203 SUVs evaluated by targeted consumers and 180,000 high-quality unrated images. Our model predicts well the appeal of new aesthetic designs-43.5% improvement relative to a uniform baseline and substantial improvement over conventional machine learning models and pretrained deep neural networks. New automotive designs are generated in a controllable manner for use by design teams. We empirically verify that automatically generated designs are (1) appealing to consumers and (2) resemble designs which were introduced to the market five years after our data were collected. We provide an additional proof-of-concept application using opensource images of dining room chairs.


title: Constraining the astrophysics and cosmology from 21cm tomography using deep learning with the SKA

abstract: Future Square Kilometre Array (SKA) surveys are expected to generate huge datasets of 21cm maps on cosmological scales from the Epoch of Reionization (EoR). We assess the viability of exploiting machine learning techniques, namely, convolutional neural networks (CNN), to simultaneously estimate the astrophysical and cosmological parameters from 21cm maps from semi-numerical simulations. We further convert the simulated 21cm maps into SKA-like mock maps using the detailed SKA antennae distribution, thermal noise and a recipe for foreground cleaning. We successfully design two CNN architectures (VGGNet-like and ResNet-like) that are both efficiently able to extract simultaneously three astrophysical parameters, namely the photon escape fraction (f$_{\rm esc}$), the ionizing emissivity power dependence on halo mass ($C_{\rm ion}$) and the ionizing emissivity redshift evolution index ($D_{\rm ion}$), and three cosmological parameters, namely the matter density parameter ($\Omega_{m}$), the dimensionless Hubble constant ($h$), and the matter fluctuation amplitude ($\sigma_{8}$), from 21cm maps at several redshifts. With the presence of noise from SKA, our designed CNNs are still able to recover these astrophysical and cosmological parameters with great accuracy ($R^{2} > 92\%$), improving to $R^{2} > 99\%$ towards low redshift and low neutral fraction values. Our results show that future 21cm observations can play a key role to break degeneracy between models and tightly constrain the astrophysical and cosmological parameters, using only few frequency channels.


title: SubTSBR to tackle high noise and outliers for data-driven discovery of differential equations

abstract: Data-driven discovery of differential equations has been an emerging research topic. We propose a novel algorithm subsampling-based threshold sparse Bayesian regression (SubTSBR) to tackle high noise and outliers. The subsampling technique is used for improving the accuracy of the Bayesian learning algorithm. It has two parameters: subsampling size and the number of subsamples. When the subsampling size increases with fixed total sample size, the accuracy of our algorithm goes up and then down. When the number of subsamples increases, the accuracy of our algorithm keeps going up. We demonstrate how to use our algorithm step by step and compare our algorithm with threshold sparse Bayesian regression (TSBR) for the discovery of differential equations. We show that our algorithm produces better results. We also discuss the merits of discovering differential equations from data and demonstrate how to discover models with random initial and boundary condition as well as models with bifurcations. The numerical examples are: (1) predator-prey model with noise, (2) shallow water equations with outliers, (3) heat diffusion with random initial and boundary condition, and (4) fish-harvesting problem with bifurcations.


title: Genome-wide Causation Studies of Complex Diseases

abstract: Despite significant progress in dissecting the genetic architecture of complex diseases by genome-wide association studies (GWAS), the signals identified by association analysis may not have specific pathological relevance to diseases so that a large fraction of disease causing genetic variants is still hidden. Association is used to measure dependence between two variables or two sets of variables. Genome-wide association studies test association between a disease and SNPs (or other genetic variants) across the genome. Association analysis may detect superficial patterns between disease and genetic variants. Association signals provide limited information on the causal mechanism of diseases. The use of association analysis as a major analytical platform for genetic studies of complex diseases is a key issue that hampers discovery of the mechanism of diseases, calling into question the ability of GWAS to identify loci underlying diseases. It is time to move beyond association analysis toward techniques enabling the discovery of the underlying causal genetic strctures of complex diseases. To achieve this, we propose a concept of a genome-wide causation studies (GWCS) as an alternative to GWAS and develop additive noise models (ANMs) for genetic causation analysis. Type I error rates and power of the ANMs to test for causation are presented. We conduct GWCS of schizophrenia. Both simulation and real data analysis show that the proportion of the overlapped association and causation signals is small. Thus, we hope that our analysis will stimulate discussion of GWAS and GWCS.


title: Topological invariance of torsion sensitive intersection homology

abstract: Torsion sensitive intersection homology was introduced to unify several versions of Poincare duality for stratified spaces into a single theorem. This unified duality theorem holds with ground coefficients in an arbitrary PID and with no local cohomology conditions on the underlying space. In this paper we consider for torsion sensitive intersection homology analogues of another important property of classical intersection homology: topological invariance. In other words, we consider to what extent the defining sheaf complexes of the theory are independent (up to quasi-isomorphism) of choice of stratification. In addition to providing torsion sensitive versions of the existing invariance theorems for classical intersection homology, our techniques provide some new results even in the classical setting.


title: Organic complexity in protostellar disk candidates

abstract: We present ALMA observations of organic molecules towards five low-mass Class 0/I protostellar disk candidates in the Serpens cluster. Three sources (Ser-emb 1, Ser-emb 8, and Ser-emb 17) present emission of CH3OH as well as CH3OCH3, CH3OCHO, and CH2CO, while NH2CHO is detected in just Ser-emb 8 and Ser-emb 17. Detecting hot corino-type chemistry in three of five sources represents a high occurrence rate given the relative sparsity of these sources in the literature, and this suggests a possible link between protostellar disk formation and hot corino formation. For sources with CH3OH detections, we derive column densities of 10^{17}-10^{18} cm^{-2} and rotational temperatures of ~200-250 K. The CH3OH-normalized column density ratios of large, oxygen-bearing COMs in the Serpens sources and other hot corinos span two orders of magnitude, demonstrating a high degree of chemical diversity at the hot corino stage. Resolved observations of a larger sample of objects are needed to understand the origins of chemical diversity in hot corinos, and the relationship between different protostellar structural elements on disk-forming scales.


title: GRIP++: Enhanced Graph-based Interaction-aware Trajectory Prediction for Autonomous Driving

abstract: Despite the advancement in the technology of autonomous driving cars, the safety of a self-driving car is still a challenging problem that has not been well studied. Motion prediction is one of the core functions of an autonomous driving car. Previously, we propose a novel scheme called GRIP which is designed to predict trajectories for traffic agents around an autonomous car efficiently. GRIP uses a graph to represent the interactions of close objects, applies several graph convolutional blocks to extract features, and subsequently uses an encoder-decoder long short-term memory (LSTM) model to make predictions. Even though our experimental results show that GRIP improves the prediction accuracy of the state-of-the-art solution by 30%, GRIP still has some limitations. GRIP uses a fixed graph to describe the relationships between different traffic agents and hence may suffer some performance degradations when it is being used in urban traffic scenarios. Hence, in this paper, we describe an improved scheme called GRIP++ where we use both fixed and dynamic graphs for trajectory predictions of different types of traffic agents. Such an improvement can help autonomous driving cars avoid many traffic accidents. Our evaluations using a recently released urban traffic dataset, namely ApolloScape showed that GRIP++ achieves better prediction accuracy than state-of-the-art schemes. GRIP++ ranked #1 on the leaderboard of the ApolloScape trajectory competition in October 2019. In addition, GRIP++ runs 21.7 times faster than a state-of-the-art scheme, CS-LSTM.


title: A window into NV center kinetics via repeated annealing and spatial tracking of thousands of individual NV centers

abstract: Knowledge of the nitrogen-vacancy center formation kinetics in diamond is critical to engineering sensors and quantum information devices based on this defect. Here we utilize the longitudinal tracking of single NV centers to elucidate NV defect kinetics during high-temperature annealing from 800-1100 $^\circ$C in high-purity chemical-vapor-deposition diamond. We observe three phenomena which can coexist: NV formation, NV quenching, and NV orientation changes. Of relevance to NV-based applications, a 6 to 24-fold enhancement in the NV density, in the absence of sample irradiation, is observed by annealing at 980 $^\circ$C, and NV orientation changes are observed at 1050 $^\circ$C. With respect to the fundamental understanding of defect kinetics in ultra-pure diamond, our results indicate a significant vacancy source can be activated for NV creation between 950-980 $^\circ$C and suggests that native hydrogen from NVH$_y$ complexes plays a dominant role in NV quenching, in agreement with recent {\it ab initio} calculations. Finally, the direct observation of orientation changes allows us to estimate an NV diffusion barrier of 5.1~eV.


title: Generating Correctness Proofs with Neural Networks

abstract: Foundational verification allows programmers to build software which has been empirically shown to have high levels of assurance in a variety of important domains. However, the cost of producing foundationally verified software remains prohibitively high for most projects,as it requires significant manual effort by highly trained experts. In this paper we present Proverbot9001,a proof search system using machine learning techniques to produce proofs of software correctness in interactive theorem provers. We demonstrate Proverbot9001 on the proof obligations from a large practical proof project,the CompCert verified C compiler,and show that it can effectively automate what were previously manual proofs,automatically producing proofs for 27.5% of theorem statements in our test dataset, when combined with solver-based tooling. Without any additional solvers,we exhibit a proof completion rate that is a 4X improvement over prior state-of-the-art machine learning models for generating proofs in Coq.


title: Efficient computation of the Jacobi symbol

abstract: The family of left-to-right GCD algorithms reduces input numbers by repeatedly subtracting the smaller number, or multiple of the smaller number, from the larger number. This paper describes how to extend any such algorithm to compute the Jacobi symbol, using a single table lookup per reduction. For both quadratic time GCD algorithms (Euclid, Lehmer) and subquadratic algorithms (Knuth, Sch\"onhage, M\"oller), the additional cost is linear, roughly one table lookup per quotient in the quotient sequence. This method was used for the 2010 rewrite of the Jacobi symbol computation in GMP.


title: An Ongoing Project to Improve the Rectilinear and the Pseudolinear Crossing Constants

abstract: A drawing of a graph in the plane is {\it pseudolinear} if the edges of the drawing can be extended to doubly-infinite curves that form an arrangement of pseudolines, that is, any pair of edges crosses precisely once. A special case are {\it rectilinear} drawings where the edges of the graph are drawn as straight line segments. The rectilinear (pseudolinear) crossing number of a graph is the minimum number of pairs of edges of the graph that cross in any of its rectilinear (pseudolinear) drawings. In this paper we describe an ongoing project to continuously obtain better asymptotic upper bounds on the rectilinear and pseudolinear crossing number of the complete graph $K_n$.


title: One-relator quotients of Partially Commutative Groups

abstract: We generalise a key result of one-relator group theory, namely Magnus's Freiheitssatz, to partially commutative groups, under sufficiently strong conditions on the relator. The main theorem shows that under our conditions, on an element $r$ of a partially commutative group $\mathbb{G}$, certain Magnus subgroups embed in the quotient $G=\mathbb{G}/N(r)$; that if $r=s^n$ has root $s$ in $\mathbb{G}$ then the order of $s$ in $G$ is $n$, and under slightly stronger conditions that the word problem of $G$ is decidable. We also give conditions under which the question of which Magnus subgroups of $\mathbb{G}$ embed in $G$ reduces to the same question in the minimal parabolic subgroup of $\mathbb{G}$ containing $r$. In many cases this allows us to characterise Magnus subgroups which embed in $G$, via a condition on $r$ and the commutation graph of $\mathbb{G}$, and to find further examples of quotients $G$ where the word and conjugacy problems are decidable. We give evidence that situations in which our main theorem applies are not uncommon, by proving that for cycle graphs with a chord $\Gamma$, almost all cyclically reduced elements of the partially commutative group $\mathbb{G}(\Gamma)$ satisfy the conditions of the theorem.


title: 3.8um Imaging of 400-600K Brown Dwarfs and Orbital Constraints for WISEP J045853.90+643452.6AB

abstract: Half of the energy emitted by late-T- and Y-type brown dwarfs emerges at 3.5 < lambda um < 5.5. We present new L' (3.43 < lambda um < 4.11) photometry obtained at the Gemini North telescope for nine late-T and Y dwarfs, and synthesize L' from spectra for an additional two dwarfs. The targets include two binary systems which were imaged at a resolution of 0.25". One of these, WISEP J045853.90+643452.6AB, shows significant motion, and we present an astrometric analysis of the binary using Hubble Space Telescope, Keck Adaptive Optics, and Gemini images. We compare lambda ~4um observations to models, and find that the model fluxes are too low for brown dwarfs cooler than ~700K. The discrepancy increases with decreasing temperature, and is a factor of ~2 at T_eff=500K and ~4 at T_eff=400K. Warming the upper layers of a model atmosphere generates a spectrum closer to what is observed. The thermal structure of cool brown dwarf atmospheres above the radiative-convective boundary may not be adequately modelled using pure radiative equilibrium; instead heat may be introduced by thermochemical instabilities (previously suggested for the L- to T-type transition) or by breaking gravity waves (previously suggested for the solar system giant planets). One-dimensional models may not capture these atmospheres, which likely have both horizontal and vertical pressure/temperature variations.


title: Confining Strings, Infinite Statistics and Integrability

abstract: We study confining strings in massive adjoint two-dimensional chromodynamics. Off-shell, as a consequence of zigzag formation, the resulting worldsheet theory provides a non-trivial dynamical realization of infinite quon statistics. Taking the high energy limit we identify a remarkably simple and novel integrable relativistic $N$-body system. Its symmetry algebra contains an additional "shadow" Poincar\'e subalgebra. This model describes the $N$-particle subsector of a $T\bar T$-deformed massless fermion.


title: The Social Sciences Interdisciplinarity for Astronomy and Astrophysics -- Lessons from the History of NASA and Related Fields

abstract: In this paper we showcase the importance of understanding and measuring interdisciplinarity and other -disciplinarity concepts for all scientists, the role social sciences have historically played in NASA research and missions, the sparsity of social science interdisciplinarity in space and planetary sciences, including astronomy and astrophysics, while there is an imperative necessity for it, and the example of interdisciplinarity between social sciences and astrobiology. Ultimately we give voice to the scientists across all fields with respect to their needs, aspirations and experiences in their interdisciplinary work with social sciences through an ad-hoc survey we conducted within the Astro2020 Decadal Survey scientific community.

categories: astro-ph.IM astro-ph.CO astro-ph.EP astro-ph.GA astro-ph.HE


title: Iterated chromatic localisation

abstract: We study a certain monoid of endofunctors of the stable homotopy category that includes localizations with respect to finite unions of Morava $K$-theories. We work in an axiomatic framework that can also be applied to analogous questions in equivariant stable homotopy theory. Our results should be helpful for the study of transchromatic phenomena, including the Chromatic Splitting Conjecture. The combinatorial parts of this work have been formalised in the Lean proof assistant.


title: Multi-Purposing Domain Adaptation Discriminators for Pseudo Labeling Confidence

abstract: Often domain adaptation is performed using a discriminator (domain classifier) to learn domain-invariant feature representations so that a classifier trained on labeled source data will generalize well to unlabeled target data. A line of research stemming from semi-supervised learning uses pseudo labeling to directly generate "pseudo labels" for the unlabeled target data and trains a classifier on the now-labeled target data, where the samples are selected or weighted based on some measure of confidence. In this paper, we propose multi-purposing the discriminator to not only aid in producing domain-invariant representations but also to provide pseudo labeling confidence.


title: Syntax and Stack Overflow: A methodology for extracting a corpus of syntax errors and fixes

abstract: One problem when studying how to find and fix syntax errors is how to get natural and representative examples of syntax errors. Most syntax error datasets are not free, open, and public, or they are extracted from novice programmers and do not represent syntax errors that the general population of developers would make. Programmers of all skill levels post questions and answers to Stack Overflow which may contain snippets of source code along with corresponding text and tags. Many snippets do not parse, thus they are ripe for forming a corpus of syntax errors and corrections. Our primary contribution is an approach for extracting natural syntax errors and their corresponding human made fixes to help syntax error research. A Python abstract syntax tree parser is used to determine preliminary errors and corrections on code blocks extracted from the SOTorrent data set. We further analyzed our code by executing the corrections in a Python interpreter. We applied our methodology to produce a public data set of 62,965 Python Stack Overflow code snippets with corresponding tags, errors, and stack traces. We found that errors made by Stack Overflow users do not match errors made by student developers or random mutations, implying there is a serious representativeness risk within the field. Finally we share our dataset openly so that future researchers can re-use and extend our syntax errors and fixes.


title: OmniNet: A unified architecture for multi-modal multi-task learning

abstract: Transformer is a popularly used neural network architecture, especially for language understanding. We introduce an extended and unified architecture that can be used for tasks involving a variety of modalities like image, text, videos, etc. We propose a spatio-temporal cache mechanism that enables learning spatial dimension of the input in addition to the hidden states corresponding to the temporal input sequence. The proposed architecture further enables a single model to support tasks with multiple input modalities as well as asynchronous multi-task learning, thus we refer to it as OmniNet. For example, a single instance of OmniNet can concurrently learn to perform the tasks of part-of-speech tagging, image captioning, visual question answering and video activity recognition. We demonstrate that training these four tasks together results in about three times compressed model while retaining the performance in comparison to training them individually. We also show that using this neural network pre-trained on some modalities assists in learning unseen tasks such as video captioning and video question answering. This illustrates the generalization capacity of the self-attention mechanism on the spatio-temporal cache present in OmniNet.


title: Novel Receiver for Superparamagnetic Iron Oxide Nanoparticles in a Molecular Communication Setting

abstract: Superparamagnetic iron oxide nanoparticles (SPIONs) have recently been introduced as information carriers in a testbed for molecular communication (MC) in duct flow. Here, a new receiver for this testbed is presented, based on the concept of a bridge circuit. The capability for a reliable transmission using the testbed and detection of the proposed receiver was evaluated by sending a text message and a 80 bit random sequence at a bit rate of 1/s, which resulted in a bit error rate of 0 %. Furthermore, the sensitivity of the device was assessed by a dilution series, which gave a limit for the detectability of peaks between 0.1 to 0.5 mg/mL. Compared to the commercial susceptometer that was previously used as receiver, the new detector provides an increased sampling rate of 100 samples/s and flexibility in the dimensions of the propagation channel. Furthermore, it allows to implement both single-ended and differential signaling in SPION-bases MC testbeds.


title: Optimization of a partial differential equation on a complex network

abstract: Differential equations on metric graphs can describe many phenomena in the physical world but also the spread of information on social media. To efficiently compute the solution is a hard task in numerical analysis. Solving a design problem, where the optimal setup for a desired state is given, is even more challenging. In this work, we focus on the task of solving an optimization problem subject to a differential equation on a metric graph with the control defined on a small set of Dirichlet nodes. We discuss the discretization by finite elements and provide rigorous error bounds as well as an efficient preconditioning strategy to deal with the large-scale case. We show in various examples that the method performs very robustly.


title: A fully 3D multi-path convolutional neural network with feature fusion and feature weighting for automatic lesion identification in brain MRI images

abstract: We propose a fully 3D multi-path convolutional network to predict stroke lesions from 3D brain MRI images. Our multi-path model has independent encoders for different modalities containing residual convolutional blocks, weighted multi-path feature fusion from different modalities, and weighted fusion modules to combine encoder and decoder features. Compared to existing 3D CNNs like DeepMedic, 3D U-Net, and AnatomyNet, our networks achieves the highest statistically significant cross-validation accuracy of 60.5% on the large ATLAS benchmark of 220 patients. We also test our model on multi-modal images from the Kessler Foundation and Medical College Wisconsin and achieve a statistically significant cross-validation accuracy of 65%, significantly outperforming the multi-modal 3D U-Net and DeepMedic. Overall our model offers a principled, extensible multi-path approach that outperforms multi-channel alternatives and achieves high Dice accuracies on existing benchmarks.


title: Lagrangian formulation for electric charge in a magnetic monopole distribution

abstract: We give a Lagrangian description of an electric charge in a field sourced by a continuous magnetic monopole distribution. The description is made possible thanks to a doubling of the configuration space. The Legendre transform of the nonrelativistic Lagrangian agrees with the Hamiltonian description given recently by Kupriyanov and Szabo. The covariant relativistic version of the Lagrangian is shown to introduce a new gauge symmetry, in addition to standard reparametrizations. The generalization of the system to open strings coupled to a magnetic monopole distribution is also given, as well as the generalization to particles in a non-Abelian gauge field which does not satisfy Bianchi identities in some region of the space-time.


title: Accounting for total variation and robustness in profiling health care providers

abstract: Monitoring outcomes of health care providers, such as patient deaths, hospitalizations and hospital readmissions, helps in assessing the quality of health care. We consider a large database on patients being treated at dialysis facilities in the United States, and the problem of identifying facilities with outcomes that are better than or worse than expected. Analyses of such data have been commonly based on random or fixed facility effects, which have shortcomings that can lead to unfair assessments. A primary issue is that they do not appropriately account for variation between providers that is outside the providers' control due, for example, to unobserved patient characteristics that vary between providers. In this article, we propose a smoothed empirical null approach that accounts for the total variation and adapts to different provider sizes. The linear model provides an illustration that extends easily to other nonlinear models for survival or binary outcomes, for example. The empirical null method is generalized to allow for some variation being due to quality of care. These methods are examined with numerical simulations and applied to the monitoring of survival in the dialysis facility data.


title: Stability selection enables robust learning of partial differential equations from limited noisy data

abstract: We present a statistical learning framework for robust identification of partial differential equations from noisy spatiotemporal data. Extending previous sparse regression approaches for inferring PDE models from simulated data, we address key issues that have thus far limited the application of these methods to noisy experimental data, namely their robustness against noise and the need for manual parameter tuning. We address both points by proposing a stability-based model selection scheme to determine the level of regularization required for reproducible recovery of the underlying PDE. This avoids manual parameter tuning and provides a principled way to improve the method's robustness against noise in the data. Our stability selection approach, termed PDE-STRIDE, can be combined with any sparsity-promoting penalized regression model and provides an interpretable criterion for model component importance. We show that in particular the combination of stability selection with the iterative hard-thresholding algorithm from compressed sensing provides a fast, parameter-free, and robust computational framework for PDE inference that outperforms previous algorithmic approaches with respect to recovery accuracy, amount of data required, and robustness to noise. We illustrate the performance of our approach on a wide range of noise-corrupted simulated benchmark problems, including 1D Burgers, 2D vorticity-transport, and 3D reaction-diffusion problems. We demonstrate the practical applicability of our method on real-world data by considering a purely data-driven re-evaluation of the advective triggering hypothesis for an embryonic polarization system in C.~elegans. Using fluorescence microscopy images of C.~elegans zygotes as input data, our framework is able to recover the PDE model for the regulatory reaction-diffusion-flow network of the associated proteins.


title: A Simpler Approach to Linear Programming

abstract: Dantzig and Eaves claimed that fundamental duality theorems of linear programming were a trivial consequence of Fourier elimination. Another property of Fourier elimination is considered here, regarding the existence of implicit equalities rather than solvability. This leads to a different interpretation of duality theory which allows us to use Gaussian elimination to decide solvability of systems of linear inequalities, for bounded systems.


title: Birational geometry for the covering of a nilpotent orbit closure

abstract: A nilpotent orbit $O$ of a complex semisimple Lie algebra $\mathfrak{g}$ has finite fundamental group. Associated with an etale cover of $O$, we have a finite cover of the closure $\bar{O}$ of $O$. In this article we consider the finite cover X associated with the universal cover of a nilpotent orbit $O$ of a classical simple Lie algebra $\mathfrak{g}$. We construct explicitly a Q-factorial terminalization of X in a group theoretic way.


title: Multi-Scale Process Modelling and Distributed Computation for Spatial Data

abstract: Recent years have seen a huge development in spatial modelling and prediction methodology, driven by the increased availability of remote-sensing data and the reduced cost of distributed-processing technology. It is well known that modelling and prediction using infinite-dimensional process models is not possible with large data sets, and that both approximate models and, often, approximate-inference methods, are needed. The problem of fitting simple global spatial models to large data sets has been solved through the likes of multi-resolution approximations and nearest-neighbour techniques. Here we tackle the next challenge, that of fitting complex, nonstationary, multi-scale models to large data sets. We propose doing this through the use of superpositions of spatial processes with increasing spatial scale and increasing degrees of nonstationarity. Computation is facilitated through the use of Gaussian Markov random fields and parallel Markov chain Monte Carlo based on graph colouring. The resulting model allows for both distributed computing and distributed data. Importantly, it provides opportunities for genuine model and data scaleability and yet is still able to borrow strength across large spatial scales. We illustrate a two-scale version on a data set of sea-surface temperature containing on the order of one million observations, and compare our approach to state-of-the-art spatial modelling and prediction methods.


title: Gon\v{c}arov Polynomials in Partition Lattices and Exponential Families

abstract: Classical Gon\v{c}arov polynomials arose in numerical analysis as a basis for the solutions of the Gon\v{c}arov interpolation problem. These polynomials provide a natural algebraic tool in the enumerative theory of parking functions. By replacing the differentiation operator with a delta operator and using the theory of finite operator calculus, Lorentz, Tringali and Yan introduced the sequence of generalized Gon\v{c}arov polynomials associated to a pair $(\Delta, Z)$ of a delta operator $\Delta$ and an interpolation grid $Z$. Generalized Gon\v{c}arov polynomials share many nice algebraic properties and have a connection with the theories of binomial enumeration and order statistics. In this paper we give a complete combinatorial interpretation for any sequence of generalized Gon\v{c}arov polynomials. First, we show that they can be realized as weight enumerators in partition lattices. Then, we give a more concrete realization in exponential families and show that these polynomials enumerate various enriched structures of vector parking functions.


title: Degrees of Randomized Computability

abstract: In this survey we discuss work of Levin and V'yugin on collections of sequences that are non-negligible in the sense that they can be computed by a probabilistic algorithm with positive probability. More precisely, Levin and V'yugin introduced an ordering on collections of sequences that are closed under Turing equivalence. Roughly speaking, given two such collections $\mathcal{A}$ and $\mathcal{B}$, $\mathcal{A}$ is below $\mathcal{B}$ in this ordering if $\mathcal{A}\setminus\mathcal{B}$ is negligible. The degree structure associated with this ordering, the Levin-V'yugin degrees (or LV-degrees), can be shown to be a Boolean algebra, and in fact a measure algebra. We demonstrate the interactions of this work with recent results in computability theory and algorithmic randomness: First, we recall the definition of the Levin-V'yugin algebra and identify connections between its properties and classical properties from computability theory. In particular, we apply results on the interactions between notions of randomness and Turing reducibility to establish new facts about specific LV-degrees, such as the LV-degree of the collection of 1-generic sequences, that of the collection of sequences of hyperimmune degree, and those collections corresponding to various notions of effective randomness. Next, we provide a detailed explanation of a complex technique developed by V'yugin that allows the construction of semi-measures into which computability-theoretic properties can be encoded. We provide two examples of the use of this technique by explicating a result of V'yugin's about the LV-degree of the collection of Martin-L\"of random sequences and extending the result to the LV-degree of the collection of sequences of DNC degree.


title: Unsupervised Task Design to Meta-Train Medical Image Classifiers

abstract: Meta-training has been empirically demonstrated to be the most effective pre-training method for few-shot learning of medical image classifiers (i.e., classifiers modeled with small training sets). However, the effectiveness of meta-training relies on the availability of a reasonable number of hand-designed classification tasks, which are costly to obtain, and consequently rarely available. In this paper, we propose a new method to unsupervisedly design a large number of classification tasks to meta-train medical image classifiers. We evaluate our method on a breast dynamically contrast enhanced magnetic resonance imaging (DCE-MRI) data set that has been used to benchmark few-shot training methods of medical image classifiers. Our results show that the proposed unsupervised task design to meta-train medical image classifiers builds a pre-trained model that, after fine-tuning, produces better classification results than other unsupervised and supervised pre-training methods, and competitive results with respect to meta-training that relies on hand-designed classification tasks.


title: Scheduling Discovery in the 2020s

abstract: The 2020s will be the most data-rich decade of astronomy in history. As the scale and complexity of our surveys increase, the problem of scheduling becomes more critical. We must develop high-quality scheduling approaches, implement them as open-source software, and begin linking the typically separate stages of observation and data analysis.


title: Decoding the Style and Bias of Song Lyrics

abstract: The central idea of this paper is to gain a deeper understanding of song lyrics computationally. We focus on two aspects: style and biases of song lyrics. All prior works to understand these two aspects are limited to manual analysis of a small corpus of song lyrics. In contrast, we analyzed more than half a million songs spread over five decades. We characterize the lyrics style in terms of vocabulary, length, repetitiveness, speed, and readability. We have observed that the style of popular songs significantly differs from other songs. We have used distributed representation methods and WEAT test to measure various gender and racial biases in the song lyrics. We have observed that biases in song lyrics correlate with prior results on human subjects. This correlation indicates that song lyrics reflect the biases that exist in society. Increasing consumption of music and the effect of lyrics on human emotions makes this analysis important.


title: Collective Heavy Top Dynamics

abstract: We construct a Poisson map $\mathbf{M}\colon T^{*}\mathbb{C}^{2} \to \mathfrak{se}(3)^{*}$ with respect to the canonical Poisson bracket on $T^{*}\mathbb{C}^{2} \cong T^{*}\mathbb{R}^{4}$ and the $(-)$-Lie--Poisson bracket on the dual $\mathfrak{se}(3)^{*}$ of the Lie algebra of the special Euclidean group $\mathsf{SE}(3)$. The essential part of this map is the momentum map associated with the cotangent lift of the natural right action of the semidirect product Lie group $\mathsf{SU}(2) \ltimes \mathbb{C}^{2}$ on $\mathbb{C}^{2}$. This Poisson map gives rise to a canonical Hamiltonian system on $T^{*}\mathbb{C}^{2}$ whose solutions are mapped by $\mathbf{M}$ to solutions of the heavy top equations. We show that the Casimirs of the heavy top dynamics and the additional conserved quantity of the Lagrange top correspond to the Noether conserved quantities associated with certain symmetries of the canonical Hamiltonian system. We also construct a Lie--Poisson integrator for the heavy top dynamics by combining the Poisson map $\mathbf{M}$ with a simple symplectic integrator, and demonstrate that the integrator exhibits either exact or near conservation of the conserved quantities of the Kovalevskaya top.


title: Pressure dependence of the superconducting transition temperature of compressed LaH$_{10}$

abstract: Two recent experiments [M. Somayazulu $et$ $al$., Phys. Rev. Lett. ${\bf 122}$, 027001 (2019) and A. P. Drozdov $et$ $al$., Nature ${\bf 569}$, 528 (2019)] reported the discovery of superconductivity in the fcc phase of LaH$_{10}$ at a critical temperature $T_{\rm c}$ between 250${\sim}$260 K under a pressure of about 170 GPa. However, the dependence of $T_{\rm c}$ on pressure showed different patterns: i.e., the former experiment observed a continuous increase of $T_{\rm c}$ up to ${\sim}$275 K on further increase of pressure to 202 GPa, while the latter one observed an abrupt decrease of $T_{\rm c}$ with increasing pressure. Here, based on first-principles calculations, we reveal that for the fcc-LaH$_{10}$ phase, softening of the low-frequency optical phonon modes of H atoms dramatically occurs as pressure decreases, giving rise to a significant increase of the electron-phonon coupling (EPC) constant. Meanwhile, the electronic band structure near the Fermi energy is insensitive to change with respect to pressure. These results indicate that the pressure-dependent phonon softening is unlikely associated with Fermi-surface nesting, but driven by effective screening with the electronic states near the Fermi energy. It is thus demonstrated that the strong variation of EPC with respect to pressure plays a dominant role in the decrease of $T_{\rm c}$ with increasing pressure, supporting the measurements of Drozdov $et$ $al$.


title: Before and beyond the Wilson-Cowan equations

abstract: The Wilson-Cowan equations represent a landmark in the history of computational neuroscience. Among the insights Wilson and Cowan offered for neuroscience, they crystallized an approach to modeling neural dynamics and brain function. Although their iconic equations are used in various guises today, the ideas that led to their formulation and the relationship to other approaches are not well known. Here, we give a little context to some of the biological and theoretical concepts that lead to the Wilson-Cowan equations and discuss how to extend beyond them.


title: Pointed Hopf actions on central simple division algebras

abstract: We examine actions of finite-dimensional pointed Hopf algebras on central simple division algebras in characteristic 0. (By a Hopf action we mean a Hopf module algebra structure.) In all examples considered, we show that the given Hopf algebra does admit a faithful action on a central simple division algebra, and we construct such a division algebra. This is in contrast to earlier work of Etingof and Walton, in which it was shown that most pointed Hopf algebras do not admit faithful actions on fields. We consider all bosonizations of Nichols algebras of finite Cartan type, small quantum groups, generalized Taft algebras with non-nilpotent skew primitive generators, and an example of non-Cartan type.


title: Remote sensing of angular scattering effect of aerosols in a North American megacity

abstract: The angle-dependent scattering effect of aerosols in the atmosphere can be used to infer their compositions, which in turn is important to understand their impacts of human health and Earth climate. The aerosol phase function, which characterizes the angular signature of scattering, has been continuously monitored from ground-based and space-borne observations. However, the range of scattering angles these instruments can sample is very limited. There is a dearth of research on the remote sensing of aerosol angular scattering effect at a city scale that analyzes diurnal variability and includes a wide range of scattering angles. Here, we quantify the aerosol angular scattering effect using measurements from a mountain-top remote sensing instrument: the California Laboratory for Atmospheric Remote Sensing Fourier Transform Spectrometer (CLARS-FTS). CLARS-FTS is located on top of the Mt. Wilson (1.67km above sea level) overlooking the Los Angeles (LA) megacity and receives reflected sunlight from targeted surface reflection points. The observational geometries of CLARS-FTS provide a wide range of scattering angles, from about 20 degrees (forward) to about 140 degrees (backward). The O2 ratio, which is the ratio of retrieved O2 Slant Column Density (SCD) to geometric O2 SCD, quantifies the aerosol transmission with a value of 1.0 represent aerosol-free and with a value closer to 0.0 represents stronger aerosol loadings. The aerosol transmission quantified by the O2 ratio from CLARS measurements provides an effective indicator of the aerosol scattering effect.


title: Derivation of holographic negativity in AdS$_3$/CFT$_2$

abstract: We present a derivation of the holographic dual of logarithmic negativity in AdS$_3$/CFT$_2$ that was recently conjectured in [Phys. Rev. D 99, 106014 (2019)]. This is given by the area of an extremal cosmic brane that terminates on the boundary of the entanglement wedge. The derivation consists of relating the recently introduced R\'enyi reflected entropy to the logarithmic negativity in holographic conformal field theories. Furthermore, we clarify previously mysterious aspects of negativity at large central charge seen in conformal blocks and comment on generalizations to generic dimensions, dynamical settings, and quantum corrections.


title: Search-Based Motion Planning for Performance Autonomous Driving

abstract: Driving on the limits of vehicle dynamics requires predictive planning of future vehicle states. In this work, a search-based motion planning is used to generate suitable reference trajectories of dynamic vehicle states with the goal to achieve the minimum lap time on slippery roads. The search-based approach enables to explicitly consider a nonlinear vehicle dynamics model as well as constraints on states and inputs so that even challenging scenarios can be achieved in a safe and optimal way. The algorithm performance is evaluated in simulated driving on a track with segments of different curvatures.


title: Comparison of Classical Machine Learning Approaches on Bangla Textual Emotion Analysis

abstract: Detecting emotions from text is an extension of simple sentiment polarity detection. Instead of considering only positive or negative sentiments, emotions are conveyed using more tangible manner; thus, they can be expressed as many shades of gray. This paper manifests the results of our experimentation for fine-grained emotion analysis on Bangla text. We gathered and annotated a text corpus consisting of user comments from several Facebook groups regarding socio-economic and political issues, and we made efforts to extract the basic emotions (sadness, happiness, disgust, surprise, fear, anger) conveyed through these comments. Finally, we compared the results of the five most popular classical machine learning techniques namely Naive Bayes, Decision Tree, k-Nearest Neighbor (k-NN), Support Vector Machine (SVM) and K-Means Clustering with several combinations of features. Our best model (SVM with a non-linear radial-basis function (RBF) kernel) achieved an overall average accuracy score of 52.98% and an F1 score (macro) of 0.3324


title: A study of multivalent q-starlike functions connected with circular domain

abstract: In the present article, our aim is to examine some useful problems including the convolution problem, sufficiency criteria, coefficient estimates and Fekete-Szego type inequalities for a new subfamily of analytic and multivalent functions associated with circular domain. In addition, we also define and study a Bernardi integral operator in its $q$-extension for multivalent functions.


title: The Sym(3) Conjecture and Alt(8)

abstract: We give an alternate computer-free proof of a result of Z. Arad, M. Muzychuk, and A. Oliver: if G is a minimal counterexample to the Sym(3) conjecture, then Soc(G)' cannot be isomorphic to Alt(8).


title: Technosignatures in the Thermal Infrared

abstract: WISE, Gaia, and JWST provide an opportunity to compute the first robust upper limits on the energy supplies of extraterrestrial civilizations, both for stars in the Galaxy (Kardashev Type II civilizations) and for other galaxies (Kardashev Type III civilizations). Together, they allow for a nearly-complete catalog of nearby stars with infrared excesses, which is valuable for both stellar astrophysics and searches for technosignatures; and of the MIR luminosities of galaxies, important for studies of galaxies' star and star-formation properties, but also for the identification of potential galaxies endemic with alien technology. JWST will provide the crucial mid-infrared spectroscopy necessary to identify the origin of these infrared excesses, advancing both traditional astronomy and searches for technosignatures. Such signatures are distinguished from dust by their lack of far-infrared emission and lack of association with star formation.


title: Technosignatures in Transit

abstract: Kepler, K2, TESS, and similar time-domain photometric projects, while designed with exoplanet detection in mind, are also well-suited projects for searches for large artificial structures orbiting other stars in the Galaxy. An effort to examine these data sets with an eye towards non-spherical or otherwise anomalous transit events, and a robust follow-up program to understand the stars and occulters that generate them, would enable the first robust upper limits on such megastructures in terms of their sizes, occurrence rates, and orbital properties. Such work also has the ancillary benefit of improving our understanding of stellar photometric variability and orbital and physical parameter estimation of exoplanets from photometric time series, and may lead to the identification of new, unexpected classes of stellar variables and exoplanets. Ultimately, searching for the most unusual and anomalous signatures benefits not only the search for technologies, but also the entire astronomical community by uncovering new mysteries to advance our understanding of the Universe.


title: Searches for Technosignatures in Astronomy and Astrophysics

abstract: The search for life beyond the Solar System-a major part of the Planetary Systems thematic area of the Astro2020 Decadal process-includes the search for technological life. Although financial support for such searches at the NSF and NASA and in past decadal surveys has been weak to nonexistent, recent advances in astrobiology, astrophysics, and advances in technical capability make searches for technosignatures a compelling theme for 2020-2030 and beyond.


title: Searches for Technosignatures: The State of the Profession

abstract: The search for life in the universe is a major theme of astronomy and astrophysics for the next decade. Searches for technosignatures are complementary to searches for biosignatures, in that they offer an alternative path to discovery, and address the question of whether complex (i.e. technological) life exists elsewhere in the Galaxy. This approach has been endorsed in prior Decadal Reviews and National Academies reports, and yet the field still receives almost no federal support in the US. Because of this lack of support, searches for technosignatures, precisely the part of the search of greatest public interest, suffers from a very small pool of trained practitioners. A major source of this issue is institutional inertia at NASA, which avoids the topic as a result of decades-past political grandstanding, conflation of the effort with non-scientific topics such as UFOs, and confusion regarding the scope of the term "SETI." The Astro2020 Decadal should address this issue by making developing the field an explicit priority for the next decade. It should recommend that NASA and the NSF support training and curricular development in the field in a way that supports equity and diversity, and make explicit calls for proposals to fund searches for technosignatures.


title: Approximating Constraint Satisfaction Problems on High-Dimensional Expanders

abstract: We consider the problem of approximately solving constraint satisfaction problems with arity $k > 2$ ($k$-CSPs) on instances satisfying certain expansion properties, when viewed as hypergraphs. Random instances of $k$-CSPs, which are also highly expanding, are well-known to be hard to approximate using known algorithmic techniques (and are widely believed to be hard to approximate in polynomial time). However, we show that this is not necessarily the case for instances where the hypergraph is a high-dimensional expander. We consider the spectral definition of high-dimensional expansion used by Dinur and Kaufman [FOCS 2017] to construct certain primitives related to PCPs. They measure the expansion in terms of a parameter $\gamma$ which is the analogue of the second singular value for expanding graphs. Extending the results by Barak, Raghavendra and Steurer [FOCS 2011] for 2-CSPs, we show that if an instance of MAX k-CSP over alphabet $[q]$ is a high-dimensional expander with parameter $\gamma$, then it is possible to approximate the maximum fraction of satisfiable constraints up to an additive error $\epsilon$ using $q^{O(k)} \cdot (k/\epsilon)^{O(1)}$ levels of the sum-of-squares SDP hierarchy, provided $\gamma \leq \epsilon^{O(1)} \cdot (1/(kq))^{O(k)}$. Based on our analysis, we also suggest a notion of threshold-rank for hypergraphs, which can be used to extend the results for approximating 2-CSPs on low threshold-rank graphs. We show that if an instance of MAX k-CSP has threshold rank $r$ for a threshold $\tau = (\epsilon/k)^{O(1)} \cdot (1/q)^{O(k)}$, then it is possible to approximately solve the instance up to additive error $\epsilon$, using $r \cdot q^{O(k)} \cdot (k/\epsilon)^{O(1)}$ levels of the sum-of-squares hierarchy. As in the case of graphs, high-dimensional expanders (with sufficiently small $\gamma$) have threshold rank 1 according to our definition.


title: Moderate deviations for the size of the giant component in a random hypergraph

abstract: We prove a moderate deviations principles for the size of the largest connected component in a random $d$-uniform hypergraph. The key tool is a version of the exploration process, that is also used to investigate the giant component of an Erd\"os-R\'enyi graph, a moderate deviations principle for the martingale associated with this exploration process, and exponential estimates.


title: EEG-Based Emotion Recognition Using Regularized Graph Neural Networks

abstract: Electroencephalography (EEG) measures the neuronal activities in different brain regions via electrodes. Many existing studies on EEG-based emotion recognition do not fully exploit the topology of EEG channels. In this paper, we propose a regularized graph neural network (RGNN) for EEG-based emotion recognition. RGNN considers the biological topology among different brain regions to capture both local and global relations among different EEG channels. Specifically, we model the inter-channel relations in EEG signals via an adjacency matrix in a graph neural network where the connection and sparseness of the adjacency matrix are inspired by neuroscience theories of human brain organization. In addition, we propose two regularizers, namely node-wise domain adversarial training (NodeDAT) and emotion-aware distribution learning (EmotionDL), to better handle cross-subject EEG variations and noisy labels, respectively. Extensive experiments on two public datasets, SEED and SEED-IV, demonstrate the superior performance of our model than state-of-the-art models in most experimental settings. Moreover, ablation studies show that the proposed adjacency matrix and two regularizers contribute consistent and significant gain to the performance of our RGNN model. Finally, investigations on the neuronal activities reveal important brain regions and inter-channel relations for EEG-based emotion recognition.


title: Multi-year Long-term Load Forecast for Area Distribution Feeders based on Selective Sequence Learning

abstract: Long-term load forecast (LTLF) for area distribution feeders is one of the most critical tasks frequently performed in electric distribution utility companies. For a specific planning area, cost-effective system upgrades can only be planned out based on accurate feeder LTLF results. In our previous research, we established a unique sequence prediction method which has the tremendous advantage of combining area top-down, feeder bottom-up and multi-year historical data all together for forecast and achieved a superior performance over various traditional methods by real-world tests. However, the previous method only focused on the forecast of the next one-year. In our current work, we significantly improved this method: the forecast can now be extended to a multi-year forecast window in the future; unsupervised learning techniques are used to group feeders by their load composition features to improve accuracy; we also propose a novel selective sequence learning mechanism which uses Gated Recurrent Unit network to not only learn how to predict sequence values but also learn to select the best-performing sequential configuration for each individual feeder. The proposed method was tested on an actual urban distribution system in West Canada. It was compared with traditional methods and our previous sequence prediction method. It demonstrates the best forecasting performance as well as the possibility of using sequence prediction models for multi-year component-level load forecast.


title: The relation between the independence number and rank of a signed graph

abstract: A signed graph $(G, \sigma)$ is a graph with a sign attached to each of its edges, where $G$ is the underlying graph of $(G, \sigma)$. Let $c(G)$, $\alpha(G)$ and $r(G, \sigma)$ be the cyclomatic number, the independence number and the rank of the adjacency matrix of $(G, \sigma)$, respectively. In this paper, we study the relation among the independence number, the rank and the cyclomatic number of a signed graph $(G, \sigma)$ with order $n$, and prove that $2n-2c(G) \leq r(G, \sigma)+2\alpha(G) \leq 2n$. Furthermore, the signed graphs that reaching the lower bound are investigated.


title: An inverse problem for a class of lacunary canonical systems with diagonal Hamiltonian

abstract: Hamiltonians are 2-by-2 positive semidefinite real symmetric matrix-valued functions satisfying certain conditions. In this paper, we solve the inverse problem for which recovers a Hamiltonian from the solution of a first-order system attached to a given Hamiltonian, consisting of ordinary differential equations parametrized by a set of complex numbers, under certain conditions for the solutions. This inverse problem is a generalization of the inverse problem for two-dimensional canonical systems.


title: Optimal strong approximation for quadrics over $\mathbb{F}_q[t]$

abstract: Suppose $q$ is a fixed odd prime power, $F(\vec{x})$ is a non-degenerate quadratic form over $\mathbb{F}_q[t]$ of discriminant $\Delta$ in $d\geq 5$ variables $\vec{x}$, and $f,g\in\mathbb{F}_q[t]$, $\boldsymbol{\lambda}\in\mathbb{F}_q[t]^d$. We show that whenever $\text{deg} f\geq (4+\varepsilon)\text{deg} g+O_{\varepsilon,F}(1)$, $\gcd(\Delta^{\infty},fg)=O(1)$, and the necessary local conditions are satisfied, we have a solution $\vec{x}\in\mathbb{F}_q[t]^d$ to $F(\vec{x})=f$ such that $\vec{x}\equiv\boldsymbol{\lambda}\bmod g$. For $d=4$, we show that the same conclusion holds if we instead have $\text{deg} f\geq (6+\varepsilon)\text{deg} g+O_{\varepsilon,F}(1)$. This gives us a new proof (independent of the Ramanujan conjecture over function fields proved by Drinfeld) that the diameter of any $k$-regular Morgenstern Ramanujan graphs $G$ is at most $(2+\varepsilon)\log_{k-1}|G|+O_{\varepsilon}(1)$. In contrast to the $d=4$ case, our result is optimal for $d\geq 5$. Our main new contributions are a stationary phase theorem over function fields for bounding oscillatory integrals, and a notion of anisotropic cones to circumvent isotropic phenomena in the function field setting.


title: Global Nonlinear Stability of Geodesic Solutions of Evolutionary Faddeev Model

abstract: In this paper, for evolutionary Faddeev model corresponding to maps from the Minkowski space $\mathbb{R}^{1+n}$ to the unit sphere $\mathbb{S}^2$, we show the global nonlinear stability of geodesic solutions, which are a kind of nontrivial and large solutions.


title: Optimal Downlink-Uplink Scheduling of Wireless Networked Control for Industrial IoT

abstract: This paper considers a wireless networked control system (WNCS) consisting of a dynamic system to be controlled (i.e., a plant), a sensor, an actuator and a remote controller for mission-critical Industrial Internet of Things (IIoT) applications. A WNCS has two types of wireless transmissions, i.e., the sensor's measurement transmission to the controller and the controller's command transmission to the actuator. In this work, we consider a practical half-duplex controller, which introduces a novel transmission-scheduling problem for WNCSs. A frequent scheduling of sensor's transmission results in a better estimation of plant states at the controller and thus a higher quality of control command, but it leads to a less frequent/timely control of the plant. Therefore, considering the overall control performance of the plant in terms of its average cost function, there exists a fundamental tradeoff between the sensor's and the controller's transmissions. We formulate a new problem to optimize the transmission-scheduling policy for minimizing the long-term average cost function. We derive the necessary and sufficient condition of the existence of a stationary and deterministic optimal policy that results in a bounded average cost in terms of the transmission reliabilities of the sensor-to-controller and controller-to-actuator channels. Also, we derive an easy-to-compute suboptimal policy, which notably reduces the average cost of the plant compared to a naive alternative-scheduling policy.


title: Discontinuous Galerkin methods for short pulse type equations via hodograph transformations

abstract: In the present paper, we consider the discontinuous Galerkin (DG) methods for solving short pulse (SP) type equations. The short pulse equation has been shown to be completely integrable, which admits the loop-soliton, cuspon-soliton solutions as well as smooth-soliton solutions. Through hodograph transformations, these nonclassical solutions can be profiled as the smooth solutions of the coupled dispersionless (CD) system or the sine-Gordon equation. Thus, DG methods can be developed for the CD system or the sine-Gordon equation to simulate the loop-soliton or cuspon-soliton solutions of the SP equation. The conservativeness or dissipation of the Hamiltonian or momentum for the semi-discrete DG schemes can be proved. Also we modify the above DG schemes and obtain an integration DG scheme. Theoretically the a-priori error estimates have been provided for the momentum conserved DG scheme and the integration DG scheme. We also propose the DG scheme and the integration DG scheme for the sine-Gordon equation, in case the SP equation can not be transformed to the CD system. All these DG schemes can be adopted to the generalized or modified SP type equations. Numerical experiments are provided to illustrate the optimal order of accuracy and capability of these DG schemes.


title: An Adaptive Approach for Anomaly Detector Selection and Fine-Tuning in Time Series

abstract: The anomaly detection of time series is a hotspot of time series data mining. The own characteristics of different anomaly detectors determine the abnormal data that they are good at. There is no detector can be optimizing in all types of anomalies. Moreover, it still has difficulties in industrial production due to problems such as a single detector can't be optimized at different time windows of the same time series. This paper proposes an adaptive model based on time series characteristics and selecting appropriate detector and run-time parameters for anomaly detection, which is called ATSDLN(Adaptive Time Series Detector Learning Network). We take the time series as the input of the model, and learn the time series representation through FCN. In order to realize the adaptive selection of detectors and run-time parameters according to the input time series, the outputs of FCN are the inputs of two sub-networks: the detector selection network and the run-time parameters selection network. In addition, the way that the variable layer width design of the parameter selection sub-network and the introduction of transfer learning make the model be with more expandability. Through experiments, it is found that ATSDLN can select appropriate anomaly detector and run-time parameters, and have strong expandability, which can quickly transfer. We investigate the performance of ATSDLN in public data sets, our methods outperform other methods in most cases with higher effect and better adaptation. We also show experimental results on public data sets to demonstrate how model structure and transfer learning affect the effectiveness.


title: Growing a Brain: Fine-Tuning by Increasing Model Capacity

abstract: CNNs have made an undeniable impact on computer vision through the ability to learn high-capacity models with large annotated training sets. One of their remarkable properties is the ability to transfer knowledge from a large source dataset to a (typically smaller) target dataset. This is usually accomplished through fine-tuning a fixed-size network on new target data. Indeed, virtually every contemporary visual recognition system makes use of fine-tuning to transfer knowledge from ImageNet. In this work, we analyze what components and parameters change during fine-tuning, and discover that increasing model capacity allows for more natural model adaptation through fine-tuning. By making an analogy to developmental learning, we demonstrate that "growing" a CNN with additional units, either by widening existing layers or deepening the overall network, significantly outperforms classic fine-tuning approaches. But in order to properly grow a network, we show that newly-added units must be appropriately normalized to allow for a pace of learning that is consistent with existing units. We empirically validate our approach on several benchmark datasets, producing state-of-the-art results.


title: Linear-semiorders and their incomparability graphs

abstract: A linear-interval order is the intersection of a linear order and an interval order. For this class of orders, several structural results have been known. This paper introduces a new subclass of linear-interval orders. We call a partial order a \emph{linear-semiorder} if it is the intersection of a linear order and a semiorder. We show a characterization and a polynomial-time recognition algorithm for linear-semiorders. We also prove that being a linear-semiorder is a comparability invariant, showing that incomparability graphs of linear-semiorders can be recognized in polynomial time.


title: Investigations and Improvement of Robustness of Reduced-Order Models of Reacting Flow

abstract: The impact of chemical reactions on the robustness and accuracy of projection-based Reduced-Order Models (ROMs) of fluid flows is investigated. Both Galerkin and Least-Squares Petrov Galerkin ROMs are shown to be less robust in reacting flows as compared to non-reacting flows. In particular, reacting flow ROMs show a strong sensitivity to \st{the} resolution and are often unstable. To identify the main underlying causes, a representative problem that contains the essential physics encountered in typical combustion dynamics problems is chosen. Comparisons with non-reacting solutions are used to assess the impact of reactions. Investigations are focused on three potential areas of significance: 1) preservation of conservation laws; 2) loss of dissipation; and 3) existence of unphysical local phenomena. Results indicate that conservation is relatively well-controlled and the global dissipation in the ROMs is actually larger than that in the underlying CFD solutions. Spurious local phenomena are, however, highly deleterious. Specifically, the steep temperature gradients that characterize combustion can lead to oscillations in local temperatures even in the absence of reactions. Representative calculations with physics-based temperature constraints verify that eliminating such excursions results in considerable improvement in both stability and future-state prediction capability.


title: Prioritized Guidance for Efficient Multi-Agent Reinforcement Learning Exploration

abstract: Exploration efficiency is a challenging problem in multi-agent reinforcement learning (MARL), as the policy learned by confederate MARL depends on the collaborative approach among multiple agents. Another important problem is the less informative reward restricts the learning speed of MARL compared with the informative label in supervised learning. In this work, we leverage on a novel communication method to guide MARL to accelerate exploration and propose a predictive network to forecast the reward of current state-action pair and use the guidance learned by the predictive network to modify the reward function. An improved prioritized experience replay is employed to better take advantage of the different knowledge learned by different agents which utilizes Time-difference (TD) error more effectively. Experimental results demonstrates that the proposed algorithm outperforms existing methods in cooperative multi-agent environments. We remark that this algorithm can be extended to supervised learning to speed up its training.


title: Game of Sloanes: Best known packings in complex projective space

abstract: It is often of interest to identify a given number of points in projective space such that the minimum distance between any two points is as large as possible. Such configurations yield representations of data that are optimally robust to noise and erasures. The minimum distance of an optimal configuration not only depends on the number of points and the dimension of the projective space, but also on whether the space is real or complex. For decades, Neil Sloane's online Table of Grassmannian Packings has been the go-to resource for putatively or provably optimal packings of points in real projective spaces. Using a variety of numerical algorithms, we have created a similar table for complex projective spaces. This paper surveys the relevant literature, explains some of the methods used to generate the table, presents some new putatively optimal packings, and invites the reader to competitively contribute improvements to this table.


title: Active galactic nuclei feedback at parsec scale

abstract: We perform simulations to study the effects of active galactic nuclei (AGNs) radiation and wind feedback on the properties of slowly rotating accretion flow at parsec scale. We find that when only radiative feedback is considered, outflows can be produced by the radiation pressure due to Thomson scattering. The mass flux of outflow is comparable to that of inflow. Although strong outflow is present, the luminosity of AGN can be easily super-Eddington. When wind feedback is also taken into account, the mass flux of outflow does not change much. Consequently, the luminosity of the central AGN can still be super-Eddington. However, observations show that the luminosity of most AGNs is sub-Eddington. Some other mechanisms are needed to reduce the AGNs luminosity. Although the mass outflow rate does not change much by wind feedback, other properties of outflow (density, temperature, velocity and kinetic power) can be significantly changed by wind feedback. In the presence of wind feedback, the density of outflow becomes significantly lower, temperature of outflow becomes significantly higher, velocity of outflow is increased by one order of magnitude, the kinetic power of outflow is increased by a factor of $40-100$.


title: Interval estimators for inequality measures using grouped data

abstract: Income inequality measures are often used as an indication of economic health. How to obtain reliable confidence intervals for these measures based on sampled data has been studied extensively in recent years. To preserve confidentiality, income data is often made available in summary form only (i.e. histograms, frequencies between quintiles, etc.). In this paper, we show that good coverage can be achieved for bootstrap and Wald-type intervals for quantile-based measures when only grouped (binned) data are available. These coverages are typically superior to those that we have been able to achieve for intervals for popular measures such as the Gini index in this grouped data setting. To facilitate the bootstrapping, we use the Generalized Lambda Distribution and also a linear interpolation approximation method to approximate the underlying density. The latter is possible when groups means are available. We also apply our methods to real data sets.


title: The props of quantum mechanics

abstract: We introduce a formalism that exploits the many-input many-output nature of nodes in quantum circuits. There is a diagrammatic and an algebraic version, the latter similar to the spinor formalism of general relativity. This allows us to work in truly basis independent ways, clarifying and simplifying many aspects of quantum state processing. The narrative is at times interrupted by antics of characters from quantum age fairy tales.