title: Constructing invariant tori using guaranteed Euler method

abstract: We show here how, using Euler's integration method and an associated function bounding the error in function of time, one can generate structures closely surrounding the invariant tori of dynamical systems. Such structures are constructed from a finite number of balls of $\mathbb{R}^n$ and encompass the deformations of the tori when small perturbations of the flow of the system occur.


title: Meta-Learning to Compositionally Generalize

abstract: Natural language is compositional; the meaning of a sentence is a function of the meaning of its parts. This property allows humans to create and interpret novel sentences, generalizing robustly outside their prior experience. Neural networks have been shown to struggle with this kind of generalization, in particular performing poorly on tasks designed to assess compositional generalization (i.e. where training and testing distributions differ in ways that would be trivial for a compositional strategy to resolve). Their poor performance on these tasks may in part be due to the nature of supervised learning which assumes training and testing data to be drawn from the same distribution. We implement a meta-learning augmented version of supervised learning whose objective directly optimizes for out-of-distribution generalization. We construct pairs of tasks for meta-learning by sub-sampling existing training data. Each pair of tasks is constructed to contain relevant examples, as determined by a similarity metric, in an effort to inhibit models from memorizing their input. Experimental results on the COGS and SCAN datasets show that our similarity-driven meta-learning can improve generalization performance.


title: A likelihood based sensitivity analysis for publication bias on summary ROC in meta-analysis of diagnostic test accuracy

abstract: In meta-analysis of diagnostic test accuracy, summary receiver operating characteristic (SROC) is a recommended method to summarize the discriminant capacity of a diagnostic test in the presence of study-specific cutoff values and the area under the SROC (SAUC) gives the aggregate measure of test accuracy. SROC or SAUC can be estimated by bivariate modelling of pairs of sensitivity and specificity over the primary diagnostic studies. However, publication bias is a major threat to the validity of estimates in meta-analysis. To address this issue, we propose to adopt sensitivity analysis to make an objective inference for the impact of publication bias on SROC or SAUC. We extend Copas likelihood based sensitivity analysis to the bivariate normal model used for meta-analysis of diagnostic test accuracy to evaluate how much SROC or SAUC would change with different selection probabilities under several selective publication mechanisms dependent on sensitivity and/or specificity. The selection probability is modelled by a selection function on $t$-type statistic for the linear combination of logit-transformed sensitivity and specificity, allowing the selective publication of each study to be influenced by the cutoff-dependent $p$-value for sensitivity, specificity, or diagnostic odds ratio. By embedding the selection function into the bivariate normal model, the conditional likelihood is proposed and the bias-corrected SROC or SAUC can be estimated by maximizing the likelihood. We illustrate the proposed sensitivity analysis by reanalyzing a meta-analysis of test accuracy for intravascular device related infection. Simulation studies are conducted to investigate the performance of proposed methods.


title: Coresets for Classification -- Simplified and Strengthened

abstract: We give relative error coresets for training linear classifiers with a broad class of loss functions, including the logistic loss and hinge loss. Our construction achieves $(1\pm \epsilon)$ relative error with $\tilde O(d \cdot \mu_y(X)^2/\epsilon^2)$ points, where $\mu_y(X)$ is a natural complexity measure of the data matrix $X \in \mathbb{R}^{n \times d}$ and label vector $y \in \{-1,1\}^n$, introduced in by Munteanu et al. 2018. Our result is based on subsampling data points with probabilities proportional to their $\ell_1$ $Lewis$ $weights$. It significantly improves on existing theoretical bounds and performs well in practice, outperforming uniform subsampling along with other importance sampling methods. Our sampling distribution does not depend on the labels, so can be used for active learning. It also does not depend on the specific loss function, so a single coreset can be used in multiple training scenarios.


title: Nonparametric Regression for 3D Point Cloud Learning

abstract: Over the past two decades, we have seen an exponentially increased amount of point clouds collected with irregular shapes in various areas. Motivated by the importance of solid modeling for point clouds, we develop a novel and efficient smoothing tool based on multivariate splines over the tetrahedral partitions to extract the underlying signal and build up a 3D solid model from the point cloud. The proposed smoothing method can denoise or deblur the point cloud effectively and provide a multi-resolution reconstruction of the actual signal. In addition, it can handle sparse and irregularly distributed point clouds and recover the underlying trajectory. The proposed smoothing and interpolation method also provides a natural way of numerosity data reduction. Furthermore, we establish the theoretical guarantees of the proposed method. Specifically, we derive the convergence rate and asymptotic normality of the proposed estimator and illustrate that the convergence rate achieves the optimal nonparametric convergence rate. Through extensive simulation studies and a real data example, we demonstrate the superiority of the proposed method over traditional smoothing methods in terms of estimation accuracy and efficiency of data reduction.


title: Production cross sections of $^{3,4}_\Lambda$H bound states in $^{3,4}$He($K^-$, $\pi^0$) reactions at 1 GeV/$c$

abstract: We investigate theoretically productions of $^{3,4}_\Lambda$H bound states in the exothermic ($K^-$, $\pi^0$) reactions on $^{3,4}$He targets at $p_{K^-}=$ 1.0 GeV/c in a distorted-wave impulse approximation with the optimal Fermi-averaging $K^-p\to\pi^0\Lambda$ $t$ matrix. We calculate angular distributions of the laboratory differential cross sections $d\sigma/d\Omega_{\rm lab}$ and the integrated cross sections $\sigma_{\rm lab}$ for a $J^P= 0^+$ ground state and a $J^P=1^+$ excited state of $^4_\Lambda$H at $\pi^0$ forward-direction angles of $\theta_{\rm lab}=$ 0$^\circ$-20$^\circ$, and those for a $J^P= 1/2^+$ ground state of $^3_\Lambda$H. The production of a $J^P=3/2^+$ excited state of $^3_\Lambda$H as a virtual state is also evaluated. The comparison in $d\sigma/d\Omega_{\rm lab}$ and $\sigma_{\rm lab}$ between $^4_\Lambda$H and $^3_\Lambda$H provides examining the mechanism of the production and structure of $^{3,4}_\Lambda$H, as well as in the endothermic ($\pi^-$, $K^0$) reactions at $p_{\pi^-}=$ 1.05 GeV/c. This investigation confirms the feasibility of lifetime measurements of $^3_\Lambda$H at the J-PARC experiments.


title: Augmenting Teleportation in Virtual Reality With Discrete Rotation Angles

abstract: Locomotion is one of the most essential interaction tasks in virtual reality (VR) with teleportation being widely accepted as the state-of-the-art locomotion technique at the time of this writing. A major draw-back of teleportation is the accompanying physical rotation that is necessary to adjust the users' orientation either before or after teleportation. This is a limiting factor for tethered head-mounted displays (HMDs) and static body postures and can induce additional simulator sickness for HMDs with three degrees-of-freedom (DOF) due to missing parallax cues. To avoid physical rotation, previous work proposed discrete rotation at fixed intervals (InPlace) as a controller-based technique with low simulator sickness, yet the impact of varying intervals on spatial disorientation, user presence and performance remains to be explored. An unevaluated technique found in commercial VR games is reorientation during the teleportation process (TeleTurn), which prevents physical rotation but potentially increases interaction time due to its continuous orientation selection. In an exploratory user study, where participants were free to apply both techniques, we evaluated the impact of rotation parameters of either technique on user performance and preference. Our results indicate that discrete InPlace rotation introduced no significant spatial disorientation, while user presence scores were increased. Discrete TeleTurn and teleportation without rotation was ranked higher and achieved a higher presence score than continuous TeleTurn, which is the current state-of-the-art found in VR games. Based on observations, that participants avoided TeleTurn rotation when discrete InPlace rotation was available, we distilled guidelines for designing teleportation without physical rotation.


title: Interpretable agent communication from scratch (with a generic visual processor emerging on the side)

abstract: As deep networks begin to be deployed as autonomous agents, the issue of how they can communicate with each other becomes important. Here, we train two deep nets from scratch to perform realistic referent identification through unsupervised emergent communication. We show that the largely interpretable emergent protocol allows the nets to successfully communicate even about object types they did not see at training time. The visual representations induced as a by-product of our training regime, moreover, show comparable quality, when re-used as generic visual features, to a recent self-supervised learning model. Our results provide concrete evidence of the viability of (interpretable) emergent deep net communication in a more realistic scenario than previously considered, as well as establishing an intriguing link between this field and self-supervised visual learning.


title: Spectro-imagery of an active tornado-like prominence: formation and evolution

abstract: The nature of flows in tornado-prominences is an open issue. While the AIA imager aboard the Solar Dynamics Observatory (SDO) allowed us to follow the global structure of a tornado-like prominence during five hours, the Interface Region Imaging Spectrograph (IRIS), and the Multi subtractive Double pass spectrograph (MSDP) permitted to obtain plasma diagnostics of its fine structures. We aim to address two questions. Is the observed plasma rotation conceptually acceptable in a flux rope magnetic support configuration with dips? How is the plasma density distributed in the tornado-like prominence? We calculated line-of-sight velocities and non-thermal line widths using Gaussian fitting for Mg II lines and bisector method for H-alpha line. We determined the electron density from Mg II line integrated intensities and profile fitting methods using 1D NLTE radiative transfer theory models. The global structure of the prominence observed in H-alpha, and Mg II h and k lines fits with a magnetic field structure configuration with dips. Coherent Dopplershifts in red- and blue-shifted areas observed in both lines were detected along rapidly-changing vertical and horizontal structures. However, the tornado at the top of the prominence consists of multiple-fine threads with opposite flows suggesting counter streaming flows rather than rotation. Surprisingly we found that the electron density at the top of the prominence could be larger (10^11 cm^{-3}) than in the inner part of the prominence. We suggest that the tornado is in a formation state with cooling of hot plasma in a first phase, and following that, a phase of leakage of the formed blobs with large transverse flows of material along long loops extended away of the UV prominence top. The existence of such long magnetic field lines on both sides of the prominence would avoid the tornado-like prominence to really turn around its axis.


title: Provably Robust Detection of Out-of-distribution Data (almost) for free

abstract: The application of machine learning in safety-critical systems requires a reliable assessment of uncertainty. However, deep neural networks are known to produce highly overconfident predictions on out-of-distribution (OOD) data. Even if trained to be non-confident on OOD data, one can still adversarially manipulate OOD data so that the classifier again assigns high confidence to the manipulated samples. We show that two previously published defenses can be broken by better adapted attacks, highlighting the importance of robustness guarantees around OOD data. Since the existing method for this task is hard to train and significantly limits accuracy, we construct a classifier that can simultaneously achieve provably adversarially robust OOD detection and high clean accuracy. Moreover, by slightly modifying the classifier's architecture our method provably avoids the asymptotic overconfidence problem of standard neural networks. We provide code for all our experiments.


title: Minimal Gaps and Additive Energy in real-valued sequences

abstract: We study the minimal gap statistic for sequences of the form $\left( \alpha x_n \right)_{n = 1}^{\infty}$ where $\left( x_n \right)_{n = 1}^{\infty}$ is a sequence of real numbers, and its connection to the additive energy of $\left( x_n \right)_{n = 1}^{\infty}$. Inspired by a recent paper of Aistleitner, El-Baz and Munsch we show conditionally on the Lindel\"{o}f Hypothesis that if the additive energy is of lowest possible order then for almost all $\alpha$, the minimal gap $\delta_{\min}^{\alpha} (N) = \min \left\{ \alpha x_m - \alpha x_n \bmod \ 1 : 1 \leq m \neq n \leq N \right\}$ is close to that of a random sequence, a result Rudnick showed for integer-valued sequences. We also show unconditional results in this direction, as well as some converse theorems about sequences with large additive energy.


title: Question Generation for Adaptive Education

abstract: Intelligent and adaptive online education systems aim to make high-quality education available for a diverse range of students. However, existing systems usually depend on a pool of hand-made questions, limiting how fine-grained and open-ended they can be in adapting to individual students. We explore targeted question generation as a controllable sequence generation task. We first show how to fine-tune pre-trained language models for deep knowledge tracing (LM-KT). This model accurately predicts the probability of a student answering a question correctly, and generalizes to questions not seen in training. We then use LM-KT to specify the objective and data for training a model to generate questions conditioned on the student and target difficulty. Our results show we succeed at generating novel, well-calibrated language translation questions for second language learners from a real online education platform.


title: On the Connection between Local Attention and Dynamic Depth-wise Convolution

abstract: Vision Transformer (ViT) attains state-of-the-art performance in visual recognition, and the variant, Local Vision Transformer, makes further improvements. The major component in Local Vision Transformer, local attention, performs the attention separately over small local windows. We rephrase local attention as a channel-wise locally-connected layer and analyze it from two network regularization manners, sparse connectivity and weight sharing, as well as weight computation. Sparse connectivity: there is no connection across channels, and each position is connected to the positions within a small local window. Weight sharing: the connection weights for one position are shared across channels or within each group of channels. Dynamic weight: the connection weights are dynamically predicted according to each image instance. We point out that local attention resembles depth-wise convolution and its dynamic version in sparse connectivity. The main difference lies in weight sharing - depth-wise convolution shares connection weights (kernel weights) across spatial positions. We empirically observe that the models based on depth-wise convolution and the dynamic variant with lower computation complexity perform on-par with or sometimes slightly better than Swin Transformer, an instance of Local Vision Transformer, for ImageNet classification, COCO object detection and ADE semantic segmentation. These observations suggest that Local Vision Transformer takes advantage of two regularization forms and dynamic weight to increase the network capacity. Code is available at https://github.com/Atten4Vis/DemystifyLocalViT.


title: DKP oscillator in the presence of a spinning cosmic string

abstract: We investigate the modification of gravitational fields generated by topological defects on a generalized Duffin-Kemmer-Petiau (DKP) oscillator for spin-0 particle under spinning cosmic string background. The generalized DKP oscillator equation under spinning cosmic string background is established, and the impact of the Cornell potential on the generalized DKP oscillator is presented. We give the influence of space-time and potential parameters on energy levels.


title: Towards Social Role-Based Interruptibility Management

abstract: Pervasive and ubiquitous computing facilitates immediate access to information in the sense of always-on. Information such as news, messages, or reminders can significantly enhance our daily routines but are rendered useless or disturbing when not being aligned with our intrinsic interruptibility preferences. Attention management systems use machine learning to identify short-term opportune moments, so that information delivery leads to fewer interruptions. Humans' intrinsic interruptibility preferences - established for and across social roles and life domains - would complement short-term attention and interruption management approaches. In this article, we present our comprehensive results towards social role-based attention and interruptibility management. Our approach combines on-device sensing and machine learning with theories from social science to form a personalized two-stage classification model. Finally, we discuss the challenges of the current and future AI-driven attention management systems concerning privacy, ethical issues, and future directions.


title: Evolution of the galaxy merger fraction in the CLAUDS+HSC-SSP deep fields

abstract: We estimate the evolution of the galaxy-galaxy merger fraction for $M_\star>10^{10.5}M_\odot$ galaxies over $0.25<z<1$ in the $\sim$18.6 deg$^2$ deep CLAUDS+HSC-SSP surveys. We do this by training a Random Forest Classifier to identify merger candidates from a host of parametric morphological features, and then visually follow-up likely merger candidates to reach a high-purity, high-completeness merger sample. Correcting for redshift-dependent detection bias, we find that the merger fraction at $z=0$ is 1.0$\pm$0.2%, that the merger fraction evolves as $(1+z)^{2.3 \pm 0.4}$, and that a typical massive galaxy has undergone $\sim$0.3 major mergers since $z=1$. This pilot study illustrates the power of very deep ground-based imaging surveys combined with machine learning to detect and study mergers through the presence of faint, low surface brightness merger features out to at least $z\sim1$.


title: Supervised Machine Learning with Plausible Deniability

abstract: We study the question of how well machine learning (ML) models trained on a certain data set provide privacy for the training data, or equivalently, whether it is possible to reverse-engineer the training data from a given ML model. While this is easy to answer negatively in the most general case, it is interesting to note that the protection extends over non-recoverability towards plausible deniability: Given an ML model $f$, we show that one can take a set of purely random training data, and from this define a suitable ``learning rule'' that will produce a ML model that is exactly $f$. Thus, any speculation about which data has been used to train $f$ is deniable upon the claim that any other data could have led to the same results. We corroborate our theoretical finding with practical examples, and open source implementations of how to find the learning rules for a chosen set of raining data.


title: Silting complexes of coherent sheaves and the Humphreys conjecture

abstract: Let $G$ be a connected reductive algebraic group over an algebraically closed field $\Bbbk$ of characteristic $p \ge 0$, and let $\mathcal{N}$ be its nilpotent cone. Under mild hypotheses, we construct for each nilpotent $G$-orbit $C$ and each indecomposable tilting vector bundle $T$ on $C$ a certain complex $S(C,T)$ of $G \times \mathbb{G}_m$-equivariant coherent sheaves on $\mathcal{N}$. We prove that these objects are (up to shift) precisely the indecomposable objects in the coheart of a certain co-$t$-structure. We then show that if $p$ is larger than the Coxeter number, then the hypercohomology $H^\bullet(S(C,T))$ is identified with the cohomology of a tilting module for $G$. This confirms a conjecture of Humphreys on the support of the cohomology of tilting modules.


title: HPRNet: Hierarchical Point Regression for Whole-Body Human Pose Estimation

abstract: In this paper, we present a new bottom-up one-stage method for whole-body pose estimation, which we call "hierarchical point regression," or HPRNet for short. In standard body pose estimation, the locations of $\sim 17$ major joints on the human body are estimated. Differently, in whole-body pose estimation, the locations of fine-grained keypoints (68 on face, 21 on each hand and 3 on each foot) are estimated as well, which creates a scale variance problem that needs to be addressed. To handle the scale variance among different body parts, we build a hierarchical point representation of body parts and jointly regress them. The relative locations of fine-grained keypoints in each part (e.g. face) are regressed in reference to the center of that part, whose location itself is estimated relative to the person center. In addition, unlike the existing two-stage methods, our method predicts whole-body pose in a constant time independent of the number of people in an image. On the COCO WholeBody dataset, HPRNet significantly outperforms all previous bottom-up methods on the keypoint detection of all whole-body parts (i.e. body, foot, face and hand); it also achieves state-of-the-art results on face (75.4 AP) and hand (50.4 AP) keypoint detection. Code and models are available at \url{https://github.com/nerminsamet/HPRNet}.


title: Conelike radiant structures

abstract: Analogues of the classical affine-projective correspondence are developed in the context of statistical manifolds compatible with a radiant vector field. These utilize a formulation of Einstein equations for special statistical structures that generalizes the usual Einstein equations for pseudo-Riemannian metrics and is of independent interest. A conelike radiant structure is a not necessarily flat affine connection equipped with a family of surfaces that behave like the intersections of the planes through the origin with a convex cone in a real vector space. A radiant structure is a torsion-free affine connection and a vector field whose covariant derivative is the identity endomorphism. A radiant structure is conelike if for every point and every two-dimensional subspace containing the radiant vector field there is a totally geodesic surface tangent to the subspace at the point. Such structures exist on the total space of any principal bundle with one-dimensional fiber and on any Lie group with a quadratic structure on its Lie algebra. The affine connection of a conelike radiant structure can be canonically normalized to have antisymmetric Ricci tensor. The compatibility of radiant and conelike structures with metrics is investigated and yields connections for which the symmetrized Ricci curvature is a constant multiple of a compatible metric that generalizes well-known constructions of Riemannian and Lorentzian Einstein-Weyl structures over K\"ahler-Einstein manifolds having nonzero scalar curvature. Einstein equations for special statistical manifolds are given that generalize the Einstein-Weyl equations and encompass these more general examples. There are constructed left-invariant conelike radiant structures on a Lie group endowed with a left-invariant nondegenerate bilinear form, and the case of three-dimensional unimodular Lie groups is described in detail.


title: Inference for Network Regression Models with Community Structure

abstract: Network regression models, where the outcome comprises the valued edge in a network and the predictors are actor or dyad-level covariates, are used extensively in the social and biological sciences. Valid inference relies on accurately modeling the residual dependencies among the relations. Frequently homogeneity assumptions are placed on the errors which are commonly incorrect and ignore critical, natural clustering of the actors. In this work, we present a novel regression modeling framework that models the errors as resulting from a community-based dependence structure and exploits the subsequent exchangeability properties of the error distribution to obtain parsimonious standard errors for regression parameters.


title: Quasi-plurisubharmonic envelopes 2: Bounds on Monge-Amp\`ere volumes

abstract: In \cite{GL21a} we have developed a new approach to $L^{\infty}$-a priori estimates for degenerate complex Monge-Amp\`ere equations, when the reference form is closed. This simplifying assumption was used to ensure the constancy of the volumes of Monge-Amp\`ere measures. We study here the way these volumes stay away from zero and infinity when the reference form is no longer closed. We establish a transcendental version of the Grauert-Riemenschneider conjecture, partially answering conjectures of Demailly-P\u{a}un \cite{DP04} and Boucksom-Demailly-P\u{a}un-Peternell \cite{BDPP13}. Our approach relies on a fine use of quasi-plurisubharmonic envelopes. The results obtained here will be used in \cite{GL21b} for solving degenerate complex Monge-Amp\`ere equations on compact Hermitian varieties.


title: Quasi-plurisubharmonic envelopes 1: Uniform estimates on K\"ahler manifolds

abstract: We develop a new approach to $L^{\infty}$-a priori estimates for degenerate complex Monge-Amp\`ere equations on complex manifolds. It only relies on compactness and envelopes properties of quasi-plurisubharmonic functions. Our method allows one to obtain new and efficient proofs of several fundamental results in K\"ahler geometry as we explain in this article. In a sequel we shall explain how this approach also applies to the hermitian setting producing new relative a priori bounds, as well as existence results.


title: 3D Human Pose Estimation Based on 2D-3D Consistency with Synchronized Adversarial Training

abstract: 3D human pose estimation from a single image is still a challenging problem despite the large amount of work that has been performed in this field. Generally, most methods directly use neural networks and ignore certain constraints (e.g., reprojection constraints, joint angle, and bone length constraints). While a few methods consider these constraints but train the network separately, they cannot effectively solve the depth ambiguity problem. In this paper, we propose a GAN-based model for 3D human pose estimation, in which a reprojection network is employed to learn the mapping of the distribution from 3D poses to 2D poses, and a discriminator is employed for 2D-3D consistency discrimination. We adopt a novel strategy to synchronously train the generator, the reprojection network and the discriminator. Furthermore, inspired by the typical kinematic chain space (KCS) matrix, we introduce a weighted KCS matrix and take it as one of the discriminator's inputs to impose joint angle and bone length constraints. The experimental results on Human3.6M show that our method significantly outperforms state-of-the-art methods in most cases.


title: Raw Waveform Encoder with Multi-Scale Globally Attentive Locally Recurrent Networks for End-to-End Speech Recognition

abstract: End-to-end speech recognition generally uses hand-engineered acoustic features as input and excludes the feature extraction module from its joint optimization. To extract learnable and adaptive features and mitigate information loss, we propose a new encoder that adopts globally attentive locally recurrent (GALR) networks and directly takes raw waveform as input. We observe improved ASR performance and robustness by applying GALR on different window lengths to aggregate fine-grain temporal information into multi-scale acoustic features. Experiments are conducted on a benchmark dataset AISHELL-2 and two large-scale Mandarin speech corpus of 5,000 hours and 21,000 hours. With faster speed and comparable model size, our proposed multi-scale GALR waveform encoder achieved consistent character error rate reductions (CERRs) from 7.9% to 28.1% relative over strong baselines, including Conformer and TDNN-Conformer. In particular, our approach demonstrated notable robustness than the traditional handcrafted features and outperformed the baseline MFCC-based TDNN-Conformer model by a 15.2% CERR on a music-mixed real-world speech test set.


title: Slug bubble growth and dissolution by solute exchange

abstract: In many environmental and industrial applications, the mass transfer of gases in liquid solvents is a fundamental process during the generation of bubbles for specific purposes or, vice versa, the removal of entrapped bubbles. We address the growth dynamics of a trapped slug bubble in a vertical glass cylinder under a water barrier. In the studied process, the ambient air atmosphere is replaced by a CO$_2$ atmosphere at the same or higher pressure. The asymmetric exchange of the gaseous solutes between the CO$_2$-rich water barrier and the air-rich bubble always results in net bubble growth. We refer to this process as solute exchange. The dominant transport of CO$_2$ across the water barrier is driven by a combination of diffusion and convective dissolution. The experimental results are compared to and explained with a simple numerical model, with which the underlying mass transport processes are quantified. Analytical solutions that accurately predict the bubble growth dynamics are subsequently derived. The effect of convective dissolution across the water layer is treated as a reduction of the effective diffusion length, in accordance with the mass transfer scaling observed in laminar or natural convection. Finally, the binary water-bubble system is extended to a ternary water-bubble-alkane system. It is found that the alkane (n-hexadecane) layer bestows a buffering (hindering) effect on bubble growth and dissolution. The resulting growth dynamics and underlying fluxes are characterised theoretically.


title: Operating Tor Relays at Universities: Experiences and Considerations

abstract: In today's digital society, the Tor network has become an indispensable tool for individuals to protect their privacy on the Internet. Operated by volunteers, relay servers constitute the core component of Tor and are used to geographically escape surveillance. It is therefore essential to have a large, yet diverse set of relays. In this work, we analyze the contribution of educational institutions to the Tor network and report on our experience of operating exit relays at a university. Taking Germany as an example (but arguing that the global situation is similar), we carry out a quantitative study and find that universities contribute negligible amounts of relays and bandwidth. Since many universities all over the world have excellent conditions that render them perfect places to host Tor (exit) relays, we encourage other interested people and institutions to join. To this end, we discuss and resolve common concerns and provide lessons learned.


title: Factorizations of almost simple unitary groups

abstract: This is the second one in a series of papers classifying the factorizations of almost simple groups with nonsolvable factors. In this paper we deal with almost simple unitary groups.


title: Staircase Attention for Recurrent Processing of Sequences

abstract: Attention mechanisms have become a standard tool for sequence modeling tasks, in particular by stacking self-attention layers over the entire input sequence as in the Transformer architecture. In this work we introduce a novel attention procedure called staircase attention that, unlike self-attention, operates across the sequence (in time) recurrently processing the input by adding another step of processing. A step in the staircase comprises of backward tokens (encoding the sequence so far seen) and forward tokens (ingesting a new part of the sequence), or an extreme Ladder version with a forward step of zero that simply repeats the Transformer on each step of the ladder, sharing the weights. We thus describe a family of such models that can trade off performance and compute, by either increasing the amount of recurrence through time, the amount of sequential processing via recurrence in depth, or both. Staircase attention is shown to be able to solve tasks that involve tracking that conventional Transformers cannot, due to this recurrence. Further, it is shown to provide improved modeling power for the same size model (number of parameters) compared to self-attentive Transformers on large language modeling and dialogue tasks, yielding significant perplexity gains.


title: Optimizing a Binary Intelligent Reflecting Surface for OFDM Communications under Mutual Coupling

abstract: An intelligent reflecting surface (IRS) can greatly improve the channel quality over a frequency-flat channel, if it is configured to reflect the incident signal as a beam towards the receiver. However, the fundamental limitations of the IRS technology become apparent over practical frequency-selective channels, where the same configuration must be used over the entire bandwidth. In this paper, we consider a wideband orthogonal frequency-division multiplexing (OFDM) system that is supported by a fairly realistic IRS setup with two unbalanced states per element and also mutual coupling. We describe the simulation setup considered in the IEEE Signal Processing Cup 2021, propose a low-complexity solution for channel estimation and IRS configuration, and evaluate it on that setup.


title: Generative adversarial network with object detector discriminator for enhanced defect detection on ultrasonic B-scans

abstract: Non-destructive testing is a set of techniques for defect detection in materials. While the set of imaging techniques are manifold, ultrasonic imaging is the one used the most. The analysis is mainly performed by human inspectors manually analyzing recorded images. The low number of defects in real ultrasonic inspections and legal issues considering data from such inspections make it difficult to obtain proper results from automatic ultrasonic image (B-scan) analysis. In this paper, we present a novel deep learning Generative Adversarial Network model for generating ultrasonic B-scans with defects in distinct locations. Furthermore, we show that generated B-scans can be used for synthetic data augmentation, and can improve the performance of deep convolutional neural object detection networks. Our novel method is demonstrated on a dataset of almost 4000 B-scans with more than 6000 annotated defects. Defect detection performance when training on real data yielded average precision of 71%. By training only on generated data the results increased to 72.1%, and by mixing generated and real data we achieve 75.7% average precision. We believe that synthetic data generation can generalize to other challenges with limited datasets and could be used for training human personnel.


title: A Structured Optimal Controller with Feed-Forward for Transportation

abstract: We study an optimal control problem for a simple transportation model on a path graph. We give a closed form solution for the optimal controller, which can also account for planned disturbances using feed-forward. The optimal controller is highly structured, which allows the controller to be implemented using only local communication, conducted through two sweeps through the graph.


title: NWT: Towards natural audio-to-video generation with representation learning

abstract: In this work we introduce NWT, an expressive speech-to-video model. Unlike approaches that use domain-specific intermediate representations such as pose keypoints, NWT learns its own latent representations, with minimal assumptions about the audio and video content. To this end, we propose a novel discrete variational autoencoder with adversarial loss, dVAE-Adv, which learns a new discrete latent representation we call Memcodes. Memcodes are straightforward to implement, require no additional loss terms, are stable to train compared with other approaches, and show evidence of interpretability. To predict on the Memcode space, we use an autoregressive encoder-decoder model conditioned on audio. Additionally, our model can control latent attributes in the generated video that are not annotated in the data. We train NWT on clips from HBO's Last Week Tonight with John Oliver. NWT consistently scores above other approaches in Mean Opinion Score (MOS) on tests of overall video naturalness, facial naturalness and expressiveness, and lipsync quality. This work sets a strong baseline for generalized audio-to-video synthesis. Samples are available at https://next-week-tonight.github.io/NWT/.


title: LLAMA: The Low-Level Abstraction For Memory Access

abstract: The performance gap between CPU and memory widens continuously. Choosing the best memory layout for each hardware architecture is increasingly important as more and more programs become memory bound. For portable codes that run across heterogeneous hardware architectures, the choice of the memory layout for data structures is ideally decoupled from the rest of a program. This can be accomplished via a zero-runtime-overhead abstraction layer, underneath which memory layouts can be freely exchanged. We present the Low-Level Abstraction of Memory Access (LLAMA), a C++ library that provides such a data structure abstraction layer with example implementations for multidimensional arrays of nested, structured data. LLAMA provides fully C++ compliant methods for defining and switching custom memory layouts for user-defined data types. The library is extensible with third-party allocators. Providing two close-to-life examples, we show that the LLAMA-generated AoS (Array of Structs) and SoA (Struct of Arrays) layouts produce identical code with the same performance characteristics as manually written data structures. Integrations into the SPEC CPU\textsuperscript{\textregistered} lbm benchmark and the particle-in-cell simulation PIConGPU demonstrate LLAMA's abilities in real-world applications. LLAMA's layout-aware copy routines can significantly speed up transfer and reshuffling of data between layouts compared with naive element-wise copying. LLAMA provides a novel tool for the development of high-performance C++ applications in a heterogeneous environment.


title: Sensitivity analysis for random measurement error using regression calibration and simulation-extrapolation

abstract: Sensitivity analysis for measurement error can be applied in the absence of validation data by means of regression calibration and simulation-extrapolation. These have not been compared for this purpose. A simulation study was conducted comparing the performance of regression calibration and simulation-extrapolation in a multivariable model. The performance of the two methods was evaluated in terms of bias, mean squared error (MSE) and confidence interval coverage, for ranging reliability of the error-prone measurement (0.2-0.9), sample size (125-1,000), number of replicates (2-10), and R-squared (0.03-0.75). It was assumed that no validation data were available about the error-free measures, while measurement error variance was correctly estimated. In various scenarios, regression calibration was unbiased while simulation-extrapolation was biased: median bias was 1.4% (interquartile range (IQR): 0.8;2%), and -12.8% (IQR: -13.2;-11.0%), respectively. A small gain in efficiency was observed for simulation-extrapolation (median MSE: 0.005, IQR: 0.004;0.006) versus regression calibration (median MSE: 0.006, IQR: 0.004;0.007). Confidence interval coverage was at the nominal level of 95% for regression calibration, and smaller than 95% for simulation-extrapolation (median coverage: 92%, IQR: 85;94%). In the absence of validation data, the use of regression calibration is recommended for sensitivity analysis for measurement error.


title: Measurement of branching fractions and search for $CP$ violation in $D^{0}\to\pi^{+}\pi^{-}\eta$, $D^{0}\to K^{+}K^{-}\eta$, and $D^{0}\to\phi\eta$ at Belle

abstract: We measure the branching fractions and $CP$ asymmetries for the singly Cabibbo-suppressed decays $D^{0}\to\pi^{+}\pi^{-}\eta$, $D^{0}\to K^{+}K^{-}\eta$, and $D^{0}\to\phi\eta$, using 980 fb$^{-1}$ of data from the Belle experiment at the KEKB $e^+e^-$ collider. We obtain \begin{eqnarray} \mathcal{B}(D^{0}\to\pi^{+}\pi^{-}\eta) & = & [1.22\pm 0.02\,({\rm stat})\pm 0.02\,({\rm syst})\pm 0.03\,(\mathcal{B}_{\rm ref})]\times 10^{-3}\,, \nonumber \\ \mathcal{B}(D^{0}\to K^{+}K^{-}\eta) & = & [1.80\,^{+0.07}_{-0.06}\,({\rm stat})\pm 0.04\,({\rm syst})\pm 0.05\,(\mathcal{B}_{\rm ref})]\times 10^{-4}\,, \nonumber \\ \mathcal{B}(D^{0}\to\phi\eta) & = & [1.84\pm 0.09\,({\rm stat})\pm 0.06\,({\rm syst})\pm 0.05\,(\mathcal{B}_{\rm ref})]\times 10^{-4}\,, \nonumber \end{eqnarray} where the third uncertainty ($\mathcal{B}_{\rm ref}$) is from the uncertainty in the branching fraction of the reference mode $D^{0}\to K^{-}\pi^{+}\eta$. The color-suppressed decay $D^{0}\to\phi\eta$ is observed for the first time, with very high significance. The results for the $CP$ asymmetries are \begin{eqnarray} A_{CP}(D^{0}\to\pi^{+}\pi^{-}\eta) & = & [0.9\pm 1.2\,({\rm stat})\pm 0.5\,({\rm syst})]\%\,, \nonumber \\ A_{CP}(D^{0}\to K^{+}K^{-}\eta) & = & [-1.4\pm 3.3\,({\rm stat})\pm 1.1\,({\rm syst})]\%\,, \nonumber \\ A_{CP}(D^{0}\to\phi\eta)&= & [-1.9\pm 4.4\,({\rm stat})\pm 0.6\,({\rm syst})]\%\,. \nonumber \end{eqnarray} The results for $D^{0}\to\pi^{+}\pi^{-}\eta$ are a significant improvement over previous results. The branching fraction and $A_{CP}$ results for $D^{0}\to K^{+}K^{-}\eta$, and the $A_{CP}$ result for $D^{0}\to\phi\eta$, are the first such measurements. No evidence for $CP$ violation is found in any of these decays.


title: Spectrally-wide acoustic frequency combs generated using oscillations of polydisperse gas bubble clusters in liquids

abstract: Acoustic frequency combs leverage unique properties of the optical frequency comb technology in high-precision measurements and innovative sensing in optically inaccessible environments such as under water, under ground or inside living organisms. Because acoustic combs with wide spectra would be required for many of these applications but techniques of their generation have not yet been developed, here we propose a new approach to the creation of spectrally-wide acoustic combs using oscillations of polydisperse gas bubble clusters in liquids. By means of numerical simulations we demonstrate that clusters consisting of bubbles with precisely controlled sizes can produce wide acoustic spectra composed of equally-spaced coherent peaks. We show that under typical experimental conditions bubble clusters remain stable over time required for a reliable recording of comb signals. We also demonstrate that the spectral composition of combs can be tuned by adjusting the number and size of bubbles in a cluster.


title: Infinitely many solutions for Schr\"{o}dinger-Newton equations

abstract: We prove the existence of infinitely many non-radial positive solutions for the Schr\"{o}dinger-Newton system $$ \left\{\begin{array}{ll} \Delta u- V(|x|)u + \Psi u=0, &x\in\mathbb{R}^3,\newline \Delta \Psi+\frac12 u^2=0, &x\in\mathbb{R}^3, \end{array}\right. $$ provided that $V(r)$ has the following behavior at infinity: $$ V(r)=V_0+\frac{a}{r^m}+O\left(\frac{1}{r^{m+\theta}}\right) \quad\mbox{ as } r\rightarrow\infty, $$ where $\frac12\le m<1$ and $a, V_0, \theta$ are some positive constants. In particular, for any $s$ large we use a reduction method to construct $s-$bump solutions lying on a circle of radius $r\sim (s\log s)^{\frac{1}{1-m}}$.


title: Morphing tree drawings in a small 3D grid

abstract: We study crossing-free grid morphs for planar tree drawings using 3D. A morph consists of morphing steps, where vertices move simultaneously along straight-line trajectories at constant speeds. A crossing-free morph is known between two drawings of an $n$-vertex planar graph $G$ with $\mathcal{O}(n)$ morphing steps and using the third dimension it can be reduced to $\mathcal{O}(\log n)$ for an $n$-vertex tree [Arseneva et al.\ 2019]. However, these morphs do not bound one practical parameter, the resolution. Can the number of steps be reduced substantially by using the third dimension while keeping the resolution bounded throughout the morph? We answer this question in an affirmative and present a 3D non-crossing morph between two planar grid drawings of an $n$-vertex tree in $\mathcal{O}(\sqrt{n} \log n)$ morphing steps. Each intermediate drawing lies in a $3D$ grid of polynomial volume.


title: No symplectic-Lipschitz structures on $S^{2n \geq 4}$

abstract: We prove that a closed manifold which admits a symplectic-Lipschitz structure has non-vanishing even-degree cohomology groups with real coefficients. In particular, spheres $S^{2n}, n \geq 2$ do not admit symplectic-Lipschitz structures.


title: Explicit descent on elliptic curves and splitting Brauer classes

abstract: We prove new results on splitting Brauer classes by genus 1 curves, settling in particular the case of degree 7 classes over global fields. Though our method is cohomological in nature, and proceeds by considering the more difficult problem of splitting $\mu_N$-gerbes, we use crucial input from the arithmetic of modular curves and explicit $N$-descent on elliptic curves.


title: Principled Hyperedge Prediction with Structural Spectral Features and Neural Networks

abstract: Hypergraph offers a framework to depict the multilateral relationships in real-world complex data. Predicting higher-order relationships, i.e hyperedge, becomes a fundamental problem for the full understanding of complicated interactions. The development of graph neural network (GNN) has greatly advanced the analysis of ordinary graphs with pair-wise relations. However, these methods could not be easily extended to the case of hypergraph. In this paper, we generalize the challenges of GNN in representing higher-order data in principle, which are edge- and node-level ambiguities. To overcome the challenges, we present SNALS that utilizes bipartite graph neural network with structural features to collectively tackle the two ambiguity issues. SNALS captures the joint interactions of a hyperedge by its local environment, which is retrieved by collecting the spectrum information of their connections. As a result, SNALS achieves nearly 30% performance increase compared with most recent GNN-based models. In addition, we applied SNALS to predict genetic higher-order interactions on 3D genome organization data. SNALS showed consistently high prediction accuracy across different chromosomes, and generated novel findings on 4-way gene interaction, which is further validated by existing literature.


title: Modeling Uplink Coverage Performance in Hybrid Satellite-Terrestrial Networks

abstract: Once deemed a far-fetched vision, emerging deployments of massive satellite constellations will soon offer true global coverage. When these constellations overlay the evolving terrestrial networks, a new hybrid continuum is formed which has the potential to provide uninterrupted coverage. In this paper, we provide an analytic framework for the uplink coverage probability in hybrid satellite-terrestrial networks. The framework extends well-developed terrestrial models that utilize tools from stochastic geometry by incorporating an additional layer that fits the emerging next generation satellite-terrestrial networks. The paper captures the impact of both (i) the constellation size, and (ii) the terrestrial base station density on the coverage of the uplink traffic which is dominant in applications relying on wireless sensor networks, such as the Internet of Things. This framework provides insights that guide the design of hybrid network infrastructures for a desired quality of service.


title: Physics graduate student employment: what we can learn from professional social media

abstract: In this study, we investigated the employment status of recent University of Ottawa physics MSc and PhD graduates, finding that 94% of graduates are either employed or pursuing further physics education one year post-graduation. Our database was populated from the public online repository of MSc and PhD theses submitted between the academic years of 2011 to 2019, with employment information collected in 2020 from the professional social media platform LinkedIn. Our results highlight that graduates primarily find employment quickly and in their field of study, with most graduates employed in either academia or physics-related industries. We also found that a significant portion of employed graduates, 20%, find employment in non-traditional physics careers, such as business management and healthcare. Graduates with careers in academia tend to have lower online connectivity compared to graduates with careers in industry or non-traditional fields, suggesting a greater importance for online networking for students interested in non-academic careers.


title: Nodal-line semimetal HMTSF-TCNQ: Anomalous orbital diamagnetism and charge density wave

abstract: This study investigates the electronic states and physical quantities of an organic charge-transfer complex HMTSF-TCNQ, which undergoes a charge-density-wave (CDW) phase transition at temperature $T_c\simeq 30$ K. A first-principles calculation is utilized to determine that the normal state is a topological semimetal with open nodal lines. Besed on the first-principles calculation, we develop a tight-binding model to investigate the electronic state in detail. Below $T_c$, the CDW phase is examined in the tight-binding scheme using the mean-field approximation. It is shown that the open nodal lines are deformed into closed ones, and their shapes are sensitive to the order parameter. Using this tight-binding model, we theoretically evaluate the temperature dependencies of two physical quantities: the spin-lattice relaxation time $T_1$ and the orbital magnetic susceptibility. In particular, an anomalous plateau is obtained at low temperatures in the orbital diamagnetism. We presume that this anomalous plateau originates owing to the conflict between the interband diamagnetism, impurity scattering, and the nodal line deformation. We also conduct an experiment to investigate the orbital magnetism, and the results are in excellent quantitative agreement with the theory.


title: On a problem with conjugation conditions for an equation of even order with a fractional derivative in the sense of Caputo

abstract: In this work, for an equation of high even order with a fractional derivative in the sense of Caputo, a problem is studied in a rectangular domain with conjugation conditions. A criterion for the uniqueness of a solution is given. The solution is constructed in the form of a Fourier series in the eigenfunctions of a one-dimensional problem.


title: Hamming graphs and concentration properties in non-quasi-reflexive Banach spaces

abstract: In this note, we study some concentration properties for Lipschitz maps defined on Hamming graphs, as well as their stability under sums of Banach spaces. As an application, we extend a result of Causey on the coarse Lipschitz structure of quasi-reflexive spaces satisfying upper $\ell_p$ tree estimates to the setting of $\ell_p$-sums of such spaces. Our result provides us with a tool for constructing the first examples of Banach spaces that are not quasi-reflexive but nevertheless admit some concentration inequality. We also give a sufficient condition for a space to be asymptotic-$c_0$ in terms of a concentration property, as well as relevant counterexamples.


title: Unsupervised Word Segmentation from Discrete Speech Units in Low-Resource Settings

abstract: Documenting languages helps to prevent the extinction of endangered dialects, many of which are otherwise expected to disappear by the end of the century. When documenting oral languages, unsupervised word segmentation (UWS) from speech is a useful, yet challenging, task. It consists in producing time-stamps for slicing utterances into smaller segments corresponding to words, being performed from phonetic transcriptions, or in the absence of these, from the output of unsupervised speech discretization models. These discretization models are trained using raw speech only, producing discrete speech units that can be applied for downstream (text-based) tasks. In this paper we compare five of these models: three Bayesian and two neural approaches, with regards to the exploitability of the produced units for UWS. For the UWS task, we experiment with two models, using as our target language the Mboshi (Bantu C25), an unwritten language from Congo-Brazzaville. Additionally, we report results for Finnish, Hungarian, Romanian and Russian in equally low-resource settings, using only 4 hours of speech. Our results suggest that neural models for speech discretization are difficult to exploit in our setting, and that it might be necessary to adapt them to limit sequence length. We obtain our best UWS results by using Bayesian models that produce high quality, yet compressed, discrete representations of the input speech signal.


title: A direct product theorem for quantum communication complexity with applications to device-independent cryptography

abstract: We give a direct product theorem for the entanglement-assisted interactive quantum communication complexity of an $l$-player predicate $\mathsf{V}$. In particular we show that for a distribution $p$ that is product across the input sets of the $l$ players, the success probability of any entanglement-assisted quantum communication protocol for computing $n$ copies of $\mathsf{V}$, whose communication is $o(\log(\mathrm{eff}^*(\mathsf{V},p))\cdot n)$, goes down exponentially in $n$. Here $\mathrm{eff}^*(\mathsf{V}, p)$ is a distributional version of the quantum efficiency or partition bound introduced by Laplante, Lerays and Roland (2014), which is a lower bound on the distributional quantum communication complexity of computing a single copy of $\mathsf{V}$ with respect to $p$. Applying our direct product theorem for small communication and techniques related to $\mathrm{eff}^*$, we show that it is possible to do device-independent (DI) quantum cryptography without the assumption that devices do not leak any information. First, we analyze the parallel DI quantum key distribution protocol given by Jain, Miller and Shi (2020), and show that when the protocol is carried out with devices that are compatible with $n$ copies of the Magic Square game, it is possible to extract $\Omega(n)$ bits of key from it, even in the presence of $O(n)$ bits of leakage. Second, we show that it is possible to do sequential versions of the Jain, Miller and Shi protocol, which give a better key rate for QKD with leakage, and let us do sequential DI randomness expansion with leakage (it is not known how to do parallel DI randomness expansion even without leakage). Third, we show that proofs of quantumness with two entangled provers are resistant to leakage, i.e., classical players who communicate $O(n)$ bits with each other cannot convince the verifier that they share entanglement.


title: A Unified Generative Framework for Aspect-Based Sentiment Analysis

abstract: Aspect-based Sentiment Analysis (ABSA) aims to identify the aspect terms, their corresponding sentiment polarities, and the opinion terms. There exist seven subtasks in ABSA. Most studies only focus on the subsets of these subtasks, which leads to various complicated ABSA models while hard to solve these subtasks in a unified framework. In this paper, we redefine every subtask target as a sequence mixed by pointer indexes and sentiment class indexes, which converts all ABSA subtasks into a unified generative formulation. Based on the unified formulation, we exploit the pre-training sequence-to-sequence model BART to solve all ABSA subtasks in an end-to-end framework. Extensive experiments on four ABSA datasets for seven subtasks demonstrate that our framework achieves substantial performance gain and provides a real unified end-to-end solution for the whole ABSA subtasks, which could benefit multiple tasks.


title: Selection rules in symmetry-broken systems by symmetries in synthetic dimensions

abstract: Selection rules are often considered a hallmark of symmetry. When a symmetry is broken, e.g., by an external perturbation, the system exhibits selection rule deviations which are often analyzed by perturbation theory. Here, we employ symmetry-breaking degrees of freedom as synthetic dimensions, to demonstrate that symmetry-broken systems systematically exhibit a new class of symmetries and selection rules. These selection rules determine the scaling of a system's observables (to all orders in the strength of the symmetry-breaking perturbation) as it transitions from symmetric to symmetry-broken. We specifically analyze periodically driven (Floquet) systems subject to two driving fields, where the first field imposes a spatio-temporal symmetry, and the second field breaks it, imposing a symmetry in synthetic dimensions. We tabulate the resulting synthetic symmetries for (2+1)D Floquet group symmetries and derive the corresponding selection rules for high harmonic generation (HHG) and above-threshold ionization (ATI). Finally, we observe experimentally HHG selection rules imposed by symmetries in synthetic dimensions. The new class of symmetries & selection rules extends the scope of existing symmetry breaking spectroscopy techniques, opening new routes for ultrafast spectroscopy of phonon-polarization, spin-orbit coupling, and more.


title: Obtaining Better Static Word Embeddings Using Contextual Embedding Models

abstract: The advent of contextual word embeddings -- representations of words which incorporate semantic and syntactic information from their context -- has led to tremendous improvements on a wide variety of NLP tasks. However, recent contextual models have prohibitively high computational cost in many use-cases and are often hard to interpret. In this work, we demonstrate that our proposed distillation method, which is a simple extension of CBOW-based training, allows to significantly improve computational efficiency of NLP applications, while outperforming the quality of existing static embeddings trained from scratch as well as those distilled from previously proposed methods. As a side-effect, our approach also allows a fair comparison of both contextual and static embeddings via standard lexical evaluation tasks.


title: On Primordial Black Holes and secondary gravitational waves generated from inflation with solo/multi-bumpy potential

abstract: It is well known that a primordial black hole (PBH) can be generated in the inflation process of the early universe, especially when the inflation field has a number of non-trivial features that could break the slow-roll condition. In this study, we investigate a toy model of inflation with bumpy potential, which has one or several bumps. We determined that the potential with multi-bump can generate power spectra with multi-peaks in small-scale region, which can in turn predict the generation of primordial black holes in various mass ranges. We also consider the two possibilities of PBH formation by spherical and elliptical collapses. Finally, we discuss the scalar-induced gravitational waves(SIGWs)generated by linear scalar perturbations at second order.


title: Methodological considerations for estimating policy effects in the context of co-occurring policies

abstract: Understanding how best to estimate state-level policy effects is important, and several unanswered questions remain, particularly about the ability of statistical models to disentangle the effects of concurrently enacted policies. In practice, many policy evaluation studies do not attempt to control for effects of co-occurring policies, and this issue has not received extensive attention in the methodological literature to date. In this study, we utilized Monte Carlo simulations to assess the impact of co-occurring policies on the performance of commonly-used statistical models in state policy evaluations. Simulation conditions varied effect sizes of the co-occurring policies and length of time between policy enactment dates, among other factors. Outcome data (annual state-specific opioid mortality rate per 100,000) were obtained from 1999-2016 National Vital Statistics System (NVSS) Multiple Cause of Death mortality files, thus yielding longitudinal annual state-level data over 18 years from 50 states. When co-occurring policies are ignored (i.e., omitted from the analytic model), our results demonstrated that high relative bias (>85%) arises, particularly when policies are enacted in rapid succession. Moreover, as expected, controlling for all co-occurring policies will effectively mitigate the threat of confounding bias; however, effect estimates may be relatively imprecise (i.e., larger variance) when policies are enacted in near succession. Our findings highlight several key methodological issues regarding co-occurring policies in the context of opioid-policy research yet also generalize more broadly to evaluation of other state-level policies, such as policies related to firearms or COVID-19, showcasing the need to think critically about co-occurring policies that are likely to influence the outcome when specifying analytic models.


title: A hybrid classical-quantum approach to solve the heat equation using quantum annealers

abstract: The numerical solution of partial differential equations by discretization techniques is ubiquitous in computational physics. In this work we benchmark this approach in the quantum realm by solving the heat equation for a square plate subject to fixed temperatures at the edges and random heat sources and sinks within the domain. The hybrid classical-quantum approach consists in the solution on a quantum computer of the coupled linear system of equations that result from the discretization step. Owing to the limitations in the number of qubits and their connectivity, we use the Gauss-Seidel method to divide the full system of linear equations into subsystems, which are solved iteratively in block fashion. Each of the linear subsystems were solved using 2000Q and Advantage quantum computers developed by D-Wave Systems Inc. By comparing classical numerical and quantum solutions, we observe that the errors and chain break fraction are, on average, greater on the 2000Q system. Unlike the classical Gauss-Seidel method, the errors of the quantum solutions level off after a few iterations of our algorithm. This is partly a result of the span of the real number line available from the mapping of the chosen size of the set of qubit states. We verified this by using techniques to progressively shrink the range mapped by the set of qubit states at each iteration (increasing floating-point accuracy). As a result, no leveling off is observed. However, an increase in qubits does not translate to an overall lower error. This is believed to be indicative of the increasing length of chains required for the mapping to real numbers and the ensuing limitations of hardware.


title: Residual Feedback Learning for Contact-Rich Manipulation Tasks with Uncertainty

abstract: While classic control theory offers state of the art solutions in many problem scenarios, it is often desired to improve beyond the structure of such solutions and surpass their limitations. To this end, residual policy learning (RPL) offers a formulation to improve existing controllers with reinforcement learning (RL) by learning an additive "residual" to the output of a given controller. However, the applicability of such an approach highly depends on the structure of the controller. Often, internal feedback signals of the controller limit an RL algorithm to adequately change the policy and, hence, learn the task. We propose a new formulation that addresses these limitations by also modifying the feedback signals to the controller with an RL policy and show superior performance of our approach on a contact-rich peg-insertion task under position and orientation uncertainty. In addition, we use a recent Cartesian impedance control architecture as the control framework which can be available to us as a black-box while assuming no knowledge about its input/output structure, and show the difficulties of standard RPL. Furthermore, we introduce an adaptive curriculum for the given task to gradually increase the task difficulty in terms of position and orientation uncertainty. A video showing the results can be found at https://youtu.be/SAZm_Krze7U .


title: Infinite-color randomly reinforced urns with dominant colors

abstract: We define and prove limit results for a class of dominant P\'olya sequences, which are randomly reinforced urn processes with color-specific random weights and unbounded number of possible colors. Under fairly mild assumptions on the expected reinforcement, we show that the predictive and the empirical distributions converge almost surely (a.s.) in total variation to the same random probability measure $\tilde{P}$; moreover, $\tilde{P}(\mathcal{D})=1$ a.s., where $\mathcal{D}$ denotes the set of dominant colors for which the expected reinforcement is maximum. In the general case, the predictive probabilities and the empirical frequencies of any $\delta$-neighborhood of $\mathcal{D}$ converge a.s. to one. That is, although non-dominant colors continue to be regularly observed, their distance to $\mathcal{D}$ converges in probability to zero. We refine the above results with rates of convergence. We further hint potential applications of dominant P\'olya sequences in randomized clinical trials and species sampling, and use our central limit results for Bayesian inference.


title: Relativistic Fock space coupled cluster study of bismuth electronic structure to extract the Bi nuclear quadrupole moment

abstract: We report the value of the electric quadrupole moment of $^{209}$Bi extracted from the atomic data. For this, we performed electronic structure calculations for the ground $^4S^o_{3/2}$ and excited $^2P^o_{3/2}$ states of atomic Bi using the Dirac-Coulomb-Breit Hamiltonian and the Fock space coupled cluster method with single, double, and full triple amplitudes for the three-particle Fock space sector. The value of the quadrupole moment of $^{209}$Bi, $Q(^{209}$Bi$)=-418(6)$~mb, derived from the resulting electric field gradient values and available atomic hyperfine splittings is in excellent agreement with molecular data. Due to the availability of the hyperfine constants for unstable isotopes of Bi, current atomic calculation allows also to correct their quadrupole moments.


title: On the $16$-rank of class groups of $\mathbb{Q}(\sqrt{-3p})$ for primes $p$ congruent to $1$ modulo $4$

abstract: For fixed $q\in\{3,7,11,19, 43,67,163\}$, we consider the density of primes $p$ congruent to $1$ modulo $4$ such that the class group of the number field $\mathbb{Q}(\sqrt{-qp})$ has order divisible by $16$. We show that this density is equal to $1/8$, in line with a more general conjecture of Gerth. Vinogradov's method is the key analytic tool for our work.


title: Principled Data Completion of Network Constraints for Day Ahead Auctions in Power Markets

abstract: Network constraints play a key role in the price finding mechanism for European Power Markets, but historical data is very sparse and usually insufficient for many quantitative applications. We reconstruct the constraints data, known as the Power Transmission Distribution Factors (PTDFs) and Remaining Available Margins (RAMs), by first recovering the underlying time dependent signals known as the Generation Shift Keys (GSKs) and Phase Angles (PAs), and the electricity grid characteristics, via a mathematical optimisation problem. This is solved by exploiting marginal convexity in certain subspaces via alternating minimisation. The GSKs and PAs are then mapped to the PTDFs and RAMs, using the grid structure. Our reconstruction achieves good in-sample and out-of-sample relative errors for the PTDFs and RAMs. We further show that our model outperforms the naive approach, and that the reconstructed GSKs and PAs recover specific structure.


title: Hyperbolic Temporal Knowledge Graph Embeddings with Relational and Time Curvatures

abstract: Knowledge Graph (KG) completion has been excessively studied with a massive number of models proposed for the Link Prediction (LP) task. The main limitation of such models is their insensitivity to time. Indeed, the temporal aspect of stored facts is often ignored. To this end, more and more works consider time as a parameter to complete KGs. In this paper, we first demonstrate that, by simply increasing the number of negative samples, the recent AttH model can achieve competitive or even better performance than the state-of-the-art on Temporal KGs (TKGs), albeit its nontemporality. We further propose Hercules, a time-aware extension of AttH model, which defines the curvature of a Riemannian manifold as the product of both relation and time. Our experiments show that both Hercules and AttH achieve competitive or new state-of-the-art performances on ICEWS04 and ICEWS05-15 datasets. Therefore, one should raise awareness when learning TKGs representations to identify whether time truly boosts performances.


title: Speech BERT Embedding For Improving Prosody in Neural TTS

abstract: This paper presents a speech BERT model to extract embedded prosody information in speech segments for improving the prosody of synthesized speech in neural text-to-speech (TTS). As a pre-trained model, it can learn prosody attributes from a large amount of speech data, which can utilize more data than the original training data used by the target TTS. The embedding is extracted from the previous segment of a fixed length in the proposed BERT. The extracted embedding is then used together with the mel-spectrogram to predict the following segment in the TTS decoder. Experimental results obtained by the Transformer TTS show that the proposed BERT can extract fine-grained, segment-level prosody, which is complementary to utterance-level prosody to improve the final prosody of the TTS speech. The objective distortions measured on a single speaker TTS are reduced between the generated speech and original recordings. Subjective listening tests also show that the proposed approach is favorably preferred over the TTS without the BERT prosody embedding module, for both in-domain and out-of-domain applications. For Microsoft professional, single/multiple speakers and the LJ Speaker in the public database, subjective preference is similarly confirmed with the new BERT prosody embedding. TTS demo audio samples are in https://judy44chen.github.io/TTSSpeechBERT/.


title: On the approximation exponents for subspaces of $\mathbb{R}^n$

abstract: This paper follows the generalisation of the classical theory of Diophantine approximation to subspaces of $\mathbb{R}^n$ established by W. M. Schmidt in 1967. Let $A$ and $B$ be two subspaces of $\mathbb{R}^n$ of respective dimensions $d$ and $e$ with $d+e\leqslant n$. The proximity between $A$ and $B$ is measured by $t=\min(d,e)$ canonical angles $0\leqslant \theta_1\leqslant \cdots\leqslant \theta_t\leqslant \pi/2$; we set $\psi_j(A,B)=\sin\theta_j$. If $B$ is a rational subspace, his complexity is measured by its height $H(B)=\mathrm{covol}(B\cap\mathbb{Z}^n)$. We denote by $\mu_n(A\vert e)_j$ the exponent of approximation defined as the upper bound (possibly equal to $+\infty$) of the set of $\beta>0$ such that the inequality $\psi_j(A,B)\leqslant H(B)^{-\beta}$ holds for infinitely many rational subspaces $B$ of dimension $e$. We are interested in the minimal value $\mathring{\mu}_n(d\vert e)_j$ taken by $\mu_n(A\vert e)_j$ when $A$ ranges through the set of subspaces of dimension $d$ of $\mathbb{R}^n$ such that for all rational subspaces $B$ of dimension $e$ one has $\dim (A\cap B)<j$. We show that $\mathring{\mu}_4(2\vert 2)_1=3$, $\mathring{\mu}_5(3\vert 2)_1\le 6$ and $\mathring{\mu}_{2d}(d\vert \ell)_1\leqslant 2d^2/(2d-\ell)$. We also prove a lower bound in the general case, which implies that $\mathring{\mu}_n(d\vert d)_d\xrightarrow[n\to+\infty]{} 1/d$.


title: A Perspective on Time towards Wireless 6G

abstract: With the advent of 5G technology, the notion of latency got a prominent role in wireless connectivity, serving as a proxy term for addressing the requirements for real-time communication. As wireless systems evolve towards 6G, the ambition to immerse the digital into the physical reality will increase. Besides making the real-time requirements more stringent, this immersion will bring the notions of time, simultaneity, presence, and causality to a new level of complexity. A growing body of research points out that latency is insufficient to parameterize all real-time requirements. Notably, one such requirement that received a significant attention is information freshness, defined through the Age of Information (AoI) and its derivatives. The objective of this article is to investigate the general notion of timing in wireless communication systems and networks and its relation to effective information generation, processing, transmission, and reconstruction at the senders and receivers. We establish a general statistical framework of timing requirements in wireless communication systems, which subsumes both latency and AoI. The framework is made by associating a timing component with the two basic statistical operations, decision and estimation. We first use the framework to present a representative sample of the existing works that deal with timing in wireless communication. Next, it is shown how the framework can be used with different communication models of increasing complexity, starting from the basic Shannon one-way communication model and arriving to communication models for consensus, distributed learning, and inference. Overall, this paper fills an important gap in the literature by providing a systematic treatment of various timing measures in wireless communication and sets the basis for design and optimization for the next-generation real-time systems.


title: Learning Riemannian Manifolds for Geodesic Motion Skills

abstract: For robots to work alongside humans and perform in unstructured environments, they must learn new motion skills and adapt them to unseen situations on the fly. This demands learning models that capture relevant motion patterns, while offering enough flexibility to adapt the encoded skills to new requirements, such as dynamic obstacle avoidance. We introduce a Riemannian manifold perspective on this problem, and propose to learn a Riemannian manifold from human demonstrations on which geodesics are natural motion skills. We realize this with a variational autoencoder (VAE) over the space of position and orientations of the robot end-effector. Geodesic motion skills let a robot plan movements from and to arbitrary points on the data manifold. They also provide a straightforward method to avoid obstacles by redefining the ambient metric in an online fashion. Moreover, geodesics naturally exploit the manifold resulting from multiple--mode tasks to design motions that were not explicitly demonstrated previously. We test our learning framework using a 7-DoF robotic manipulator, where the robot satisfactorily learns and reproduces realistic skills featuring elaborated motion patterns, avoids previously unseen obstacles, and generates novel movements in multiple-mode settings.


title: Exploration and preference satisfaction trade-off in reward-free learning

abstract: Biological agents have meaningful interactions with their environment despite the absence of immediate reward signals. In such instances, the agent can learn preferred modes of behaviour that lead to predictable states -- necessary for survival. In this paper, we pursue the notion that this learnt behaviour can be a consequence of reward-free preference learning that ensures an appropriate trade-off between exploration and preference satisfaction. For this, we introduce a model-based Bayesian agent equipped with a preference learning mechanism (pepper) using conjugate priors. These conjugate priors are used to augment the expected free energy planner for learning preferences over states (or outcomes) across time. Importantly, our approach enables the agent to learn preferences that encourage adaptive behaviour at test time. We illustrate this in the OpenAI Gym FrozenLake and the 3D mini-world environments -- with and without volatility. Given a constant environment, these agents learn confident (i.e., precise) preferences and act to satisfy them. Conversely, in a volatile setting, perpetual preference uncertainty maintains exploratory behaviour. Our experiments suggest that learnable (reward-free) preferences entail a trade-off between exploration and preference satisfaction. Pepper offers a straightforward framework suitable for designing adaptive agents when reward functions cannot be predefined as in real environments.


title: Equations of state of new boron-rich selenides B$_6$Se and B$_{12}$Se

abstract: Two novel of boron-rich selenides, orthorhombic B$_6$Se and rhombohedral B$_{12}$Se, have been recently synthesized at high pressure - high temperature conditions. Room-temperature compressibilities of these phases were studied in a diamond anvil cell using synchrotron powder X-ray diffraction. A fit of experimental p-V data by third-order Birch-Murnaghan equation of state yielded the bulk moduli of 155(2) GPa for B$_{12}$Se and 144(3) GPa for B$_6$Se. No pressure-induced phase transitions have been observed in the studied pressure range, i.e., up to 35 GPa.


title: Observation of mutual extinction and transparency in light scattering

abstract: Interference of scattered waves is fundamental for modern light-scattering techniques, such as optical wavefront shaping. Recently, a new type of wavefront shaping was introduced where the extinction is manipulated instead of the scattered intensity. The underlying idea is that upon changing the phases or the amplitudes of incident beams, the total extinction will change due to interference described by the cross terms between different incident beams. Here, we experimentally demonstrate the mutual extinction and transparency effects in scattering media, in particular, a human hair and a silicon bar. To this end, we send two light beams with a variable mutual angle on the sample. Depending on the relative phase of the incident beams we observe either nearly zero extinction, mutual transparency, or almost twice the single-beam extinction, mutual extinction, in agreement with theory. We use an analytical approximation for the scattering amplitude, starting from a completely opaque object and we discuss the limitations of our approximation. We discuss the applications of the mutual extinction and transparency effects in various fields such as non-line-of-sight communications, microscopy, and biomedical imaging.


title: Breaking the Limits of Message Passing Graph Neural Networks

abstract: Since the Message Passing (Graph) Neural Networks (MPNNs) have a linear complexity with respect to the number of nodes when applied to sparse graphs, they have been widely implemented and still raise a lot of interest even though their theoretical expressive power is limited to the first order Weisfeiler-Lehman test (1-WL). In this paper, we show that if the graph convolution supports are designed in spectral-domain by a non-linear custom function of eigenvalues and masked with an arbitrary large receptive field, the MPNN is theoretically more powerful than the 1-WL test and experimentally as powerful as a 3-WL existing models, while remaining spatially localized. Moreover, by designing custom filter functions, outputs can have various frequency components that allow the convolution process to learn different relationships between a given input graph signal and its associated properties. So far, the best 3-WL equivalent graph neural networks have a computational complexity in $\mathcal{O}(n^3)$ with memory usage in $\mathcal{O}(n^2)$, consider non-local update mechanism and do not provide the spectral richness of output profile. The proposed method overcomes all these aforementioned problems and reaches state-of-the-art results in many downstream tasks.


title: Solar Cycle Evolution of Filaments over a Century: Investigations with the Meudon and McIntosh Hand-drawn Archives

abstract: Hand-drawn synoptic maps from the Meudon Observatory (1919 onwards) and the McIntosh archive (1967 onwards) are two important sources of long-term, manually recorded filament observations. In this study, we calibrate the Meudon maps and subsequently identify filaments through an automated method. We extract physical parameters from this filament database and perform a comparative study of their long-term evolution focusing on the cotemporal period of McIntosh and Meudon observations. The spatio-temporal evolution of filaments manifests in the form of a filament butterfly diagram, indicating further that they are intimately related to the large-scale solar cycle. Physical descriptors such as the number and length of filaments, which are tracers of solar surface magnetic field, have cycles which are phase-locked with the 11 year sunspot cycle. The tilt angle distribution of filaments - both near or distant from active region locations - indicates that their origin is due to either large-scale surface magnetic field or inter-active region field evolution. This study paves the way for constructing a composite series of hand-drawn filament data with minimal gaps stretching the time span of solar filament observations to a century. On the one hand, this would serve as useful constraints for models of magnetic field emergence and evolution on the Sun's surface over multiple solar cycles, and on the other hand, this filament database may be used to guide the reconstruction of filament-prominence associated eruptive events before the space age


title: $\mathbf{{}^{12}{C} + {}^{12}{C}}$ Fusion $\boldsymbol{S^*}$-factor from a Full-microscopic Nuclear Model

abstract: The ${}^{12}\mathrm{C} + {}^{12}\mathrm{C}$ fusion reaction plays a vital role in the explosive phenomena of the universe. The resonances in the Gamow window rule its reaction rate and products. Hence, the determination of the resonance parameters by nuclear models is indispensable as the direct measurement is not feasible. Here, for the first time, we report the resonances in the ${}^{12}\mathrm{C} + {}^{12}\mathrm{C}$ fusion reaction described by a full-microscopic nuclear model. The model plausibly reproduces the measured low-energy astrophysical $S$-factors and predicts the resonances in the Gamow window. Contradictory to the hindrance model, we conclude that there is no low-energy suppression of the $S$-factor.


title: Rediscussion of eclipsing binaries. Paper V. The triple system V455 Aurigae

abstract: V455 Aur is a detached eclipsing binary containing two F-stars in a 3.15-d orbit with a small eccentricity. Its eclipses were discovered in data from the Hipparcos satellite and a spectroscopic orbit was obtained by Griffin (2001, 2013). Griffin found a long-term variation of the systemic velocity of the eclipsing system due to a third body in a highly eccentric orbit (e = 0.73) with a period of 4200 d. We have used these data, the light curve of V455 Aur from the TESS satellite, and the Gaia EDR3 parallax to determine the physical properties of the components of the system to high precision. We find the eclipsing stars to have masses of 1.289 +/- 0.006 Msun and 1.232 +/- 0.005 Msun, radii of 1.389 +/- 0.011 Rsun and 1.318 +/- 0.014 Rsun and effective temperatures of 6500 +/- 200 and 6400 +/- 200 K. Light from the tertiary component is directly detected for the first time, in the form of a third light of l_3 = 0.028 +/- 0.002 in the solution of the TESS light curve. From this l_3, theoretical spectra and empirical calibrations we estimate the star to have a mass of 0.72 +/- 0.05 Msun, a radius of 0.74 +/- 0.05 Rsun and a temperature of 4300 +/- 300 K. The inclination of the outer orbit is 53 +/- 3 degrees, so the two orbits in the system are not coplanar. We show that a measured spectroscopic light ratio of the two eclipsing stars could lower the uncertainties in radius from 1% to 0.25%. A detailed spectroscopic analysis could also yield precise temperatures and chemical abundances of the system, thus making V455 Aur one of the most precisely measured eclipsing systems known.


title: Rediscussion of eclipsing binaries. Paper VI. The F-type system V505 Persei

abstract: V505 Per is a detached eclipsing binary containing two F5 V stars in a 4.22-d circular orbit. We use a light curve from the TESS satellite and published radial velocity measurements to establish the properties of the system to high precision. The masses of the stars are 1.275 +/- 0.004 Msun and 1.258 +/- 0.003 Msun, and their radii are 1.294 +/- 0.002 Rsun and 1.264 +/- 0.002 Rsun. Adding published effective temperature estimates, we precisely measure the luminosities and absolute bolometric magnitudes of the stars, and the distance to the system. The distance is slightly shorter than that obtained from the Gaia EDR3 parallax, a discrepancy most easily explained by uncertainty in the 2MASS K-band apparent magnitude. We reanalyse existing light and radial velocity curves from three previous studies of this system and conclude that, in this case, formal errors are reliable for the spectroscopic orbits but not light curves, that errorbars from a residual-permutation algorithm are suitable for light curves but not spectroscopic orbits, and that published results are not always reproducible. The precisions in the measured properties of V505 Per are high and among the best ever obtained for a detached eclipsing binary system.


title: Contrastive Representation Learning for Hand Shape Estimation

abstract: This work presents improvements in monocular hand shape estimation by building on top of recent advances in unsupervised learning. We extend momentum contrastive learning and contribute a structured collection of hand images, well suited for visual representation learning, which we call HanCo. We find that the representation learned by established contrastive learning methods can be improved significantly by exploiting advanced background removal techniques and multi-view information. These allow us to generate more diverse instance pairs than those obtained by augmentations commonly used in exemplar based approaches. Our method leads to a more suitable representation for the hand shape estimation task and shows a 4.7% reduction in mesh error and a 3.6% improvement in F-score compared to an ImageNet pretrained baseline. We make our benchmark dataset publicly available, to encourage further research into this direction.


title: Non-Newtonian rheology in a capillary tube with varying radius

abstract: The flow through a capillary tube with non-constant radius and where bubbles of yield stress fluid are injected is strongly non-linear. In particular below a finite yield pressure drop, $P_y$, flow is absent, while a singular behaviour is expected above it. In this paper we compute the yield pressure drop statistics and the mean flow rate in two cases: (i) when a single bubble is injected, (ii) when many bubbles are randomly injected in the fluid.


title: One-directional polarization transport in electron/nuclear spin chains with loss and gain

abstract: Understanding the joint dynamics of electron and nuclear spins is central to core concepts in solid-state magnetic resonance - such as spin-lattice relaxation and dynamic nuclear polarization - but a generalization that capitalizes on competing polarization loss and gain channels is still lacking. Here, we theoretically study the non-Hermitian dynamics of hybrid electron/nuclear spin systems in the simultaneous presence of electron spin pumping and spin-lattice relaxation. Focusing on periodic, one-dimensional chains, we find that by adjusting the electron spin pumping to a critical level, it is possible to steer the flow of nuclear polarization to create site-dependent distributions where either end of the array polarizes in opposite ways, irrespective of the initial state. By contrast, we show that ring-like patterns - where the limit nuclear polarization is uniform - exhibit a non-decaying, externally-driven nuclear spin current. Interestingly, cyclic magnetic field modulation can render these processes largely robust to defects in the chain, a response featuring some interesting similarities - and differences - with recent findings in other non-Hermitian physical platforms.


title: WISDOM Project -- IX Giant Molecular Clouds in the Lenticular Galaxy NGC4429: Effects of Shear and Tidal Forces on Clouds

abstract: We present high spatial resolution (12pc) Atacama Large Millimeter/sub-millimeter Array CO(J=3-2) observations of the nearby lenticular galaxy NGC4429. We identify 217 giant molecular clouds within the 450pc radius molecular gas disc. The clouds generally have smaller sizes and masses but higher surface densities and observed linewidths than those of Milky Way disc clouds. An unusually steep size - line width relation and large cloud internal velocity gradients (0.05 - 0.91 km s^-1 pc^-1) and observed Virial parameters (alpha_obs,vir = 4.0) are found, that appear due to internal rotation driven by the background galactic gravitational potential. Removing this rotation, an internal Virial equilibrium appears to be established between the self-gravitational (Usg) and turbulent kinetic (Eturb) energies of each cloud, i.e. alpha_sg,vir=Usg/Eturb = 1.3. However, to properly account for both self and external gravity (shear and tidal forces), we formulate a modified Virial theorem and define an effective Virial parameter alpha_eff,vir = alpha_sg,vir + Usg/Eext (and associated effective velocity dispersion). The NGC4429 clouds then appear to be in a critical state in which the self-gravitational energy and the contribution of external gravity to the cloud's energy budget (Eext) are approximately equal, i.e. Eext/Usg~1. As such, alpha_eff,vir = 2.2 and most clouds are not virialised but remain marginally gravitationally bound. We show this is consistent with the clouds having sizes similar to their tidal radii and being generally radially elongated. External gravity is thus as important as self-gravity to regulate the clouds of NGC4429.


title: The sensitivity of the Deffuant model to measurement error

abstract: Opinion dynamics models have an enormous potential for studying current phenomena such as vaccine hesitancy. Unfortunately, to date, most of the models have little to no empirical validation. One major problem in testing these models against real-world data relates to the difficulties in measuring opinions in ways that map directly to representations in models. Indeed, this kind of measurement is complex in nature and presents more types of measurement error than just classical random noise. Thus, it is crucial to know how these different error types may affect the model's predictions. In this work, we analyze this relationship in the Deffuant model. Starting from the psychometrics literature, we first discuss how opinion measurements are affected by three types of errors: random noise, binning, and distortions (i.e. uneven intervals between scale points). While the first two are known to most of the scientific community, the third one is mostly unknown outside psychometrics. Because of that, we highlight the nature and peculiarities of each of these measurement errors. By simulating these types of error, we show that the Deffuant model is robust to binning but not to noise and distortions. Indeed, if a scale has 4 or more points (like most self-report scales), binning has almost no effect on the final predictions. However, prediction error increases almost linearly with random noise, up to a maximum error of 40%. After reaching this value, increasing the amount of noise does not worsen the prediction. Distortions are most problematic, reaching a maximum prediction error of 80%. Up to now most of the research focused on the properties of the models without analyzing the types of data they may be used with. Here we show that when studying a model, we should also analyze its robustness to these types of measurement error.


title: Global Solutions to an initial boundary problem for the compressible 3-D MHD equations with Navier-slip and perfectly conducting boundary conditions in exterior domains

abstract: An initial boundary value problem for compressible Magnetohydrodynamics (MHD) is considered on an exterior domain (with the first Betti number vanishes) in $R^3$ in this paper. The global existence of smooth solutions near a given constant state for compressible MHD with the boundary conditions of Navier-slip for the velocity filed and perfect conduction for the magnetic field is established. Moreover the explicit decay rate is given. In particular, the results obtained in this paper also imply the global existence of classical solutions for the full compressible Navier-Stokes equations with Navier-slip boundary conditions on exterior domains in three dimensions, which is not available in literature, to the best of knowledge of the authors'.


title: Weighted Sparse Subspace Representation: A Unified Framework for Subspace Clustering, Constrained Clustering, and Active Learning

abstract: Spectral-based subspace clustering methods have proved successful in many challenging applications such as gene sequencing, image recognition, and motion segmentation. In this work, we first propose a novel spectral-based subspace clustering algorithm that seeks to represent each point as a sparse convex combination of a few nearby points. We then extend the algorithm to constrained clustering and active learning settings. Our motivation for developing such a framework stems from the fact that typically either a small amount of labelled data is available in advance; or it is possible to label some points at a cost. The latter scenario is typically encountered in the process of validating a cluster assignment. Extensive experiments on simulated and real data sets show that the proposed approach is effective and competitive with state-of-the-art methods.


title: Generic properties of free boundary problems in plasma physics

abstract: We are concerned with the global bifurcation analysis of positive solutions to free boundary problems arising in plasma physics. We show that in general, in the sense of domain variations, the following alternative holds: either the shape of the branch of solutions resembles the monotone one of the model case of the two-dimensional disk, or it is a continuous simple curve without bifurcation points which ends up at a point where the boundary density vanishes. On the other hand, we deduce a general criterion ensuring the existence of a free boundary in the interior of the domain. Application to a classic nonlinear eigenvalue problem is also discussed.


title: Progressive Spatio-Temporal Bilinear Network with Monte Carlo Dropout for Landmark-based Facial Expression Recognition with Uncertainty Estimation

abstract: Deep neural networks have been widely used for feature learning in facial expression recognition systems. However, small datasets and large intra-class variability can lead to overfitting. In this paper, we propose a method which learns an optimized compact network topology for real-time facial expression recognition utilizing localized facial landmark features. Our method employs a spatio-temporal bilinear layer as backbone to capture the motion of facial landmarks during the execution of a facial expression effectively. Besides, it takes advantage of Monte Carlo Dropout to capture the model's uncertainty which is of great importance to analyze and treat uncertain cases. The performance of our method is evaluated on three widely used datasets and it is comparable to that of video-based state-of-the-art methods while it has much less complexity.


title: Minimax and adaptive tests for detecting abrupt and possibly transitory changes in a Poisson process

abstract: Motivated by applications in cybersecurity and epidemiology, we consider the problem of detecting an abrupt change in the intensity of a Poisson process, characterised by a jump (non transitory change) or a bump (transitory change) from constant. We propose a complete study from the nonasymptotic minimax testing point of view, when the constant baseline intensity is known or unknown. The question of minimax adaptation with respect to each parameter (height, location, length) of the change is tackled, leading to a comprehensive overview of the various minimax separation rate regimes. We exhibit three such regimes and identify the factors of the two phase transitions, by giving the cost of adaptation to each parameter. For each alternative hypothesis, depending on the knowledge or not of each change parameter, we propose minimax or minimax adaptive tests based on linear statistics, close to CUSUM statistics, or quadratic statistics more adapted to the L 2-distance considered in our minimax criteria and typically more powerful in practice, as our simulation study shows. When the change location or length is unknown, our adaptive tests are constructed from a scan aggregation principle combined with Bonferroni or min-p level correction, and a conditioning trick when the baseline intensity is unknown.