title: The effect of biomimetic mineralization of 3D-printed mesoporous bioglass scaffolds on physical properties and in vitro osteogenicity

abstract: Three-dimensional Mesoporous bioactive glasses (MBGs) scaffolds has been widely considered for bone regeneration purposes and additive manufacturing enables the fabrication of highly bioactive patient-specific constructs for bone defects. Commonly, this process is performed with the addition of polymeric binders that facilitate the printability of scaffolds. However, these additives cover the MBG particles resulting in the reduction of their osteogenic potential. The present work investigates a simple yet effective phosphate-buffered saline immersion method for achieving polyvinyl alcohol binder removal while enables the maintenance of the mesoporous structure of MBG 3D-printed scaffolds. This resulted in significantly modifying the surface of the scaffold via the spontaneous formation of a biomimetic mineralized layer which positively affected the physical and biological properties of the scaffold. The extensive surface remodeling induced by the deposition of the apatite-like layer lead to a 3-fold increase in surface area, a 5-fold increase in the roughness, and 4-fold increase in the hardness of the PBS-immersed scaffolds when compared to the as-printed counterpart. The biomimetic mineralization also occurred throughout the bulk of the scaffold connecting the MBGs particles and was responsible for the maintenance of structural integrity. In vitro assays using MC3T3-E1 pre-osteoblast like cells demonstrated a significant upregulation of osteogenic-related genes for the scaffolds previously immersed in PBS when compared to the as-printed PVA-containing scaffolds. Although the pre-immersion scaffolds performed equally towards osteogenic cell differentiation, our data suggest that a short immersion in PBS of MBG scaffolds is beneficial for the osteogenic properties and might accelerate bone formation after implantation.


title: Josephson dynamical simulation using the electronic circuit simulator APLAC: a tutorial

abstract: Analysis Program for Linear Active Circuits (APLAC) is a general-purpose electronics circuit simulator, which has included a built-in model of the Josephson junction (JJ) since late 80's. It capabilities in simulating eg. noisy Superconducting Quantum Interference Devices (SQUIDs), Rapid Single Flux Quantum (RSFQ) logic circuits, or Superconducting Transition Edge Sensors (TESes) are relatively unknown within the superconducting electronics community. Here we give a brief step-to-step tutorial for APLAC users to unleash those capabilities.


title: Evidence for local spots of viscous electron flow in graphene at moderate mobility

abstract: Dominating electron-electron scattering enables viscous electron flow exhibiting hydrodynamic current density patterns such as Poiseuille profiles or vortices. The viscous regime has recently been observed in graphene by non-local transport experiments and mapping of the Poiseuille profile. Here, we probe the current-induced surface potential maps of graphene field effect transistors with moderate mobility using scanning probe microscopy at room temperature. We discover micron-sized large areas appearing close to charge neutrality that show current induced electric fields opposing the externally applied field. By estimating the local scattering lengths from the gate dependence of local in-plane electric fields, we find that electron-electron scattering dominates in these areas as expected for viscous flow. Moreover, we suppress the inverted fields by artificially decreasing the electron-disorder scattering length via mild ion bombardment. These results imply that viscous electron flow is omnipresent in graphene devices, even at moderate mobility.

categories: cond-mat.mes-hall cond-mat.dis-nn cond-mat.mtrl-sci cond-mat.other


title: Robust Cell-Load Learning with a Small Sample Set

abstract: Learning of the cell-load in radio access networks (RANs) has to be performed within a short time period. Therefore, we propose a learning framework that is robust against uncertainties resulting from the need for learning based on a relatively small training sample set. To this end, we incorporate prior knowledge about the cell-load in the learning framework. For example, an inherent property of the cell-load is that it is monotonic in downlink (data) rates. To obtain additional prior knowledge we first study the feasible rate region, i.e., the set of all vectors of user rates that can be supported by the network. We prove that the feasible rate region is compact. Moreover, we show the existence of a Lipschitz function that maps feasible rate vectors to cell-load vectors. With these results in hand, we present a learning technique that guarantees a minimum approximation error in the worst-case scenario by using prior knowledge and a small training sample set. Simulations in the network simulator NS3 demonstrate that the proposed method exhibits better robustness and accuracy than standard multivariate learning techniques, especially for small training sample sets.


title: Learning Multi-Scene Absolute Pose Regression with Transformers

abstract: Absolute camera pose regressors estimate the position and orientation of a camera from the captured image alone. Typically, a convolutional backbone with a multi-layer perceptron head is trained with images and pose labels to embed a single reference scene at a time. Recently, this scheme was extended for learning multiple scenes by replacing the MLP head with a set of fully connected layers. In this work, we propose to learn multi-scene absolute camera pose regression with Transformers, where encoders are used to aggregate activation maps with self-attention and decoders transform latent features and scenes encoding into candidate pose predictions. This mechanism allows our model to focus on general features that are informative for localization while embedding multiple scenes in parallel. We evaluate our method on commonly benchmarked indoor and outdoor datasets and show that it surpasses both multi-scene and state-of-the-art single-scene absolute pose regressors. We make our code publicly available from https://github.com/yolish/multi-scene-pose-transformer.


title: Fairmandering: A column generation heuristic for fairness-optimized political districting

abstract: The American winner-take-all congressional district system empowers politicians to engineer electoral outcomes by manipulating district boundaries. Existing computational solutions mostly focus on drawing unbiased maps by ignoring political and demographic input, and instead simply optimize for compactness. We claim that this is a flawed approach because compactness and fairness are orthogonal qualities, and introduce a scalable two-stage method to explicitly optimize for arbitrary piecewise-linear definitions of fairness. The first stage is a randomized divide-and-conquer column generation heuristic which produces an exponential number of distinct district plans by exploiting the compositional structure of graph partitioning problems. This district ensemble forms the input to a master selection problem to choose the districts to include in the final plan. Our decoupled design allows for unprecedented flexibility in defining fairness-aligned objective functions. The pipeline is arbitrarily parallelizable, is flexible to support additional redistricting constraints, and can be applied to a wide array of other regionalization problems. In the largest ever ensemble study of congressional districts, we use our method to understand the range of possible expected outcomes and the implications of this range on potential definitions of fairness.


title: NeBula: Quest for Robotic Autonomy in Challenging Environments; TEAM CoSTAR at the DARPA Subterranean Challenge

abstract: This paper presents and discusses algorithms, hardware, and software architecture developed by the TEAM CoSTAR (Collaborative SubTerranean Autonomous Robots), competing in the DARPA Subterranean Challenge. Specifically, it presents the techniques utilized within the Tunnel (2019) and Urban (2020) competitions, where CoSTAR achieved 2nd and 1st place, respectively. We also discuss CoSTAR's demonstrations in Martian-analog surface and subsurface (lava tubes) exploration. The paper introduces our autonomy solution, referred to as NeBula (Networked Belief-aware Perceptual Autonomy). NeBula is an uncertainty-aware framework that aims at enabling resilient and modular autonomy solutions by performing reasoning and decision making in the belief space (space of probability distributions over the robot and world states). We discuss various components of the NeBula framework, including: (i) geometric and semantic environment mapping; (ii) a multi-modal positioning system; (iii) traversability analysis and local planning; (iv) global motion planning and exploration behavior; (i) risk-aware mission planning; (vi) networking and decentralized reasoning; and (vii) learning-enabled adaptation. We discuss the performance of NeBula on several robot types (e.g. wheeled, legged, flying), in various environments. We discuss the specific results and lessons learned from fielding this solution in the challenging courses of the DARPA Subterranean Challenge competition.


title: Conditional Generative Adversarial Networks for Speed Control in Trajectory Simulation

abstract: Motion behaviour is driven by several factors -- goals, presence and actions of neighbouring agents, social relations, physical and social norms, the environment with its variable characteristics, and further. Most factors are not directly observable and must be modelled from context. Trajectory prediction, is thus a hard problem, and has seen increasing attention from researchers in the recent years. Prediction of motion, in application, must be realistic, diverse and controllable. In spite of increasing focus on multimodal trajectory generation, most methods still lack means for explicitly controlling different modes of the data generation. Further, most endeavours invest heavily in designing special mechanisms to learn the interactions in latent space. We present Conditional Speed GAN (CSG), that allows controlled generation of diverse and socially acceptable trajectories, based on user controlled speed. During prediction, CSG forecasts future speed from latent space and conditions its generation based on it. CSG is comparable to state-of-the-art GAN methods in terms of the benchmark distance metrics, while being simple and useful for simulation and data augmentation for different contexts such as fast or slow paced environments. Additionally, we compare the effect of different aggregation mechanisms and show that a naive approach of concatenation works comparable to its attention and pooling alternatives.


title: 3-Lie Algebras, Ternary Nambu-Lie algebras and the Yang-Baxter equation

abstract: We construct ternary self-distributive (TSD) objects from compositions of binary Lie algebras, $3$-Lie algebras and, in particular, ternary Nambu-Lie algebras. We show that the structures obtained satisfy an invertibility property resembling that of racks. We prove that these structures give rise to Yang-Baxter operators in the tensor product of the base vector space and, upon defining suitable twisting isomorphisms, we obtain representations of the infinite (framed) braid group. We consider examples for low-dimensional Lie algebras, where the ternary bracket is defined by composition of the binary ones, along with simple $3$-Lie algebras. We show that the Yang-Baxter operators obtained are not gauge equivalent to the transposition operator, and we consider the problem of deforming the operators to obtain new solutions to the Yang-Baxter equation. We discuss the applications of this deformation procedure to the construction of (framed) link invariants.


title: Construction of non-trivial relativistic quantum fields in arbitrary space-time dimension via superposition of free fields

abstract: We construct Schwinger functions as the superposition of Schwinger functions which correspond to those of free fields with sharp masses $ m $. We prove that all axioms of Osterwalder and Schrader are satisfied. This construction works independently of the space-time dimension $ d $. The Schwinger functions under consideration are the moments of a non-Gaussian measure on the space of tempered distributions $ S'(\mathbb R^d) $.


title: #PraCegoVer: A Large Dataset for Image Captioning in Portuguese

abstract: Automatically describing images using natural sentences is an important task to support visually impaired people's inclusion onto the Internet. It is still a big challenge that requires understanding the relation of the objects present in the image and their attributes and actions they are involved in. Then, visual interpretation methods are needed, but linguistic models are also necessary to verbally describe the semantic relations. This problem is known as Image Captioning. Although many datasets were proposed in the literature, the majority contains only English captions, whereas datasets with captions described in other languages are scarce. Recently, a movement called PraCegoVer arose on the Internet, stimulating users from social media to publish images, tag #PraCegoVer and add a short description of their content. Thus, inspired by this movement, we have proposed the #PraCegoVer, a multi-modal dataset with Portuguese captions based on posts from Instagram. It is the first large dataset for image captioning in Portuguese with freely annotated images. Further, the captions in our dataset bring additional challenges to the problem: first, in contrast to popular datasets such as MS COCO Captions, #PraCegoVer has only one reference to each image; also, both mean and variance of our reference sentence length are significantly greater than those in the MS COCO Captions. These two characteristics contribute to making our dataset interesting due to the linguistic aspect and the challenges that it introduces to the image captioning problem. We publicly-share the dataset at https://github.com/gabrielsantosrv/PraCegoVer.


title: Implementable coupling of L\'evy process and Brownian motion

abstract: We provide a simple algorithm for construction of Brownian paths approximating those of a L\'evy process on a finite time interval. It requires knowledge of the L\'evy process trajectory on a chosen regular grid and the law of its endpoint, or the ability to simulate from that. This algorithm is based on reordering of Brownian increments, and it can be applied in a recursive manner. We establish an upper bound on the mean squared maximal distance between the paths and determine a suitable mesh size in various asymptotic regimes. The analysis proceeds by reduction to the comonotonic coupling of increments. Applications to model risk and multilevel Monte Carlo are discussed in detail, and numerical examples are provided.


title: Under pressure: Hydrogel swelling in a granular medium

abstract: Hydrogels hold promise in agriculture as reservoirs of water in dry soil, potentially alleviating the burden of irrigation. However, confinement in soil can markedly reduce the ability of hydrogels to absorb water and swell, limiting their widespread adoption. Unfortunately, the underlying reason remains unknown. By directly visualizing the swelling of hydrogels confined in three-dimensional granular media, we demonstrate that the extent of hydrogel swelling is determined by the competition between the force exerted by the hydrogel due to osmotic swelling and the confining force transmitted by the surrounding grains. Furthermore, the medium can itself be restructured by hydrogel swelling, as set by the balance between the osmotic swelling force, the confining force, and intergrain friction. Together, our results provide quantitative principles to predict how hydrogels behave in confinement, potentially improving their use in agriculture as well as informing other applications such as oil recovery, construction, mechanobiology, and filtration.

categories: cond-mat.soft cond-mat.dis-nn cond-mat.mtrl-sci physics.geo-ph


title: Paying Attention to Activation Maps in Camera Pose Regression

abstract: Camera pose regression methods apply a single forward pass to the query image to estimate the camera pose. As such, they offer a fast and light-weight alternative to traditional localization schemes based on image retrieval. Pose regression approaches simultaneously learn two regression tasks, aiming to jointly estimate the camera position and orientation using a single embedding vector computed by a convolutional backbone. We propose an attention-based approach for pose regression, where the convolutional activation maps are used as sequential inputs. Transformers are applied to encode the sequential activation maps as latent vectors, used for camera pose regression. This allows us to pay attention to spatially-varying deep features. Using two Transformer heads, we separately focus on the features for camera position and orientation, based on how informative they are per task. Our proposed approach is shown to compare favorably to contemporary pose regressors schemes and achieves state-of-the-art accuracy across multiple outdoor and indoor benchmarks. In particular, to the best of our knowledge, our approach is the only method to attain sub-meter average accuracy across outdoor scenes. We make our code publicly available from here.


title: Multiscale dynamics of colloidal deposition and erosion in porous media

abstract: Diverse processes -- e.g., environmental pollution, groundwater remediation, oil recovery, filtration, and drug delivery -- involve the transport of colloidal particles in porous media. Using confocal microscopy, we directly visualize this process in situ and thereby identify the fundamental mechanisms by which particles are distributed throughout a medium. At high injection pressures, hydrodynamic stresses cause particles to be continually deposited on and eroded from the solid matrix -- notably, forcing them to be distributed throughout the entire medium. By contrast, at low injection pressures, the relative influence of erosion is suppressed, causing particles to localize near the inlet of the medium. Unexpectedly, these macroscopic distribution behaviors depend on imposed pressure in similar ways for particles of different charges, although the pore-scale distribution of deposition is sensitive to particle charge. These results reveal how the multiscale interactions between fluid, particles, and the solid matrix control how colloids are distributed in a porous medium.

categories: cond-mat.soft cond-mat.dis-nn cond-mat.mes-hall cond-mat.mtrl-sci physics.flu-dyn


title: Spectral control of high order harmonics through non-linear propagation effects

abstract: High harmonic generation (HHG) in crystals has revealed a wealth of perspectives such as all-optical mapping of the electronic band structure, ultrafast quantum information and the creation of novel all-solid-state attosecond sources. Significant efforts have been made to understand the microscopic aspects of HHG in crystals, whereas the macroscopic effects, such as non-linear propagation effects of the driving pulse inside the dense solid media and its impact on the HHG process is often overlooked. In this work, we study macroscopic effects by comparing two materials with distinct optical properties, silicon (Si) and zinc oxide (ZnO). By scanning the focal position of 85 fs, 2.123 $\mu$m wavelength pulses inside the crystals (Z-scan) we reveal spectral shifts in the generated harmonics. We interpret the overall blueshift of the emitted harmonic spectrum as an imprint of the driving field spectral modulation occurring during the propagation inside the crystal. This is supported with numerical simulations. This study demonstrates that through manipulation of the fundamental driving field through non-linear propagation effects, precise control of the emitted HHG spectrum in solids can be realised. This method could offer a robust way to tailor HHG spectra for a range of spectroscopic applications.


title: Monadic Intuitionistic and Modal Logics Admitting Provability Interpretations

abstract: The G\"odel translation provides an embedding of the intuitionistic logic $\mathsf{IPC}$ into the modal logic $\mathsf{Grz}$, which then embeds into the modal logic $\mathsf{GL}$ via the splitting translation. Combined with Solovay's theorem that $\mathsf{GL}$ is the modal logic of the provability predicate of Peano Arithmetic $\mathsf{PA}$, both $\mathsf{IPC}$ and $\mathsf{Grz}$ admit arithmetical interpretations. When attempting to 'lift' these results to the monadic extensions $\mathsf{MIPC}$, $\mathsf{MGrz}$, and $\mathsf{MGL}$ of these logics, the same techniques no longer work. Following a conjecture made by Esakia, we add an appropriate version of Casari's formula to these monadic extensions (denoted by a '+'), obtaining that the G\"odel translation embeds $\mathsf{M^{+}IPC}$ into $\mathsf{M^{+}Grz}$ and the splitting translation embeds $\mathsf{M^{+}Grz}$ into $\mathsf{MGL}$. As proven by Japaridze, Solovay's result extends to the monadic system $\mathsf{MGL}$, which leads us to an arithmetical interpretation of both $\mathsf{M^{+}IPC}$ and $\mathsf{M^{+}Grz}$.


title: How do we Evaluate Self-adaptive Software Systems?

abstract: With the increase of research in self-adaptive systems, there is a need to better understand the way research contributions are evaluated. Such insights will support researchers to better compare new findings when developing new knowledge for the community. However, so far there is no clear overview of how evaluations are performed in self-adaptive systems. To address this gap, we conduct a mapping study. The study focuses on experimental evaluations published in the last decade at the prime venue of research in software engineering for self-adaptive systems -- the International Symposium on Software Engineering for Adaptive and Self-Managing Systems (SEAMS). Results point out that specifics of self-adaptive systems require special attention in the experimental process, including the distinction of the managing system (i.e., the target of evaluation) and the managed system, the presence of uncertainties that affect the system behavior and hence need to be taken into account in data analysis, and the potential of managed systems to be reused across experiments, beyond replications. To conclude, we offer a set of suggestions derived from our study that can be used as input to enhance future experiments in self-adaptive systems.


title: Heat kernel bounds for parabolic equations with singular (form-bounded) vector fields

abstract: We consider Kolmogorov operator $-\nabla \cdot a \cdot \nabla + b \cdot \nabla$ with measurable uniformly elliptic matrix $a$ and prove Gaussian lower and upper bounds on its heat kernel under minimal assumptions on the vector field $b$ and its divergence ${\rm div\,}b$. More precisely, we prove: (1) Gaussian lower bound, provided that ${\rm div\,}b \geq 0$, and $b$ is in the class of form-bounded vector fields (containing e.g.\,the class $L^d$, the weak $L^d$ class, as well as some vector fields that are not even in $L_{\rm loc}^{2+\varepsilon}$, $\varepsilon>0$); in these assumptions, the Gaussian upper bound is in general invalid; (2) Gaussian upper bound, provided that $b$ is form-bounded, and the positive part of ${\rm div\,}b$ is in the Kato class; in these assumptions, the Gaussian lower bound is in general invalid; (3) Gaussian upper and lower bounds, provided that $b$ is form-bounded, ${\rm div\,}b$ is in the Kato class; (4) A priori Gaussian upper and lower bounds, provided that $b$ is in a large class containing the class of form-bounded vector fields, ${\rm div\,}b$ is in the Kato class.


title: Galactic internet made possible by star gravitational lensing

abstract: Galactic Internet may already exist, if all stars are exploited as gravitational lenses. In fact, the gravitational lens of the Sun is a well-known astrophysical phenomenon predicted by Einstein's general theory of relativity. It implies that, if we can send a probe along any radial direction away from the Sun up to the minimal distance of 550 AU and beyond, the Sun's mass will act as a huge magnifying lens, letting us "see" detailed radio maps of whatever may lie on the other side of the Sun even at very large distances. The 2009 book by this author, ref. [1], studies such future FOCAL space missions to 550 AU and beyond. In this paper, however, we want to study another possibility yet: how to create the future interstellar radio links between the solar system and any future interstellar probe by utilizing the gravitational lens of the Sun as a huge antenna. In particular, we study the Bit Error Rate (BER) across interstellar distances with and without using the gravitational lens effect of the Sun (ref. [2]). The conclusion is that only when we will exploit the Sun as a gravitational lens we will be able to communicate with our own probes (or with nearby Aliens) across the distances of even the nearest stars to us in the Galaxy, and that at a reasonable Bit Error Rate. We also study the radio bridge between the Sun and any other Star that is made up by the two gravitational lenses of both the Sun and that Star. The alignment for this radio bridge to work is very strict, but the power-saving is enormous, due to the huge contributions of the two stars' lenses to the overall antenna gain of the system. We study a few cases in detail. Finally, we find the information channel capacity for each of those radio bridges, putting thus a physical constraint to the amount of information transfer that will be possible even by exploiting the stars as gravitational lenses


title: S-matrix bootstrap in 3+1 dimensions: regularization and dual convex problem

abstract: The S-matrix bootstrap maps out the space of S-matrices allowed by analyticity, crossing, unitarity, and other constraints. For the $2\rightarrow 2$ scattering matrix $S_{2\rightarrow 2}$ such space is an infinite dimensional convex space whose boundary can be determined by maximizing linear functionals. On the boundary interesting theories can be found, many times at vertices of the space. Here we consider $3+1$ dimensional theories and focus on the equivalent dual convex minimization problem that provides strict upper bounds for the regularized primal problem and has interesting practical and physical advantages over the primal problem. Its variables are dual partial waves $k_\ell(s)$ that are free variables, namely they do not have to obey any crossing, unitarity or other constraints. Nevertheless they are directly related to the partial waves $f_\ell(s)$, for which all crossing, unitarity and symmetry properties result from the minimization. Numerically, it requires only a few dual partial waves, much as one wants to possibly match experimental results. We consider the case of scalar fields which is related to pion physics.


title: Occupancy-Driven Stochastic Decision Framework for Ranking Commercial Building Loads

abstract: For effective integration of building operations into the evolving demand response programs of the power grid, real-time decisions concerning the use of building appliances for grid services must excel on multiple criteria, ranging from the added value to occupants' comfort to the quality of the grid services. In this paper, we present a data-driven decision-support framework to dynamically rank load control alternatives in a commercial building, addressing the needs of multiple decision criteria (e.g. occupant comfort, grid service quality) under uncertainties in occupancy patterns. We adopt a stochastic multi-criteria decision algorithm recently applied to prioritize residential on/off loads, and extend it to i) complex load control decisions (e.g. dimming of lights, changing zone temperature set-points) in a commercial building; and ii) systematic integration of zonal occupancy patterns to accurately identify short-term grid service opportunities. We evaluate the performance of the framework for curtailment of air-conditioning, lighting, and plug-loads in a multi-zone office building for a range of design choices. With the help of a prototype system that integrates an interactive \textit{Data Analytics and Visualization} frontend, we demonstrate a way for the building operators to monitor the flexibility in energy consumption and to develop trust in the decision recommendations by interpreting the rationale behind the ranking.


title: Experiences on Managing Technical Debt with Code Smells and AntiPatterns

abstract: Technical debt has become a common metaphor for the accumulation of software design and implementation choices that seek fast initial gains but that are under par and counterproductive in the long run. However, as a metaphor, technical debt does not offer actionable advice on how to get rid of it. To get to a practical level in solving problems, more focused mechanisms are needed. Commonly used approaches for this include identifying code smells as quick indications of possible problems in the codebase and detecting the presence of AntiPatterns that refer to overt, recurring problems in design. There are known remedies for both code smells and AntiPatterns. In paper, our goal is to show how to effectively use common tools and the existing body of knowledge on code smells and AntiPatterns to detect technical debt and pay it back. We present two main results: (i) How a combination of static code analysis and manual inspection was used to detect code smells in a codebase leading to the discovery of AntiPatterns; and (ii) How AntiPatterns were used to identify, characterize, and fix problems in the software. The experiences stem from a private company and its long-lasting software product development effort.


title: On the (in)stability of sunspots

abstract: The stability of sunspots is one of the long-standing unsolved puzzles in the field of solar magnetism. We study the effects that destabilise and stabilise the flux tube of a simulated sunspot in the upper convection zone. The depth-varying effects of fluting instability, buoyancy forces, and timescales on the flux tube are analysed. The simulation was calculated with the MURaM code. The domain has a lateral extension of 98 Mm x 98 Mm and extends almost 18 Mm below the solar surface. The analysed data set of 30 hours shows a stable sunspot at the solar surface. We studied the evolution of the flux tube at horizontal layers by means of the relative change in perimeter and area with a linear stability analysis. We find a corrugation along the perimeter of the flux tube that proceeds fastest at a depth of about 8 Mm below the surface. Towards the surface and towards deeper layers, the decrease in compactness is damped. From the stability analysis, we find that above a depth of 2 Mm, the sunspot is stabilised by buoyancy forces. The spot is least stable at a depth of about 3 Mm because of fluting instability. The stability of the sunspot at the surface affects the behaviour of the field lines in deeper layers by magnetic tension. Therefore the fluting instability is damped at depths of about 3 Mm, and the decrease in compactness is strongest at a depth of about 8 Mm. The more vertical orientation of the magnetic field and the longer convective timescale slow down the corrugation process in layers deeper than 10 Mm. The formation of large intrusions of field-free plasma below the surface destabilises the sunspot, and eventually lead to the disruption and decay of the sunspot. This process is not visible at the surface, where the sunspot is stabilised by buoyancy forces. The onset of sunspot decay occurs in deeper layers, while the sunspot still appears stable in the photosphere.


title: The Discovery of Dynamics via Linear Multistep Methods and Deep Learning: Error Estimation

abstract: Identifying hidden dynamics from observed data is a significant and challenging task in a wide range of applications. Recently, the combination of linear multistep methods (LMMs) and deep learning has been successfully employed to discover dynamics, whereas a complete convergence analysis of this approach is still under development. In this work, we consider the deep network-based LMMs for the discovery of dynamics. We put forward error estimates for these methods using the approximation property of deep networks. It indicates, for certain families of LMMs, that the $\ell^2$ grid error is bounded by the sum of $O(h^p)$ and the network approximation error, where $h$ is the time step size and $p$ is the local truncation error order. Numerical results of several physically relevant examples are provided to demonstrate our theory.


title: UCB-based Algorithms for Multinomial Logistic Regression Bandits

abstract: Out of the rich family of generalized linear bandits, perhaps the most well studied ones are logisitc bandits that are used in problems with binary rewards: for instance, when the learner/agent tries to maximize the profit over a user that can select one of two possible outcomes (e.g., `click' vs `no-click'). Despite remarkable recent progress and improved algorithms for logistic bandits, existing works do not address practical situations where the number of outcomes that can be selected by the user is larger than two (e.g., `click', `show me later', `never show again', `no click'). In this paper, we study such an extension. We use multinomial logit (MNL) to model the probability of each one of $K+1\geq 2$ possible outcomes (+1 stands for the `not click' outcome): we assume that for a learner's action $\mathbf{x}_t$, the user selects one of $K+1\geq 2$ outcomes, say outcome $i$, with a multinomial logit (MNL) probabilistic model with corresponding unknown parameter $\bar{\boldsymbol\theta}_{\ast i}$. Each outcome $i$ is also associated with a revenue parameter $\rho_i$ and the goal is to maximize the expected revenue. For this problem, we present MNL-UCB, an upper confidence bound (UCB)-based algorithm, that achieves regret $\tilde{\mathcal{O}}(dK\sqrt{T})$ with small dependency on problem-dependent constants that can otherwise be arbitrarily large and lead to loose regret bounds. We present numerical simulations that corroborate our theoretical results.


title: Investigating the electronic origins of the repulsion between substitutional and interstitial solutes in hcp Ti

abstract: The high solubility of oxygen in Ti, Zr and Hf makes it difficult to stabilize the protective oxide scales on their surfaces as the subsurface regions can serve as boundless sinks that continuously dissolve oxygen. Alloying elements are crucial to reduce the oxygen solubility and diffusivity within early transition metals. Past studies have shown that substitutional alloying additions to titanium repel interstitial oxygen. Here we use first-principles calculations to show that this repulsion is short ranged and identify a variety of factors that are likely responsible for the repulsive interaction. We identify a unique hybridization phenomenon between dissolved substitutional elements and interstitial oxygen within hcp Ti that leads to a repulsive interaction at short distances, similar to that between closed-shell atoms. Calculations of Bader charges also suggest the existence of short-range Coulomb interactions due to the accumulation of charge on the substitutional solute and interstitial oxygen that is drawn from the Ti host.


title: Potential Gap: Using Reactive Policies to Guarantee Safe Navigation

abstract: This paper considers the integration of gap-based local navigation methods with artificial potential field (APF) methods to derive a local planning module for hierarchical navigation systems that has provable collision-free properties. Given that APF theory applies to idealized robot models, the provable properties are lost when applied to more realistic models. We describe a set of algorithm modifications that correct for these errors and enhance robustness to non-ideal models. Central to the construction of the local planner is the use of sensory-derived local free-space models that detect gaps and use them for the synthesis of the APF. Modifications are given for a nonholonomic robot model. Integration of the local planner, called potential gap, into a hierarchical navigation system provides the local goals and trajectories needed for collision-free navigation through unknown environments. Monte Carlo experiments in benchmark worlds confirm the asserted safety and robustness properties by testing under various robot models.


title: Kolyvagin classes versus non-cristalline diagonal classes

abstract: Let $E/\mathbb{Q}$ be an elliptic curve having multiplicative reduction at a prime $p$. Let $(g,h)$ be a pair of eigenforms of weight $1$ arising as the theta series of an imaginary quadratic field $K$, and assume that the triple-product $L$-function $L(f,g,h,s)$ is self-dual and does not vanish at the central critical point $s=1$. The main result of this article is a formula expressing the $p$-adic iterated integrals introduced in [DLR] to the Kolyvagin classes associated by Bertolini and Darmon to a system of Heegner points on $E$.


title: On Efficient Design of Pilot Experiment for Generalized Linear Models

abstract: The experimental design for a generalized linear model (GLM) is important but challenging since the design criterion often depends on model specification including the link function, the linear predictor, and the unknown regression coefficients. Prior to constructing locally or globally optimal designs, a pilot experiment is usually conducted to provide some insights on the model specifications. In pilot experiments, little information on the model specification of GLM is available. Surprisingly, there is very limited research on the design of pilot experiments for GLMs. In this work, we obtain some theoretical understanding of the design efficiency in pilot experiments for GLMs. Guided by the theory, we propose to adopt a low-discrepancy design with respect to some target distribution for pilot experiments. The performance of the proposed design is assessed through several numerical examples.


title: Production of GeV-scale heavy neutral leptons in 3-body decays. Comparison with the PYTHIA approach

abstract: Despite the undeniable success of the Standard Model of particle physics (SM), there are some phenomena that the SM can't explain. These phenomena indicate that the SM has to be modified. One of the possible ways to extend the SM is to introduce heavy neutral leptons (HNLs). To search for HNLs in intensity frontier experiments, one has to consider HNL production both in 2-body and 3-body decays of some mesons. We verified the possibility of using the parton level PYTHIA default matrix elements (without the form-factor formalism) to calculate HNL production in 3-body semileptonic decays of $B$ and $D$ mesons in the experimentally interesting mass range of the produced HNLs. We conclude that this approach is quite suitable for the estimation of the sensitivity region for HNLs in the intensity frontier experiments, provided one uses suitable parton level PYTHIA default matrix elements. Our study was driven by the usage of such an approximation by the SHiP collaboration. We conclude that in this case the parton level PYTHIA default matrix elements could have been chosen more appropriately.


title: Light curve properties of SN 2017fgc and HV SNe Ia

abstract: Photometric and spectroscopic observations of type Ia supernova (SN) 2017fgc which cover the period from $-$12 to +137 days since the $B$-band maximum are presented. SN 2017fgc is a photometrically normal SN Ia with the luminosity decline rate, $ \Delta m_{15} (B)_{true} $= 1.10 $ \pm $ 0.10 mag. Spectroscopically, it belongs to the High Velocity (HV) SNe Ia group, with the Si II $\lambda$6355 velocity near the $B$-band maximum estimated to be 15,200 $ \pm $ 480 km $s^{-1}$. At the epochs around the near-infrared secondary peak, the $R$ and $I$ bands show an excess of $\sim$0.2 mag level compared to the light curves of the normal velocity (NV) SNe Ia. Further inspection of the samples of HV and NV SNe Ia indicates that the excess is a generic feature among HV SNe Ia, different from NV SNe Ia. There is also a hint that the excess is seen in the V band, both in SN 2017fgc and other HV SNe Ia, which behaves like a less prominent shoulder in the light curve. The excess is not obvious in the B band (and unknown in the U band), and the color is consistent with the fiducial SN color. This might indicate the excess is attributed to the bolometric luminosity, not in the color. This excess is less likely caused by external effects, like an echo or change in reddening but could be due to an ionization effect, which reflects an intrinsic, either distinct or continuous, difference in the ejecta properties between HV and NV SNe Ia.


title: Effective protection of quantum coherence by non-Hermitian driving

abstract: In this letter, we investigate the effects of non-Hermitian driving on quantum coherence in a bipartite system. The results that the dynamical localization destroyed by the Hermitian interaction revives are an evidence of the restoration of quantum coherence by non-Hermitian driving. Besides, the entanglement between the two subsystems also decays with the boosting of non-hermitian driving strength, which provides another evidence that non-Hermitian driving will protect quantum coherence. The physics behind this phenomenon is the domination of the quasieigenstate with maximum imaginary value of the quasieigenvalue on the dynamics of the non-Hermitian system. Our discovery establishes a restoration mechanism of quantum coherence in interacting and dissipative quantum systems, which is highly relevant to experiments in diverse fields from many-body physics to quantum information.


title: Dynamical Phase Transitions in models of Collective Neutrino Oscillations

abstract: Collective neutrino oscillations can potentially play an important role in transporting lepton flavor in astrophysical scenarios where the neutrino density is large, typical examples are the early universe and supernova explosions. It has been argued in the past that simple models of the neutrino Hamiltonian designed to describe forward scattering can support substantial flavor evolution on very short time scales $t\approx\log(N)/(G_F\rho_\nu)$, with $N$ the number of neutrinos, $G_F$ the Fermi constant and $\rho_\nu$ the neutrino density. This finding is in tension with results for similar but exactly solvable models for which $t\approx\sqrt{N}/(G_F\rho_\nu)$ instead. In this work we provide a coherent explanation of this tension in terms of Dynamical Phase Transitions (DPT) and study the possible impact that a DPT could have in more realistic models of neutrino oscillations and their mean-field approximation.


title: Spinodal-assisted nucleation in the two-dimensional $q-$state Potts model with short-to-long range interactions

abstract: We study homogeneous nucleation in the two-dimensional $q-$state Potts model for $q=3,5,10,20$ and ferromagnetic couplings $J_{ij} \propto \Theta (R - |i-j|)$, by means of Monte Carlo simulations employing heat bath dynamics. Metastability is induced in the low temperature phase through an instantaneous quench of the magnetic field coupled to one of the $q$ spin states. The quench depth is adjusted, depending on the value of temperature $T$, interaction range $R$, and number of states $q$, in such a way that a constant nucleation time is always obtained. In this setup we analyze the crossover between the classical compact droplet regime occurring in presence of short range interactions $R \sim 1$, and the long-range regime $R\gg 1$ where the properties of nucleation are influenced by the presence of a mean-field spinodal singularity. We evaluate the metastable susceptibility of the order parameter as well as various critical droplet properties, which along with the evolution of the quench depth as a function of $q,T$ and $R$, are then compared with the field theoretical predictions valid in the large $R$ limit in order to find the onset of spinodal-assisted nucleation. We find that, with a mild dependence on the values of $q$ and $T$ considered, spinodal scaling holds for interaction ranges $R\gtrsim 8-10$, and that signatures of the presence of a pseudo-spinodal are already visible for remarkably small interaction ranges $R\sim 4-5$. The influence of spinodal singularities on the occurrence of multi-step nucleation is also discussed.


title: Sum of squares generalizations for conic sets

abstract: In polynomial optimization problems, nonnegativity constraints are typically handled using the sum of squares condition. This can be efficiently enforced using semidefinite programming formulations, or as more recently proposed by Papp and Yildiz [18], using the sum of squares cone directly in a nonsymmetric interior point algorithm. Beyond nonnegativity, more complicated polynomial constraints (in particular, generalizations of the positive semidefinite, second order and $\ell_1$-norm cones) can also be modeled through structured sum of squares programs. We take a different approach and propose using more specialized polynomial cones instead. This can result in lower dimensional formulations, more efficient oracles for interior point methods, or self-concordant barriers with smaller parameters. In most cases, these algorithmic advantages also translate to faster solving times in practice.


title: Sinusoidal Parameter Estimation from Signed Measurements via Majorization-Minimization Based RELAX

abstract: We consider the problem of sinusoidal parameter estimation using signed observations obtained via one-bit sampling with fixed as well as time-varying thresholds. In a previous paper, a relaxation-based algorithm, referred to as 1bRELAX, has been proposed to iteratively maximize the likelihood function. However, the exhaustive search procedure used in each iteration of 1bRELAX is time-consuming. In this paper, we present a majorization-minimization (MM) based 1bRELAX algorithm, referred to as 1bMMRELAX, to enhance the computational efficiency of 1bRELAX. Using the MM technique, 1bMMRELAX maximizes the likelihood function iteratively using simple FFT operations instead of the more computationally intensive search used by 1bRELAX. Both simulated and experimental results are presented to show that 1bMMRELAX can significantly reduce the computational cost of 1bRELAX while maintaining its excellent estimation accuracy.


title: Quantum Phases of a Weakly Disordered Josephson Ladder

abstract: The interplay of interactions and disorder in low-dimensional superconductors supports the formation of multiple quantum phases, as possible instabilities of the Superconductor-Insulator Transition (SIT) at a singular quantum critical point. We explore a one-dimensional model which exhibits such variety of phases in the strongly quantum fluctuations regime. Specifically, we study the effect of weak disorder on a two-leg Josephson ladder with comparable Josephson and charging energies ($E_J$~$E_C$). An additional key feature of our model is the requirement of perfect $\mathbb{Z}_2$-symmetry, respected by all parameters including the disorder. Using a perturbative renormalization-group (RG) analysis, we derive the phase diagram and identify at least one intermediate phase between a full-fledged superconductor and a disorder-dominated insulator. Most prominently, for repulsive interactions on the rungs we identify two distinct mixed phases: in both of them the longitudinal charge mode is a gapless superconductor, however one phase exhibits a dipolar charge density order on the rungs, while the other is disordered. This latter phase is characterized by coexisting superconducting (phase-locked) and charge-ordered rungs, and encompasses the potential of evolving into a Griffith's phase characteristic of the random-field Ising model in the strong disorder limit.


title: Charged Throats in the Ho\v{r}ava--Lifshitz Theory

abstract: A spherically symmetric solution of the field equations of the Ho\v{r}ava--Lifshitz gravity--gauge vector interaction theory is obtained and analyzed. It describes a charged throat. The solution exists provided a restriction on the relation between the mass and charge is satisfied. The restriction reduces to the Reissner--Nordstr\"{o}m one in the limit in which the coupling constants tend to the relativistic values of General Relativity. We introduce the correct charts to describe the solution across the entire manifold, including the throat connecting an asymptotic Minkowski space--time with a singular 3+1 dimensional manifold. The solution external to the throat on the asymptotically flat side tends to the Reissner--Nordstr\"{o}m space--time at the limit when the coupling parameter, associated with the term in the low energy Hamiltonian that manifestly breaks the relativistic symmetry, tends to zero. Also, when the electric charge is taken to be zero the solution becomes the spherically symmetric and static solution of the Ho\v{r}ava--Lifshitz gravity.


title: Formation and fate of the born-again planetary nebula HuBi 1

abstract: We present the first 3D radiation-hydrodynamic simulations on the formation and evolution of born-again planetary nebulae (PNe), with particular emphasis to the case of HuBi1, the inside-out PN. We use the extensively-tested GUACHO code to simulate the formation of HuBi1 adopting mass-loss and stellar wind terminal velocity estimates obtained from observations presented by our group. We found that, if the inner shell of HuBi1 was formed by an explosive very late thermal pulse (VLTP) ejecting material with velocities of $\sim$300 km s$^{-1}$, the age of this structure is consistent with that of $\simeq$200 yr derived from multi-epoch narrow-band imaging. Our simulations predict that, as a consequence of the dramatic reduction of the stellar wind velocity and photon ionizing flux during the VLTP, the velocity and pressure structure of the outer H-rich nebula are affected creating turbulent ionized structures surrounding the inner shell. These are indeed detected in Gran Telescopio Canarias MEGARA optical observations. Furthermore, we demonstrate that the current relatively low ionizing photon flux from the central star of HuBi1 is not able to completely ionize the inner shell, which favors previous suggestions that its excitation is dominated by shocks. Our simulations suggest that the kinetic energy of the H-poor ejecta of HuBi1 is at least 30 times that of the clumps and filaments in the evolved born-again PNe A30 and A78, making it a truly unique VLTP event.


title: Purchase history and product personalization

abstract: Product personalization opens the door to price discrimination. A rich product line allows firms to better tailor products to consumers' tastes, but the mere choice of a product carries valuable information about consumers that can be leveraged for price discrimination. We study this trade-off in an upstream-downstream model, where a consumer buys a good of variable quality upstream, followed by an indivisible good downstream. The downstream firm's use of the consumer's purchase history for price discrimination introduces a novel distortion: The upstream firm offers a subset of the products that it would offer if, instead, it could jointly design its product line and downstream pricing. By controlling the degree of product personalization the upstream firm curbs ratcheting forces that result from the consumer facing downstream price discrimination.


title: Policy-Guided Heuristic Search with Guarantees

abstract: The use of a policy and a heuristic function for guiding search can be quite effective in adversarial problems, as demonstrated by AlphaGo and its successors, which are based on the PUCT search algorithm. While PUCT can also be used to solve single-agent deterministic problems, it lacks guarantees on its search effort and it can be computationally inefficient in practice. Combining the A* algorithm with a learned heuristic function tends to work better in these domains, but A* and its variants do not use a policy. Moreover, the purpose of using A* is to find solutions of minimum cost, while we seek instead to minimize the search loss (e.g., the number of search steps). LevinTS is guided by a policy and provides guarantees on the number of search steps that relate to the quality of the policy, but it does not make use of a heuristic function. In this work we introduce Policy-guided Heuristic Search (PHS), a novel search algorithm that uses both a heuristic function and a policy and has theoretical guarantees on the search loss that relates to both the quality of the heuristic and of the policy. We show empirically on the sliding-tile puzzle, Sokoban, and a puzzle from the commercial game `The Witness' that PHS enables the rapid learning of both a policy and a heuristic function and compares favorably with A*, Weighted A*, Greedy Best-First Search, LevinTS, and PUCT in terms of number of problems solved and search time in all three domains tested.


title: QCD Analysis of Near-Threshold Photon-Proton Production of Heavy Quarkonium

abstract: The near threshold photo or electroproduction of heavy vector quarkonium off the proton is studied in quantum chromodynamics. Similar to the high-energy limit, the production amplitude can be factorized in terms of gluonic Generalized Parton Distributions and the quarkonium distribution amplitude. At the threshold, the threshold kinematics has a large skewness parameter $\xi$, leading to the dominance of the spin-2 contribution over higher-spin twist-2 operators. Thus threshold production data are useful to extract the gluonic gravitational form factors, allowing studying the gluonic contributions to the quantum anomalous energy, mass radius, spin and mechanical pressure in the proton. We use the recent GlueX data on the $J/\psi$ photoproduction to illustrate the potential physics impact from the high-precision data from future JLab 12 GeV and EIC physics program.


title: Intrinsic tame filling functions are equivalent to intrinsic diameter functions

abstract: Intrinsic tame filling functions are quasi-isometry invariants that are refinements of the intrinsic diameter function of a group. The main purpose of this paper is to show that every finite presentation of a group has an intrinsic tame filling function that is equivalent to its intrinsic diameter function. We also introduce some alternative filling functions--based on concepts similar to those used to define intrinsic tame filling functions--that are potential proper refinements of the intrinsic diameter function.


title: The G\'alvez-Kock-Tonks conjecture for locally discrete decomposition spaces

abstract: G\'alvez-Carrillo, Kock, and Tonks constructed a decomposition space $U$ of all M\"obius intervals, as a recipient of Lawvere's interval construction for M\"obius categories, and conjectured that $U$ enjoys a certain universal property: for every M\"obius decomposition space $X$, the space of culf functors from $X$ to $U$ is contractible. In this paper, we work at the level of homotopy 1-types to prove the first case of the conjecture, namely for locally discrete decomposition spaces. This provides also the first substantial evidence for the general conjecture. This case is general enough to cover all locally finite posets, Cartier--Foata monoids, M\"obius categories and strict (directed) restriction species. The proof is 2-categorical. First, we construct a local strict model of $U$, which is then used to show by hand that the Lawvere interval construction, considered as a natural transformation, does not admit other self-modifications than the identity.


title: Scatter Correction in X-ray CT by Physics-Inspired Deep Learning

abstract: Scatter due to interaction of photons with the imaged object is a fundamental problem in X-ray Computed Tomography (CT). It manifests as various artifacts in the reconstruction, making its abatement or correction critical for image quality. Despite success in specific settings, hardware-based methods require modification in the hardware, or increase in the scan time or dose. This accounts for the great interest in software-based methods, including Monte-Carlo based scatter estimation, analytical-numerical, and kernel-based methods, with data-driven learning-based approaches demonstrated recently. In this work, two novel physics-inspired deep-learning-based methods, PhILSCAT and OV-PhILSCAT, are proposed. The methods estimate and correct for the scatter in the acquired projection measurements. Different from previous works, they incorporate both an initial reconstruction of the object of interest and the scatter-corrupted measurements related to it, and use a deep neural network architecture and cost function, both specifically tailored to the problem. Numerical experiments with data generated by Monte-Carlo simulations of the imaging of phantoms reveal consistent improvement over a recent purely projection-domain deep neural network scatter correction method.


title: A combinatorial realization of Kirillov-Reshetikhin crystals for type E arising from translations

abstract: The main purpose of this paper is to give a combinatorial realization of Kirillov-Reshetikhin (KR simply) crystals $B^{r, s}$ for type $\text{E}_n^{(1)}$ with a minuscule node $r$ and $s \ge 1$. To do this, we describe explicitly the crystal of the quantum nilpotent subalgebra associated with the translation by the negative of the $r$-th fundamental weight. Then the crystal can be extended as an affine crystal, in which a certain subcrystal characterized by the $\varepsilon_r^*$-statistic is isomorphic to $B^{r,s}$ as an affine crystal, where $\varepsilon_r^*$ is also realized precisely in terms of triple and quadruple paths.


title: MoViNets: Mobile Video Networks for Efficient Video Recognition

abstract: We present Mobile Video Networks (MoViNets), a family of computation and memory efficient video networks that can operate on streaming video for online inference. 3D convolutional neural networks (CNNs) are accurate at video recognition but require large computation and memory budgets and do not support online inference, making them difficult to work on mobile devices. We propose a three-step approach to improve computational efficiency while substantially reducing the peak memory usage of 3D CNNs. First, we design a video network search space and employ neural architecture search to generate efficient and diverse 3D CNN architectures. Second, we introduce the Stream Buffer technique that decouples memory from video clip duration, allowing 3D CNNs to embed arbitrary-length streaming video sequences for both training and inference with a small constant memory footprint. Third, we propose a simple ensembling technique to improve accuracy further without sacrificing efficiency. These three progressive techniques allow MoViNets to achieve state-of-the-art accuracy and efficiency on the Kinetics, Moments in Time, and Charades video action recognition datasets. For instance, MoViNet-A5-Stream achieves the same accuracy as X3D-XL on Kinetics 600 while requiring 80% fewer FLOPs and 65% less memory. Code will be made available at https://github.com/tensorflow/models/tree/master/official/vision.


title: Robust Multi-Modal Policies for Industrial Assembly via Reinforcement Learning and Demonstrations: A Large-Scale Study

abstract: Over the past several years there has been a considerable research investment into learning-based approaches to industrial assembly, but despite significant progress these techniques have yet to be adopted by industry. We argue that it is the prohibitively large design space for Deep Reinforcement Learning (DRL), rather than algorithmic limitations per se, that are truly responsible for this lack of adoption. Pushing these techniques into the industrial mainstream requires an industry-oriented paradigm which differs significantly from the academic mindset. In this paper we define criteria for industry-oriented DRL, and perform a thorough comparison according to these criteria of one family of learning approaches, DRL from demonstration, against a professional industrial integrator on the recently established NIST assembly benchmark. We explain the design choices, representing several years of investigation, which enabled our DRL system to consistently outperform the integrator baseline in terms of both speed and reliability. Finally, we conclude with a competition between our DRL system and a human on a challenge task of insertion into a randomly moving target. This study suggests that DRL is capable of outperforming not only established engineered approaches, but the human motor system as well, and that there remains significant room for improvement. Videos can be found on our project website: https://sites.google.com/view/shield-nist.


title: Maximal $L^1$ regularity for solutions to inhomogeneous incompressible Navier-Stokes equations

abstract: This paper is devoted to the maximal $L^1$ regularity and asymptotic behavior for solutions to the inhomogeneous incompressible Navier-Stokes equations under a scaling-invariant smallness assumption on the initial velocity. We obtain a new global $L^1$-in-time estimate for the Lipschitz seminorm of the velocity field without any smallness assumption on the initial density fluctuation. In the derivation of this estimate, we study the maximal $L^1$ regularity for a linear Stokes system with variable coefficients. The analysis tools are a use of the semigroup generated by a generalized Stokes operator to characterize some Besov norms and a new gradient estimate for a class of second-order elliptic equations of divergence form. Our method might be used to study some other issues arising from incompressible or compressible viscous fluids.


title: Higher Siegel--Weil formula for unitary groups: the non-singular terms

abstract: We construct special cycles on the moduli stack of unitary shtukas. We prove an identity between (1) the r-th central derivative of non-singular Fourier coefficients of a normalized Siegel--Eisenstein series, and (2) the degree of special cycles of "virtual dimension 0" on the moduli stack of unitary shtukas with r legs. This may be viewed as a function-field analogue of the Kudla-Rapoport Conjecture, that has the additional feature of encompassing all higher derivatives of the Eisenstein series.


title: A curated dataset for data-driven turbulence modelling

abstract: The recent surge in machine learning augmented turbulence modelling is a promising approach for addressing the limitations of Reynolds-averaged Navier-Stokes (RANS) models. This work presents the development of the first open-source dataset, curated and structured for immediate use in machine learning augmented turbulence closure modelling. The dataset features a variety of RANS simulations with matching direct numerical simulation (DNS) and large-eddy simulation (LES) data. Four turbulence models are selected to form the initial dataset: $k$-$\varepsilon$, $k$-$\varepsilon$-$\phi_t$-$f$, $k$-$\omega$, and $k$-$\omega$ SST. The dataset consists of 29 cases per turbulence model, for several parametrically sweeping reference DNS/LES cases: periodic hills, square duct, parametric bumps, converging-diverging channel, and a curved backward-facing step. At each of the 895,640 points, various RANS features with DNS/LES labels are available. The feature set includes quantities used in current state-of-the-art models, and additional fields which enable the generation of new feature sets. The dataset reduces effort required to train, test, and benchmark new models. The dataset is available at https://doi.org/10.34740/kaggle/dsv/2044393 .


title: Homophily Outlier Detection in Non-IID Categorical Data

abstract: Most of existing outlier detection methods assume that the outlier factors (i.e., outlierness scoring measures) of data entities (e.g., feature values and data objects) are Independent and Identically Distributed (IID). This assumption does not hold in real-world applications where the outlierness of different entities is dependent on each other and/or taken from different probability distributions (non-IID). This may lead to the failure of detecting important outliers that are too subtle to be identified without considering the non-IID nature. The issue is even intensified in more challenging contexts, e.g., high-dimensional data with many noisy features. This work introduces a novel outlier detection framework and its two instances to identify outliers in categorical data by capturing non-IID outlier factors. Our approach first defines and incorporates distribution-sensitive outlier factors and their interdependence into a value-value graph-based representation. It then models an outlierness propagation process in the value graph to learn the outlierness of feature values. The learned value outlierness allows for either direct outlier detection or outlying feature selection. The graph representation and mining approach is employed here to well capture the rich non-IID characteristics. Our empirical results on 15 real-world data sets with different levels of data complexities show that (i) the proposed outlier detection methods significantly outperform five state-of-the-art methods at the 95%/99% confidence level, achieving 10%-28% AUC improvement on the 10 most complex data sets; and (ii) the proposed feature selection methods significantly outperform three competing methods in enabling subsequent outlier detection of two different existing detectors.


title: Dual Monte Carlo Tree Search

abstract: AlphaZero, using a combination of Deep Neural Networks and Monte Carlo Tree Search (MCTS), has successfully trained reinforcement learning agents in a tabula-rasa way. The neural MCTS algorithm has been successful in finding near-optimal strategies for games through self-play. However, the AlphaZero algorithm has a significant drawback; it takes a long time to converge and requires high computational power due to complex neural networks for solving games like Chess, Go, Shogi, etc. Owing to this, it is very difficult to pursue neural MCTS research without cutting-edge hardware, which is a roadblock for many aspiring neural MCTS researchers. In this paper, we propose a new neural MCTS algorithm, called Dual MCTS, which helps overcome these drawbacks. Dual MCTS uses two different search trees, a single deep neural network, and a new update technique for the search trees using a combination of the PUCB, a sliding-window, and the epsilon-greedy algorithm. This technique is applicable to any MCTS based algorithm to reduce the number of updates to the tree. We show that Dual MCTS performs better than one of the most widely used neural MCTS algorithms, AlphaZero, for various symmetric and asymmetric games.


title: An Empirical Study of OSS-Fuzz Bugs

abstract: Continuous fuzzing is an increasingly popular technique for automated quality and security assurance. Google maintains OSS-Fuzz: a continuous fuzzing service for open source software. We conduct the first empirical study of OSS-Fuzz, analyzing 23,907 bugs found in 316 projects. We examine the characteristics of fuzzer-found faults, the lifecycles of such faults, and the evolution of fuzzing campaigns over time. We find that OSS-Fuzz is often effective at quickly finding bugs, and developers are often quick to patch them. However, flaky bugs, timeouts, and out of memory errors are problematic, people rarely file CVEs for security vulnerabilities, and fuzzing campaigns often exhibit punctuated equilibria, where developers might be surprised by large spikes in bugs found. Our findings have implications on future fuzzing research and practice.


title: Towards Improving the Trustworthiness of Hardware based Malware Detector using Online Uncertainty Estimation

abstract: Hardware-based Malware Detectors (HMDs) using Machine Learning (ML) models have shown promise in detecting malicious workloads. However, the conventional black-box based machine learning (ML) approach used in these HMDs fail to address the uncertain predictions, including those made on zero-day malware. The ML models used in HMDs are agnostic to the uncertainty that determines whether the model "knows what it knows," severely undermining its trustworthiness. We propose an ensemble-based approach that quantifies uncertainty in predictions made by ML models of an HMD, when it encounters an unknown workload than the ones it was trained on. We test our approach on two different HMDs that have been proposed in the literature. We show that the proposed uncertainty estimator can detect >90% of unknown workloads for the Power-management based HMD, and conclude that the overlapping benign and malware classes undermine the trustworthiness of the Performance Counter-based HMD.


title: Unsupervised and self-adaptative techniques for cross-domain person re-identification

abstract: Person Re-Identification (ReID) across non-overlapping cameras is a challenging task and, for this reason, most works in the prior art rely on supervised feature learning from a labeled dataset to match the same person in different views. However, it demands the time-consuming task of labeling the acquired data, prohibiting its fast deployment, specially in forensic scenarios. Unsupervised Domain Adaptation (UDA) emerges as a promising alternative, as it performs feature-learning adaptation from a model trained on a source to a target domain without identity-label annotation. However, most UDA-based algorithms rely upon a complex loss function with several hyper-parameters, which hinders the generalization to different scenarios. Moreover, as UDA depends on the translation between domains, it is important to select the most reliable data from the unseen domain, thus avoiding error propagation caused by noisy examples on the target data -- an often overlooked problem. In this sense, we propose a novel UDA-based ReID method that optimizes a simple loss function with only one hyper-parameter and that takes advantage of triplets of samples created by a new offline strategy based on the diversity of cameras within a cluster. This new strategy adapts the model and also regularizes it, avoiding overfitting on the target domain. We also introduce a new self-ensembling strategy, in which weights from different iterations are aggregated to create a final model combining knowledge from distinct moments of the adaptation. For evaluation, we consider three well-known deep learning architectures and combine them for final decision-making. The proposed method does not use person re-ranking nor any label on the target domain, and outperforms the state of the art, with a much simpler setup, on the Market to Duke, the challenging Market1501 to MSMT17, and Duke to MSMT17 adaptation scenarios.


title: Conditional Frechet Inception Distance

abstract: We consider distance functions between conditional distributions. We focus on the Wasserstein metric and its Gaussian case known as the Frechet Inception Distance (FID). We develop conditional versions of these metrics, analyze their relations and provide a closed form solution to the conditional FID (CFID) metric. We numerically compare the metrics in the context of performance evaluation of modern conditional generative models. Our results show the advantages of CFID compared to the classical FID and mean squared error (MSE) measures. In contrast to FID, CFID is useful in identifying failures where realistic outputs which are not related to their inputs are generated. On the other hand, compared to MSE, CFID is useful in identifying failures where a single realistic output is generated even though there is a diverse set of equally probable outputs.


title: Multi-directional Bicycle Robot for Steel Structure Inspection

abstract: This paper presents a novel design of a multi-directional bicycle robot, which targets inspecting general ferromagnetic structures including complex-shaped structures. The locomotion concept is based on arranging two magnetic wheels in a bicycle-like configuration with two independent steering actuators. This configuration allows the robot to possess multi-directional mobility. An additional free joint helps the robot naturally adapt to non-flat and complex surfaces of steel structures. The robot has the biggest advantage to be mechanically simple with high mobility. Besides, the robot is equipped with sensing tools for structure health monitoring. We demonstrate the deployment of our robot to perform steel rust detection on steel bridges. The final inspection results are visualized as 3D models of the bridges together with marked locations of detected rusty areas.


title: Lipschitz normally embedded set and tangent cones at infinity

abstract: We prove that any analytic set in $\C^n$ with a unique tangent cone at infinity is an algebraic set. We prove that the degree of a complex algebraic set in $\C^n$, which is Lipschitz normally embedded at infinity, is equal to the degree of its tangent cone at infinity.


title: Local Cohomology of Module of Differentials of integral extensions

abstract: The main focus of this paper is on determining the highest non-vanishing local cohomology modules of \Omega_(B/R), \Omega_(B/V)(\Omega_(B/k)) where R is either a complete regular local ring or a complete local normal domain with coefficient ring V (field k) and B is its integral closure in an algebraic extension of Q(R). Similar problem is also studied over a normal domain R containing a field k of characteristic 0. In this connection new observations on the direct summand property for integral extensions are also presented.


title: hep_tables: Heterogeneous Array Programming for HEP

abstract: Array operations are one of the most concise ways of expressing common filtering and simple aggregation operations that is the hallmark of the first step of a particle physics analysis: selection, filtering, basic vector operations, and filling histograms. The High Luminosity run of the Large Hadron Collider (HL-LHC), scheduled to start in 2026, will require physicists to regularly skim datasets that are over a PB in size, and repeatedly run over datasets that are 100's of TB's - too big to fit in memory. Declarative programming techniques are a way of separating the intent of the physicist from the mechanics of finding the data, processing the data, and using distributed computing to process it efficiently that is required to extract the plot or data desired in a timely fashion. This paper describes a prototype library that provides a framework for different sub-systems to cooperate in producing this data, using an array-programming declarative interface. This prototype has a ServiceX data-delivery sub-system and an awkward array sub-system cooperating to generate requested data. The ServiceX system runs against ATLAS xAOD data.


title: ExAD: An Ensemble Approach for Explanation-based Adversarial Detection

abstract: Recent research has shown Deep Neural Networks (DNNs) to be vulnerable to adversarial examples that induce desired misclassifications in the models. Such risks impede the application of machine learning in security-sensitive domains. Several defense methods have been proposed against adversarial attacks to detect adversarial examples at test time or to make machine learning models more robust. However, while existing methods are quite effective under blackbox threat model, where the attacker is not aware of the defense, they are relatively ineffective under whitebox threat model, where the attacker has full knowledge of the defense. In this paper, we propose ExAD, a framework to detect adversarial examples using an ensemble of explanation techniques. Each explanation technique in ExAD produces an explanation map identifying the relevance of input variables for the model's classification. For every class in a dataset, the system includes a detector network, corresponding to each explanation technique, which is trained to distinguish between normal and abnormal explanation maps. At test time, if the explanation map of an input is detected as abnormal by any detector model of the classified class, then we consider the input to be an adversarial example. We evaluate our approach using six state-of-the-art adversarial attacks on three image datasets. Our extensive evaluation shows that our mechanism can effectively detect these attacks under blackbox threat model with limited false-positives. Furthermore, we find that our approach achieves promising results in limiting the success rate of whitebox attacks.


title: Quadratic Gravity and Restricted Weyl Symmetry

abstract: We investigate the relationship between quadratic gravity and a restricted Weyl symmetry where a gauge parameter $\Omega(x)$ of Weyl transformation satisfies a constraint $\Box \Omega = 0$ in a curved space-time. First, we briefly review a model with a restricted gauge symmetry on the basis of QED where a $U(1)$ gauge parameter $\theta(x)$ obeys a similar constraint $\Box \theta = 0$ in a flat Minkowski space-time, and explain that the restricted gauge symmetry removes one on-shell mode of gauge field, which together with the Feynman gauge leaves only two transverse polarizations as physical states. Next, it is shown that the restricted Weyl symmetry also eliminates one component of a dipole field in quadratic gravity around a flat Minkowski background, leaving only a single scalar state. Finally, we show that the restricted Weyl symmetry cannot remove any dynamical degrees of freedom in static background metrics by using the zero-energy theorem of quadratic gravity. This fact also holds for the Euclidean background metrics without imposing the static condition.


title: A Large-scale Dataset for Hate Speech Detection on Vietnamese Social Media Texts

abstract: In recent years, Vietnam witnesses the mass development of social network users on different social platforms such as Facebook, Youtube, Instagram, and Tiktok. On social medias, hate speech has become a critical problem for social network users. To solve this problem, we introduce the ViHSD - a human-annotated dataset for automatically detecting hate speech on the social network. This dataset contains over 30,000 comments, each comment in the dataset has one of three labels: CLEAN, OFFENSIVE, or HATE. Besides, we introduce the data creation process for annotating and evaluating the quality of the dataset. Finally, we evaluated the dataset by deep learning models and transformer models.


title: DCT and DST Filtering with Sparse Graph Operators

abstract: Graph filtering is a fundamental tool in graph signal processing. Polynomial graph filters (PGFs), defined as polynomials of a fundamental graph operator, can be implemented in the vertex domain, and usually have a lower complexity than frequency domain filter implementations. In this paper, we focus on the design of filters for graphs with graph Fourier transform (GFT) corresponding to a discrete trigonometric transform (DTT), i.e., one of 8 types of discrete cosine transforms (DCT) and 8 discrete sine transforms (DST). In this case, we show that multiple sparse graph operators can be identified, which allows us to propose a generalization of PGF design: multivariate polynomial graph filter (MPGF). First, for the widely used DCT-II (type-2 DCT), we characterize a set of sparse graph operators that share the DCT-II matrix as their common eigenvector matrix. This set contains the well-known connected line graph. These sparse operators can be viewed as graph filters operating in the DCT domain, which allows us to approximate any DCT graph filter by a MPGF, leading to a design with more degrees of freedom than the conventional PGF approach. Then, we extend those results to all of the 16 DTTs as well as their 2D versions, and show how their associated sets of multiple graph operators can be determined. We demonstrate experimentally that ideal low-pass and exponential DCT/DST filters can be approximated with higher accuracy with similar runtime complexity. Finally, we apply our method to transform-type selection in a video codec, AV1, where we demonstrate significant encoding time savings, with a negligible compression loss.


title: Defending Quantum Objectivity

abstract: A recent argument, attributed to Masanes, is claimed to show that the assumption that quantum measurements have definite, objective outcomes, is incompatible with quantum predictions. In this work, a detailed examination of the argument shows that it has a much narrower field of application than previously recognized. In particular, it is found: i) that the argument only applies to hidden-variable models with a particular feature; and ii) that such a feature is not present in most hidden-variable models, including pilot-wave theory. It is concluded that the argument does not succeed in calling into question the objectivity of quantum measurements.


title: Lyapunov exponents for random maps

abstract: It has been recently realized that for abundant dynamical systems on a compact manifold, the set of points for which Lyapunov exponents fail to exist, called the Lyapunov irregular set, has positive Lebesgue measure. In the present paper, we show that under any physical noise, the Lyapunov irregular set has zero Lebesgue measure and the number of such Lyapunov exponents is finite. This result is a Lyapunov exponent version of Ara\'{u}jo's theorem on the existence and finitude of time averages. Furthermore, we numerically compute the Lyapunov exponents for a surface flow with an attracting heteroclinic connection, which enjoys the Lyapunov irregular set of positive Lebesgue measure, under a physical noise. This paper also contains the proof of the disappearance of Lyapunov irregular behavior on a positive Lebesgue measure set for a surface flow with an attracting homoclinic/heteroclinic connection under a non-physical noise.


title: Faithfullness of geometric action of skein algebras

abstract: We show that the action of the Kauffman bracket skein algebra of a surface $\Sigma$ on the skein module of the handlebody bounded by $\Sigma$ is faithful if and only if the quantum parameter is not a root of 1.


title: Metastability of Blume-Capel model with zero chemical potential and zero external field

abstract: In this study, we investigate the metastable behavior of Metropolis-type Glauber dynamics associated with the Blume-Capel model with zero chemical potential and zero external field at very low temperatures. The corresponding analyses for the same model with zero chemical potential and positive small external field were performed in [Cirillo and Nardi, Journal of Statistical Physics, 150: 1080-1114, 2013] and [Landim and Lemire, Journal of Statistical Physics, 164: 346-376, 2016]. We obtain both large deviation-type and potential-theoretic results on the metastable behavior in our setting. To this end, we perform highly thorough investigation on the energy landscape, where it is revealed that no critical configurations exist and alternatively a massive flat plateau of saddle configurations resides therein.


title: Leveraging Models to Reduce Test Cases in Software Repositories

abstract: Given a failing test case, test case reduction yields a smaller test case that reproduces the failure. This process can be time consuming due to repeated trial and error with smaller test cases. Current techniques speed up reduction by only exploring syntactically valid candidates, but they still spend significant effort on semantically invalid candidates. In this paper, we propose a model-guided approach to speed up test case reduction. The approach trains a model of semantic properties driven by syntactic test case properties. By using this model, we can skip testing even syntactically valid test case candidates that are unlikely to succeed. We evaluate this model-guided reduction on a suite of 14 large fuzzer-generated C test cases from the bug repositories of two well-known C compilers, GCC and Clang. Our results show that with an average precision of 77%, we can decrease the number of removal trials by 14% to 61%. We observe a 30% geomean improvement in reduction time over the state of the art technique while preserving similar reduction power.


title: On derivatives of Kato's Euler system and the Mazur-Tate Conjecture

abstract: We provide a new interpretation of the Mazur-Tate Conjecture and then use it to obtain the first (unconditional) theoretical evidence in support of the conjecture for elliptic curves of strictly positive rank.


title: General condition of quantum teleportation by one-dimensional quantum walks

abstract: We extend the scheme of quantum teleportation by quantum walks introduced by Wang et al. (2017). First, we introduce the mathematical definition of the accomplishment of quantum teleportation by this extended scheme. Secondly, we show a useful necessary and sufficient condition that the quantum teleportation is accomplished rigorously. Our result classifies the parameters of the setting for the accomplishment of the quantum teleportation.


title: How to Design Sample and Computationally Efficient VQA Models

abstract: In multi-modal reasoning tasks, such as visual question answering (VQA), there have been many modeling and training paradigms tested. Previous models propose different methods for the vision and language tasks, but which ones perform the best while being sample and computationally efficient? Based on our experiments, we find that representing the text as probabilistic programs and images as object-level scene graphs best satisfy these desiderata. We extend existing models to leverage these soft programs and scene graphs to train on question answer pairs in an end-to-end manner. Empirical results demonstrate that this differentiable end-to-end program executor is able to maintain state-of-the-art accuracy while being sample and computationally efficient.


title: An averaging formula for the cohomology of PEL-type Rapoport--Zink spaces

abstract: We prove under certain assumptions a formula for the cohomology of PEL-type Rapoport--Zink spaces that "averages" over the Kottwitz set. Our formula generalizes that of Shin's beyond the EL-type case and is proven by combining Mantovan's formula with descriptions of the cohomology of Shimura and Igusa varieties. We then use this averaging formula to derive a conjectural description of the cohomology of Rapoport--Zink spaces, generalizing our earlier work for EL-type spaces. Along the way, we give a description of the cohomology of Igusa varieties in terms of the Langlands correspondence, generalizing work in Shin's thesis.


title: Interpretable, predictive spatio-temporal models via enhanced Pairwise Directions Estimation

abstract: This article concerns the predictive modeling for spatio-temporal data as well as model interpretation using data information in space and time. We develop a novel approach based on supervised dimension reduction for such data in order to capture nonlinear mean structures without requiring a prespecified parametric model. In addition to prediction as a common interest, this approach emphasizes the exploration of geometric information from the data. The method of Pairwise Directions Estimation (PDE; Lue, 2019) is implemented in our approach as a data-driven function searching for spatial patterns and temporal trends. The benefit of using geometric information from the method of PDE is highlighted, which aids effectively in exploring data structures. We further enhance PDE, referring to it as PDE+, by incorporating kriging to estimate the random effects not explained in the mean functions. Our proposal can not only increase prediction accuracy, but also improve the interpretation for modeling. Two simulation examples are conducted and comparisons are made with four existing methods. The results demonstrate that the proposed PDE+ method is very useful for exploring and interpreting the patterns and trends for spatio-temporal data. Illustrative applications to two real datasets are also presented.


title: Search for the charged lepton flavor violating decay $J/\psi\to e\tau$

abstract: A search for the charged lepton flavor violating decay $J/\psi\to e^{\pm}\tau^{\mp}$ with $\tau^{\mp} \to \pi^{\mp}\pi^0\nu_{\tau}$ is performed with about $10$ billion $J/ \psi$ events collected with the BESIII detector at the BEPCII. No significant signal is observed, and an upper limit is set on the branching fraction $\mathcal{B}(J/\psi\to e^{\pm}\tau^{\mp})<7.5\times10^{-8}$ at the 90$\%$ confidence level. This improves the previously published limit by two orders of magnitude.


title: Fluence Adaptation for Task-based Dose Optimization in X-ray Phase-Contrast Imaging

abstract: Purpose: Grating-based imaging (GBI) and edge-illumination (EI) are two promising types of XPCI as the conventional x-ray sources can be directly utilized. For GBI and EI systems, the phase-stepping acquisition with multiple exposures at a constant fluence is usually adopted in the literature. This work, however, attempts to challenge such a constant fluence concept during the phase-stepping process and proposes a fluence adaptation mechanism for dose reduction. Method: Recently, analytic multi-order moment analysis has been proposed to improve the computing efficiency. In these algorithms, multiple contrasts can be calculated by summing together the weighted phase-stepping curves (PSCs) with some kernel functions, which suggests us that the raw data at different steps have different contributions for the noise in retrieved contrasts. Based on analytic retrieval formulas and the Gaussian noise model for detected signals, we derived an optimal adaptive fluence distribution, which is proportional to the absolute weighting kernel functions and the root of original sample PSCs acquired under the constant fluence. Results: To validate our analyses, simulations and experiments are conducted for GBI and EI systems. Simulated results demonstrate that the dose reduction ratio between our proposed fluence distributions and the typical constant one can be about 20% for the phase contrast, which is consistent with our theoretical predictions. Although the experimental noise reduction ratios are a little smaller than the theoretical ones, synthetic and real experiments both observe better noise performance by our proposed method. Our simulated results also give out the effective ranges of the parameters of the PSCs, such as the visibility in GBI, the standard deviation and the mean value in EI, providing a guidance for the use of our proposed approach in practice.


title: Smart Scheduling based on Deep Reinforcement Learning for Cellular Networks

abstract: To improve the system performance towards the Shannon limit, advanced radio resource management mechanisms play a fundamental role. In particular, scheduling should receive much attention, because it allocates radio resources among different users in terms of their channel conditions and QoS requirements. The difficulties of scheduling algorithms are the tradeoffs need to be made among multiple objectives, such as throughput, fairness and packet drop rate. We propose a smart scheduling scheme based on deep reinforcement learning (DRL). We not only verify the performance gain achieved, but also provide implementation-friend designs, i.e., a scalable neural network design for the agent and a virtual environment training framework. With the scalable neural network design, the DRL agent can easily handle the cases when the number of active users is time-varying without the need to redesign and retrain the DRL agent. Training the DRL agent in a virtual environment offline first and using it as the initial version in the practical usage helps to prevent the system from suffering from performance and robustness degradation due to the time-consuming training. Through both simulations and field tests, we show that the DRL-based smart scheduling outperforms the conventional scheduling method and can be adopted in practical systems.


title: A theorem on meromorphic descent and the specialization of the pro-\'etale fundamental group

abstract: Given a Noetherian formal scheme $\hat X$ over ${\rm Spf}(R)$, where $R$ is a complete DVR, we first prove a theorem of meromorphic descent along a possibly infinite cover of $\hat{X}$. Using this we construct a specialization functor from the category of continuous representations of the pro-\'etale fundamental group of the special fiber to the category of $F$-divided sheaves on the generic fiber. This specialization functor partially recovers the specialization functor of the \'etale fundamental groups. We also express the pro-\'etale fundamental group of a connected scheme $X$ of finite type over a field as coproducts and quotients of the free group and the \'etale fundamental groups of the normalizations of the irreducible components of $X$ and those of its singular loci.


title: Optimal convergence rate of modified Milstein scheme for SDEs with rough fractional diffusions

abstract: We combine the rough path theory and stochastic backward error analysis to develop a new framework for error analysis on numerical schemes. Based on our approach, we prove that the almost sure convergence rate of the modified Milstein scheme for stochastic differential equations driven by multiplicative multidimensional fractional Brownian motion with Hurst parameter $H\in(\frac14,\frac12)$ is $(2H-\frac12)^-$ for sufficiently smooth coefficients, which is optimal in the sense that it is consistent with the result of the corresponding implementable approximation of the L\'evy area of fractional Brownian motion. Our result gives a positive answer to the conjecture proposed in [11] for the case $H\in(\frac13,\frac12)$, and reveals for the first time that numerical schemes constructed by a second-order Taylor expansion do converge for the case $H\in(\frac14,\frac13]$.


title: On entire function $e^{p(z)}\int_0^{z}\beta(t)e^{-p(t)}dt$ with applications to Tumura--Clunie equations and complex dynamics

abstract: Let $p(z)$ be a nonconstant polynomial and $\beta(z)$ be a small entire function of $e^{p(z)}$ in the sense of Nevanlinna. We first describe the growth behavior of the entire function $H(z):=e^{p(z)}\int_0^{z}\beta(t)e^{-p(t)}dt$ on the complex plane $\mathbb{C}$. As an application, we solve entire solutions of Tumura--Clunie type differential equation $f(z)^n+P(z,f)=b_1(z)e^{p_1(z)}+b_2(z)e^{p_2(z)}$, where $b_1(z)$ and $b_2(z)$ are nonzero polynomials, $p_1(z)$ and $p_2(z)$ are two polynomials of the same degree~$k\geq 1$ and $P(z,f)$ is a differential polynomial in $f$ of degree $\leq n-1$ with meromorphic functions of order~$<k$ as coefficients. These results allow us to determine all solutions with relatively few zeros of the second-order differential equation $f''-[b_1(z)e^{p_1(z)}+b_2(z)e^{p_2(z)}+b_3(z)]f=0$, where $b_3(z)$ is a polynomial. We also prove a theorem on certain first-order linear differential equation related to complex dynamics.


title: Surface measures and related functional inequalities on configuration spaces

abstract: Using finite difference operators, we define a notion of boundary and surface measure for configuration sets under Poisson measures. A Margulis-Russo type identity and a co-area formula are stated with applications to deviation inequalities and functional inequalities, and bounds are obtained on the associated isoperimetric constants.


title: Grand challenges and emergent modes of convergence science

abstract: To address complex problems, scholars are increasingly faced with challenges of integrating diverse knowledge domains. We analyzed the evolution of this convergence paradigm in the broad ecosystem of brain science, which provides a real-time testbed for evaluating two modes of cross-domain integration - subject area exploration via expansive learning and cross-disciplinary collaboration among domain experts. We show that research involving both modes features a 16% citation premium relative to a mono-disciplinary baseline. Further comparison of research integrating neighboring versus distant research domains shows that the cross-disciplinary mode is essential for integrating across relatively large disciplinary distances. Yet we find research utilizing cross-domain subject area exploration alone - a convergence shortcut - to be growing in prevalence at roughly 3% per year, significantly faster than the alternative cross-disciplinary mode, despite being less effective at integrating domains and markedly less impactful. By measuring shifts in the prevalence and impact of different convergence modes in the 5-year intervals before and after 2013, our results indicate that these counterproductive patterns may relate to competitive pressures associated with global Human Brain flagship funding initiatives. Without additional policy guidance, such Grand Challenge flagships may unintentionally incentivize such convergence shortcuts, thereby undercutting the advantages of cross-disciplinary teams in tackling challenges calling on convergence.

categories: cs.DL cs.SI econ.GN physics.soc-ph q-fin.EC


title: Secure list decoding and its application to bit-string commitment

abstract: We propose a new concept of secure list decoding, which is related to bit-string commitment. While the conventional list decoding requires that the list contains the transmitted message, secure list decoding requires the following additional security conditions to work as a modification of bit-string commitment. The first additional security condition is the receiver's uncertainty for the transmitted message, which is stronger than the impossibility of the correct decoding, even though the transmitted message is contained in the list. The other additional security condition is the impossibility for the sender to estimate another element of the decoded list except for the transmitted message. The first condition is evaluated by the equivocation rate. The asymptotic property is evaluated by three parameters, the rates of the message and list sizes, and the equivocation rate. We derive the capacity region of this problem. We show that the combination of hash function and secure list decoding yields the conventional bit-string commitment. Our results hold even when the input and output systems are general probability spaces including continuous systems. When the input system is a general probability space, we formulate the abilities of the honest sender and the dishonest sender in a different way.