title: GraphFlow: Exploiting Conversation Flow with Graph Neural Networks for Conversational Machine Comprehension

abstract: Conversational machine comprehension (MC) has proven significantly more challenging compared to traditional MC since it requires better utilization of conversation history. However, most existing approaches do not effectively capture conversation history and thus have trouble handling questions involving coreference or ellipsis. Moreover, when reasoning over passage text, most of them simply treat it as a word sequence without exploring rich semantic relationships among words. In this paper, we first propose a simple yet effective graph structure learning technique to dynamically construct a question and conversation history aware context graph at each conversation turn. Then we propose a novel Recurrent Graph Neural Network, and based on that, we introduce a flow mechanism to model the temporal dependencies in a sequence of context graphs. The proposed GraphFlow model can effectively capture conversational flow in a dialog, and shows competitive performance compared to existing state-of-the-art methods on CoQA, QuAC and DoQA benchmarks. In addition, visualization experiments show that our proposed model can offer good interpretability for the reasoning process.


title: Novel Method for the Precise Determination of the QCD Running Coupling from Event Shapes Distributions in Electron-Positron Annihilation

abstract: We present a novel method for precisely determining the running QCD coupling constant $\alpha_s(Q^2)$ over a wide range of $Q^2$ from event shapes for electron-positron annihilation measured at a single annihilation energy $\sqrt{s}$. The renormalization scale $Q^2$ of the running coupling depends dynamically on the virtuality of the underlying quark and gluon subprocess and thus the specific kinematics of each event. The determination of the renormalization scale for event shape distributions is obtained by using the Principle of Maximum Conformality (PMC), a rigorous scale-setting method for gauge theories which satisfies all the requirements of Renormalization Group Invariance, including renormalization-scheme independence and consistency with Abelian theory in the $N_C \to 0$ limit. In this paper we apply the PMC to two classic event shapes measured in $e^+ e^-$ annihilation: the thrust ($T$) and C-parameter ($C$). The PMC renormalization scale depends differentially on the values of $T$ and $C$. The application of PMC scale-setting determines the running coupling $\alpha_s(Q^2)$ to high precision over a wide range of $Q^2$ from $10$ to $250~\rm{GeV}^2$ from measurements of the event shape distributions at the $Z^0$ peak. The extrapolation of the running coupling using pQCD evolution gives the value $\alpha_s(M^2_Z)=0.1185\pm0.0012$ from the thrust, and $\alpha_s(M^2_Z)=0.1193^{+0.0021}_{-0.0019}$ from the C-parameter in the $\bar{\rm{MS}}$ scheme. These determinations of $\alpha_s(M^2_Z)$ are consistent with the world average and are more precise than the values obtained from analyses of event shapes currently used in the world average. The highly-consistent results for the $T$ and $C$ event-shape distributions provide an additional verification of the applicability of the PMC to pQCD.


title: An Empirical Study of Batch Normalization and Group Normalization in Conditional Computation

abstract: Batch normalization has been widely used to improve optimization in deep neural networks. While the uncertainty in batch statistics can act as a regularizer, using these dataset statistics specific to the training set impairs generalization in certain tasks. Recently, alternative methods for normalizing feature activations in neural networks have been proposed. Among them, group normalization has been shown to yield similar, in some domains even superior performance to batch normalization. All these methods utilize a learned affine transformation after the normalization operation to increase representational power. Methods used in conditional computation define the parameters of these transformations as learnable functions of conditioning information. In this work, we study whether and where the conditional formulation of group normalization can improve generalization compared to conditional batch normalization. We evaluate performances on the tasks of visual question answering, few-shot learning, and conditional image generation.


title: M\"untz Sturm-Liouville Problems: Theory and Numerical Experiments

abstract: This paper presents two new classes of M\"untz functions which are called Jacobi-M\"untz functions of the first and second types. These newly generated functions satisfy in two self-adjoint fractional Sturm-Liouville problems and thus they have some spectral properties such as: orthogonality, completeness, three-term recurrence relations and so on. With respect to these functions two new orthogonal projections and their error bounds are derived. Also, two new M\"untz type quadrature rules are introduced. As two applications of these basis functions some fractional ordinary and partial differential equations are considered and numerical results are given.


title: Intrinsic Interleaving Distance for Merge Trees

abstract: Merge trees are a type of graph-based topological summary that tracks the evolution of connected components in the sublevel sets of scalar functions. They enjoy widespread applications in data analysis and scientific visualization. In this paper, we consider the problem of comparing two merge trees via the notion of interleaving distance in the metric space setting. We investigate various theoretical properties of such a metric. In particular, we show that the interleaving distance is intrinsic on the space of labeled merge trees and provide an algorithm to construct metric 1-centers for collections of labeled merge trees. We further prove that the intrinsic property of the interleaving distance also holds for the space of unlabeled merge trees. Our results are a first step toward performing statistics on graph-based topological summaries.


title: Normal completions of toric varieties over rank one valuation rings and completions of $\Gamma$-admissible fans

abstract: We show that any normal toric variety over a rank one valuation ring admits an equivariant open embedding in a normal toric variety which is proper over the valuation ring, after a base-change by a finite extension of valuation rings. If the value group $\Gamma$ is discrete or divisible then no base-change is needed. We give explicit examples which show that existing methods do not produce such normal equivariant completions. Our approach is combinatorial and proceeds by showing that $\Gamma$-admissible fans admit $\Gamma$-admissible completions. In order to show this we prove a combinatorial analog of noetherian reduction which we believe will be of independent interest.


title: Impact of new physics on $B+L$ violation at colliders

abstract: In the Standard Model, chiral electroweak anomalies predict nonperturbative interactions that violate baryon ($B$) plus lepton number ($L$). The potential observability of these processes at colliders has been amply discussed in the literature, mostly focusing on the impact of the accompanying boson emission, which contributes to the cross sections through an exponential function of the center-of-mass energy. We focus instead on the impact of exotic fermions charged under $SU(2)_L$, which not only can be emitted in these processes, but also affect the non-exponential contributions to the cross-sections. Estimating the latter using instanton techniques, we find sizable effects that suggest that if $B+L$-violating processes are ever seen at colliders, they may involve physics beyond the Standard Model.


title: Equilibrium states for non-uniformly hyperbolic systems: statistical properties and analyticity

abstract: We consider a wide family of non-uniformly expanding maps and hyperbolic H\"older continuous potentials. We prove that the unique equilibrium state associated to each element of this family is given by the eigenfunction of the transfer operator and the eigenmeasure of the dual operator (both having the spectral radius as eigenvalue). We show that the transfer operator has the spectral gap property in some space of H\"older continuous observables and from this we obtain an exponential decay of correlations and a central limit theorem for the equilibrium state. Moreover, we establish the analyticity with respect to the potential of the equilibrium state as well as that of other thermodynamic quantities. Furthermore, we derive similar results for the equilibrium state associated to a family of non-uniformly hyperbolic skew products and hyperbolic H\"older continuous potentials.


title: Experimental biases on the heliospheric contribution to the observed TeV cosmic ray anisotropy

abstract: The arrival direction distribution of cosmic ray particles observed on Earth is shaped by the cumulative effects of their galactic source locations and of trajectory bending in the turbulent interstellar magnetic field. Coherent magnetic structures are expected to disrupt particle trajectories and their observed distribution, as well. The heliosphere, the large magnetic bubble generated by the sweeping effect of solar wind on the local interstellar plasma, strongly affects the TeV cosmic ray particles detected on Earth. By unfolding the heliospheric influence on the observed anisotropy, it is possible to determine the pitch angle distribution of cosmic rays in the interstellar medium. This information makes it possible to study in detail the global diffusion properties of TeV cosmic rays in the Galaxy. However, observational blindness to key features of the cosmic ray arrival direction distribution may lead to biases in the determination of the heliospheric influence and of the interstellar CR distribution. Any inference of galactic TeV CR diffusion properties must carefully account for local propagation phenomena and observational limitations.


title: From Logistic Growth to Exponential Growth in a Population Dynamical Model

abstract: Dynamics among central sources (hubs) providing a resource and large number of components enjoying and contributing to this resource describes many real life situations. Modeling, controlling, and balancing this dynamics is a general problem that arises in many scientific disciplines. We analyze a stochastic dynamical system exhibiting this dynamics with a multiplicative noise. We show that this model can be solved exactly by passing to variables that describe the mass ratio between the components and the hub. We derive a deterministic equation for the average mass ratio. This equation describes logistic growth. We derive the full phase diagram of the model and identify three regimes by calculating the sample and moment Lyapunov exponent of the system. The first regime describes full balance between the non-hub components and the hub, in the second regime the entire resource is concentrated mainly in the hub, and in the third regime the resource is localized on a few non-hub components and the hub. Surprisingly, in the limit of large number of components the transition values do not depend on the amount of resource given by the hub. This model has interesting application in the context of analysis of porous media using Magnetic Resonance (MR) techniques.


title: Simultaneous Iris and Periocular Region Detection Using Coarse Annotations

abstract: In this work, we propose to detect the iris and periocular regions simultaneously using coarse annotations and two well-known object detectors: YOLOv2 and Faster R-CNN. We believe coarse annotations can be used in recognition systems based on the iris and periocular regions, given the much smaller engineering effort required to manually annotate the training images. We manually made coarse annotations of the iris and periocular regions (122K images from the visible (VIS) spectrum and 38K images from the near-infrared (NIR) spectrum). The iris annotations in the NIR databases were generated semi-automatically by first applying an iris segmentation CNN and then performing a manual inspection. These annotations were made for 11 well-known public databases (3 NIR and 8 VIS) designed for the iris-based recognition problem and are publicly available to the research community. Experimenting our proposal on these databases, we highlight two results. First, the Faster R-CNN + Feature Pyramid Network (FPN) model reported an Intersection over Union (IoU) higher than YOLOv2 (91.86% vs 85.30%). Second, the detection of the iris and periocular regions being performed simultaneously is as accurate as performed separately, but with a lower computational cost, i.e., two tasks were carried out at the cost of one.


title: Third- and fourth-order virial coefficients of harmonically trapped fermions in a semiclassical approximation

abstract: Using a leading-order semiclassical approximation, we calculate the third- and fourth-order virial coefficients of nonrelativistic spin-1/2 fermions in a harmonic trapping potential in arbitrary spatial dimensions, and as functions of temperature, trapping frequency and coupling strength. Our simple, analytic results for the interaction-induced changes $\Delta b_3$ and $\Delta b_4$ agree qualitatively, and in some regimes quantitatively, with previous numerical calculations for the unitary limit of three-dimensional Fermi gases.


title: Session-Based Hotel Recommendations: Challenges and Future Directions

abstract: In the year 2019, the Recommender Systems Challenge deals with a real-world task from the area of e-tourism for the first time, namely the recommendation of hotels in booking sessions. In this context, this article aims at identifying and investigating what we believe are important domain-specific challenges recommendation systems research in hotel search is facing, from both academic and industry perspectives. We focus on three main challenges, namely dealing with (1) multiple stakeholders and value-awareness in recommendations, (2) sparsity of user data and the extensive cold-start problem, and (3) dynamic input data and computational requirements. To this end, we review the state of the art toward solving these challenges and discuss shortcomings. We detail possible future directions and visions we contemplate for the further evolution of the field. This article should, therefore, serve two purposes: giving the interested reader an overview of current challenges in the field and inspiring new approaches for the ACM Recommender Systems Challenge 2019 and beyond.


title: Discovery of A Mg II Changing-look AGN and its Implications for a Unification Sequence of Changing-look AGNs

abstract: Changing-Look (CL) is a rare phenomenon of Active Galactic Nuclei (AGNs) that exhibit emerging or disappearing broad lines accompanied by continuum variations on astrophysically short timescales ($\lesssim$ 1 yr to a few decades). While previous studies have found Balmer-line (broad H$\alpha$ and/or H$\beta$) CL AGNs, the broad Mg II line is persistent even in dim states. No unambiguous Mg II CL AGN has been reported to date. We perform a systematic search of Mg II CL AGNs using multi-epoch spectra of a special population of Mg II-emitters (characterized by strong broad Mg II emission with little evidence for AGN from other normal indicators such as broad H$\alpha$ and H$\beta$ or blue power-law continua) from the Fourteenth Data Release of the Sloan Digital Sky Survey. We present the discovery of the first unambiguous case of an Mg II CL AGN, SDSS J152533.60+292012.1 (at redshift $z$ = 0.449), which is turning off within rest-frame 286 days. The dramatic diminishing of Mg II equivalent width (from 110 $\pm$ 26 Angstrom to being consistent with zero), together with little optical continuum variation ($\Delta V_{\rm max-min}$ $=$ 0.17 $\pm$ 0.05 mag) coevally over $\sim$ 10 years, rules out dust extinction or a tidal disruption event. Combined with previously known H$\beta$ CL AGNs, we construct a sequence that represents different temporal stages of CL AGNs. This CL sequence is best explained by the photoionization model of Guo et al. (2019). In addition, we present two candidate turn-on Mg II CL AGNs and a sample of 361 Mg II-emitters for future Mg II CL AGN searches.


title: Biased Average Position Estimates in Line and Bar Graphs: Underestimation, Overestimation, and Perceptual Pull

abstract: In visual depictions of data, position (i.e., the vertical height of a line or a bar) is believed to be the most precise way to encode information compared to other encodings (e.g., hue). Not only are other encodings less precise than position, but they can also be prone to systematic biases (e.g., color category boundaries can distort perceived differences between hues). By comparison, position's high level of precision may seem to protect it from such biases. In contrast, across three empirical studies, we show that while position may be a precise form of data encoding, it can also produce systematic biases in how values are visually encoded, at least for reports of average position across a short delay. In displays with a single line or a single set of bars, reports of average positions were significantly biased, such that line positions were underestimated and bar positions were overestimated. In displays with multiple data series (i.e., multiple lines and/or sets of bars), this systematic bias still persisted. We also observed an effect of "perceptual pull", where the average position estimate for each series was 'pulled' toward the other. These findings suggest that, although position may still be the most precise form of visual data encoding, it can also be systematically biased.


title: A worldsheet supersymmetric Newton-Cartan string

abstract: We construct a (locally) supersymmetric worldsheet action for a string in a non-relativistic Newton-Cartan background. We do this using a doubled string action, which describes the target space geometry in an $O(D,D)$ covariant manner using a doubled metric and doubled vielbeins. By adopting different parametrisations of these doubled background fields, we can describe both relativistic and non-relativistic geometries. We focus on the torsional Newton-Cartan geometry which can be obtained by null duality/reduction (such null duality is particularly simple for us to implement). The doubled action we use gives the Hamiltonian form of the supersymmetric Newton-Cartan string action automatically, from which we then obtain the equivalent Lagrangian. We extract geometric quantities of interest from the worldsheet couplings and write down the supersymmetry transformations. Our general results should apply to other non-relativistic backgrounds. We comment on the usefulness of the doubled approach as a tool for studying non-relativistic string theory.


title: Keplerian orbits through the Conley-Zehnder index

abstract: It was discovered by Gordon in 1977 that Keplerian ellipses in the plane are minimizers of the Lagrangian action and spectrally stable as periodic points of the associated Hamiltonian flow. The aim of this paper is to give a homotopy theoretical proof of these results through a self-contained, explicit and simple computation of the Conley-Zehnder index. The techniques developed in this paper can be used to investigate the higher dimensional case of Keplerian ellipses, where the classical variational proof no longer applies.


title: Complexity of magnetic fields on red dwarfs

abstract: Magnetic fields in cool stars can be investigated by measuring Zeeman line broadening and polarization in atomic and molecular lines. Similar to the Sun, these fields are complex and height-dependent. Many molecular lines dominating M-dwarf spectra (e.g., FeH, CaH, MgH, and TiO) are temperature -- and Zeeman -- sensitive and form at different atmospheric heights, which makes them excellent probes of magnetic fields on M dwarfs. Our goal is to analyze the complexity of magnetic fields in M dwarfs. We investigate how magnetic fields vary with the stellar temperature and how "surface" inhomogeneities are distributed in height -- the dimension that is usually neglected in stellar magnetic studies. We have determined effective temperatures of the photosphere and of magnetic features, magnetic field strengths and filling factors for nine M dwarfs (M1-M7). Our chi^2 analysis is based on a comparison of observed and synthetic intensity and circular polarization profiles. Stokes profiles were calculated by solving polarized radiative transfer equations. Properties of magnetic structures depend on the analyzed atomic or molecular species and their formation heights. Two types of magnetic features similar to those on the Sun have been found: a cooler (starspots) and a hotter (network) one. The magnetic field strength in both starspots and network is within 3 kG to 6 kG, on average it is 5 kG. These fields occupy a large fraction of M dwarf atmospheres at all heights, up to 100%. The plasma beta is less than one, implying highly magnetized stars. A combination of molecular and atomic species and a simultaneous analysis of intensity and circular polarization spectra have allowed us to better decipher the complexity of magnetic fields on M dwarfs, including their dependence on the atmospheric height. This work provides an opportunity to investigate a larger sample of M dwarfs and L-type brown dwarfs.


title: Learning in brain-computer interface control evidenced by joint decomposition of brain and behavior

abstract: Motor imagery-based brain-computer interfaces (BCIs) use an individuals ability to volitionally modulate localized brain activity as a therapy for motor dysfunction or to probe causal relations between brain activity and behavior. However, many individuals cannot learn to successfully modulate their brain activity, greatly limiting the efficacy of BCI for therapy and for basic scientific inquiry. Previous research suggests that coherent activity across diverse cognitive systems is a hallmark of individuals who can successfully learn to control the BCI. However, little is known about how these distributed networks interact through time to support learning. Here, we address this gap in knowledge by constructing and applying a multimodal network approach to decipher brain-behavior relations in motor imagery-based brain-computer interface learning using MEG. Specifically, we employ a minimally constrained matrix decomposition method (non-negative matrix factorization) to simultaneously identify regularized, covarying subgraphs of functional connectivity, to assess their similarity to task performance, and to detect their time-varying expression. Individuals also displayed marked variation in the spatial properties of subgraphs such as the connectivity between the frontal lobe and the rest of the brain, and in the temporal properties of subgraphs such as the stage of learning at which they reached maximum expression. From these observations, we posit a conceptual model in which certain subgraphs support learning by modulating brain activity in regions important for sustaining attention. To test this model, we use tools that stipulate regional dynamics on a networked system (network control theory), and find that good learners display a single subgraph whose temporal expression tracked performance and whose architecture supports easy modulation of brain regions important for attention.


title: GRB lensing parallax: Closing primordial black hole dark matter mass window

abstract: The primordial black hole (PBH) comprising full dark matter (DM) abundance is currently allowed if its mass lies between $10^{-16}M_{\odot} \lesssim M \lesssim 10^{-11} M_{\odot}$. This lightest mass range is hard to be probed by ongoing gravitational lensing observations. In this paper, we advocate that an old idea of the lensing parallax of Gamma-ray bursts (GRBs), observed simultaneously by spatially separated detectors, can probe the unconstrained mass range; and that of nearby stars can probe a heavier mass range. In addition to various good properties of GRBs, astrophysical separations achievable around us --- $r_\oplus \text{--}$ AU --- is just large enough to resolve the GRB lensing by lightest PBH DM.


title: Chiral active matter: microscopic `torque dipoles' have more than one hydrodynamic description

abstract: Many biological systems, such as bacterial suspensions and actomyosin networks, form polar liquid crystals. These systems are `active' or far-from-equilibrium, due to local forcing of the solvent by the constituent particles. In many cases the source of activity is chiral; since forcing is internally generated, some sort of `torque dipole' is then present locally. But it is not obvious how `torque dipoles' should be encoded in the hydrodynamic equations that describe the system at continuum level: different authors have arrived at contradictory conclusions on this issue. In this work, we resolve the paradox by presenting a careful derivation, from linear irreversible thermodynamics, of the general equations of motion of a single-component chiral active fluid with spin degrees of freedom. We find that there is no unique hydrodynamic description for such a fluid in the presence of torque dipoles of a given strength. Instead, at least three different hydrodynamic descriptions emerge, depending on whether we decompose each torque dipole as two point torques, two force pairs, or one point torque and one force pair -- where point torques create internal angular momenta of the chiral bodies (spin), whereas force pairs impart centre of mass motion that contributes to fluid velocity. By considering a general expansion of the Onsager coefficients, we also derive a new shear-elongation parameter and cross-coupling viscosity, which can lead to unpredicted phenomena even in passive polar liquid crystals. Finally, elimination of the angular variables gives an effective polar hydrodynamics with renormalized active stresses, viscosities and kinetic coefficients. Remarkably, this can include a direct contribution of chiral activity to the equation of motion for the polar order parameter, which survives even in `dry' active systems where the fluid velocity is set to zero.

categories: cond-mat.soft cond-mat.stat-mech physics.bio-ph physics.flu-dyn


title: Machine Learning at the Network Edge: A Survey

abstract: Resource-constrained IoT devices, such as sensors and actuators, have become ubiquitous in recent years. This has led to the generation of large quantities of data in real-time, which is an appealing target for AI systems. However, deploying machine learning models on such end-devices is nearly impossible. A typical solution involves offloading data to external computing systems (such as cloud servers) for further processing but this worsens latency, leads to increased communication costs, and adds to privacy concerns. To address this issue, efforts have been made to place additional computing devices at the edge of the network, i.e close to the IoT devices where the data is generated. Deploying machine learning systems on such edge computing devices alleviates the above issues by allowing computations to be performed close to the data sources. This survey describes major research efforts where machine learning systems have been deployed at the edge of computer networks, focusing on the operational aspects including compression techniques, tools, frameworks, and hardware used in successful applications of intelligent edge systems.


title: Direct and inverse results on restricted signed sumsets in integers

abstract: Let $G$ be an additive abelian group. Let $A=\{a_{0}, a_{1},\ldots, a_{k-1}\}$ be a nonempty finite subset of $G$. For a positive integer $h$ satisfying $1\leq h\leq k$, we let \[h\hat{}_{\underline{+}}A:=\{\Sigma_{i=0}^{k-1}\lambda_{i} a_{i}: (\lambda_{0},\lambda_{1}, \ldots, \lambda_{k-1}) \in \{-1,0,1\}^{k},~\Sigma_{i=0}^{k-1}|\lambda_{i}|=h \},\] be the restricted signed sumset of $A$. The direct problem for the restricted signed sumset $h\hat{}_{\underline{+}}A$ is to find the minimum number of elements in $h\hat{}_{\underline{+}}A$ in terms of $|A|$. The inverse problem for $h\hat{}_{\underline{+}}A$ is to determine the structure of the finite set $A$ for which $|h\hat{}_{\underline{+}}A|$ is minimal. In this article, we solve some cases of both direct and inverse problems for $h\hat{}_{\underline{+}}A$, when $A$ is a finite set of integers. In this connection, we also pose some questions as conjectures in the remaining cases.


title: A refinement of Khovanov homology

abstract: We refine Khovanov homology in the presence of an involution on the link. This refinement takes the form of a triply-graded theory, arising from a pair of filtrations. We focus primarily on strongly invertible knots and show, for instance, that this refinement is able to detect mutation.


title: Cyclic sieving, skew Macdonald polynomials and Schur positivity

abstract: When $\lambda$ is a partition, the specialized non-symmetric Macdonald polynomial $E_{\lambda}(x;q;0)$ is symmetric and related to a modified Hall--Littlewood polynomial. We show that whenever all parts of the integer partition $\lambda$ is a multiple of $n$, the underlying set of fillings exhibit the cyclic sieving phenomenon (CSP) under a cyclic shift of the columns. The corresponding CSP polynomial is given by $E_{\lambda}(x;q;0)$. In addition, we prove a refined cyclic sieving phenomenon where the content of the fillings is fixed. This refinement is closely related to an earlier result by B.~Rhoades. We also introduce a skew version of $E_{\lambda}(x;q;0)$. We show that these are symmetric and Schur-positive via a variant of the Robinson--Schenstedt--Knuth correspondence and we also describe crystal raising- and lowering operators for the underlying fillings. Moreover, we show that the skew specialized non-symmetric Macdonald polynomials are in some cases vertical-strip LLT polynomials. As a consequence, we get a combinatorial Schur expansion of a new family of LLT polynomials.


title: Dynamic Information Design with Diminishing Sensitivity Over News

abstract: A Bayesian agent experiences gain-loss utility each period over changes in belief about future consumption ("news utility"), with diminishing sensitivity over the magnitude of news. Diminishing sensitivity induces a preference over news skewness: gradual bad news, one-shot good news is worse than one-shot resolution, which is in turn worse than gradual good news, one-shot bad news. So, the agent's preference between gradual information and one-shot resolution can depend on his consumption ranking of different states. In a dynamic cheap-talk framework where a benevolent sender communicates the state over multiple periods, the babbling equilibrium is essentially unique without loss aversion. More loss-averse agents may enjoy higher news utility in equilibrium, contrary to the commitment case. We characterize the family of gradual good news equilibria that exist with high enough loss aversion, and find the sender conveys progressively larger pieces of good news. We discuss applications to media competition and game shows.


title: Why Does My Model Fail? Contrastive Local Explanations for Retail Forecasting

abstract: In various business settings, there is an interest in using more complex machine learning techniques for sales forecasting. It is difficult to convince analysts, along with their superiors, to adopt these techniques since the models are considered to be "black boxes," even if they perform better than current models in use. We examine the impact of contrastive explanations about large errors on users' attitudes towards a "black-box'" model. We propose an algorithm, Monte Carlo Bounds for Reasonable Predictions. Given a large error, MC-BRP determines (1) feature values that would result in a reasonable prediction, and (2) general trends between each feature and the target, both based on Monte Carlo simulations. We evaluate on a real dataset with real users by conducting a user study with 75 participants to determine if explanations generated by MC-BRP help users understand why a prediction results in a large error, and if this promotes trust in an automatically-learned model. Our study shows that users are able to answer objective questions about the model's predictions with overall 81.1% accuracy when provided with these contrastive explanations. We show that users who saw MC-BRP explanations understand why the model makes large errors in predictions significantly more than users in the control group. We also conduct an in-depth analysis on the difference in attitudes between Practitioners and Researchers, and confirm that our results hold when conditioning on the users' background.


title: Learning Effective Embeddings From Crowdsourced Labels: An Educational Case Study

abstract: Learning representation has been proven to be helpful in numerous machine learning tasks. The success of the majority of existing representation learning approaches often requires a large amount of consistent and noise-free labels. However, labels are not accessible in many real-world scenarios and they are usually annotated by the crowds. In practice, the crowdsourced labels are usually inconsistent among crowd workers given their diverse expertise and the number of crowdsourced labels is very limited. Thus, directly adopting crowdsourced labels for existing representation learning algorithms is inappropriate and suboptimal. In this paper, we investigate the above problem and propose a novel framework of \textbf{R}epresentation \textbf{L}earning with crowdsourced \textbf{L}abels, i.e., "RLL", which learns representation of data with crowdsourced labels by jointly and coherently solving the challenges introduced by limited and inconsistent labels. The proposed representation learning framework is evaluated in two real-world education applications. The experimental results demonstrate the benefits of our approach on learning representation from limited labeled data from the crowds, and show RLL is able to outperform state-of-the-art baselines. Moreover, detailed experiments are conducted on RLL to fully understand its key components and the corresponding performance.


title: explAIner: A Visual Analytics Framework for Interactive and Explainable Machine Learning

abstract: We propose a framework for interactive and explainable machine learning that enables users to (1) understand machine learning models; (2) diagnose model limitations using different explainable AI methods; as well as (3) refine and optimize the models. Our framework combines an iterative XAI pipeline with eight global monitoring and steering mechanisms, including quality monitoring, provenance tracking, model comparison, and trust building. To operationalize the framework, we present explAIner, a visual analytics system for interactive and explainable machine learning that instantiates all phases of the suggested pipeline within the commonly used TensorBoard environment. We performed a user-study with nine participants across different expertise levels to examine their perception of our workflow and to collect suggestions to fill the gap between our system and framework. The evaluation confirms that our tightly integrated system leads to an informed machine learning process while disclosing opportunities for further extensions.


title: Advanced material system for the design of an intermediate band solar cell: type-II CdTe quantum dots in a ZnCdSe matrix

abstract: We explore CdTe fractional monolayer quantum dots (QDs) in a ZnCdSe host matrix for potential application in an intermediate band solar cell device. Careful consideration has been taken during the initiation of the growth process of QDs by migration enhanced epitaxy, in order to avoid the formation of undesirable interfacial layers that can form due to the lack of common anion between the two materials. A superlattice structure of 100 periods of alternating QD and spacer layers is analyzed by high resolution X-ray diffraction (XRD) and photoluminescent (PL) spectroscopy. Simple arguments are used following continuum elastic theory to deduce the size of the dots and the strain within the superlattice from XRD data. This is further verified using PL and used in the energy calculations that yield the values of the intermediate band energy. The results suggest that the optimized materials are highly suitable for these high efficiency solar cells.


title: A Model of Random Industrial SAT

abstract: One of the most studied models of SAT is random SAT. In this model, instances are composed from clauses chosen uniformly randomly and independently of each other. This model may be unsatisfactory in that it fails to describe various features of SAT instances, arising in real-world applications. Various modifications have been suggested to define models of industrial SAT. Here, we focus mainly on the aspect of community structure. Namely, here the set of variables consists of a number of disjoint communities, and clauses tend to consist of variables from the same community. Thus, we suggest a model of random industrial SAT, in which the central generalization with respect to random SAT is the additional community structure. There has been a lot of work on the satisfiability threshold of random $k$-SAT, starting with the calculation of the threshold of $2$-SAT, up to the recent result that the threshold exists for sufficiently large $k$. In this paper, we endeavor to study the satisfiability threshold for the proposed model of random industrial SAT. Our main result is that the threshold in this model tends to be smaller than its counterpart for random SAT. Moreover, under some conditions, this threshold even vanishes.


title: An Optimal Linear Dynamic Detection Method for Replay Attack in Cyber-Physical Systems

abstract: The problem of detecting replay attack to the linear stochastic system with Kalman filer state estimator and LQG controller is addressed. To this end, a dynamic attack detector method is proposed which is coupled with the dynamics of the system. While preserving stability of the main system, conditions on parameters of the attack detector dynamics are obtained such that the attack can be revealed by destabilization of a residual trajectory which is the difference between the estimated and measured output of the system. Using this method, system operator can adjust the detection rate based on the proposed scheme by changing the design parameters. Nevertheless, since the exogenous attack detector signal affects the performance of the closed loop control system, we propose an optimization problem to determine such a detector with minimum loss effect. In the simulation results, the proposed dynamical attack detector approach is compared with the well-known additive white noise watermarking method and the results confirm the superiority of the new scheme.


title: Triple product p-adic L-functions for Shimura curves over totally real number fields

abstract: Let F be a totally real number field. Using a recent geometric approach developed by Andreatta and Iovita we construct several variables p-adic families of finite slope quaternionic automorphic forms over F. It is achieved by interpolating the modular sheaves defined over some auxiliary unitary Shimura curves. Secondly, we attach p-adic L-functions to triples of ordinary p-adic families of quaternionic automorphic eigenforms. This is done by relating trilinear periods to some trilinear products over unitary Shimura curves which can be interpolated adapting the work of Liu-Zhang-Zhang to our families.


title: Hypergraph based Berge hypergraphs

abstract: Fix a hypergraph $\mathcal{F}$. A hypergraph $\mathcal{H}$ is called a {\it Berge copy of $\mathcal{F}$} or {\it Berge-$\mathcal{F}$} if we can choose a subset of each hyperedge of $\mathcal{H}$ to obtain a copy of $\mathcal{F}$. A hypergraph $\mathcal{H}$ is {\it Berge-$\mathcal{F}$-free} if it does not contain a subhypergraph which is Berge copy of $\mathcal{F}$. This is a generalization of the usual, graph based Berge hypergraphs, where $\mathcal{F}$ is a graph. In this paper, we study extremal properties of hypergraph based Berge hypergraphs and generalize several results from the graph based setting. In particular, we show that for any $r$-uniform hypregraph $\mathcal{F}$, the sum of the sizes of the hyperedges of a (not necessarily uniform) Berge-$\mathcal{F}$-free hypergraph $\mathcal{H}$ on $n$ vertices is $o(n^r)$ when all the hyperedges of $\mathcal{H}$ are large enough. We also give a connection between hypergraph based Berge hypergraphs and generalized hypergraph Tur\'an problems.


title: Aquarium: Cassiopea and Alewife Languages

abstract: This technical report describes two of the domain specific languages used in the Aquarium kernel code synthesis project. It presents the language cores in terms of abstract syntax. Cassiopea is a machine description language for describing the semantics of processor instruction sets. Alewife is a specification language that can be used to write machine-independent specifications for assembly-level instruction blocks. An Alewife specification can be used to verify and synthesize code for any machine described in Cassiopea, given a machine-specific translation for abstractions used in the specification. This article does not include an introduction to either the Aquarium system or the use of the languages. In addition to this version of the article being a draft, the Aquarium project and the languages are works in progress. This article cannot currently be considered either final or complete.


title: Applications of Kinetic Tools to Inverse Transport Problems

abstract: We show that the inverse problems for a class of kinetic equations can be solved by classical tools in PDE analysis including energy estimates and the celebrated averaging lemma. Using these tools, we give a unified framework for the reconstruction of the absorption coefficient for transport equations in the subcritical and critical regimes. Moreover, we apply this framework to obtain, to the best of our knowledge, the first result in a nonlinear setting. We also extend the result of recovering the scattering coefficient in [14] from 3D to 2D convex domains.


title: Scintillation Arc Brightness and Electron Density for an Analytical Noodle Model

abstract: We show that narrow filaments or sheets of over- or under-dense plasma, or "noodles," with fluctuations of scattering phase of less than a radian, can form the scintillation arcs seen for many pulsars. The required local fluctuations of electron density are indefinitely small. We assume a cosine profile for the electron column and find the scattered field by analytic Kirchhoff integration. For a large electron column, corresponding to large amplitude of phase variation, the stationary-phase approximation is accurate; we call this regime "ray optics". For smaller-amplitude phase variation, the stationary-phase approximation is inaccurate or inapplicable; we call this regime "wave optics". We show that scattering is most efficient when the width of the strip equals that of one pair of Fresnel zones, and in the wave-optics regime. We show that the resolution of present observations is about 100 Fresnel zones on the scattering screen. Incoherent superposition of strips within a resolution element tends to increase the scattered field. We find that observations match a single noodle per resolution element with phase of up to 12 radians; or many noodles per resolution element with arbitrarily small phase variation each, for net phase of less than a radian. Observations suggest a minimum radius for noodles of about 650 km, comparable to the ion inertial scale or the ion cyclotron radius in the scattering plasma.


title: Adversarial Robustness Curves

abstract: The existence of adversarial examples has led to considerable uncertainty regarding the trust one can justifiably put in predictions produced by automated systems. This uncertainty has, in turn, lead to considerable research effort in understanding adversarial robustness. In this work, we take first steps towards separating robustness analysis from the choice of robustness threshold and norm. We propose robustness curves as a more general view of the robustness behavior of a model and investigate under which circumstances they can qualitatively depend on the chosen norm.


title: Heavy-hole band splitting observed in mobility spectrum of p-type InAs grown on GaAs substrate

abstract: High quality berylium doped InAs layer grown by MBE on GaAs substrate has been examined via magnetotransport measurements and high resolution quantitative mobility spectrum analysis in the range from 5 to 300 K and up to 15 T magnetic field. The layer homogenity and dopant concentration has been proofed via HR-SIMS. The results shew four channel conductivity and essential splitting of the most populated holelike channel below 55 K. The multilayer model concluded from the QMSA results has been compared with nextnano simulation.


title: On equations and first-order theory of one-relator monoids

abstract: We investigate systems of equations and the first-order theory of one-relator monoids. We describe a family $\mathcal{F}$ of one-relator monoids of the form $\langle A\mid w=1\rangle$ where for each monoid $M$ in $\mathcal{F}$, the longstanding open problem of decidability of word equations with length constraints reduces to the Diophantine problem (i.e.\ decidability of systems of equations) in $M$. We achieve this result by finding an interpretation in $M$ of a free monoid, using only systems of equations together with length relations. It follows that each monoid in $\mathcal{F}$ has undecidable positive AE-theory, hence in particular it has undecidable first-order theory. The family $\mathcal{F}$ includes many one-relator monoids with torsion $\langle A\mid w^n = 1\rangle$ ($n>1$). In contrast, all one-relator groups with torsion are hyperbolic, and all hyperbolic groups are known to have decidable Diophantine problem. We further describe a different class of one-relator monoids with decidable Diophantine problem.


title: Testing for Externalities in Network Formation Using Simulation

abstract: We discuss a simplified version of the testing problem considered by Pelican and Graham (2019): testing for interdependencies in preferences over links among N (possibly heterogeneous) agents in a network. We describe an exact test which conditions on a sufficient statistic for the nuisance parameter characterizing any agent-level heterogeneity. Employing an algorithm due to Blitzstein and Diaconis (2011), we show how to simulate the null distribution of the test statistic in order to estimate critical values and/or p-values. We illustrate our methods using the Nyakatoke risk-sharing network. We find that the transitivity of the Nyakatoke network far exceeds what can be explained by degree heterogeneity across households alone.


title: The Rate of Short Duration Gamma-Ray Bursts in the local Universe

abstract: The binary neutron star merger responsible for the gravitational wave event, GW170817, strengthened the merger association with short-duration gamma-ray bursts (SGRBs) following the detection of the SGRB counterpart, GRB 170817A. Here we consider the constraints on a population of local gamma-ray bursts with moderately short duration ($T_{90}<4 \textrm{ s}$) and within $d < 200 \textrm{ Mpc}$, that may have originated from similar compact binary mergers. Using well localised gamma-ray bursts from $\sim14.5 \textrm{ years}$ of Swift/Burst Alert Telescope monitoring, we find no events with high likelihood of being in this distance range, and place an upper limit for the all-sky rate of such events of $<4\ \textrm{y}^{-1}$. For Fermi/Gamma-ray Burst Monitor (GBM) and CGRO/Burst And Transient Source Experiment (BATSE) detected bursts, where the localisation has considerably larger uncertainties, we cross-correlated with 2MASS Redshift Survey galaxies at $d<100 \textrm{ Mpc}$, obtaining a weaker constraint of $< 12\ \textrm{y}^{-1}$. A separate correlation search from the GBM and BATSE bursts for giant flares originating from soft gamma-ray repeaters in nearby galaxies ($d < 11 \textrm{ Mpc} $ ) yields an upper limit of $<3\ \textrm{y}^{-1}$.


title: Perspective: Heusler interfaces -- opportunities beyond spintronics?

abstract: Heusler compounds, in both cubic and hexagonal polymorphs, exhibit a remarkable range of electronic, magnetic, elastic, and topological properties, rivaling that of the transition metal oxides. To date, research on these quantum materials has focused primarily on bulk magnetic and thermoelectric properties or on applications in spintronics. More broadly, however, Heuslers provide a platform for discovery and manipulation of emergent properties at well-defined crystalline interfaces. Here, motivated by advances in the epitaxial growth of layered Heusler heterostructures, I present a vision for Heusler interfaces, focusing on the frontiers and challenges that lie beyond spintronics. The ability to grow these materials epitaxially on technologically important semiconductor substrates, such as GaAs, Ge, and Si, provides a direct path for their integration with modern electronics. Further advances will require new methods to control the stoichiometry and defects to "electronic grade" quality, and to control the interface abruptness and ordering at the atomic scale.


title: OCT Fingerprints: Resilience to Presentation Attacks

abstract: Optical coherent tomography (OCT) fingerprint technology provides rich depth information, including internal fingerprint (papillary junction) and sweat (eccrine) glands, in addition to imaging any fake layers (presentation attacks) placed over finger skin. Unlike 2D surface fingerprint scans, additional depth information provided by the cross-sectional OCT depth profile scans are purported to thwart fingerprint presentation attacks. We develop and evaluate a presentation attack detector (PAD) based on deep convolutional neural network (CNN). Input data to CNN are local patches extracted from the cross-sectional OCT depth profile scans captured using THORLabs Telesto series spectral-domain fingerprint reader. The proposed approach achieves a TDR of 99.73% @ FDR of 0.2% on a database of 3,413 bonafide and 357 PA OCT scans, fabricated using 8 different PA materials. By employing a visualization technique, known as CNN-Fixations, we are able to identify the regions in the OCT scan patches that are crucial for fingerprint PAD detection.


title: Solvable lattice models for metals with Z2 topological order

abstract: We present quantum dimer models in two dimensions which realize metallic ground states with Z2 topological order. Our models are generalizations of a dimer model introduced in [PNAS 112,9552-9557 (2015)] to provide an effective description of unconventional metallic states in hole-doped Mott insulators. We construct exact ground state wave functions in a specific parameter regime and show that the ground state realizes a fractionalized Fermi liquid. Due to the presence of Z2 topological order the Luttinger count is modified and the volume enclosed by the Fermi surface is proportional to the density of doped holes away from half filling. We also comment on possible applications to magic-angle twisted bilayer graphene.


title: Evaluation of the Implementation of an Abstract Interpretation Algorithm using Tabled CLP

abstract: CiaoPP is an analyzer and optimizer for logic programs, part of the Ciao Prolog system. It includes PLAI, a fixpoint algorithm for the abstract interpretation of logic programs which we adapt to use tabled constraint logic programming. In this adaptation, the tabling engine drives the fixpoint computation, while the constraint solver handles the LUB of the abstract substitutions of different clauses. That simplifies the code and improves performance, since termination, dependencies, and some crucial operations (e.g., branch switching and resumption) are directly handled by the tabling engine. Determining whether the fixpoint has been reached uses semantic equivalence, which can decide that two syntactically different abstract substitutions represent the same element in the abstract domain. Therefore, the tabling analyzer can reuse answers in more cases than an analyzer using syntactical equality. This helps achieve better performance, even taking into account the additional cost associated to these checks. Our implementation is based on the TCLP framework available in Ciao Prolog and is one-third the size of the initial fixpoint implementation in CiaoPP. Its performance has been evaluated by analyzing several programs using different abstract domains. This paper is under consideration for publication in Theory and Practice of Logic Programming (TPLP).


title: Conditional independence testing: a predictive perspective

abstract: Conditional independence testing is a key problem required by many machine learning and statistics tools. In particular, it is one way of evaluating the usefulness of some features on a supervised prediction problem. We propose a novel conditional independence test in a predictive setting, and show that it achieves better power than competing approaches in several settings. Our approach consists in deriving a p-value using a permutation test where the predictive power using the unpermuted dataset is compared with the predictive power of using dataset where the feature(s) of interest are permuted. We conclude that the method achives sensible results on simulated and real datasets.


title: On (unitary) perfect polynomials over $\mathbb{F}_2$ with only Mersenne primes as odd divisors

abstract: The only (unitary) perfect polynomials over $\mathbb{F}_2$ that are products of $x$, $x+1$ and Mersenne primes are precisely the nine (resp. nine "classes") known ones. This follows from a new result about the factorization of $M^{2h+1} +1$, for a Mersenne prime $M$ and for a positive integer $h$. Other consequences of such a factorization are new results about odd perfect polynomials.


title: Single-timescale distributed GNE seeking for aggregative games over networks via forward-backward operator splitting

abstract: We consider aggregative games with affine coupling constraints, where agents have partial information on the aggregate value and can only communicate with neighbouring agents. We propose a single-layer distributed algorithm that reaches a variational generalized Nash equilibrium, under constant step sizes. The algorithm works on a single timescale, i.e., does not require multiple communication rounds between agents before updating their action. The convergence proof leverages an invariance property of the aggregate estimates and relies on a forward-backward splitting for two preconditioned operators and their restricted (strong) monotonicity properties on the consensus subspace.


title: Strongly correlated quantum droplets in quasi-1D dipolar Bose gas

abstract: We exploit a few- to many-body approach to study strongly interacting dipolar bosons in the quasi-one-dimensional system. The dipoles attract each other while the short range interactions are repulsive. Solving numerically the multi-atom Schr\"{o}dinger equation, we discover that such systems can exhibit not only the well known bright soliton solutions but also novel quantum droplets for a strongly coupled case. For larger systems, basing on microscopic properties of the found few-body solution, we propose a new equation for a density amplitude of atoms. It accounts for fermionization for strongly repelling bosons by incorporating the Lieb-Liniger energy in a local density approximation and approaches the standard Gross-Pitaevskii equation (GPE) in the weakly interacting limit. Not only does such a framework provide an alternative mechanism of the droplet stability, but it also introduces means to further analyze this previously unexplored quantum phase. In the limiting strong repulsion case, yet another simple multi-atom model is proposed. We stress that the celebrated Lee-Huang-Yang term in the GPE is not applicable in this case.


title: General rogue waves in the Boussinesq equation

abstract: We derive general rogue wave solutions of arbitrary orders in the Boussinesq equation by the bilinear Kadomtsev-Petviashvili (KP) reduction method. These rogue solutions are given as Gram determinants with $2N-2$ free irreducible real parameters, where $N$ is the order of the rogue wave. Tuning these free parameters, rogue waves of various patterns are obtained, many of which have not been seen before. Compared to rogue waves in other integrable equations, a new feature of rogue waves in the Boussinesq equation is that the rogue wave of maximum amplitude at each order is generally asymmetric in space. On the technical aspect, our contribution to the bilinear KP-reduction method for rogue waves is a new judicious choice of differential operators in the procedure, which drastically simplifies the dimension reduction calculation as well as the analytical expressions of rogue wave solutions.


title: Towards an in situ, full-power gauge of the focal-volume intensity of petawatt-class lasers

abstract: About 50 years ago, Sarachick and Schappert [Phys. Rev. D. 1, 2738 (1970)] showed that relativistic Thomson scattering leads to wavelength shifts that are proportional to the laser intensity. About 28 years later Chen et al. [Nature 396, 653 (1998)] used these shifts to estimate their laser intensity near 10^18 W/cm^2. More recently there have been several theoretical studies aimed at exploiting nonlinear Thomson scattering as a tool for direct measurement of intensities well into the relativistic regime. We present the first quantitative study of this approach for intensities between 10^18 and 10^19 W/cm^2. We show that the spectral shifts are in reasonable agreement with estimates of the peak intensity extracted from images of the focal area obtained at reduced power. Finally, we discuss the viability of the approach, its range of usefulness and how it might be extended to gauge intensities well in excess of 10^19 W/cm^2.


title: Few-Shot Meta-Denoising

abstract: We study the problem of few-shot learning-based denoising where the training set contains just a handful of clean and noisy samples. A solution to mitigate the small training set issue is to pre-train a denoising model with small training sets containing pairs of clean and synthesized noisy signals, produced from empirical noise priors, and fine-tune on the available small training set. While such transfer learning seems effective, it may not generalize well because of the limited amount of training data. In this work, we propose a new meta-learning training approach for few-shot learning-based denoising problems. Our model is meta-trained using known synthetic noise models, and then fine-tuned with the small training set, with the real noise, as a few-shot learning task. Meta-learning from small training sets of synthetically generated data during meta-training enables us to not only generate an infinite number of training tasks, but also train a model to learn with small training sets -- both advantages have the potential to improve the generalisation of the denoising model. Our approach is empirically shown to produce more accurate denoising results than supervised learning and transfer learning in three denoising evaluations for images and 1-D signals. Interestingly, our study provides strong indications that meta-learning has the potential to become the main learning algorithm for denoising.


title: Domain-Independent Cost-Optimal Planning in ASP

abstract: We investigate the problem of cost-optimal planning in ASP. Current ASP planners can be trivially extended to a cost-optimal one by adding weak constraints, but only for a given makespan (number of steps). It is desirable to have a planner that guarantees global optimality. In this paper, we present two approaches to addressing this problem. First, we show how to engineer a cost-optimal planner composed of two ASP programs running in parallel. Using lessons learned from this, we then develop an entirely new approach to cost-optimal planning, stepless planning, which is completely free of makespan. Experiments to compare the two approaches with the only known cost-optimal planner in SAT reveal good potentials for stepless planning in ASP. The paper is under consideration for acceptance in TPLP.


title: A Structural Average of Labeled Merge Trees for Uncertainty Visualization

abstract: Physical phenomena in science and engineering are frequently modeled using scalar fields. In scalar field topology, graph-based topological descriptors such as merge trees, contour trees, and Reeb graphs are commonly used to characterize topological changes in the (sub)level sets of scalar fields. One of the biggest challenges and opportunities to advance topology-based visualization is to understand and incorporate uncertainty into such topological descriptors to effectively reason about their underlying data. In this paper, we study a structural average of a set of labeled merge trees and use it to encode uncertainty in data. Specifically, we compute a 1-center tree that minimizes its maximum distance to any other tree in the set under a well-defined metric called the interleaving distance. We provide heuristic strategies that compute structural averages of merge trees whose labels do not fully agree. We further provide an interactive visualization system that resembles a numerical calculator that takes as input a set of merge trees and outputs a tree as their structural average. We also highlight structural similarities between the input and the average and incorporate uncertainty information for visual exploration. We develop a novel measure of uncertainty, referred to as consistency, via a metric-space view of the input trees. Finally, we demonstrate an application of our framework through merge trees that arise from ensembles of scalar fields. Our work is the first to employ interleaving distances and consistency to study a global, mathematically rigorous, structural average of merge trees in the context of uncertainty visualization.


title: 3D Virtual Garment Modeling from RGB Images

abstract: We present a novel approach that constructs 3D virtual garment models from photos. Unlike previous methods that require photos of a garment on a human model or a mannequin, our approach can work with various states of the garment: on a model, on a mannequin, or on a flat surface. To construct a complete 3D virtual model, our approach only requires two images as input, one front view and one back view. We first apply a multi-task learning network called JFNet that jointly predicts fashion landmarks and parses a garment image into semantic parts. The predicted landmarks are used for estimating sizing information of the garment. Then, a template garment mesh is deformed based on the sizing information to generate the final 3D model. The semantic parts are utilized for extracting color textures from input images. The results of our approach can be used in various Virtual Reality and Mixed Reality applications.


title: Nonlinear scalar perturbations of extremal Reissner-Nordstr\"{o}m spacetimes

abstract: We present the first rigorous study of nonlinear wave equations on extremal black hole spacetimes without any symmetry assumptions on the solution. Specifically, we prove global existence with asymptotic blow-up for solutions to nonlinear wave equations satisfying the null condition on extremal Reissner-Nordstr\"{o}m backgrounds. This result shows that the extremal horizon instability persists in model nonlinear theories. Our proof crucially relies on a new vector field method that allows us to obtain almost sharp decay estimates.


title: Selected Topics in Numerical Methods for Cosmology

abstract: The large amount of cosmological data already available (and in the near future) makes necessary the development of efficient numerical codes. Many software products have been implemented to perform cosmological analyses considering one or few probes. The need of multi-task software is rapidly increasing, in order to combine numerous cosmological probes along with their specificity (e.g., astrophysical descriptions and systematic errors). In this work we mention some of these libraries, bringing out some challenges they will face in the few-percent error era (on the cosmological parameters). We review some concepts of the standard cosmological model, and examine some specific topics on their implementation, bringing, for example, the discussion on how some quantities are numerically defined in different codes. We also consider implementation differences between public codes, mentioning their advantages/disadvantages.


title: Global properties of the growth index of matter inhomogeneities in the universe

abstract: We perform here a global analysis of the growth index $\gamma$ behaviour from deep in the matter era till the far future. For a given cosmological model in GR or in modified gravity, the value of $\gamma(\Omega_{m})$ is unique when the decaying mode of scalar perturbations is negligible. However, $\gamma_{\infty}$, the value of $\gamma$ in the asymptotic future, is unique even in the presence of a nonnegligible decaying mode today. Moreover $\gamma$ becomes arbitrarily large deep in the matter era. Only in the limit of a vanishing decaying mode do we get a finite $\gamma$, from the past to the future in this case. We find further a condition for $\gamma(\Omega_{m})$ to be monotonically decreasing (or increasing). This condition can be violated inside general relativity (GR) for varying $w_{DE}$ though generically $\gamma(\Omega_{m})$ will be monotonically decreasing (like $\Lambda$CDM), except in the far future and past. A bump or a dip in $G_{\rm eff}$ can also lead to a significant and rapid change in the slope $\frac{d\gamma}{d\Omega_{m}}$. On a $\Lambda$CDM background, a $\gamma$ substantially lower (higher) than $0.55$ with a negative (positive) slope reflects the opposite evolution of $G_{\rm eff}$. In DGP models, $\gamma(\Omega_{m})$ is monotonically increasing except in the far future. While DGP gravity becomes weaker than GR in the future and $w^{DGP}\to -1$, we still get $\gamma_{\infty}^{DGP}= \gamma_{\infty}^{\Lambda CDM}=\frac{2}{3}$. In contrast, despite $G^{DGP}_{\rm eff}\to G$ in the past, $\gamma$ does not tend to its value in GR because $\frac{dG^{DGP}_{\rm eff}}{d\Omega_{m}}\Big|_{-\infty}\ne 0$.


title: Invariants of long knots

abstract: By using the notion of a rigid R-matrix in a monoidal category and the Reshetikhin--Turaev functor on the category of tangles, we review the definition of the associated invariant of long knots. In the framework of the monoidal categories of relations and spans over sets, by introducing racks associated with pointed groups, we illustrate the construction and the importance of consideration of long knots. Else, by using the restricted dual of algebras and Drinfeld's quantum double construction, we show that to any Hopf algebra $H$ with invertible antipode, one can associate a universal long knot invariant $Z_H(K)$ taking its values in the convolution algebra $((D(H))^o)^*$ of the restricted dual Hopf algebra $(D(H))^o$ of the quantum double $D(H)$ of $H$. That extends the known constructions of universal invariants previously considered mostly either in the case of finite dimensional Hopf algebras or by using some topological completions.


title: Gate Leakage Suppression and Breakdown Voltage Enhancement in p-GaN HEMTs using Metal/Graphene Gates

abstract: In this work, single-layer intrinsic and fluorinated graphene were investigated as gate insertion layers in normally-OFF p-GaN gate HEMTs, which wraps around the bottom of the gate forming Ti/graphene/p-GaN at the bottom and Ti/graphene/ SiNx on the two sides. Compared to the Au/Ti/p-GaN HEMTs without graphene, the insertion of graphene can increase the ION/IOFF ratios by a factor of 50, increase the VTH by 0.30 V and reduce the off-state gate leakage by 50 times. Additionally, this novel gate structure has better thermal stability. After thermal annealing at 350 {\deg}C, gate breakdown voltage holds at 12.1 V, which is first reported for Schottky gate p-GaN HEMTs. This is considered to be a result of the 0.24 eV increase in Schottky barrier height and the better quality of the Ti/graphene/p-GaN and Ti/graphene/SiNx interfaces. This approach is very effective in improving the Ion/Iff ratio and gate BV of normally-OFF GaN HEMTs.


title: ShapeCaptioner: Generative Caption Network for 3D Shapes by Learning a Mapping from Parts Detected in Multiple Views to Sentences

abstract: 3D shape captioning is a challenging application in 3D shape understanding. Captions from recent multi-view based methods reveal that they cannot capture part-level characteristics of 3D shapes. This leads to a lack of detailed part-level description in captions, which human tend to focus on. To resolve this issue, we propose ShapeCaptioner, a generative caption network, to perform 3D shape captioning from semantic parts detected in multiple views. Our novelty lies in learning the knowledge of part detection in multiple views from 3D shape segmentations and transferring this knowledge to facilitate learning the mapping from 3D shapes to sentences. Specifically, ShapeCaptioner aggregates the parts detected in multiple colored views using our novel part class specific aggregation to represent a 3D shape, and then, employs a sequence to sequence model to generate the caption. Our outperforming results show that ShapeCaptioner can learn 3D shape features with more detailed part characteristics to facilitate better 3D shape captioning than previous work.


title: Continued fractions over non-Euclidean imaginary quadratic rings

abstract: We propose and study a generalized continued fraction algorithm that can be executed in an arbitrary imaginary quadratic field, the novelty being a non-restriction to the five Euclidean cases. Many hallmark properties of classical continued fractions are shown to be retained, including exponential convergence, best-of-the-second-kind approximation quality (up to a constant), periodicity of quadratic irrational expansions, and polynomial time complexity.


title: Correlation Tensor Magnetic Resonance Imaging

abstract: Diffusional Kurtosis Imaging (DKI) is a sensitive biomarker for microstructure in health and disease. However, DKI is not specific to any microstructural property since it may emerge from several different sources. Q-space trajectory encoding has been proposed for decoupling isotropic from anisotropic kurtosis. Still, this method assume that the system is comprised of multiple Gaussian diffusion components. Here, we develop a more general framework for resolving the underlying kurtosis sources. We introduce Correlation Tensor MRI (CTI) - an approach harnessing the versatility of double diffusion encoding (DDE) and capable of explicitly decoupling isotropic and anisotropic kurtosis components from intra-compartmental kurtosis effects arising from restricted diffusion. Additionally, CTI provides an index that is potentially sensitive to intra-compartmental kurtosis. The theoretical foundations of CTI, as well as the first proof-of-concept CTI ex vivo experiments in mouse brains at a field of 16.4T, are presented. We find that anisotropic and isotropic kurtosis can decouple microscopic anisotropy from substantial partial volume effects between tissue and free water. Our intra-compartmental kurtosis index exhibited positive values in both white and grey matter tissues. Simulations in different microenvironments show, however, that our current CTI protocol for estimating intra-compartmental kurtosis is limited by higher-order terms that were not considered in this study. CTI measurements were then extended to in vivo settings and used to map healthy rat brains at 9.4T. These in vivo CTI results were consistent with our ex vivo findings. Although future studies are still required to mitigate the higher-order effects on the intra-compartmental kurtosis index, results show that CTI's more general estimates of anisotropic and isotropic kurtosis contributions are already ripe for future studies.


title: Temporal orders and causal vector for physiological data analysis

abstract: In addition to the global parameter- and time-series-based approaches, physiological analyses should constitute a local temporal one, particularly when analyzing data within protocol segments. Hence, we introduce the R package implementing the estimation of temporal orders with a causal vector (CV). It may use linear modeling or time series distance. The algorithm was tested on cardiorespiratory data comprising tidal volume and tachogram curves, obtained from elite athletes (supine and standing, in static conditions) and a control group (different rates and depths of breathing, while supine). We checked the relation between CV and body position or breathing style. The rate of breathing had a greater impact on the CV than does the depth. The tachogram curve preceded the tidal volume relatively more when breathing was slower. Clinical relevance - The method can assess (1) relationships between two signals, having one of the two time-shifted in a given range, (2) curves of most optimal inter-signal shift, and (c) their stability across time; also for other physiological studies.


title: Oxygen-based digital etching of AlGaN/GaN structures with AlN as etch-stop layers

abstract: O2-plamsa-based digital etching of Al0.25Ga0.75N with a 0.8 nm AlN spacer on GaN was investigated. At 40 W RF bias power and 40 sccm oxygen flow, the etch depth of Al0.25Ga0.75N was 5.7 nm per cycle. The 0.8 nm AlN spacer layer acted as an etch-stop layer in 3 cycles. The surface roughness improved to 0.33 nm after 7 digital etch cycles. Compared to the dry etch only approach, this technique causes less damages. It was shown to be effective in precisely controlling the AlGaN etch depth required for recessed-AlGaN HEMTs.


title: Threshold Voltage Improvement and Leakage Reduction of AlGaN/GaN HEMTs Using Dual-Layer SiNx Stressors

abstract: In this work, AlGaN/GaN HEMTs with dual-layer SiNx stressors (composed of a low-stress layer and a high-stress layer) were investigated. The low-stress padding layer solved the surface damage problem caused during the deposition of the high-stress SiNx, and provided a good passivated interface. The HEMTs with dual-layer stressors showed a 1 V increase in the threshold voltage (Vth) with comparable on-current and RF current gain to those without stressors. Moreover, the off-current (I_off) was shown to be reduced by one to three orders of magnitude in the strained devices as a result of the lower electric field in AlGaN, which suppressed the gate injection current. The dual-layer stressor scheme supports strain engineering as an effective approach in the pursuit of the normally-off operation of AlGaN/GaN HEMTs.


title: Tensor network approach to real-time path integral

abstract: We present a tensor network representation of the path integral for the one-component real scalar field theory in 1+1 dimensional Minkowski space-time. It is numerically verified by comparing with the exact result in the non-interacting case.


title: Self-avoiding walks and polygons on hyperbolic graphs

abstract: We prove that for the $d$-regular tessellations of the hyperbolic plane by $k$-gons, there are exponentially more self-avoiding walks of length $n$ than there are self-avoiding polygons of length $n$. We then prove that this property implies that the self-avoiding walk is ballistic, even on an arbitrary vertex-transitive graph. Moreover, for every fixed $k$, we show that the connective constant for self-avoiding walks satisfies the asymptotic expansion $d-1-O(1/d)$ as $d\to \infty$; on the other hand, the connective constant for self-avoiding polygons remains bounded. Finally, we show for all but two tessellations that the number of self-avoiding walks of length $n$ is comparable to the $n$th power of their connective constant. Some of these results were previously obtained by Madras and Wu \cite{MaWuSAW} for all but finitely many regular tessellations of the hyperbolic plane.


title: A Method to Determine Broadband Complex Permittivity of Thin Film Dielectric Materials up to 60 GHz

abstract: This paper describes a method to determine the complex permittivity of a thin dielectric film from finite element analysis and microstrip line measurements. Two transmission line equivalent circuit models were used for the cases of an air-filled line and a lossless line, whose distributed elements can be calculated from full wave finite element simulations. With these calculated distributed elements and microstrip line measurements, the complex permittivity was extracted. The technique utilizes a simple way to separate the dielectric loss from the measured total loss and the complex permittivity was extracted using the measured propagation constant. A rational dielectric model was employed to fit the extracted complex permittivity, which ensures causality of the final solution. The best fitting results obtained through this procedure are considered as the final permittivity results, which have shown excellent match to the data sheet values. Moreover, simulations using the fitted permittivity exhibit good agreement with the experimental propagation constant data of microstrip lines up to 60 GHz. The proposed method was demonstrated on polyimide, and it can be applied to other thin film materials.


title: On the geometry of lattices and finiteness of Picard groups

abstract: Let $(K,\mathcal O, k)$ be a $p$-modular system with $k$ algebraically closed and $\mathcal O$ unramified, and let $\Lambda$ be an $\mathcal O$-order in a separable $K$-algebra. We call a $\Lambda$-lattice $L$ rigid if ${\rm Ext}^1_{\Lambda}(L,L)=0$, in analogy with the definition of rigid modules over a finite-dimensional algebra. By partitioning the $\Lambda$-lattices of a given dimension into "varieties of lattices", we show that there are only finitely many rigid $\Lambda$-lattices $L$ of any given dimension. As a consequence we show that if the first Hochschild cohomology of $\Lambda$ vanishes, then the Picard group and the outer automorphism group of $\Lambda$ are finite. In particular the Picard groups of blocks of finite groups defined over $\mathcal O$ are always finite.


title: Online terrain estimation for autonomous vehicles on deformable terrains

abstract: In this work, a terrain estimation framework is developed for autonomous vehicles operating on deformable terrains. Previous work in this area usually relies on steady state tire operation, linearized classical terramechanics models, or on computationally expensive algorithms that are not suitable for real-time estimation. To address these shortcomings, this work develops a reduced-order nonlinear terramechanics model as a surrogate of the Soil Contact Model (SCM) through extending a state-of-the-art Bekker model to account for additional dynamic effects. It is shown that this reduced-order surrogate model is able to accurately replicate the forces predicted by the SCM while reducing the computation cost by an order of magnitude. This surrogate model is then utilized in a unscented Kalman filter to estimate the sinkage exponent. Simulations suggest this parameter can be estimated within 4% of its true value for clay and sandy loam terrains. It is also shown that utilizing this estimated parameter can reduce the prediction errors of the future vehicle states by orders of magnitude, which could assist with achieving more robust model-predictive autonomous navigation strategies.


title: Complexity of Proximal augmented Lagrangian for nonconvex optimization with nonlinear equality constraints

abstract: We analyze worst-case complexity of a Proximal augmented Lagrangian (Proximal AL) framework for nonconvex optimization with nonlinear equality constraints. When an approximate first-order (second-order) optimal point is obtained in the subproblem, an $\epsilon$ first-order (second-order) optimal point for the original problem can be guaranteed within $\mathcal{O}(1/ \epsilon^{2 - \eta})$ outer iterations (where $\eta$ is a user-defined parameter with $\eta\in[0,2]$ for the first-order result and $\eta \in [1,2]$ for the second-order result) when the proximal term coefficient $\beta$ and penalty parameter $\rho$ satisfy $\beta = \mathcal{O}(\epsilon^\eta)$ and $\rho = \Omega (1/\epsilon^\eta)$, respectively. We also investigate the total iteration complexity and operation complexity when a Newton-conjugate-gradient algorithm is used to solve the subproblems. Finally, we discuss an adaptive scheme for determining a value of the parameter $\rho$ that satisfies the requirements of the analysis.


title: High-Resolution Near Infrared Spectroscopy of HD~100546: IV. Orbiting Companion Disappears on Schedule

abstract: HD~100546 is a Herbig Ae/Be star surrounded by a disk with a large central region that is cleared of gas and dust (i.e., an inner hole). High-resolution near-infrared spectroscopy reveals a rich emission spectrum of fundamental ro-vibrational CO emission lines whose time variable properties point to the presence of an orbiting companion within the hole. The Doppler shift and spectroastrometric signal of the CO v=1-0 P26 line, observed from 2003 to 2013, are consistent with a source of excess CO emission that orbits the star near the inner rim of the disk. The properties of the excess emission are consistent with those of a circumplanetary disk. In this paper, we report follow up observations that confirm our earlier prediction that the orbiting source of excess emission would disappear behind the near side of the inner rim of the outer disk in 2017. We find that while the hotband CO lines remained unchanged in 2017, the v=1-0 P26 line and its spectroastrometric signal returned to the profile observed in 2003. With these new observations, we further constrain the origin of the emission and discuss possible ways of confirming the presence of an orbiting planetary companion in the inner disk.


title: Modelling and Simulation of of high efficiency GaAs PIN-Solar Cell

abstract: Solar energy is the most convenient and reliable energy source among all renewable energy resources and an efficient photovoltaic device is required to convert this energy into utilizable energy. Different types of solar cells (SC) are commercially available. However, various parameters need to be optimized to get maximum efficiency from a SC. In this study we have presented a SC model in which dependence of quantum efficiency (QE) on various parameters has been investigated. The mobility of the carriers has been varied with wide range along with the carrier life time (LT). Results show that maximum efficiencies can be achieved up to 11.10% and 10.81% keeping the electron and hole mobility to be 1500 cm2V-1s-1 and 300 cm2V-1s-1 respectively with electron and hole carrier LT to be 3ns and 7ns respectively. The effect of surface recombination velocity (SRV) has also been brought under observation and the maximum efficiency is found to be 13.75% at electron and hole SRV equal to be 103ms-1. Results shows that the higher photovoltaic efficiencies can be achieved by increasing the mobility and carrier LT while decreasing the surface recombination velocities.


title: Effects of Hole-Transport Layer Homogeneity in Organic Solar Cells, A Multi-Length Scale Study

abstract: Irreproducibility is a serious issue in thin film organic photovoltaic (OPV) devices, as smallest local inhomogeneities can change the entire behaviour of identically built devices without showing obvious failure. Inhomogeneities can occur at various steps of device preparation and appear in all layers with different length scales and impact. The hole-transport interlayer (HTL) in OPV devices blocks unwanted electron diffusion to the anode and corrects energetic mismatch between oxide electrode and organic semiconductor. Most commonly used is commercial ink based on poly(3,4-ethylenedioxythiophene) poly(styrenesulfonate) (PEDOT PSS) colloidal particles. However, exactly these are suspected to cause microscopic inhomogeneities, causing known irreproducibility of device characteristics. Considering PEDOT:PSS acidity-caused electrode corrosion, it is questionable how much impact colloids have on device homogeneity. In this report, we give proof that a colloidal HTL does not necessarily cause device inhomogeneity and decreased efficiency, by comparing OPV devices with different HTLs, namely from commercial PEDOT:PSS ink and from MoO3, obtained from two liquid precursors, leading to quasi-continuous or colloidal layers. With a combination of X-ray diffraction, atomic force and Kelvin probe microscopy, photoelectron and ambient air photoemission spectroscopy, we discuss the layers properties from nano- to macroscale and demonstrate their impact upon implementation into OPV devices, via spatially-resolved characterization.


title: Influence of Environmentally Affected Hole Transport Layers on Spatial Homogeneity and Charge Transport Dynamics of Organic Solar Cells

abstract: After the efficiency of organic photovoltaic (OPV) cells achieved more than 10%, the control of stability and degradation mechanisms of solar cells became a prominent task. The improvement of device efficiency due to incorporation of a hole-transport layer (HTL) in bulk-heterojunction solar cells has been extensively reported. However, the most widely used HTL material, PEDOT:PSS is frequently suspected to be the dominating source for devices instability under environmental conditions. Thereby effects like photooxidation and electrode corrosion are often reported to shorten device lifetime. However, often in environmental device studies, the source of degradation, whether being from the HTL, the active layer or the metal cathode are rather difficult to distinguish, because the external diffusion of oxygen and water affects all components. In this study, different HTLs, namely prepared from traditional PEDOT:PSS and also two types of molybdenum trioxide (MoO3), are exposed to different environments such as oxygen, light or humidity, prior to device finalization under inert conditions. This allows investigating any effects within the HTL and from reactions at its interface to the indium-tin-oxide electrode or the active layer. The surface and bulk chemistry of the exposed HTL has been monitored and discussed in context to the observed device physics, dynamic charge transport and spatial performance homogeneity of the according OPV device. The results show that merely humidity-exposure of the HTL leads to decreased device performance for PEDOT:PSS, but also for one type of the tested MoO3. The losses are related to the amount of absorbed water in the HTL, inducing loss of active area in terms of interfacial contact. The device with PEDOT:PSS HTL after humid air exposure showed seriously decreased photocurrent by micro-delamination of swelling/shrinkage of the hygroscopic layer.


abstract: A covariant simultaneous action for branes in an arbitrary curved background spacetime is considered. The action depends on a pair of independent field variables, the brane embedding functions, through the canonical momentum of a reparametrization invariant geometric model for the brane, and an auxiliary vector field. The form of the action is analogous to a symplectic potential. Extremization of the simultaneous action produces at once the equations of motion and the Jacobi equations for the brane geometric model, and it also provides a convenient shortcut towards its second variation. In this note, we consider geometric models depending only on the intrinsic geometry of the brane worldvolume, and discuss briefly the generalization to extrinsic geometry dependent models. The approach is illustrated for Dirac-NambuGoto [DNG] branes. For a relativistic particle, a simultaneous action was introduced by Bazanski, that served as an inspiration for the present work.


title: Excitation of bound states in the continuum via second harmonic generations

abstract: A bound state in the continuum (BIC) on a periodic structure sandwiched between two homogeneous media is a guided mode with a frequency and a wavenumber such that propagating plane waves with the same frequency and wavenumber exist in the homogeneous media. Optical BICs are of significant current interest, since they have applications in lasing, sensing, filtering, switching, and many light emission processes, but they cannot be excited by incident plane waves when the structure consists of linear materials. In this paper, we study the diffraction of a plane wave by a periodic structure with a second order nonlinearity, assuming the structure has a BIC and the frequency and wavenumber of the incident wave are one half of those of the BIC. Based on a scaling analysis and a perturbation theory, we show that the incident wave may induce a very strong second harmonic wave dominated by the BIC, and also a fourth harmonic wave that cannot be ignored. The perturbation theory reveals that the amplitude of the BIC is inversely proportional to a small parameter depending on the amplitude of the incident wave and the nonlinear coefficient. In addition, a system of four nonlinearly coupled Helmholtz equations (the four-wave model) is proposed to model the nonlinear process. Numerical solutions of the four-wave model are presented for a periodic array of circular cylinders and used to validate the perturbation results.


title: Does the spacelike character of the Minkowski four-momentum show up in analog gravity?

abstract: In analog gravity the recent experiment of Drori {\it et al.} [Phys. Rev. Lett. {\bf 122}, 010404 (2019)] is impressive, as it shows how the emission of two Hawking quanta emitted in opposite directions lead to measurable consequences in the medium's rest system in a straightforward way. This result raises however the following problem: how can this experiment be explained in terms of classical electrodynamics? There must necessarily exist such an explanation (the experiment is after all classical); otherwise classical electrodynamics would be an incomplete theory. This is the main topic of the present paper. We propose that the measured effect is a demonstration of the spacelike character of the Minkowski four-momentum. Moreover, we extend the discussion by analyzing a Gedanken experiment (making use of the Kerr effect as a formal agency), to illustrate the transition from subluminal to superluminal phenomena in a straightforward way. Finally, we emphasize the close relationship that exists between the spacelike Minkowski momentum and the anomalous Doppler effect.


title: Modeling Kepler Eclipsing Binaries: Homogeneous Inference of Orbital & Stellar Properties

abstract: We report on the properties of eclipsing binaries from the Kepler mission with a newly developed photometric modeling code, which uses the light curve, spectral energy distribution of each binary, and stellar evolution models to infer stellar masses without the need for radial velocity measurements. We present solutions and posteriors to orbital and stellar parameters for 728 systems, forming the largest homogeneous catalogue of full Kepler binary parameter estimates to date. Using comparisons to published radial velocity measurements, we demonstrate that the inferred properties (e.g., masses) are reliable for well-detached main-sequence binaries, which make up the majority of our sample. The fidelity of our inferred parameters degrades for a subset of systems not well described by input isochrones, such as short-period binaries that have undergone interactions, or binaries with post-main sequence components. Additionally, we identify 35 new systems which show evidence of eclipse timing variations, perhaps from apsidal motion due to binary tides or tertiary companions. We plan to subsequently use these models to search for and constrain the presence of circumbinary planets in Kepler eclipsing binary systems.


title: "Sliced" Subwindow Search: a Sublinear-complexity Solution to the Maximum Rectangle Problem

abstract: Considering a 2D matrix of positive and negative numbers, how might one draw a rectangle within it whose contents sum higher than all other rectangles'? This fundamental problem, commonly known the maximum rectangle problem or subwindow search, spans many computational domains. Yet, the problem has not been solved without demanding computational resources at least linearly proportional to the size of the matrix. In this work, we present a new approach to the problem which achieves sublinear time and memory complexities by interpolating between a small amount of equidistant sections of the matrix. Applied to natural images, our solution outperforms the state-of-the-art by achieving an 11x increase in speed and memory efficiency at 99% comparative accuracy. In general, our solution outperforms existing solutions when matrices are sufficiently large and a marginal decrease in accuracy is acceptable, such as in many problems involving natural images. As such, it is well-suited for real-time application and in a variety of computationally hard instances of the maximum rectangle problem.


title: Projection pursuit based generalized betas accounting for higher order co-moment effects in financial market analysis

abstract: Betas are possibly the most frequently applied tool to analyze how securities relate to the market. While in very widespread use, betas only express dynamics derived from second moment statistics. Financial returns data often deviate from normal assumptions in the sense that they have significant third and fourth order moments and contain outliers. This paper targets to introduce a way to calculate generalized betas that also account for higher order moment effects, while maintaining the conceptual simplicity and interpretability of betas. Thereunto, the co-moment analysis projection index (CAPI) is introduced. When applied as a projection index in the projection pursuit (PP) framework, generalized betas are obtained as the directions optimizing the CAPI objective. A version of CAPI based on trimmed means is introduced as well, which is more stable in the presence of outliers. Simulation results underpin the statistical properties of all projections and a small, yet highly illustrative example is presented.


title: An L0-Norm Constrained Non-Negative Matrix Factorization Algorithm for the Simultaneous Disaggregation of Fixed and Shiftable Loads

abstract: Energy disaggregation refers to the decomposition of energy use time series data into its constituent loads. This paper decomposes daily use data of a household unit into fixed loads and one or more classes of shiftable loads. The latter is characterized by ON OFF duty cycles. A novel algorithm based on nonnegative matrix factorization NMF for energy disaggregation is proposed, where fixed loads are represented in terms of real-valued basis vectors, whereas shiftable loads are divided into binary signals. This binary decomposition approach directly applies L0 norm constraints on individual shiftable loads. The new approach obviates the need for more computationally intensive methods e.g. spectral decomposition or mean field annealing that have been used in earlier research for these constraints. A probabilistic framework for the proposed approach has been addressed. The proposed approach s effectiveness has been demonstrated with real consumer energy data.


title: Bilinear embedding for Schr\"odinger-type operators with complex coefficients

abstract: We prove a variant of the so-called bilinear embedding theorem for operators in divergence form with complex coefficients and with nonnegative locally integrable potentials, subject to mixed boundary conditions, and acting on arbitrary open subsets of $\mathbb{R}^{d}$.


title: Massive MIMO Channel Estimation with an Untrained Deep Neural Network

abstract: This paper proposes a deep learning-based channel estimation method for multi-cell interference-limited massive MIMO systems, in which base stations equipped with a large number of antennas serve multiple single-antenna users. The proposed estimator employs a specially designed deep neural network (DNN) to first denoise the received signal, followed by a conventional least-squares (LS) estimation. We analytically prove that our LS-type deep channel estimator can approach minimum mean square error (MMSE) estimator performance for high-dimensional signals, while avoiding MMSE's requirement for complex channel inversions and knowledge of the channel covariance matrix. This analytical result, while asymptotic, is observed in simulations to be operational for just 64 antennas and 64 subcarriers per OFDM symbol. The proposed method also does not require any training and utilizes several orders of magnitude fewer parameters than conventional DNNs. The proposed deep channel estimator is also robust to pilot contamination and can even completely eliminate it under certain conditions.


title: Interplay between boundary conditions and Wilson's mass in Dirac-like Hamiltonians

abstract: Dirac-like Hamiltonians, linear in momentum $k$, describe the low-energy physics of a large set of novel materials, including graphene, topological insulators, and Weyl fermions. We show here that the inclusion of a minimal $k^2$ Wilson's mass correction improves the models and allows for systematic derivations of appropriate boundary conditions for the envelope functions on finite systems. Considering only Wilson's masses allowed by symmetry, we show that the $k^2$ corrections are equivalent to Berry-Mondragon's discontinuous boundary conditions. This allows for simple numerical implementations of regularized Dirac models on a lattice, while properly accounting for the desired boundary condition. We apply our results on graphene nanoribbons (zigzag and armchair), and on a PbSe monolayer (topological crystalline insulator). For graphene, we find generalized Brey-Fertig boundary conditions, which correctly describe the small gap seen on \textit{ab initio} data for the metallic armchair nanoribbon. On PbSe, we show how our approach can be used to find spin-orbital-coupled boundary conditions. Overall, our discussions are set on a generic model that can be easily generalized for any Dirac-like Hamiltonian.