title: Transition to turbulence in nonuniform coronal loops driven by torsional Alfv\'en waves. II. Extended analysis and effect of magnetic twist

abstract: It has been shown in a previous work that torsional Alfv\'en waves can drive turbulence in nonuniform coronal loops with a purely axial magnetic field. Here we explore the role of the magnetic twist. We model a coronal loop as a transversely nonuniform straight flux tube, anchored in the photosphere, and embedded in a uniform coronal environment. We consider that the magnetic field is twisted and control the strength of magnetic twist by a free parameter of the model. We excite the longitudinally fundamental mode of standing torsional Alfv\'en waves, whose temporal evolution is obtained by means of high-resolution three-dimensional ideal magnetohydrodynamic numerical simulations. We find that phase mixing of torsional Alfv\'en waves creates velocity shear in the direction perpendicular to the magnetic field lines. The velocity shear eventually triggers the Kelvin-Helmholtz instability (KHi). In weakly twisted magnetic tubes, the KHi is able to grow nonlinearly and, subsequently, turbulence is driven in the coronal loop in a similar manner as in the untwisted case. Provided that magnetic twist remains weak, the effect of magnetic twist is to delay the onset of the KHi and to slow down the development of turbulence. In contrast, magnetic tension can suppress the nonlinear growth of the KHi when magnetic twist is strong enough, even if the KHi has locally been excited by the phase-mixing shear. Thus, turbulence is not generated in strongly twisted loops


title: Majorization-minimization for Sparse Nonnegative Matrix Factorization with the $\beta$-divergence

abstract: This article introduces new multiplicative updates for nonnegative matrix factorization with the $\beta$-divergence and sparse regularization of one of the two factors (say, the activation matrix). It is well known that the norm of the other factor (the dictionary matrix) needs to be controlled in order to avoid an ill-posed formulation. Standard practice consists in constraining the columns of the dictionary to have unit norm, which leads to a nontrivial optimization problem. Our approach leverages a reparametrization of the original problem into the optimization of an equivalent scale-invariant objective function. From there, we derive block-descent majorization-minimization algorithms that result in simple multiplicative updates for either $\ell_{1}$-regularization or the more "aggressive" log-regularization. In contrast with other state-of-the-art methods, our algorithms are universal in the sense that they can be applied to any $\beta$-divergence (i.e., any value of $\beta$) and that they come with convergence guarantees. We report numerical comparisons with existing heuristic and Lagrangian methods using various datasets: face images, an audio spectrogram, hyperspectral data, and song play counts. We show that our methods obtain solutions of similar quality at convergence (similar objective values) but with significantly reduced CPU times.


title: Massive galaxy clusters like "El Gordo'' hint at primordial quantum diffusion

abstract: It is generally assumed within the standard cosmological model that initial density perturbations are Gaussian at all scales. However, primordial quantum diffusion unavoidably generates non-Gaussian, exponential tails in the distribution of inflationary perturbations. These exponential tails have direct consequences for the formation of collapsed structures in the universe, as has been studied in the context of primordial black holes. We show that these tails also affect the very-large-scale structures, making heavy clusters like "El Gordo", or large voids like the one associated with the cosmic microwave background cold spot, more probable. We compute the halo mass function and cluster abundance as a function of redshift in the presence of exponential tails. We find that quantum diffusion generically enlarges the number of heavy clusters and depletes subhalos, an effect that cannot be captured by the famed $f_{\mathrm{NL}}$ corrections. These late-universe signatures could thus be fingerprints of quantum dynamics during inflation that should be incorporated in $N$-body simulations and checked against astrophysical data.


title: Dynamic Car Dispatching and Pricing: Revenue and Fairness for Ridesharing Platforms

abstract: A major challenge for ridesharing platforms is to guarantee profit and fairness simultaneously, especially in the presence of misaligned incentives of drivers and riders. We focus on the dispatching-pricing problem to maximize the total revenue while keeping both drivers and riders satisfied. We study the computational complexity of the problem, provide a novel two-phased pricing solution with revenue and fairness guarantees, extend it to stochastic settings and develop a dynamic (a.k.a., learning-while-doing) algorithm that actively collects data to learn the demand distribution during the scheduling process. We also conduct extensive experiments to demonstrate the effectiveness of our algorithms.


title: Stable Centres of Iwahori-Hecke Algebras of type A

abstract: A celebrated result of Farahat and Higman constructs an algebra $\mathrm{FH}$ which "interpolates" the centres $Z(\mathbb{Z}S_n)$ of group algebras of the symmetric groups $S_n$. We extend these results from symmetric group algebras to type $A$ Iwahori-Hecke algebras, $H_n(q)$. In particular, we explain how to construct an algebra $\mathrm{FH}_q$ "interpolating" the centres $Z(H_n(q))$. We prove that $\mathrm{FH}_q$ is isomorphic to $\mathcal{R}[q,q^{-1}] \otimes_{\mathbb{Z}} \Lambda$ (where $\mathcal{R}$ is the ring of integer-valued polynomials, and $\Lambda$ is the ring of symmetric functions). The isomorphism can be described as "evaluation at Jucys-Murphy elements", leading to a proof of a conjecture of Francis and Wang. This yields character formulae for the Geck-Rouquier basis of $Z(H_n(q))$ when acting on Specht modules.


title: Linear regression with unmatched data: a deconvolution perspective

abstract: Consider the regression problem where the response $Y\in\mathbb{R}$ and the covariate $X\in\mathbb{R}^d$ for $d\geq 1$ are \textit{unmatched}. Under this scenario, we do not have access to pairs of observations from the distribution of $(X, Y)$, but instead, we have separate datasets $\{Y_i\}_{i=1}^n$ and $\{X_j\}_{j=1}^m$, possibly collected from different sources. We study this problem assuming that the regression function is linear and the noise distribution is known or can be estimated. We introduce an estimator of the regression vector based on deconvolution and demonstrate its consistency and asymptotic normality under an identifiability assumption. In the general case, we show that our estimator (DLSE: Deconvolution Least Squared Estimator) is consistent in terms of an extended $\ell_2$ norm. Using this observation, we devise a method for semi-supervised learning, i.e., when we have access to a small sample of matched pairs $(X_k, Y_k)$. Several applications with synthetic and real datasets are considered to illustrate the theory.


title: Braids, links, cobordisms and formal groups

abstract: V.V. Sharko in his papers and books has investigated functions on manifolds and cobordism. Braids intimately connect with functions on manifolds. These connections are represented by mapping class groups of corresponding discs, by fundamental groups of corresponding punctured discs, and by some other topological or algebraic structures. This paper presents selected algebraic methods and results of braids, links, cobordism connect with investigations by V.V. Sharko. These includes group theoretic results on braids and links, infinitesimal braid group relations and connections as well as connections on coherent sheaves on smooth schemes, a sketch of our algorithm for constructing of Lazard`s one dimensional universal commutative formal group and selected results on applications of commutative formal groups to cobordism theory.


title: Bayesian hierarchical modelling of the $\mathrm{M_{\star}}$-SFR relation from 1<z<6 in ASTRODEEP

abstract: The Hubble Frontier Fields represent the opportunity to probe the high-redshift evolution of the main sequence of star-forming galaxies to lower masses than possible in blank fields thanks to foreground lensing of massive galaxy clusters. We use the BEAGLE SED-fitting code to derive stellar masses, $\mathrm{M_{\star}}=\log(M/\mathrm{M_{\odot}})$, SFRs, $\Psi=\log(\psi/\mathrm{M_{\odot}}\,\mathrm{yr}^{-1})$ and redshifts from galaxies within the ASTRODEEP catalogue. We fit a fully Bayesian hierarchical model of the main sequence over $1.25<z<6$ of the form $\Psi = \alpha_\mathrm{9.7}(z) + \beta(\mathrm{M_{\star}}-9.7) + \mathcal{N}(0,\sigma^2)$ while explicitly modelling the outlier distribution. The redshift-dependent intercept at $\mathrm{M_{\star}}=9.7$ is parametrized as $\alpha_\mathrm{9.7}(z) = \log[N (1+z)^{\gamma}] + 0.7$. Our results agree with an increase in normalization of the main sequence to high redshifts that follows the redshift-dependent rate of accretion of gas onto dark matter halos with $\gamma=2.40^{+0.18}_{-0.18}$. We measure a slope and intrinsic scatter of $\beta=0.79^{+0.03}_{-0.04}$ and $\sigma=0.26^{+0.02}_{-0.02}$. We find that the sampling of the SED provided by the combination of filters (Hubble + ground-based Ks-band + Spitzer 3.6 and 4.5 $\mathrm{\mu m}$) is insufficient to constrain $\mathrm{M_{\star}}$ and $\Psi$ over the full dynamic range of the observed main sequence, even at the lowest redshifts studied. While this filter set represents the best current sampling of high-redshift galaxy SEDs out to $z>3$, measurements of the main sequence to low masses and high redshifts still strongly depend on priors employed in SED fitting (as well as other fitting assumptions). Future data-sets with JWST should improve this.


title: Microscopic dynamics and Bose-Einstein condensation in liquid helium

abstract: We review fundamental problems involved in liquid theory including both classical and quantum liquids. Understanding classical liquids involves exploring details of their microscopic dynamics and its consequences. Here, we apply the same general idea to quantum liquids. We discuss momentum condensation in liquid helium which is consistent with microscopic dynamics in liquids and high mobility of liquid atoms. We propose that mobile transit atoms accumulate in the finite-energy state where the transit speed is close to the speed of sound. In this state, the transit energy is close to the zero-point energy. In momentum space, the accumulation operates on a sphere with the radius set by interatomic spacing and corresponds to zero net momentum. We show that this picture is supported by experiments, including the measured kinetic energy of helium atoms below the superfluid transition and sharp peaks of scattered intensity at predicted energy. We discuss the implications of this picture including the macroscopic wave function and superfluidity.


title: PointNorm: Dual Normalization is All You Need for Point Cloud Analysis

abstract: Point cloud analysis is challenging due to the irregularity of the point cloud data structure. Existing works typically employ the ad-hoc sampling-grouping operation of PointNet++, followed by sophisticated local and/or global feature extractors for leveraging the 3D geometry of the point cloud. Unfortunately, the sampling-grouping operations do not address the point cloud's irregularity, whereas the intricate local and/or global feature extractors led to poor computational efficiency. In this paper, we introduce a novel DualNorm module after the sampling-grouping operation to effectively and efficiently address the irregularity issue. The DualNorm module consists of Point Normalization, which normalizes the grouped points to the sampled points, and Reverse Point Normalization, which normalizes the sampled points to the grouped points. The proposed framework, PointNorm, utilizes local mean and global standard deviation to benefit from both local and global features while maintaining a faithful inference speed. Experiments show that we achieved excellent accuracy and efficiency on ModelNet40 classification, ScanObjectNN classification, ShapeNetPart Part Segmentation, and S3DIS Semantic Segmentation. Code is available at https://github.com/ShenZheng2000/PointNorm-for-Point-Cloud-Analysis.


title: Non-Myopic Multifidelity Bayesian Optimization

abstract: Bayesian optimization is a popular framework for the optimization of black box functions. Multifidelity methods allows to accelerate Bayesian optimization by exploiting low-fidelity representations of expensive objective functions. Popular multifidelity Bayesian strategies rely on sampling policies that account for the immediate reward obtained evaluating the objective function at a specific input, precluding greater informative gains that might be obtained looking ahead more steps. This paper proposes a non-myopic multifidelity Bayesian framework to grasp the long-term reward from future steps of the optimization. Our computational strategy comes with a two-step lookahead multifidelity acquisition function that maximizes the cumulative reward obtained measuring the improvement in the solution over two steps ahead. We demonstrate that the proposed algorithm outperforms a standard multifidelity Bayesian framework on popular benchmark optimization problems.


title: Search for Lorentz-Invariance Violation with the first KATRIN data

abstract: Some extensions of the Standard Model of Particle Physics allow for Lorentz invariance and Charge-Parity-Time (CPT)-invariance violations. In the neutrino sector strong constraints have been set by neutrino-oscillation and time-of-flight experiments. However, some Lorentz-invariance-violating parameters are not accessible via these probes. In this work, we focus on the parameters $(a_{\text{of}}^{(3)})_{00}$, $(a_{\text{of}}^{(3)})_{10}$ and $(a_{\text{of}}^{(3)})_{11}$ which would manifest themselves in a non-isotropic beta-decaying source as a sidereal oscillation and an overall shift of the spectral endpoint. Based on the data of the first scientific run of the KATRIN experiment, we set the first limit on $\left|(a_{\text{of}}^{(3)})_{11}\right|$ of $< 3.7\cdot10^{-6}$ GeV at 90\% confidence level. Moreover, we derive new constraints on $(a_{\text{of}}^{(3)})_{00}$ and $(a_{\text{of}}^{(3)})_{10}$.


title: Event-triggered Control of Port-Hamiltonian Systems under Time-delay Communication

abstract: We study the problem of periodic event-triggered control of interconnected port-Hamiltonian systems subject to time-varying delays in their communication. In particular, we design a threshold parameter for the event-triggering condition, a sampling period, and a maximum allowable delay such that interconnected port-Hamiltonian control systems with periodic event-triggering mechanism under a time-delayed communication are able to achieve asymptotically stable behaviour. Simulation results are presented to validate the theory.


title: Modal properties of dielectric bowtie cavities with deep sub-wavelength confinement

abstract: We present a design for an optical dielectric bowtie cavity which features deep sub-wavelength confinement of light. The cavity is derived via simplification of a complex geometry identified through inverse design by topology optimization, and it successfully retains the extreme properties of the original structure, including an effective mode volume of $V_\text{eff}= 0.083 \pm 0.001 \; (\lambda_\text{c}/2n_\text{Si})^3$ at its center. Based on this design, we present a modal analysis to show that the Purcell factor can be well described by a single quasinormal mode in a wide bandwidth of interest. Owing to the small mode volume, moreover, the cavity exhibits a remarkable sensitivity to local shape deformations, which we show to be well described by perturbation theory. The intuitive simplification approach to inverse design geometries coupled with the quasinormal mode analysis demonstrated in this work provides a powerful modeling framework for the emerging field of dielectric cavities with deep sub-wavelength confinement.


title: GAN with an Auxiliary Regressor for the Fast Simulation of the Electromagnetic Calorimeter Response

abstract: High energy physics experiments essentially rely on simulated data for physics analyses. However, running detailed simulation models requires a tremendous amount of computation resources. New approaches to speed up detector simulation are therefore needed. The generation of calorimeter responses is often the most expensive component of the simulation chain for HEP experiments. It was shown that deep learning techniques, especially Generative Adversarial Networks, may be used to reproduce the calorimeter response. However, those applications are challenging, as the generated responses need evaluation not only in terms of image consistency: different physics-based quality metrics should be also taken into consideration. In our work, we develop a multitask GAN-based framework with the goal to speed up the response generation of the Electromagnetic Calorimeter (ECAL) of the LHCb detector at LHC. We introduce the Auxiliary Regressor as a second task to evaluate a proxy metric of the given input that is used by the Discriminator of the GAN. We show that this approach improves the stability of GAN and the model produces samples with better physics distributions.


title: Left Ventricle Contouring of Apical Three-Chamber Views on 2D Echocardiography

abstract: We propose a new method to automatically contour the left ventricle on 2D echocardiographic images. Unlike most existing segmentation methods, which are based on predicting segmentation masks, we focus at predicting the endocardial contour and the key landmark points within this contour (basal points and apex). This provides a representation that is closer to how experts perform manual annotations and hence produce results that are physiologically more plausible. Our proposed method uses a two-headed network based on the U-Net architecture. One head predicts the 7 contour points, and the other head predicts a distance map to the contour. This approach was compared to the U-Net and to a point based approach, achieving performance gains of up to 30\% in terms of landmark localisation (<4.5mm) and distance to the ground truth contour (<3.5mm).


title: Two-dimensional delta-Bose gas: skew-product relative motion and exact non-Gaussianity for the stochastic heat equation

abstract: For the two-dimensional delta-Bose gas, we prove the first Feynman-Kac type diffusion representation of the Schr\"odinger semigroup generated by the relative motion of two particles. The representation builds on a unique skew-product diffusion, in the sense of Erickson [22], such that the Markov process ${\rm BES}(0,\beta{\downarrow})$ from Donati-Martin and Yor [20] defines the radial part. Moreover, the law of this skew-product diffusion is singular to the law of planar Brownian motion. As an application, we determine some exact forms of non-Gaussianity for the stochastic heat equation. We give two different proofs of the representation: one applies the original characterization of ${\rm BES}(0,\beta{\downarrow})$ by excursion theory, and the other extends [20] by approximations of SDEs and resolvents. In the latter, we establish a new characterization of ${\rm BES}(0,\beta{\downarrow})$ as a strongly well-posed SDE. The approximations of resolvents construct the two-dimensional delta potentials from local times of lower-dimensional Bessel processes, which may be of independent interest for being an alternative to the long-standing approach of using mollifiers.


title: Symmetry-Aware Transformer-based Mirror Detection

abstract: Mirror detection aims to identify the mirror regions in the given input image. Existing works mainly focus on integrating the semantic features and structural features to mine specific relations between mirror and non-mirror regions, or introducing mirror properties like depth or chirality to help analyze the existence of mirrors. In this work, we observe that a real object typically forms a loose symmetry relationship with its corresponding reflection in the mirror, which is beneficial in distinguishing mirrors from real objects. Based on this observation, we propose a dual-path Symmetry-Aware Transformer-based mirror detection Network (SATNet), which includes two novel modules: Symmetry-Aware Attention Module (SAAM) and Contrast and Fusion Decoder Module (CFDM). Specifically, we first adopt a transformer backbone to model global information aggregation in images, extracting multi-scale features in two paths. We then feed the high-level dual-path features to SAAMs to capture the symmetry relations. Finally, we fuse the dual-path features and refine our prediction maps progressively with CFDMs to obtain the final mirror mask. Experimental results show that SATNet outperforms both RGB and RGB-D mirror detection methods on all available mirror detection datasets. Codes and trained models are available at: https://github.com/tyhuang0428/SATNet.


title: 6D Camera Relocalization in Visually Ambiguous Extreme Environments

abstract: We propose a novel method to reliably estimate the pose of a camera given a sequence of images acquired in extreme environments such as deep seas or extraterrestrial terrains. Data acquired under these challenging conditions are corrupted by textureless surfaces, image degradation, and presence of repetitive and highly ambiguous structures. When naively deployed, the state-of-the-art methods can fail in those scenarios as confirmed by our empirical analysis. In this paper, we attempt to make camera relocalization work in these extreme situations. To this end, we propose: (i) a hierarchical localization system, where we leverage temporal information and (ii) a novel environment-aware image enhancement method to boost the robustness and accuracy. Our extensive experimental results demonstrate superior performance in favor of our method under two extreme settings: localizing an autonomous underwater vehicle and localizing a planetary rover in a Mars-like desert. In addition, our method achieves comparable performance with state-of-the-art methods on the indoor benchmark (7-Scenes dataset) using only 20% training data.


title: A nonstandard proof of continuity of affine varieties

abstract: Extending the classical result that the roots of a polynomial with coefficients in $\mathbf{C}$ are continuous functions of the coefficients of the polynomial, nonstandard analysis is used to prove that if $\mathcal{F} = \{f_{\lambda} :\lambda \in \Lambda\}$ is a set of polynomials in $\mathbf{C}[t_1,\ldots, t_n]$ and if $^*\mathcal{G} = \{g_{\lambda} :\lambda \in \Lambda\}$ is a set of polynomials in $^*\mathbf{C}_0[t_1,\ldots, t_n]$ such that $g_{\lambda}$ is an infinitesimal deformation of $f_{\lambda}$ for all $\lambda \in \Lambda$, then the nonstandard affine variety $^*V_0(\mathcal{G})$ is an infinitesimal deformation of the affine variety $V(\mathcal{F})$.


title: Persistence and the Sheaf-Function Correspondence

abstract: The sheaf-function correspondence identifies the group of constructible functions on a real analytic manifold $M$ with the Grothendieck group of constructible sheaves on $M$. When $M$ is a finite dimensional real vector space, Kashiwara-Schapira have recently introduced the convolution distance between sheaves of $k$-vector spaces on $M$. In this paper, we characterize distances on the group of constructible functions on a real finite dimensional vector space that can be controlled by the convolution distance through the sheaf-function correspondence. Our main result asserts that such distances are almost trivial: they vanish as soon as two constructible functions have the same Euler integral. We formulate consequences of our result for Topological Data Analysis: there cannot exists non-trivial additive invariants of persistence modules that are continuous for the interleaving distance.


title: QT-Routenet: Improved GNN generalization to larger 5G networks by fine-tuning predictions from queueing theory

abstract: In order to promote the use of machine learning in 5G, the International Telecommunication Union (ITU) proposed in 2021 the second edition of the ITU AI/ML in 5G challenge, with over 1600 participants from 82 countries. This work details the second place solution overall, which is also the winning solution of the Graph Neural Networking Challenge 2021. We tackle the problem of generalization when applying a model to a 5G network that may have longer paths and larger link capacities than the ones observed in training. To achieve this, we propose to first extract robust features related to Queueing Theory (QT), and then fine-tune the analytical baseline prediction using a modification of the Routenet Graph Neural Network (GNN) model. The proposed solution generalizes much better than simply using Routenet, and manages to reduce the analytical baseline's 10.42 mean absolute percent error to 1.45 (1.27 with an ensemble). This suggests that making small changes to an approximate model that is known to be robust can be an effective way to improve accuracy without compromising generalization.


title: Search for keV-scale Sterile Neutrinos with first KATRIN Data

abstract: In this work we present a keV-scale sterile-neutrino search with the first tritium data of the KATRIN experiment, acquired in the commissioning run in 2018. KATRIN performs a spectroscopic measurement of the tritium $\beta$-decay spectrum with the main goal of directly determining the effective electron anti-neutrino mass. During this commissioning phase a lower tritium activity facilitated the search for sterile neutrinos with a mass of up to $1.6\, \mathrm{keV}$. We do not find a signal and set an exclusion limit on the sterile-to-active mixing amplitude of down to $\sin^2\theta < 5\cdot10^{-4}$ ($95\,\%$ C.L.), improving current laboratory-based bounds in the sterile-neutrino mass range between 0.1 and $1.0\, \mathrm{keV}$.


title: Memristive and tunneling effects in 3D interconnected silver nanowires

abstract: Due to their memristive properties nanowire networks are very promising for neuromorphic computing applications. Indeed, the resistance of such systems can evolve with the input voltage or current as it confers a synaptic behaviour to the device. Here, we propose a network of silver nanowires (Ag-NWs) which are grown in a nanopourous membrane with interconnected nanopores by electrodeposition. This bottom-up approach fabrication method gives a conducting network with a 3D architecture and a high density of Ag-NWs. The resulting 3D interconnected Ag-NW network exhibits a high initial resistance as well as a memristive behavior. It is expected to arise from the creation and the destruction of conducting silver filaments inside the Ag-NW network. Moreover, after several cycles of measurement, the resistance of the network switches from a high resistance regime, in the GOhm range, with a tunnel conduction to a low resistance regime, in the kOhm range.


title: Learning robust marking policies for adaptive mesh refinement

abstract: In this work, we revisit the marking decisions made in the standard adaptive finite element method (AFEM). Experience shows that a na\"{i}ve marking policy leads to inefficient use of computational resources for adaptive mesh refinement (AMR). Consequently, using AFEM in practice often involves ad-hoc or time-consuming offline parameter tuning to set appropriate parameters for the marking subroutine. To address these practical concerns, we recast AMR as a Markov decision process in which refinement parameters can be selected on-the-fly at run time, without the need for pre-tuning by expert users. In this new paradigm, the refinement parameters are also chosen adaptively via a marking policy that can be optimized using methods from reinforcement learning. We use the Poisson equation to demonstrate our techniques on $h$- and $hp$-refinement benchmark problems, and our experiments suggest that superior marking policies remain undiscovered for many classical AFEM applications. Furthermore, an unexpected observation from this work is that marking policies trained on one family of PDEs are sometimes robust enough to perform well on problems far outside the training family. For illustration, we show that a simple $hp$-refinement policy trained on 2D domains with only a single re-entrant corner can be deployed on far more complicated 2D domains, and even 3D domains, without significant performance loss. For reproduction and broader adoption, we accompany this work with an open-source implementation of our methods.


title: Quasi Normal Modes of Ayon-Beato Garcia Regular Black Holes For Scalar Field

abstract: In this paper, we compute quasi-normal modes of ABG black holes (which has a non-linear electrodynamical source) using the WKB methods and AIM. A comparison between the spectrum of QNMs calculated by both methods is made. We analyse how the spectrum of QNMs depends on the black hole parameters, multipole number and overtone number and establish that the ABG black hole is stable against the scalar field.


title: Adaptive Edge Content Delivery Networks for Web-Scale File Systems

abstract: The InterPlanetary File System (IPFS) is an hyper-media distribution protocol, addressed by content and identities. It aims to make the web faster, safer, and more open. The JavaScript implementation of IPFS runs on the browser, benefiting from the mass adoption potential that it yields. Startrail takes advantage of the IPFS ecosystem and strives to further evolve it, making it more scalable and performant through the implementation of an adaptive network caching mechanism. Our solution aims to add resilience to IPFS and improve its overall scalability, by avoiding overloading the nodes providing highly popular content, particularly during flash-crowd-like conditions where popularity and demand grow suddenly. We add a novel crucial key component to enable an IPFS-based decentralized Content Distribution Network (CDN). Following a peer-to-peer architecture, it runs on a scalable, highly available network of untrusted nodes that distribute immutable and authenticated objects which are cached progressively towards the sources of requests.


title: Efficient error and variance estimation for randomized matrix computations

abstract: Randomized matrix algorithms have become workhorse tools in scientific computing and machine learning. To use these algorithms safely in applications, they should be coupled with posterior error estimates to assess the quality of the output. To meet this need, this paper proposes two diagnostics: a leave-one-out error estimator for randomized low-rank approximations and a jackknife resampling method to estimate the variance of the output of a randomized matrix computation. Both of these diagnostics are rapid to compute for randomized low-rank approximation algorithms such as the randomized SVD and randomized Nystr\"om approximation, and they provide useful information that can be used to assess the quality of the computed output and guide algorithmic parameter choices.


title: TCT: Convexifying Federated Learning using Bootstrapped Neural Tangent Kernels

abstract: State-of-the-art federated learning methods can perform far worse than their centralized counterparts when clients have dissimilar data distributions. For neural networks, even when centralized SGD easily finds a solution that is simultaneously performant for all clients, current federated optimization methods fail to converge to a comparable solution. We show that this performance disparity can largely be attributed to optimization challenges presented by nonconvexity. Specifically, we find that the early layers of the network do learn useful features, but the final layers fail to make use of them. That is, federated optimization applied to this non-convex problem distorts the learning of the final layers. Leveraging this observation, we propose a Train-Convexify-Train (TCT) procedure to sidestep this issue: first, learn features using off-the-shelf methods (e.g., FedAvg); then, optimize a convexified problem obtained from the network's empirical neural tangent kernel approximation. Our technique yields accuracy improvements of up to +36% on FMNIST and +37% on CIFAR10 when clients have dissimilar data.


title: Pinching Azumaya algebras

abstract: We show, that for a morphism of schemes from X to Y, that is a finite modification in finitely many closed points, a cohomological Brauer class on Y is represented by an Azumaya algebra if its pullback to X is represented by an Azumaya algebra. Part of the proof uses an extension of a result by Ferrand, on pinching of finite locally free sheaves, to Azumaya algebras.


title: You Only Align Once: Bidirectional Interaction for Spatial-Temporal Video Super-Resolution

abstract: Spatial-Temporal Video Super-Resolution (ST-VSR) technology generates high-quality videos with higher resolution and higher frame rates. Existing advanced methods accomplish ST-VSR tasks through the association of Spatial and Temporal video super-resolution (S-VSR and T-VSR). These methods require two alignments and fusions in S-VSR and T-VSR, which is obviously redundant and fails to sufficiently explore the information flow of consecutive spatial LR frames. Although bidirectional learning (future-to-past and past-to-future) was introduced to cover all input frames, the direct fusion of final predictions fails to sufficiently exploit intrinsic correlations of bidirectional motion learning and spatial information from all frames. We propose an effective yet efficient recurrent network with bidirectional interaction for ST-VSR, where only one alignment and fusion is needed. Specifically, it first performs backward inference from future to past, and then follows forward inference to super-resolve intermediate frames. The backward and forward inferences are assigned to learn structures and details to simplify the learning task with joint optimizations. Furthermore, a Hybrid Fusion Module (HFM) is designed to aggregate and distill information to refine spatial information and reconstruct high-quality video frames. Extensive experiments on two public datasets demonstrate that our method outperforms state-of-the-art methods in efficiency, and reduces calculation cost by about 22%.


title: A guideline proposal for minimizing cybersickness in VR-based serious games and applications

abstract: Head-mounted displays (HMDs) are popular immersive tools in general, not limited to entertainment but also for education, military, and serious games for health. While these displays have strong popularity, they still have user experience issues, triggering possible symptoms of discomfort to users. This condition is known as cybersickness (CS) and is one of the most popular research topics tied to virtual reality (VR) issues. We first present the main strategies focused on minimizing cybersickness problems in virtual reality. Following this, we propose a guideline framework based on CS causes such as locomotion, acceleration, the field of view, depth of field, degree of freedom, exposition use time, latency-lag, static rest frame, and camera rotation. Additionally, serious games applications and broader categories of games can also adopt it. Additionally, we categorized the imminent challenges for CS minimization into four different items. Conclusively, this work contributes as a consulting reference to enable VR developers and designers to optimize their VR users' experience and VR serious games.


title: Giant orbital Hall effect and orbital-to-spin conversion in 3d, 5d, and 4f metallic heterostructures

abstract: The orbital Hall effect provides an alternative means to the spin Hall effect to convert a charge current into a flow of angular momentum. Recently, compelling signatures of orbital Hall effects have been identified in 3d transition metals. Here, we report a systematic study of the generation, transmission, and conversion of orbital currents in heterostructures comprising 3d, 5d, and 4f metals. We show that the orbital Hall conductivity of Cr reaches giant values of the order of 5*10^5 Ohm^{-1} m^{-1} and that Pt presents a strong orbital Hall effect in addition to the spin Hall effect. Measurements performed as a function of thickness of nonmagnetic Cr, Mn, and Pt layers and ferromagnetic Co and Ni layers reveal how the orbital and spin currents compete or assist each other in determining the spin-orbit torques acting on the magnetic layer. We further show how this interplay can be drastically modulated by introducing 4 f spacers between the nonmagnetic and magnetic layers. Gd and Tb act as very efficient orbital-to-spin current converters, boosting the spin-orbit torques generated by Cr by a factor of 4 and reversing the sign of the torques generated by Pt. To interpret our results, we present a generalized drift-diffusion model that includes both spin and orbital Hall effects and describes their interconversion mediated by spin-orbit coupling.


title: Toda systems for Takiff algebras

abstract: We study completely integrable systems attached to Takiff algebras $\mathfrak{g}_N$, extending open Toda systems of split simple Lie algebras $\mathfrak{g}$. With respect to Darboux coordinates on coadjoint orbits $\mathcal{O}$, the potentials of the hamiltonians are products of polynomial and exponential functions. General solutions for equations of motion for $\mathfrak{g}_N$ are obtained using differential operators called jet transformations. These results are applied to a $3$-body problem based on $\mathfrak{sl}(2)$, and to an extension of soliton solutions for $A_\infty$ to associated Takiff algebras. The new classical integrable systems are then lifted to families of commuting operators in an enveloping algebra, solving a Vinberg problem and quantizing the Poisson algebra of functions on $\mathcal{O}$.


title: Polyphonic sound event detection for highly dense birdsong scenes

abstract: One hour before sunrise, one can experience the dawn chorus where birds from different species sing together. In this scenario, high levels of polyphony, as in the number of overlapping sound sources, are prone to happen resulting in a complex acoustic outcome. Sound Event Detection (SED) tasks analyze acoustic scenarios in order to identify the occurring events and their respective temporal information. However, highly dense scenarios can be hard to process and have not been studied in depth. Here we show, using a Convolutional Recurrent Neural Network (CRNN), how birdsong polyphonic scenarios can be detected when dealing with higher polyphony and how effectively this type of model can face a very dense scene with up to 10 overlapping birds. We found that models trained with denser examples (i.e., higher polyphony) learn at a similar rate as models that used simpler samples in their training set. Additionally, the model trained with the densest samples maintained a consistent score for all polyphonies, while the model trained with the least dense samples degraded as the polyphony increased. Our results demonstrate that highly dense acoustic scenarios can be dealt with using CRNNs. We expect that this study serves as a starting point for working on highly populated bird scenarios such as dawn chorus or other dense acoustic problems.


title: A sharpened energy-Strichartz inequality for the wave equation

abstract: We consider the sharp Strichartz estimate for the wave equation on $\mathbb R^{1+5}$ in the energy space, due to Bez and Rogers. We show that it can be refined by adding a term proportional to the distance from the set of maximisers, in the spirit of the classical sharpened Sobolev estimate of Bianchi and Egnell.


title: Joint Prediction of Monocular Depth and Structure using Planar and Parallax Geometry

abstract: Supervised learning depth estimation methods can achieve good performance when trained on high-quality ground-truth, like LiDAR data. However, LiDAR can only generate sparse 3D maps which causes losing information. Obtaining high-quality ground-truth depth data per pixel is difficult to acquire. In order to overcome this limitation, we propose a novel approach combining structure information from a promising Plane and Parallax geometry pipeline with depth information into a U-Net supervised learning network, which results in quantitative and qualitative improvement compared to existing popular learning-based methods. In particular, the model is evaluated on two large-scale and challenging datasets: KITTI Vision Benchmark and Cityscapes dataset and achieve the best performance in terms of relative error. Compared with pure depth supervision models, our model has impressive performance on depth prediction of thin objects and edges, and compared to structure prediction baseline, our model performs more robustly.


title: Influence of Solution Parameters on Phase Formation and Morphology of Electrospun Poly(vinylidene fluoride) Nanofiber

abstract: Poly(vinylidene fluoride), PVDF nanofibers were prepared using electrospinning method and the influence of electrospinning parameters such as PVDF concentration, DMF:Acetone ratio on the formation of different phases ({\alpha}, \b{eta}, and {\gamma}) and desired morphology were investigated. A detail analysis of X-ray diffraction (XRD) and Fourier-transform infrared (FTIR) spectrometer data by deconvoluting the peaks were carried out for assigning and accurately identifying the peaks associated to a particular PVDF phase. SEM micrographs further helped in understanding the mechanism of desired \b{eta}-phase as well as in obtaining a condition for bead-free nanofibers. Interestingly, the mechanism of formation of phases are found mostly governed by the balance between surface tension and viscosity which is controlled by PVDF concentration and acetone content in the solution. A PVDF concentration of 20% (w/v) and DMF:Acetone ratio 1:1 was found suitable for maximum \b{eta}-phase in a completely bead-free nanofibers. The methodology and understanding of underlying mechanism of desired phase and morphology suggest a guideline for optimization of \b{eta}-phase and bead-free fibers in PVDF-based nanocomposite as well.


title: Massey products in the \'etale cohomology of number fields

abstract: We give formulas for 3-fold Massey products in the \'etale cohomology of the ring of integers of a number field and use these to find the first known examples of imaginary quadratic fields with class group of $p$-rank two possessing an infinite $p$-class field tower, where $p$ is an odd prime. Furthermore, a necessary and sufficient condition, in terms of class groups of $p$-extensions, for the vanishing of 3-fold Massey products is given. As a consequence, we give an elementary and sufficient condition for the infinitude of class field towers of imaginary quadratic fields. We also disprove McLeman's $(3,3)$-conjecture. Lastly, we relate the vanishing of Massey products to the existence of Galois representations of $G_{\mathbb{Q},S}$ which realize an unexpectedly large class group for certain extensions of a quadratic imaginary number field.


title: The Geometrically Frustrated Spin Glass (Fe1-pGap)2TiO5

abstract: The unusual anisotropy of the spin glass transition in the pseudobrookite system Fe$_2$TiO$_5$ has been interpreted as arising from an induced, van der Waals-like, interaction among magnetic clusters. Here we present susceptibility ($\chi$) and specific heat data (C) for Fe2TiO5 diluted with non-magnetic Ga, (Fe$_{1-p}$Ga$_p$)$_2$TiO$_5$, for disorder parameter p = 0, 0.11, and 0.42, and elastic neutron scattering data for p = 0.20. A uniform suppression of T{_g} is observed upon increasing p, along with a value of $\chi (T_g)$ that increases as T$_g$ decreases, i.e. $d\chi(T_g)/dT_g< 0$. We also observe C(T) $\propto$ T$^2$ in the low temperature limit. The observed behavior places (Fe$_{1-p}$Ga$_p$)$_2$TiO$_5$ in the category of a strongly geometrically frustrated spin glass.


title: Contextual Decision Trees

abstract: Focusing on Random Forests, we propose a multi-armed contextual bandit recommendation framework for feature-based selection of a single shallow tree of the learned ensemble. The trained system, which works on top of the Random Forest, dynamically identifies a base predictor that is responsible for providing the final output. In this way, we obtain local interpretations by observing the rules of the recommended tree. The carried out experiments reveal that our dynamic method is superior to an independent fitted CART decision tree and comparable to the whole black-box Random Forest in terms of predictive performances.


title: Deep Transformer Model with Pre-Layer Normalization for COVID-19 Growth Prediction

abstract: Coronavirus disease or COVID-19 is an infectious disease caused by the SARS-CoV-2 virus. The first confirmed case caused by this virus was found at the end of December 2019 in Wuhan City, China. This case then spread throughout the world, including Indonesia. Therefore, the COVID-19 case was designated as a global pandemic by WHO. The growth of COVID-19 cases, especially in Indonesia, can be predicted using several approaches, such as the Deep Neural Network (DNN). One of the DNN models that can be used is Deep Transformer which can predict time series. The model is trained with several test scenarios to get the best model. The evaluation is finding the best hyperparameters. Then, further evaluation was carried out using the best hyperparameters setting of the number of prediction days, the optimizer, the number of features, and comparison with the former models of the Long Short-Term Memory (LSTM) and Recurrent Neural Network (RNN). All evaluations used metric of the Mean Absolute Percentage Error (MAPE). Based on the results of the evaluations, Deep Transformer produces the best results when using the Pre-Layer Normalization and predicting one day ahead with a MAPE value of 18.83. Furthermore, the model trained with the Adamax optimizer obtains the best performance among other tested optimizers. The performance of the Deep Transformer also exceeds other test models, which are LSTM and RNN.


title: Shrinkage Estimation of Higher Order Bochner Integrals

abstract: We consider shrinkage estimation of higher order Hilbert space valued Bochner integrals in a non-parametric setting. We propose estimators that shrink the $U$-statistic estimator of the Bochner integral towards a pre-specified target element in the Hilbert space. Depending on the degeneracy of the kernel of the $U$-statistic, we construct consistent shrinkage estimators with fast rates of convergence, and develop oracle inequalities comparing the risks of the the $U$-statistic estimator and its shrinkage version. Surprisingly, we show that the shrinkage estimator designed by assuming complete degeneracy of the kernel of the $U$-statistic is a consistent estimator even when the kernel is not complete degenerate. This work subsumes and improves upon Krikamol et al., 2016, JMLR and Zhou et al., 2019, JMVA, which only handle mean element and covariance operator estimation in a reproducing kernel Hilbert space. We also specialize our results to normal mean estimation and show that for $d\ge 3$, the proposed estimator strictly improves upon the sample mean in terms of the mean squared error.


title: Smooth Anonymity for Sparse Graphs

abstract: When working with user data providing well-defined privacy guarantees is paramount. In this work, we aim to manipulate and share an entire sparse dataset with a third party privately. In fact, differential privacy has emerged as the gold standard of privacy, however, when it comes to sharing sparse datasets, e.g. sparse networks, as one of our main results, we prove that \emph{any} differentially private mechanism that maintains a reasonable similarity with the initial dataset is doomed to have a very weak privacy guarantee. In such situations, we need to look into other privacy notions such as $k$-anonymity. In this work, we consider a variation of $k$-anonymity, which we call smooth-$k$-anonymity, and design simple large-scale algorithms that efficiently provide smooth-$k$-anonymity. We further perform an empirical evaluation to back our theoretical guarantees and show that our algorithm improves the performance in downstream machine learning tasks on anonymized data.


title: Novel Superposed Kinklike and Pulselike Solutions for Several Nonlocal Nonlinear Equations

abstract: We show that a number of nonlocal nonlinear equations including the Ablowitz-Musslimani and the Yang variant of the nonlocal nonlinear Schr\"od-inger (NLS) equation, nonlocal modified Korteweg de Vries (mKdV) equation as well as the nonlocal Hirota equation admit novel kinklike and pulselike superposed periodic solutions. Further, we show that the nonlocal mKdV equation, in addition also admits the superposed (hyperbolic) kink-antikink solution. Besides, we show that while the nonlocal Ablowitz-Musslimani variant of the NLS admits complex (parity-time reversal or) PT-invariant kink and pulse solutions, neither the local NLS nor the Yang variant of the nonlocal NLS admits such solutions. Finally, except for the Yang variant of the nonlocal NLS, we show that the other three nonlocal equations admit both the kink and pulse solutions in the same model.


title: Developing an NLP-based Recommender System for the Ethical, Legal, and Social Implications of Synthetic Biology

abstract: Synthetic biology is an emerging field that involves the engineering and re-design of organisms for purposes such as food security, health, and environmental protection. As such, it poses numerous ethical, legal, and social implications (ELSI) for researchers and policy makers. Various efforts to ensure socially responsible synthetic biology are underway. Policy making is one regulatory avenue, and other initiatives have sought to embed social scientists and ethicists on synthetic biology projects. However, given the nascency of synthetic biology, the number of heterogeneous domains it spans, and the open nature of many ethical questions, it has proven challenging to establish widespread concrete policies, and including social scientists and ethicists on synthetic biology teams has met with mixed success. This text proposes a different approach, asking instead is it possible to develop a well-performing recommender model based upon natural language processing (NLP) to connect synthetic biologists with information on the ELSI of their specific research? This recommender was developed as part of a larger project building a Synthetic Biology Knowledge System (SBKS) to accelerate discovery and exploration of the synthetic biology design space. Our approach aims to distill for synthetic biologists relevant ethical and social scientific information and embed it into synthetic biology research workflows.


title: Contact behaviour of simulated rough spheres generated with spherical harmonics

abstract: Normal contact behaviour between non adhesive fractal rough particles is studied using a finite element method (FEM). A series of spherical grain surfaces with distinguished roughness features are generated by means of Spherical Harmonics. These surfaces are described by two roughness descriptors, namely, relative roughness (Rr) and fractal dimension (FD). The contact behaviour of rough spheres with a rigid flat surface is simulated using FEM to quantify the influences of surface structure and sphere morphology by focusing on contact stiffness and true contact area. The dependence of normal contact stiffness (k) on applied normal force (F) is found to follow a power law over four orders of magnitude, with both alpha and beta being highly correlated with Rr and FD. With increasing load, the power exponent converges to that of Hertzian contact, e.g., 1/3, independent of Rr. Regions of true contact evolve through the formation of new microcontacts and their progressive merging, meanwhile the area distributions of contact island induced by various forces tend to obey similar Weibull distributions due to fractal nature in their surfaces. Contacts with larger values of Rr are found to produce contact contours with higher fractal dimension as calculated by a 2D box-counting method. Our results suggest that the correlation between radial lengths in a quasi spherical particle should be considered in studying contact behaviour.


title: Iterative Linear Quadratic Optimization for Nonlinear Control: Differentiable Programming Algorithmic Templates

abstract: We present the implementation of nonlinear control algorithms based on linear and quadratic approximations of the objective from a functional viewpoint. We present a gradient descent, a Gauss-Newton method, a Newton method, differential dynamic programming approaches with linear quadratic or quadratic approximations, various line-search strategies, and regularized variants of these algorithms. We derive the computational complexities of all algorithms in a differentiable programming framework and present sufficient optimality conditions. We compare the algorithms on several benchmarks, such as autonomous car racing using a bicycle model of a car. The algorithms are coded in a differentiable programming language in a publicly available package.


title: Capacity Results for the Wiretapped Oblivious Transfer

abstract: In this paper, we study the problem of the 1-of-2 string oblivious transfer (OT) between Alice and Bob in the presence of a passive eavesdropper Eve. The eavesdropper Eve is not allowed to get any information about the private data of Alice or Bob. When Alice and Bob are honest-but-curious users, we propose a protocol that satisfies $1$-private (neither Alice nor Bob colludes with Eve) OT requirements for the binary erasure symmetric broadcast channel, in which the channel provides dependent erasure patterns to Bob and Eve. We find that when the erasure probabilities satisfy certain conditions, the derived lower and upper bounds on the wiretapped OT capacity meet. Our results generalize and improve upon the results on wiretapped OT capacity by Mishra et al. Finally, we propose a protocol for a larger class of wiretapped channels and derive a lower bound on the wiretapped OT capacity.


title: BR-SNIS: Bias Reduced Self-Normalized Importance Sampling

abstract: Importance Sampling (IS) is a method for approximating expectations under a target distribution using independent samples from a proposal distribution and the associated importance weights. In many applications, the target distribution is known only up to a normalization constant, in which case self-normalized IS (SNIS) can be used. While the use of self-normalization can have a positive effect on the dispersion of the estimator, it introduces bias. In this work, we propose a new method, BR-SNIS, whose complexity is essentially the same as that of SNIS and which significantly reduces bias without increasing the variance. This method is a wrapper in the sense that it uses the same proposal samples and importance weights as SNIS, but makes clever use of iterated sampling--importance resampling (ISIR) to form a bias-reduced version of the estimator. We furnish the proposed algorithm with rigorous theoretical results, including new bias, variance and high-probability bounds, and these are illustrated by numerical examples.


title: Unexpected Biases between Congruence Classes for Parts in k-indivisible Partitions

abstract: For integers $k,t \geq 2$, and $1\leq r \leq t$ let $D_k^\times(r,t;n)$ be the number of parts among all $k$-indivisible partitions of $n$ (i.e., partitions where all parts are not divisible by $k$) of $n$ that are congruent to $r$ modulo $t$. Using Wright's circle method, we derive an asymptotic for $D_k^\times(r,t;n)$ as $n \to \infty$ when $k,t$ are coprime. The main term of this asymptotic does not depend on $r$, and so, in a weak asymptotic sense, the parts are equidistributed among congruence classes. However, inspection of the lower order terms indicates a bias towards different congruence classes modulo $t$. This induces an ordering on the congruence classes modulo $t$, which we call the $k$-indivisible ordering. We prove that for $k \geq \frac{6(t^2-1)}{\pi^2}$ the $k$-indivisible ordering matches the natural ordering. We also explore the properties of these orderings when $k < \frac{6(t^2-1)}{\pi^2}$.


title: N-Grammer: Augmenting Transformers with latent n-grams

abstract: Transformer models have recently emerged as one of the foundational models in natural language processing, and as a byproduct, there is significant recent interest and investment in scaling these models. However, the training and inference costs of these large Transformer language models are prohibitive, thus necessitating more research in identifying more efficient variants. In this work, we propose a simple yet effective modification to the Transformer architecture inspired by the literature in statistical language modeling, by augmenting the model with n-grams that are constructed from a discrete latent representation of the text sequence. We evaluate our model, the N-Grammer on language modeling on the C4 data-set as well as text classification on the SuperGLUE data-set, and find that it outperforms several strong baselines such as the Transformer and the Primer. We open-source our model for reproducibility purposes in Jax.


title: Configurational entropy of generalized sine-Gordon-type models

abstract: A family of deformed models of the sine-Gordon-type can be generated by twisting the sine-Gordon model. As a particular case, the 3-sine-Gordon model is here addressed, whose differential configurational entropy and the differential configurational complexity of three topological sectors are discussed, using two complementary approaches. Stability aspects are also discussed.


title: On The Babenko-Bechner-Type Inequality associated with the Weinstein Operator

abstract: In this paper, we study the Babenko-Bechner-type inequality for the Fourier Weinstein transform associated with the Weinstein operator. We use this inequality to establish a new version of Young's type inequality.


title: SmartPubSub: Content-based Pub-Sub on IPFS

abstract: The InterPlanetary File System (IPFS) is a hypermedia distribution protocol enabling the creation of completely distributed applications. One of the most efficient and effective ways to distribute information is through notifications, with a producer of content (publisher) sharing content with other interested parts (subscribers). IPFS already implements topic-based publish-subscribe systems under an experimental flag. The goal of this work is to advance on that, by developing a content-based pub-sub system (with subscriptions as predicates about event content) to disseminate information on top of IPFS in an efficient and decentralized way, leveraging its infrastructure. We design two protocols: ScoutSubs that is completely decentralized; FastDelivery that is centered in the publisher. With these two approaches, we show the different advantages of having each of these protocols simultaneously by comparing ScoutSubs full decentralization, and FastDelivery centralization at data sources.


title: A note on existence of an optimal set for a bonnesen type quantitative isoperimetric ratio in the plane

abstract: In this note we prove the existence of a set $E_0\subset\mathbb{R}^2$, different from a ball, which minimizes, among the convex sets that satisfy a suitable interior cone condition, the ratio \begin{equation} \label{eq:0} \frac{D(E)}{\lambda_\mathcal{H}^2(E)}, \end{equation} where $D$ is the isoperimetric deficit and $\lambda_\mathcal{H}$ the deviation from the spherical shape of a set $E\subset \mathbb{R}^2$.


title: Markovian Foundations for Quasi-Stochastic Approximation with Applications to Extremum Seeking Control

abstract: This paper concerns quasi-stochastic approximation (QSA) to solve root finding problems commonly found in applications to optimization and reinforcement learning. The general constant gain algorithm may be expressed as the time-inhomogeneous ODE $ \frac{d}{dt}\Theta_t=\alpha f_t (\Theta_t)$, with state process $\Theta$ evolving on $\mathbb{R}^d$. Theory is based on an almost periodic vector field, so that in particular the time average of $f_t(\theta)$ defines the time-homogeneous mean vector field $\bar{f} \colon \mathbb{R}^d \to \mathbb{R}^d$ with $\bar{f}(\theta^*)=0$. Under smoothness assumptions on the functions involved, the following exact representation is obtained: \[\frac{d}{dt}\Theta_t=\alpha[\bar{f}(\Theta_t)-\alpha\bar\Upsilon_t+\alpha^2\mathcal{W}_t^0+\alpha\frac{d}{dt}\mathcal{W}_t^1+\frac{d^2}{dt^2}\mathcal{W}_t^2]\] along with formulae for the smooth signals $\{\bar \Upsilon_t , \mathcal{W}_t^i : i=0, 1, 2\}$. This new representation, combined with new conditions for ultimate boundedness, has many applications for furthering the theory of QSA and its applications, including the following implications that are developed in this paper: (i) A proof that the estimation error $\|\Theta_t-\theta^*\|$ is of order $O(\alpha)$, but can be reduced to $O(\alpha^2)$ using a second order linear filter. (ii) In application to extremum seeking control, it is found that the results do not apply because the standard algorithms are not Lipschitz continuous. A new approach is presented to ensure that the required Lipschitz bounds hold, and from this we obtain stability, transient bounds, and asymptotic bias of order $O(\alpha^2)$, and asymptotic variance of order $O(\alpha^4)$. (iii) It is in general possible to obtain better than $O(\alpha)$ bounds on error in traditional stochastic approximation when there is Markovian noise.


title: The Monitoring Logging and Alarm System of the ASTRI Mini-Array gamma-ray air-Cherenkov experiment at the Observatorio del Teide

abstract: The ASTRI Mini-Array is a project for the Cherenkov astronomy in the TeV energy range. ASTRI Mini-Array consists of nine Imaging Atmospheric Cherenkov telescopes located at the Teide Observatory (Canarias Islands). Large volumes of monitoring and logging data result from the operation of a large-scale astrophysical observatory. In the last few years, several "Big Data" technologies have been developed to deal with such volumes of data, especially in the Internet of Things (IoT) framework. We present the Monitoring, Logging, and Alarm (MLA) system for the ASTRI Mini-Array aimed at supporting the analysis of scientific data and improving the operational activities of the telescope facility. The MLA system was designed and built considering the latest software tools and concepts coming from Big Data and IoT to respond to the challenges posed by the operation of the array. A particular relevance has been given to satisfying the reliability, availability, and maintainability requirements towards all the array sub-systems and auxiliary devices. The system architecture has been designed to scale up with the number of devices to be monitored and with the number of software components to be considered in the distributed logging system.


title: A fast method for solving the linear optimization p roblem subjected to simplified Dombi fuzzy relational equations

abstract: In this paper, an optimization model with a linear objective function subjected to a system of fuzzy relation equations (FRE) is studied where the feasible region is defined by the Dombi t-norm. Dombi family of t-norms includes a parametric family of continuous strict t-norms, whose members are increasing functions of the parameter. This family of t-norms covers the whole spectrum of t-norms when the parameter is changed from zero to infinity. Since the feasible solutions set of FREs is non-convex and the finding of all minimal solutions is an NP-hard problem, designing an efficient solution procedure for solving such problems is not a trivial job. Firstly, the feasible domain is characterized and then, based on some theoretical properties of the problem, a modified branch-and-bound solution technique is presented, which solves the problem by considering a few number of feasible paths. After presenting our solution procedure, a concrete example is included for illustration purpose.


title: Grassmannian packings: Trust-region stochastic tuning for matrix incoherence

abstract: We provide a new numerical procedure for constructing low coherence matrices, Trust-Region Stochastic Tuning for Matrix Incoherence (TRSTMI) and detail the results of experiments with a CPU/GPU parallelized implementation of this method. These trials suggest the superiority of this approach over other existing methods when the size of the matrix is large. We also present new conjectures on optimal complex matrices motivated and guided by the experimental results.


title: Affine Fractional Sobolev and Isoperimetric Inequalities

abstract: Sharp affine fractional Sobolev inequalities for functions on $\mathbb R^n$ are established. For each $0<s<1$, the new inequalities are significantly stronger than (and directly imply) the sharp fractional Sobolev inequalities of Almgren and Lieb. In the limit as $s\to 1^-$, the new inequalities imply the sharp affine Sobolev inequality of Gaoyong Zhang. As a consequence, fractional Petty projection inequalities are obtained that are stronger than the fractional Euclidean isoperimetric inequalities and a natural conjecture for radial mean bodies is proved.


title: Shock-avoiding slicing conditions: tests and calibrations

abstract: While the 1+log slicing condition has been extremely successful in numerous numerical relativity simulations, it is also known to develop "gauge-shocks" in some examples. Alternative "shock-avoiding" slicing conditions suggested by Alcubierre prevent these pathologies in those examples, but have not yet been explored and tested very broadly. In this paper we compare the performance of shock-avoiding slicing conditions with those of 1+log slicing for a number of "text-book" problems, including black holes and relativistic stars. While, in some simulations, the shock-avoiding slicing conditions feature some unusual properties and lead to more "gauge-dynamics" than the 1+log slicing condition, we find that they perform quite similarly in terms of stability and accuracy, and hence provide a very viable alternative to 1+log slicing.


title: Tilt-then-Blur or Blur-then-Tilt? Clarifying the Atmospheric Turbulence Model

abstract: Imaging at a long distance often requires advanced image restoration algorithms to compensate for the distortions caused by atmospheric turbulence. However, unlike many standard restoration problems such as deconvolution, the forward image formation model of the atmospheric turbulence does not have a simple expression. Thanks to the Zernike representation of the phase, one can show that the forward model is a combination of tilt (pixel shifting due to the linear phase terms) and blur (image smoothing due to the high order aberrations). Confusions then arise between the ordering of the two operators. Should the model be tilt-then-blur, or blur-then-tilt? Some papers in the literature say that the model is tilt-then-blur, whereas more papers say that it is blur-then-tilt. This paper clarifies the differences between the two and discusses why the tilt-then-blur is the correct model. Recommendations are given to the research community.


title: A fast algorithm for the linear programming problem constrained with the Weighted power mean-Fuzzy Relational Equalities (WPM-FRE)

abstract: In this paper, a linear programming problem is investigated in which the feasible region is formed as a special type of fuzzy relational equalities (FRE). In this type of FRE, fuzzy composition is considered as the weighted power mean operator (WPM). Some theoretical properties of the feasible region are derived and some necessary and sufficient conditions are also presented to determine the feasibility of the problem. Moreover, two procedures are proposed for simplifying the problem. Based on some structural properties of the problem, an algorithm is presented to find the optimal solutions and finally, an example is described to illustrate the algorithm.


title: Radar Image Reconstruction from Raw ADC Data using Parametric Variational Autoencoder with Domain Adaptation

abstract: This paper presents a parametric variational autoencoder-based human target detection and localization framework working directly with the raw analog-to-digital converter data from the frequency modulated continous wave radar. We propose a parametrically constrained variational autoencoder, with residual and skip connections, capable of generating the clustered and localized target detections on the range-angle image. Furthermore, to circumvent the problem of training the proposed neural network on all possible scenarios using real radar data, we propose domain adaptation strategies whereby we first train the neural network using ray tracing based model data and then adapt the network to work on real sensor data. This strategy ensures better generalization and scalability of the proposed neural network even though it is trained with limited radar data. We demonstrate the superior detection and localization performance of our proposed solution compared to the conventional signal processing pipeline and earlier state-of-art deep U-Net architecture with range-doppler images as inputs


title: Singularities of generic linkage via Frobenius powers

abstract: Let $I$ be an equidimensional ideal of a ring polynomial $R$ over $\mathbb{C}$ and let $J$ be its generic linkage. We prove that there is a uniform bound of the difference between the F-pure thresholds of $I_p$ and $J_p$ via the generalized Frobenius powers of ideals. This provides evidence that the F-pure threshold of an equidimensional ideal $I$ is less than that of its generic linkage. As a corollary we recover a result on log canonical thresholds of generic linkage by Niu.


title: Optimization of the storage database for the Monitoring system of the CTA

abstract: We present preliminary test results for the correct sizing of the bare metal hardware that will host the database of the Monitoring system (MON) for the Cherenkov Telescope Array (CTA). The MON is the subsystem of the Array Control and Data Acquisition System (ACADA) that is responsible for monitoring and logging the overall CTA array. It acquires and stores monitoring points and logging information from the array elements, at each of the CTA sites. MON is designed and built in order to deal with big data time series, and exploits some of the currently most advanced technologies in the fields of databases and Internet of Things (IoT). To dimension the bare metal hardware required by the monitoring system (excluding the logging), we performed the test campaign that is discussed in this paper. We discuss here the best set of parameters and the optimized configuration to maximize the database data writing in terms of the number of updated rows per second. We also demonstrate the feasibility of our approach in the frame of the CTA requirements.


title: On the generalised Springer correspondence for groups of type $E_8$

abstract: We complete the determination of the generalised Springer correspondence for connected reductive algebraic groups, by proving a conjecture of Lusztig on the last open cases which occur for groups of type $E_8$.


title: Exploration of 3D wavelet scattering transform coefficients for line-intensity mapping measurements

abstract: The wavelet scattering transform (WST) has recently gained attention in the context of large-scale structure studies, being a possible generator of summary statistics encapsulating non-Gaussianities beyond the reach of the conventional power spectrum. This work examines the three-dimensional solid harmonic WST in the context of a three-dimensional line-intensity mapping measurement to be undertaken by current and proposed phases of the CO Mapping Array Project (COMAP). The WST coefficients demonstrate interpretable behaviour in the context of noiseless CO line-intensity simulations. The contribution of the cosmological $z\sim3$ signal to these coefficients is also detectable in principle even in the Pathfinder phase of COMAP. Using the peak-patch method to generate large numbers of simulations and incorporating observational noise, we numerically estimate covariance matrices and show that careful choices of WST hyperparameters and rescaled or reduced coefficient sets are both necessary to keep covariances well-conditioned. Fisher forecasts show that even a reduced `shapeless' set of $\ell$-averaged WST coefficients show constraining power that can exceed that of the power spectrum alone even with similar detection significance. The full WST could improve parameter constraints even over the combination of the power spectrum and the voxel intensity distribution, showing that it uniquely encapsulates shape information about the line-intensity field. However, practical applications urgently require further understanding of the WST in key contexts like covariances and cross-correlations.


title: Trees with many leaves in tournaments

abstract: Sumner's universal tournament conjecture states that every $(2n-2)$-vertex tournament should contain a copy of every $n$-vertex oriented tree. If we know the number of leaves of an oriented tree, or its maximum degree, can we guarantee a copy of the tree with fewer vertices in the tournament? Due to work initiated by H\"aggkvist and Thomason (for number of leaves) and K\"uhn, Mycroft and Osthus (for maximum degree), it is known that improvements can be made over Sumner's conjecture in some cases, and indeed sometimes an $(n+o(n))$-vertex tournament may be sufficient. In this paper, we give new results on these problems. Specifically, we show i) for every $\alpha>0$, there exists $n_0\in\mathbb{N}$ such that, whenever $n\geqslant n_0$, every $((1+\alpha)n+k)$-vertex tournament contains a copy of every $n$-vertex oriented tree with $k$ leaves, and ii) for every $\alpha>0$, there exists $c>0$ and $n_0\in\mathbb{N}$ such that, whenever $n\geqslant n_0$, every $(1+\alpha)n$-vertex tournament contains a copy of every $n$-vertex oriented tree with maximum degree $\Delta(T)\leqslant cn$. Our first result gives an asymptotic form of a conjecture by Havet and Thomass\'e, while the second improves a result of Mycroft and Naia which applies to trees with polylogarithmic maximum degree.


title: Unit Reducible Fields and Perfect Unary Forms

abstract: In this paper, we introduce the notion of unit reducibility for number fields, that is, number fields in which all positive unary forms attain their nonzero minimum at a unit. Furthermore, we investigate the link between unit reducibility and the number of homothety classes of perfect unary forms for a given number field, and prove an open conjecture about the number of classes of perfect unary forms in real quadratic fields, stated by D. Yasaki.


title: Efficient geodesics in the curve complex and their dot graphs

abstract: For the complex of curves of a closed orientable surface of genus $g$, $\mathcal{C}(S_{g>1})$, the notion of efficient geodesic in was introduced in arXiv:1408.4133. There it was established that there always exists (finitely many) efficient geodesics between any two vertices, $ v_{\alpha} , v_{\beta} \in \mathcal{C}(S_g)$, representing homotopy classes of simple closed curves, $\alpha , \beta \subset S_g$. The main tool for used in establishing the existence of efficient geodesic was a dot graph, a booking scheme for recording the intersection pattern of a reference arc, $\gamma \subset S_g$, with the simple closed curves associated with the vertices of geodesic path in the zero skeleton, $\mathcal{C}^0(S_g)$. In particular, for an efficient geodesic between $v_\alpha $ and $v_\beta$ of length $d \geq 3$, it was shown that any curve corresponding to the vertex that is distance one from $v_\alpha$ intersects any $\gamma$ at most $d -2$ times. In this note we make a more expansive study of the characterizing "shape" of the dot graphs over the entire set of vertices in an efficient geodesic edge-path. The key take away of this study is that the shape of a dot graph for any efficient geodesic is contained within a spindle shape region. Since the Nielson-Thurston coordinates of any curve on $S_g$ are directly derived from its intersection number with finitely many reference arcs, spindle shaped dot graphs control the coordinate behavior of curves associated with the vertices of an efficient geodesic.


title: A General Framework for Partial to Full Image Registration

abstract: Image registration is a research field in which images must be compared and aligned independently of the point of view or camera characteristics. In some applications (such as forensic biometrics, satellite photography or outdoor scene identification) classical image registration systems fail due to one of the images compared represents a tiny piece of the other image. For instance, in forensics palmprint recognition, it is usual to find only a small piece of the palmprint, but in the database, the whole palmprint has been enrolled. The main reason of the poor behaviour of classical image registration methods is the gap between the amounts of salient points of both images, which is related to the number of points to be considered as outliers. Usually, the difficulty of finding a good match increases when the image that represents the tiny part of the scene has been drastically rotated. Again, in the case of palmprint forensics, it is difficult to decide a priori the orientation of the found tiny palmprint image. We present a rotation invariant registration method that explicitly considers that the image to be matched is a small piece of a larger image. We have experimentally validated our method in two different scenarios; palmprint identification and outdoor image registration.


title: River Surface Patch-wise Detector Using Mixture Augmentation for Scum-cover-index

abstract: Urban rivers provide a water environment that influences residential living. River surface monitoring has become crucial for making decisions about where to prioritize cleaning and when to automatically start the cleaning treatment. We focus on the organic mud, or "scum", that accumulates on the river's surface and contributes to the river's odor and has external economic effects on the landscape. Because of its feature of a sparsely distributed and unstable pattern of organic shape, automating the monitoring process has proved difficult. We propose a patch-wise classification pipeline to detect scum features on the river surface using mixture image augmentation to increase the diversity between the scum floating on the river and the entangled background on the river surface reflected by nearby structures like buildings, bridges, poles, and barriers. Furthermore, we propose a scum-index cover on rivers to help monitor worse grade online, collect floating scum, and decide on chemical treatment policies. Finally, we demonstrate the application of our method on a time series dataset with frames every ten minutes recording river scum events over several days. We discuss the significance of our pipeline and its experimental findings.


title: ProDiff: Progressive Fast Diffusion Model For High-Quality Text-to-Speech

abstract: Denoising diffusion probabilistic models (DDPMs) have recently achieved leading performances in many generative tasks. However, the inherited iterative sampling process costs hinder their applications to text-to-speech deployment. Through the preliminary study on diffusion model parameterization, we find that previous gradient-based TTS models require hundreds or thousands of iterations to guarantee high sample quality, which poses a challenge for accelerating sampling. In this work, we propose ProDiff, on progressive fast diffusion model for high-quality text-to-speech. Unlike previous work estimating the gradient for data density, ProDiff parameterizes the denoising model by directly predicting clean data to avoid distinct quality degradation in accelerating sampling. To tackle the model convergence challenge with decreased diffusion iterations, ProDiff reduces the data variance in the target site via knowledge distillation. Specifically, the denoising model uses the generated mel-spectrogram from an N-step DDIM teacher as the training target and distills the behavior into a new model with N/2 steps. As such, it allows the TTS model to make sharp predictions and further reduces the sampling time by orders of magnitude. Our evaluation demonstrates that ProDiff needs only 2 iterations to synthesize high-fidelity mel-spectrograms, while it maintains sample quality and diversity competitive with state-of-the-art models using hundreds of steps. ProDiff enables a sampling speed of 24x faster than real-time on a single NVIDIA 2080Ti GPU, making diffusion models practically applicable to text-to-speech synthesis deployment for the first time. Our extensive ablation studies demonstrate that each design in ProDiff is effective, and we further show that ProDiff can be easily extended to the multi-speaker setting. Audio samples are available at \url{https://ProDiff.github.io/.}


title: Dynamic Selection of Perception Models for Robotic Control

abstract: Robotic perception models, such as Deep Neural Networks (DNNs), are becoming more computationally intensive and there are several models being trained with accuracy and latency trade-offs. However, modern latency accuracy trade-offs largely report mean accuracy for single-step vision tasks, but there is little work showing which model to invoke for multi-step control tasks in robotics. The key challenge in a multi-step decision making is to make use of the right models at right times to accomplish the given task. That is, the accomplishment of the task with a minimum control cost and minimum perception time is a desideratum; this is known as the model selection problem. In this work, we precisely address this problem of invoking the correct sequence of perception models for multi-step control. In other words, we provide a provably optimal solution to the model selection problem by casting it as a multi-objective optimization problem balancing the control cost and perception time. The key insight obtained from our solution is how the variance of the perception models matters (not just the mean accuracy) for multi-step decision making, and to show how to use diverse perception models as a primitive for energy-efficient robotics. Further, we demonstrate our approach on a photo-realistic drone landing simulation using visual navigation in AirSim. Using our proposed policy, we achieved 38.04% lower control cost with 79.1% less perception time than other competing benchmarks.


title: Geometric vertex decomposition and liaison for toric ideals of graphs

abstract: The geometric vertex decomposability property for polynomial ideals is an ideal-theoretic generalization of the vertex decomposability property for simplicial complexes. Indeed, a homogeneous geometrically vertex decomposable ideal is radical and Cohen-Macaulay, and is in the Gorenstein liaison class of a complete intersection (glicci). In this paper, we initiate an investigation into when the toric ideal $I_G$ of a finite simple graph $G$ is geometrically vertex decomposable. We first show how geometric vertex decomposability behaves under tensor products, which allows us to restrict to connected graphs. We then describe a graph operation that preserves geometric vertex decomposability, thus allowing us to build many graphs whose corresponding toric ideals are geometrically vertex decomposable. Using work of Constantinescu and Gorla, we prove that toric ideals of bipartite graphs are geometrically vertex decomposable. We also propose a conjecture that all toric ideals of graphs with a square-free degeneration with respect to a lexicographic order are geometrically vertex decomposable. As evidence, we prove the conjecture in the case that the universal Gr\"obner basis of $I_G$ is a set of quadratic binomials. We also prove that some other families of graphs have the property that $I_G$ is glicci.


title: Relationship Design for Socially-Aware Behavior in Static Games

abstract: Autonomous agents can adopt socially-aware behaviors to reduce social costs, mimicking the way animals interact in nature and humans in society. We present a new approach to model socially-aware decision-making that includes two key elements: bounded rationality and inter-agent relationships. We capture the interagent relationships by introducing a novel model called a relationship game and encode agents' bounded rationality using quantal response equilibria. For each relationship game, we define a social cost function and formulate a mechanism design problem to optimize weights for relationships that minimize social cost at the equilibrium. We address the multiplicity of equilibria by presenting the problem in two forms: Min-Max and Min-Min, aimed respectively at minimization of the highest and lowest social costs in the equilibria. We compute the quantal response equilibrium by solving a least-squares problem defined with its Karush-Kuhn-Tucker conditions, and propose two projected gradient descent algorithms to solve the mechanism design problems. Numerical results, including two-lane congestion and congestion with an ambulance, confirm that these algorithms consistently reach the equilibrium with the intended social costs.


title: Fraisse Structures with SDAP+, Part I: Indivisibility

abstract: This is Part I of a two-part series regarding Ramsey properties of Fraisse structures satisfying a property called SDAP+, which strengthens the Disjoint Amalgamation Property. We prove that every Fraisse structure in a finite relational language with relation symbols of any finite arity satisfying this property is indivisible. Novelties include a new formulation of coding trees in terms of 1-types over initial segments of the Fraisse structure, and a direct proof of indivisibility which uses the method of forcing to conduct unbounded searches for finite sets. In Part II, we prove that every Fraisse structure in a finite relational language with relation symbols of arity at most two having this property has finite big Ramsey degrees which have a simple characterization. It follows that any such Fraisse structure admits a big Ramsey structure. Part II utilizes a theorem from Part I as a pigeonhole principle for induction arguments. This work offers a streamlined and unifying approach to Ramsey theory on some seemingly disparate classes of Fraisse structures.


title: Confining complex ghost degrees of freedom

abstract: We show a theorem proving that a non-local bosonic field upon a covariant interaction with a confining gauge field undergoes the confinement of its degrees of freedom present in the free theory changing completely the physical mass spectrum following Kugo-Ojima criterion. This is applicable to an infinite number of excitations of the bosonic field including ghosts whereas we pay special attention to the modes with the complex conjugate masses, states appearing in the string field theory motivated infinite-derivative models. The same recipe will obviously work for the Lee-Wick models.


title: Rough homogenization for Langevin dynamics on fluctuating Helfrich surfaces

abstract: In this paper, we study different scaling rough path limit regimes in space and time for the Langevin dynamics on a quasi-planar fluctuating Helfrich surfaces. The convergence results of the processes were already proven in the work by Duncan, Elliott, Pavliotis and Stuart (2015). We extend this work by proving the convergence of the It\^o and Stratonovich rough path lift. For the rough path limit, there appears, typically, an area correction term to the It\^o iterated integrals, and in certain regimes to the Stratonovich iterated integrals. This yields additional information on the homogenization limit and enables to conclude on homogenization results for diffusions driven by the Brownian motion on the membrane using the continuity of the It\^o-Lyons map in rough paths topology.


title: On Market Clearing of Day Ahead Auctions for European Power Markets: Cost Minimisation versus Social Welfare Maximisation

abstract: For the case of inflexible demand and considering network constraints, we introduce a Cost Minimisation (CM) based market clearing mechanism, and a model representing the standard Social Welfare Maximisation mechanism used in European Day Ahead Electricity Markets. Since the CM model corresponds to a more challenging optimisation problem, we propose four numerical algorithms that leverage the problem structure, each with different trade-offs between computational cost and convergence guarantees. These algorithms are evaluated on synthetic data to provide some intuition of their performance. We also provide strong (but partial) analytical results to facilitate efficient solution of the CM problem, which call for the introduction of a new concept: optimal zonal stack curves, and these results are used to devise one of the four solution algorithms. An evaluation of the CM and SWM models and their comparison is performed, under the assumption of truthful bidding, on the real world data of Central Western European Day Ahead Power Market during the period of 2019-2020. We show that the SWM model we introduce gives a good representation of the historical time series of the real prices. Further, the CM reduces the market power of producers, as generally this results in decreased zonal prices and always decreases the total cost of electricity procurement when compared to the currently employed SWM.


title: Quantum state tomography with tensor train cross approximation

abstract: It has been recently shown that a state generated by a one-dimensional noisy quantum computer is well approximated by a matrix product operator with a finite bond dimension independent of the number of qubits. We show that full quantum state tomography can be performed for such a state with a minimal number of measurement settings using a method known as tensor train cross approximation. The method works for reconstructing full rank density matrices and only requires measuring local operators, which are routinely performed in state-of-art experimental quantum platforms. Our method requires exponentially fewer state copies than the best known tomography method for unstructured states and local measurements. The fidelity of our reconstructed state can be further improved via supervised machine learning, without demanding more experimental data. Scalable tomography is achieved if the full state can be reconstructed from local reductions.


title: A numerical study of stellar discs in galactic nuclei

abstract: We explore the dynamics of stellar discs in the close vicinity of a supermassive black hole (SMBH) by means of direct $N$-body simulations. We show that an isolated nuclear stellar disc exhibits anisotropic mass segregation meaning that massive stars settle to lower orbital inclinations and more circular orbits than the light stars. However, in systems in which the stellar disc is embedded in a much more massive isotropic stellar cluster, anisotropic mass segregation tends to be suppressed. In both cases, an initially thin stellar disc becomes thicker, especially in the inner parts due to the fluctuating anisotropy in the spherical component. We find that vector resonant relaxation is quenched in the disc by nodal precession, but it is still the most efficient relaxation process around SMBHs of mass $10^6M_\odot$ and above. Two body relaxation may dominate for less massive SMBHs found in dwarf galaxies. Stellar discs embedded in massive isotropic stellar clusters ultimately tend to become isotropic on the local two-body relaxation time-scale. Our simulations show that the dynamics of young stars at the centre of the Milky Way is mostly driven by vector resonant relaxation leading to an anticorrelation between the scatter of orbital inclinations and distance from the SMBH. If the $S$-stars formed in a disc less than 10 Myr ago, they may coexist with a cusp of stellar mass black holes or an intermediate mass black hole with mass up to $1000M_\odot$ to reproduce the observed scatter of angular momenta.


title: Pattern recognition in the nucleation kinetics of non-equilibrium self-assembly

abstract: Inspired by biology's most sophisticated computer, the brain, neural networks constitute a profound reformulation of computational principles. Remarkably, analogous high-dimensional, highly-interconnected computational architectures also arise within information-processing molecular systems inside living cells, such as signal transduction cascades and genetic regulatory networks. Might neuromorphic collective modes be found more broadly in other physical and chemical processes, even those that ostensibly play non-information-processing roles such as protein synthesis, metabolism, or structural self-assembly? Here we examine nucleation during self-assembly of multicomponent structures, showing that high-dimensional patterns of concentrations can be discriminated and classified in a manner similar to neural network computation. Specifically, we design a set of 917 DNA tiles that can self-assemble in three alternative ways such that competitive nucleation depends sensitively on the extent of co-localization of high-concentration tiles within the three structures. The system was trained in-silico to classify a set of 18 grayscale 30 x 30 pixel images into three categories. Experimentally, fluorescence and atomic force microscopy monitoring during and after a 150-hour anneal established that all trained images were correctly classified, while a test set of image variations probed the robustness of the results. While slow compared to prior biochemical neural networks, our approach is surprisingly compact, robust, and scalable. This success suggests that ubiquitous physical phenomena, such as nucleation, may hold powerful information processing capabilities when scaled up as high-dimensional multicomponent systems.


title: PyMAF-X: Towards Well-aligned Full-body Model Regression from Monocular Images

abstract: We present PyMAF-X, a regression-based approach to recovering parametric full-body models from monocular images. This task is very challenging since minor parametric deviation may lead to noticeable misalignment between the estimated mesh and the input image. Moreover, when integrating part-specific estimations into the full-body model, existing solutions tend to either degrade the alignment or produce unnatural wrist poses. To address these issues, we propose a Pyramidal Mesh Alignment Feedback (PyMAF) loop in our regression network for well-aligned human mesh recovery and extend it as PyMAF-X for the recovery of expressive full-body models. The core idea of PyMAF is to leverage a feature pyramid and rectify the predicted parameters explicitly based on the mesh-image alignment status. Specifically, given the currently predicted parameters, mesh-aligned evidence will be extracted from finer-resolution features accordingly and fed back for parameter rectification. To enhance the alignment perception, an auxiliary dense supervision is employed to provide mesh-image correspondence guidance while spatial alignment attention is introduced to enable the awareness of the global contexts for our network. When extending PyMAF for full-body mesh recovery, an adaptive integration strategy is proposed in PyMAF-X to produce natural wrist poses while maintaining the well-aligned performance of the part-specific estimations. The efficacy of our approach is validated on several benchmark datasets for body, hand, face, and full-body mesh recovery, where PyMAF and PyMAF-X effectively improve the mesh-image alignment and achieve new state-of-the-art results. The project page with code and video results can be found at https://www.liuyebin.com/pymaf-x.


title: Mechanical overtone frequency combs

abstract: Mechanical frequency combs are poised to bring the applications and utility of optical frequency combs into the mechanical domain. So far, their main challenge has been strict requirements on drive frequencies and power, which complicate operation. We demonstrate a straightforward mechanism to create a frequency comb consisting of mechanical overtones (integer multiples) of a single eigenfrequency, by monolithically integrating a suspended dielectric membrane with a counter-propagating optical trap. The periodic optical field modulates the dielectrophoretic force on the membrane at the overtones of a membrane's motion. These overtones share a fixed frequency and phase relation, and constitute a mechanical frequency comb. The periodic optical field also creates an optothermal parametric drive that requires no additional power or external frequency reference. This combination of effects results in an easy-to-use mechanical frequency comb platform that requires no precise alignment, no additional feedback or control electronics, and only uses a single, mW continuous wave laser beam. This highlights the overtone frequency comb as the straightforward future for applications in sensing, metrology and quantum acoustics.


title: How to make and trap a pseudo-vesicle with a micropipette

abstract: We present a simple method to produce giant lipid pseudo-vesicles (vesicles with an oily cap on the top), trapped in an agarose gel. The method can be implemented using only a regular micropipette and relies on the formation of a water/oil/water double droplet in liquid agarose. We characterize the produced vesicle with fluorescence imaging and establish the presence and integrity of the lipid bilayer by the successful insertion of {\alpha}-Hemolysin transmembrane proteins. Finally, we show that the vesicle can be easily mechanically deformed, non-intrusively, by indenting the surface of the gel.