title: Van der Waals Ferromagnetic Josephson Junctions

abstract: Superconductor-ferromagnet (S-F) interfaces in two-dimensional (2D) heterostructures present a unique opportunity to study the interplay between superconductivity and ferromagnetism. The realization of such nanoscale heterostructures in van der Waals (vdW) crystals remains largely unexplored due to the challenge of making an atomically-sharp interface from their layered structures. Here, we build a vdW ferromagnetic Josephson junction (JJ) by inserting a few-layer ferromagnetic insulator Cr2Ge2Te6 into two layers of superconductor NbSe2. Owing to the remanent magnetic moment of the barrier, the critical current and the corresponding junction resistance exhibit a hysteretic and oscillatory behavior against in-plane magnetic fields, manifesting itself as a strong Josephson coupling state. Through the control of this hysteresis, we can effectively trace the magnetic properties of atomic Cr2Ge2Te6 in response to the external magnetic field. Also, we observe a central minimum of critical current in some thick JJ devices, evidencing the coexistence of 0 and {\pi} phase coupling in the junction region. Our study paves the way to exploring the sensitive probes of weak magnetism and multifunctional building blocks for phase-related superconducting circuits with the use of vdW heterostructures.


title: Perfect matching and distance spectral radius in graphs and bipartite graphs

abstract: A perfect matching in a graph $G$ is a set of nonadjacent edges covering every vertex of $G$. Motivated by recent progress on the relations between the eigenvalues and the matching number of a graph, in this paper, we aim to present a distance spectral radius condition to guarantee the existence of a perfect matching. Let $G$ be an $n$-vertex connected graph where $n$ is even and $\lambda_{1}(D(G))$ be the distance spectral radius of $G$. Then the following statements are true. \noindent$\rm{I)}$ If $4\le n\le10$ and ${\lambda }_{1} (D\left(G\right))\le {\lambda }_{1} (D(S_{n,{\frac{n}{2}}-1}))$, then $G$ contains a perfect matching unless $G\cong S_{n,{\frac{n}{2}-1}}$ where $S_{n,{\frac{n}{2}-1}}\cong K_{{\frac{n}{2}-1}}\vee ({\frac{n}{2}+1})K_1$. \noindent$\rm{II)}$ If $n\ge 12$ and ${\lambda }_{1} (D\left(G\right))\le {\lambda }_{1} (D(G^*))$, then $G$ contains a perfect matching unless $G\cong G^*$ where $G^*\cong K_1\vee (K_{n-3}\cup2K_1)$. Moreover, if $G$ is a connected $2n$-vertex balanced bipartite graph with $\lambda_{1}(D(G))\le \lambda_{1}(D(B_{n-1,n-2})) $, then $G$ contains a perfect matching, unless $G\cong B_{n-1,n-2}$ where $B_{n-1,n-2}$ is obtained from $K_{n,n-2}$ by attaching two pendent vertices to a vertex in the $n$-vertex part.


title: Testing High-latitude Curvature Effect of Gamma-Ray Bursts with {\it Fermi} Data: Evidence of Bulk Acceleration in Prompt Emission

abstract: When a gamma-ray burst (GRB) emitter stops emission abruptly, the observer receives rapidly fading emission from high latitudes with respect to the line of sight, known as the ``curvature effect''. Identifying such emission from GRB prompt-emission lightcurves would constrain the radius of prompt emission from the central engine and the composition of GRB jets. We perform a dedicated search of high-latitude emission (HLE) through spectral and temporal analyses of a sample of single-pulse bursts detected by the Gamma-ray Burst Monitor on board the {\it Fermi} satellite. We identify HLE from a subsample of bursts and constrain the emission radius to be $R_{\rm GRB} \sim (10^{15}-10^{16})$ cm from the central engine. Some bursts have the HLE decay faster than predicted by a constant Lorentz factor jet, suggesting that the emission region is undergoing acceleration during prompt emission. This supports the Poynting-flux-dominated jet composition for these bursts. The conclusion is consistent with previous results drawn from spectral-lag modeling of prompt emission and HLE analysis of X-ray flares.


title: A basis of algebraic de Rham cohomology of complete intersections over a characteristic zero field

abstract: Let $\Bbbk$ be a field of characteristic 0. Let $X$ be a smooth complete intersection over $k$ of dimension $n-k$ in the projective space $\mathbf{P}^n_{k}$, for given positive integers $n$ and $k$. When $k=\mathbb{C}$, Terasoma (\cite{Ter90}) and Konno (\cite{Ko91}) provided an explicit representative (in terms of differential forms) of a basis for the primitive middle-dimensional algebraic de Rham cohomology $H_{dR,\operatorname{prim}}^{n-k}(X;\mathbb{C})$. Later Dimca constructed another explicit representative of a basis of $H_{dR,\operatorname{prim}}^{n-k}(X;\mathbb{C})$ in \cite{Dim95}. Moreover, he proved that his representative gives the same cohomology class as the previous representative of Terasoma and Konno. The goal of this article is to examine the above two different approaches without assuming that $k=\mathbb{C}$ and provide a similar comparison result for any field $k$. Dimca's argument depends heavily on the condition $k=\mathbb{C}$ and our idea is to find appropriate Cech-de Rham complexes and spectral sequences corresponding to those two approaches, which work without restrictions on $k$.


title: Security Analysis and Improvement of Source Independent Quantum Random Number Generators with Imperfect Devices

abstract: A quantum random number generator (QRNG) as a genuine source of randomness is essential in many applications, such as number simulation and cryptography. Recently, a source-independent quantum random number generator (SI-QRNG), which can generate secure random numbers with untrusted sources, has been realized. However, the measurement loopholes of the trusted but imperfect devices used in SI-QRNGs have not yet been fully explored, which will cause security problems, especially in high-speed systems. Here, we point out and evaluate the security loopholes of practical imperfect measurement devices in SI-QRNGs. We also provide corresponding countermeasures to prevent these information leakages by recalculating the conditional minimum entropy and adding a monitor. Furthermore, by taking into account the finite-size effect,we show that the influence of the afterpulse can exceed that of the finite-size effect with the large number of sampled rounds. Our protocol is simple and effective, and it promotes the security of SI-QRNG in practice as well as the compatibility with high-speed measurement devices, thus paving the way for constructing ultrafast and security-certified commercial SI-QRNG systems.


title: Neural News Recommendation with Negative Feedback

abstract: News recommendation is important for online news services. Precise user interest modeling is critical for personalized news recommendation. Existing news recommendation methods usually rely on the implicit feedback of users like news clicks to model user interest. However, news click may not necessarily reflect user interests because users may click a news due to the attraction of its title but feel disappointed at its content. The dwell time of news reading is an important clue for user interest modeling, since short reading dwell time usually indicates low and even negative interest. Thus, incorporating the negative feedback inferred from the dwell time of news reading can improve the quality of user modeling. In this paper, we propose a neural news recommendation approach which can incorporate the implicit negative user feedback. We propose to distinguish positive and negative news clicks according to their reading dwell time, and respectively learn user representations from positive and negative news clicks via a combination of Transformer and additive attention network. In addition, we propose to compute a positive click score and a negative click score based on the relevance between candidate news representations and the user representations learned from the positive and negative news clicks. The final click score is a combination of positive and negative click scores. Besides, we propose an interactive news modeling method to consider the relatedness between title and body in news modeling. Extensive experiments on real-world dataset validate that our approach can achieve more accurate user interest modeling for news recommendation.


title: Non-Bayesian Parametric Missing-Mass Estimation

abstract: We consider the classical problem of missing-mass estimation, which deals with estimating the total probability of unseen elements in a sample. The missing-mass estimation problem has various applications in machine learning, statistics, language processing, ecology, sensor networks, and others. The naive, constrained maximum likelihood (CML) estimator is inappropriate for this problem since it tends to overestimate the probability of the observed elements. Similarly, the conventional constrained Cramer-Rao bound (CCRB), which is a lower bound on the mean-squared-error (MSE) of unbiased estimators, does not provide a relevant bound on the performance for this problem. In this paper, we introduce a frequentist, non-Bayesian parametric model of the problem of missing-mass estimation. We introduce the concept of missing-mass unbiasedness by using the Lehmann unbiasedness definition. We derive a non-Bayesian CCRB-type lower bound on the missing-mass MSE (mmMSE), named the missing-mass CCRB (mmCCRB), based on the missing-mass unbiasedness. The missing-mass unbiasedness and the proposed mmCCRB can be used to evaluate the performance of existing estimators. Based on the new mmCCRB, we propose a new method to improve existing estimators by an iterative missing-mass Fisher scoring method. Finally, we demonstrate via numerical simulations that the proposed mmCCRB is a valid and informative lower bound on the mmMSE of state-of-the-art estimators for this problem: the CML, the Good-Turing, and Laplace estimators. We also show that the performance of the Laplace estimator is improved by using the new Fisher-scoring method.


title: Explosive synchronization in interlayer phase-shifted Kuramoto oscillators on multiplex networks

abstract: We show that an introduction of a phase parameter ($\alpha$), with $0 \le \alpha \le \pi/2$, in the interlayer coupling terms of multiplex networks of Kuramoto oscillators can induce explosive synchronization (ES) in the multiplexed layers. Along with the {\alpha} values, the hysteresis width is determined by the interlayer coupling strength and the frequency mismatch between the mirror (inter-connected) nodes. A mean-field analysis is performed to support the numerical results. Similar to the earlier works, we find that the suppression of synchronization is accountable for the origin of ES. The robustness of ES against changes in the network topology and frequency distribution is tested. Finally, taking a suggestion from the synchronized state of the multiplex networks, we extend the results to the classical concept of the single-layer networks in which some specific links are assigned a phase-shifted coupling. Different methods have been introduced in the past years to incite ES in coupled oscillators; our results indicate that a phase-shifted coupling can also be one such method to achieve ES.


title: Equifocal submanifolds with non-flat section and topological Tits buildings

abstract: From the Lytchak's result for polar foliations on an irreducible simply connected symmetric space $G/K$ of compact type and rank greater than one, we can derive that there exists no equifocal submanifold with non-flat section whose codimension is greater than two in the symmetric space $G/K$. In the first-half part of this paper, we give a new proof of this non-existence theorem. The recipi of our new proof is as follows. Suppose that there exists an equifocal submanifold $M$ with non-flat section whose codimension is greater than two in an irreducible symmertric space $G/K$ of compact type and rank greater than one. We introduce the notion of a slice topology of $G/K$ associated to $M$. We consider the universal covering $\pi:\widehat{G/K}\to G/K$ of the slice topological space $G/K$ and give $\widehat{G/K}$ the manifold structure and the Riemannian metric such that $\pi$ is a Riemannian submersion onto the symmetric space $G/K$. First we show that a simplicial decomposition of the Riemannian manifold $\widehat{G/K}$ gives an irreducible topological Tits building of spherical type and rank greater than two. By applying Burns-Spatzier's theorem to this topological Tits building, we show that the Riemannian manifold $\widehat{G/K}$ is homothetic to the unit sphere. Furthermore, from this fact, we show that $G/K$ is isometric to a sphere, a complex projective space or a quaternionic projective space. This contradicts that $G/K$ is of rank greater than one. This is the recipi of our proof. In the second-half part, we estimate the codimension of $M$ from above by using the multiplicities of the roots of the root system of $G/K$. As its result, we can show that there exists no equifocal submanifold with non-flat section in some irreducible simply connected symmetric spaces of compact type.


title: The Lattice Sine-Gordon Equation as a Superposition Formula for an NLS-Type System

abstract: We treat the lattice sine-Gordon equation and two of its generalised symmetries as a compatible system. Elimination of shifts from the two symmetries of the lattice sine-Gordon equation yields an integrable NLS-type system. An auto-B\"acklund transformation and a superposition formula for the NLS-type system is obtained by elimination of shifts from the lattice sine-Gordon equation and its down-shifted version. We use the obtained formulae to calculate a superposition of two and three elementary solutions.


title: Seed Stocking Via Multi-Task Learning

abstract: Sellers of crop seeds need to plan for the variety and quantity of seeds to stock at least a year in advance. There are a large number of seed varieties of one crop, and each can perform best under different growing conditions. Given the unpredictability of weather, farmers need to make decisions that balance high yield and low risk. A seed vendor needs to be able to anticipate the needs of farmers and have them ready. In this study, we propose an analytical framework for estimating seed demand with three major steps. First, we will estimate the yield and risk of each variety as if they were planted at each location. Since past experiments performed with different seed varieties are highly unbalanced across varieties, and the combination of growing conditions is sparse, we employ multi-task learning to borrow information from similar varieties. Second, we will determine the best mix of seeds for each location by seeking a tradeoff between yield and risk. Third, we will aggregate such mix and pick the top five varieties to re-balance the yield and risk for each growing location. We find that multi-task learning provides a viable solution for yield prediction, and our overall analytical framework has resulted in a good performance.


title: Change-point detection using spectral PCA for multivariate time series

abstract: We propose a two-stage approach Spec PC-CP to identify change points in multivariate time series. In the first stage, we obtain a low-dimensional summary of the high-dimensional time series by Spectral Principal Component Analysis (Spec-PCA). In the second stage, we apply cumulative sum-type test on the Spectral PCA component using a binary segmentation algorithm. Compared with existing approaches, the proposed method is able to capture the lead-lag relationship in time series. Our simulations demonstrate that the Spec PC-CP method performs significantly better than competing methods for detecting change points in high-dimensional time series. The results on epileptic seizure EEG data and stock data also indicate that our new method can efficiently {detect} change points corresponding to the onset of the underlying events.


title: Panorama: A Framework to Support Collaborative Context Monitoring on Co-Located Mobile Devices

abstract: A key challenge in wide adoption of sophisticated context-aware applications is the requirement of continuous sensing and context computing. This paper presents Panorama, a middleware that identifies collaboration opportunities to offload context computing tasks to nearby mobile devices as well as cloudlets/cloud. At the heart of Panorama is a multi-objective optimizer that takes into account different constraints such as access cost, computation capability, access latency, energy consumption and data privacy, and efficiently computes a collaboration plan optimized simultaneously for different objectives such as minimizing cost, energy and/or execution time. Panorama provides support for discovering nearby devices and cloudlets/cloud, computing an optimal collaboration plan, distributing computation to participating devices, and getting the results back. The paper provides an extensive evaluation of Panorama via two representative context monitoring applications over a set of Android devices and a cloudlet/cloud under different constraints.


title: Fluid pressurisation and earthquake propagation in the Hikurangi subduction zone

abstract: In subduction zones, seismic slip at shallow crustal depths can lead to the generation of tsunamis. Large slip displacements during tsunamogenic earthquakes are attributed to the low coseismic shear strength of the fluid-saturated and non-lithified clay-rich fault rocks. However, because of experimental challenges in confining these materials, the physical processes responsible of the coseismic reduction in fault shear strength are poorly understood. Using a novel experimental setup, we measured pore fluid pressure during simulated seismic slip in clay-rich materials sampled from the deep oceanic drilling of the P\=apaku thrust (Hikurangi subduction zone, New Zealand). Here we show that at seismic velocity, shear-induced dilatancy is followed by pressurisation of fluids. The thermal and mechanical pressurisation of fluids, enhanced by the low permeability of the fault, reduces the energy required to propagate earthquake rupture. We suggest that fluid-saturated clay-rich sediments, occurring at shallow depth in subduction zones, can promote earthquake rupture propagation and slip because of their low permeability and tendency to pressurise when sheared at seismic slip velocities.


title: Blind Modulation Classification via Combined Machine Learning and Signal Feature Extraction

abstract: In this study, an algorithm to blind and automatic modulation classification has been proposed. It well benefits combined machine leaning and signal feature extraction to recognize diverse range of modulation in low signal power to noise ratio (SNR). The presented algorithm contains four. First, it advantages spectrum analyzing to branching modulated signal based on regular and irregular spectrum character. Seconds, a nonlinear soft margin support vector (NS SVM) problem is applied to received signal, and its symbols are classified to correct and incorrect (support vectors) symbols. The NS SVM employment leads to discounting in physical layer noise effect on modulated signal. After that, a k-center clustering can find center of each class. finally, in correlation function estimation of scatter diagram is correlated with pre-saved ideal scatter diagram of modulations. The correlation outcome is classification result. For more evaluation, success rate, performance, and complexity in compare to many published methods are provided. The simulation prove that the proposed algorithm can classified the modulated signal in less SNR. For example, it can recognize 4-QAM in SNR=-4.2 dB, and 4-FSK in SNR=2.1 dB with %99 success rate. Moreover, due to using of kernel function in dual problem of NS SVM and feature base function, the proposed algorithm has low complexity and simple implementation in practical issues.


title: A Survey of Privacy-Preserving Techniques for Encrypted Traffic Inspection over Network Middleboxes

abstract: Middleboxes in a computer network system inspect and analyse network traffic to detect malicious communications, monitor system performance and provide operational services. However, encrypted traffic hinders the ability of middleboxes to perform such services. A common practice in addressing this issue is by employing a "Man-in-the-Middle" (MitM) approach, wherein an encrypted traffic flow between two endpoints is interrupted, decrypted and analysed by the middleboxes. The MitM approach is straightforward and is used by many organisations, but there are both practical and privacy concerns. Due to the cost of the MitM appliances and the latency incurred in the encrypt-decrypt processes, enterprises continue to seek solutions that are less costly. There were discussion on the many efforts required to configure MitM. Besides, MitM violates end-to-end privacy guarantee, raising privacy concerns and issues on compliance especially with the rising awareness on user privacy. Furthermore, some of the MitM implementations were found to be flawed. Consequently, new practical and privacy-preserving techniques for inspection over encrypted traffic were proposed. We examine them to compare their advantages, limitations and challenges. We categorise them into four main categories by defining a framework that consist of system architectures, use cases, trust and threat models. These are searchable encryption, access control, machine learning and trusted hardware. We first discuss the man-in-the-middle approach as a baseline, then discuss in details each of them, and provide an in-depth comparisons of their advantages and limitations. By doing so we describe practical constraints, advantages and pitfalls towards adopting the techniques. We also give insights on the gaps between research work and industrial deployment, which leads us to the discussion on the challenges and research directions.


title: Locality Sensitive Hashing for Efficient Similar Polygon Retrieval

abstract: Locality Sensitive Hashing (LSH) is an effective method of indexing a set of items to support efficient nearest neighbors queries in high-dimensional spaces. The basic idea of LSH is that similar items should produce hash collisions with higher probability than dissimilar items. We study LSH for (not necessarily convex) polygons, and use it to give efficient data structures for similar shape retrieval. Arkin et al. represent polygons by their "turning function" - a function which follows the angle between the polygon's tangent and the $ x $-axis while traversing the perimeter of the polygon. They define the distance between polygons to be variations of the $ L_p $ (for $p=1,2$) distance between their turning functions. This metric is invariant under translation, rotation and scaling (and the selection of the initial point on the perimeter) and therefore models well the intuitive notion of shape resemblance. We develop and analyze LSH near neighbor data structures for several variations of the $ L_p $ distance for functions (for $p=1,2$). By applying our schemes to the turning functions of a collection of polygons we obtain efficient near neighbor LSH-based structures for polygons. To tune our structures to turning functions of polygons, we prove some new properties of these turning functions that may be of independent interest. As part of our analysis, we address the following problem which is of independent interest. Find the vertical translation of a function $ f $ that is closest in $ L_1 $ distance to a function $ g $. We prove tight bounds on the approximation guarantee obtained by the translation which is equal to the difference between the averages of $ g $ and $ f $.


title: Take More Positives: An Empirical Study of Contrastive Learing in Unsupervised Person Re-Identification

abstract: Unsupervised person re-identification (re-ID) aims at closing the performance gap to supervised methods. These methods build reliable relationship between data points while learning representations. However, we empirically show that the reason why they are successful is not only their label generation mechanisms, but also their unexplored designs. By studying two unsupervised person re-ID methods in a cross-method way, we point out a hard negative problem is handled implicitly by their designs of data augmentations and PK sampler respectively. In this paper, we find another simple solution for the problem, i.e., taking more positives during training, by which we generate pseudo-labels and update models in an iterative manner. Based on our findings, we propose a contrastive learning method without a memory back for unsupervised person re-ID. Our method works well on benchmark datasets and outperforms the state-of-the-art methods. Code will be made available.


title: Tunable plasmons in large area WTe2 thin films

abstract: The observation of the electrically tunable and highly confined plasmons in graphene has stimulated the exploration of interesting properties of plasmons in other two dimensional materials. Recently, hyperbolic plasmon resonance modes are observed in exfoliated WTe2 films, a type-II Weyl semimetal with layered structure, providing a platform for the assembly of plasmons with hyperbolicity and exotic topological properties. However, the plasmon modes were observed in relatively thick and small-area films, which restrict the tunability and application for plasmons. Here, large-area (~ cm) WTe2 films with different thickness are grown by chemical vapor deposition method, in which plasmon resonance modes are observed in films with different thickness down to about 8 nm. Hybridization of plasmon and surface polar phonons of the substrate is revealed by mapping the plasmon dispersion. The plasmon frequency is demonstrated to be tunable by changing the temperature and film thickness. Our results facilitate the development of a tunable and scalable WTe2 plasmonic system for revealing topological properties and towards various applications in sensing, imaging and light modulation.


title: Mixup Without Hesitation

abstract: Mixup linearly interpolates pairs of examples to form new samples, which is easy to implement and has been shown to be effective in image classification tasks. However, there are two drawbacks in mixup: one is that more training epochs are needed to obtain a well-trained model; the other is that mixup requires tuning a hyper-parameter to gain appropriate capacity but that is a difficult task. In this paper, we find that mixup constantly explores the representation space, and inspired by the exploration-exploitation dilemma in reinforcement learning, we propose mixup Without hesitation (mWh), a concise, effective, and easy-to-use training algorithm. We show that mWh strikes a good balance between exploration and exploitation by gradually replacing mixup with basic data augmentation. It can achieve a strong baseline with less training time than original mixup and without searching for optimal hyper-parameter, i.e., mWh acts as mixup without hesitation. mWh can also transfer to CutMix, and gain consistent improvement on other machine learning and computer vision tasks such as object detection. Our code is open-source and available at https://github.com/yuhao318/mwh


title: Parallelized POD-based Suboptimal Economic Model Predictive Control of a State-Constrained Boussinesq approximation

abstract: Motivated by an energy efficient building application, we want to optimize a quadratic cost functional subject to the Boussinesq approximation of the Navier-Stokes equations and to bilateral state and control constraints. Since the computation of such an optimal solution is numerically costly, we design an efficient strategy to compute a sub-optimal (but applicationally acceptable) solution with significantly reduced computational effort. We employ an economic Model Predictive Control (MPC) strategy to obtain a feedback control. The MPC sub-problems are based on a linear-quadratic optimal control problem subjected to mixed control and state constraints and a convection-diffusion equation, reduced with proper orthogonal decomposition. To solve each sub-problem, we apply a primal-dual active set strategy. The method can be fully parallelized, which enables the solution of large problems with real-world parameters.


title: Invertibility In the Sense of Ehrenpreis

abstract: We consider zero sets of entire functions belonging to the Schwartz algebra. This algebra is defined as the Fourier-Laplace transform image of the space of all distributions compactly supported on the real line. We study the conditions under which given complex sequence forms zero set of some invertible in the sense of Ehrenpreis element of the Schwartz algebra (slowly decreasing function).


title: From Control to Mathematics-Part I: Controllability-Based Design for Iterative Methods in Solving Linear Equations

abstract: In the interaction between control and mathematics, mathematical tools are fundamental for all the control methods, but it is unclear how control impacts mathematics. This is the first part of our paper that attempts to give an answer with focus on solving linear algebraic equations (LAEs) from the perspective of systems and control, where it mainly introduces the controllability-based design results. By proposing an iterative method that integrates a learning control mechanism, a class of tracking problems for iterative learning control (ILC) is explored for the problem solving of LAEs. A trackability property of ILC is newly developed, by which analysis and synthesis results are established to disclose the equivalence between the solvability of LAEs and the controllability of discrete control systems. Hence, LAEs can be solved by equivalently achieving the perfect tracking tasks of resulting ILC systems via the classic state feedback-based design and analysis methods. It is shown that the solutions for any solvable LAE can all be calculated with different selections of the initial input. Moreover, the presented ILC method is applicable to determining all the least squares solutions of any unsolvable LAE. In particular, a deadbeat design is incorporated to ILC such that the solving of LAEs can be completed within finite iteration steps. The trackability property is also generalized to conventional two-dimensional ILC systems, which creates feedback-based methods, instead of the common used contraction mapping-based methods, for the design and convergence analysis of ILC.


title: Empirical Decomposition of the IV-OLS Gap with Heterogeneous and Nonlinear Effects

abstract: This study proposes an econometric framework to interpret and empirically decompose the difference between IV and OLS estimates given by a linear regression model when the true causal effects of the treatment are nonlinear in treatment levels and heterogeneous across covariates. I show that the IV-OLS coefficient gap consists of three estimable components: the difference in weights on the covariates, the difference in weights on the treatment levels, and the difference in identified marginal effects that arises from endogeneity bias. Applications of this framework to return-to-schooling estimates demonstrate the empirical relevance of this distinction in properly interpreting the IV-OLS gap.


title: Proceedings of the NeurIPS 2020 Workshop on Machine Learning for the Developing World: Improving Resilience

abstract: These are the proceedings of the 4th workshop on Machine Learning for the Developing World (ML4D), held as part of the Thirty-fourth Conference on Neural Information Processing Systems (NeurIPS) on Saturday, December 12th 2020.


title: Phase Retrieval using Expectation Consistent Signal Recovery Algorithm based on Hypernetwork

abstract: Phase retrieval (PR) is an important component in modern computational imaging systems. Many algorithms have been developed over the past half-century. Recent advances in deep learning have introduced new possibilities for a robust and fast PR. An emerging technique called deep unfolding provides a systematic connection between conventional model-based iterative algorithms and modern data-based deep learning. Unfolded algorithms, which are powered by data learning, have shown remarkable performance and convergence speed improvement over original algorithms. Despite their potential, most existing unfolded algorithms are strictly confined to a fixed number of iterations when layer-dependent parameters are used. In this study, we develop a novel framework for deep unfolding to overcome existing limitations. Our development is based on an unfolded generalized expectation consistent signal recovery (GEC-SR) algorithm, wherein damping factors are left for data-driven learning. In particular, we introduce a hypernetwork to generate the damping factors for GEC-SR. Instead of learning a set of optimal damping factors directly, the hypernetwork learns how to generate the optimal damping factors according to the clinical settings, thereby ensuring its adaptivity to different scenarios. To enable the hypernetwork to adapt to varying layer numbers, we use a recurrent architecture to develop a dynamic hypernetwork that generates a damping factor that can vary online across layers. We also exploit a self-attention mechanism to enhance the robustness of the hypernetwork. Extensive experiments show that the proposed algorithm outperforms existing ones in terms of convergence speed and accuracy and still works well under very harsh settings, even under which many classical PR algorithms are unstable.


title: A continuous generalization of domination-like invariants

abstract: In this paper, we define a new domination-like invariant of graphs. Let $\mathbb{R}^{+}$ be the set of non-negative numbers. Let $c\in \mathbb{R}^{+}-\{0\}$ be a number, and let $G$ be a graph. A function $f:V(G)\rightarrow \mathbb{R}^{+}$ is a $c$-self-dominating function of $G$ if for every $u\in V(G)$, $f(u)\geq c$ or $\max\{f(v):v\in N_{G}(u)\}\geq 1$. The $c$-self-domination number $\gamma ^{c}(G)$ of $G$ is defined as $\gamma ^{c}(G):=\min\{\sum_{u\in V(G)}f(u):f$ is a $c$-self-dominating function of $G\}$. Then $\gamma ^{1}(G)$, $\gamma ^{\infty }(G)$ and $\gamma ^{\frac{1}{2}}(G)$ are equal to the domination number, the total domination number and the half of the Roman domination number of $G$, respectively. Our main aim is to continuously fill in the gaps among such three invariants. In this paper, we give a sharp upper bound of the $c$-self-domination number for all $c\geq \frac{1}{2}$.


title: Automated Detection of Patellofemoral Osteoarthritis from Knee Lateral View Radiographs Using Deep Learning: Data from the Multicenter Osteoarthritis Study (MOST)

abstract: Objective: To assess the ability of imaging-based deep learning to predict radiographic patellofemoral osteoarthritis (PFOA) from knee lateral view radiographs. Design: Knee lateral view radiographs were extracted from The Multicenter Osteoarthritis Study (MOST) (n = 18,436 knees). Patellar region-of-interest (ROI) was first automatically detected, and subsequently, end-to-end deep convolutional neural networks (CNNs) were trained and validated to detect the status of patellofemoral OA. Patellar ROI was detected using deep-learning-based object detection method. Manual PFOA status assessment provided in the MOST dataset was used as a classification outcome for the CNNs. Performance of prediction models was assessed using the area under the receiver operating characteristic curve (ROC AUC) and the average precision (AP) obtained from the precision-recall (PR) curve in the stratified 5-fold cross validation setting. Results: Of the 18,436 knees, 3,425 (19%) had PFOA. AUC and AP for the reference model including age, sex, body mass index (BMI), the total Western Ontario and McMaster Universities Arthritis Index (WOMAC) score, and tibiofemoral Kellgren-Lawrence (KL) grade to predict PFOA were 0.806 and 0.478, respectively. The CNN model that used only image data significantly improved the prediction of PFOA status (ROC AUC= 0.958, AP= 0.862). Conclusion: We present the first machine learning based automatic PFOA detection method. Furthermore, our deep learning based model trained on patella region from knee lateral view radiographs performs better at predicting PFOA than models based on patient characteristics and clinical assessments.


title: The Beta-Mixture Shrinkage Prior for Sparse Covariances with Posterior Minimax Rates

abstract: Statistical inference for sparse covariance matrices is crucial to reveal dependence structure of large multivariate data sets, but lacks scalable and theoretically supported Bayesian methods. In this paper, we propose beta-mixture shrinkage prior, computationally more efficient than the spike and slab prior, for sparse covariance matrices and establish its minimax optimality in high-dimensional settings. The proposed prior consists of beta-mixture shrinkage and gamma priors for off-diagonal and diagonal entries, respectively. To ensure positive definiteness of the resulting covariance matrix, we further restrict the support of the prior to a subspace of positive definite matrices. We obtain the posterior convergence rate of the induced posterior under the Frobenius norm and establish a minimax lower bound for sparse covariance matrices. The class of sparse covariance matrices for the minimax lower bound considered in this paper is controlled by the number of nonzero off-diagonal elements and has more intuitive appeal than those appeared in the literature. The obtained posterior convergence rate coincides with the minimax lower bound unless the true covariance matrix is extremely sparse. In the simulation study, we show that the proposed method is computationally more efficient than competitors, while achieving comparable performance. Advantages of the shrinkage prior are demonstrated based on two real data sets.


title: The free energy of spherical pure $p$-spin models -- computation from the TAP approach

abstract: We compute the free energy at all temperatures for the spherical pure $p$-spin models from the generalized Thouless-Anderson-Palmer representation. This is the first example of a mixed $p$-spin model for which the free energy is computed in the whole replica symmetry breaking phase, without appealing to the famous Parisi formula.


title: Exponential convergence of distributed optimization for heterogeneous linear multi-agent systems

abstract: In this work we study a distributed optimal output consensus problem for heterogeneous linear multi-agent systems where the agents aim to reach consensus with the purpose of minimizing the sum of private convex costs. Based on output feedback, a fully distributed control law is proposed by using the proportional-integral (PI) control technique. For strongly convex cost functions with Lipschitz gradients, the designed controller can achieve convergence exponentially in an undirected and connected network. Furthermore, to remove the requirement of continuous communications, the proposed control law is then extended to periodic and event-triggered communication schemes, which also achieve convergence exponentially. Two simulation examples are given to verify the proposed control algorithms.


title: Activation Density based Mixed-Precision Quantization for Energy Efficient Neural Networks

abstract: As neural networks gain widespread adoption in embedded devices, there is a need for model compression techniques to facilitate deployment in resource-constrained environments. Quantization is one of the go-to methods yielding state-of-the-art model compression. Most approaches take a fully trained model, apply different heuristics to determine the optimal bit-precision for different layers of the network, and retrain the network to regain any drop in accuracy. Based on Activation Density (AD)-the proportion of non-zero activations in a layer-we propose an in-training quantization method. Our method calculates bit-width for each layer during training yielding a mixed precision model with competitive accuracy. Since we train lower precision models during training, our approach yields the final quantized model at lower training complexity and also eliminates the need for re-training. We run experiments on benchmark datasets like CIFAR-10, CIFAR-100, TinyImagenet on VGG19/ResNet18 architectures and report the accuracy and energy estimates for the same. We achieve ~4.5x benefit in terms of estimated multiply-and-accumulate (MAC) reduction while reducing the training complexity by 50% in our experiments. To further evaluate the energy benefits of our proposed method, we develop a mixed-precision scalable Process In Memory (PIM) hardware accelerator platform. The hardware platform incorporates shift-add functionality for handling multi-bit precision neural network models. Evaluating the quantized models obtained with our proposed method on the PIM platform yields ~5x energy reduction compared to 16-bit models. Additionally, we find that integrating AD based quantization with AD based pruning (both conducted during training) yields up to ~198x and ~44x energy reductions for VGG19 and ResNet18 architectures respectively on PIM platform compared to baseline 16-bit precision, unpruned models.


title: Neural Contract Element Extraction Revisited: Letters from Sesame Street

abstract: We investigate contract element extraction. We show that LSTM-based encoders perform better than dilated CNNs, Transformers, and BERT in this task. We also find that domain-specific WORD2VEC embeddings outperform generic pre-trained GLOVE embeddings. Morpho-syntactic features in the form of POS tag and token shape embeddings, as well as context-aware ELMO embeddings do not improve performance. Several of these observations contradict choices or findings of previous work on contract element extraction and generic sequence labeling tasks, indicating that contract element extraction requires careful task-specific choices. Analyzing the results of (i) plain TRANSFORMER-based and (ii) BERT-based models, we find that in the examined task, where the entities are highly context-sensitive, the lack of recurrency in TRANSFORMERs greatly affects their performance.


title: On the Calibration and Uncertainty of Neural Learning to Rank Models

abstract: According to the Probability Ranking Principle (PRP), ranking documents in decreasing order of their probability of relevance leads to an optimal document ranking for ad-hoc retrieval. The PRP holds when two conditions are met: [C1] the models are well calibrated, and, [C2] the probabilities of relevance are reported with certainty. We know however that deep neural networks (DNNs) are often not well calibrated and have several sources of uncertainty, and thus [C1] and [C2] might not be satisfied by neural rankers. Given the success of neural Learning to Rank (L2R) approaches-and here, especially BERT-based approaches-we first analyze under which circumstances deterministic, i.e. outputs point estimates, neural rankers are calibrated. Then, motivated by our findings we use two techniques to model the uncertainty of neural rankers leading to the proposed stochastic rankers, which output a predictive distribution of relevance as opposed to point estimates. Our experimental results on the ad-hoc retrieval task of conversation response ranking reveal that (i) BERT-based rankers are not robustly calibrated and that stochastic BERT-based rankers yield better calibration; and (ii) uncertainty estimation is beneficial for both risk-aware neural ranking, i.e.taking into account the uncertainty when ranking documents, and for predicting unanswerable conversational contexts.


title: On the Differential Private Data Market: Endogenous Evolution, Dynamic Pricing, and Incentive Compatibility

abstract: Privacy is an essential issue in data trading markets. This work uses a mechanism design approach to study the optimal market model to economize the value of privacy of personal data, using differential privacy. The buyer uses a finite number of randomized algorithms to get access to the owners' data in a sequential-composition manner, in which each randomized algorithm is differentially private. Each usage of a randomized algorithm is referred to as a period. Due to the composability of differential privacy, there are inevitable privacy losses accumulated over periods. Hence, we allow the owners to leave the market at the end of any period by making stopping decisions. We define an instrumental kernel function to capture the instrumentalness of owners' preferences and model the formation of each owner's (both intrinsic and instrumental) preference over periods by taking into consideration the composability of differential privacy and the time-varying nature of privacy concerns. Our desideratum is to study the buyer's design regime of optimal market models in a dynamic environment when each owner makes coupled decisions of stopping and reporting of their preferences. The buyer seeks to design a privacy allocation rule that dynamically specifies the degree of privacy protections and a payment rule to compensate for the privacy losses of the owners. The buyer additionally chooses a payment rule which is independent of owners' reports of their preferences to influence the owners' stopping decisions. We characterize the dynamic incentive compatibility and provide a design principle to construct the payment rules in terms of the privacy allocation rule. Further, we relax the buyer's market design problem and provide a sufficient condition for an approximated dynamic incentive-compatible market model.


title: Alignment of electron optical beam shaping elements using a convolutional neural network

abstract: A convolutional neural network is used to align an orbital angular momentum sorter in a transmission electron microscope. The method is demonstrated using simulations and experimentally. As a result of its accuracy and speed, it offers the possibility of real-time tuning of other electron optical devices and electron beam shaping configurations.


title: Dirac Metamaterial Assembled by Pyrene Derivative and its Topological Photonics

abstract: Over the past decade, topology has garnered great attention in a wide area of physics. In particular, it has exerted influence on photonics because carefully engineered photonic crystals and metamaterials can help explore the non-trivial state of materials. In this regard, all dielectric metamaterials with large anisotropy, and dipole and multipole Mie resonators have played an increasingly important role in topological photonics. Advantages of Mie resonators make it possible to quest for non-trivial states in three dimensions and theoretical calculation supports its potential. However, it is very difficult to demonstrate this experimentally because it is hard to make the metacrystal by anisotropic meta-atoms despite much effort. Here we report a Dirac metamaterial for 3D topological photonics. It is implemented by a metacrystal self-assembled by a molecule, HYLION-12 which has both anisotropic polarizability and ring current. As its peculiar properties, it has an exotic optical constant that can be used for the electric and magnetic hyperbolic metamaterial, and the double hyperbolic metamaterial in the ultraviolet region. It also showed 142% of reflectance at 242nm as an amplified reflector and asymmetric transmittance up to 30% through the opaque substrate as a Huygens source under 300nm. Furthermore, it demonstrated various phenomena of topological photonics such as Pancharatnam-Berry and waveguide phase merging, wavefront shaping and waveguide on edges as a 3D topological photonic material. The new strategy using polyaromatic hydrocarbons (PAHs) is expected to be an effective way to realize 3D topological photonics.

categories: physics.optics cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.other


title: Asymptotic Scattering by Poissonian Thermostats

abstract: We consider an infinite chain of coupled harmonic oscillators with a Poisson thermostat at the origin. In the high frequency limit, we establish the reflection-transmission-scattering coefficients for the wave energy scattered off the thermostat. Unlike the case of the Langevin thermostat [5], in the macroscopic limit the Poissonian thermostat scattering generates a continuous cloud of waves of frequencies different from that of the incident wave.


title: Statistical-physics-inspired model for intrinsic fluctuations driving supply and demand in markets

abstract: We propose a simple statistical-physics-inspired model for the effect of intrinsic fluctuations on supply and demand in markets. The model consists of agents that trade in two types of goods of which the total number is separately conserved. The relative preference of an individual agent for the two types of goods is determined by a utility that is identical for all agents. Market supply and demand curves are computed and compared for various motivated choices of the distribution of goods over the agents. In particular, we compare the "mean-field" case, in which all agents have the same number of goods and that is akin to the economics textbook case, to the case of Boltzmann-Gibbs distributed goods, in which agents have a fluctuating number of goods. We find that the resulting equilibrium prices are not equal for these two approaches, especially when a large fraction of the agents can neither buy nor sell.


title: Density matrix renormalization group approach to the low temperature thermodynamics of correlated 1D fermionic models

abstract: The low temperature thermodynamics of correlated 1D fermionic models with spin and charge degrees of freedom is obtained by exact diagonalization (ED) of small systems and followed by density matrix renormalization group (DMRG) calculations that target the lowest hundreds of states $\{E(N)\}$ at system size $N$ instead of the ground state. Progressively larger $N$ reaches $T < 0.05t$ in correlated models with electron transfer $t$ between first neighbors and bandwidth $4t$. The size dependence of the many-fermion basis is explicitly included for arbitrary interactions by scaling the partition function. The remaining size dependence is then entirely due to the energy spectrum $\{E(N)\}$ of the model. The ED/DMRG method is applied to Hubbard and extended Hubbard models, both gapped and gapless, with $N_e = N$ or $N/2$ electrons and is validated against exact results for the magnetic susceptibility $\chi(T)$ and entropy $S(T)$ per site. Some limitations of the method are noted. Special attention is given to the bond-order-wave phase of the extended Hubbard model with competing interactions and low $T$ thermodynamics sensitive to small gaps.


title: The Daily Life of Software Engineers during the COVID-19 Pandemic

abstract: Following the onset of the COVID-19 pandemic and subsequent lockdowns, software engineers' daily life was disrupted and abruptly forced into remote working from home. This change deeply impacted typical working routines, affecting both well-being and productivity. Moreover, this pandemic will have long-lasting effects in the software industry, with several tech companies allowing their employees to work from home indefinitely if they wish to do so. Therefore, it is crucial to analyze and understand how a typical working day looks like when working from home and how individual activities affect software developers' well-being and productivity. We performed a two-wave longitudinal study involving almost 200 globally carefully selected software professionals, inferring daily activities with perceived well-being, productivity, and other relevant psychological and social variables. Results suggest that the time software engineers spent doing specific activities from home was similar when working in the office. However, we also found some significant mean differences. The amount of time developers spent on each activity was unrelated to their well-being, perceived productivity, and other variables. We conclude that working remotely is not per se a challenge for organizations or developers.


title: Impact of non-gravitational effects on chaotic properties of retrograde orbits

abstract: Dynamical studies of asteroid populations in retrograde orbits, that is with orbital inclinations greater than 90 degrees, are interesting because the origin of such orbits is still unexplained. Generally, the population of retrograde asteroids includes mostly Centaurs and transneptunian objects (TNOs). A special case is the near-Earth object (343158) 2009 HC82 from the Apollo group. Another interesting object is the comet 333P/LINEAR, which for several years was considered the second retrograde object approaching Earth. Another comet in retrograde orbit, 161P Hartley/IRAS appears to be an object of similar type. Thanks to the large amount of observational data for these two comets, we tested various models of cometary non-gravitational forces applied to their dynamics. The goal was to estimate which of non-gravitational perturbations could affect the stability of retrograde bodies. In principle, we study the local stability by measuring the divergence of nearby orbits. We numerically determined Lyapunov characteristic indicators (LCI) and the associated Lyapunov times (LT). This time, our calculations were extended by more advanced models of non-gravitational perturbations (i.e. Yarkovsky drift and in selected cases cometary forces). This allowed us to estimate chaos in the Lyapunov sense. We found that the Yarkovsky effect for obliquities of $\gamma=0^{\circ}$ and $\gamma=180^{\circ}$ can change the LT substantially. In most cases, for the prograde rotation, we received more stable solutions. Moreover, we confirmed the role of retrograde resonances in this process. Additionally, the studied cometary effects also significantly influence the long-term behaviour of the selected comets. The LT can reach values from 100 to over 1000 years. Conclusions. All of our results indicate that the use of models with non-gravitational effects for retrograde bodies is clearly justified.


title: Event-Driven Source Traffic Prediction in Machine-Type Communications Using LSTM Networks

abstract: Source traffic prediction is one of the main challenges of enabling predictive resource allocation in machine type communications (MTC). In this paper, a Long Short-Term Memory (LSTM) based deep learning approach is proposed for event-driven source traffic prediction. The source traffic prediction problem can be formulated as a sequence generation task where the main focus is predicting the transmission states of machine-type devices (MTDs) based on their past transmission data. This is done by restructuring the transmission data in a way that the LSTM network can identify the causal relationship between the devices. Knowledge of such a causal relationship can enable event-driven traffic prediction. The performance of the proposed approach is studied using data regarding events from MTDs with different ranges of entropy. Our model outperforms existing baseline solutions in saving resources and accuracy with a margin of around 9%. Reduction in Random Access (RA) requests by our model is also analyzed to demonstrate the low amount of signaling required as a result of our proposed LSTM based source traffic prediction approach.


title: Atomic Layer Deposition of Aluminum (111) Thin Film by Dimethylethylaminealane Precursor

abstract: We report the growth of aluminum (111) thin film by atomic layer deposition (ALD) technique with dimethylethylaminealane (DMEAA) as a precursor. It is found that the metallic underlayer is essential to grow uniform aluminum films by DMEAA precursor. As a titanium thin film is used as the underlayer, grown aluminum thin film shows (111) orientation irrespective of substrates. The lattice constant and superconducting transition temperature of the aluminum thin films are the same as the bulk one. These findings suggest that ALD technique provides high quality of the aluminum thin films and have potential for the applications of superconducting devices. We discuss ALD technique with DMEAA precursor is the promising method for fabricating vertical small Josephson tunnel junctions, which can be used as the superconducting quantum bits.


title: Sensitivity-conditioning: Beyond Singular Perturbation for Control Design on Multiple Time Scales

abstract: A classical approach to design controllers for interconnected systems is to assume that the different subsystems operate at different time scales, then design simpler controllers within each time scale, and finally certify stability of the interconnected system via singular perturbation analysis. In this work, we propose an alternative approach that also allows to design the controllers of the individual subsystems separately. However, instead of requiring a sufficiently large time-scale separation, our approach consists of adding a feed-forward term to modify the dynamics of faster systems in order to anticipate the dynamics of slower ones. We present several examples in bilevel optimization and cascade control design, where our approach improves the performance of currently available methods.


title: Riemannian manifolds with entire Grauert tube are rationally elliptic

abstract: It was conjectured by Bott-Grove-Halperin that a compact simply connected Riemannian manifold $M$ with nonnegative sectional curvature is rationally elliptic. We confirm this conjecture under the stronger assumption that $M$ has entire Grauert tube,i.e.,$M$ is a real analytic Riemannian manifold that has a unique adapted complex structure defined on the whole tangent bundle $TM$.


title: Wake fields and impedance of space charge in Cartesian coordinate system

abstract: In Cartesian coordinate system, the fields generated by a point charge moving parallel to the axis of a rectangular vacuum chamber can be formulated in terms of eigenfunctions of the rectangular waveguide using the mode expansion method. In combination with the conventional impedance theory, the Green-function forms of the wake functions and impedance for space-charge effects can be obtained, and are found to be functions of the positions of both the source and test particles. Using the Green's functions calculated for a point charge, the wake fields and impedance of a beam with various distributions can be calculated, and should be useful to model the three-dimensional space-charge effects. This paper summarizes our findings and also show comparisons to the existing theories.


title: Micropillar compression of single crystal tungsten carbide, Part 2: Lattice rotation axis to identify deformation slip mechanisms

abstract: The plastic deformation mechanisms of tungsten carbide at room and elevated temperatures influence the wear and fracture properties of WC-Co hardmetal composite materials. The relationship between residual defect structures, including glissile and sessile dislocations and stacking faults, and the slip deformation activity, which produce slip traces, is not clear. Part 1 of this study showed that {10-10} was the primary slip plane at all measured temperatures and orientations, but secondary slip on the basal plane was activated at 600 {\deg}C, which suggests that <a> dislocations can cross-slip onto the basal plane at 600 {\deg}C. In the present work, Part 2, lattice rotation axis analysis of deformed WC micropillar mid-sections was used to discriminate <a> prismatic slip from multiple <c+a> prismatic slip in WC, which enabled the dislocation types contributing to plastic slip to be distinguished, independently of TEM residual defect analysis. Prismatic-oriented micropillars deformed primarily by multiple <c+a> prismatic slip at room temperature, but by <a> prismatic slip at 600 {\deg}C. Deformation in the near-basal oriented pillar at 600 {\deg}C can be modelled as prismatic slip along <c> constrained by the indenter face and pillar base. Secondary <a> basal slip, which was observed near the top of the pillar, was activated to maintain deformation compatibility with the indenter face. The lattice rotations, buckled pillar shape, mechanical data, and slip traces observed in the pillar are all consistent with this model.


title: Data-driven peakon and periodic peakon travelling wave solutions of some nonlinear dispersive equations via deep learning

abstract: In the field of mathematical physics, there exist many physically interesting nonlinear dispersive equations with peakon solutions, which are solitary waves with discontinuous first-order derivative at the wave peak. In this paper, we apply the multi-layer physics-informed neural networks (PINNs) deep learning to successfully study the data-driven peakon and periodic peakon solutions of some well-known nonlinear dispersion equations with initial-boundary value conditions such as the Camassa-Holm (CH) equation, Degasperis-Procesi equation, modified CH equation with cubic nonlinearity, Novikov equation with cubic nonlinearity, mCH-Novikov equation, b-family equation with quartic nonlinearity, generalized modified CH equation with quintic nonlinearity, and etc. These results will be useful to further study the peakon solutions and corresponding experimental design of nonlinear dispersive equations.

categories: nlin.PS cs.LG math-ph math.MP nlin.SI physics.flu-dyn


title: Thermalization and Prethermalization in Periodically Kicked Quantum Spin Chains

abstract: We study the dynamics of periodically-kicked many-body systems away from the high-frequency regime, and discuss a family of Floquet systems where the notion of prethermalization can be naturally extended to intermediate and low driving frequencies. We investigate numerically the dynamics of both integrable and nonintegrable systems, and report on the formation of a long-lived prethermal plateau, akin to the high-frequency limit, where the system thermalizes with respect to an effective Hamiltonian captured by the inverse-frequency expansion (IFE). Unlike the high-frequency regime, we find that the relevant heating times are model dependent: we analyze the stability of the prethermal plateau to small perturbations in the drive period, and show that, in a spin chain whose IFE is intractable, the plateau duration is insensitive to the perturbation strength, in contrast to a chain where the IFE admits the resummation of an entire subseries. Infinitesimal perturbations are enough to restore the ergodic properties of the system, and decrease residual finite-size effects. Although the regime where the Floquet system leaves the prethermal plateau and starts heating up to infinite temperature is not captured by the IFE, we provide evidence that the evolved subsystem is described well by a thermal state w.r.t.~the IFE Hamiltonian, with a gradually changing temperature, in accord with the Eigenstate Thermalization Hypothesis.


title: Quotient maps of 2,3-uniform tilings of the plane on the torus

abstract: A 2-uniform tiling is an edge-to-edge tiling by regular polygons having $2$ distinct transitivity classes of vertices. There are 20 distinct 2-uniform tilings (these are of $14$ different types) on the plane, and since the plane is the universal cover of the torus, it is natural to explore maps on the torus that correspond to the 2-uniform tilings. In this article, we discuss that if a map is the quotient of a plane's $2$-uniform lattice then what would be the bounds of the number of vertex orbits. A 3-uniform tiling is an edge-to-edge tiling by regular polygons having $3$ distinct transitivity classes of vertices. There are $61$ distinct $3$-uniform tilings on the plane. In this article, we discuss that if a map is the quotient of a plane's $3$-uniform lattice then what would be the bounds of the number of vertex orbits.


title: Micropillar compression of single crystal tungsten carbide, Part 1: temperature and orientation dependence of deformation behaviour

abstract: Tungsten carbide cobalt hardmetals are commonly used as cutting tools subject to high operation temperature and pressures, where the mechanical performance of the tungsten carbide phase affects the wear and lifetime of the material. In this study, the mechanical behaviour of the isolated tungsten carbide (WC) phase was investigated using single crystal micropillar compression. Micropillars in two crystal orientations, 1-5 ${\mu}$m in diameter, were fabricated using focused ion beam (FIB) machining and subsequently compressed between room temperature and 600 {\deg}C. The activated plastic deformation mechanisms were strongly anisotropic and weakly temperature dependent. The flow stresses of basal-oriented pillars were about three times higher than the prismatic pillars, and pillars of both orientations soften slightly with increasing temperature. The basal pillars tended to deform by either unstable cracking or unstable yield, whereas the prismatic pillars deformed by slip-mediated cracking. However, the active deformation mechanisms were also sensitive to pillar size and shape. Slip trace analysis of the deformed pillars showed that {10-10} prismatic planes were the dominant slip plane in WC. Basal slip was also activated as a secondary slip system at high temperatures.


title: Towards Stratified Space Learning: Linearly Embedded Graphs

abstract: In this paper, we consider the simplest class of stratified spaces -- linearly embedded graphs. We present an algorithm that learns the abstract structure of an embedded graph and models the specific embedding from a point cloud sampled from it. We use tools and inspiration from computational geometry, algebraic topology, and topological data analysis and prove the correctness of the identified abstract structure under assumptions on the embedding. The algorithm is implemented in the Julia package http://github.com/yossibokor/Skyler.jl , which we used for the numerical simulations in this paper.


title: Fingerprints of optical absorption in the perovskite LaInO$_{3}$: Insight from many-body theory and experiment

abstract: We provide a combined theoretical and experimental study of the electronic structure and the optical absorption edge of the orthorhombic perovskite LaInO$_{3}$. Employing density-functional theory and many-body perturbation theory, we predict a direct electronic quasiparticle band gap of about 5 eV and an effective electron (hole) mass of 0.31 (0.48) m$_{0}$. We find the lowest-energy excitation at 0.2 eV below the fundamental gap, reflecting a sizeable electron-hole attraction. Since the transition from the valence band maximum (VBM, $\Gamma$ point) is, however, dipole forbidden the onset is characterized by weak excitations from transitions around it. The first intense excitation appears about 0.32 eV above. Interestingly, this value coincides with an experimental value obtained by ellipsometry (4.80 eV) which is higher than the onset from optical absorption spectroscopy (4.35 eV). The latter discrepancy is attributed to the fact that the weak transitions that define the optical gap are not resolved by the ellipsometry measurement. The absorption edge shows a strong dependency on the light polarization, reflecting the character of the involved valence states. Temperature-dependent measurements show a redshift of the optical gap by about 120 meV by increasing the temperature from 5 to 300 K. Renormalization due to zero-point vibrations is extrapolated from the latter measurement to amount to 150 meV. By adding the excitonic binding energy of 0.2 eV obtained theoretically to the experimental optical absorption onset, we determine the fundamental band gap at room temperature to be 4.55 eV.


title: CAnet: Uplink-aided Downlink Channel Acquisition in FDD Massive MIMO using Deep Learning

abstract: In frequency-division duplexing systems, the downlink channel state information (CSI) acquisition scheme leads to high training and feedback overheads. In this paper, we propose an uplink-aided downlink channel acquisition framework using deep learning to reduce these overheads. Unlike most existing works that focus only on channel estimation or feedback modules, to the best of our knowledge, this is the first study that considers the entire downlink CSI acquisition process, including downlink pilot design, channel estimation, and feedback. First, we propose an adaptive pilot design module by exploiting the correlation in magnitude among bidirectional channels in the angular domain to improve channel estimation. Next, to avoid the bit allocation problem during the feedback module, we concatenate the complex channel and embed the uplink channel magnitude to the channel reconstruction at the base station. Lastly, we combine the above two modules and compare two popular downlink channel acquisition frameworks. The former framework estimates and feeds back the channel at the user equipment subsequently. The user equipment in the latter one directly feeds back the received pilot signals to the base station. Our results reveal that, with the help of uplink, directly feeding back the pilot signals can save approximately 20% of feedback bits, which provides a guideline for future research.


title: Rethinking Interactive Image Segmentation: Feature Space Annotation

abstract: Despite the progress of interactive image segmentation methods, high-quality pixel-level annotation is still time-consuming and laborious - a bottleneck for several deep learning applications. We take a step back to propose interactive and simultaneous segment annotation from multiple images guided by feature space projection. This strategy is in stark contrast to existing interactive segmentation methodologies, which perform annotation in the image domain. We show that feature space annotation achieves competitive results with state-of-the-art methods in foreground segmentation datasets: iCoSeg, DAVIS, and Rooftop. Moreover, in the semantic segmentation context, it achieves 91.5% accuracy in the Cityscapes dataset, being 74.75 times faster than the original annotation procedure. Further, our contribution sheds light on a novel direction for interactive image annotation that can be integrated with existing methodologies. The supplementary material presents video demonstrations. Code available at https://github.com/LIDS-UNICAMP/rethinking-interactive-image-segmentation.


title: Rigorous derivation of dark energy and inflation as geometry effects in Covariant Canonical Gauge Gravity

abstract: The cosmological implications of the Covariant Canonical Gauge Theory of Gravity (CCGG) are investigated. CCGG is a Palatini theory derived from first principles using the canonical transformation formalism in the covariant Hamiltonian formulation. The Einstein-Hilbert theory is thereby extended by a quadratic Riemann-Cartan term in the Lagrangian. Moreover, the requirement of covariant conservation of the stress-energy tensor leads to necessary presence of torsion. In the Friedman universe that promotes the cosmological constant to a time-dependent function, and gives rise to a geometrical correction with the EOS of dark radiation. The resulting cosmology, compatible with the $\Lambda$CDM parameter set, encompasses bounce and bang scenarios with graceful exits into the late dark energy era. Testing those scenarios against low-z observations shows that CCGG is a viable theory.


title: Adaptation of the tapered double cantilever beam test for the measurement of fracture energy and its variations with crack speed

abstract: In this work we present the design of a new test geometry inspired by the Tapered Double Cantilever Beam (TDCB) specimen that is shown to provide an improved characterization of the fracture properties of brittle solids. First, we show that our new design results in an exponential increase of the specimen compliance with crack length, leading to an extremely stable crack growth during the test. We determine an analytical description of this behavior, which provides a simple procedure to extract the fracture energy without depending on finite element calculations. Validation tests are done on polymethylmethacrylate (PMMA) specimens. We use both finite element simulations and our analytical model to interpret the data. We find a very good agreement between the toughness determined by both methods. The stable nature of crack growth in our improved TDCB specimens results in a precise control of the crack speed. This feature is employed to go one step further and characterize the variations of toughness with crack speed. We propose an original optimization procedure for the determination of the material parameters characterizing the kinetic law describing the toughness rate dependency. Overall, the approach proposed together with the newly designed test geometry offer unprecedented possibilities for the full and accurate characterization of the fracture behavior of brittle materials such as rocks, sandstone, mortar etc.


title: "It is like egg": Paul Lorenzen and the collapse of proofs of consistency

abstract: Paul Lorenzen, mathematician and philosopher of the 20th century, mentions October 1947 as the date of a crisis in his mathematical and philosophical investigations. An autograph dated 15 October 1947 documents this crisis. This article proposes a translation and a commentary of it and sketches the circumstances of its writing on the base of his correspondence with Paul Bernays. A letter from Lorenzen to Carl Friedrich Gethmann dated 14 January 1988 carves out the story of this crisis by showing how he soaks up the indications of his correspondents and transmutes them into an absolutely original research. Paul Lorenzen, math{\'e}maticien et philosophe du 20e si{\`e}cle, {\'e}voque le mois d'octobre 1947 comme la date d'une c{\'e}sure dans ses investigations tant math{\'e}matiques que philosophiques. Un {\'e}crit autographe dat{\'e} du 15 octobre 1947 marque cette c{\'e}sure. Cet article en propose une traduction et un commentaire ligne {\`a} ligne, puis esquisse les circonstances de son {\'e}criture sur la base de la correspondance avec Paul Bernays, figure tut{\'e}laire de la logique math{\'e}matique et l'un de ses principaux interlocuteurs {\`a} cette {\'e}poque, {\`a} l'instar de Heinrich Scholz et Oskar Becker. Une lettre de Lorenzen {\`a} Carl Friedrich Gethmann dat{\'e}e du 14 janvier 1988 {\'e}toffe le r{\'e}cit de cette c{\'e}sure en montrant comment il s'impr{\`e}gne des indications que lui fournissent ses correspondants pour les transmuter en une recherche absolument originale.


title: Homogenization of the Stokes system in a non-periodically perforated domain

abstract: In our recent work [8], we have studied the homogenization of the Poisson equation in a class of non periodically perforated domains. In this paper, we examine the case of the Stokes system. We consider a porous medium in which the characteristic distance between two holes, denoted by $\epsilon$, is proportional to the characteristic size of the holes. It is well known (see [1],[17] and [19]) that, when the holes are periodically distributed in space, the velocity converges to a limit given by the Darcy's law when the size of the holes tends to zero. We generalize these results to the setting of [8]. The non-periodic domains are defined as a local perturbation of a periodic distribution of holes. We obtain classical results of the homogenization theory in perforated domains (existence of correctors and regularity estimates uniform in $\epsilon$) and we prove $H^2$--convergence estimates for particular force fields.


title: Interpretable discovery of new semiconductors with machine learning

abstract: Machine learning models of materials$^{1-5}$ accelerate discovery compared to ab initio methods: deep learning models now reproduce density functional theory (DFT)-calculated results at one hundred thousandths of the cost of DFT$^{6}$. To provide guidance in experimental materials synthesis, these need to be coupled with an accurate yet effective search algorithm and training data consistent with experimental observations. Here we report an evolutionary algorithm powered search which uses machine-learned surrogate models trained on high-throughput hybrid functional DFT data benchmarked against experimental bandgaps: Deep Adaptive Regressive Weighted Intelligent Network (DARWIN). The strategy enables efficient search over the materials space of ~10$^8$ ternaries and 10$^{11}$ quaternaries$^{7}$ for candidates with target properties. It provides interpretable design rules, such as our finding that the difference in the electronegativity between the halide and B-site cation being a strong predictor of ternary structural stability. As an example, when we seek UV emission, DARWIN predicts K$_2$CuX$_3$ (X = Cl, Br) as a promising materials family, based on its electronegativity difference. We synthesized and found these materials to be stable, direct bandgap UV emitters. The approach also allows knowledge distillation for use by humans.


title: Modeling Heavy-Gas Dispersion in Air with Two-Layer Shallow Water Equations

abstract: Computation of gas dispersal in urban places or hilly grounds requires a large amount of computer time when addressed with conventional multidimensional models. Those are usually based on two-phase flow or Navier-Stokes equations. Different classes of simplified models exist. Among them, two-layer shallow water models are interesting to address large-scale dispersion. Indeed, compared to conventional multidimensional approaches, 2D simulations are carried out to mimic 3D effects. The computational gain in CPU time is consequently expected to be tremendous. However, such models involve at least three fundamental difficulties. The first one is related to the lack of hyperbolicity of most existing formulations, yielding serious consequences regarding wave propagation. The second is related to the non-conservative terms in the momentum equations. Those terms account for interactions between fluid layers. Recently, these two difficulties have been addressed in Chiapolino and Saurel (2018) and an unconditional hyperbolic model has been proposed along with a Harten-Lax-van Leer (HLL) type Riemann solver dealing with the non-conservative terms. In the same reference, numerical experiments showed robustness and accuracy of the formulation. In the present paper, a third difficulty is addressed. It consists of the determination of appropriate drag effect formulation. Such effects also account for interactions between fluid layers and become of particular importance when dealing with heavy-gas dispersion. With this aim, the model is compared to laboratory experiments in the context of heavy gas dispersal in quiescent air. It is shown that the model accurately reproduces experimental results thanks to an appropriate drag force correlation. This function expresses drag effects between the heavy and light gas layers. It is determined thanks to various experimental configurations of dam-break test problems.


title: Structure and Morphology of Crystalline Metal Nanoparticles: Polyhedral Cubic Particles

abstract: We examine nanoparticles (NPs) forming polyhedral sections of the ideal cubic lattice, simple (sc), body centered (bcc), and face centered (fcc) cubic, which are confined by facets characterized by densest and second densest {h k l} monolayers of the lattice. Together with the constraint that the NPs exhibit the same point symmetry as the ideal cubic lattice, i.e. Oh, different types of generic NPs serve for the definition of general polyhedral cubic NPs. Their structural properties, such as shape, size, and surface facets, are discussed in analytical and numerical detail with visualization of characteristic examples. The geometric relationships of the model particles expressed by corresponding formulas and numerical tables can be used to estimate shapes and sizes of real compact metal nanoparticles observed by experiment.


title: Using uncertainty estimation to reduce false positives in liver lesion detection

abstract: Despite the successes of deep learning techniques at detecting objects in medical images, false positive detections occur which may hinder an accurate diagnosis. We propose a technique to reduce false positive detections made by a neural network using an SVM classifier trained with features derived from the uncertainty map of the neural network prediction. We demonstrate the effectiveness of this method for the detection of liver lesions on a dataset of abdominal MR images. We find that the use of a dropout rate of 0.5 produces the least number of false positives in the neural network predictions and the trained classifier filters out approximately 90% of these false positives detections in the test-set.


title: Ancillary services in Great Britain during the COVID-19 lockdown: a glimpse of the carbon-free future

abstract: The COVID-19 pandemic led to partial or total lockdowns in several countries during the first half of 2020, which in turn caused a depressed electricity demand. In Great Britain (GB), this low demand combined with large renewable output at times, created conditions that were not expected until renewable capacity increases to meet emissions targets in coming years. The GB system experienced periods of very high instantaneous penetration of non-synchronous renewables, compromising system stability due to the lack of inertia in the grid. In this paper, a detailed analysis of the consequences of the lockdown on the GB electricity system is provided, focusing on the ancillary services procured to guarantee stability. Ancillary-services costs increased by {\pounds}200m in the months of May to July 2020 compared to the same period in 2019 (a threefold increase), highlighting the importance of ancillary services in low-carbon systems. Furthermore, a frequency-secured scheduling model is used in the present paper to showcase the future trends that GB is expected to experience, as penetration of renewables increases on the road to net-zero emissions by 2050. Several sensitivities are considered, demonstrating that the share of total operating costs represented by ancillary services could reach 35%.


title: Dynamic Spectrum Access using Stochastic Multi-User Bandits

abstract: A stochastic multi-user multi-armed bandit framework is used to develop algorithms for uncoordinated spectrum access. In contrast to prior work, it is assumed that rewards can be non-zero even under collisions, thus allowing for the number of users to be greater than the number of channels. The proposed algorithm consists of an estimation phase and an allocation phase. It is shown that if every user adopts the algorithm, the system wide regret is order-optimal of order $O(\log T)$ over a time-horizon of duration $T$. The regret guarantees hold for both the cases where the number of users is greater than or less than the number of channels. The algorithm is extended to the dynamic case where the number of users in the system evolves over time, and is shown to lead to sub-linear regret.


title: Virgo Filaments I: Processing of gas in cosmological filaments around the Virgo cluster

abstract: Galaxies have different morphology, gas content and star formation rate (SFR) in dense environments like galaxy clusters. The impact of environmental density extends to several virial radii, and galaxies are pre-processed in filaments and groups, before falling into the cluster. Our goal is to quantify this pre-processing, in terms of gas content and SFR, as a function of density in cosmic filaments. We have observed the two first CO transitions in 163 galaxies with the IRAM-30m telescope, and added 82 more measurements from the literature, for a sample of 245 galaxies in the filaments around Virgo. We gathered HI-21cm measurements from the literature, and observed 69 galaxies with the Nan\c{c}ay telescope, to complete our sample. We compare our filament galaxies with comparable samples from the Virgo cluster and with the isolated galaxies of the AMIGA sample. We find a progression from field galaxies to filament and cluster ones for decreasing SFR, increasing fraction of galaxies in the quenching phase, increasing proportion of early-type galaxies and decreasing gas content. Galaxies in the quenching phase, defined as having SFR below 1/3 of the main sequence rate, are only between 0-20\% in the isolated sample, while they are 20-60\% in the filaments and 30-80\% in the Virgo cluster. Processes that lead to star formation quenching are already at play in filaments. They depend mostly on the local galaxy density, while the distance to the filament spine is a secondary parameter. While the HI to stellar mass ratio decreases with local density by ~1 dex in the filaments, and ~2 dex in Virgo with respect to the field, the decrease is much less for the H$_2$ to stellar mass ratio. As the environmental density increases, the gas depletion time decreases, since the gas content decreases faster than the SFR. This suggests that gas depletion significantly precedes star formation quenching.


title: New Bias Calibration for Robust Estimation in Small Areas

abstract: Using sample surveys as a cost effective tool to provide estimates for characteristics of interest at population and sub-populations (area/domain) level has a long tradition in "small area estimation". However, the existence of outliers in the sample data can significantly affect the estimation for areas in which they occur, especially where the domain-sample size is small. Based on existing robust estimators for small area estimation we propose two novel approaches for bias calibration. A series of simulations shows that our methods lead to more efficient estimators in comparison with other existing bias-calibration methods. As a real data example we apply our estimators to obtain \textit{Gini} coefficients in labour market areas of the Tuscany region of Italy, where our sources of information are the EU-SILC survey and the Italian census. This analysis shows that the new methods reveal a different picture than existing methods. We extend our ideas to predictions for non-sampled areas.


title: Spatially-resolved decoherence of donor spins in silicon strained by a metallic electrode

abstract: Electron spins are amongst the most coherent solid-state systems known, however, to be used in devices for quantum sensing and information processing applications, they must be typically placed near interfaces. Understanding and mitigating the impacts of such interfaces on the coherence and spectral properties of electron spins is critical to realize such applications, but is also challenging: inferring such data from single-spin studies requires many measurements to obtain meaningful results, while ensemble measurements typically give averaged results that hide critical information. Here, we report a comprehensive study of the coherence of near-surface bismuth donor spins in 28-silicon at millikelvin temperatures. In particular, we use strain-induced frequency shifts caused by a metallic electrode to make spatial maps of spin coherence as a function of depth and position relative to the electrode. By measuring magnetic-field-insensitive clock transitions we separate magnetic noise caused by surface spins from charge noise. Our results include quantitative models of the strain-split spin resonance spectra and extraction of paramagnetic impurity concentrations at the silicon surface. The interplay of these decoherence mechanisms for such near-surface electron spins is critical for their application in quantum technologies, while the combination of the strain splitting and clock transition extends the coherence lifetimes by up to two orders of magnitude, reaching up to 300 ms at a mean depth of only 100nm. The technique we introduce here to spatially map coherence in near-surface ensembles is directly applicable to other spin systems of active interest, such as defects in diamond, silicon carbide, and rare earth ions in optical crystals.


title: On the set of robust sustainable thresholds for uncertain control systems

abstract: In natural resource management, or more generally in the study of sustainability issues, often the objective is to maintain the state of a given system within a desirable configuration, typically established in terms of standards or thresholds. For instance, in fisheries management, the procedure of designing policies may include keeping the spawning stock biomass over a critical threshold and also ensuring minimal catches. Given a controlled dynamical system in discrete-time, representing the evolution of some natural resources under the action of controls and uncertainties, and an initial endowment of the resources, the aim of this paper is to characterize the set of robust sustainable thresholds, that is, the thresholds for which there exists some control path, along with its corresponding state trajectory, satisfying for all possible uncertainty scenarios, prescribed mixed constraints parametrized by such thresholds. This set provides useful information to users and decision-makers, illustrating the trade-offs between constraints and it is strongly related to the robust viability, one of the key concepts in viability theory, discipline that study the consistency between a controlled dynamical system and given constraints. Specifically, we are concerned with characterizing the weak and strong Pareto fronts of the set of robust sustainable thresholds, providing a practical method for computing such objects based on optimal control theory and a level-set approach. A numerical example, relying on renewable resource management, is shown to demonstrate the proposed method.


title: Detecting and quantifying palaeoseasonality in stalagmites using geochemical and modelling approaches

abstract: Stalagmites are an extraordinarily powerful resource for the reconstruction of climatological palaeoseasonality. Here, we provide a comprehensive review of different types of seasonality preserved by stalagmites and methods for extracting this information. A new drip classification scheme is introduced, which facilitates the identification of stalagmites fed by seasonally responsive drips and which highlights the wide variability in drip types feeding stalagmites. This hydrological variability, combined with seasonality in Earth atmospheric processes, meteoric precipitation, biological processes within the soil, and cave atmosphere composition means that every stalagmite retains a different and distinct (but correct) record of environmental conditions. Replication of a record is extremely useful but should not be expected unless comparing stalagmites affected by the same processes in the same proportion. A short overview of common microanalytical techniques is presented, and suggested best practice discussed. In addition to geochemical methods, a new modelling technique for extracting meteoric precipitation and temperature palaeoseasonality from stalagmite d18O data is discussed and tested with both synthetic and real-world datasets. Finally, world maps of temperature, meteoric precipitation amount, and meteoric precipitation oxygen isotope ratio seasonality are presented and discussed, with an aim of helping to identify regions most sensitive to shifts in seasonality.


title: Combining magneto-hydrostatic constraints with Stokes profile inversions. II. Application to Hinode/SP observations

abstract: Inversion techniques applied to the radiative transfer equation for polarized light are capable of inferring the physical parameters in the solar atmosphere (temperature $T$, magnetic field ${\bf B}$, and line-of-sight velocity $v_{\rm los}$) from observations of the Stokes vector (i.e., spectropolarimetric observations) in spectral lines. Inferences are usually performed in the $(x,y,\tau_c)$ domain, where $\tau_c$ refers to the optical-depth scale. Generally, their determination in the $(x,y,z)$ volume is not possible due to the lack of a reliable estimation of the gas pressure, particularly in regions of the solar surface harboring strong magnetic fields. We aim to develop a new inversion code capable of reliably inferring the physical parameters in the $(x,y,z)$ domain. We combine, in a self-consistent way, an inverse solver for the radiative transfer equation (Firtez-DZ) with a solver for the magneto-hydrostatic (MHS) equilibrium, which derives realistic values of the gas pressure by taking the magnetic pressure and tension into account. We test the correct behavior of the newly developed code with spectropolarimetric observations of two sunspots recorded with the spectropolarimeter (SP) instrument on board the Hinode spacecraft, and we show how the physical parameters are inferred in the $(x,y,z)$ domain, with the Wilson depression of the sunspots arising as a natural consequence of the force balance. In particular, our approach significantly improves upon previous determinations that were based on semiempirical models. Our results open the door for the possibility of calculating reliable electric currents in three dimensions, ${\bf j}(x,y,z)$, in the solar photosphere. Further consistency checks would include a comparison with other methods that have recently been proposed and which achieve similar goals.


title: Symbolic Loop Compilation for Tightly Coupled Processor Arrays

abstract: Loop compilation for Tightly Coupled Processor Arrays (TCPAs), a class of massively parallel loop accelerators, entails solving NP-hard problems, yet depends on the loop bounds and number of available processing elements (PEs), parameters known only at runtime because of dynamic resource management and input sizes. Therefore, this article proposes a two-phase approach called symbolic loop compilation: At compile time, the necessary NP-complete problems are solved and the solutions compiled into a space-efficient symbolic configuration. At runtime, a concrete configuration is generated from the symbolic configuration according to the parameters values. We show that the latter phase, called instantiation, runs in polynomial time with its most complex step, program instantiation, not depending on the number of PEs. As validation, we performed symbolic loop compilation on real-world loops and measured time and space requirements. Our experiments confirm that a symbolic configuration is space-efficient and suited for systems with little memory -- often, a symbolic configuration is smaller than a single concrete configuration -- and that program instantiation scales well with the number of PEs -- for example, when instantiating a symbolic configuration of a matrix-matrix multiplication, the execution time is similar for $4\times 4$ and $32\times 32$ PEs.


title: Numerical radius in Hilbert $C^*$-modules

abstract: Utilizing the linking algebra of a Hilbert $C^*$-module $\big(\mathscr{V}, {\|\!\cdot\!\|}\big)$, we introduce $\Omega(x)$ as a definition of numerical radius for an element $x\in\mathscr{V}$ and then show that $\Omega(\cdot)$ is a norm on $\mathscr{V}$ such that $\frac{1}{2}{\|x\|} \leq \Omega(x) \leq {\|x\|}$. In addition, we obtain an equivalent condition for $\Omega(x) = \frac{1}{2}{\|x\|}$. Moreover, we present a refinement of the triangle inequality for the norm $\Omega(\cdot)$. Some other related results are also discussed.


title: On the number of isolate dominating sets of certain graphs

abstract: Let $G=(V,E)$ be a simple graph. A dominating set of $G$ is a subset $S\subseteq V$ such that every vertex not in $S$ is adjacent to at least one vertex in $S$. The cardinality of a smallest dominating set of $G$, denoted by $\gamma(G)$, is the domination number of $G$. A dominating set $S$ is an isolate dominating set of $G$, if the induced subgraph $G[S]$ has at least one isolated vertex. The isolate domination number, $\gamma_0(G)$, is the minimum cardinality of an isolate dominating set of $G$. In this paper, we count the number of isolate dominating sets of some specific graphs.


title: On the distribution of prime divisors in Krull monoid algebras

abstract: In the present work, we prove that every class of the divisor class group of a Krull monoid algebra contains infinitely many prime divisors. Several attempts to this result have been made in the literature so far, unfortunately with open gaps. We present a complete proof of this fact.


title: Low thermal conductivity of the superfast rotator (499998) 2011 PT

abstract: Context: Asteroids with a diameter of up to a few dozen meters may spin very fast and complete an entire rotation within a few minutes. These small and fast-rotating bodies are thought to be monolithic objects because the gravitational force due to their small size is not strong enough to counteract the strong centripetal force caused by the fast rotation. Additionally, it is not clear whether the fast spin prevents dust and small particles (regolith) from being kept on their surface. Aims: We develop a model for constraining the thermal conductivity of the surface of the small, fast-rotating near-Earth asteroids. This model may suggest whether regolith is likely present on these objects. Methods: Our approach is based on the comparison of the measured Yarkovsky drift and a predicted value using a theoretical model that depends on the orbital, physical and thermal parameters of the object. The necessary parameters are either deduced from statistical distribution derived for near-Earth asteroids population or determined from observations with associated uncertainty. With this information, we performed Monte Carlo simulations and produced a probability density distribution for the thermal conductivity. Results: Applying our model to the superfast rotator asteroid (499998) 2011 PT, we find that the measured Yarkovsky drift can only be achieved when the thermal conductivity $K$ of the surface is low. The resulting probability density function for the conductivity is bimodal, with two most likely values being around 0.0001 and 0.005 W m$^{-1}$ K$^{-1}$. Based on this, we find that the probability that $K$ is lower than 0.1 W m$^{-1}$ K$^{-1}$ is at least 95\%. This low thermal conductivity might indicate that the surface of 2011 PT is covered with a thermal insulating layer, composed of a regolith-like material similar to lunar dust.