title: Evolutionary Variational Optimization of Generative Models

abstract: We combine two popular optimization approaches to derive learning algorithms for generative models: variational optimization and evolutionary algorithms. The combination is realized for generative models with discrete latents by using truncated posteriors as the family of variational distributions. The variational parameters of truncated posteriors are sets of latent states. By interpreting these states as genomes of individuals and by using the variational lower bound to define a fitness, we can apply evolutionary algorithms to realize the variational loop. The used variational distributions are very flexible and we show that evolutionary algorithms can effectively and efficiently optimize the variational bound. Furthermore, the variational loop is generally applicable ("black box") with no analytical derivations required. To show general applicability, we apply the approach to three generative models (we use noisy-OR Bayes Nets, Binary Sparse Coding, and Spike-and-Slab Sparse Coding). To demonstrate effectiveness and efficiency of the novel variational approach, we use the standard competitive benchmarks of image denoising and inpainting. The benchmarks allow quantitative comparisons to a wide range of methods including probabilistic approaches, deep deterministic and generative networks, and non-local image processing methods. In the category of "zero-shot" learning (when only the corrupted image is used for training), we observed the evolutionary variational algorithm to significantly improve the state-of-the-art in many benchmark settings. For one well-known inpainting benchmark, we also observed state-of-the-art performance across all categories of algorithms although we only train on the corrupted image. In general, our investigations highlight the importance of research on optimization methods for generative models to achieve performance improvements.


title: Modulation spaces associated to tensor products of amalgam spaces

abstract: We identify the modulation spaces associated to tensor products of amalgam spaces having a large class of Banach spaces as their local component. As consequences of the main results, we describe the modulation spaces associated to tensor products of various $L^p$ spaces.


title: Semi-analytical model for planetary resonances: application to planets around single and binary stars

abstract: In spite of planetary resonances being a common dynamical mechanism acting on planetary systems, no general model exists for describing their properties, particularly for commensurabilities of any order and arbitrary values of the eccentricities and inclinations. The present work presents a semi-analytical model that describes the resonance strength, width, location and stability of fixed points, as well as periods of small-amplitude librations. The model is valid to any two gravitationally interacting massive bodies, and thus applicable to planets around single or binary stars. Using a theoretical framework in Poincar\'e and Jacobi reference system we develop a semi-analytical method that employs a numerical evaluation of the averaged resonant disturbing function. Validations of the model are presented comparing its predictions with dynamical maps for real and fictitious systems. The model is shown to describe very well many dynamical features of planetary resonances. Notwithstanding the good agreements found in all cases, a small deviation is noted in the location of the resonance centers for circumbinary systems. As a consequence of its application to the HD31527 system we have found that the updated best-fit solution leads to a high-eccentricity stable libration between the middle and outer planets inside the 16/3 mean-motion resonance. This is the first planetary system whose long-term dynamics appears dominated by such a high-order commensurability. In the case of circumbinary planets, the overlap of N/1 mean-motion resonances coincides very well with the size of the global chaotic region close to the binary, as well as its dependence with the mutual inclination.


title: Comparative Analysis of the Solar Wind: Modeling Charge State Distributions in the Heliosphere

abstract: Non-equilibrium ionization (NEI) is a key process often times neglected when modeling astrophysical plasmas with thermodynamical timescales much shorter than the timescales for ionization and recombination. In this paper, we perform NEI modeling on a magnetohydrodynamic (MHD) simulation of the solar wind during the Whole Sun Month (Carrington Rotation 1913 from 1996 August 22 to September 18), and compare the resulting charge state distributions with in situ measurements made with the Solar Wind Ion Composition Spectrometer (SWICS) on Ulysses. We trace the wind trajectory back to 20$R_\odot$ using the velocity measured by Ulysses at about 4 AU, and obtain the plasma flow trajectory within 20$R_\odot$ using the MHD simulation data performed by Predictive Science Inc. We assume that the wind started from the solar surface and was initially in ionization equilibrium, and analyze the time-dependent ionization state of solar wind along each of the trajectories. In our analysis we: (1) obtain charge state densities and ratios for slow and fast winds based on the NEI model, and compare them with in situ observations; and (2) measure the "freeze-in" distance for several ions observed by SWICS to determine a possible correlation between when the ionization states become fixed and the electron density and outward velocity of the plasma at the freeze-in height. This study provides a stringent test on comparing charge state distributions from the outer corona predicted from simulations with in situ measurements made in the far heliosphere. This showcases the challenges in matching plasma conditions near the corona to those observed in interplanetary space.


title: Zeros of Gaussian Weyl-Heisenberg functions and hyperuniformity of charge

abstract: We study Gaussian random functions on the complex plane whose stochastics are invariant under the Weyl-Heisenberg group (twisted stationarity). The theory is modeled on translation invariant Gaussian entire functions, but allows for non-analytic examples, in which case winding numbers can be either positive or negative. We calculate the first intensity of zero sets of such functions, both when considered as points on the plane, or as charges according to their phase winding. In the latter case, charges are shown to be in a certain average equilibrium independently of the particular covariance structure (universal screening). We investigate the corresponding fluctuations, and show that in many cases they are suppressed at large scales (hyperuniformity). This means that universal screening is empirically observable at large scales. We also derive an asymptotic expression for the charge variance. As a main application, we obtain statistics for the zero sets of the short-time Fourier transform of complex white noise with general windows, and also prove the following uncertainty principle: the expected number of zeros per unit area is minimized, among all window functions, exactly by generalized Gaussians. Further applications include poly-entire functions such as covariant derivatives of Gaussian entire functions.


title: CFHTLenS: Galaxy bias as function of scale, stellar mass, and colour. Conflicts with predictions by semi-analytic models

abstract: Galaxy models predict a tight relation between the clustering of galaxies and dark matter on cosmological scales, but predictions differ notably in the details. We used this opportunity and tested two semi-analytic models by the Munich and Durham groups with data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). For the test we measured the scale-dependent galaxy bias factor $b(k)$ and correlation factor $r(k)$ from linear to non-linear scales of $k\approx10\,h\,\rm Mpc^{-1}$ at two redshifts $\bar{z}=0.35,0.51$ for galaxies with stellar mass between $5\times10^9$ and $3\times10^{11}\,h_{\rm 70}^{-2}\,{\rm M}_\odot$. Our improved gravitational lensing technique accounts for the intrinsic alignment of sources and the magnification of lens galaxies for better constraints for the galaxy-matter correlation $r(k)$. Galaxy bias in CFHTLenS increases with $k$ and stellar mass, it is colour-dependent, revealing the individual footprints of galaxy types. Despite a reasonable model agreement for the relative change with both scale and galaxy properties, there is a clear conflict for $b(k)$ with no model preference: the model galaxies are too weakly clustered. This may flag a model problem at $z\gtrsim0.3$ for all stellar masses. As in the models, however, there is a high correlation $r(k)$ between matter and galaxy density on all scales, and galaxy bias is typically consistent with a deterministic bias on linear scales. Only our blue and low-mass galaxies of about $7\times10^9\,h_{\rm 70}^{-2}\,{\rm M}_\odot$ at $\bar{z}=0.51$ show, contrary to the models, a weak tendency towards a stochastic bias on linear scales where $r_{\rm ls}=0.75\pm0.14\,{\rm(stat.)}\pm0.06\,{\rm(sys.)}$. This result is of interest for cosmological probes, such as $E_{\rm G}$, that rely on a deterministic galaxy bias.


title: Finite size effects around pseudo-transition in one-dimensional models with nearest neighbor interaction

abstract: Recently gigantic peaks in thermodynamic response functions have been observed at finite temperature for one-dimensional models with short-range coupling, closely resembling a second-order phase transition. Thus, we will analyze the finite temperature pseudo-transition property observed in some one-dimensional models and its relationship with finite size effect. In particular, we consider two chain models to study the finite size effects; these are the Ising-Heisenberg tetrahedral chain and an Ising-Heisenberg-type ladder model. Although the anomalous peaks of these one-dimensional models have already been studied in the thermodynamic limit, here we will discuss the finite size effects of the chain and why the peaks do not diverge in the thermodynamic limit. So, we discuss the dependence of the finite size effects, for moderately and sufficiently large systems, in which the specific heat and magnetic susceptibility exhibit peculiar rounded towering peaks for a given temperature. This behavior is quite similar to a continuous phase transition, but there is no singularity. For moderately large systems, the peaks narrow and increase in height as the number of unit cells is increased, and the location of peak shifts slightly. Hence, one can naively induce that the sharp peak should lead to a divergence in the thermodynamic limit. However, for a rather large system, the height of a peak goes asymptotically to a finite value. Our result confirms the dependence of the peak height with the number of unit cells at the pseudo-critical temperature. We also provide an alternative empirical function that satisfactorily fits specific heat and magnetic susceptibility at pseudo-critical temperature. Certainly, our result is crucial to understand the finite size correction behavior in quantum spin models, which in general are only numerically tractable within the framework of the finite size analysis.


title: Mueller matrix maps of dichroic filters reveal polarization aberrations

abstract: Dichroic filters are used by instrument designers to split a field of view into different optical paths for simultaneous measurement of different spectral bands. Quantifying the polarization aberrations of a dichroic is relevant for predicting the incident polarization states downstream, which could affect the performance of diffraction limited systems. One important application is the fore-optics of exoplanet imaging coronagraphs. In this work, the polarization properties of the Edmund #69-205 650 nm roll-off dichroic are measured using a rotating retarder Mueller matrix imaging polarimeter. The polarization properties of this commercial dichroic are compared at normal and 45$^{\circ}$ angle of incidence. The normal incidence measurements verify the instrument calibration since no polarization aberrations were observed. Transmission measurements at 680 nm and 45$^{\circ}$ yield a 2.9 rad magnitude of retardance and 0.95 diattenuation. Effectively, at 630 nm the dichroic is a $\lambda$/4 waveplate with a horizontal fast-axis.


title: Flexible deep transfer learning by separate feature embeddings and manifold alignment

abstract: Object recognition is a key enabler across industry and defense. As technology changes, algorithms must keep pace with new requirements and data. New modalities and higher resolution sensors should allow for increased algorithm robustness. Unfortunately, algorithms trained on existing labeled datasets do not directly generalize to new data because the data distributions do not match. Transfer learning (TL) or domain adaptation (DA) methods have established the groundwork for transferring knowledge from existing labeled source data to new unlabeled target datasets. However, current DA approaches assume similar source and target feature spaces and suffer in the case of massive domain shifts or changes in the feature space. Existing methods assume the data are either the same modality, or can be aligned to a common feature space. Therefore, most methods are not designed to support a fundamental domain change such as visual to auditory data. We propose a novel deep learning framework that overcomes this limitation by learning separate feature extractions for each domain while minimizing the distance between the domains in a latent lower-dimensional space. The alignment is achieved by considering the data manifold along with an adversarial training procedure. We demonstrate the effectiveness of the approach versus traditional methods with several ablation experiments on synthetic, measured, and satellite image datasets. We also provide practical guidelines for training the network while overcoming vanishing gradients which inhibit learning in some adversarial training settings.


title: Exact Christoffel-Darboux Expansions: A New, Multidimensional, Algebraic, Eigenenergy Bounding Method

abstract: Although the Christoffel-Darboux representation (CDR) plays an important role within the theory of orthogonal polynomials, and many important bosonic and fermionic multidimensional Schrodinger equation systems can be transformed into a moment equation representation (MER), the union of the two into an effective, algebraic, eigenenergy bounding method has been overlooked. This particular fusion of the two representations, suitable for bounding bosonic or fermionic systems, defines the Orthonormal Polynomial Projection Quantization - Bounding Method (OPPQ-BM), as developed here. We use it to analyze several one dimensional and two dimensional systems, including the quadratic Zeeman effect for strong-superstrong magnetic fields. For this problem, we match or surpass the excellent, but intricate, results of Kravchenko et al (1996 Phys. Rev. A 54287) for a broad range of magnetic fields, without the need for any truncations or approximations.


title: Boosted Top Quarks in the Peak Region with N$^3$LL Resummation

abstract: We present results for the 2-jettiness differential distribution for boosted top quark pairs produced in $e^+e^-$ collisions in the peak region accounting for QCD large-logarithm resummation at next-to-next-to-next-to-leading logarithmic (N$^3$LL) order and fixed-order corrections to matrix elements at next-to-next-to-leading order (NNLO) calculated in the framework of soft-collinear effective theory and boosted heavy quark effective theory. Electroweak and finite-width effects are included at leading order. We study the perturbative convergence of the cross section in the pole and MSR mass schemes, with and without soft gap subtractions. We find that there is a partial cancellation between the pole mass and soft function renormalons. When renormalon subtractions concerning the top mass and the soft function are implemented, the perturbative uncertainties are, however, systematically smaller and an improvement in the stability of the peak position is observed. We find that the top MSR mass may be determined with perturbative uncertainties well below $100$\,MeV from the peak position of the 2-jettiness distribution. This result has important applications for Monte Carlo top quark mass calibrations.


title: Confronting Abusive Language Online: A Survey from the Ethical and Human Rights Perspective

abstract: The pervasiveness of abusive content on the internet can lead to severe psychological and physical harm. Significant effort in Natural Language Processing (NLP) research has been devoted to addressing this problem through abusive content detection and related sub-areas, such as the detection of hate speech, toxicity, cyberbullying, etc. Although current technologies achieve high classification performance in research studies, it has been observed that the real-life application of this technology can cause unintended harms, such as the silencing of under-represented groups. We review a large body of NLP research on automatic abuse detection with a new focus on ethical challenges, organized around eight established ethical principles: privacy, accountability, safety and security, transparency and explainability, fairness and non-discrimination, human control of technology, professional responsibility, and promotion of human values. In many cases, these principles relate not only to situational ethical codes, which may be context-dependent, but are in fact connected to universal human rights, such as the right to privacy, freedom from discrimination, and freedom of expression. We highlight the need to examine the broad social impacts of this technology, and to bring ethical and human rights considerations to every stage of the application life-cycle, from task formulation and dataset design, to model training and evaluation, to application deployment. Guided by these principles, we identify several opportunities for rights-respecting, socio-technical solutions to detect and confront online abuse, including `nudging', `quarantining', value sensitive design, counter-narratives, style transfer, and AI-driven public education applications.


title: Pattern Recognition Scheme for Large-Scale Cloud Detection over Landmarks

abstract: Landmark recognition and matching is a critical step in many Image Navigation and Registration (INR) models for geostationary satellite services, as well as to maintain the geometric quality assessment (GQA) in the instrument data processing chain of Earth observation satellites. Matching the landmark accurately is of paramount relevance, and the process can be strongly impacted by the cloud contamination of a given landmark. This paper introduces a complete pattern recognition methodology able to detect the presence of clouds over landmarks using Meteosat Second Generation (MSG) data. The methodology is based on the ensemble combination of dedicated support vector machines (SVMs) dependent on the particular landmark and illumination conditions. This divide-and-conquer strategy is motivated by the data complexity and follows a physically-based strategy that considers variability both in seasonality and illumination conditions along the day to split observations. In addition, it allows training the classification scheme with millions of samples at an affordable computational costs. The image archive was composed of 200 landmark test sites with near 7 million multispectral images that correspond to MSG acquisitions during 2010. Results are analyzed in terms of cloud detection accuracy and computational cost. We provide illustrative source code and a portion of the huge training data to the community.


title: Nonlinear Cook distance for Anomalous Change Detection

abstract: In this work we propose a method to find anomalous changes in remote sensing images based on the chronochrome approach. A regressor between images is used to discover the most {\em influential points} in the observed data. Typically, the pixels with largest residuals are decided to be anomalous changes. In order to find the anomalous pixels we consider the Cook distance and propose its nonlinear extension using random Fourier features as an efficient nonlinear measure of impact. Good empirical performance is shown over different multispectral images both visually and quantitatively evaluated with ROC curves.


title: Randomized RX for target detection

abstract: This work tackles the target detection problem through the well-known global RX method. The RX method models the clutter as a multivariate Gaussian distribution, and has been extended to nonlinear distributions using kernel methods. While the kernel RX can cope with complex clutters, it requires a considerable amount of computational resources as the number of clutter pixels gets larger. Here we propose random Fourier features to approximate the Gaussian kernel in kernel RX and consequently our development keep the accuracy of the nonlinearity while reducing the computational cost which is now controlled by an hyperparameter. Results over both synthetic and real-world image target detection problems show space and time efficiency of the proposed method while providing high detection performance.


title: Influence Maximization Under Generic Threshold-based Non-submodular Model

abstract: As a widely observable social effect, influence diffusion refers to a process where innovations, trends, awareness, etc. spread across the network via the social impact among individuals. Motivated by such social effect, the concept of influence maximization is coined, where the goal is to select a bounded number of the most influential nodes (seed nodes) from a social network so that they can jointly trigger the maximal influence diffusion. A rich body of research in this area is performed under statistical diffusion models with provable submodularity, which essentially simplifies the problem as the optimal result can be approximated by the simple greedy search. When the diffusion models are non-submodular, however, the research community mostly focuses on how to bound/approximate them by tractable submodular functions so as to estimate the optimal result. In other words, there is still a lack of efficient methods that can directly resolve non-submodular influence maximization problems. In this regard, we fill the gap by proposing seed selection strategies using network graphical properties in a generalized threshold-based model, called influence barricade model, which is non-submodular. Specifically, under this model, we first establish theories to reveal graphical conditions that ensure the network generated by node removals has the same optimal seed set as that in the original network. We then exploit these theoretical conditions to develop efficient algorithms by strategically removing less-important nodes and selecting seeds only in the remaining network. To the best of our knowledge, this is the first graph-based approach that directly tackles non-submodular influence maximization.


title: Mixture Model Framework for Traumatic Brain Injury Prognosis Using Heterogeneous Clinical and Outcome Data

abstract: Prognoses of Traumatic Brain Injury (TBI) outcomes are neither easily nor accurately determined from clinical indicators. This is due in part to the heterogeneity of damage inflicted to the brain, ultimately resulting in diverse and complex outcomes. Using a data-driven approach on many distinct data elements may be necessary to describe this large set of outcomes and thereby robustly depict the nuanced differences among TBI patients' recovery. In this work, we develop a method for modeling large heterogeneous data types relevant to TBI. Our approach is geared toward the probabilistic representation of mixed continuous and discrete variables with missing values. The model is trained on a dataset encompassing a variety of data types, including demographics, blood-based biomarkers, and imaging findings. In addition, it includes a set of clinical outcome assessments at 3, 6, and 12 months post-injury. The model is used to stratify patients into distinct groups in an unsupervised learning setting. We use the model to infer outcomes using input data, and show that the collection of input data reduces uncertainty of outcomes over a baseline approach. In addition, we quantify the performance of a likelihood scoring technique that can be used to self-evaluate the extrapolation risk of prognosis on unseen patients.


title: Influencer Videos: Unboxing the Mystique

abstract: Influencer marketing has become a very popular tool to reach customers. Despite the rapid growth in influencer videos, there has been little research on the effectiveness of their constituent features in explaining video engagement. We study YouTube influencers and analyze their unstructured video data across text, audio and images using an "interpretable deep learning" framework that accomplishes both goals of prediction and interpretation. Our prediction-based approach analyzes unstructured data and finds that "what is said" in words (text) is more influential than "how it is said" in imagery (images) or acoustics (audio). Our novel interpretation-based approach is implemented after completion of model prediction by analyzing the same source of unstructured data to measure importance attributed to the video features. We eliminate several spurious relationships in two steps, identifying a subset of relationships which are confirmed using theory. We uncover novel findings that establish distinct associations for measures of shallow and deep engagement based on the dual-system framework of human thinking. Our approach is validated using simulated data, and we discuss the learnings from our findings for influencers and brands.


title: Predictions of astrometric jitter for Sun-like stars. I. The model and its application to the Sun as seen from the ecliptic

abstract: The advent of Gaia, capable of measuring stellar wobbles caused by orbiting planets, raised an interest to the astrometric detection of exoplanets. Another source of such wobbles (often also called jitter) is stellar magnetic activity. A quantitative assessment of the stellar astrometric jitter is important for a more reliable astrometric detection and characterisation of exoplanets. We calculate the displacement of the solar photocentre due to the magnetic activity for an almost 16-year period (February 2, 1999 - August 1, 2014). We also investigate how the displacement depends on the spectral passband chosen for observations, including the wavelength range to be covered by the upcoming Small-JASMINE mission of JAXA. This is done by extending the SATIRE-S model for solar irradiance variability to calculating the displacement of the solar photocentre caused by the magnetic features on the surface of the Sun. We found that the peak to peak amplitude of the solar photocentre displacement would reach 0.5 mas if the Sun were located 10 pc away from the observer and observed in the Gaia G filter. This is by far too small to be detected by the Gaia mission. However, the Sun is a relatively inactive star so that one can expect significantly larger signals for younger, and, consequently, more active stars. The model developed in this study can be combined with the simulations of emergence and surface transport of magnetic flux which have recently became available to model the astrometric jitter over the broad range of magnetic activities.


title: Charge sensitivity of a cavity-embedded Cooper pair transistor limited by single-photon shot noise

abstract: Using an operator scattering approach, we analyze the quantum dynamics of an ultrasensitive electrometer -- a Cooper pair transistor embedded in a quarter-wave microwave cavity (cCPT). While the cCPT is inherently a tunable, strongly nonlinear system affording a diverse range of functionalities, we restrict our present analysis to a necessary first investigation of its linear charge sensing capabilities, limiting to low pump powers corresponding to an average cavity photon number $\lesssim 1$. Assuming realizable cCPT parameters (B. L. Brock et al., Phys. Rev. Applied, 044009), and not including noise from the subsequent amplifier chain, we predict the fundamental, photon shot noise-limited charge sensitivity to be $0.12\, \mu{\mathrm{e}}/\sqrt{\mathrm{Hz}}$ when the pumped cavity has an average of one photon and the cCPT is operated close to charge degeneracy. This is to be compared with a first reported charge sensitivity value $14\, \mu{\mathrm{e}}/\sqrt{\mathrm{Hz}}$ in the single-photon regime (B. L. Brock et al., arXiv:2102.05362).


title: Hierarchical Recurrent Attention Networks for Structured Online Maps

abstract: In this paper, we tackle the problem of online road network extraction from sparse 3D point clouds. Our method is inspired by how an annotator builds a lane graph, by first identifying how many lanes there are and then drawing each one in turn. We develop a hierarchical recurrent network that attends to initial regions of a lane boundary and traces them out completely by outputting a structured polyline. We also propose a novel differentiable loss function that measures the deviation of the edges of the ground truth polylines and their predictions. This is more suitable than distances on vertices, as there exists many ways to draw equivalent polylines. We demonstrate the effectiveness of our method on a 90 km stretch of highway, and show that we can recover the right topology 92\% of the time.


title: Gamma attenuation parameters of some borate glasses

abstract: In this study we have reported the borate-based glasses in the form of $4(75 - x)B-2O-3 - xBi_2O_3 - 10Na_2O - 10CaO - 5Al_2O_3$ ($0 \le x \le 25$ mol%) for gamma-ray and neutron shielding applications, and for this purpose mass attenuation coefficients of these glasses have been investigated using Geant4 code and XCOM program in the energy range between 1 keV and 10 MeV. A good agreement between the simulation results and computed values was observed. The obtained values were then used to calculate the other related shielding parameters such as half value layer (HVL), mean free path (MFP) and effective atomic numbers ($Z_{\rm eff}$). Exposure buildup factors of the glasses were also evaluated by using the geometrical progression (G-P) method. Additionally, neutron effective removal cross-sections of each selected glass have been computed. The results indicate that the Bi2O3 content in the glass sample affects the radiation shielding properties. It can be concluded that the glass coded as F offers superior gamma-ray shielding properties, while the glass coded as D offers better neutron shielding features among the other samples.


title: Fluctuations of the log-gamma polymer free energy with general parameters and slopes

abstract: We prove that the free energy of the log-gamma polymer between lattice points $(1,1)$ and $(M,N)$ converges to the GUE Tracy-Widom distribution in the $M^{1/3}$ scaling, provided that $N/M$ remains bounded away from zero and infinity. We prove this result for the model with inverse gamma weights of any shape parameter $\theta>0$ and furthermore establish a moderate deviation estimate for the upper tail of the free energy in this case. Finally, we consider a non i.i.d. setting where the weights on finitely many rows and columns have different parameters, and we show that when these parameters are critically scaled the limiting free energy fluctuations are governed by a generalization of the Baik-Ben Arous-P\'ech\'e distribution from spiked random matrices with two sets of parameters.


title: On Holder Continuity and Equivalent Formulation of Intrinsic Harnack Estimates for an Anisotropic Parabolic Degenerate Prototype Equation

abstract: We give a proof of H\"older continuity for bounded local weak solutions to the equation $u_t= \sum_{i=1}^N (|u_{x_i}|^{p_i-2} u_{x_i})_{x_i}$, in $\Omega \times [0,T]$, with $\Omega \subset \subset \mathbb{R}^N$, under the condition $ 2<p_i<\bar{p}(1+2/N)$ for each $i=1,..,N$, being $\bar{p}$ the harmonic mean of the $p_i$s, via recently discovered intrinsic Harnack estimates. Moreover we establish equivalent forms of these Harnack estimates within the proper intrinsic geometry.


title: The Keck Planet Imager and Characterizer: Phase I fiber injection unit early performance and commissioning

abstract: The Keck Planet Imager and Characterizer (KPIC) is an upgrade to the Keck II adaptive optics system and instrument suite with the goal of improving direct imaging and high-resolution spectroscopic characterization capabilities for giant exoplanets. KPIC Phase I includes a fiber injection unit (FIU) downstream of a new pyramid wavefront sensor, coupling planet light to a single mode fiber fed into NIRSPEC, Keck's high-resolution infrared spectrograph. This enables high-dispersion spectroscopy (HDS) of directly imaged exoplanets at smaller separation and higher contrast, improving our spectral characterization capabilities for these objects. Here, we report performance results from the KPIC Phase I FIU commissioning, including analysis of throughput, stability, and sensitivity of the instrument.


title: Letter to the Editor and Comments on: Intravital dynamic and correlative imaging reveals diffusion-dominated canalicular and flow-augmented ductular bile flux

abstract: Bile, the central metabolic product of the liver, is secreted by hepatocytes into bile canaliculi (BC), tubular subcellular structures of 0.5-2 $\mu$m diameter which are formed by the apical membranes of juxtaposed hepatocytes. BC interconnect to build a highly ramified 3D network that collects and transports bile towards larger interlobular bile ducts (IBD). The transport mechanism of bile is of fundamental interest and the current text-book model of "osmotic canalicular bile flow", i.e. enforced by osmotic water influx, has recently been quantitatively verified based on flux measurements and BC contractility (Meyer et al., 2017) but challenged in the article entitled "Intravital dynamic and correlative imaging reveals diffusion-dominated canalicular and flow-augmented ductular bile flux" (Vartak et al., 2020). We feel compelled to share a series of arguments that question the key conclusions of this study (Vartak et al., 2020), particularly because it risks to be misinterpreted and misleading for the field in view of potential clinical applications. We summarized the arguments in the following 8 points.


title: SU(2/1) superchiral self-duality: a new quantum, algebraic and geometric paradigm to describe the electroweak interactions

abstract: We propose an extension of the Yang-Mills paradigm from Lie algebras to internal chiral superalgebras. We replace the Lie algebra-valued connection one-form $A$, by a superalgebra-valued polyform $\widetilde{A}$ mixing exterior-forms of all degrees and satisfying the chiral self-duality condition $\widetilde{A} = {}^*\widetilde{A} \,\chi$, where $\chi$ denotes the superalgebra grading operator. This superconnection contains Yang-Mills vectors valued in the even Lie subalgebra, together with scalars and self-dual tensors valued in the odd module, all coupling only to the charge parity CP-positive Fermions. The Fermion quantum loops then induce the usual Yang-Mills-scalar Lagrangian, the self-dual Avdeev-Chizhov propagator of the tensors, plus a new vector-scalar-tensor vertex and several quartic terms which match the geometric definition of the supercurvature. Applied to the $SU(2/1)$ Lie-Kac simple superalgebra, which naturally classifies all the elementary particles, the resulting quantum field theory is anomaly-free and the interactions are governed by the super-Killing metric and by the structure constants of the superalgebra.


title: Quantum Request-Answer Game with Buffer Model for Online Algorithms

abstract: We consider online algorithms as a request-answer game. An adversary that generates input requests, and an online algorithm answers. We consider a generalized version of the game that has a buffer of limited size. The adversary loads data to the buffer, and the algorithm has random access to elements of the buffer. We consider quantum and classical (deterministic or randomized) algorithms for the model. In the paper, we provide a specific problem (The Most Frequent Keyword Problem) and a quantum algorithm that works better than any classical (deterministic or randomized) algorithm in terms of competitive ratio. At the same time, for the problem, classical online algorithms in the standard model are equivalent to the classical algorithms in the request-answer game with buffer model.


title: Third components with elliptical orbits in the eclipsing binaries: BD and, SV Cam, V0836 Cyg & XZ CMi

abstract: This is our regular research in the field of cyclic O-C changes and third components. In this article the results of four new stars investigation is described. These stars are: BD And, SV Cam, V0836 Cyg and XZ CMi. All of them have cyclic O-C curve with linear (BD And, SV Cam and XZ CMi) or parabolic trend (V0836 Cyg). We computed the mass transfer rate, minima possible mass of the third component and their errors for each of the researched stars.


title: Summability of formal solutions for a family of generalized moment integro-differential equations

abstract: Generalized summability results are obtained regarding formal solutions of certain families of linear moment integro-differential equations with time variable coefficients. The main result leans on the knowledge of the behavior of the moment derivatives of the elements involved in the problem. A refinement of the main result is also provided giving rise to more accurate results which remain valid in wide families of problems of high interest in practice, such as fractional integro-differential equations.


title: Do We Need Improved Code Quality Metrics?

abstract: The software development community has been using code quality metrics for the last five decades. Despite their wide adoption, code quality metrics have attracted a fair share of criticism. In this paper, first, we carry out a qualitative exploration by surveying software developers to gauge their opinions about current practices and potential gaps with the present set of metrics. We identify deficiencies including lack of soundness, i.e., the ability of a metric to capture a notion accurately as promised by the metric, lack of support for assessing software architecture quality, and insufficient support for assessing software testing and infrastructure. In the second part of the paper, we focus on one specific code quality metric-LCOM as a case study to explore opportunities towards improved metrics. We evaluate existing LCOM algorithms qualitatively and quantitatively to observe how closely they represent the concept of cohesion. In this pursuit, we first create eight diverse cases that any LCOM algorithm must cover and obtain their cohesion levels by a set of experienced developers and consider them as a ground truth. We show that the present set of LCOM algorithms do poorly w.r.t. these cases. To bridge the identified gap, we propose a new approach to compute LCOM and evaluate the new approach with the ground truth. We also show, using a quantitative analysis using more than 90 thousand types belonging to 261 high-quality Java repositories, the present set of methods paint a very inaccurate and misleading picture of class cohesion. We conclude that the current code quality metrics in use suffer from various deficiencies, presenting ample opportunities for the research community to address the gaps.


title: Drug-Target Interaction Prediction via an Ensemble of Weighted Nearest Neighbors with Interaction Recovery

abstract: Predicting drug-target interactions (DTI) via reliable computational methods is an effective and efficient way to mitigate the enormous costs and time of the drug discovery process. Structure-based drug similarities and sequence-based target protein similarities are the commonly used information for DTI prediction. Among numerous computational methods, neighborhood-based chemogenomic approaches that leverage drug and target similarities to perform predictions directly are simple but promising ones. However, existing similarity-based methods need to be re-trained to predict interactions for any new drugs or targets and cannot directly perform predictions for both new drugs, new targets, and new drug-target pairs. Furthermore, a large amount of missing (undetected) interactions in current DTI datasets hinders most DTI prediction methods. To address these issues, we propose a new method denoted as Weighted k-Nearest Neighbor with Interaction Recovery (WkNNIR). Not only can WkNNIR estimate interactions of any new drugs and/or new targets without any need of re-training, but it can also recover missing interactions (false negatives). In addition, WkNNIR exploits local imbalance to promote the influence of more reliable similarities on the interaction recovery and prediction processes. We also propose a series of ensemble methods that employ diverse sampling strategies and could be coupled with WkNNIR as well as any other DTI prediction method to improve performance. Experimental results over five benchmark datasets demonstrate the effectiveness of our approaches in predicting drug-target interactions. Lastly, we confirm the practical prediction ability of proposed methods to discover reliable interactions that were not reported in the original benchmark datasets.


title: Electric Forces in the Charged SYK Model

abstract: The connection between gravitational force and operator growth, reported in earlier papers, is generalized to include electromagnetic forces. It is shown how in the U(1)-invariant SYK system electric forces emerge through the same mechanisms -- the momentum-size correspondence, and operator growth -- that give rise to gravitational force. The unit of electric charge implied by the U(1)-SYK theory is consistent with a Kaluza-Klein radius inverse to the characteristic energy scale of the SYK theory.


title: An Introduction to Barenblatt Solutions for Anisotropic $p$-Laplace Equations

abstract: We introduce Fundamental solutions of Barenblatt type for the equation $u_t=\sum_{i=1}^N \bigg( |u_{x_i}|^{p_i-2}u_{x_i} \bigg)_{x_i}$, $p_i >2 \quad \forall i=1,..,N$, on $\Sigma_T=\mathbb{R}^N \times[0,T]$, and we prove their importance for the regularity properties of the solutions.


title: Explicit strong boundedness for higher rank symplectic groups

abstract: This paper gives an explicit argument to show strong boundedness for ${\rm Sp}_{2n}(R)$ for $R$ a ring of S-algebraic integers or a semi-local ring. This gives a quantitative version of a related abstract result in a previous paper of the author. The results presented further generalize older results regarding strong boundedness by Kedra, Libman and Martin and Morris from ${\rm SL}_n$ to ${\rm Sp}_{2n}$. Further, the presented results solve the question of the asymptotic of strong boundedness for ${\rm Sp}_{2n}(R)$ for $R$ semi-local case with an argument that immediately generalizes to all other split Chevalley groups.


title: Hybrid RIS-Empowered Reflection and Decode-and-Forward Relaying for Coverage Extension

abstract: In this letter, we introduce two hybrid transmission schemes combining a passive reconfigurable intelligent surface (RIS) with decode-and-forward relaying in a synergistic manner. The proposed schemes offer a flexible as well as cost- and power-efficient solution for coverage extension in future generation wireless networks. We present closed-form expressions for the end-to-end signal-to-noise ratio of both schemes and a sequential optimization algorithm for the power allocation and the RIS phase configurations. Our computer simulations and theoretical analysis demonstrate that the RIS and relaying technologies enhance the achievable rate and error performance remarkably when working complementary to each other, rather than being considered as competing technologies.


title: Self-assembled monolayers of oligophenylenes stiffer than steel and silicon, possibly even stiffer than Si$_3$N$_4$

abstract: To quantify charge transport through molecular junctions fabricated using the conducting probe atomic force microscopy (CP-AFM) platform, information on the number of molecules $N$ per junction is absolutely necessary. $N$ can be currently obtained only via contact mechanics, and the Young's modulus $E$ of the self-assembled monolayer (SAM) utilized in the key quantity for this approach. The experimental determination of $E$ for SAMs of CP-AFM junctions fabricated using oligophenylene dithiols (OPDn, $1 \leq n \leq 4$) and gold electrodes turned out to be too challenging. Recent measurements (Z. Xie et al, J. Am. Chem. Soc. 139 (2017) 5696) merely succeeded to provide a low bound estimate ($E \approx 58\,$GPa). It is this state of affairs that motivated the present theoretical investigation. Our microscopic calculations yield values $E \approx 240 \pm 6\,$GPa for the OPDn SAMs of the aforementioned experimental study, which are larger than those of steel ($ E \approx 180 - 200\,$GPa) and silicon ($E \approx 130 - 185\,$GPa). The fact that the presently computed $E$ is much larger than the aforementioned experimental lower bound explain why experimentally measuring $E$ of OPDn SAM's is so challenging. Having $E \approx 337 \pm 8\,$GPa, OPDn SAMs with herringbone arrangement adsorbed on fcc (111)Au are even stiffer than Si$_3$N$_4$ ($E \approx 160 - 290\,$GPa).


title: Real-Time Vehicular Wireless System-Level Simulation

abstract: Future automation and control units for advanced driver assistance systems (ADAS) will exchange sensor and kinematic data with nearby vehicles using wireless communication links to improve traffic safety. In this paper we present an accurate real-time system-level simulation for multi-vehicle communication scenarios to support the development and test of connected ADAS systems. The physical and data-link layer are abstracted and provide the frame error rate (FER) to a network simulator. The FER is strongly affected by the non-stationary doubly dispersive fading process of the vehicular radio communication channel. We use a geometry-based stochastic channel model (GSCM) to enable a simplified but still accurate representation of the non-stationary vehicular fading process. The propagation path parameters of the GSCM are used to efficiently compute the time-variant condensed radio channel parameters per stationarity region of each communication link during run-time. Five condensed radio channel parameters mainly determine the FER forming a parameter vector: path loss, root mean square delay spread, Doppler bandwidth, $K$-factor, and line-of-sight Doppler shift. We measure the FER for a pre-defined set of discrete grid points of the parameter vector using a channel emulator and a given transmitter-receiver modem pair. The FER data is stored in a table and looked up during run-time of the real-time system-level simulation. We validate our methodology using empirical measurement data from a street crossing scenarios demonstrating a close match in terms of FER between simulation and measurement.


title: Ultraholomorphic sectorial extensions of Beurling type

abstract: We prove sectorial extension theorems for ultraholomorphic function classes of Beurling type defined by weight functions with a controlled loss of regularity. The proofs are based on a reduction lemma, due to the second author, which allows to extract the Beurling from the Roumieu case, which was treated recently by Jim\'{e}nez-Garrido, Sanz, and the third author. In order to have control on the opening of the sectors, where the extensions exist, we use the (mixed) growth index and the order of quasianalyticity of weight functions. As a consequence we obtain corresponding extension results for classes defined by weight sequences. Additionally, we give information on the existence of continuous linear extension operators.


title: Dynamical heterogeneities in non-entangled polystyrene and poly(ethylene oxide) star melts

abstract: Star polymers can exhibit a heterogeneous dynamical behavior due to their internal structure. In this work we employ atomistic molecular dynamics simulations to study translational motion in non-entangled polystyrene and poly(ethylene oxide) star-shaped melts. We focus on the local heterogeneous dynamics originating from the multi-arm star-like architecture and quantify the intramolecular dynamical gradient. By examining the translational motion at length scales of the order of the Kuhn length, we aim to find common features for both studied chemistries and to provide a critical and direct comparison with theoretical models of polymer dynamics. We discuss the observed tendencies with respect to the continuous Rouse model adjusted for the star-like architectures. Two versions of the Rouse model are examined: one assuming uniform friction on every Rouse bead and another one considering larger branch point friction. Apart from chain connectivity between neighboring beads, both versions disregard interactions between the chains. Despite the tolerable description of the simulation data, neither model appears to reflect the mobility gradient accurately. The detailed quantitative atomistic models employed here bridge the gap between the theoretical and general, coarse-granined models of star-like polymers which lack the indispensable chemical details.


title: Deep Unsupervised Image Hashing by Maximizing Bit Entropy

abstract: Unsupervised hashing is important for indexing huge image or video collections without having expensive annotations available. Hashing aims to learn short binary codes for compact storage and efficient semantic retrieval. We propose an unsupervised deep hashing layer called Bi-half Net that maximizes entropy of the binary codes. Entropy is maximal when both possible values of the bit are uniformly (half-half) distributed. To maximize bit entropy, we do not add a term to the loss function as this is difficult to optimize and tune. Instead, we design a new parameter-free network layer to explicitly force continuous image features to approximate the optimal half-half bit distribution. This layer is shown to minimize a penalized term of the Wasserstein distance between the learned continuous image features and the optimal half-half bit distribution. Experimental results on the image datasets Flickr25k, Nus-wide, Cifar-10, Mscoco, Mnist and the video datasets Ucf-101 and Hmdb-51 show that our approach leads to compact codes and compares favorably to the current state-of-the-art.


title: Goal Reasoning by Selecting Subgoals with Deep Q-Learning

abstract: In this work we propose a goal reasoning method which learns to select subgoals with Deep Q-Learning in order to decrease the load of a planner when faced with scenarios with tight time restrictions, such as online execution systems. We have designed a CNN-based goal selection module and trained it on a standard video game environment, testing it on different games (planning domains) and levels (planning problems) to measure its generalization abilities. When comparing its performance with a satisfying planner, the results obtained show both approaches are able to find plans of good quality, but our method greatly decreases planning time. We conclude our approach can be successfully applied to different types of domains (games), and shows good generalization properties when evaluated on new levels (problems) of the same game (domain).


title: An Approach to Deploy Interactive Robotic Simulators on the Web for HRI Experiments: Results in Social Robot Navigation

abstract: Evaluation of social robot navigation inherently requires human input due to its qualitative nature. Motivated by the need to scale human evaluation, we propose a general method for deploying interactive, rich-client robotic simulations on the web. Prior approaches implement specific web-compatible simulators or provide tools to build a simulator for a specific study. Instead, our approach builds on standard Linux tools to share a graphical desktop with remote users. We leverage these tools to deploy simulators on the web that would typically be constrained to desktop computing environments. As an example implementation of our approach, we introduce the SEAN Experimental Platform (SEAN-EP). With SEAN-EP, remote users can virtually interact with a mobile robot in the Social Environment for Autonomous Navigation, without installing any software on their computer or needing specialized hardware. We validated that SEAN-EP could quickly scale the collection of human feedback and its usability through an online survey. In addition, we compared human feedback from participants that interacted with a robot using SEAN-EP with feedback obtained through a more traditional video survey. Our results suggest that human perceptions of robots may differ based on whether they interact with the robots in simulation or observe them in videos. Also, they suggest that people perceive the surveys with interactive simulations as less mentally demanding than video surveys.


title: Spying on the prior of the number of data clusters and the partition distribution in Bayesian cluster analysis

abstract: Cluster analysis aims at partitioning data into groups or clusters. In applications, it is common to deal with problems where the number of clusters is unknown. Bayesian mixture models employed in such applications usually specify a flexible prior that takes into account the uncertainty with respect to the number of clusters. However, a major empirical challenge involving the use of these models is in the characterisation of the induced prior on the partitions. This work introduces an approach to compute descriptive statistics of the prior on the partitions for three selected Bayesian mixture models developed in the areas of Bayesian finite mixtures and Bayesian nonparametrics. The proposed methodology involves computationally efficient enumeration of the prior on the number of clusters in-sample (termed as ``data clusters'') and determining the first two prior moments of symmetric additive statistics characterising the partitions. The accompanying reference implementation is made available in the R package 'fipp'. Finally, we illustrate the proposed methodology through comparisons and also discuss the implications for prior elicitation in applications.


title: One-particle and two-particle visibilities in bipartite entangled Gaussian states

abstract: Complementarity between one- and two-particle visibility in discrete systems can be extended to bipartite quantum-entangled Gaussian states. The meaning of the two-particle visibility originally defined by Jaeger, Horne, Shimony, and Vaidman with the use of an indirect method that first corrects the two-particle probability distribution by adding and subtracting other distributions with varying degree of entanglement, however, deserves further analysis. Furthermore, the origin of complementarity between one-particle visibility and two-particle visibility is somewhat elusive and it is not entirely clear what is the best way to associate particular two-particle quantum observables with the two-particle visibility. Here, we develop a direct method for quantifying the two-particle visibility based on measurement of a pair of two-particle observables that are compatible with the measured pair of single-particle observables. For each of the two-particle observables the corresponding visibility is computed, after which the absolute difference of the latter pair of visibilities is considered as a redefinition of the two-particle visibility. Our approach reveals a mathematical symmetry as it treats the two pairs of one-particle or two-particle observables on equal footing by formally identifying all four observable distributions as rotated marginal distributions of the original two-particle probability distribution. The complementarity relation between one-particle visibility and two-particle visibility obtained with the direct method is exact in the limit of infinite Gaussian precision where the entangled Gaussian state approaches an ideal EPR state. The presented results demonstrate the theoretical utility of rotated marginal distributions for elucidating the nature of two-particle visibility and provide tools for the development of quantum applications employing continuous variables.


title: Arithmetic subsequences in a random ordering of an additive set

abstract: For a finite set $A$ of size $n$, an ordering is an injection from $\{1,2,\ldots,n\}$ to $A$. We present results concerning the asymptotic properties of the length $L_n$ of the longest arithmetic subsequence in a random ordering of an additive set $A$. In the torsion-free case where $A = [1,n]^d\subseteq {\bf Z}^d$, we prove that $L_n\sim 2d\log n/\log\log n$. We show that the case $A = {\bf Z}/n{\bf Z}$ behaves asymptotically like the torsion-free case with $d=1$, and then use this fact to compute the expected length of the longest arithmetic subsequence in a random ordering of an arbitrary finite abelian group. We also prove that the number of orderings of ${\bf Z}/n{\bf Z}$ without any arithmetic subsequence of length $3$ is $2^{n-1}$ when $n\geq 2$ is a power of $2$, and zero otherwise. We conclude with a concrete application to elementary $p$-groups and a discussion of possible noncommutative generalisations.


title: NMR Investigation on Honeycomb Iridate Ag$_3$LiIr$_2$O$_6$

abstract: We investigate the structural and magnetic properties of a Kitaev spin liquid candidate material Ag$_3$LiIr$_2$O$_6$ based on $^7$Li nuclear magnetic resonance line shape, Knight shift and spin-lattice relaxation rate $1/T_1$. The first sample A shows signatures of magnetically ordered spins, and exhibits one sharp $^7$Li peak with FWHM increasing significantly below 14~K. $1/T_1^{stretch}$ of this sample displays a broad local maximum at 40~K, followed by a very sharp peak at $T_N = 9\pm1$~K due to critical slowing down of Ir spin fluctuations, a typical signature of magnetic long range order. In order to shed light on the position-by-position variation of $1/T_1$ throughout the sample, we use Inverse Laplace Transform $T_1$ analysis based on Tikhonov regularization to deduce the density distribution function $P(1/T_1)$. We demonstrate that $\sim 60\%$ of Ir spins are statically ordered at the NMR measurement timescale but the rest of the sample volume remains paramagnetic even at 4.2~K, presumably because of structural disorder induced primarily by stacking faults. In order to further investigate the influence of structural disorder, we compare these NMR results with those of a second sample B, which has been shown by transmission electron microscope to have domains with unwanted Ag inclusion at Li and Ir sites within the Ir honeycomb planes. The sample B displays an additional NMR peak with relative intensity of $\sim 17\%$. The small Knight shift and $1/T_1$ of these defect-induced $^7$Li sites and the enhancement of bulk susceptibility at low temperatures suggest that these defects generate domains of only weakly magnetic Ir spins accompanied by free spins, leading to a lack of clear signatures of long-range order. The apparent lack of long-range order could be easily misinterpreted as evidence for the realization of a spin liquid ground state in highly disordered Kitaev lattice.


title: Confinement and Graded Partition Functions for $\mathcal{N}=4$ SYM

abstract: Gauge theories with confining phases at low temperatures tend to deconfine at high temperatures. In some cases, for example in supersymmetric theories, confinement can persist for all temperatures provided the partition function includes a grading by $(-1)^F$. When it is possible to define partition functions which smoothly interpolate between no grading and $(-1)^F$ grading, it is natural to ask if there are other choices of grading that have the same effect as $(-1)^F$ on confinement. We explore how this works for $\mathcal{N}=4$ SYM on $S^1\times S^3$ in the large $N$ limit at both small and large coupling. We find evidence for a continuous range of grading parameters that preserve confinement for all temperatures at large coupling, while at small coupling only a discrete set of gradings preserves confinement.


title: Inverse classical scattering using fractional derivative

abstract: The fractional calculus framework will be used to invert the potential energy function from the classical scattering angle, which will be related to Riemann-Liouville fractional integral. Numerical solution of this fractional order problem will be applied to the inverse Rutherford scattering and to the inverse scattering of Xe--Rn atoms, in which the potential is given by Lennard-Jones function. Proofs of existence will be presented for more clarity and completness of the present work. In the two cases considered, the potential energy function can be retrieved with a desired precision. The present method gives a clear understanding of the inverse fractional problem framework.


title: Vacuum Polarization Energy of the Kinks in the Sinh-Deformed Models

abstract: We compute the one-loop quantum corrections to the kink energies of the sinh-deformed $\phi^{4}$ and $\varphi^{6}$ models in one space and one time dimensions. These models are constructed from the well-known polynomial $\phi^{4}$ and $\varphi^{6}$ models by a deformation procedure. We also compute the vacuum polarization energy to the non-polynomial function $U(\phi)=\frac{1}{4}(1-\sinh^{2}\phi)^{2}$. This potential approaches the $\phi^{4}$ model in the limit of small values of the scalar function. These energies are extracted from scattering data for fluctuations about the kink solutions. We show that for certain topological sectors with non-equivalent vacua the kink solutions of the sinh-deformed models are destabilized.


title: Thermoelectric properties of InAs nanowires from full-band atomistic simulations

abstract: In this work we theoretically explore the effect of dimensionality on the thermoelectric power factor of InAs nanowires by coupling atomistic tight-binding calculations to the Linearized Boltzmann transport formalism. We consider nanowires with diameters from 40nm (bulk-like) down to 3nm (1D), which allows for the proper exploration of the power factor within a unified large-scale atomistic description across a large diameter range. We find that as the diameter of the nanowires is reduced below d < 10 nm, the Seebeck coefficient increases substantially, a consequence of strong subband quantization. Under phonon-limited scattering conditions, a considerable improvement of ~6x in the power factor is observed around d = 10 nm. The introduction of surface roughness scattering in the calculation reduces this power factor improvement to ~2x. As the diameter is decreased down to d = 3 nm, the power factor is diminished. Our results show that, although low effective mass materials such as InAs can reach low-dimensional behavior at larger diameters and demonstrate significant thermoelectric power factor improvements, surface roughness is also stronger at larger diameters, which takes most of the anticipated power factor advantages away. However, the power factor improvement that can be observed around d = 10 nm, could prove to be beneficial as both the Lorenz number and the phonon thermal conductivity are reduced at that diameter. Thus, this work, by using large-scale full-band simulations that span the corresponding length scales, clarifies properly the reasons behind power factor improvements (or degradations) in low-dimensional materials. The elaborate computational method presented can serve as a platform to develop similar schemes for 2D and 3D material electronic structures.


title: COVID-19 and other viruses: holding back its expansion by massive testing

abstract: The experience of Singapur and South Korea makes it clear that under certain circumstances massive testing is an effective way for containing the advance of the COVID-19. In this paper, we propose a modified SEIR model which takes into account tracing and massive testing, proving theoretically that more tracing and testing implies a reduction of the total number of infected people in the long run. We apply this model to the spread of the first wave of the disease in Spain, obtaining numerical results. After that, we introduce a heuristic approach in order to minimize the COVID-19 spreading by planning effective test distributions among the populations of a region over a period of time. As an application, the impact of distributing tests among the counties of New York according to this method is computed in terms of the number of saved infected individuals.


title: A machine learning solver for high-dimensional integrals: Solving Kolmogorov PDEs by stochastic weighted minimization and stochastic gradient descent through a high-order weak approximation scheme of SDEs with Malliavin weights

abstract: The paper introduces a very simple and fast computation method for high-dimensional integrals to solve high-dimensional Kolmogorov partial differential equations (PDEs). The new machine learning-based method is obtained by solving a stochastic weighted minimization with stochastic gradient descent which is inspired by a high-order weak approximation scheme for stochastic differential equations (SDEs) with Malliavin weights. Then solutions to high-dimensional Kolmogorov PDEs or expectations of functionals of solutions to high-dimensional SDEs are accurately approximated without suffering from the curse of dimensionality. Numerical examples for PDEs and SDEs up to 100 dimensions are shown by using second and third-order discretization schemes in order to demonstrate the effectiveness of our method.


title: Beating Random Assignment for Approximating Quantum 2-Local Hamiltonian Problems

abstract: The quantum k-Local Hamiltonian problem is a natural generalization of classical constraint satisfaction problems (k-CSP) and is complete for QMA, a quantum analog of NP. Although the complexity of k-Local Hamiltonian problems has been well studied, only a handful of approximation results are known. For Max 2-Local Hamiltonian where each term is a rank 3 projector, a natural quantum generalization of classical Max 2-SAT, the best known approximation algorithm was the trivial random assignment, yielding a 0.75-approximation. We present the first approximation algorithm beating this bound, a classical polynomial-time 0.764-approximation. For strictly quadratic instances, which are maximally entangled instances, we provide a 0.801 approximation algorithm, and numerically demonstrate that our algorithm is likely a 0.821-approximation. We conjecture these are the hardest instances to approximate. We also give improved approximations for quantum generalizations of other related classical 2-CSPs. Finally, we exploit quantum connections to a generalization of the Grothendieck problem to obtain a classical constant-factor approximation for the physically relevant special case of strictly quadratic traceless 2-Local Hamiltonians on bipartite interaction graphs, where a inverse logarithmic approximation was the best previously known (for general interaction graphs). Our work employs recently developed techniques for analyzing classical approximations of CSPs and is intended to be accessible to both quantum information scientists and classical computer scientists.


title: An overview on deep learning-based approximation methods for partial differential equations

abstract: It is one of the most challenging problems in applied mathematics to approximatively solve high-dimensional partial differential equations (PDEs). Recently, several deep learning-based approximation algorithms for attacking this problem have been proposed and tested numerically on a number of examples of high-dimensional PDEs. This has given rise to a lively field of research in which deep learning-based methods and related Monte Carlo methods are applied to the approximation of high-dimensional PDEs. In this article we offer an introduction to this field of research by revisiting selected mathematical results related to deep learning approximation methods for PDEs and reviewing the main ideas of their proofs. We also provide a short overview of the recent literature in this area of research.


title: A study of measure-theoretic area formulas

abstract: We present a complete study of measure-theoretic area formulas in metric spaces, providing different measurability conditions.


title: ActionBert: Leveraging User Actions for Semantic Understanding of User Interfaces

abstract: As mobile devices are becoming ubiquitous, regularly interacting with a variety of user interfaces (UIs) is a common aspect of daily life for many people. To improve the accessibility of these devices and to enable their usage in a variety of settings, building models that can assist users and accomplish tasks through the UI is vitally important. However, there are several challenges to achieve this. First, UI components of similar appearance can have different functionalities, making understanding their function more important than just analyzing their appearance. Second, domain-specific features like Document Object Model (DOM) in web pages and View Hierarchy (VH) in mobile applications provide important signals about the semantics of UI elements, but these features are not in a natural language format. Third, owing to a large diversity in UIs and absence of standard DOM or VH representations, building a UI understanding model with high coverage requires large amounts of training data. Inspired by the success of pre-training based approaches in NLP for tackling a variety of problems in a data-efficient way, we introduce a new pre-trained UI representation model called ActionBert. Our methodology is designed to leverage visual, linguistic and domain-specific features in user interaction traces to pre-train generic feature representations of UIs and their components. Our key intuition is that user actions, e.g., a sequence of clicks on different UI components, reveals important information about their functionality. We evaluate the proposed model on a wide variety of downstream tasks, ranging from icon classification to UI component retrieval based on its natural language description. Experiments show that the proposed ActionBert model outperforms multi-modal baselines across all downstream tasks by up to 15.5%.


title: Is the brain macroscopically linear? A system identification of resting state dynamics

abstract: A central challenge in the computational modeling of neural dynamics is the trade-off between accuracy and simplicity. At the level of individual neurons, nonlinear dynamics are both experimentally established and essential for neuronal functioning. An implicit assumption has thus formed that an accurate computational model of whole-brain dynamics must also be highly nonlinear, whereas linear models may provide a first-order approximation. Here, we provide a rigorous and data-driven investigation of this hypothesis at the level of whole-brain blood-oxygen-level-dependent (BOLD) and macroscopic field potential dynamics by leveraging the theory of system identification. Using functional MRI (fMRI) and intracranial EEG (iEEG), we model the resting state activity of 700 subjects in the Human Connectome Project (HCP) and 122 subjects from the Restoring Active Memory (RAM) project using state-of-the-art linear and nonlinear model families. We assess relative model fit using predictive power, computational complexity, and the extent of residual dynamics unexplained by the model. Contrary to our expectations, linear auto-regressive models achieve the best measures across all three metrics, eliminating the trade-off between accuracy and simplicity. To understand and explain this linearity, we highlight four properties of macroscopic neurodynamics which can counteract or mask microscopic nonlinear dynamics: averaging over space, averaging over time, observation noise, and limited data samples. Whereas the latter two are technological limitations and can improve in the future, the former two are inherent to aggregated macroscopic brain activity. Our results, together with the unparalleled interpretability of linear models, can greatly facilitate our understanding of macroscopic neural dynamics and the principled design of model-based interventions for the treatment of neuropsychiatric disorders.


title: Seeing past words: Testing the cross-modal capabilities of pretrained V&L models on counting tasks

abstract: We investigate the reasoning ability of pretrained vision and language (V&L) models in two tasks that require multimodal integration: (1) discriminating a correct image-sentence pair from an incorrect one, and (2) counting entities in an image. We evaluate three pretrained V&L models on these tasks: ViLBERT, ViLBERT 12-in-1 and LXMERT, in zero-shot and finetuned settings. Our results show that models solve task (1) very well, as expected, since all models are pretrained on task (1). However, none of the pretrained V&L models is able to adequately solve task (2), our counting probe, and they cannot generalise to out-of-distribution quantities. We propose a number of explanations for these findings: LXMERT (and to some extent ViLBERT 12-in-1) show some evidence of catastrophic forgetting on task (1). Concerning our results on the counting probe, we find evidence that all models are impacted by dataset bias, and also fail to individuate entities in the visual input. While a selling point of pretrained V&L models is their ability to solve complex tasks, our findings suggest that understanding their reasoning and grounding capabilities requires more targeted investigations on specific phenomena.


title: On the structure of certain $\Gamma$-difference modules

abstract: This is a largely expository paper, providing a self-contained account on the results of [Sch-Si1, Sch-Si2], in the cases denoted there 2Q and 2M. These papers of Sch\"afke and Singer supplied new proofs to the main theorems of [Bez-Bou, Ad-Be], on the rationality of power series satisfying a pair of independent q-difference, or Mahler, equations. We emphasize the language of $\Gamma$-difference modules, instead of difference equations or systems. Although in the two cases mentioned above this is only a semantic change, we also treat a new case, which may be labeled 1M1Q. Here the group $\Gamma$ is generalized dihedral rather than abelian, and the language of equations is inadequate. In the last section we explain how to generalize the main theorems in case 2Q to finite characteristic.


title: Injection locking of a levitated optomechanical oscillator for precision force sensing

abstract: We report on the injection locking of an optically levitated nanomechanical oscillator (a silica nanosphere) to resonant intensity modulations of an external optical signal. We explore the characteristic features of injection locking in this system, e.g. the phase pull-in effect and the injection-induced reduction of the oscillation linewidth. Our measurements are in good agreement with theoretical predictions and deepen the analogy of injection locking in levitated optomechanical systems to that in optical systems (lasers). By measuring the force noise of our feedback cooled free-running oscillator, we attain a force sensitivity of $\sim23~\rm{zN}/\sqrt{\rm{Hz}}$. This can readily allow, in fairly short integration times, for tests of violations of Newtonian gravity and searching for new small-scale forces. As a proof of concept, we show that the injection locking can be exploited to measure the forces optically induced on levitated nanoparticles, with potential applications in explorations of optical binding and entanglement between optically coupled nanomechanical oscillators.


title: What can we learn from the conformal noninvariance of the Klein-Gordon equation?

abstract: It is well known that the Klein-Gordon equation in curved spacetime is conformally noninvariant, both with and without a mass term. We show that such a noninvariance provides nontrivial physical insights at different levels, first within the fully relativistic regime, then in the nonrelativistic regime leading to the Schr\"odinger equation, and then within the de Broglie-Bohm causal interpretation of quantum mechanics. The conformal noninvariance of the Klein-Gordon equation coupled to a vector potential is confronted with the conformal invariance of Maxwell's equations in the presence of a charged current. The conformal invariance of the non-minimally coupled Klein-Gordon equation to gravity is then examined in light of the conformal invariance of Maxwell's equations. Finally, the consequence of the noninvariance of the equation on the Aharonov-Bohm effect in curved spacetime is discussed.


title: Unbiased Subdata Selection for Fair Classification: A Unified Framework and Scalable Algorithms

abstract: As an important problem in modern data analytics, classification has witnessed varieties of applications from different domains. Different from conventional classification approaches, fair classification concerns the issues of unintentional biases against the sensitive features (e.g., gender, race). Due to high nonconvexity of fairness measures, existing methods are often unable to model exact fairness, which can cause inferior fair classification outcomes. This paper fills the gap by developing a novel unified framework to jointly optimize accuracy and fairness. The proposed framework is versatile and can incorporate different fairness measures studied in literature precisely as well as can be applicable to many classifiers including deep classification models. Specifically, in this paper, we first prove Fisher consistency of the proposed framework. We then show that many classification models within this framework can be recast as mixed-integer convex programs, which can be solved effectively by off-the-shelf solvers when the instance sizes are moderate and can be used as benchmarks to compare the efficiency of approximation algorithms. We prove that in the proposed framework, when the classification outcomes are known, the resulting problem, termed "unbiased subdata selection," is strongly polynomial-solvable and can be used to enhance the classification fairness by selecting more representative data points. This motivates us to develop an iterative refining strategy (IRS) to solve the large-scale instances, where we improve the classification accuracy and conduct the unbiased subdata selection in an alternating fashion. We study the convergence property of IRS and derive its approximation bound. More broadly, this framework can be leveraged to improve classification models with unbalanced data by taking F1 score into consideration.


title: Persistence properties of a Camassa-Holm type equation with $(n+1)-$order non-linearities

abstract: Lower order conservation laws and symmetries of a family of hyperbolic equations having the Camassa-Holm equation as a particular member are obtained. We show that the equation has two conservation laws with zeroth order characteristics and that its symmetries are generated by translations in the independent variables and a certain scaling, as well as some invariant solutions are studied. Next, we consider persistence and asymptotic properties for the solutions of the equation considered. In particular, we analyse the behaviour of the solutions of the equation for large values of the spatial variable. We show that if the initial data has a certain asymptotic exponential decaying, then such property persists for any time as long as the solution exists. Moreover, depending on the behaviour of the initial data for large values of the spatial variable and if for some further time the solution shares the same behaviour, then it necessarily vanishes identically. Finally, we prove unique continuation results for the solutions of the equation.


title: A new unified model of dark matter and dark energy in 5-dimensional $f(R,\phi)$ gravity

abstract: We propose a new unified model that describes~dark energy and dark matter in the context of $f(R,\phi )$ gravity using a massive scalar field in five dimensions. The scalar field is considered in the bulk that surrounds the 3-brane in branworld model. We show that for a specific choice of the $% f(R,\phi )$ function, we can describe the Einstein gravitation in 4-dimensional space-time. We obtain a relationship between the speed of the universe's expansion and the speed of the bulk's expansion. We also propose that the dark matter is represented by the scalar field mass and that the dark energy is a kinetic energy of this field. Finally, we show that, according to conditions, one can obtain the percentages of density\ of dark matter and the density of ordinary matter.


title: Topological K-theory for discrete groups and Index theory

abstract: We give a complete solution, for discrete countable groups, to the problem of defining and computing a geometric pairing between the left hand side of the Baum-Connes assembly map, given in terms of geometric cycles associated to proper actions on manifolds, and cyclic periodic cohomology of the group algebra. Indeed, for any such group $\Gamma$ (without any further assumptions on it) we construct an explicit morphism from the Left-Hand side of the Baum-Connes assembly map to the periodic cyclic homology of the group algebra. This morphism, called here the Chern-Baum-Connes assembly map, allows to give a proper and explicit formulation for a Chern-Connes pairing with the periodic cyclic cohomology of the group algebra. Several theorems are needed to formulate the Chern-Baum-Connes assembly map. In particular we establish a delocalised Riemann-Roch theorem, the wrong way functoriality for periodic delocalised cohomology for $\Gamma$-proper actions, the construction of a Chern morphism between the Left-Hand side of Baum-Connes and a delocalised cohomology group associated to $\Gamma$ which is an isomorphism once tensoring with $\mathbb{C}$, and the construction of an explicit cohomological assembly map between the delocalised cohomology group associated to $\Gamma$ and the homology group $H_*(\Gamma,F\Gamma)$. We then give an index theoretical formula for the above mentioned pairing (for any $\Gamma$) in terms of pairings of invariant forms, associated to geometric cycles and given in terms of delocalized Chern and Todd classes, and currents naturally associated to group cocycles using Burghelea's computation. As part of our results we prove that left-Hand side group used in this paper is isomorphic to the usual analytic model for the left-hand side of the assembly map.


title: A Structure-Aware Method for Direct Pose Estimation

abstract: Estimating camera pose from a single image is a fundamental problem in computer vision. Existing methods for solving this task fall into two distinct categories, which we refer to as direct and indirect. Direct methods, such as PoseNet, regress pose from the image as a fixed function, for example using a feed-forward convolutional network. Such methods are desirable because they are deterministic and run in constant time. Indirect methods for pose regression are often non-deterministic, with various external dependencies such as image retrieval and hypothesis sampling. We propose a direct method that takes inspiration from structure-based approaches to incorporate explicit 3D constraints into the network. Our approach maintains the desirable qualities of other direct methods while achieving much lower error in general.


title: Properties of RG interfaces for 2D boundary flows

abstract: We consider RG interfaces for boundary RG flows in two-dimensional QFTs. Such interfaces are particular boundary condition changing operators linking the UV and IR conformal boundary conditions. We refer to them as RG operators. In this paper we study their general properties putting forward a number of conjectures. We conjecture that an RG operator is always a conformal primary such that the OPE of this operator with its conjugate must contain the perturbing UV operator when taken in one order and the leading irrelevant operator (when it exists) along which the flow enters the IR fixed point, when taken in the other order. We support our conjectures by perturbative calculations for flows between nearby fixed points, by a non-perturbative variational method inspired by the variational method proposed by J.~Cardy for massive RG flows, and by numerical results obtained using boundary TCSA. The variational method has a merit of its own as it can be used as a first approximation in charting the global structure of the space of boundary RG flows. We also discuss the role of the RG operators in the transport of states and local operators. Some of our considerations can be generalised to two-dimensional bulk flows, clarifying some conceptual issues related to the RG interface put forward by D.~Gaiotto for bulk $\phi_{1,3}$ flows.


title: Designing an Adaptive Bandwidth Management for Higher Education Institutions

abstract: Purpose: This study proposes an adaptive bandwidth management system which can be explicitly used by educational institutions. The primary goal of the system is to increase the bandwidth of the users who access more on educational websites. Through this proposed bandwidth management, the users of the campus networks is encouraged to utilize the internet for educational purposes. Method: The weblog from a university's pfSense proxy server was utilized and undergo Web Usage Mining (WUM) to determine the number of educational and non-educational websites accessed by the users. Certain formulas were used in the computation of the bandwidth which was dynamically assigned to the users. A prototyping technique was applied in developing adaptive bandwidth management system. The prototype was simulated and evaluated by experts in compliance with ISO/IEC 14598-6 and ISO/IEC 9126-1 standards. Results: This study found that the prototype is capable of adjusting the bandwidth of the network users dynamically. The users who browsed more on educational websites or contents were assigned with higher bandwidth compared to those who are not. Further, the evaluated prototype met the software standards of ISO. Conclusion: The proposed adaptive bandwidth management can contribute to the continuous development in the area of computer networking, especially in designing and managing campus networks. It also helps the network administrators or IT managers in allocating bandwidth with minimal effort.


title: The Circlet Inequalities: A New, Circulant-Based Facet-Defining Inequality for the TSP

abstract: Facet-defining inequalities of the symmetric Traveling Salesman Problem (TSP) polytope play a prominent role in both polyhedral TSP research and state-of-the-art TSP solvers. In this paper, we introduce a new class of facet-defining inequalities, the \emph{circlet inequalities}. These inequalities were first conjectured in Gutekunst and Williamson \cite{Gut19b} when studying Circulant TSP, and they provide a bridge between polyhedral TSP research and number-theoretic investigations of Hamiltonian cycles stemming from a conjecture due to Marco Buratti in 2017. The circlet inequalities exhibit circulant symmetry by placing the same weight on all edges of a given length; our main proof exploits this symmetry to prove the validity of the circlet inequalities. We then show that the circlet inequalities are facet-defining and compute their strength following Goemans \cite{Goe95b}; they achieve the same worst-case strength as the similarly circulant crown inequalities of Naddef and Rinaldi \cite{Nad92}, but are generally stronger.


title: Effect of inter-system coupling on heat transport in a microscopic collision model

abstract: In this paper we consider a bipartite system composed of two subsystems each coupled to its own thermal environment. Based on a collision model, we mainly study whether the approximation (i.e., the inter-system interaction is ignored when modeling the system-environment coupling) is valid or not. We also address the problem of heat transport unitedly for both conventional energy-preserving system-environment interactions and non-energy preserving system-environment interactions. For the former interaction, as the inter-system interaction strength increases, at first this approximation gets worse as expected, but then counterintuitively gets better even for a stronger inter-system coupling. For the latter interaction with asymmetry, this approximation gets progressively worse. In this case we realize a perfect thermal rectification, and we can not find apparent rectification effect for the former interaction. Finally and more importantly, our results show that whether this approximation is valid or not is closely related to the quantum correlations between the subsystems, i.e., the weaker the quantum correlations, the more justified the approximation and vice versa.


title: Deep decarbonization and the Supergrid, prospects for electricity transmission between Europe and China

abstract: Long distance transmission within continents has been shown to be one of the most important variation management strategies in renewable energy systems, where allowing for transmission expansion will reduce system cost by around 20%. In this paper, we test whether the system cost further decreases when transmission is extended to intercontinental connections. We analyze a Eurasian interconnection between China, Mid-Asia and Europe, using a capacity expansion model with hourly time resolution. The model is constrained by an increasingly tighter global cap on CO2 emissions in order to investigate the effect of different levels of reliance on variable sources. Our results show that a supergrid option decreases total system cost by a maximum of 5%, compared to continental grid integration. This maximum effect is achieved when (i) the generation is constrained to be made up almost entirely by renewables, (ii) the land available for VRE farms is relatively limited and the demand is relatively high and (iii) the cost for solar PV and storage is high. The importance of these two factors is explained by that a super grid allows for harnessing of remote wind-, solar- and hydro resources as well as management of variations, both of which are consequential only in cases where dispatchable resources are limited or very costly. As for the importance of the cost for storage, it represents a competing variation management option, and when it has low cost, it substitutes part of the role of the supergrid, which is to manage variations through long-distance trade. The cost decrease from a Eurasian supergrid was found to be between 0% and 5%, compared to the cost in the case of continental-scale grids. We conclude that the benefits of a supergrid from a techno-economic perspective are in most cases negligible, or modest at best.17 pages, f


title: Multi-Head Self-Attention with Role-Guided Masks

abstract: The state of the art in learning meaningful semantic representations of words is the Transformer model and its attention mechanisms. Simply put, the attention mechanisms learn to attend to specific parts of the input dispensing recurrence and convolutions. While some of the learned attention heads have been found to play linguistically interpretable roles, they can be redundant or prone to errors. We propose a method to guide the attention heads towards roles identified in prior work as important. We do this by defining role-specific masks to constrain the heads to attend to specific parts of the input, such that different heads are designed to play different roles. Experiments on text classification and machine translation using 7 different datasets show that our method outperforms competitive attention-based, CNN, and RNN baselines.


title: Unbiased Gradient Estimation for Distributionally Robust Learning

abstract: Seeking to improve model generalization, we consider a new approach based on distributionally robust learning (DRL) that applies stochastic gradient descent to the outer minimization problem. Our algorithm efficiently estimates the gradient of the inner maximization problem through multi-level Monte Carlo randomization. Leveraging theoretical results that shed light on why standard gradient estimators fail, we establish the optimal parameterization of the gradient estimators of our approach that balances a fundamental tradeoff between computation time and statistical variance. Numerical experiments demonstrate that our DRL approach yields significant benefits over previous work.


title: Understanding and Increasing Efficiency of Frank-Wolfe Adversarial Training

abstract: Deep neural networks are easily fooled by small perturbations known as adversarial attacks. Adversarial Training (AT) is a technique that approximately solves a robust optimization problem to minimize the worst-case loss and is widely regarded as the most effective defense. Due to the high computation time for generating strong adversarial examples in the AT process, single-step approaches have been proposed to reduce training time. However, these methods suffer from catastrophic overfitting where adversarial accuracy drops during training, and although improvements have been proposed, they increase training time and robustness is far from that of multi-step AT. We develop a theoretical framework for adversarial training with FW optimization (FW-AT) that reveals a geometric connection between the loss landscape and the $\ell_2$ distortion of $\ell_\infty$ FW attacks. We analytically show that high distortion of FW attacks is equivalent to small gradient variation along the attack path. It is then experimentally demonstrated on various deep neural network architectures that $\ell_\infty$ attacks against robust models achieve near maximal distortion, while standard networks have lower distortion. It is experimentally shown that catastrophic overfitting is strongly correlated with low distortion of FW attacks. This mathematical transparency differentiates FW from Projected Gradient Descent (PGD) optimization. To demonstrate the utility of our theoretical framework we develop FW-AT-Adapt, a novel adversarial training algorithm which uses a simple distortion measure to adapt the number of attack steps during training to increase efficiency without compromising robustness. FW-AT-Adapt provides training time on par with single-step fast AT methods and closes the gap between fast AT methods and multi-step PGD-AT with minimal loss in adversarial accuracy in white-box and black-box settings.


title: Integer Division by Constants: Optimal Bounds

abstract: The integer division of a numerator n by a divisor d gives a quotient q and a remainder r. Optimizing compilers accelerate software by replacing the division of n by d with the division of c * n (or c * n + c) by m for convenient integers c and m chosen so that they approximate the reciprocal: c/m ~= 1/d. Such techniques are especially advantageous when m is chosen to be a power of two and when d is a constant so that c and m can be precomputed. The literature contains many bounds on the distance between c/m and the divisor d. Some of these bounds are optimally tight, while others are not. We present optimally tight bounds for quotient and remainder computations.


title: Regularity and finite element approximation for two-dimensional elliptic equations with line Dirac sources

abstract: We study the elliptic equation with a line Dirac delta function as the source term subject to the Dirichlet boundary condition in a two-dimensional domain. Such a line Dirac measure causes different types of solution singularities in the neighborhood of the line fracture. We establish new regularity results for the solution in a class of weighted Sobolev spaces and propose finite element algorithms that approximate the singular solution at the optimal convergence rate. Numerical tests are presented to justify the theoretical findings.


title: How Discrete Spectrum and Resonances Influence the Asymptotics of the Toda Shock Wave

abstract: We rigorously derive the long-time asymptotics of the Toda shock wave in a middle region where the solution is asymptotically finite gap. In particular, we describe the influence of the discrete spectrum in the spectral gap on the shift of the phase in the theta-function representation for this solution. We also study the effect of possible resonances at the endpoints of the gap on this phase. This paper is a continuation of research started in [arXiv:2001.05184].


title: Out-distribution aware Self-training in an Open World Setting

abstract: Deep Learning heavily depends on large labeled datasets which limits further improvements. While unlabeled data is available in large amounts, in particular in image recognition, it does not fulfill the closed world assumption of semi-supervised learning that all unlabeled data are task-related. The goal of this paper is to leverage unlabeled data in an open world setting to further improve prediction performance. For this purpose, we introduce out-distribution aware self-training, which includes a careful sample selection strategy based on the confidence of the classifier. While normal self-training deteriorates prediction performance, our iterative scheme improves using up to 15 times the amount of originally labeled data. Moreover, our classifiers are by design out-distribution aware and can thus distinguish task-related inputs from unrelated ones.


title: The Life and Death of SSDs and HDDs: Similarities, Differences, and Prediction Models

abstract: Data center downtime typically centers around IT equipment failure. Storage devices are the most frequently failing components in data centers. We present a comparative study of hard disk drives (HDDs) and solid state drives (SSDs) that constitute the typical storage in data centers. Using a six-year field data of 100,000 HDDs of different models from the same manufacturer from the BackBlaze dataset and a six-year field data of 30,000 SSDs of three models from a Google data center, we characterize the workload conditions that lead to failures and illustrate that their root causes differ from common expectation but remain difficult to discern. For the case of HDDs we observe that young and old drives do not present many differences in their failures. Instead, failures may be distinguished by discriminating drives based on the time spent for head positioning. For SSDs, we observe high levels of infant mortality and characterize the differences between infant and non-infant failures. We develop several machine learning failure prediction models that are shown to be surprisingly accurate, achieving high recall and low false positive rates. These models are used beyond simple prediction as they aid us to untangle the complex interaction of workload characteristics that lead to failures and identify failure root causes from monitored symptoms.