title: Identifying the Development and Application of Artificial Intelligence in Scientific Text

abstract: We describe a strategy for identifying the universe of research publications relevant to the application and development of artificial intelligence. The approach leverages the arXiv corpus of scientific preprints, in which authors choose subject tags for their papers from a set defined by editors. We compose a functional definition of AI relevance by learning these subjects from paper metadata, and then inferring the arXiv-subject labels of papers in larger corpora: Clarivate Web of Science, Digital Science Dimensions, and Microsoft Academic Graph. This yields predictive classification $F_1$ scores between .75 and .86 for Natural Language Processing (cs.CL), Computer Vision (cs.CV), and Robotics (cs.RO). For a single model that learns these and four other AI-relevant subjects (cs.AI, cs.LG, stat.ML, and cs.MA), we see precision of .83 and recall of .85. We evaluate the out-of-domain performance of our classifiers against other sources of topic information and predictions from alternative methods. We find that a supervised solution can generalize to identify publications that belong to the high-level fields of study represented on arXiv. This offers a method for identifying AI-relevant publications that updates at the pace of research output, without reliance on subject-matter experts for query development or labeling.


title: Quantifying Noninvertibility in Discrete Dynamical Systems

abstract: Given a finite set $X$ and a function $f:X\to X$, we define the degree of noninvertibility of $f$ to be $\displaystyle\text{deg}(f)=\frac{1}{|X|}\sum_{x\in X}|f^{-1}(f(x))|$. This is a natural measure of how far the function $f$ is from being bijective. We compute the degrees of noninvertibility of some specific discrete dynamical systems, including the Carolina solitaire map, iterates of the bubble sort map acting on permutations, bubble sort acting on multiset permutations, and a map that we call "nibble sort." We also obtain estimates for the degrees of noninvertibility of West's stack-sorting map and the Bulgarian solitaire map. We then turn our attention to arbitrary functions and their iterates. In order to compare the degree of noninvertibility of an arbitrary function $f:X\to X$ with that of its iterate $f^k$, we prove that \[\max_{\substack{f:X\to X\\ |X|=n}}\frac{\text{deg}(f^k)}{\text{deg}(f)^\gamma}=\Theta(n^{1-1/2^{k-1}})\] for every real number $\gamma\geq 2-1/2^{k-1}$. We end with several conjectures and open problems.


title: Dark Matter bound-state formation at higher order: a non-equilibrium quantum field theory approach

abstract: The formation of meta-stable dark matter bound states in coannihilating scenarios could efficiently occur through the scattering with a variety of Standard Model bath particles, where light bosons during the electroweak cross over or even massless photons and gluons are exchanged in the t-channel. The amplitudes for those higher-order processes, however, are divergent in the collinear direction of the in- and out-going bath particles if the mediator is massless. To address the issue of collinear divergences, we derive the bound-state formation collision term in the framework of non-equilibrium quantum field theory. The main result is an expression for a more general cross section, which allows to compute higher-order bound-state formation processes inside the primordial plasma background in a comprehensive manner. Based on this result, we show that next-to-leading order contributions, including the bath-particle scattering, are i) collinear finite and ii) generically dominate over the on-shell emission for temperatures larger than the absolute value of the binding energy. Based on a simplified model, we demonstrate that the impact of these new effects on the thermal relic abundance is significant enough to make it worthwhile to study more realistic coannihilation scenarios.


title: Boundaries of the Amplituhedron with amplituhedronBoundaries

abstract: Positive geometries provide a modern approach for computing scattering amplitudes in a variety of physical models. In order to facilitate the exploration of these new geometric methods, we introduce a Mathematica package called ``amplituhedronBoundaries'' for calculating the boundary structures of three positive geometries: the amplituhedron $\mathcal{A}_{n,k}^{(m)}$, the momentum amplituhedron $\mathcal{M}_{n,k}^{(m)}$ and the hypersimplex $\Delta_{k,n}$. The first two geometries are relevant for scattering amplitudes in planar $\mathcal{N}=4$ SYM, while the last one is a well-studied polytope appearing in many contexts in mathematics, and is closely related to $\mathcal{M}_{n,k}^{(2)}$. The package includes an array of useful tools for the study of these positive geometries, including their boundary stratifications, drawing their boundary posets, and additional tools for manipulating combinatorial structures useful for positive Grassmannians.


title: Fair Prediction with Endogenous Behavior

abstract: There is increasing regulatory interest in whether machine learning algorithms deployed in consequential domains (e.g. in criminal justice) treat different demographic groups "fairly." However, there are several proposed notions of fairness, typically mutually incompatible. Using criminal justice as an example, we study a model in which society chooses an incarceration rule. Agents of different demographic groups differ in their outside options (e.g. opportunity for legal employment) and decide whether to commit crimes. We show that equalizing type I and type II errors across groups is consistent with the goal of minimizing the overall crime rate; other popular notions of fairness are not.


title: Geometrically Nonlinear Response of a Fractional-Order Nonlocal Model of Elasticity

abstract: This study presents the analytical and finite element formulation of a geometrically nonlinear and fractional-order nonlocal model of an Euler-Bernoulli beam. The finite nonlocal strains in the Euler-Bernoulli beam are obtained from a frame-invariant and dimensionally consistent fractional-order (nonlocal) continuum formulation. The finite fractional strain theory provides a positive definite formulation that results in a mathematically well-posed formulation which is consistent across loading and boundary conditions. The governing equations and the corresponding boundary conditions of the geometrically nonlinear and nonlocal Euler-Bernoulli beam are obtained using variational principles. Further, a nonlinear finite element model for the fractional-order system is developed in order to achieve the numerical solution of the integro-differential nonlinear governing equations. Following a thorough validation with benchmark problems, the fractional finite element model (f-FEM) is used to study the geometrically nonlinear response of a nonlocal beam subject to various loading and boundary conditions. Although presented in the context of a 1D beam, this nonlinear f-FEM formulation can be extended to higher dimensional fractional-order boundary value problems.


title: Co-adjoint orbits and time-optimal problems for two-step free-nilpotent Lie groups

abstract: We describe co-adjoint orbits and Casimir functions for two-step free-nilpotent Lie algebras. The symplectic foliation consists of affine subspaces of the Lie coalgebra of different dimensions. Further, we consider left-invariant time-optimal problems on two-step Carnot groups, for which the set of admissible velocities is a strictly convex compactum in the first layer of the Lie algebra containing the origin in its interior. We describe integrals for the vertical subsystem of the Hamiltonian system of Pontryagin maximum principle. Further, we describe constancy and periodicity of solutions to this subsystem and controls, and characterize its flow, for two-dimensional co-adjoint orbits.


title: Asymptotic Behavior of Spanning Forests and Connected Spanning Subgraphs on Two-Dimensional Lattices

abstract: We calculate exponential growth constants $\phi$ and $\sigma$ describing the asymptotic behavior of spanning forests and connected spanning subgraphs on strip graphs, with arbitrarily great length, of several two-dimensional lattices, including square, triangular, honeycomb, and certain heteropolygonal Archimedean lattices. By studying the limiting values as the strip widths get large, we infer lower and upper bounds on these exponential growth constants for the respective infinite lattices. Since our lower and upper bounds are quite close to each other, we can infer very accurate approximate values for these exponential growth constants, with fractional uncertainties ranging from $O(10^{-4})$ to $O(10^{-2})$. We show that $\phi$ and $\sigma$, are monotonically increasing functions of vertex degree for these lattices.


title: On the complexity of finding tensor ranks

abstract: The purpose of this note is to give a linear algebra algorithm to find out if a rank of a given tensor over a field $\F$ is at most $k$ over the algebraic closure of $\F$, where $k$ is a given positive integer. We estimate the arithmetic complexity of our algorithm.


title: Weighted Additive Spanners

abstract: A \emph{spanner} of a graph $G$ is a subgraph $H$ that approximately preserves shortest path distances in $G$. Spanners are commonly applied to compress computation on metric spaces corresponding to weighted input graphs. Classic spanner constructions can seamlessly handle edge weights, so long as error is measured \emph{multiplicatively}. In this work, we investigate whether one can similarly extend constructions of spanners with purely \emph{additive} error to weighted graphs. These extensions are not immediate, due to a key lemma about the size of shortest path neighborhoods that fails for weighted graphs. Despite this, we recover a suitable amortized version, which lets us prove direct extensions of classic $+2$ and $+4$ unweighted spanners (both all-pairs and pairwise) to $+2W$ and $+4W$ weighted spanners, where $W$ is the maximum edge weight. Specifically, we show that a weighted graph $G$ contains all-pairs (pairwise) $+2W$ and $+4W$ weighted spanners of size $O(n^{3/2})$ and $\widetilde{O}(n^{7/5})$ ($O(np^{1/3})$ and $O(np^{2/7})$) respectively. For a technical reason, the $+6$ unweighted spanner becomes a $+8W$ weighted spanner; closing this error gap is an interesting remaining open problem. That is, we show that $G$ contains all-pairs (pairwise) $+8W$ weighted spanners of size $O(n^{4/3})$ ($O(np^{1/4})$).


title: Cover Combinatorial Filters and their Minimization Problem (Extended Version)

abstract: Recent research has examined algorithms to minimize robots' resource footprints. The class of combinatorial filters (discrete variants of widely-used probabilistic estimators) has been studied and methods for reducing their space requirements introduced. This paper extends existing combinatorial filters by introducing a natural generalization that we dub cover combinatorial filters. In addressing the new -- but still NP-complete -- problem of minimization of cover filters, this paper shows that multiple concepts previously believed to be true about combinatorial filters (and actually conjectured, claimed, or assumed to be) are in fact false. For instance, minimization does not induce an equivalence relation. We give an exact algorithm for the cover filter minimization problem. Unlike prior work (based on graph coloring) we consider a type of clique-cover problem, involving a new conditional constraint, from which we can find more general relations. In addition to solving the more general problem, the algorithm also corrects flaws present in all prior filter reduction methods. In employing SAT, the algorithm provides a promising basis for future practical development.


title: Forward-backward algorithms with different inertial terms for structured non-convex minimization problems

abstract: We investigate two inertial forward-backward algorithms in connection with the minimization of the sum of a non-smooth and possibly non-convex and a non-convex differentiable function. The algorithms are formulated in the spirit of the famous FISTA method, however the setting is non-convex and we allow different inertial terms. Moreover, the inertial parameters in our algorithms can take negative values too. We also treat the case when the non-smooth function is convex and we show that in this case a better step size can be allowed. We prove some abstract convergence results which applied to our numerical schemes allow us to show that the generated sequences converge to a critical point of the objective function, provided a regularization of the objective function satisfies the Kurdyka-Lojasiewicz property. Further, we obtain a general result that applied to our numerical schemes ensures convergence rates for the generated sequences and for the objective function values formulated in terms of the KL exponent of a regularization of the objective function. Finally, we apply our results to image restoration.


title: Enabling Low-Power OFDM for IoT by Exploiting Asymmetric Clock Rates

abstract: The conventional high-speed Wi-Fi has recently become a contender for low-power Internet-of-Things (IoT) communications. OFDM continues its adoption in the new IoT Wi-Fi standard due to its spectrum efficiency that can support the demand of massive IoT connectivity. While the IoT Wi-Fi standard offers many new features to improve power and spectrum efficiency, the basic physical layer (PHY) structure of transceiver design still conforms to its conventional design rationale where access points (AP) and clients employ the same OFDM PHY. In this paper, we argue that current Wi-Fi PHY design does not take full advantage of the inherent asymmetry between AP and IoT. To fill the gap, we propose an asymmetric design where IoT devices transmit uplink packets using the lowest power while pushing all the decoding burdens to the AP side. Such a design utilizes the sufficient power and computational resources at AP to trade for the transmission (TX) power of IoT devices. The core technique enabling this asymmetric design is that the AP takes full power of its high clock rate to boost the decoding ability. We provide an implementation of our design and show that it can reduce up to 88% of the IoT's TX power when the AP sets $8\times$ clock rate.


title: Using anisotropic 3D Minkowski functionals for trabecular bone characterization and biomechanical strength prediction in proximal femur specimens

abstract: The ability of Anisotropic Minkowski Functionals (AMFs) to capture local anisotropy while evaluating topological properties of the underlying gray-level structures has been previously demonstrated. We evaluate the ability of this approach to characterize local structure properties of trabecular bone micro-architecture in ex vivo proximal femur specimens, as visualized on multi-detector CT, for purposes of biomechanical bone strength prediction. To this end, volumetric AMFs were computed locally for each voxel of volumes of interest (VOI) extracted from the femoral head of 146 specimens. The local anisotropy captured by such AMFs was quantified using a fractional anisotropy measure; the magnitude and direction of anisotropy at every pixel was stored in histograms that served as a feature vectors that characterized the VOIs. A linear multi-regression analysis algorithm was used to predict the failure load (FL) from the feature sets; the predicted FL was compared to the true FL determined through biomechanical testing. The prediction performance was measured by the root mean square error (RMSE) for each feature set. The best prediction performance was obtained from the fractional anisotropy histogram of AMF Euler Characteristic (RMSE = 1.01 +- 0.13), which was significantly better than MDCT-derived mean BMD (RMSE = 1.12 +- 0.16, p<0.05). We conclude that such anisotropic Minkowski Functionals can capture valuable information regarding regional trabecular bone quality and contribute to improved bone strength prediction, which is important for improving the clinical assessment of osteoporotic fracture risk.


title: Neutron scattering off Weyl semimetals

abstract: We present how to detect type-$1$ Weyl nodes in a material by inelastic neutron scattering. Such an experiment first of all allows one to determine the dispersion of the Weyl fermions. We extend the reasoning to produce a quantitative test of the Weyl equation taking into account realistic anisotropic properties. These anisotropies are mostly contained in the form of the emergent magnetic moment of the excitations, which determines how they couple to the neutron. Although there are many material parameters, we find several quantitative predictions that are universal and demonstrate that the excitations are described by solutions to the Weyl equation. The realistic, anisotropic coupling between electrons and neutrons implies that even fully unpolarized neutrons can reveal the spin-momentum locking of the Weyl fermions because the neutrons will couple to some components of the Weyl fermion pseudospin more strongly. On the other hand, in an experiment with polarized neutrons, the scattered neutron beam remains fully polarized in a direction that varies as a function of momentum transfer (within the range of validity of the Weyl equation). This allows measurement of the chirality of Weyl fermions for inversion symmetric nodes. Furthermore, we estimate that the scattering rate may be large enough for such experiments to be practical; in particular, the magnetic moment may be larger than the ordinary Bohr magneton, compensating for a small density of states.


title: H\"older conditions and $\tau$-spikes for analytic Lipschitz functions

abstract: Let $U$ be an open subset of $\mathbb{C}$ with boundary point $x_0$ and let $A_{\alpha}(U)$ be the space of functions analytic on $U$ that belong to lip$\alpha(U)$, the "little Lipschitz class". We consider the condition $S= \displaystyle \sum_{n=1}^{\infty}2^{(t+\lambda+1)n}M_*^{1+\alpha}(A_n \setminus U)< \infty,$ where $t$ is a non-negative integer, $0<\lambda<1$, $M_*^{1+\alpha}$ is the lower $1+\alpha$ dimensional Hausdorff content, and $A_n = \{z: 2^{-n-1}<|z-x_0|<2^{-n}\}$. This is similar to a necessary and sufficient condition for bounded point derivations on $A_{\alpha}(U)$ at $x_0$. We show that $S= \infty$ implies that $x_0$ is a $(t+\lambda)$-spike for $A_{\alpha}(U)$ and that if $S<\infty$ and $U$ satisfies a cone condition, then the $t$-th derivatives of functions in $A_{\alpha}(U)$ satisfy a H\"older condition at $x_0$ for a non-tangential approach.


title: Genuine NMSSM deviations in the 125 GeV Higgs boson decays

abstract: In Supersymmetry deviations from the SM signal strengths of the 125 GeV Higgs boson can occur, because of possible SUSY contributions in diagrams with loops, which leads to different deviations in signal strengths for processes with and without loop diagrams. In the Next-to Minimal Supersymmetric Standard Model (NMSSM) additional deviations may occur, because of the mixing with the additional singlet-like Higgs boson and/or additional decays into pairs of light particles, like neutralinos, pseudo-scalar Higgs bosons or singlet-like Higgs bosons. In this paper we study these "genuine" NMSSM deviations in detail to check not only their possible size, but also look for correlations or anti-correlations with respect to other channels in the hope to find specific patterns in the deviations as for the SUSY contributions, which occur predominantly in processes including loops in the diagrams. The whole NMSSM parameter space is sampled in a deterministic way. The novel scanning method is discussed in detail in the Appendix. We found three different regions with "genuine" NMSSM deviations, which are largely independent of the production mode. What was surprising: some regions show negativ correlations between final states with fermions and bosons meaning if the signal strengths for fermions decrease the signal strengths for bosonic final states increase, while other regions show positive correlations. The sign of the correlations is a strong function of the mass difference between the observed 125 GeV Higgs boson and the Higgs singlet, so looking for correlations could give useful hints about the existence and mass of the singlet Higgs boson. The features of these "genuine" NMSSM effects have been investigated in representative benchmark points for each of the three regions where a single effect is dominant. These benchmark points have been detailed in the Appendix.


title: Lugares geom\'etricos asociados a dos puntos del plano

abstract: Both the ellipse and the hyperbola are geometric places that can be defined by establishing a relationship between points $P$ of the plane and two fixed points $A$ and $B$ (which are its foci $F'=A$ and $F=B$). Given two points $A$ and $B$ of the plan (which we no longer call the foci $F$ and $F'$), we are going to present three geometric places associated with $A$ and $B$ other than the ellipse and the hyperbola.


title: Surrogate-free machine learning-based organ dose reconstruction for pediatric abdominal radiotherapy

abstract: To study radiotherapy-related adverse effects, detailed dose information (3D distribution) is needed for accurate dose-effect modeling. For childhood cancer survivors who underwent radiotherapy in the pre-CT era, only 2D radiographs were acquired, thus 3D dose distributions must be reconstructed from limited information. State-of-the-art methods achieve this by using 3D surrogate anatomies. These can lack personalization and lead to coarse reconstructions. We present and validate a surrogate-free dose reconstruction method based on Machine Learning (ML). Abdominal planning CTs ($n$=142) of recently-treated childhood cancer patients were gathered, their organs at risk were segmented, and 300 artificial Wilms' tumor plans were sampled automatically. Each artificial plan was automatically emulated on the 142 CTs, resulting in 42,600 3D dose distributions from which dose-volume metrics were derived. Anatomical features were extracted from digitally reconstructed radiographs simulated from the CTs to resemble historical radiographs. Further, patient and radiotherapy plan features typically available from historical treatment records were collected. An evolutionary ML algorithm was then used to link features to dose-volume metrics. Besides 5-fold cross-validation, a further evaluation was done on an independent dataset of five CTs each associated with two clinical plans. Cross-validation resulted in Mean Absolute Errors (MAEs) $\leq$0.6 Gy for organs completely inside or outside the field. For organs positioned at the edge of the field, MAEs $\leq$1.7 Gy for D$_{mean}$, $\leq$2.9 Gy for D$_{2cc}$, and $\leq$13% for V$_{5Gy}$ and V$_{10Gy}$, were obtained, without systematic bias. Similar results were found for the independent dataset. Our novel, ML-based organ dose reconstruction method is not only accurate but also efficient, as the setup of a surrogate is no longer needed.


title: AIBench: An Agile Domain-specific Benchmarking Methodology and an AI Benchmark Suite

abstract: Domain-specific software and hardware co-design is encouraging as it is much easier to achieve efficiency for fewer tasks. Agile domain-specific benchmarking speeds up the process as it provides not only relevant design inputs but also relevant metrics, and tools. Unfortunately, modern workloads like Big data, AI, and Internet services dwarf the traditional one in terms of code size, deployment scale, and execution path, and hence raise serious benchmarking challenges. This paper proposes an agile domain-specific benchmarking methodology. Together with seventeen industry partners, we identify ten important end-to-end application scenarios, among which sixteen representative AI tasks are distilled as the AI component benchmarks. We propose the permutations of essential AI and non-AI component benchmarks as end-to-end benchmarks. An end-to-end benchmark is a distillation of the essential attributes of an industry-scale application. We design and implement a highly extensible, configurable, and flexible benchmark framework, on the basis of which, we propose the guideline for building end-to-end benchmarks, and present the first end-to-end Internet service AI benchmark. The preliminary evaluation shows the value of our benchmark suite---AIBench against MLPerf and TailBench for hardware and software designers, micro-architectural researchers, and code developers. The specifications, source code, testbed, and results are publicly available from the web site \url{http://www.benchcouncil.org/AIBench/index.html}.


title: Satellite reveals age and extent of oil palm plantations in Southeast Asia

abstract: In recent decades, global oil palm production has shown an abrupt increase, with almost 90% produced in Southeast Asia alone. Monitoring oil palm is largely based on national surveys and inventories or one-off mapping studies. However, they do not provide detailed spatial extent or timely updates and trends in oil palm expansion or age. Palm oil yields vary significantly with plantation age, which is critical for landscape-level planning. Here we show the extent and age of oil palm plantations for the year 2017 across Southeast Asia using remote sensing. Satellites reveal a total of 11.66 (+/- 2.10) million hectares (Mha) of plantations with more than 45% located in Sumatra. Plantation age varies from ~7 years in Kalimantan to ~13 in Insular Malaysia. More than half the plantations on Kalimantan are young (<7 years) and not yet in full production compared to Insular Malaysia where 45% of plantations are older than 15 years, with declining yields. For the first time, these results provide a consistent, independent, and transparent record of oil palm plantation extent and age structure, which are complementary to national statistics.


title: Study of the $h\gamma Z$ coupling at the ILC

abstract: We study the $h \gamma Z$ coupling, which is a loop induced coupling in the Standard Model (SM), to probe new physics. In a global fit based on the SM Effective Field Theory, measurement of the SM $h \gamma Z$ coupling can provide a very useful constraint, in particular for the precise determination of $hZZ$ and $hWW$ couplings. At the International Linear Collider (ILC), there are two direct ways to study the $h \gamma Z$ coupling: one is to measure the branching ratio of the $h \to \gamma Z$ decay and the other to measure the cross section for the $e^+e^- \to h \gamma$ process. We have performed a full simulation study of the $e^+e^- \to h \gamma$ process at the 250 GeV ILC, assuming 2 ab$^{-1}$ data collected by the International Large Detector (ILD). The expected 1$\sigma$ bound on the effective $h\gamma Z$ coupling ($\zeta_{AZ}$) combining measurements of the cross section for $e^+e^- \to h \gamma$ followed by $h \to b \bar{b}$ and the $h \to \gamma Z$ branching ratio is $-0.0015<\zeta_{AZ}<0.0015$. The expected significance for the signal cross section in the fully hadronic $h \to WW^*$ channel is 0.09 $\sigma$ for beam polarizations of $P(e^-,e^+)=(-80\%,+30\%)$.


title: Bloch-Siegert oscillations in the Rabi model with an amplitude-modulated driving field

abstract: We study the coherent dynamics of a qubit excited by an amplitude-modulated electromagnetic field under the Rabi resonance when the frequency of the low-frequency modulation field matches the Rabi frequency in the high-frequency field. Due to destructive interference of multiple photon processes at the ultrastrong coupling between the qubit and the low-frequency driving field, Rabi oscillations result exclusively from the Bloch-Siegert effect. It is directly observed in the time-resolved coherent dynamics as the Bloch-Siegert oscillation. In this case, triplets in Fourier spectra of the coherent response are transformed into doublets with the splitting between the lines equal to twice the Bloch-Siegert shift. These unusual properties are demonstrated in conditions of experiments with a nitrogen vacancy center in diamond.


title: Local Spectrum, Local Spectral radius,and Growth Conditions

abstract: In this paper we present complete description of the elements of Banach space with one-point spectrum. Some applications of these results are also given.


title: Globally generated vector bundles with $c_1 = 5$ on $\mathbb{P}^n$, $n \geq 4$

abstract: We complete the classification of globally generated vector bundles with small $c_1$ on projective spaces by treating the case $c_1 = 5$ on $\mathbb{P}^n$, $n \geq 4$ (the case $c_1 \leq 3$ has been considered by Sierra and Ugaglia, while the cases $c_1 = 4$ on any projective space and $c_1 = 5$ on $\mathbb{P}^2$ and $\mathbb{P}^3$ have been studied in two of our previous papers). It turns out that there are very few indecomposable bundles of this kind: besides some obvious examples there are, roughly speaking, only the (first twist of the) rank 5 vector bundle which is the middle term of the monad defining the Horrocks bundle of rank 3 on $\mathbb{P}^5$, and its restriction to $\mathbb{P}^4$. We recall, in an appendix, from our preprint [arXiv:1805.11336], the main results allowing the classification of globally generated vector bundles with $c_1 = 5$ on $\mathbb{P}^3$. Since there are many such bundles, a large part of the main body of the paper is occupied with the proof of the fact that, except for the simplest ones, they do not extend to $\mathbb{P}^4$ as globally generated vector bundles.


title: Density of Binary Disc Packings: The 9 Compact Packings

abstract: A disc packing in the plane is compact if its contact graph is a triangulation. There are $9$ values of $r$ such that a compact packing by discs of radii $1$ and $r$ exists. We prove, for each of these $9$ values, that the maximal density over all the packings by discs of radii $1$ and $r$ is reached for a compact packing (we give it as well as its density).


title: On commutative homogeneous vector bundles attached to nilmanifolds

abstract: The notion of Gelfand pair (G, K) can be generalized if we consider homogeneous vector bundles over G/K instead of the homogeneous space G/K and matrix-valued functions instead of scalar-valued functions. This gives the definition of commutative homogeneous vector bundles. Being a Gelfand pair is a necessary condition of being a commutative homogeneous vector bundle. For the case in which G/K is a nilmanifold having square-integrable representations, in a previous article we determined a big family of commutative homogeneous vector bundles. In this paper, we complete that classification.


title: The automorphism group of Rauzy diagrams

abstract: We give a description of the automorphism group of a Rauzy diagram as a subgroup of the symmetric group. This is based on an example that appear in some personnal notes of Yoccoz that are to be published in the project ''Yoccoz archives''.


title: Kalman meets Bellman: Improving Policy Evaluation through Value Tracking

abstract: Policy evaluation is a key process in Reinforcement Learning (RL). It assesses a given policy by estimating the corresponding value function. When using parameterized value functions, common approaches minimize the sum of squared Bellman temporal-difference errors and receive a point-estimate for the parameters. Kalman-based and Gaussian-processes based frameworks were suggested to evaluate the policy by treating the value as a random variable. These frameworks can learn uncertainties over the value parameters and exploit them for policy exploration. When adopting these frameworks to solve deep RL tasks, several limitations are revealed: excessive computations in each optimization step, difficulty with handling batches of samples which slows training and the effect of memory in stochastic environments which prevents off-policy learning. In this work, we discuss these limitations and propose to overcome them by an alternative general framework, based on the extended Kalman filter. We devise an optimization method, called Kalman Optimization for Value Approximation (KOVA) that can be incorporated as a policy evaluation component in policy optimization algorithms. KOVA minimizes a regularized objective function that concerns both parameter and noisy return uncertainties. We analyze the properties of KOVA and present its performance on deep RL control tasks.


title: Shape Optimization of Actuators over Banach Spaces for Nonlinear Systems

abstract: In this paper, optimal actuator shape for nonlinear parabolic systems is discussed. The system under study is an abstract differential equation with a locally Lipschitz nonlinear part. A quadratic cost on the state and input of the system is considered. The existence of an optimal actuator shape has been established in the literature. This paper focuses on driving the optimality conditions for actuator shapes belonging to a Banach space. The application of the theory to the optimal actuator shape design for railway track model is considered.


title: How Impersonators Exploit Instagram to Generate Fake Engagement?

abstract: Impersonators on Online Social Networks such as Instagram are playing an important role in the propagation of the content. These entities are the type of nefarious fake accounts that intend to disguise a legitimate account by making similar profiles. In addition to having impersonated profiles, we observed a considerable engagement from these entities to the published posts of verified accounts. Toward that end, we concentrate on the engagement of impersonators in terms of active and passive engagements which is studied in three major communities including ``Politician'', ``News agency'', and ``Sports star'' on Instagram. Inside each community, four verified accounts have been selected. Based on the implemented approach in our previous studies, we have collected 4.8K comments, and 2.6K likes across 566 posts created from 3.8K impersonators during 7 months. Our study shed light into this interesting phenomena and provides a surprising observation that can help us to understand better how impersonators engaging themselves inside Instagram in terms of writing Comments and leaving Likes.


title: The importance of finding the upper bounds for prime gaps in order to solve the twin primes conjecture and the Goldbach conjecture

abstract: ABSTRACT. In this article we present a point of view that highlights the importance of finding the upper bounds for prime gaps, in order to solve the twin primes conjecture and the Goldbach conjecture. For this purpose, we present a procedure for the determination of the upper bounds for prime gaps different from the most famous and known approaches. The proposed method analyzes the distribution of prime numbers using the set of relative numbers. Using negative numbers too, it becomes intuitive to understand that that the arrangement of 2P+1 consecutive numbers that goes -P to P, is the only arrangement that minimizes the distance between two powers having the same absolute value of the base D, with |D|<=P. This arrangement is considered important because by increasing the number of powers of the prime numbers within a range of consecutive numbers, it is presumed to decrease the overlap between the prime numbers considered. Consequently, by reducing these overlaps, we suppose to obtain an arrangement, in which the prime numbers less than and equal to P and their multiples occupy the greatest possible number of positions within a range of 2P+1 consecutive numbers. If this result could be demonstrated, would imply not only the resolution of the Legendre conjecture, but also a step forward in the resolution of the twin primes conjecture and the Goldbach conjecture.


title: A Survey on Feasibility and Suitability of Blockchain Techniques for the E-Voting Systems

abstract: In the second decade of the 21st century, blockchain definitely became one of the most trending computational technologies. This research aims to question the feasibility and suitability of using blockchain technology within e-voting systems, regarding both technical and non-technical aspects. In today's world, although the course of this spreading is considerably slow, several countries already use means of e-voting due to many social and economic reasons, which we further investigated. Nevertheless, the number of countries offering various e-government solutions, apart from e-voting, is significantly high. E-voting systems, naturally, require much more attention and assurance regarding potential security and anonymity issues, since voting is one of the few extremely critical governmental processes. Nevertheless, e-voting is not purely a governmental service, but many companies and nonprofit organizations would benefit the cost-efficiency, scalability, remote accessibility, and ease of use that it provides. Blockchain technology is claimed to be able to address some, obviously not all, important security concerns, including anonymity, confidentiality, integrity, and non-repudiation. The analysis results presented in this article mostly confirm these claims.


title: Transitive nilpotent Lie algebras of vector fields and their Tanaka prolongations

abstract: Transitive local Lie algebras of vector fields can be easily constructed from dilations of $\mathbb{R}^n$ associating with coordinates positive weights (give me a sequence of $n$ positive integers and I will give you a transitive nilpotent Lie algebra of vector fields on $\mathbb{R}^n$). It is interesting that all transitive nilpotent local Lie algebra of vector fields can be obtained as subalgebras of nilpotent algebras of this kind. Starting with a graded nilpotent Lie algebra one constructs graded parts of its Tanaka prolongations inductively as `derivations of degree 0, 1, etc. Of course, vector fields of weight $k$ with respect to the dilation define automatically derivations of weight $k$, so the Tanaka prolongation is in this case never finite. Are they all such derivations given by vector fields or there are additional `strange'ones? We answer this question. Except for special cases, derivations of degree 0 are given by vector fields of degree 0 and the Tanaka prolongation recovers the whole algebra of polynomial vectors defined by the dilation. However, in some particular cases of dilations we can find `strange' derivations which we describe in detail


title: Limit of Gaussian and normal curvatures of surfaces in Riemannian approximation scheme for sub-Riemannian three dimensional manifolds and Gauss-Bonnet theorem

abstract: The authors Balogh-Tyson-Vecchi in arXiv:1604.00180 utilize the Riemannian approximations scheme $(\mathbb H^1,<,>_L)$, in the Heisenberg group, introduced by Gromov, to calculate the limits of Gaussian and normal curvatures defined on surfaces of $\mathbb H^1$ when $L\rightarrow\infty$. They show that these limits exist (unlike the limit of Riemannian surface area form or length form), and they obtain Gauss-Bonnet theorem in $\mathbb H^1$ as limit of Gauss-Bonnet theorems in $(\mathbb H^1,<,>_L)$ when $L$ goes to infinity. This construction was extended by Wang-Wei in arXiv:1912.00302 to the affine group and the group of rigid motions of the Minkowski plane. We generalize constructions of both papers to surfaces in sub-Riemannian three dimensional manifolds following the approach of arXiv:1909.13341, and prove analogous Gauss-Bonnet theorem.


title: Energy-Efficiency Routing algorithms in Wireless Sensor Networks: a Survey

abstract: A Wireless Sensor Network (WSN) is a collection of tiny nodes that have low energy levels and have become an essential component of the modern communication infrastructure and very important in industry and academia. Energy is crucial in WSN, and thus the design of WSN in the research community is based on energy efficiency, and node energy consumption is a great challenge to enhance WSN's lifetime. It may be costly or even impossible to charge or replace consumed batteries because of the difficult environment. Many energy efficiency methods are introduced in this article to decrease energy consumption, improve network performance and increase network lifetime.


title: ScamPy -- A sub-halo clustering & abundance matching based Python interface for painting galaxies on the dark matter halo/sub-halo hierarchy

abstract: We present a computational framework for "painting" galaxies on top of the Dark Matter Halo/Sub-Halo hierarchy obtained from N-body simulations. The method we use is based on the sub-halo clustering and abundance matching (SCAM) scheme which requires observations of the 1- and 2-point statistics of the target (observed) population we want to reproduce. This method is particularly tailored for high redshift studies and thereby relies on the observed high-redshift galaxy luminosity functions and correlation properties. The core functionalities are written in c++ and exploit Object Oriented Programming, with a wide use of polymorphism, to achieve flexibility and high computational efficiency. In order to have an easily accessible interface, all the libraries are wrapped in python and provided with an extensive documentation. We validate our results and provide a simple and quantitative application to reionization, with an investigation of physical quantities related to the galaxy population, ionization fraction and bubble size distribution.


title: The ALMA Spectroscopic Survey in the HUDF: A model to explain observed 1.1 and 0.85 millimeter dust continuum number counts

abstract: We present a new semi-empirical model for the dust continuum number counts of galaxies at 1.1 millimeter and 850 \micron. Our approach couples an observationally motivated model for the stellar mass and SFR distribution of galaxies with empirical scaling relations to predict the dust continuum flux density of these galaxies. Without a need to tweak the IMF, the model reproduces the currently available observations of the 1.1 millimeter and 850 \micron number counts, including the observed flattening in the 1.1 millimeter number counts below 0.3 mJy \citep{Gonzalez2019numbercounts} and the number counts in discrete bins of different galaxy properties. Predictions of our work include : (1) the galaxies that dominate the number counts at flux densities below 1 mJy (3 mJy) at 1.1 millimeter (850 $\mu$m) have redshifts between $z=1$ and $z=2$, stellar masses of $\sim 5\times10^{10}~\rm{M}_\odot$, and dust masses of $\sim 10^{8}~\rm{M}_\odot$; (2) the flattening in the observed 1.1 millimeter number counts corresponds to the knee of the 1.1 millimeter luminosity function. A similar flattening is predicted for the number counts at 850 $\mu$m; (3) the model reproduces the redshift distribution of current 1.1 millimeter detections; (4) to efficiently detect large numbers of galaxies through their dust continuum, future surveys should scan large areas once reaching a 1.1 millimeter flux density of 0.1 mJy rather than integrating to fainter fluxes. Our modeling framework also suggests that the amount of information on galaxy physics that can be extracted from the 1.1 millimeter and 850 $\mu$m number counts is almost exhausted.


title: Linear resistivity and Sachdev-Ye-Kitaev (SYK) spin liquid behavior in a quantum critical metal with spin-$1/2$ fermions

abstract: `Strange metals' with resistivity depending linearly on temperature $T$ down to low-$T$ have been a long-standing puzzle in condensed matter physics. Here, we consider a model of itinerant spin-$1/2$ fermions interacting via on-site Hubbard interaction and random infinite-ranged spin-spin interaction. We show that the quantum critical point associated with the melting of the spin-glass phase by charge fluctuations displays non-Fermi liquid behaviour, with local spin dynamics identical to that of the Sachdev-Ye-Kitaev family of models. This extends the quantum spin liquid dynamics previously established in the large-$M$ limit of $SU(M)$ symmetric models, to models with physical $SU(2)$ spin-$1/2$ electrons. Remarkably, the quantum critical regime also features a Planckian linear-$T$ resistivity associated with a $T$-linear scattering rate and a frequency dependence of the electronic self-energy consistent with the Marginal Fermi Liquid phenomenology.


title: Sterile neutrinos and neutrinoless double beta decay in effective field theory

abstract: We investigate neutrinoless double beta decay ($0\nu\beta\beta$) in the presence of sterile neutrinos with Majorana mass terms. These gauge-singlet fields are allowed to interact with Standard-Model (SM) fields via renormalizable Yukawa couplings as well as higher-dimensional gauge-invariant operators up to dimension seven in the Standard Model Effective Field Theory extended with sterile neutrinos. At the GeV scale, we use Chiral effective field theory involving sterile neutrinos to connect the operators at the level of quarks and gluons to hadronic interactions involving pions and nucleons. This allows us to derive an expression for $0\nu\beta\beta$ rates for various isotopes in terms of phase-space factors, hadronic low-energy constants, nuclear matrix elements, the neutrino masses, and the Wilson coefficients of higher-dimensional operators. The needed hadronic low-energy constants and nuclear matrix elements depend on the neutrino masses, for which we obtain interpolation formulae grounded in QCD and chiral perturbation theory that improve existing formulae that are only valid in a small regime of neutrino masses. The resulting framework can be used directly to assess the impact of $0\nu\beta\beta$ experiments on scenarios with light sterile neutrinos and should prove useful in global analyses of sterile-neutrino searches. We perform several phenomenological studies of $0\nu\beta\beta$ in the presence of sterile neutrinos with and without higher-dimensional operators. We find that non-standard interactions involving sterile neutrinos have a dramatic impact on $0\nu\beta\beta$ phenomenology, and next-generation experiments can probe such interactions up to scales of $\mathcal O(100)$ TeV.


title: Tracing the assembly of the Milky Way's disk through abundance clustering

abstract: A major goal in the field of galaxy formation is to understand the formation of the Milky Way's disk. The first step toward doing this is to empirically describe its present state. We use the new high-dimensional dataset of 19 abundances from 27,135 red clump APOGEE stars to examine the distribution of clusters defined using abundances. We explore different dimensionality reduction techniques and implement a non-parametric agglomerate hierarchical clustering method. We see that groups defined using abundances are spatially separated, as a function of age. Furthermore, the abundance groups represent different distributions in the [Fe/H]-age plane. Ordering our clusters by age reveals patterns suggestive of the sequence of chemical enrichment in the disk over time. Our results indicate that a promising avenue to trace the details of the disk's assembly is via a full interpretation of the empirical connections we report.


title: Beta decays as sensitive probes of lepton flavor universality

abstract: Nuclear $\beta$ decays as well as the decay of the neutron are well-established low-energy probes of physics beyond the Standard Model (SM). In particular, with the axial-vector coupling of the nucleon $g_A$ determined from lattice QCD, the comparison between experiment and SM prediction is commonly used to derive constraints on right-handed currents. Further, in addition to the CKM element $V_{us}$ from kaon decays, $V_{ud}$ from $\beta$ decays is a critical input for the test of CKM unitarity. Here, we point out that the available information on $\beta$ decays can be re-interpreted as a stringent test of lepton flavor universality (LFU). In fact, we find that the ratio of $V_{us}$ from kaon decays over $V_{us}$ from $\beta$ decays (assuming CKM unitarity) is extremely sensitive to LFU violation (LFUV) in $W$-$\mu$-$\nu$ couplings thanks to a CKM enhancement by $(V_{ud}/V_{us})^2\sim 20$. From this perspective, recent hints for the violation of CKM unitarity can be viewed as further evidence for LFUV, fitting into the existing picture exhibited by semi-leptonic $B$ decays and the anomalous magnetic moments of muon and electron. Finally, we comment on the future sensitivity that can be reached with this LFU violating observable and discuss complementary probes of LFU that may reach a similar level of precision, such as $\Gamma(\pi\to\mu\nu)/\Gamma(\pi\to e\nu)$ at the PEN and PiENu experiments or even direct measurements of $W\to\mu\nu$ at an FCC-ee.


title: On the coexistence of the streaming instability and the vertical shear instability in protoplanetary disks

abstract: The streaming instability is a leading candidate mechanism to explain the formation of planetesimals. Yet, the role of this instability in the driving of turbulence in protoplanetary disks, given its fundamental nature as a linear hydrodynamical instability, has so far not been investigated in detail. We study the turbulence that is induced by the streaming instability as well as its interaction with the vertical shear instability. For this purpose, we employ the FLASH Code to conduct two-dimensional axisymmetric global disk simulations spanning radii from $1$ au to $100$ au, including the mutual drag between gas and dust as well as the radial and vertical stellar gravity. If the streaming instability and the vertical shear instability start their growth at the same time, we find the turbulence in the dust mid-plane layer to be primarily driven by the streaming instability. It gives rise to vertical gas motions with a Mach number of up to ${\sim}10^{-2}$. The dust scale height is set in a self-regulatory manner to about $1\%$ of the gas scale height. In contrast, if the vertical shear instability is allowed to saturate before the dust is introduced into our simulations, then it continues to be the main source of the turbulence in the dust layer. The vertical shear instability induces turbulence with a Mach number of ${\sim}10^{-1}$ and thus impedes dust sedimentation. Nonetheless, we find the vertical shear instability and the streaming instability in combination to lead to radial dust concentration in long-lived accumulations which are significantly denser than those formed by the streaming instability alone. Thus, the vertical shear instability may promote planetesimal formation by creating weak overdensities that act as seeds for the streaming instability.


title: Smoke on the water: CO and H2O in a circumnuclear disc around a quasar at redshift 2.64

abstract: We present an analysis of observations with the Atacama Large (sub-)Millimetre Array (ALMA) of the 380 GHz water emission line and CO (11-10) emission line from MG J0414+0534, a gravitationally lensed dusty star-forming galaxy that hosts a type 1 quasar. We also present observations at 1.6 GHz with global very long baseline interferometry (VLBI) of the radio source. We confirm the previously reported detection of the 380 GHz water line and that the flux density ratio between the two merging lensed images is reversed with respect to the radio/mm continuum. We further find tentative evidence of variability in the integrated line intensity on timescales of days. We show that the 380 GHz water line has two components of emission: a disc that is around 35 pc in diameter around the quasar, and another component of emission that is offset $\sim$600 pc perpendicular to the disc that lies close to the lensing caustic. With lens modelling of the multi-wavelength data sets, we construct a model for the quasar system consisting of a circumnuclear disc of molecular gas with a size of about 60 pc bisected by radio jets extending to a distance of about 200 pc from the radio core. Our findings suggest that observations with ALMA of high-excitation molecular lines from strongly lensed quasars could allow detailed studies of AGN accretion and feedback at the cosmic peak of black hole and galaxy growth.


title: Hyper-Eddington accretion flows onto black holes accompanied by powerful outflows

abstract: We perform two-dimensional radiation hydrodynamical simulations of accretion flows onto black holes (BHs) at the nuclei of protogalaxies, and study the impact of mechanical and radiative feedback on rapid growth of BHs. The outflows deposit mass, momentum and energy into the surrounding medium and prevent mass accretion onto the BH, resulting in the reduction of radiative output. We find that when the BH is embedded in a dense gas core, ionizing radiation attenuated by inefficient BH feeding owing to mechanical feedback hardly affects the gas dynamics at the BH gravitational sphere of influence, from which intense inflows of neutral gas occur at rates substantially exceeding the Eddington limit without impeded by photoionization and heating. Since mechanical power of outflows driven by the rapidly accreting BH is sufficiently strong, bipolar outflows completely evacuate the surrounding gas in the polar region but mass inflows through the equatorial region maintain the BH accretion rate as high as $\sim 300-10^3~\dot{M}_{\rm Edd}$, which is reduced by one order of magnitude from those with radiative feedback alone. Furthermore, we find that the critical gas density required for rapid accretion is lower by a factor of $\sim 3$ nearly independently of BH mass, when mechanical feedback is considered. By studying the dependence on outflow model parameters (e.g., opening angle, mass-loading degree into outflows, velocity), we conclude that contrary to naive expectation, the existence of stronger outflow leads to the transition to rapid accretion phases more efficiently. Rapidly growing BHs inject mechanical power with $\sim 0.1-1\%$ of the radiative luminosity into their host galaxy scales, which is used for cosmological simulations.


title: How efficient is the streaming instability in viscous protoplanetary disks?

abstract: The streaming instability is a popular candidate for planetesimal formation by concentrating dust particles to trigger gravitational collapse. However, its robustness against physical conditions expected in protoplanetary disks is unclear. In particular, particle stirring by turbulence may impede the instability. To quantify this effect, we develop the linear theory of the streaming instability with external turbulence modelled by gas viscosity and particle diffusion. We find the streaming instability is sensitive to turbulence, with growth rates becoming negligible for alpha-viscosity parameters $\alpha \gtrsim \mathrm{St} ^{1.5}$, where $\mathrm{St}$ is the particle Stokes number. We explore the effect of non-linear drag laws, which may be applicable to porous dust particles, and find growth rates are modestly reduced. We also find that gas compressibility increase growth rates by reducing the effect of diffusion. We then apply linear theory to global models of viscous protoplanetary disks. For minimum-mass Solar nebula disk models, we find the streaming instability only grows within disk lifetimes beyond $\sim 10$s of AU, even for cm-sized particles and weak turbulence ($\alpha\sim 10^{-4}$). Our results suggest it is rather difficult to trigger the streaming instability in non-laminar protoplanetary disks, especially for small particles.


title: On the (Lack of) Evolution of the Stellar Mass Function of Massive Galaxies from $z$=1.5 to 0.4

abstract: We study the evolution in the number density of the highest mass galaxies over $0.4<z<1.5$ (covering 9 Gyr). We use the Spitzer/HETDEX Exploratory Large-Area (SHELA) Survey, which covers 17.5 $\mathrm{deg}^2$ with eight photometric bands spanning 0.3-4.5 $\mu$m within the SDSS Stripe 82 field. This size produces the lowest counting uncertainties and cosmic variance yet for massive galaxies at $z\sim1.0$. We study the stellar mass function (SMF) for galaxies with $\log(M_\ast/M_\odot)>10.3$ using a forward-modeling method that fully accounts for statistical and systematic uncertainties on stellar mass. From $z$=0.4 to 1.5 the massive end of the SMF shows minimal evolution in its shape: the characteristic mass ($M^\ast$) evolves by less than 0.1 dex ($\pm$0.05 dex); the number density of galaxies with $\log (M_\ast/M_\odot) >11$ stays roughly constant at $\log (n/\mathrm{Mpc}^{-3})$ $\simeq$ $-$3.4 ($\pm$0.05), then declines to $\log (n/\mathrm{Mpc}^{-3})$=$-$3.7 ($\pm$0.05) at $z$=1.5. We discuss the uncertainties in the SMF, which are dominated by assumptions in the star formation history and details of stellar population synthesis models for stellar mass estimations. For quiescent galaxies, the data are consistent with no (or slight) evolution ($\lesssim0.1$ dex) in the characteristic mass nor number density from $z\sim 1.5$ to the present. This implies that any mass growth (presumably through "dry" mergers) of the quiescent massive galaxy population must balance the rate of mass losses from late-stage stellar evolution and the formation of quenching galaxies from the star-forming population. We provide a limit on this mass growth from $z=1.0$ to 0.4 of $\Delta M_\ast/M_\ast\leq$ 45% (i.e., $\simeq0.16$ dex) for quiescent galaxies more massive than $10^{11}$ $M_\odot$.


title: Enhancement of magnetization plateaus in low dimensional spin systems

abstract: We study the low-energy properties and, in particular, the magnetization process of a spin-1/2 Heisenberg $J_1-J_2$ sawtooth and frustrated chain (also known as zig-zag ladder) with a spatially anisotropic $g$-factor. We treat the problem both analytically and numerically while keeping the $J_2/J_1$ ratio generic. Numerically, we use complete and Lanczos diagonalization as well as the infinite time-evolving block decimation (iTEBD) method. Analytically we employ (non-)Abelian bosonization. Additionally for the sawtooth chain, we provide an analytical description in terms of flat bands and localized magnons. By considering a specific pattern for the $g$-factor anisotropy for both models, we show that a small anisotropy significantly enhances a magnetization plateau at half saturation. For the magnetization of the frustrated chain, we show the destruction of the $1/3$ of the full saturation plateau in favor of the creation of a plateau at half-saturation. For large anisotropies, the existence of an additional plateau at zero magnetization is possible. Here and at higher magnetic fields, the system is locked in the half-saturation plateau, never reaching full saturation.


title: MNELLS: The MUSE Nearby Early-Type Galaxy Lens Locator Survey

abstract: Low-redshift strong-lensing galaxies can provide robust measurements of the stellar mass-to-light ratios in early-type galaxies (ETG), and hence constrain variations in the stellar initial mass function (IMF). At present, only a few such systems are known. Here, we report the first results from a blind search for gravitationally-lensed emission line sources behind 52 massive $z$ $<$ 0.07 ETGs with MUSE integral field spectroscopy. For 16 galaxies, new observations were acquired, whilst the other 36 were analysed from archival data. This project has previously yielded one confirmed galaxy-scale strong lens (J0403-0239) which we report in an earlier paper. J0403-0239 has since received follow-up observations, presented here, which indicate support for our earlier IMF results. Three cluster-scale, and hence dark-matter-dominated, lensing systems were also discovered (central galaxies of A4059, A2052 and AS555). For nine further galaxies, we detect a singly-imaged but closely-projected source within 6 arcsec (including one candidate with sources at three different redshifts); such cases can be exploited to derive upper limits on the IMF mass-excess factor, $\alpha$. Combining the new lens and new upper limits, with the previously-discovered systems, we infer an average $\langle \alpha \rangle$ = 1.06 $\pm$ 0.08 (marginalised over the intrinsic scatter), which is inconsistent with a Salpeter-like IMF ($\alpha$ = 1.55) at the 6$\sigma$ level. We test the detection threshold in these short-exposure MUSE observations with the injection and recovery of simulated sources, and predict that one in twenty-five observations is expected to yield a new strong-lens system. Our observational results are consistent with this expected yield.


title: New constraints on the origin of medium-energy neutrinos observed by IceCube

abstract: The recent IceCube publication claims the observation of cosmic neutrinos with energies down to $\sim 10$ TeV, reinforcing the growing evidence that the neutrino flux in the 10-100 TeV range is unexpectedly large. Any conceivable source of these neutrinos must also produce a $\gamma$-ray flux which degrades in energy en route to the Earth and contributes to the extragalactic $\gamma$-ray background measured by the Fermi satellite. In a quantitative multimessenger analysis, featuring minimalistic assumptions, we find a $\geq 3\sigma$ tension in the data, reaching $\sim 5\sigma$ for cosmic neutrinos extended down to $\sim 1$ TeV, interpreted as evidence for a population of hidden cosmic-ray accelerators.


title: Room-Temperature Photonic Logical Qubits via Second-Order Nonlinearities

abstract: Recent progress in nonlinear optical materials and microresonators has brought quantum computing with bulk optical nonlinearities into the realm of possibility. This platform is of great interest, not only because photonics is an obvious choice for quantum networks, but also because it may be the only feasible route to quantum information processing at room temperature. We introduce a paradigm for room-temperature photonic quantum logic that significantly simplifies the realization of various quantum circuits, and in particular, of error correction. It uses only the strongest available bulk nonlinearity, namely the $\chi^{(2)}$ nonlinear susceptibility. The key element is a three-mode resonator that implements programmable bosonic quantum logic gates. We show that just two of these elements suffice for a complete, compact error-correction circuit on a bosonic code, without the need for measurement or feed-forward control. An extrapolation of current progress in nonlinear optical materials and photonic circuits indicates that such circuitry should be achievable within the next decade.


title: Holographic renormalization group flows in two-dimensional gravity and $AdS$ black holes

abstract: We look into the $AdS$ black holes from two-dimensional gravity perspective. In this work, we extend the previous results of holographic renormalization group flows to dimensions two. By introducing a superpotential, we derive the flow equations in two-dimensional dilaton gravity. We find a quantity which monotonically decreases along flows and give some comments on holographic $c$-theorem. As examples, we show that recently studied supersymmetric $AdS$ black hole solutions generically dimensionally reduce to two-dimensional dilaton gravity, and obtain the flow equations for black hole solutions.


title: Constraining the intergalactic medium at $z\approx$ 9.1 using LOFAR Epoch of Reionization observations

abstract: We derive constraints on the thermal and ionization states of the intergalactic medium (IGM) at redshift $\approx$ 9.1 using new upper limits on the 21-cm power spectrum measured by the LOFAR radio-telescope and a prior on the ionized fraction at that redshift estimated from recent cosmic microwave background (CMB) observations. We have used results from the reionization simulation code GRIZZLY and a Bayesian inference framework to constrain the parameters which describe the physical state of the IGM. We find that, if the gas heating remains negligible, an IGM with ionized fraction $\gtrsim 0.13$ and a distribution of the ionized regions with a characteristic size $\gtrsim 8 ~h^{-1}$ comoving megaparsec (Mpc) and a full width at the half maximum (FWHM) $\gtrsim 16 ~h^{-1}$ Mpc is ruled out. For an IGM with a uniform spin temperature $T_{\rm S} \gtrsim 3$ K, no constraints on the ionized component can be computed. If the large-scale fluctuations of the signal are driven by spin temperature fluctuations, an IGM with a volume fraction $\lesssim 0.34$ of heated regions with a temperature larger than CMB, average gas temperature 7-160 K and a distribution of the heated regions with characteristic size 3.5-70 $h^{-1}$ Mpc and FWHM of $\lesssim 110$ $h^{-1}$ Mpc is ruled out. These constraints are within the 95 per cent credible intervals. With more stringent future upper limits from LOFAR at multiple redshifts, the constraints will become tighter and will exclude an increasingly large region of the parameter space.


title: Improved upper limits on the 21-cm signal power spectrum of neutral hydrogen at $\boldsymbol{z \approx 9.1}$ from LOFAR

abstract: A new upper limit on the 21-cm signal power spectrum at a redshift of $z \approx 9.1$ is presented, based on 141 hours of data obtained with the Low-Frequency Array (LOFAR). The analysis includes significant improvements in spectrally-smooth gain-calibration, Gaussian Process Regression (GPR) foreground mitigation and optimally-weighted power spectrum inference. Previously seen `excess power' due to spectral structure in the gain solutions has markedly reduced but some excess power still remains with a spectral correlation distinct from thermal noise. This excess has a spectral coherence scale of $0.25 - 0.45$\,MHz and is partially correlated between nights, especially in the foreground wedge region. The correlation is stronger between nights covering similar local sidereal times. A best 2-$\sigma$ upper limit of $\Delta^2_{21} < (73)^2\,\mathrm{mK^2}$ at $k = 0.075\,\mathrm{h\,cMpc^{-1}}$ is found, an improvement by a factor $\approx 8$ in power compared to the previously reported upper limit. The remaining excess power could be due to residual foreground emission from sources or diffuse emission far away from the phase centre, polarization leakage, chromatic calibration errors, ionosphere, or low-level radio-frequency interference. We discuss future improvements to the signal processing chain that can further reduce or even eliminate these causes of excess power.


title: Quantum Phase Transition induced by Topological Frustration

abstract: In quantum many-body systems with local interactions, the effects of boundary conditions are considered to be negligible, at least for sufficiently large systems. Here we show an example of the opposite. We consider a spin chain with two competing interactions, set on a ring with an odd number of sites. When only the dominant interaction is antiferromagnetic, and thus induces topological frustration, the standard antiferromagnetic order (expressed by the magnetization) is destroyed. When also the second interaction turns from ferro to antiferro, an antiferromagnetic order characterized by a site-dependent magnetization which varies in space with an incommensurate pattern, emerges. This modulation results from a ground state degeneracy, which allows to break the translational invariance. The transition between the two cases is signaled by a discontinuity in the first derivative of the ground state energy and represents a quantum phase transition induced by a special choice of boundary conditions.

categories: cond-mat.stat-mech cond-mat.str-el hep-th quant-ph


title: The Chandra High Resolution X-ray Spectrum of Quiescent Emission from Sgr A*

abstract: In quiescence, Sgr A* is surprisingly dim, shining 100,000 times less than expected for its environment. This problem has motivated a host of theoretical models to explain radiatively inefficient accretion flows (RIAFs). The Chandra Galactic Center (GC) X-ray Visionary Program obtained approximately 3 Ms (one month) of Chandra HETG data, offering the only opportunity to examine the quiescent X-ray emission of Sgr A* with high resolution spectroscopy. Utilizing custom background regions and filters for removing overlapping point sources, this work provides the first ever look at stacked HETG spectra of Sgr A*. We model the background datasets with a cubic spline and fit the unbinned Sgr A* spectra with a simple parametric model of a power law plus Gaussian lines under the effects of interstellar extinction. We detect a strong 6.7 keV iron emission line in the HEG spectra and a 3.1 keV emission line in the MEG spectra. In all cases, the line centroids and equivalent widths are consistent with those measured from low-resolution CCD spectra. An examination of the unbinned, stacked HEG+/-1 spectrum reveals fine structure in the iron line complex. In addition to resolving the resonant and forbidden lines from He-like iron, there are apparent emission features arising with higher statistical significance at lower energy, potentially associated with FeXX-XXIV ions in a ~1 keV plasma arising near the Bondi radius of Sgr A*. With this work, we release the cleaned and stacked Sgr A* and background HETG spectra to the public as a special legacy dataset.


title: The ALMA Spectroscopic Survey in the HUDF: Deep 1.2 mm continuum number counts

abstract: We present the results from the 1.2 mm continuum image obtained as part of the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS). The 1.2 mm continuum image has a size of 2.9 (4.2) arcmin$^2$ within a primary beam response of 50% (10%) and a rms value of $9.3\thinspace{\rm\mu Jy\thinspace beam^{-1}}$. We detect 35 sources at high significance (Fidelity $\geq0.5$), 32 of these have well characterized near-infrared HST counterparts. We estimate the 1.2 mm number counts to flux levels of $<30\thinspace{\rm\mu Jy}$ in two different ways: we first use the detected sources to constrain the number counts and find a significant flattening of the counts below $S_\nu \sim 0.1$ mJy. In a second approach, we constrain the number counts by using a probability of deflection statistics (P(D)) analysis. For this latter approach, we describe new methods to accurately measure the noise in interferometric imaging (employing jack-knifing in the cube and in the visibility plane). This independent measurement confirms the flattening of the number counts. Our analysis of the differential number counts shows that we are detecting $\sim$93% ($\sim$100% if we include the lower fidelity detections) of the total continuum dust emission associated to galaxies in the HUDF. The ancillary data allows us to study the dependence of the 1.2 mm number counts on redshift ($z=0-4$), galaxy dust mass (${\rm M}_{\rm dust}=10^{7}-10^{9}{\rm M}_{\odot}$), stellar mass (${\rm M}_{*}=10^{9}-10^{12}{\rm M}_{\odot}$), and star-formation rate (${\rm SFR}=1-1000\thinspace{\rm M}_{\odot}\thinspace{\rm yr^{-1}}$). In an accompanying paper we show that the number counts are crucial to constrain galaxy evolution models and the understanding of star-forming galaxies at high redshift.


title: Holographic vector meson melting in a thermal gravity-dilaton background related to QCD

abstract: A holographic model of probe vector mesons (quarkonia) is presented, where the dynamical gravity-dilaton background is adjusted to the thermodynamics of 2 +1 flavor QCD with physical quark masses. The vector meson action is modified to account for various quark masses. We focus on the $\Phi$, $J/\psi$ and $\Upsilon$ meson melting in agreement with hadron phenomenology in heavy-ion collisions at LHC, that is the formation of hadrons at the observed freeze-out temperature of 155 MeV.


title: Stochastic evolution of scalar fields with continuous symmetries during inflation

abstract: During inflation, scalar fields with masses less than the Hubble scale acquire vacuum expectation values (vevs) via stochastic processes driven by quantum fluctuations. For nearly massless spectator scalars transforming nontrivially under a continuous symmetry group, we demonstrate that the evolution of the vev depends on the dimensionality of the scalar field space. Fields in larger representations both attain larger vacuum expectation values and converge more rapidly to equilibrium. We present an argument demonstrating how this higher-dimensional evolution can be obtained in unitary gauge for fields transforming under local symmetries with a mass gap that is small compared to the Hubble scale. Finally, we show that accounting for the full number of degrees of freedom in the Standard Model Higgs multiplet tightens Higgs stability constraints on the inflationary scale at the percent level and has more dramatic consequences for both the vev and the energy stored in the Higgs field after inflation.


title: Strings of diquark-quark (QQ)Q baryon before phase transition

abstract: We explore the limit at which the effective baryonic Y-string model of the junction approaches the mesonic stringlike behavior. We calculate and compare the numerical values of the static potential and energy-density correlators of diquark-quark and quark-antiquark configurations. The gauge model is pure Yang-Mills $SU(3)$ lattice gauge theory at coupling $\beta=6.0$ and finite temperature. The diquark setup is approximated as two quarks confined within a sphere of radius $0.1$ fm. The lattice data of the potential and energy show that the string binding the diquark-quark configuration displays an identical behavior to the quark-antiquark confining string. However, with the temperature increase to a small enough neighborhood of the critical point $T_{c}$, the gluonic similarities between the two systems do not manifest neither at short nor intermediate distance scales $R<1.0$ fm. The comparison between the potential and the second moment of the action-density correlators for both systems shows significant splitting. This suggests that subsisted baryonic decoupled states overlap with the mesonic spectrum. The baryonic junction's model for the potential and the profile returns a good fit to the numerical lattice data of the diquark-quark arrangement. However, near the critical point, the mesonic string displays large deviations compared to fits of the corresponding quark-antiquark data.


title: Multilinear Compressive Learning with Prior Knowledge

abstract: The recently proposed Multilinear Compressive Learning (MCL) framework combines Multilinear Compressive Sensing and Machine Learning into an end-to-end system that takes into account the multidimensional structure of the signals when designing the sensing and feature synthesis components. The key idea behind MCL is the assumption of the existence of a tensor subspace which can capture the essential features from the signal for the downstream learning task. Thus, the ability to find such a discriminative tensor subspace and optimize the system to project the signals onto that data manifold plays an important role in Multilinear Compressive Learning. In this paper, we propose a novel solution to address both of the aforementioned requirements, i.e., How to find those tensor subspaces in which the signals of interest are highly separable? and How to optimize the sensing and feature synthesis components to transform the original signals to the data manifold found in the first question? In our proposal, the discovery of a high-quality data manifold is conducted by training a nonlinear compressive learning system on the inference task. Its knowledge of the data manifold of interest is then progressively transferred to the MCL components via multi-stage supervised training with the supervisory information encoding how the compressed measurements, the synthesized features, and the predictions should be like. The proposed knowledge transfer algorithm also comes with a semi-supervised adaption that enables compressive learning models to utilize unlabeled data effectively. Extensive experiments demonstrate that the proposed knowledge transfer method can effectively train MCL models to compressively sense and synthesize better features for the learning tasks with improved performances, especially when the complexity of the learning task increases.


title: A theoretical investigation of the Humphreys-Davidson limit at high and low metallicity

abstract: Current massive star evolution grids are not able to simultaneously reproduce the empirical upper luminosity limit of red supergiants, the Humphrey-Davidson (HD) limit at high and low metallicity. In this study, we provide a better understanding of what drives massive star evolution to blue and red supergiant phases, with the ultimate aim of reproducing the HD limit at varied metallicities (Z). For solar, LMC, and SMC Z, we develop eight grids of MESA models for the mass range 20-60M to probe the effect of semiconvection and overshooting. We compare rotating and non-rotating models with efficient (alpha_semi = 100) and inefficient semi-convection (alpha_semi = 0.1), with high and low core overshooting (alpha_ov of 0.1 or 0.5). The red and blue supergiant evolutionary phases are investigated by comparing the fraction of core He-burning lifetimes spent in each phase. We find that the extension of the convective core by overshooting alpha_ov = 0.5 has an effect on the post-MS evolution which can disable semiconvection leading to more RSGs, but a lack of BSGs. We therefore implement alpha_ov = 0.1 which switches on semiconvective mixing, though for standard alpha_semi = 1, would result in an HD limit which is higher than observed at low Z. Therefore, we need to implement very efficient semiconvection of alpha_semi = 100 which reproduces the HD limit at log L ~ 5.5 for the Magellanic Clouds while simultaneously reproducing the Galactic HD limit of log L ~ 5.8 naturally. The effect of semiconvection is not active at high Z due to the depletion of the envelope structure by strong mass loss such that semiconvective regions could not form. Z-dependent mass loss plays an indirect, yet decisive role in setting the HD limit as a function of Z. For a combination of efficient semiconvection and low overshooting with standard Z-dependent mass loss, we find a natural HD limit at all metallicities.


title: Lipschitz Extensions and Approximations

abstract: The classical Hahn-Banach theorem is based on a successive point-by-point procedure of extending bounded linear functionals. In the setting of a general metric domain, the conditions are less restrictive and the extension is only required to be Lipschitz with the same Lipschitz constant. In this case, the successive procedure can be replaced by a much simpler one which was done by McShane and Whitney in the 1930s. Using virtually the same construction, Czipszer and Geh\'er showed a similar extension property for pointwise Lipschitz functions. In the present paper, we relate this construction to another classical result obtained previously by Hausdorff and dealing with pointwise Lipschitz approximations of semi-continuous functions. Moreover, we furnish complementary extension-approximation results for locally Lipschitz functions which fit naturally in this framework.


title: Ripple Walk Training: A Subgraph-based training framework for Large and Deep Graph Neural Network

abstract: Graph neural networks (GNNs) have achieved outstanding performance in learning graph-structured data and various tasks. However, many current GNNs suffer from three common problems when facing large-size graphs or using a deeper structure: neighbors explosion, node dependence, and oversmoothing. Such problems attribute to the data structures of the graph itself or the designing of the multi-layers GNNs framework, and can lead to low training efficiency and high space complexity. To deal with these problems, in this paper, we propose a general subgraph-based training framework, namely Ripple Walk Training (RWT), for deep and large graph neural networks. RWT samples subgraphs from the full graph to constitute a mini-batch, and the full GNN is updated based on the mini-batch gradient. We analyze the high-quality subgraphs to train GNNs in a theoretical way. A novel sampling method Ripple Walk Sampler works for sampling these high-quality subgraphs to constitute the mini-batch, which considers both the randomness and connectivity of the graph-structured data. Extensive experiments on different sizes of graphs demonstrate the effectiveness and efficiency of RWT in training various GNNs (GCN & GAT).


title: Semi-dynamic Algorithms for Strongly Chordal Graphs

abstract: There is an extensive literature on dynamic algorithms for a large number of graph theoretic problems, particularly for all varieties of shortest path problems. Germane to this paper are a number fully dynamic algorithms that are known for chordal graphs. However, to the best of our knowledge no study has been done for the problem of dynamic algorithms for strongly chordal graphs. To address this gap, in this paper, we propose a semi-dynamic algorithm for edge-deletions and a semi-dynamic algorithm for edge-insertions in a strongly chordal graph, $G = (V, E)$, on $n$ vertices and $m$ edges. The query complexity of an edge-deletion is $O(d_u^2d_v^2 (n + m))$, where $d_u$ and $d_v$ are the degrees of the vertices $u$ and $v$ of the candidate edge $\{u, v\}$, while the query-complexity of an edge-insertion is $O(n^2)$.


title: Optimal Error Pseudodistributions for Read-Once Branching Programs

abstract: In a seminal work, Nisan (Combinatorica'92) constructed a pseudorandom generator for length $n$ and width $w$ read-once branching programs with seed length $O(\log n\cdot \log(nw)+\log n\cdot\log(1/\varepsilon))$ and error $\varepsilon$. It remains a central question to reduce the seed length to $O(\log (nw/\varepsilon))$, which would prove that $\mathbf{BPL}=\mathbf{L}$. However, there has been no improvement on Nisan's construction for the case $n=w$, which is most relevant to space-bounded derandomization. Recently, in a beautiful work, Braverman, Cohen and Garg (STOC'18) introduced the notion of a pseudorandom pseudo-distribution (PRPD) and gave an explicit construction of a PRPD with seed length $\tilde{O}(\log n\cdot \log(nw)+\log(1/\varepsilon))$. A PRPD is a relaxation of a pseudorandom generator, which suffices for derandomizing $\mathbf{BPL}$ and also implies a hitting set. Unfortunately, their construction is quite involved and complicated. Hoza and Zuckerman (FOCS'18) later constructed a much simpler hitting set generator with seed length $O(\log n\cdot \log(nw)+\log(1/\varepsilon))$, but their techniques are restricted to hitting sets. In this work, we construct a PRPD with seed length $$O(\log n\cdot \log (nw)\cdot \log\log(nw)+\log(1/\varepsilon)).$$ This improves upon the construction in [BCG18] by a $O(\log\log(1/\varepsilon))$ factor, and is optimal in the small error regime. In addition, we believe our construction and analysis to be simpler than the work of Braverman, Cohen and Garg.


title: Enhanced topological protection in planar quasi-one-dimensional channels with periodically-modulated width

abstract: We study one dimensional (1D) and quasi-1D periodic structures as possible platforms for the emergence of Majorana bound states with enhanced robustness against disorder and system inhomogeneity. First, using a simple 1D model, we analytically derive the effective parameters characterizing the minibands generated by the periodic potential. We show that, for strong enough periodic potentials, the higher energy minibands hosting Majorana bound states have significant advantages compared to their counterparts in uniform systems, including increased topological gaps, enhanced robustness against disorder, and enlarged parameter space regions consistent with the presence of topological superconductivity. We identify the problem of engineering a strong enough periodic potential as a key roadblock to realizing efficient periodic 1D structures. To address this challenge, we propose an efficient implementation of the periodic potential based on quasi-1D channels realized in 2D semiconductor heterostructures proximity coupled to superconductor strips of periodically modulated width. Our numerical study of the modulated channel device shows excellent agreement with the simple 1D model, reveals a topological phase diagram that is quite insensitive to the details of the confining potential associated with screening by the superconductor, and demonstrates that engineering patterned 2D structures represents a powerful and versatile approach to realizing robust Majorana bound states.


title: Positive Hermitian curvature flow on complex 2-step nilpotent Lie groups

abstract: We study the positive Hermitian curvature flow of left-invariant metrics on complex 2-step nilpotent Lie groups. In this setting we completely characterize the long-time behaviour of the flow, showing that normalized solutions to the flow subconverge to a non-flat algebraic soliton, in Cheeger- Gromov topology. We also exhibit a uniqueness result for algebraic solitons on such Lie groups.


title: Spectrum preserving short cycle removal on regular graphs

abstract: We describe a new method to remove short cycles on regular graphs while maintaining spectral bounds (the nontrivial eigenvalues of the adjacency matrix), as long as the graphs have certain combinatorial properties. These combinatorial properties are related to the number and distance between short cycles and are known to happen with high probability in uniformly random regular graphs. Using this method we can show two results involving high girth spectral expander graphs. First, we show that given $d \geq 3$ and $n$, there exists an explicit distribution of $d$-regular $\Theta(n)$-vertex graphs where with high probability its samples have girth $\Omega(\log_{d - 1} n)$ and are $\epsilon$-near-Ramanujan; i.e., its eigenvalues are bounded in magnitude by $2\sqrt{d - 1} + \epsilon$ (excluding the single trivial eigenvalue of $d$). Then, for every constant $d \geq 3$ and $\epsilon > 0$, we give a deterministic poly$(n)$-time algorithm that outputs a $d$-regular graph on $\Theta(n)$-vertices that is $\epsilon$-near-Ramanujan and has girth $\Omega(\sqrt{\log n})$, based on the work of arXiv:1909.06988 .


title: Computing Classical Modular Forms for Arbitrary Congruence Subgroups

abstract: In this paper, we prove the existence of an efficient algorithm for the computation of $q$-expansions of modular forms of weight $k$ and level $\Gamma$, where $\Gamma \subseteq SL_{2}({\mathbb{Z}})$ is an arbitrary congruence subgroup. We also discuss some practical aspects and provide the necessary theoretical background.


title: Extreme Nonuniqueness of End-Sum

abstract: We give explicit examples of pairs of one-ended, open 4-manifolds whose end-sums yield uncountably many manifolds with distinct proper homotopy types. This answers strongly in the affirmative a conjecture of Siebenmann regarding the nonuniqueness of end-sums. In addition to the construction of these examples, we provide a detailed discussion of the tools used to distinguish them; most importantly, the end-cohomology algebra. Key to our Main Theorem is an understanding of this algebra for an end-sum in terms of the algebras of the summands together with ray-fundamental classes determined by the rays used to perform the end-sum. Differing ray-fundamental classes allow us to distinguish the various examples, but only through the subtle theory of infinitely generated abelian groups. An appendix is included which contains the necessary background from that area.


title: Robust Stochastic Bandit Algorithms under Probabilistic Unbounded Adversarial Attack

abstract: The multi-armed bandit formalism has been extensively studied under various attack models, in which an adversary can modify the reward revealed to the player. Previous studies focused on scenarios where the attack value either is bounded at each round or has a vanishing probability of occurrence. These models do not capture powerful adversaries that can catastrophically perturb the revealed reward. This paper investigates the attack model where an adversary attacks with a certain probability at each round, and its attack value can be arbitrary and unbounded if it attacks. Furthermore, the attack value does not necessarily follow a statistical distribution. We propose a novel sample median-based and exploration-aided UCB algorithm (called med-E-UCB) and a median-based $\epsilon$-greedy algorithm (called med-$\epsilon$-greedy). Both of these algorithms are provably robust to the aforementioned attack model. More specifically we show that both algorithms achieve $\mathcal{O}(\log T)$ pseudo-regret (i.e., the optimal regret without attacks). We also provide a high probability guarantee of $\mathcal{O}(\log T)$ regret with respect to random rewards and random occurrence of attacks. These bounds are achieved under arbitrary and unbounded reward perturbation as long as the attack probability does not exceed a certain constant threshold. We provide multiple synthetic simulations of the proposed algorithms to verify these claims and showcase the inability of existing techniques to achieve sublinear regret. We also provide experimental results of the algorithm operating in a cognitive radio setting using multiple software-defined radios.


title: STANNIS: Low-Power Acceleration of Deep Neural Network Training Using Computational Storage

abstract: This paper proposes a framework for distributed, in-storage training of neural networks on clusters of computational storage devices. Such devices not only contain hardware accelerators but also eliminate data movement between the host and storage, resulting in both improved performance and power savings. More importantly, this in-storage processing style of training ensures that private data never leaves the storage while fully controlling the sharing of public data. Experimental results show up to 2.7x speedup and 69% reduction in energy consumption and no significant loss in accuracy.


title: Three-point Functions in $\mathcal{N}=4$ SYM at Finite $N_c$ and Background Independence

abstract: We compute non-extremal three-point functions of scalar operators in $\mathcal{N}=4$ super Yang-Mills at tree-level in $g_{YM}$ and at finite $N_c$, using the operator basis of the restricted Schur characters. We make use of the diagrammatic methods called quiver calculus to simplify the three-point functions. The results involve an invariant product of the generalized Racah-Wigner tensors ($6j$ symbols). Assuming that the invariant product is written by the Littlewood-Richardson coefficients, we show that the non-extremal three-point functions satisfy the large $N_c$ background independence; correspondence between the string excitations on $AdS_5 \times S^5$ and those in the LLM geometry.


title: Decision-Making with Auto-Encoding Variational Bayes

abstract: To make decisions based on a model fit with auto-encoding variational Bayes (AEVB), practitioners often let the variational distribution serve as a surrogate for the posterior distribution. This approach yields biased estimates of the expected risk, and therefore leads to poor decisions for two reasons. First, the model fit with AEVB may not equal the underlying data distribution. Second, the variational distribution may not equal the posterior distribution under the fitted model. We explore how fitting the variational distribution based on several objective functions other than the ELBO, while continuing to fit the generative model based on the ELBO, affects the quality of downstream decisions. For the probabilistic principal component analysis model, we investigate how importance sampling error, as well as the bias of the model parameter estimates, varies across several approximate posteriors when used as proposal distributions. Our theoretical results suggest that a posterior approximation distinct from the variational distribution should be used for making decisions. Motivated by these theoretical results, we propose learning several approximate proposals for the best model and combining them using multiple importance sampling for decision-making. In addition to toy examples, we present a full-fledged case study of single-cell RNA sequencing. In this challenging instance of multiple hypothesis testing, our proposed approach surpasses the current state of the art.


title: On Higher-Order Cryptography (Long Version)

abstract: Type-two constructions abound in cryptography: adversaries for encryption and authentication schemes, if active, are modeled as algorithms having access to oracles, i.e. as second-order algorithms. But how about making cryptographic schemes themselves higher-order? This paper gives an answer to this question, by first describing why higher-order cryptography is interesting as an object of study, then showing how the concept of probabilistic polynomial time algorithm can be generalized so as to encompass algorithms of order strictly higher than two, and finally proving some positive and negative results about the existence of higher-order cryptographic primitives, namely authentication schemes and pseudorandom functions.


title: The quasar main sequence and its potential for cosmology

abstract: The main sequence offers a method for the systematization of quasar spectral properties. Extreme FeII emitters (or extreme Population A, xA) are believed to be sources accreting matter at very high rates. They are easily identifiable along the quasar main sequence, in large spectroscopic surveys over a broad redshift range. The very high accretion rate makes it possible that massive black holes hosted in xA quasars radiate at a stable, extreme luminosity-to-mass ratio. After reviewing the basic interpretation of the main sequence, we report on the possibility of identifying virial broadening estimators from low-ionization line widths, and provide evidence of the conceptual validity of redshift-independent luminosities based on virial broadening for a known luminosity-to-mass ratio.


title: Compression stiffening of fibrous networks with stiff inclusions

abstract: Tissues commonly consist of cells embedded within a fibrous biopolymer network. Whereas cell-free reconstituted biopolymer networks typically soften under applied uniaxial compression, various tissues, including liver, brain, and fat, have been observed to instead stiffen when compressed. The mechanism for this compression stiffening effect is not yet clear. Here, we demonstrate that when a material composed of stiff inclusions embedded in a fibrous network is compressed, heterogeneous rearrangement of the inclusions can induce tension within the interstitial network, leading to a macroscopic crossover from an initial bending-dominated softening regime to a stretching-dominated stiffening regime, which occurs before and independently of jamming of the inclusions. Using a coarse-grained particle-network model, we first establish a phase diagram for compression-driven, stretching-dominated stress propagation and jamming in uniaxially compressed 2- and 3-dimensional systems. Then, we demonstrate that a more detailed computational model of stiff inclusions in a subisostatic semiflexible fiber network exhibits quantitative agreement with the predictions of our coarse-grained model as well as qualitative agreement with experiments.


title: Convolutional Support Vector Machine

abstract: The support vector machine (SVM) and deep learning (e.g., convolutional neural networks (CNNs)) are the two most famous algorithms in small and big data, respectively. Nonetheless, smaller datasets may be very important, costly, and not easy to obtain in a short time. This paper proposes a novel convolutional SVM (CSVM) that has both advantages of CNN and SVM to improve the accuracy and effectiveness of mining smaller datasets. The proposed CSVM adapts the convolution product from CNN to learn new information hidden deeply in the datasets. In addition, it uses a modified simplified swarm optimization (SSO) to help train the CSVM to update classifiers, and then the traditional SVM is implemented as the fitness for the SSO to estimate the accuracy. To evaluate the performance of the proposed CSVM, experiments were conducted to test five well-known benchmark databases for the classification problem. Numerical experiments compared favorably with those obtained using SVM, 3-layer artificial NN (ANN), and 4-layer ANN. The results of these experiments verify that the proposed CSVM with the proposed SSO can effectively increase classification accuracy.


title: A path to confine gluons and fermions through complex gauge theory

abstract: We introduced here the study of a QCD based on a complex group. Our aim is to show that a gauge theory with a complex symmetry develops some of the features required for the description of a confined phase. This theory leads to gluons with propagations characteristic of i-particles, fundamental constituents of propagators in the Gribov scenario, together with a gluon propagating with the mode predicted by 't Hooft for the quark confinement. In this way, we present in this same theory a possible gluon condensate, together with a confining interquark potential after the coupling to fermions. This is a novelty, as, up to now, all theories of the GZ kind do not generate quark confinement.


title: Intelligent Methods for Accurately Detecting Phishing Websites

abstract: With increasing technology developments, there is a massive number of websites with varying purposes. But a particular type exists within this large collection, the so-called phishing sites which aim to deceive their users. The main challenge in detecting phishing websites is discovering the techniques that have been used. Where phishers are continually improving their strategies and creating web pages that can protect themselves against many forms of detection methods. Therefore, it is very necessary to develop reliable, active and contemporary methods of phishing detection to combat the adaptive techniques used by phishers. In this paper, different phishing detection approaches are reviewed by classifying them into three main groups. Then, the proposed model is presented in two stages. In the first stage, different machine learning algorithms are applied to validate the chosen dataset and applying features selection methods on it. Thus, the best accuracy was achieved by utilizing only 20 features out of 48 features combined with Random Forest is 98.11%. While in the second stage, the same dataset is applied to various fuzzy logic algorithms. As well the experimental results from the application of Fuzzy logic algorithms were incredible. Where in applying the FURIA algorithm with only five features the accuracy rate was 99.98%. Finally, comparison and discussion of the results between applying machine learning algorithms and fuzzy logic algorithms is done. Where the performance of using fuzzy logic algorithms exceeds the use of machine learning algorithms.