title: Leveraging Planning Landmarks for Hybrid Online Goal Recognition

abstract: Goal recognition is an important problem in many application domains (e.g., pervasive computing, intrusion detection, computer games, etc.). In many application scenarios it is important that goal recognition algorithms can recognize goals of an observed agent as fast as possible and with minimal domain knowledge. Hence, in this paper, we propose a hybrid method for online goal recognition that combines a symbolic planning landmark based approach and a data-driven goal recognition approach and evaluate it in a real-world cooking scenario. The empirical results show that the proposed method is not only significantly more efficient in terms of computation time than the state-of-the-art but also improves goal recognition performance. Furthermore, we show that the utilized planning landmark based approach, which was so far only evaluated on artificial benchmark domains, achieves also good recognition performance when applied to a real-world cooking scenario.


title: The Hubble constant tension: current status and future perspectives through new cosmological probes

abstract: The Hubble constant ($H_0$) tension is one of the major open problems in modern cosmology. This tension is the discrepancy, ranging from 4 to 6 $\sigma$, between the $H_0$ value estimated locally with the combination of Supernovae Ia (SNe Ia) + Cepheids and the cosmological $H_0$ obtained through the study of the Cosmic Microwave Background (CMB) radiation. The approaches adopted in Dainotti et al. 2021 (ApJ) and Dainotti et al. 2022 (Galaxies) are introduced. Through a binning division of the Pantheon sample of SNe Ia (Scolnic et al. 2018), the value of $H_0$ has been estimated in each of the redshift-ordered bins and fitted with a function lowering with the redshift. The results show a decreasing trend of $H_0$ with redshift. If this is not due to astrophysical biases or residual redshift evolution of the SNe Ia parameters, it can be explained in light of modified gravity theories, e.g., the $f(R)$ scenarios. We also briefly describe the possible impact of high-$z$ probes on the Hubble constant tension, such as Gamma-ray bursts (GRBs) and Quasars (QSOs), reported in Dainotti et al. 2022 (Galaxies) and Lenart et al. 2022 (ApJ), respectively.


title: A bridge between convexity and quasiconvexity

abstract: We introduce a notion of convexity with respect to a one-dimensional operator and with this notion find a one-parameter family of different convexities that interpolates between classical convexity and quasiconvexity. We show that, for this interpolation family, the convex envelope of a continuous boundary datum in a strictly convex domain is continuous up to the boundary and is characterized as being the unique viscosity solution to the Dirichlet problem in the domain for a certain fully nonlinear partial differential equation that involves the associated operator. In addition we prove that the convex envelopes of a boundary datum constitute a one-parameter curve of functions that goes from the quasiconvex envelope to the convex envelope being continuous with respect to uniform convergence. Finally, we also show some regularity results for the convex envelopes proving that there is an analogous to a supporting hyperplane at every point and that convex envelopes are $C^1$ if the boundary data satisfies in particular $NV$-condition we introduce.


title: DIFFER: Decomposing Individual Reward for Fair Experience Replay in Multi-Agent Reinforcement Learning

abstract: Cooperative multi-agent reinforcement learning (MARL) is a challenging task, as agents must learn complex and diverse individual strategies from a shared team reward. However, existing methods struggle to distinguish and exploit important individual experiences, as they lack an effective way to decompose the team reward into individual rewards. To address this challenge, we propose DIFFER, a powerful theoretical framework for decomposing individual rewards to enable fair experience replay in MARL. By enforcing the invariance of network gradients, we establish a partial differential equation whose solution yields the underlying individual reward function. The individual TD-error can then be computed from the solved closed-form individual rewards, indicating the importance of each piece of experience in the learning task and guiding the training process. Our method elegantly achieves an equivalence to the original learning framework when individual experiences are homogeneous, while also adapting to achieve more muscular efficiency and fairness when diversity is observed.Our extensive experiments on popular benchmarks validate the effectiveness of our theory and method, demonstrating significant improvements in learning efficiency and fairness.


title: Trainable Loss Weights in Super-Resolution

abstract: In recent years, limited research has discussed the loss function in the super-resolution process. The majority of those studies have only used perceptual similarity conventionally. This is while the development of appropriate loss can improve the quality of other methods as well. In this article, a new weighting method for pixel-wise loss is proposed. With the help of this method, it is possible to use trainable weights based on the general structure of the image and its perceptual features while maintaining the advantages of pixel-wise loss. Also, a criterion for comparing weights of loss is introduced so that the weights can be estimated directly by a convolutional neural network. In addition, in this article, the expectation-maximization method is used for the simultaneous estimation super-resolution network and weighting network. In addition, a new activation function, called "FixedSum", is introduced which can keep the sum of all components of vector constants while keeping the output components between zero and one. As experimental results shows, weighted loss by the proposed method leads to better results than the unweighted loss and weighted loss based on uncertainty in both signal-to-noise and perceptual similarity senses on the state-of-the-art networks. Code is available online.


title: A Study on FGSM Adversarial Training for Neural Retrieval

abstract: Neural retrieval models have acquired significant effectiveness gains over the last few years compared to term-based methods. Nevertheless, those models may be brittle when faced to typos, distribution shifts or vulnerable to malicious attacks. For instance, several recent papers demonstrated that such variations severely impacted models performances, and then tried to train more resilient models. Usual approaches include synonyms replacements or typos injections -- as data-augmentation -- and the use of more robust tokenizers (characterBERT, BPE-dropout). To further complement the literature, we investigate in this paper adversarial training as another possible solution to this robustness issue. Our comparison includes the two main families of BERT-based neural retrievers, i.e. dense and sparse, with and without distillation techniques. We then demonstrate that one of the most simple adversarial training techniques -- the Fast Gradient Sign Method (FGSM) -- can improve first stage rankers robustness and effectiveness. In particular, FGSM increases models performances on both in-domain and out-of-domain distributions, and also on queries with typos, for multiple neural retrievers.


title: ARDIAS: AI-Enhanced Research Management, Discovery, and Advisory System

abstract: In this work, we present ARDIAS, a web-based application that aims to provide researchers with a full suite of discovery and collaboration tools. ARDIAS currently allows searching for authors and articles by name and gaining insights into the research topics of a particular researcher. With the aid of AI-based tools, ARDIAS aims to recommend potential collaborators and topics to researchers. In the near future, we aim to add tools that allow researchers to communicate with each other and start new projects.


title: Strongly proper connected coloring of graphs

abstract: We study a new variant of \emph{connected coloring} of graphs based on the concept of \emph{strong} edge coloring (every color class forms an \emph{induced} matching). In particular, an edge-colored path is \emph{strongly proper} if its color sequence does not contain identical terms within a distance of at most two. A \emph{strong proper connected} coloring of $G$ is the one in which every pair of vertices is joined by at least one strongly proper path. Let spc($G$) denote the least number of colors needed for such coloring of a graph $G$. We prove that the upper bound spc($G$)$\leq${5} holds for any $2$-connected graph $G$. On the other hand, we demonstrate that there are $2$-connected graphs with arbitrarily large girth satisfying spc($G$)$\geq${4}. Additionally, we prove that graphs whose cycle lengths are divisible by $3$ satisfy spc($G$)$\leq{3}$. We also consider briefly other connected colorings defined by various restrictions on color sequences of connecting paths. For instance, in a \emph{nonrepetitive connected coloring} of $G$, every pair of vertices should be joined by a path whose color sequence is \emph{nonrepetitive}, that is, it does not contain two adjacent identical blocks. We demonstrate that $2$-connected graphs are $15$-colorable while $4$-connected graphs are $6$-colorable, in the connected nonrepetitive sense. A similar conclusion with a finite upper bound on the number of colors holds for a much wider variety of connected colorings corresponding to fairly general properties of sequences. We end the paper with some open problems of concrete and general nature.


title: Fokker-Planck equation for the particle size distribution function in KJMA transformations

abstract: The Fokker-Planck (FP) equation has been derived for describing the temporal evolution of the particle size probability density function (PDF) for KJMA (Kolmogorov-Johnson-Mehl-Avrami) transformations. The classical case of transformations with constant rates of both nucleation and growth, in 3D space, has been considered. Integration of the equation shows that the PDF is given by the superposition of one-parameter Gamma distributions with time dependent mean size given by the KJMA theory. The asymptotic behavior of the FP solution offers a demonstration of the conjecture, previously proposed by Pineda et al [E. Pineda, P. Bruna, D. Crespo, Phys. Rev. E 70 (2004) 066119], according to which the set of nuclei formed at the same time are Gamma-distributed, with parameter depending on nucleus birth time. Computer simulations of the transformation with constant nucleation and growth rates, do show that the temporal evolution of the PDF, in the volume domain, is in good agreement with the Johnson-Mehl PDF. The approach based on the FP equations provides a particle size PDF that exhibits such behavior.


title: A Mathematical Programming approach to Overlapping community detection by cluster's Shapley value

abstract: We propose a new model to detect the overlapping communities of a network that is based on cooperative games and mathematical programming. More specifically, communities are defined as stable coalitions of a weighted graph community game and they are revealed as the optimal solution of a mixed-integer linear programming problem. Exact optimal solutions are obtained for small and medium sized instances and and it is shown that they provide useful information about the network structure, improving on previous contributions. Next, a heuristic algorithm is developed to solve the largest instances and used to compare two variations of the objective function.


title: A Naive Bayes Classifier for identifying Class II YSOs

abstract: A naive Bayes classifier for identifying Class II YSOs has been constructed and applied to a region of the Northern Galactic Plane containing 8 million sources with good quality Gaia EDR3 parallaxes. The classifier uses the five features: Gaia $G$-band variability, WISE mid-infrared excess, UKIDSS and 2MASS near-infrared excess, IGAPS H$\alpha$ excess and overluminosity with respect to the main sequence. A list of candidate Class II YSOs is obtained by choosing a posterior threshold appropriate to the task at hand, balancing the competing demands of completeness and purity. At a threshold posterior greater than 0.5 our classifier identifies 6504 candidate Class II YSOs. At this threshold we find a false positive rate around 0.02 per cent and a true positive rate of approximately 87 per cent for identifying Class II YSOs. The ROC curve rises rapidly to almost one with an area under the curve around 0.998 or better, indicating the classifier is efficient at identifying candidate Class II YSOs. Our map of these candidates shows what are potentially three previously undiscovered clusters or associations. When comparing our results to published catalogues from other young star classifiers, we find between one quarter and three quarters of high probability candidates are unique to each classifier, telling us no single classifier is finding all young stars.


title: Mechanism to induce geometric constriction on kinks and domain walls

abstract: We investigate scalar field theories in the multifield scenario, focusing mainly on the possibility to smoothly build internal structure and asymmetry for kinks and domain walls. The procedure requires the inclusion of an extra field which is associated to a function that modifies the dynamics of the other fields. We investigate minimum energy configurations, which support first order equations compatible with the equations of motion. The extra field allows a transition which is guided by a parameter that connects the standard solution to another one, geometrically constrained, mimicking the effects of geometrical constrictions in magnetic materials.


title: An Efficient Approximate Method for Online Convolutional Dictionary Learning

abstract: Most existing convolutional dictionary learning (CDL) algorithms are based on batch learning, where the dictionary filters and the convolutional sparse representations are optimized in an alternating manner using a training dataset. When large training datasets are used, batch CDL algorithms become prohibitively memory-intensive. An online-learning technique is used to reduce the memory requirements of CDL by optimizing the dictionary incrementally after finding the sparse representations of each training sample. Nevertheless, learning large dictionaries using the existing online CDL (OCDL) algorithms remains highly computationally expensive. In this paper, we present a novel approximate OCDL method that incorporates sparse decomposition of the training samples. The resulting optimization problems are addressed using the alternating direction method of multipliers. Extensive experimental evaluations using several image datasets show that the proposed method substantially reduces computational costs while preserving the effectiveness of the state-of-the-art OCDL algorithms.


title: A Method For Eliminating Contour Errors In Self-Encoder Reconstructed Images

abstract: In this paper, we propose a self-supervised twin network approach based on this a priori. The method of generating the approximate10 edge information of an image and then differentially eliminating the edge errors11 in the reconstructed image with a dilate algorithm. This is used to improve the12 accuracy of the reconstructed image and to separate foreign matter and noise from13 the original image, so that it can be visualized in a more practical scene


title: Evaluation of the syllables pronunciation quality in speech rehabilitation through the solution of the classification problem

abstract: The solution of the problem of assessing the quality of the pronunciation of syllables during speech rehabilitation after surgical treatment of oncological diseases of the organs of the speech-forming tract is considered in the work. The assessment is carried out by solving the problem of classifying syllables into two classes: before and immediately after surgical treatment. A classifier is built on the basis of the LSTM neural network and trained on the records before the operation and immediately after it, before the start of speech rehabilitation. The measure of assessing the quality of syllables pronunciation in the process of rehabilitation is the metric of belonging to the class before the operation. A study is being made of the influence of taking into account problematic phonemes, the gender of the patient, his individual characteristics on the resulting estimates of the quality of pronunciation. A comparison with existing types of syllable pronunciation quality assessments is carried out, recommendations are given for the practical application of the resulting new class of pronunciation quality assessments.


title: Transition from wave turbulence to acousticlike shock-wave regime

abstract: We report on the experimental observation of a transition from a dispersive wave turbulence regime to a nondispersive regime involving shock waves on the surface of a fluid. We use a magnetic fluid in a canal subjected to an external horizontal magnetic field to tune the dispersivity of the system. For a low magnetic field, gravity-capillary wave turbulence is observed, whereas for a high enough field, random steep coherent structures arise which are found to be shock waves. These shock waves create singularities in the second-order difference of the surface elevation, leading to an $\omega^{-4}$ frequency power spectrum. This spectrum is also found to be controlled by the number and amplitude of the shocks and is well captured by a model based on a random Dirac-$\delta$ distribution (Kuznetsov-like spectrum). Finally, the shock-amplitude statistics exhibits a power-law distribution with an exponent close to the predictions of the one-dimensional random-forced Burgers equation. This shock-wave regime, discovered here for surface waves, thus paves the way to better explore their properties.


title: On Batching Variable Size Inputs for Training End-to-End Speech Enhancement Systems

abstract: The performance of neural network-based speech enhancement systems is primarily influenced by the model architecture, whereas training times and computational resource utilization are primarily affected by training parameters such as the batch size. Since noisy and reverberant speech mixtures can have different duration, a batching strategy is required to handle variable size inputs during training, in particular for state-of-the-art end-to-end systems. Such strategies usually strive for a compromise between zero-padding and data randomization, and can be combined with a dynamic batch size for a more consistent amount of data in each batch. However, the effect of these strategies on resource utilization and more importantly network performance is not well documented. This paper systematically investigates the effect of different batching strategies and batch sizes on the training statistics and speech enhancement performance of a Conv-TasNet, evaluated in both matched and mismatched conditions. We find that using a small batch size during training improves performance in both conditions for all batching strategies. Moreover, using sorted or bucket batching with a dynamic batch size allows for reduced training time and GPU memory usage while achieving similar performance compared to random batching with a fixed batch size.


title: A DPG method for the quad-curl problem

abstract: We derive an ultraweak variational formulation of the quad-curl problem in two and three dimensions. We present a discontinuous Petrov-Galerkin (DPG) method for its approximation and prove its quasi-optimal convergence. We illustrate how this method can be applied to the Stokes problem in two dimensions, after an application of the curl operator to eliminate the pressure variable. In this way, DPG techniques known from Kirchhoff-Love plates can be used. We present an a priori error estimate that improves a previous approximation result for effective shear forces by using a less restrictive regularity assumption. Numerical experiments illustrate our findings.


title: Hubbard bands, Mott transition and deconfinement in strongly correlated systems

abstract: The problem of deconfinement phases in strongly correlated systems is discussed. In space-time dimension $d=3+1$, a competition of confinement and Coulomb phases occurs, but in $d=2+1$ the confining phase dominates owing to monopole proliferation, but gapless fermion excitations can change the situation. Combining the Kotliar-Ruckenstein representation and fractionalized spin-liquid deconfinement picture, the Mott transition and Hubbard subbands are treated, general expressions in the case of an arbitrary bare band spectrum being obtained. The transition into a metallic state is determined by condensation of a gapless boson mode. The spectrum picture in the insulating state is considerably influenced by the spinon spin-liquid spectrum and hidden Fermi surface.


title: Quantum optics measurement scheme for quantum geometry and topological invariants

abstract: We show how a quantum optical measurement scheme based on heterodyne detection can be used to explore geometrical and topological properties of condensed matter systems. Considering a 2D material placed in a cavity with a coupling to the environment, we compute correlation functions of the photons exiting the cavity and relate them to the hybrid light-matter state within the cavity. Different polarizations of the intracavity field give access to all components of the quantum geometric tensor on contours in the Brillouin zone defined by the transition energy. Combining recent results based on the metric-curvature correspondence with the measured quantum metric allows us to characterize the topological phase of the material. Moreover, in systems where $S_z$ is a good quantum number, the procedure also allows us to extract the spin Chern number. As an interesting application, we consider a minimal model for twisted bilayer graphene at the magic angle, and discuss the feasibility of extracting the Euler number.


title: Interstellar Extinction in Galactic Cirri in SDSS Stripe 82

abstract: We have applied the method of star counts with Wolf diagrams to determine the interstellar extinction in five Galactic cirri in Sloan Digital Sky Survey (SDSS) Stripe 82. For this purpose, we have used the photometry of stars in the GALEX NUV filter and the photometry of red dwarfs in five SDSS bands and four SkyMapper Southern Sky Survey DR2 bands. We have identified the cirri as sky regions with an enhanced infrared emission from the Schlegel+1998 map. The extinction in them has been calculated relative to the nearby comparison regions with a reduced emission. The results for different filters agree well, giving the range of distances and the extinction law for each cirrus. The distances in the range 140--415 pc found are consistent with the 3D reddening maps. In the range between the $B$ and $V$ filters the extinctions found are consistent with the estimates from Schlegel+1998 for the Cardelli+1989 extinction law with $R_\mathrm{V}=3.1$. However, the extinctions found for all of the filters are best described not by the Cardelli+1989 extinction law with some $R_\mathrm{V}=3.1$, but by the inverse proportionality of the extinction and wavelength with its own coefficient for each cirrus. In one of the cirri our results suggest a very slight decrease in extinction with wavelength, i.e., a large contribution of gray extinction. In the remaining cirri a manifestation of gray extinction is not ruled out either. This is consistent with the previous measurements of the extinction law far from the Galactic midplane.


title: Hierarchical Regularizers for Reverse Unrestricted Mixed Data Sampling Regressions

abstract: Reverse Unrestricted MIxed DAta Sampling (RU-MIDAS) regressions are used to model high-frequency responses by means of low-frequency variables. However, due to the periodic structure of RU-MIDAS regressions, the dimensionality grows quickly if the frequency mismatch between the high- and low-frequency variables is large. Additionally the number of high-frequency observations available for estimation decreases. We propose to counteract this reduction in sample size by pooling the high-frequency coefficients and further reduce the dimensionality through a sparsity-inducing convex regularizer that accounts for the temporal ordering among the different lags. To this end, the regularizer prioritizes the inclusion of lagged coefficients according to the recency of the information they contain. We demonstrate the proposed method on an empirical application for daily realized volatility forecasting where we explore whether modeling high-frequency volatility data in terms of low-frequency macroeconomic data pays off.


title: Faster DAN: Multi-target Queries with Document Positional Encoding for End-to-end Handwritten Document Recognition

abstract: Recent advances in handwritten text recognition enabled to recognize whole documents in an end-to-end way: the Document Attention Network (DAN) recognizes the characters one after the other through an attention-based prediction process until reaching the end of the document. However, this autoregressive process leads to inference that cannot benefit from any parallelization optimization. In this paper, we propose Faster DAN, a two-step strategy to speed up the recognition process at prediction time: the model predicts the first character of each text line in the document, and then completes all the text lines in parallel through multi-target queries and a specific document positional encoding scheme. Faster DAN reaches competitive results compared to standard DAN, while being at least 4 times faster on whole single-page and double-page images of the RIMES 2009, READ 2016 and MAURDOR datasets. Source code and trained model weights are available at https://github.com/FactoDeepLearning/FasterDAN.


title: A Discussion on Nonlinear Quadratic Control and Sontag's Formula

abstract: The quadratic optimal state feedback (LQR) is one of the most popular designs for linear systems and succeeds via the solution of the algebraic Riccati equation. The situation is different in the case of non-linear systems: the Riccati equation is then replaced by the Hamilton Jacobi Bellman equation (HJB), the solution of which is generally difficult. A compromise can be the so-called Inverse Optimal Control, a form of which is Sontag's formula [1]; here the minimized cost function follows from the feedback law chosen, not vice versa. Using Sontag's formula in the variant according to Freeman and Primbs [2, 9], the actually minimized cost function is given in the following sections, including cases when it reduces to the quadratic cost. Also some remarks and thoughts are presented for discussion.


title: Distilling Text into Circuits

abstract: This paper concerns the structure of meanings within natural language. Earlier, a framework named DisCoCirc was sketched that (1) is compositional and distributional (a.k.a. vectorial); (2) applies to general text; (3) captures linguistic `connections' between meanings (cf. grammar) (4) updates word meanings as text progresses; (5) structures sentence types; (6) accommodates ambiguity. Here, we realise DisCoCirc for a substantial fragment of English. When passing to DisCoCirc's text circuits, some `grammatical bureaucracy' is eliminated, that is, DisCoCirc displays a significant degree of (7) inter- and intra-language independence. That is, e.g., independence from word-order conventions that differ across languages, and independence from choices like many short sentences vs. few long sentences. This inter-language independence means our text circuits should carry over to other languages, unlike the language-specific typings of categorial grammars. Hence, text circuits are a lean structure for the `actual substance of text', that is, the inner-workings of meanings within text across several layers of expressiveness (cf. words, sentences, text), and may capture that what is truly universal beneath grammar. The elimination of grammatical bureaucracy also explains why DisCoCirc: (8) applies beyond language, e.g. to spatial, visual and other cognitive modes. While humans could not verbally communicate in terms of text circuits, machines can. We first define a `hybrid grammar' for a fragment of English, i.e. a purpose-built, minimal grammatical formalism needed to obtain text circuits. We then detail a translation process such that all text generated by this grammar yields a text circuit. Conversely, for any text circuit obtained by freely composing the generators, there exists a text (with hybrid grammar) that gives rise to it. Hence: (9) text circuits are generative for text.


title: Vanishing Theorem for Hodge Ideals on smooth hypersurfaces

abstract: We use a Koszul-type resolution to prove a weak version of Bott's vanishing theorem for smooth hypersurfaces in $\mathbb{P}^n$ and use this result to prove a vanishing theorem for Hodge ideals associated with an effective Cartier divisor on a hypersurface. This extends an earlier result of Mustata and Popa.


title: Breaking Bad: Quantifying the Addiction of Web Elements to JavaScript

abstract: While JavaScript established itself as a cornerstone of the modern web, it also constitutes a major tracking and security vector, thus raising critical privacy and security concerns. In this context, some browser extensions propose to systematically block scripts reported by crowdsourced trackers lists. However, this solution heavily depends on the quality of these built-in lists, which may be deprecated or incomplete, thus exposing the visitor to unknown trackers. In this paper, we explore a different strategy, by investigating the benefits of disabling JavaScript in the browser. More specifically, by adopting such a strict policy, we aim to quantify the JavaScript addiction of web elements composing a web page, through the observation of web breakages. As there is no standard mechanism for detecting such breakages, we introduce a framework to inspect several page features when blocking JavaScript, that we deploy to analyze 6,384 pages, including landing and internal web pages. We discover that 43% of web pages are not strictly dependent on JavaScript and that more than 67% of pages are likely to be usable as long as the visitor only requires the content from the main section of the page, for which the user most likely reached the page, while reducing the number of tracking requests by 85% on average. Finally, we discuss the viability of currently browsing the web without JavaScript and detail multiple incentives for websites to be kept usable without JavaScript.


title: A note on Hamiltonian stability for sheaves

abstract: We show a strong Hamiltonian stability result for a simpler and larger distance on the Tamarkin category. We also give a stability result with support conditions.


title: AnisoTag: 3D Printed Tag on 2D Surface via Reflection Anisotropy

abstract: In the past few years, the widespread use of 3D printing technology enables the growth of the market of 3D printed products. On Esty, a website focused on handmade items, hundreds of individual entrepreneurs are selling their 3D printed products. Inspired by the positive effects of machine-readable tags, like barcodes, on daily product marketing, we propose AnisoTag, a novel tagging method to encode data on the 2D surface of 3D printed objects based on reflection anisotropy. AnisoTag has an unobtrusive appearance and much lower extraction computational complexity, contributing to a lightweight low-cost tagging system for individual entrepreneurs. On AnisoTag, data are encoded by the proposed tool as reflective anisotropic microstructures, which would reflect distinct illumination patterns when irradiating by collimated laser. Based on it, we implement a real-time detection prototype with inexpensive hardware to determine the reflected illumination pattern and decode data according to their mapping. We evaluate AnisoTag with various 3D printer brands, filaments, and printing parameters, demonstrating its superior usability, accessibility, and reliability for practical usage.


title: Efficiency in local differential privacy

abstract: We develop a theory of asymptotic efficiency in regular parametric models when data confidentiality is ensured by local differential privacy (LDP). Even though efficient parameter estimation is a classical and well-studied problem in mathematical statistics, it leads to several non-trivial obstacles that need to be tackled when dealing with the LDP case. Starting from a standard parametric model $\mathcal P=(P_\theta)_{\theta\in\Theta}$, $\Theta\subseteq\mathbb R^p$, for the iid unobserved sensitive data $X_1,\dots, X_n$, we establish local asymptotic mixed normality (along subsequences) of the model $$Q^{(n)}\mathcal P=(Q^{(n)}P_\theta^n)_{\theta\in\Theta}$$ generating the sanitized observations $Z_1,\dots, Z_n$, where $Q^{(n)}$ is an arbitrary sequence of sequentially interactive privacy mechanisms. This result readily implies convolution and local asymptotic minimax theorems. In case $p=1$, the optimal asymptotic variance is found to be the inverse of the supremal Fisher-Information $\sup_{Q\in\mathcal Q_\alpha} I_\theta(Q\mathcal P)\in\mathbb R$, where the supremum runs over all $\alpha$-differentially private (marginal) Markov kernels. We present an algorithm for finding a (nearly) optimal privacy mechanism $\hat{Q}$ and an estimator $\hat{\theta}_n(Z_1,\dots, Z_n)$ based on the corresponding sanitized data that achieves this asymptotically optimal variance.


title: Rapid Reconstruction of 3-D Membrane Pore Structure Using a Single 2-D Micrograph

abstract: Conventional 2-D scanning electron microscopy (SEM) is commonly used to rapidly and qualitatively evaluate membrane pore structure. Quantitative 2-D analyses of pore sizes can be extracted from SEM, but without information about 3-D spatial arrangement and connectivity, which are crucial to the understanding of membrane pore structure. Meanwhile, experimental 3-D reconstruction via tomography is complex, expensive, and not easily accessible. Here, we employ data-science tools to demonstrate a proof-of-principle reconstruction of the 3-D structure of a membrane using a single 2-D image pulled from a 3-D tomographic data set. The reconstructed and experimental 3-D structures were then directly compared, with important properties such as mean pore radius, mean throat radius, coordination number and tortuosity differing by less than 15%. The developed algorithm will dramatically improve the ability of the membrane community to characterize membranes, accelerating the design and synthesis of membranes with desired structural and transport properties.


title: DreamWaQ: Learning Robust Quadrupedal Locomotion With Implicit Terrain Imagination via Deep Reinforcement Learning

abstract: Quadrupedal robots resemble the physical ability of legged animals to walk through unstructured terrains. However, designing a controller for quadrupedal robots poses a significant challenge due to their functional complexity and requires adaptation to various terrains. Recently, deep reinforcement learning, inspired by how legged animals learn to walk from their experiences, has been utilized to synthesize natural quadrupedal locomotion. However, state-of-the-art methods strongly depend on a complex and reliable sensing framework. Furthermore, prior works that rely only on proprioception have shown a limited demonstration for overcoming challenging terrains, especially for a long distance. This work proposes a novel quadrupedal locomotion learning framework that allows quadrupedal robots to walk through challenging terrains, even with limited sensing modalities. The proposed framework was validated in real-world outdoor environments with varying conditions within a single run for a long distance.


title: Stochastic particle transport by deep-water irregular breaking waves

abstract: Correct prediction of particle transport by surface waves is crucial in many practical applications such as search and rescue or salvage operations and pollution tracking and clean-up efforts. Recent results have indicated transport by deep-water breaking waves is enhanced compared to non-breaking waves. To model particle transport in irregular waves, some of which break, we develop a stochastic differential equation describing both mean particle transport and its uncertainty. The equation combines a Brownian motion, which captures non-breaking drift-diffusion effects, and a compound Poisson process, which captures jumps in particle positions due to breaking.We corroborate these predictions with new experiments, in which we track large numbers of particles in irregular breaking waves. For breaking and non-breaking wave fields, our experiments confirm that the variance of the particle position grows linearly with time, in accordance with Taylor's single-particle dispersion theory. For wave fields that include breaking, the compound Poisson process increases the linear growth rate of the mean and variance and introduces a finite skewness of the particle position distribution.


title: Automated multilingual detection of Pro-Kremlin propaganda in newspapers and Telegram posts

abstract: The full-scale conflict between the Russian Federation and Ukraine generated an unprecedented amount of news articles and social media data reflecting opposing ideologies and narratives. These polarized campaigns have led to mutual accusations of misinformation and fake news, shaping an atmosphere of confusion and mistrust for readers worldwide. This study analyses how the media affected and mirrored public opinion during the first month of the war using news articles and Telegram news channels in Ukrainian, Russian, Romanian and English. We propose and compare two methods of multilingual automated pro-Kremlin propaganda identification, based on Transformers and linguistic features. We analyse the advantages and disadvantages of both methods, their adaptability to new genres and languages, and ethical considerations of their usage for content moderation. With this work, we aim to lay the foundation for further development of moderation tools tailored to the current conflict.


title: A note on Hausdorff-Young inequalities in function spaces

abstract: The classical Hausdorff-Young inequalities for the Fourier transform acting between appropriate $L_p$ spaces are cornerstones of Fourier analysis. Here we extend it to weighted spaces of Besov or Sobolev type where the weight has the form $w(x)=(1+|x|^2)^{\alpha/2}$. This note is not a paper or draft but a sketchy complement to some earlier results where we dealt with mapping properties of the Fourier transform.


title: A Holistic Cascade System, benchmark, and Human Evaluation Protocol for Expressive Speech-to-Speech Translation

abstract: Expressive speech-to-speech translation (S2ST) aims to transfer prosodic attributes of source speech to target speech while maintaining translation accuracy. Existing research in expressive S2ST is limited, typically focusing on a single expressivity aspect at a time. Likewise, this research area lacks standard evaluation protocols and well-curated benchmark datasets. In this work, we propose a holistic cascade system for expressive S2ST, combining multiple prosody transfer techniques previously considered only in isolation. We curate a benchmark expressivity test set in the TV series domain and explored a second dataset in the audiobook domain. Finally, we present a human evaluation protocol to assess multiple expressive dimensions across speech pairs. Experimental results indicate that bi-lingual annotators can assess the quality of expressive preservation in S2ST systems, and the holistic modeling approach outperforms single-aspect systems. Audio samples can be accessed through our demo webpage: https://facebookresearch.github.io/speech_translation/cascade_expressive_s2st.


title: The quark propagator and quark-gluon vertex from lattice QCD at finite temperature

abstract: The quark-gluon vertex is an important object of QCD. Studies have shown that this quantity is relevant for the dynamical chiral symmetry breaking pattern in the vacuum. The goal of our project is to obtain the quark-gluon vertex at finite temperature around the deconfinement/chiral transition using the tools provided by lattice QCD. It will be the first time that the quark-gluon vertex at finite temperature is determined using lattice QCD. The propagators, which are a by-product of this project, are also of interest in themselves. The configurations used were generated by the FASTSUM collaboration. In this contribution, we describe our motivations and goals, some technical details of the determination and report on the status of the calculation.


title: Connecting metrics for shape-texture knowledge in computer vision

abstract: Modern artificial neural networks, including convolutional neural networks and vision transformers, have mastered several computer vision tasks, including object recognition. However, there are many significant differences between the behavior and robustness of these systems and of the human visual system. Deep neural networks remain brittle and susceptible to many changes in the image that do not cause humans to misclassify images. Part of this different behavior may be explained by the type of features humans and deep neural networks use in vision tasks. Humans tend to classify objects according to their shape while deep neural networks seem to rely mostly on texture. Exploring this question is relevant, since it may lead to better performing neural network architectures and to a better understanding of the workings of the vision system of primates. In this work, we advance the state of the art in our understanding of this phenomenon, by extending previous analyses to a much larger set of deep neural network architectures. We found that the performance of models in image classification tasks is highly correlated with their shape bias measured at the output and penultimate layer. Furthermore, our results showed that the number of neurons that represent shape and texture are strongly anti-correlated, thus providing evidence that there is competition between these two types of features. Finally, we observed that while in general there is a correlation between performance and shape bias, there are significant variations between architecture families.


title: Phase diagram of the Ashkin-Teller model

abstract: The Ashkin-Teller model is a pair of interacting Ising models and has two parameters: $J$ is a coupling constant in the Ising models and $U$ describes the strength of the interaction between them. In the ferromagnetic case $J,U>0$ on the square lattice, we establish a complete phase diagram conjectured in physics in 1970s (by Kadanoff and Wegner, Wu and Lin, Baxter and others): when $J<U$, the transitions for the Ising spins and their products occur at two distinct curves that are dual to each other; when $J\geq U$, both transitions occur at the self-dual curve. All transitions are shown to be sharp using the OSSS inequality. We use a finite-criterion argument and continuity to extend the result of Peled and the third author \cite{GlaPel19} from a self-dual point to its neighborhood. Our proofs go through the random-cluster representation of the Ashkin-Teller model introduced by Chayes-Machta and Pfister-Velenik and we rely on couplings to FK-percolation.


title: Forty Thousand Kilometers Under Quantum Protection

abstract: Quantum key distribution (QKD) is a revolutionary cryptography response to the rapidly growing cyberattacks threat posed by quantum computing. Yet, the roadblock limiting the vast expanse of secure quantum communication is the exponential decay of the transmitted quantum signal with the distance. Today's quantum cryptography is trying to solve this problem by focusing on quantum repeaters. However, efficient and secure quantum repetition at sufficient distances is still far beyond modern technology. Here, we shift the paradigm and build the long-distance security of the QKD upon the quantum foundations of the Second Law of Thermodynamics and end-to-end physical oversight over the transmitted optical quantum states. Our approach enables us to realize quantum states' repetition by optical amplifiers keeping states' wave properties and phase coherence. The unprecedented secure distance range attainable through our approach opens the door for the development of scalable quantum-resistant communication networks of the future.


title: Discriminator-free Unsupervised Domain Adaptation for Multi-label Image Classification

abstract: In this paper, a discriminator-free adversarial-based Unsupervised Domain Adaptation (UDA) for Multi-Label Image Classification (MLIC) referred to as DDA-MLIC is proposed. Recently, some attempts have been made for introducing adversarial-based UDA methods in the context of MLIC. However, these methods which rely on an additional discriminator subnet present one major shortcoming. The learning of domain-invariant features may harm their task-specific discriminative power, since the classification and discrimination tasks are decoupled. Herein, we propose to overcome this issue by introducing a novel adversarial critic that is directly deduced from the task-specific classifier. Specifically, a two-component Gaussian Mixture Model (GMM) is fitted on the source and target predictions in order to distinguish between two clusters. This allows extracting a Gaussian distribution for each component. The resulting Gaussian distributions are then used for formulating an adversarial loss based on a Frechet distance. The proposed method is evaluated on several multi-label image datasets covering three different types of domain shift. The obtained results demonstrate that DDA-MLIC outperforms existing state-of-the-art methods in terms of precision while requiring a lower number of parameters. The code is publicly available at github.com/cvi2snt/DDA-MLIC.


title: Consensus Algorithm For Calculation Of Protein Binding Affinity Using Multiple Models

abstract: The major histocompatibility complex (MHC) molecules, which bind peptides for presentation on the cell surface, play an important role in cell-mediated immunity. In light of developing databases and technologies over the years, significant progress has been made in research on peptide binding affinity calculation. Several in techniques have been developed to predict peptide binding to MHC class I. Most of the research on MHC Class I due to its nature brings better performance and more. Considering the use of different methods and different technologies, and the approach of similar methods on different proteins, a classification was created according to the binding affinity of protein peptides.


title: Hubble Space Telescope imaging of the compact elliptical galaxy M32 reveals a dearth of carbon stars

abstract: We present new Hubble Space Telescope WFC3/IR medium-band photometry of the compact elliptical galaxy M32, chemically resolving its thermally pulsating asymptotic giant branch stars. We find 2829 M-type stars and 57 C stars. The carbon stars are likely contaminants from M31. If carbon stars are present in M32 they are so in very low numbers. The uncorrected C/M ratio is 0.020 $\pm$ 0.003; this drops to less than 0.007 after taking into account contamination from M31. As the mean metallicity of M32 is just below solar, this low ratio of C to M stars is unlikely due to a metallicity ceiling for the formation of carbon stars. Instead, the age of the AGB population is likely to be the primary factor. The ratio of AGB to RGB stars in M32 is similar to that of the inner disc of M31 which contain stars that formed 1.5-4 Gyr ago. If the M32 population is at the older end of this age then its lack of C-stars may be consistent with a narrow mass range for carbon star formation predicted by some stellar evolution models. Applying our chemical classifications to the dusty variable stars identified with {\em Spitzer}, we find that the x-AGB candidates identified with Spitzer are predominately M-type stars. This substantially increases the lower limit to the cumulative dust-production rate in M32 to $>$ 1.20 $\times 10^{-5}$ ${\rm M}_{\odot} \, {\rm yr}^{-1}$.


title: Optical and near-infrared stellar activity characterization of the early M dwarf Gl~205 with SOPHIE and SPIRou

abstract: The stellar activity of M dwarfs is the main limitation for discovering and characterizing exoplanets orbiting them since it induces quasi-periodic RV variations. We aim to characterize the magnetic field and stellar activity of the early, moderately active, M dwarf Gl205 in the optical and nIR domains. We obtained high-precision quasi-simultaneous spectra in the optical and nIR with the SOPHIE spectrograph and SPIRou spectropolarimeter between 2019 and 2022. We computed the RVs from both instruments and the SPIRou Stokes V profiles. We used ZDI to map the large-scale magnetic field over the time span of the observations. We studied the temporal behavior of optical and nIR RVs and activity indicators with the Lomb-Scargle periodogram and a quasi-periodic GP regression. In the nIR, we studied the equivalent width of Al I, Ti I, K I, Fe I, and He I. We modeled the activity-induced RV jitter using a multi-dimensional GP regression with activity indicators as ancillary time series. The optical and nIR RVs have similar scatter but nIR shows a more complex temporal evolution. We observe an evolution of the magnetic field topology from a poloidal dipolar field in 2019 to a dominantly toroidal field in 2022. We measured a stellar rotation period of Prot=34.4$\pm$0.5 d in the longitudinal magnetic field. Using ZDI we measure the amount of latitudinal differential rotation (DR) shearing the stellar surface yielding rotation periods of Peq=32.0$\pm$1.8 d at the stellar equator and Ppol=45.5$\pm$0.3 d at the poles. We observed inconsistencies in the activity indicators' periodicities that could be explained by these DR values. The multi-dimensional GP modeling yields an RMS of the RV residuals down to the noise level of 3 m/s for both instruments, using as ancillary time series H$\alpha$ and the BIS in the optical, and the FWHM in the nIR.


title: Quakes of Compact Stars

abstract: Glitches are commonly observed for pulsars, which are explained by various mechanisms. One hypothesis attributes the glitch effect to the instantaneous moment of inertia change of the whole star caused by a starquake, which is similar to earthquakes caused by fast dislocation occurring on planar faults for the static stress, though the quake-induced dynamics responsible for glitch (superfluid vortex vs. pure starquake) remains still unknown. However, a theoretical model to quantitatively explain the stress loading, types of starquakes, and co-seismic change of moment of inertia is rarely discussed. In this study, we incorporate elastic deformation theories of earthquakes into the starquake problems. We compute the field of stress loading associated with rotation deceleration and determine the optimal type of starquakes at various locations. Two types of pulsar structure models, i.e. neutron and strangeon star models, are included in the computation and their differences are notable. Our calculation shows that the observed glitch amplitude can be explained by the starquakes in the strangeon star model, though the required scaled starquake magnitude is much larger than that occurred on the Earth. We further discuss the possibility to compute the energy budget and other glitch phenomena using the starquake model in the elastic medium framework.


title: Prediction of COVID-19 by Its Variants using Multivariate Data-driven Deep Learning Models

abstract: The Coronavirus Disease 2019 or the COVID-19 pandemic has swept almost all parts of the world since the first case was found in Wuhan, China, in December 2019. With the increasing number of COVID-19 cases in the world, SARS-CoV-2 has mutated into various variants. Given the increasingly dangerous conditions of the pandemic, it is crucial to know when the pandemic will stop by predicting confirmed cases of COVID-19. Therefore, many studies have raised COVID-19 as a case study to overcome the ongoing pandemic using the Deep Learning method, namely LSTM, with reasonably accurate results and small error values. LSTM training is used to predict confirmed cases of COVID-19 based on variants that have been identified using ECDC's COVID-19 dataset containing confirmed cases of COVID-19 that have been identified from 30 countries in Europe. Tests were conducted using the LSTM and BiLSTM models with the addition of RNN as comparisons on hidden size and layer size. The obtained result showed that in testing hidden sizes 25, 50, 75 to 100, the RNN model provided better results, with the minimum MSE value of 0.01 and the RMSE value of 0.012 for B.1.427/B.1.429 variant with hidden size 100. In further testing of layer sizes 2, 3, 4, and 5, the result shows that the BiLSTM model provided better results, with minimum MSE value of 0.01 and the RMSE of 0.01 for the B.1.427/B.1.429 variant with hidden size 100 and layer size 2.


title: Extended Nullstellensatz proof systems

abstract: For a finite set $\cal F$ of polynomials over fixed finite prime field of size $p$ containing all polynomials $x^2 - x$ a Nullstellensatz proof of the unsolvability of the system $$ f = 0\ ,\ \mbox{ all } f \in {\cal F} $$ in the field is a linear combination $\sum_{f \in {\cal F}} \ h_f \cdot f$ that equals to $1$ in the ring of polynomails. The measure of complexity of such a proof is its degree: $\max_f deg(h_f f)$. We study the problem to establish degree lower bounds for some {\em extended} NS proof systems: these systems prove the unsolvability of $\cal F$ by proving the unsolvability of a bigger set ${\cal F}\cup {\cal E}$, where set $\cal E$ may use new variables $r$ and contains all polynomials $r^p - r$, and satisfies the following soundness condition: -- - Any $0,1$-assignment $\overline a$ to variables $\overline x$ can be appended by an assignment $\overline b$ to variables $\overline r$ such that for all $g \in {\cal E}$ it holds that $g(\overline a, \overline b) = 0$. We define a notion of pseudo-solutions of $\cal F$ and prove that the existence of pseudo-solutions with suitable parameters implies lower bounds for two extended NS proof systems ENS and UENS defined in Buss et al. (1996/97). Further we give a combinatorial example of $\cal F$ and candidate pseudo-solutions based on the pigeonhole principle.


title: Clueless: A Tool Characterising Values Leaking as Addresses

abstract: Clueless is a binary instrumentation tool that characterises explicit cache side channel vulnerabilities of programs. It detects the transformation of data values into addresses by tracking dynamic instruction dependencies. Clueless tags data values in memory if it discovers that they are used in address calculations to further access other data. Clueless can report on the amount of data that are used as addresses at each point during execution. It can also be specifically instructed to track certain data in memory (e.g., a password) to see if they are turned into addresses at any point during execution. It returns a trace on how the tracked data are turned into addresses, if they do. We demonstrate Clueless on SPEC 2006 and characterise, for the first time, the amount of data values that are turned into addresses in these programs. We further demonstrate Clueless on a micro benchmark and on a case study. The case study is the different implementations of AES in OpenSSL: T-table, Vector Permutation AES (VPAES), and Intel Advanced Encryption Standard New Instructions (AES-NI). Clueless shows how the encryption key is transformed into addresses in the T-table implementation, while explicit cache side channel vulnerabilities are note detected in the other implementations.


title: Simultaneous Transmitting and Reflecting (STAR)-RIS for Harmonious Millimeter Wave Spectrum Sharing

abstract: The opening of the millimeter wave (mmWave) spectrum bands for 5G communications has motivated the need for novel spectrum sharing solutions at these high frequencies. In fact, reconfigurable intelligent surfaces (RISs) have recently emerged to enable spectrum sharing while enhancing the incumbents' quality-of-service (QoS). Nonetheless, co-existence over mmWave bands remains persistently challenging due to their unfavorable propagation characteristics. Hence, initiating mmWave spectrum sharing requires the RIS to further assist in improving the QoS over mmWave bands without jeopardizing spectrum sharing demands. In this paper, a novel simultaneous transmitting and reflecting RIS (STAR-RIS)-aided solution to enable mmWave spectrum sharing is proposed. In particular, the transmitting and reflecting abilities of the STAR-RIS are leveraged to tackle the mmWave spectrum sharing and QoS requirements separately. The STAR-RIS-enabled spectrum sharing problem between a primary network (e.g. a radar transmit-receive pair) and a secondary network is formulated as an optimization problem whose goal is to maximize the downlink sum-rate over a secondary multiple-input-single-output (MISO) network, while limiting interference over a primary network. Moreover, the STAR-RIS response coefficients and beamforming matrix in the secondary network are jointly optimized. To solve this non-convex problem, an alternating iterative algorithm is employed, where the STAR-RIS response coefficients and beamforming matrix are obtained using the successive convex approximation method. Simulation results show that the proposed solution outperforms conventional RIS schemes for mmWave spectrum sharing by achieving a 14.57% spectral efficiency gain.


title: Absolute continuity of self-similar measures on the plane

abstract: Consider an iterated function system consisting of similarities on the complex plane of the form $g_{i}(z) = \lambda_i z + t_i,\ \lambda_i, t_i \in \mathbb{C},\ |\lambda_i|<1, i=1,\ldots, k$. We prove that for almost every choice of $(\lambda_1, \ldots, \lambda_k)$ in the super-critical region (with fixed translations and probabilities), the corresponding self-similar measure is absolutely continuous. This extends results of Shmerkin-Solomyak (in the homogenous case) and Saglietti-Shmerkin-Solomyak (in the one-dimensional non-homogeneous case). As the main steps of the proof, we obtain results on the dimension and power Fourier decay of random self-similar measures on the plane, which may be of independent interest.


title: A signed count of 2-torsion points on real abelian varieties

abstract: We prove that a natural signed count of the $2$-torsion points on a real principally polarized abelian variety $A$ always equals to $2^{g}$ where $g$ is the dimension of $A$. When $A$ is the Jacobian of a real curve we derive signed counts of real odd theta characteristics. These can be interpreted in terms of the extrinsic geometry of contact hyperplanes to the canonical embedding of the curve. We also formulate a conjectural generalization to arbitrary fields in terms of $\mathbb{A}^1$-enumerative geometry.


title: An Approximate Algorithm for Maximum Inner Product Search over Streaming Sparse Vectors

abstract: Maximum Inner Product Search or top-k retrieval on sparse vectors is well-understood in information retrieval, with a number of mature algorithms that solve it exactly. However, all existing algorithms are tailored to text and frequency-based similarity measures. To achieve optimal memory footprint and query latency, they rely on the near stationarity of documents and on laws governing natural languages. We consider, instead, a setup in which collections are streaming -- necessitating dynamic indexing -- and where indexing and retrieval must work with arbitrarily distributed real-valued vectors. As we show, existing algorithms are no longer competitive in this setup, even against naive solutions. We investigate this gap and present a novel approximate solution, called Sinnamon, that can efficiently retrieve the top-k results for sparse real valued vectors drawn from arbitrary distributions. Notably, Sinnamon offers levers to trade-off memory consumption, latency, and accuracy, making the algorithm suitable for constrained applications and systems. We give theoretical results on the error introduced by the approximate nature of the algorithm, and present an empirical evaluation of its performance on two hardware platforms and synthetic and real-valued datasets. We conclude by laying out concrete directions for future research on this general top-k retrieval problem over sparse vectors.


title: Fourier decay of equilibrium states for bunched attractors

abstract: Let $M$ be a closed manifold, and let $f:M \rightarrow M$ be a $C^{2+\alpha}$ Axiom A diffeomorphism. Suppose that $f$ has an attractor $\Omega$ with codimension 1 stable lamination. Under a generic nonlinearity condition and a suitable bunching condition, we prove polynomial Fourier decay in the unstable direction for a large class of invariant measures on $\Omega$. Our result applies in particular for the measure of maximal entropy. We construct in the appendix an explicit solenoid that satisfies the nonlinearity and bunching assumption.


title: Energy-Delay Tradeoff in Helper-Assisted NOMA-MEC Systems: A Four-Sided Matching Algorithm

abstract: This paper designs a helper-assisted resource allocation strategy in non-orthogonal multiple access (NOMA)-enabled mobile edge computing (MEC) systems, in order to guarantee the quality of service (QoS) of the energy/delay-sensitive user equipments (UEs). To achieve a tradeoff between the energy consumption and the delay, we introduce a novel performance metric, called \emph{energy-delay tradeoff}, which is defined as the weighted sum of energy consumption and delay. The joint optimization of user association, resource block (RB) assignment, power allocation, task assignment, and computation resource allocation is formulated as a mixed-integer nonlinear programming problem with the aim of minimizing the maximal energy-delay tradeoff. Due to the non-convexity of the formulated problem with coupled and 0-1 variables, this problem cannot be directly solved with polynomial complexity. To tackle this challenge, we first decouple the formulated problem into a power allocation, task assignment and computation resource allocation (PATACRA) subproblem. Then, with the solution obtained from the PATACRA subproblem, we equivalently reformulate the original problem as a discrete user association and RB assignment (DUARA) problem. For the PATACRA subproblem, an iterative parametric convex approximation (IPCA) algorithm is proposed. Then, based on the solution obtained from the PATACRA subproblem, we first model the DUARA problem as a four-sided matching problem, and then propose a low-complexity four-sided UE-RB-helper-server matching (FS-URHSM) algorithm. Theoretical analysis demonstrates that the proposed algorithms are guaranteed to converge to stable solutions with polynomial complexity. Finally, simulation results are provided to show the superior performance of our proposed algorithm in terms of the energy consumption and the delay.


title: Navigating the Pitfalls of Active Learning Evaluation: A Systematic Framework for Meaningful Performance Assessment

abstract: Active Learning (AL) aims to reduce the labeling burden by interactively selecting the most informative samples from a pool of unlabeled data. While there has been extensive research on improving AL query methods in recent years, some studies have questioned the effectiveness of AL compared to emerging paradigms such as semi-supervised (Semi-SL) and self-supervised learning (Self-SL), or a simple optimization of classifier configurations. Thus, today's AL literature presents an inconsistent and contradictory landscape, leaving practitioners uncertain about whether and how to use AL in their tasks. In this work, we make the case that this inconsistency arises from a lack of systematic and realistic evaluation of AL methods. Specifically, we identify five key pitfalls in the current literature that reflect the delicate considerations required for AL evaluation. Further, we present an evaluation framework that overcomes these pitfalls and thus enables meaningful statements about the performance of AL methods. To demonstrate the relevance of our protocol, we present a large-scale empirical study and benchmark for image classification spanning various data sets, query methods, AL settings, and training paradigms. Our findings clarify the inconsistent picture in the literature and enable us to give hands-on recommendations for practitioners. The benchmark is hosted at https://github.com/IML-DKFZ/realistic-al .


title: Predicting chaotic statistics with unstable invariant tori

abstract: It has recently been speculated that statistical properties of chaos may be captured by weighted sums over unstable invariant tori embedded in the chaotic attractor of hyperchaotic dissipative systems; analogous to sums over periodic orbits formalized within periodic orbit theory. Using a novel numerical method for converging unstable invariant 2-tori in a chaotic PDE, we identify many quasiperiodic, unstable, invariant 2-torus solutions of a modified Kuramoto-Sivashinsky equation exhibiting hyperchaotic dynamics with two positive Lyapunov exponents. The set of tori covers significant parts of the chaotic attractor and weighted averages of the properties of the tori -- with weights computed based on their respective stability eigenvalues -- approximate statistics for the chaotic dynamics. These results are a step towards including higher-dimensional invariant sets in a generalized periodic orbit theory for hyperchaotic systems including spatio-temporally chaotic PDEs.


title: MV polytopes and reduced double Bruhat cells

abstract: When $G$ is a complex reductive algebraic group, MV polytopes are in bijection with the non-negative tropical points of the unipotent group of $G$. By fixing $w$ from the Weyl group, we can define MV polytopes whose highest vertex is labelled by $w$. We show that these polytopes are in bijection with the non-negative tropical points of the reduced double Bruhat cell labelled by $w^{-1}$. To do this, we define a collection of generalized minor functions $\Delta_\gamma^\text{new}$ which tropicalize on the reduced Bruhat cell to the BZ data of an MV polytope of highest vertex $w$. We also describe the combinatorial structure of MV polytopes of highest vertex $w$. We explicitly describe the map from the Weyl group to the subset of elements bounded by $w$ in the Bruhat order which sends $u \mapsto v$ if the vertex labelled by $u$ coincides with the vertex labelled by $v$ for every MV polytope of highest vertex $w$. As a consequence of this map, we prove that these polytopes have vertices labelled by Weyl group elements less than $w$ in the Bruhat order.


title: Incentive-weighted Anomaly Detection for False Data Injection Attacks Against Smart Meter Load Profiles

abstract: Spot pricing is often suggested as a method of increasing demand-side flexibility in electrical power load. However, few works have considered the vulnerability of spot pricing to financial fraud via false data injection (FDI) style attacks. In this paper, we consider attacks which aim to alter the consumer load profile to exploit intraday price dips. We examine an anomaly detection protocol for cyber-attacks that seek to leverage spot prices for financial gain. In this way we outline a methodology for detecting attacks on industrial load smart meters. We first create a feature clustering model of the underlying business, segregated by business type. We then use these clusters to create an incentive-weighted anomaly detection protocol for false data attacks against load profiles. This clustering-based methodology incorporates both the load profile and spot pricing considerations for the detection of injected load profiles. To reduce false positives, we model incentive-based detection, which includes knowledge of spot prices, into the anomaly tracking, enabling the methodology to account for changes in the load profile which are unlikely to be attacks.


title: Evaluating Deception and Moving Target Defense with Network Attack Simulation

abstract: In the field of network security, with the ongoing arms race between attackers, seeking new vulnerabilities to bypass defense mechanisms and defenders reinforcing their prevention, detection and response strategies, the novel concept of cyber deception has emerged. Starting from the well-known example of honeypots, many other deception strategies have been developed such as honeytokens and moving target defense, all sharing the objective of creating uncertainty for attackers and increasing the chance for the attacker of making mistakes. In this paper a methodology to evaluate the effectiveness of honeypots and moving target defense in a network is presented. This methodology allows to quantitatively measure the effectiveness in a simulation environment, allowing to make recommendations on how many honeypots to deploy and on how quickly network addresses have to be mutated to effectively disrupt an attack in multiple network and attacker configurations. With this optimum, attacks can be detected and slowed down with a minimal resource and configuration overhead. With the provided methodology, the optimal number of honeypots to be deployed and the optimal network address mutation interval can be determined. Furthermore, this work provides guidance on how to optimally deploy and configure them with respect to the attacker model and several network parameters.


title: Inference in Marginal Structural Models by Automatic Targeted Bayesian and Minimum Loss-Based Estimation

abstract: Two of the principle tasks of causal inference are to define and estimate the effect of a treatment on an outcome of interest. Formally, such treatment effects are defined as a possibly functional summary of the data generating distribution, and are referred to as target parameters. Estimation of the target parameter can be difficult, especially when it is high-dimensional. Marginal Structural Models (MSMs) provide a way to summarize such target parameters in terms of a lower dimensional working model. We introduce the semi-parametric efficiency bound for estimating MSM parameters in a general setting. We then present a frequentist estimator that achieves this bound based on Targeted Minimum Loss-Based Estimation. Our results are derived in a general context, and can be easily adapted to specific data structures and target parameters. We then describe a novel targeted Bayesian estimator and provide a Bernstein von-Mises type result analyzing its asymptotic behavior. We propose a universal algorithm that uses automatic differentiation to put the estimator into practice for arbitrary choice of working model. The frequentist and Bayesian estimators have been implemented in the Julia software package TargetedMSM.jl. Finally, we illustrate our proposed methods by investigating the effect of interventions on family planning behavior using data from a randomized field experiment conducted in Malawi.


title: Investigation of Planckian behavior in a high-conductivity oxide: PdCrO$_2$

abstract: The layered delafossite metal PdCrO$_2$ is a natural heterostructure of highly conductive Pd layers Kondo coupled to localized spins in the adjacent Mott insulating CrO$_2$ layers. At high temperatures $T$ it has a $T$-linear resistivity which is not seen in the isostructural but non-magnetic PdCoO$_2$. The strength of the Kondo coupling is known, as-grown crystals are extremely high purity and the Fermi surface is both very simple and experimentally known. It is therefore an ideal material platform in which to investigate 'Planckian metal' physics. We do this by means of controlled introduction of point disorder, measurement of the thermal conductivity and Lorenz ratio and studying the sources of its high temperature entropy. The $T$-linear resistivity is seen to be due mainly to elastic scattering and to arise from a sum of several scattering mechanisms. Remarkably, this sum leads to a scattering rate within 10$\%$ of the Planckian value of $k_BT/$$\hbar$.


title: EFX Exists for Four Agents with Three Types of Valuations

abstract: In this paper, we address the problem of determining an envy-free allocation of indivisible goods among multiple agents. EFX, which stands for envy-free up to any good, is a well-studied problem that has been shown to exist for specific scenarios, such as when there are only three agents with MMS valuations, as demonstrated by Chaudhury et al(2020), and for any number of agents when there are only two types of valuations as shown by Mahara(2020). Our contribution is to extend these results by showing that EFX exists for four agents with three distinct valuations. We further generalize this to show the existance of EFX allocations for n agents when n-2 of them have identical valuations.


title: PGD reduced-order modeling for structural dynamics applications

abstract: We propose in this paper a Proper Generalized Decomposition (PGD) approach for the solution of problems in linear elastodynamics. The novelty of the work lies in the development of weak formulations of the PGD problems based on the Lagrangian and Hamiltonian Mechanics, the main objective being to devise numerical methods that are numerically stable and energy conservative. We show that the methodology allows one to consider the Galerkin-based version of the PGD and numerically demonstrate that the PGD solver based on the Hamiltonian formulation offers better stability and energy conservation properties than the Lagrangian formulation. The performance of the two formulations is illustrated and compared on several numerical examples describing the dynamical behavior of a one-dimensional bar.


title: Freezing transition and moments of moments of the Riemann zeta function

abstract: Moments of moments of the Riemann zeta function, defined by \[ \text{MoM}_T (k,\beta) = \frac{1}{T} \int_T^{2T} \left( \int_{ |h|\leq (\log T)^\theta}|\zeta(\tfrac{1}{2} + i t + ih)|^{2\beta} dh \right)^k dt \] where $k,\beta \geq 0$ and $\theta > -1$, were introduced by Fyodorov and Keating when comparing extreme values of zeta in short intervals to those of characteristic polynomials of random unitary matrices. We study the $k = 2$ case as $T \rightarrow \infty$ and obtain sharp upper bounds for $\text{MoM}_T(2,\beta)$ for all real $0\leq \beta \leq 1$ as well as lower bounds of the conjectured order for all $\beta \geq 0$. In particular, we show that the second moment of moments undergoes a freezing phase transition with critical exponent $\beta = \tfrac{1}{\sqrt{2}}$.


title: Towards Mobility Management with Multi-Objective Bayesian Optimization

abstract: One of the consequences of network densification is more frequent handovers (HO). HO failures have a direct impact on the quality of service and are undesirable, especially in scenarios with strict latency, reliability, and robustness constraints. In traditional networks, HO-related parameters are usually tuned by the network operator, and automated techniques are still based on past experience. In this paper, we propose an approach for optimizing HO thresholds using Bayesian Optimization (BO). We formulate a multi-objective optimization problem for selecting the HO thresholds that minimize HOs too early and too late in indoor factory scenarios, and we use multi-objective BO (MOBO) for finding the optimal values. Our results show that MOBO reaches Pareto optimal solutions with few samples and ensures service continuation through safe exploration of new data points.


title: Equations for formal toric degenerations

abstract: Let $R$ be a complete equicharacteristic noetherian local domain and $\nu$ a valuation of its field of fractions whose valuation ring dominates $R$ with trivial residue field extension. The semigroup of values of $\nu$ on $R\setminus \{0\}$ is not finitely generated in general. We produce equations in an appropriate generalized power series ring for the algebra encoding the degeneration of $R$ to the toric graded algebra ${\rm gr}_\nu R$ associated to the filtration defined by $\nu$. We apply this to represent $\nu$ as the limit of a sequence of Abhyankar semivaluations (valuations on quotients) of $R$ with finitely generated semigroups.


title: Complexity of Geometric programming in the Turing model and application to nonnegative tensors

abstract: We consider a geometric programming problem consisting in minimizing a function given by the supremum of finitely many log-Laplace transforms of discrete nonnegative measures on a Euclidean space. Under a coerciveness assumption, we show that a $\varepsilon$-minimizer can be computed in a time that is polynomial in the input size and in $|\log\varepsilon|$. This is obtained by establishing bit-size estimates on approximate minimizers and by applying the ellipsoid method. We also derive polynomial iteration complexity bounds for the interior point method applied to the same class of problems. We deduce that the spectral radius of a partially symmetric, weakly irreducible nonnegative tensor can be approximated within $\varepsilon$ error in poly-time. For strongly irreducible tensors, we also show that the logarithm of the positive eigenvector is poly-time computable. Our results also yield that the the maximum of a nonnegative homogeneous $d$-form in the unit ball with respect to $d$-H\"older norm can be approximated in poly-time. In particular, the spectral radius of uniform weighted hypergraphs and some known upper bounds for the clique number of uniform hypergraphs are poly-time computable.


title: MLPGradientFlow: going with the flow of multilayer perceptrons (and finding minima fast and accurately)

abstract: MLPGradientFlow is a software package to solve numerically the gradient flow differential equation $\dot \theta = -\nabla \mathcal L(\theta; \mathcal D)$, where $\theta$ are the parameters of a multi-layer perceptron, $\mathcal D$ is some data set, and $\nabla \mathcal L$ is the gradient of a loss function. We show numerically that adaptive first- or higher-order integration methods based on Runge-Kutta schemes have better accuracy and convergence speed than gradient descent with the Adam optimizer. However, we find Newton's method and approximations like BFGS preferable to find fixed points (local and global minima of $\mathcal L$) efficiently and accurately. For small networks and data sets, gradients are usually computed faster than in pytorch and Hessian are computed at least $5\times$ faster. Additionally, the package features an integrator for a teacher-student setup with bias-free, two-layer networks trained with standard Gaussian input in the limit of infinite data. The code is accessible at https://github.com/jbrea/MLPGradientFlow.jl.


title: Low regularity error estimates for the time integration of 2D NLS

abstract: A filtered Lie splitting scheme is proposed for the time integration of the cubic nonlinear Schr\"odinger equation on the two-dimensional torus $\mathbb{T}^2$. The scheme is analyzed in a framework of discrete Bourgain spaces, which allows us to consider initial data with low regularity; more precisely initial data in $H^s(\mathbb{T}^2)$ with $s>0$. In this way, the usual stability restriction to smooth Sobolev spaces with index $s>1$ is overcome. Rates of convergence of order $\tau^{s/2}$ in $L^2(\mathbb{T}^2)$ at this regularity level are proved. Numerical examples illustrate that these convergence results are sharp.


title: Adaptive enrichment trial designs using joint modeling of longitudinal and time-to-event data

abstract: Adaptive enrichment allows for pre-defined patient subgroups of interest to be investigated throughout the course of a clinical trial. Many trials which measure a long-term time-to-event endpoint often also routinely collect repeated measures on biomarkers which may be predictive of the primary endpoint. Although these data may not be leveraged directly to support subgroup selection decisions and early stopping decisions, we aim to make greater use of these data to increase efficiency and improve interim decision making. In this work, we present a joint model for longitudinal and time-to-event data and two methods for creating standardised statistics based on this joint model. We can use the estimates to define enrichment rules and efficacy and futility early stopping rules for a flexible efficient clinical trial with possible enrichment. Under this framework, we show asymptotically that the familywise error rate is protected in the strong sense. To assess the results, we consider a trial for the treatment of metastatic breast cancer where repeated ctDNA measurements are available and the subgroup criteria is defined by patients' ER and HER2 status. Using simulation, we show that incorporating biomarker information leads to accurate subgroup identification and increases in power.


title: The upcoming spectroscopic powerhouses at the Isaac Newton Group of Telescopes

abstract: The Isaac Newton Group of Telescopes is completing a strategic change for the scientific use of its two telescopes, the 4.2-m William Herschel Telescope (WHT) and the 2.5-m Isaac Newton Telescope (INT). After more than 30 years operating as multi-purpose telescopes, the telescopes will soon complete their shift to nearly-single instrument operation dominated by large surveys. At the WHT, the WEAVE multi-fibre spectrograph is being commissioned in late 2022. Science surveys are expected to launch in 2023. 30% of the available time will be offered in open time. For the INT, construction of HARPS-3, a high-resolution ultra-stable spectrograph for extra-solar planet studies, is underway, with deployment planned for late 2024. The INT itself is being modernised and will operate as a robotic telescope. An average of 40% of the time will be offered as open time. The ING will maintain its student programme. Plans call for moving student work from the INT to the WHT once the INT starts operating robotically.


title: Group fairness in dynamic refugee assignment

abstract: Ensuring that refugees and asylum seekers thrive (e.g., find employment) in their host countries is a profound humanitarian goal, and a primary driver of employment is the geographic location within a host country to which the refugee or asylum seeker is assigned. Recent research has proposed and implemented algorithms that assign refugees and asylum seekers to geographic locations in a manner that maximizes the average employment across all arriving refugees. While these algorithms can have substantial overall positive impact, using data from two industry collaborators we show that the impact of these algorithms can vary widely across key subgroups based on country of origin, age, or educational background. Thus motivated, we develop a simple and interpretable framework for incorporating group fairness into the dynamic refugee assignment problem. In particular, the framework can flexibly incorporate many existing and future definitions of group fairness from the literature (e.g., maxmin, randomized, and proportionally-optimized within-group). Equipped with our framework, we propose two bid-price algorithms that maximize overall employment while simultaneously yielding provable group fairness guarantees. Through extensive numerical experiments using various definitions of group fairness and real-world data from the U.S. and the Netherlands, we show that our algorithms can yield substantial improvements in group fairness compared to an offline benchmark fairness constraints, with only small relative decreases ($\approx$ 1%-5%) in global performance.


title: Automatic Locally Robust Estimation with Generated Regressors

abstract: Many economic and causal parameters of interest depend on generated regressors. Examples include structural parameters in models with endogenous variables estimated by control functions and in models with sample selection, treatment effect estimation with propensity score matching, and marginal treatment effects. Inference with generated regressors is complicated by the very complex expression for influence functions and asymptotic variances. To address this problem, we propose Automatic Locally Robust/debiased GMM estimators in a general setting with generated regressors. Importantly, we allow for the generated regressors to be generated from machine learners, such as Random Forest, Neural Nets, Boosting, and many others. We use our results to construct novel Doubly Robust and Locally Robust estimators for the Counterfactual Average Structural Function and Average Partial Effects in models with endogeneity and sample selection, respectively. We provide sufficient conditions for the asymptotic normality of our debiased GMM estimators and investigate their finite sample performance through Monte Carlo simulations.


title: SALT2 Versus SALT3: Updated Model Surfaces and Their Impacts on Type Ia Supernova Cosmology

abstract: For the past decade, SALT2 has been the most common model used to fit Type Ia supernova (SN Ia) light curves for dark energy analyses. Recently, the SALT3 model was released, which upgraded a number of model features but has not yet been used for measurements of dark energy. Here, we evaluate the impact of switching from SALT2 to SALT3 for a SN cosmology analysis. We train SALT2 and SALT3 on an identical training sample of 1083 well-calibrated Type Ia supernovae, ensuring that any differences found come from the underlying model framework. We publicly release the results of this training (the SALT "surfaces"). We then run a cosmology analysis on the public Dark Energy Survey 3-Year Supernova data sample (DES-SN3YR), and on realistic simulations of those data. We provide the first estimate of the SN+CMB systematic uncertainty arising from the choice of SALT model framework (i.e. SALT2 versus SALT3), $\Delta w = +0.001 \pm 0.005$ -- a negligible effect at the current level of dark energy analyses. We also find that the updated surfaces are less sensitive to photometric calibration uncertainties than previous SALT2 surfaces, with the average spectral energy density dispersion reduced by a factor of two over optical wavelengths. This offers an opportunity to reduce the contribution of calibration errors to SN cosmology uncertainty budgets.


title: Accurate evaluation of self-heterodyne laser linewidth measurements using Wiener filters

abstract: Self-heterodyne beat note measurements are widely used for the experimental characterization of the frequency noise power spectral density (FN-PSD) and the spectral linewidth of lasers. The measured data, however, must be corrected for the transfer function of the experimental setup in a post-processing routine. The standard approach disregards the detector noise and thereby induces reconstruction artifacts, i.e., spurious spikes, in the reconstructed FN-PSD. We introduce an improved post-processing routine based on a parametric Wiener filter that is free from reconstruction artifacts, provided a good estimate of the signal-to-noise ratio is supplied. Building on this potentially exact reconstruction, we develop a new method for intrinsic laser linewidth estimation that is aimed at deliberate suppression of unphysical reconstruction artifacts. Our method yields excellent results even in the presence of strong detector noise, where the intrinsic linewidth plateau is not even visible using the standard method. The approach is demonstrated for simulated time series from a stochastic laser model including 1/f-type noise.


title: Constrained Expectation-Maximisation for inference of social graphs explaining online user-user interactions

abstract: Current network inference algorithms fail to generate graphs with edges that can explain whole sequences of node interactions in a given dataset or trace. To quantify how well an inferred graph can explain a trace, we introduce feasibility, a novel quality criterion, and suggest that it is linked to the result's accuracy. In addition, we propose CEM-*, a network inference method that guarantees 100% feasibility given online social media traces, which is a non-trivial extension of the Expectation-Maximization algorithm developed by Newman (2018). We propose a set of linear optimization updates that incorporate a set of auxiliary variables and a set of feasibility constraints; the latter takes into consideration all the hidden paths that are possible between users based on their timestamps of interaction and guide the inference toward feasibility. We provide two CEM-* variations, that assume either an Erdos Renyi (ER) or a Stochastic Block Model (SBM) prior for the underlying graph's unknown distribution. Extensive experiments on one synthetic and one real-world Twitter dataset show that for both priors CEM-* can generate a posterior distribution of graphs that explains the whole trace while being closer to the ground truth. As an additional benefit, the use of the SBM prior infers and clusters users simultaneously during optimization. CEM-* outperforms baseline and state-of-the-art methods in terms of feasibility, run-time, and precision of the inferred graph and communities. Finally, we propose a heuristic to adapt the inference to lower feasibility requirements and show how it can affect the precision of the result.


title: Finite alphabet phase retrieval

abstract: We consider the finite alphabet phase retrieval problem: recovering a signal whose entries lie in a small alphabet of possible values from its Fourier magnitudes. This problem arises in the celebrated technology of X-ray crystallography to determine the atomic structure of biological molecules. Our main result states that for generic values of the alphabet, two signals have the same Fourier magnitudes if and only if several partitions have the same difference sets. Thus, the finite alphabet phase retrieval problem reduces to the combinatorial problem of determining a signal from those difference sets. Notably, this result holds true when one of the letters of the alphabet is zero, namely, for sparse signals with finite alphabet, which is the situation in X-ray crystallography.


title: Endpoint mixed weak type extrapolation

abstract: The purpose of this note is to extend the extrapolation result by by Cruz-Uribe Martell and P\'erez as follows. Given a family $\mathcal{F}$ of pairs of functions suppose that for some $0<p<\infty$ and for every $w\in A_{\infty}$ \begin{equation} \int f^{p}w\leq c_{w}\int g^{p}w\qquad(f,g)\in\mathcal{F}\label{eq:Hip-1} \end{equation} provided the left-hand side of the estimate is finite. If we have that $A(t)=t\log(e+t)^{\rho}$ for some $\rho>0$, then, for every $u\in A_{1}$ and every $v\in A_{\infty}$ we have that \[ \left\Vert \frac{f}{v}\right\Vert_{L^{A,\infty}(uv)}\lesssim\left\Vert \frac{g}{v}\right\Vert_{L^{A,\infty}(uv)}, \] where \[ L^{A,\infty}(uv)=\inf\left\{ \lambda>0:\sup_{t>0}A(t)w\left(\left\{ x\in\mathbb{R}:|f(x)|>\lambda t\right\} \right)\leq1\right\} \] is the weak Orlicz type introduced by Iaffei. As a corollary of this extrapolation result we derive a mixed weak type inequality for Coifman-Rochberg-Weiss commutators.


title: On Creating a Comprehensive Food Database

abstract: Studies with the primary aim of addressing eating disorders focus on assessing the nutrient content of food items with an exclusive focus on caloric intake. There are two primary impediments that can be noted in these studies. The first of these relates to the fact that caloric intake of each food item is calculated from an existing database. The second concerns the scientific significance of caloric intake used as the single measure of nutrient content. By requiring an existing database, researchers are forced to find some source of a comprehensive set of food items as well as their respective nutrients. This search alone is a difficult task, and if completed often leads to the requirement of a paid API service. These services are expensive and non-customizable, taking away funding that could be aimed at other parts of the study only to give an unwieldy database that can not be modified or contributed to. In this work, we introduce a new rendition of the USDA's food database that includes both foods found in grocery stores and those found in restaurants or fast food places. At the moment, we have accumulated roughly 1.5 million food entries consisting of approximately 18,000 brands and 100 restaurants in the United States. These foods also have an abundance of nutrient data associated with them, from the caloric amount to saturated fat levels. The data is stored in MySQL format and is spread among five major tables. We have also procured images for theses foods entries when available, and have included all of our data and program scripts in an open source repository.


title: Hydrodynamic Collapse of the Leidenfrost Vapor Layer

abstract: During the Leidenfrost effect, a stable vapor film can separate a hot solid from an evaporating liquid. Eventually, after formation and upon cooling, the vapor layer cannot be sustained and undergoes a violent collapse evidenced by explosive boiling. Computationally, modeling this instability involves an interplay between hydrodynamics, thermodynamics, rapid evaporation, and length-scales from $\mu$m to cm. Selective assumptions, made to reduce computational costs, have limited most previous studies to steady-state investigations. Here, we combine two-phase laminar flow, heat transfer, and evaporation in a finite-element simulation to examine the failure of Leidenfrost vapor layers during cooling. During periods of quiescence, the geometry of the vapor layer agrees well with steady-state lubrication theory. In the simulations, we report the local temperature of the solid at failure, $T_-$, which provides a lower bound for recent experimental work using the same geometric and material conditions. Surprisingly, we find that inertial forces, which are typically ignored in theoretical treatments of the vapor layer, are responsible for initiating the instability leading to failure.


title: Overcoming Prior Misspecification in Online Learning to Rank

abstract: The recent literature on online learning to rank (LTR) has established the utility of prior knowledge to Bayesian ranking bandit algorithms. However, a major limitation of existing work is the requirement for the prior used by the algorithm to match the true prior. In this paper, we propose and analyze adaptive algorithms that address this issue and additionally extend these results to the linear and generalized linear models. We also consider scalar relevance feedback on top of click feedback. Moreover, we demonstrate the efficacy of our algorithms using both synthetic and real-world experiments.


title: Non-contact acoustic micro-tapping optical coherence elastography for quantification of corneal anisotropic elasticity: in vivo rabbit study

abstract: Purpose. To demonstrate accurate measurement of corneal elastic moduli in vivo with non-contact and non-invasive optical coherence elastography. Methods. Elastic properties (in-plane Young's modulus E and both in-plane, u, and out-of-plane, G, shear moduli) of rabbit cornea were quantified in vivo using non-contact dynamic Acoustic micro-Tapping Optical Coherence Elastography (AuT-OCE). The IOP-dependence of measured mechanical properties was explored in extracted whole globes following in vivo measurement. A nearly-incompressible transverse isotropic (NITI) model was used to reconstruct moduli from AuT-OCE data. Independently, cornea elastic moduli were also measured ex vivo with traditional, destructive mechanical tests (tensile extensometry and shear rheometry). Results. Our study demonstrates strong anisotropy of corneal elasticity in rabbits. The in-plane Young's modulus, computer as E=3u, was in the range of 20-44 MPa, whereas the out-of-plane shear modulus was in the range of 34-261 kPa. Both pressure-dependent ex vivo OCE and destructive mechanical tests performed on the same samples within an hour of euthanasia strongly support the results of AuT-OCE measurements. Conclusions. Non-contact AuT-OCE can non-invasively quantify cornea anisotropic elastic properties in vivo. Translational Relevance. As OCT is broadly accepted in Ophthalmology, these results suggest the potential for rapid translation of AuT-OCE into clinical practice. In addition, AuT-OCE can likely improve diagnostic criteria of ectatic corneal diseases, leading to early diagnosis, reduced complications, customized surgical treatment, and personalized biomechanical models of the eye.


title: Polar magneto-optic Kerr and Faraday effects in finite periodic \texorpdfstring{$\mathcal{P}\mathcal{T}$}{PT}-symmetric systems

abstract: We discuss the anomalous behavior of the Faraday (transmission) and polar Kerr (reflection) rotation angles of the propagating light, in finite periodic parity-time ($\mathcal{P}\mathcal{T}$) symmetric structures, consisting of $N$ cells. The unit cell potential is two complex $\delta$-potentials placed on both boundaries of the ordinary dielectric slab. It is shown that, for a given set of parameters describing the system, a phase transition-like anomalous behavior of Faraday and Kerr rotation angles in a parity-time symmetric systems can take place. In the anomalous phase the value of one of the Faraday and Kerr rotation angles can become negative, and both angles suffer from spectral singularities and give a strong enhancement near the singularities. We also shown that the real part of the complex angle of KR, $\theta^{R}_1$, is always equal to the $\theta^{T}_1$ of FR, no matter what phase the system is in due to the symmetry constraints. The imaginary part of KR angles $\theta^{R^{r/l}}_2$ are related to the $\theta^{T}_2$ of FR by parity-time symmetry. Calculations based on the approach of the generalized nonperturbative characteristic determinant, which is valid for a layered system with randomly distributed delta potentials, show that the Faraday and Kerr rotation spectrum in such structures has several resonant peaks. Some of them coincide with transmission peaks, providing simultaneous large Faraday and Kerr rotations enhanced by an order one or two of magnitude. We provide a recipe for funding a one-to-one relation in between KR and FR.