title: Combining cytotoxic agents with continuous dose levels in seamless phase I-II clinical trials

abstract: Phase I-II cancer clinical trial designs are intended to accelerate drug development. In cases where efficacy cannot be ascertained in a short period of time, it is common to divide the study in two stages: i) a first stage in which dose is escalated based only on toxicity data and we look for the maximum tolerated dose (MTD) set and ii) a second stage in which we search for the most efficacious dose within the MTD set. Current available approaches in the area of continuous dose levels involve fixing the MTD after stage I and discarding all collected stage I efficacy data. However, this methodology is clearly inefficient when there is a unique patient population present across stages. In this article, we propose a two-stage design for the combination of two cytotoxic agents assuming a single patient population across the entire study. In stage I, conditional escalation with overdose control (EWOC) is used to allocate successive cohorts of patients. In stage II, we employ an adaptive randomization approach to allocate patients to drug combinations along the estimated MTD curve, which is constantly updated. The proposed methodology is assessed with extensive simulations in the context of a real case study.


title: Transition probability and total crossing events in the multi-species asymmetric exclusion process

abstract: We present explicit formulas for total crossing events in the multi-species asymmetric exclusion process ($r$-ASEP) with underlying $U_q(\widehat{\mathfrak{sl}}_{r+1})$ symmetry. In the case of the two-species TASEP these can be derived using an explicit expression for the general transition probability on $\mathbb{Z}$ in terms of a multiple contour integral derived from a nested Bethe ansatz approach. For the general $r$-ASEP we employ a vertex model approach within which the probability of total crossing can be derived from partial symmetrization of an explicit high rank rainbow partition function. In the case of $r$-TASEP, the total crossing probability can be show to reduce to a multiple integral over the product of $r$ determinants. For $2$-TASEP we additionally derive convenient formulas for cumulative total crossing probabilities using Bernoulli-step initial conditions for particles of type 2 and type 1 respectively.

categories: math-ph cond-mat.stat-mech math.CO math.MP math.PR


title: A Next Basket Recommendation Reality Check

abstract: The goal of a next basket recommendation (NBR) system is to recommend items for the next basket for a user, based on the sequence of their prior baskets. Recently, a number of methods with complex modules have been proposed that claim state-of-the-art performance. They rarely look into the predicted basket and just provide intuitive reasons for the observed improvements, e.g., better representation, capturing intentions or relations, etc. We provide a novel angle on the evaluation of next basket recommendation methods, centered on the distinction between repetition and exploration: the next basket is typically composed of previously consumed items (i.e., repeat items) and new items (i.e, explore items). We propose a set of metrics that measure the repeat/explore ratio and performance of NBR models. Using these new metrics, we analyze state-of-the-art NBR models. The results of our analysis help to clarify the extent of the actual progress achieved by existing NBR methods as well as the underlying reasons for the improvements. Overall, our work sheds light on the evaluation problem of NBR and provides useful insights into the model design for this task.


title: Apertif, Phased Array Feeds for the Westerbork Synthesis Radio Telescope

abstract: We describe the APERture Tile In Focus (Apertif) system, a phased array feed (PAF) upgrade of the Westerbork Synthesis Radio Telescope which has transformed this telescope into a high-sensitivity, wide field-of-view L-band imaging and transient survey instrument. Using novel PAF technology, up to 40 partially overlapping beams can be formed on the sky simultaneously, significantly increasing the survey speed of the telescope. With this upgraded instrument, an imaging survey covering an area of 2300 deg2 is being performed which will deliver both continuum and spectral line data sets, of which the first data has been publicly released. In addition, a time domain transient and pulsar survey covering 15,000 deg2 is in progress. An overview of the Apertif science drivers, hardware and software of the upgraded telescope is presented, along with its key performance characteristics.


title: Error rate control for classification rules in multiclass mixture models

abstract: In the context of finite mixture models one considers the problem of classifying as many observations as possible in the classes of interest while controlling the classification error rate in these same classes. Similar to what is done in the framework of statistical test theory, different type I and type II-like classification error rates can be defined, along with their associated optimal rules, where optimality is defined as minimizing type II error rate while controlling type I error rate at some nominal level. It is first shown that finding an optimal classification rule boils down to searching an optimal region in the observation space where to apply the classical Maximum A Posteriori (MAP) rule. Depending on the misclassification rate to be controlled, the shape of the optimal region is provided, along with a heuristic to compute the optimal classification rule in practice. In particular, a multiclass FDR-like optimal rule is defined and compared to the thresholded MAP rules that is used in most applications. It is shown on both simulated and real datasets that the FDR-like optimal rule may be significantly less conservative than the thresholded MAP rule.


title: LightSecAgg: a Lightweight and Versatile Design for Secure Aggregation in Federated Learning

abstract: Secure model aggregation is a key component of federated learning (FL) that aims at protecting the privacy of each user's individual model while allowing for their global aggregation. It can be applied to any aggregation-based FL approach for training a global or personalized model. Model aggregation needs to also be resilient against likely user dropouts in FL systems, making its design substantially more complex. State-of-the-art secure aggregation protocols rely on secret sharing of the random-seeds used for mask generations at the users to enable the reconstruction and cancellation of those belonging to the dropped users. The complexity of such approaches, however, grows substantially with the number of dropped users. We propose a new approach, named LightSecAgg, to overcome this bottleneck by changing the design from "random-seed reconstruction of the dropped users" to "one-shot aggregate-mask reconstruction of the active users via mask encoding/decoding". We show that LightSecAgg achieves the same privacy and dropout-resiliency guarantees as the state-of-the-art protocols while significantly reducing the overhead for resiliency against dropped users. We also demonstrate that, unlike existing schemes, LightSecAgg can be applied to secure aggregation in the asynchronous FL setting. Furthermore, we provide a modular system design and optimized on-device parallelization for scalable implementation, by enabling computational overlapping between model training and on-device encoding, as well as improving the speed of concurrent receiving and sending of chunked masks. We evaluate LightSecAgg via extensive experiments for training diverse models on various datasets in a realistic FL system with large number of users and demonstrate that LightSecAgg significantly reduces the total training time.


title: Insights into the nature of the $P_{cs}(4459)$

abstract: We study the nature of the recently observed $P_{cs}(4459)$ by the LHCb collaboration by employing three methods based on the elastic effective-range expansion and the resulting size of the effective-range, the saturation of the compositeness relation and width of the resonance, and a direct fit to data involving the channels $J/\psi\Lambda$, $\Xi'_c\bar{D}$ and $\Xi_c\bar{D}^*$. We have also considered the addition of a CDD pole but this scenario can be discarded. Our different analyses clearly indicate the molecular nature of the $P_{cs}(4459)$ with a clear $\Xi_c\bar{D}^*$ dominant component. In relation with heavy-quark-spin symmetry our results also favor the actual existence of two resonances with $J=1/2$ (the lighter one) and $3/2$ (the heavier one) in the energy region of the $P_{cs}(4459)$. In the scenario of two-resonance for the $P_{cs}(4459)$, the inclusion of the $\Xi_c'\bar{D}$ channel is required for the their mass splitting and it allows one to determine the spin structures of the two resonances.


title: Phase transitions between confinement and higgs phases in ${\cal N}=1\,\, SU(N_c)$ SQCD with $1\leq N_F\leq N_c-1$ quark flavors

abstract: Considered is 4-dimensional ${\cal N}=1$ supersymmetric $SU(N_c)$ QCD (SQCD) with $1\leq N_F\leq N_c-1$ equal mass quark flavors in the fundamental representation. The gauge invariant order parameter $\rho$ is introduced distinguishing confinement (with $\rho=0$) and higgs (with $\rho\neq 0$) phases. Using a number of independent arguments for different variants of transition between the confinement and higgs regimes, it is shown that transitions between these regimes are not crossovers but the phase transitions. Besides, it is argued that these phase transitions are of the first order. This is opposite to the conclusion of the paper of E. Fradkin and S.H. Shenker [10] that the transition between the confinement and higgs regimes is the crossover, not the phase transition. And although the theories considered in this paper and in [10] are different, an experience shows that there is a widely spread opinion that the conclusion of [10] is applicable to all QCD-like theories: both lattice and continuum, both not supersymmetric and supersymmetric. This opinion is in contradiction with the results of this paper.


title: Spectral flow inside essential spectrum V: on absorbing points of coupling resonances

abstract: Let $H_0$ and $V$ be self-adjoint operators, such that $V$ admits a factorisation $V = F^*JF$ with bounded self-adjoint $J$ and $|H_0|^{1/2}$-compact $F.$ Coupling resonance functions, $r_j(z),$ of the pair $H_0$ and $V$ can be defined as $r_j(z) = - \sigma_j(z) ^{-1},$ where $\sigma_j(z)$ are eigenvalues of the compact-operator valued holomorphic function $F(H_0-z)^{-1} F^*J.$ Taken together, the functions $r_j(z)$ form an infinite-valued holomorphic function on the resolvent set of~$H_0.$ These functions contain a lot of information about the pair $H_0, V$ (this is well-known in the case of rank one $V$). A point $z_0$ of the resolvent set we call \emph{absorbing} if some $r_j(z)$ goes to $\infty$ as $z \to z_0$ along some half-interval. In this note I present some partial results concerning absorbing points of coupling resonances.


title: $N_C$ evolution of the light meson resonances

abstract: The variation of the resonance poles, i.e., the mass and width, originating from changing the values of $N_C$, the number of the QCD colors, provide an intuitive theoretical method to discriminate the interior constituents of hadrons. Many different approaches, including unitarized chiral amplitudes, dispersive methods, quark models, linear-sigma-model-like framework, etc., have been widely used to determine the $N_C$ trajectories of the various resonance pole positions. We focus on the discussions of the light-flavor meson resonances in this work.


title: Potential reconstruction from general power spectrum in single-field inflation

abstract: We suggest a new method to reconstruct, within canonical single-field inflation, the inflaton potential directly from the primordial power spectrum which may deviate significantly from near scale-invariance. Our approach relies on a more generalized slow-roll approximation than the standard one, and can probe the properties of the inflaton potential reliably. We give a few examples for reconstructing potential and discuss the validity of our method.


title: Investigation of meson-meson interaction

abstract: In the framework of the quark delocalization color screening model, we investigate the meson-meson interaction in the fully light four-quark system . The calculation of the effective potentials of all the $S-$wave states shows that for most states the interaction between two vector mesons is attractive; the one between a pseudoscalar meson and a vector meson is repulsive or weakly attractive; and the one between two pseudoscalar mesons is always repulsive. However, there is still some exception. The interaction of the $IJ=00$ $\pi\pi$ channel is attractive, while the one of the $IJ=02$ $\phi\phi$ channel is repulsive. So it is difficult to use the $S-$wave $\phi\phi$ state to explain the $X(2239)$ at present calculation. The $S-$wave $\rho\rho$ states are more likely to be resonance states, which are worthy of investigating in future work. The study of the contribution of each interaction term shows that both the one-gluon exchange and the kinetic energy interaction are important in the interaction between two mesons. The study of the variation of the delocalization parameter indicates that the contribution of the kinetic energy relates to the intermediate-range attraction mechanism in QDCSM, which is achieved by the quark delocalization.


title: DN-ADMM: Distributed Newton ADMM for Multi-agent Optimization

abstract: In a multi-agent network, we consider the problem of minimizing an objective function that is expressed as the sum of private convex and smooth functions, and a (possibly) non-differentiable convex regularizer. We propose a novel distributed second-order method based on the framework of Alternating Direction Method of Multipliers (ADMM), by invoking approximate Newton iterations to the primal update corresponding to the differentiable part. In order to achieve a distributed implementation, the total Hessian matrix is split into a diagonal component (locally computable) and an off-diagonal component (that requires communication between neighboring agents). Subsequently, the Hessian inverse is approximated by a truncation of the Taylor expansion to $K$ terms: this amounts to fully distributed updates entailing $K$ distributed communication rounds. We establish global linear convergence to the primal-dual optimal solution under the assumption that the private functions are strongly convex and have Lipschitz continuous gradient. Numerical experiments demonstrate the merits of the approach comparatively with state-of-the-art methods.


title: Error-mitigated photonic variational quantum eigensolver using a single-photon ququart

abstract: We report the experimental resource-efficient implementation of the variational quantum eigensolver (VQE) using four-dimensional photonic quantum states of single-photons. The four-dimensional quantum states are implemented by utilizing polarization and path degrees of freedom of a single-photon. Our photonic VQE is equipped with the quantum error mitigation (QEM) protocol that efficiently reduces the effects of Pauli noise in the quantum processing unit. We apply our photonic VQE to estimate the ground state energy of He--H$^{+}$ cation. The simulation and experimental results demonstrate that our resource-efficient photonic VQE can accurately estimate the bond dissociation curve, even in the presence of large noise in the quantum processing unit.


title: Effect of quantum deformed black hole on BH shadow in two-dimensional Dilaton gravity

abstract: In recent years, the study of quantum effects near the event horizon of black hole (BH) has attracted extensive attention. It has become one of the important methods to explore BH quantum properties by using the related properties of the quantum deformed black hole. In this work, we study the effect of quantum deformed black hole on BH shadow in two-dimensional Dilaton gravity. In this model, quantum effects are reflected on the quantum correction parameter m. By calculation, we find that: (1) the shape of the shadow boundary of a rotating black hole is determined by the BH spin $a$, the quantum correction parameter $m$ and the BH type parameter $n$; (2) when the spin $a=0$, the shape of the BH shadow is a perfect circle; when $a\neq 0$, the shape is distorted; if the quantum correction parameter $m=0$, their shapes reduce to the cases of Schwarzschild BH and Kerr BH respectively; (3) the degree of distortion of the BH shadow is different for various quantum correction parameters $m$; with the increase of the values of $m$, the shadow will become more and more obvious; (4) the results of different BH type parameter $n$ differ greatly. Since the value of $m$ in actual physics should be very small, the current observations of EHT cannot distinguish quantum effect from BH shadow, and can only constrain the upper limit of $m$. In future BH shadow measurements, it will be possible to distinguish quantum deformed black holes, which will help to better understand the quantum effects of BHs.


title: Localization for one-dimensional Anderson-Dirac models

abstract: We prove spectral and dynamical localization for a one-dimensional Dirac operator to which is added an ergodic random potential, with a discussion on the different types of potential. We use scattering properties to prove the positivity of the Lyapunov exponent through F\"urstenberg theorem. We get then the H\"older regularity of the integrated density of states through a new version of Thouless formula, and thus the Wegner estimate necessary for the multiscale analysis.


title: Training Feedback Spiking Neural Networks by Implicit Differentiation on the Equilibrium State

abstract: Spiking neural networks (SNNs) are brain-inspired models that enable energy-efficient implementation on neuromorphic hardware. However, the supervised training of SNNs remains a hard problem due to the discontinuity of the spiking neuron model. Most existing methods imitate the backpropagation framework and feedforward architectures for artificial neural networks, and use surrogate derivatives or compute gradients with respect to the spiking time to deal with the problem. These approaches either accumulate approximation errors or only propagate information limitedly through existing spikes, and usually require information propagation along time steps with large memory costs and biological implausibility. In this work, we consider feedback spiking neural networks, which are more brain-like, and propose a novel training method that does not rely on the exact reverse of the forward computation. First, we show that the average firing rates of SNNs with feedback connections would gradually evolve to an equilibrium state along time, which follows a fixed-point equation. Then by viewing the forward computation of feedback SNNs as a black-box solver for this equation, and leveraging the implicit differentiation on the equation, we can compute the gradient for parameters without considering the exact forward procedure. In this way, the forward and backward procedures are decoupled and therefore the problem of non-differentiable spiking functions is avoided. We also briefly discuss the biological plausibility of implicit differentiation, which only requires computing another equilibrium. Extensive experiments on MNIST, Fashion-MNIST, N-MNIST, CIFAR-10, and CIFAR-100 demonstrate the superior performance of our method for feedback models with fewer neurons and parameters in a small number of time steps. Our code is avaiable at https://github.com/pkuxmq/IDE-FSNN.


title: EBSD Grain Knowledge Graph Representation Learning for Material Structure-Property Prediction

abstract: The microstructure is an essential part of materials, storing the genes of materials and having a decisive influence on materials' physical and chemical properties. The material genetic engineering program aims to establish the relationship between material composition/process, organization, and performance to realize the reverse design of materials, thereby accelerating the research and development of new materials. However, tissue analysis methods of materials science, such as metallographic analysis, XRD analysis, and EBSD analysis, cannot directly establish a complete quantitative relationship between tissue structure and performance. Therefore, this paper proposes a novel data-knowledge-driven organization representation and performance prediction method to obtain a quantitative structure-performance relationship. First, a knowledge graph based on EBSD is constructed to describe the material's mesoscopic microstructure. Then a graph representation learning network based on graph attention is constructed, and the EBSD organizational knowledge graph is input into the network to obtain graph-level feature embedding. Finally, the graph-level feature embedding is input to a graph feature mapping network to obtain the material's mechanical properties. The experimental results show that our method is superior to traditional machine learning and machine vision methods.


title: Multipath CNN with alpha matte inference for knee tissue segmentation from MRI

abstract: Precise segmentation of knee tissues from magnetic resonance imaging (MRI) is critical in quantitative imaging and diagnosis. Convolutional neural networks (CNNs), which are state of the art, have limitations owing to the lack of image-specific adaptation, such as low tissue contrasts and structural inhomogeneities, thereby leading to incomplete segmentation results. This paper presents a deep learning based automatic segmentation framework for knee tissue segmentation. A novel multipath CNN-based method is proposed, which consists of an encoder decoder-based segmentation network in combination with a low rank tensor-reconstructed segmentation network. Low rank reconstruction in MRI tensor sub-blocks is introduced to exploit the structural and morphological variations in knee tissues. To further improve the segmentation from CNNs, trimap generation, which effectively utilizes superimposed regions, is proposed for defining high, medium and low confidence regions from the multipath CNNs. The secondary path with low rank reconstructed input mitigates the conditions in which the primary segmentation network can potentially fail and overlook the boundary regions. The outcome of the segmentation is solved as an alpha matting problem by blending the trimap with the source input. Experiments on Osteoarthritis Initiative (OAI) datasets and a self prepared scan validate the effectiveness of the proposed method. We specifically demonstrate the application of the proposed method in a cartilage segmentation based thickness map for diagnosis purposes.


title: BLEU, METEOR, BERTScore: Evaluation of Metrics Performance in Assessing Critical Translation Errors in Sentiment-oriented Text

abstract: Social media companies as well as authorities make extensive use of artificial intelligence (AI) tools to monitor postings of hate speech, celebrations of violence or profanity. Since AI software requires massive volumes of data to train computers, Machine Translation (MT) of the online content is commonly used to process posts written in several languages and hence augment the data needed for training. However, MT mistakes are a regular occurrence when translating sentiment-oriented user-generated content (UGC), especially when a low-resource language is involved. The adequacy of the whole process relies on the assumption that the evaluation metrics used give a reliable indication of the quality of the translation. In this paper, we assess the ability of automatic quality metrics to detect critical machine translation errors which can cause serious misunderstanding of the affect message. We compare the performance of three canonical metrics on meaningless translations where the semantic content is seriously impaired as compared to meaningful translations with a critical error which exclusively distorts the sentiment of the source text. We conclude that there is a need for fine-tuning of automatic metrics to make them more robust in detecting sentiment critical errors.


title: Road Network Guided Fine-Grained Urban Traffic Flow Inference

abstract: Accurate inference of fine-grained traffic flow from coarse-grained one is an emerging yet crucial problem, which can help greatly reduce the number of the required traffic monitoring sensors for cost savings. In this work, we notice that traffic flow has a high correlation with road network, which was either completely ignored or simply treated as an external factor in previous works. To facilitate this problem, we propose a novel Road-Aware Traffic Flow Magnifier (RATFM) that explicitly exploits the prior knowledge of road networks to fully learn the road-aware spatial distribution of fine-grained traffic flow. Specifically, a multi-directional 1D convolutional layer is first introduced to extract the semantic feature of the road network. Subsequently, we incorporate the road network feature and coarse-grained flow feature to regularize the short-range spatial distribution modeling of road-relative traffic flow. Furthermore, we take the road network feature as a query to capture the long-range spatial distribution of traffic flow with a transformer architecture. Benefiting from the road-aware inference mechanism, our method can generate high-quality fine-grained traffic flow maps. Extensive experiments on three real-world datasets show that the proposed RATFM outperforms state-of-the-art models under various scenarios. Our code and datasets are released at {\url{https://github.com/luimoli/RATFM}}.


title: A Comptonized Fireball Bubble: Physical Origin of Magnetar Giant Flares

abstract: Magnetar giant flares (MGFs) have been long proposed to contribute at least a sub-sample of the observed short gamma-ray bursts (GRBs). The recent discovery of the short GRB 200415A in the nearby galaxy NGC 253 established a textbook-version connection between these two phenomena. Unlike previous observations of the Galactic MGFs, the unsaturated instrument spectra of GRB 200415A provide for the first time an opportunity to test the theoretical models with the observed $\gamma$-ray photons. This paper proposed a new readily fit-able model for the MGFs, which invokes an expanding fireball Comptonized by the relativistic magnetar wind at photosphere radius. In this model, a large amount of energy is released from the magnetar crust due to the magnetic reconnection or the starquakes of the star surface and is injected into confined field lines, forming a trapped fireball bubble. After breaking through the shackles and expanding to the photospheric radius, the thermal photons of the fireball are eventually Comptonized by the relativistic $e^{\pm}$ pairs in the magnetar wind region, which produces additional higher-energy gamma-ray emission. The model predicts a modified thermal-like spectrum characterized by a low-energy component in the Rayleigh-Jeans regime, a smooth component affected by coherent Compton scattering (CC) in the intermediate energy range, and a high-energy tail due to the inverse Compton process. By performing a Monte-Carlo fit to the observational spectra of GRB 200415A, we found that the observation of the burst is entirely consistent with our model predictions.


title: Density dependent local structures in InTe phase-change materials

abstract: Chalcogenide phase-change materials (PCMs) based random access memory (PCRAM) is one of the leading candidates for the development of non-volatile memory and neuro-inspired computing technologies. Recent work shows Indium to be an important alloying element for PCRAM, while a thorough understanding of the parent compound InTe, in particular, its amorphous phase, is still lacking. In this work, we carry out ab initio simulations and chemical bonding analyses on amorphous and various crystalline polymorphs of InTe. We reveal that the local geometries are highly density dependent in amorphous structures, forming In-centered tetrahedral motifs under ambient conditions but defective octahedral motifs under pressure, which stems from the bonding characters of its crystalline polymorphs. In addition, our ab initio molecular dynamics simulations predict rapid crystallization capability of InTe under pressure. At last, we make a suggestion for better use of Indium and propose an "active" device design to utilize both thermal and mechanical effects for phase-change applications.


title: Non-local parameter degeneracy in the intrinsic space of gravitational-wave signals from extreme-mass-ratio inspirals

abstract: Extreme-mass-ratio inspirals will be prized sources for the upcoming space-based gravitational-wave observatory LISA. The hunt for these is beset by many open theoretical and computational problems in both source modeling and data analysis. We draw attention here to one of the most poorly understood: the phenomenon of non-local correlations in the space of extreme-mass-ratio-inspiral signals. Such correlations are ubiquitous in the continuum of possible signals (degeneracy), and severely hinder the search for actual signals in LISA data. However, they are unlikely to manifest in a realistic set of putative signals (confusion). We develop an inventory of new analysis tools in order to conduct an extensive qualitative study of degeneracy - its nature, causes, and implications. Previously proposed search strategies for extreme-mass-ratio inspirals are reviewed in the light of our results, and additional guidelines are suggested for the scientific analysis of such sources.


title: Functional inequalities and applications to doubly nonlinear diffusion equations

abstract: We study weighted inequalities of Hardy and Hardy-Poincar\'e type and find necessary and sufficient conditions on the weights so that the considered inequalities hold. Examples with the optimal constants are shown. Such inequalities are then used to quantify the convergence rate of solutions to doubly nonlinear fast diffusion equation towards the Barenblatt profile.


title: Lang--Trotter Conjecture for CM Elliptic Curves

abstract: Given an elliptic curve $E$ over $\mathbb{Q}$ and non-zero integer $r$, the Lang--Trotter conjecture predicts a striking asymptotic formula for the number of good primes $p\leqslant x$, denoted by $\pi_{E,r}(x)$, such that the Frobenius trace of $E$ at $p$ is equal to the given integer $r$. We focus on the CM case in this memoir, and show how to realize the following two goals: (1) to give an unconditional estimate for $\pi_{E,r}(x)$, which confirms the upper bound part of the conjecture up to a constant multiple; (2) to give a conditional explicit asymptotic formula for $\pi_{E,r}(x)$ based on the Hardy--Littlewood conjecture on primes represented by quadratic polynomials. For completeness, we also summarize classical results on quadratic, cubic and quartic residues, as well as the corresponding reciprocity laws. This part should be of independent interests and could provide useful materials for more junior readers. We also highlight some possible extensions of the arguments in this memoir that may work for other statistical problems of CM elliptic curves.


title: Adaptive-Resolution Field Mapping Using Gaussian Process Fusion with Integral Kernels

abstract: Unmanned aerial vehicles are rapidly gaining popularity in a variety of environmental monitoring tasks. A key requirement for their autonomous operation is the ability to perform efficient environmental mapping online, given limited onboard resources constraining operation time, travel distance, and computational capacity. To address this, we present an online adaptive-resolution approach for mapping terrain based on Gaussian Process fusion. A key aspect of our approach is an integral kernel encoding spatial correlation over the areas of grid cells, which enables modifying map resolution while maintaining correlations in a theoretically sound fashion. This way, we can retain details in areas of interest at higher map resolutions while compressing information in uninteresting areas at coarser resolutions to achieve a compact map representation of the environment. We evaluate the performance of our approach on both synthetic and real-world data. Results show that our method is more efficient in terms of mapping time and memory consumption without compromising on map quality. Finally, we integrate our mapping strategy into an adaptive path planning framework to show that it facilitates information gathering efficiency in online settings.


title: Energy momentum tensor of a non-minimally coupled scalar from the equivalence of the Einstein and Jordan frames

abstract: Unlike the minimally coupled gravity theory where matter is coupled with gravity in such a manner so that one can differentiate the matter and gravity sector uniquely, the non-minimally coupled theories (NMCT) are distinguished by the intermingling of two. As a consequence of this the calculation of the energy momentum tensor (EMT) in NMCT is beset with an arbitrariness. In this paper we provide an algorithm based on the well known equivalence between Jordan frame and Einstein frame formulations which enables us to construct the EMT for NMCT in a unique way.


title: Hierarchical Character Tagger for Short Text Spelling Error Correction

abstract: State-of-the-art approaches to spelling error correction problem include Transformer-based Seq2Seq models, which require large training sets and suffer from slow inference time; and sequence labeling models based on Transformer encoders like BERT, which involve token-level label space and therefore a large pre-defined vocabulary dictionary. In this paper we present a Hierarchical Character Tagger model, or HCTagger, for short text spelling error correction. We use a pre-trained language model at the character level as a text encoder, and then predict character-level edits to transform the original text into its error-free form with a much smaller label space. For decoding, we propose a hierarchical multi-task approach to alleviate the issue of long-tail label distribution without introducing extra model parameters. Experiments on two public misspelling correction datasets demonstrate that HCTagger is an accurate and much faster approach than many existing models.


title: Combinatorial Contracts

abstract: We introduce a new model of combinatorial contracts in which a principal delegates the execution of a costly task to an agent. To complete the task, the agent can take any subset of a given set of unobservable actions, each of which has an associated cost. The cost of a set of actions is the sum of the costs of the individual actions, and the principal's reward as a function of the chosen actions satisfies some form of diminishing returns. The principal incentivizes the agents through a contract, based on the observed outcome. Our main results are for the case where the task delegated to the agent is a project, which can be successful or not. We show that if the success probability as a function of the set of actions is gross substitutes, then an optimal contract can be computed with polynomially many value queries, whereas if it is submodular, the optimal contract is NP-hard. All our results extend to linear contracts for higher-dimensional outcome spaces, which we show to be robustly optimal given first moment constraints. Our analysis uncovers a new property of gross substitutes functions, and reveals many interesting connections between combinatorial contracts and combinatorial auctions, where gross substitutes is known to be the frontier for efficient computation.


title: Radiatively driven, time dependent bipolar outflows

abstract: We study the radiatively driven fluid jets around a non-rotating black hole. The radiation arising from the inner compact corona and outer sub-Keplerian part of the disc accelerates the jets. We obtain the steady state, semi-analytical, radiatively driven outflow solutions. The thermodynamics of the outflow is described by a variable adiabatic index equation of state. We develop a TVD routine to investigate the time dependent behaviour of the radiatively driven bipolar outflow. We inject with flow variables from the steady state outflow solutions in the TVD code and allow the code to settle into steady state and match the numerical results with the steady state solution. The radiation arising out of the accretion disc can provide a wide range of jet solutions, depending upon parameters like the intensity of disc, location of the inner corona etc. We induce the time dependence of the radiation field by inducing oscillation of the inner corona of the accretion disc. The radiation field then makes the bipolar outflow time dependent. We show that a non-steady radiation field arising out of disc oscillations can generate the internal shocks closer to the jet base. Depending on the disc geometry, there might be transient shocks in the jet and there might be multiple non-stationary shocks in the jet, which are of much interest in jet physics.


title: Deep-learning-driven event reconstruction applied to simulated data from a single Large-Sized Telescope of CTA

abstract: When very-high-energy gamma rays interact high in the Earth's atmosphere, they produce cascades of particles that induce flashes of Cherenkov light. Imaging Atmospheric Cherenkov Telescopes (IACTs) detect these flashes and convert them into shower images that can be analyzed to extract the properties of the primary gamma ray. The dominant background for IACTs is comprised of air shower images produced by cosmic hadrons, with typical noise-to-signal ratios of several orders of magnitude. The standard technique adopted to differentiate between images initiated by gamma rays and those initiated by hadrons is based on classical machine learning algorithms, such as Random Forests, that operate on a set of handcrafted parameters extracted from the images. Likewise, the inference of the energy and the arrival direction of the primary gamma ray is performed using those parameters. State-of-the-art deep learning techniques based on convolutional neural networks (CNNs) have the potential to enhance the event reconstruction performance, since they are able to autonomously extract features from raw images, exploiting the pixel-wise information washed out during the parametrization process. Here we present the results obtained by applying deep learning techniques to the reconstruction of Monte Carlo simulated events from a single, next-generation IACT, the Large-Sized Telescope (LST) of the Cherenkov Telescope Array (CTA). We use CNNs to separate the gamma-ray-induced events from hadronic events and to reconstruct the properties of the former, comparing their performance to the standard reconstruction technique. Three independent implementations of CNN-based event reconstruction models have been utilized in this work, producing consistent results.


title: Cooperative Task Offloading and Block Mining in Blockchain-based Edge Computing with Multi-agent Deep Reinforcement Learning

abstract: The convergence of mobile edge computing (MEC) and blockchain is transforming the current computing services in mobile networks, by offering task offloading solutions with security enhancement empowered by blockchain mining. Nevertheless, these important enabling technologies have been studied separately in most existing works. This article proposes a novel cooperative task offloading and block mining (TOBM) scheme for a blockchain-based MEC system where each edge device not only handles data tasks but also deals with block mining for improving the system utility. To address the latency issues caused by the blockchain operation in MEC, we develop a new Proof-of-Reputation consensus mechanism based on a lightweight block verification strategy. A multi-objective function is then formulated to maximize the system utility of the blockchain-based MEC system, by jointly optimizing offloading decision, channel selection, transmit power allocation, and computational resource allocation. We propose a novel distributed deep reinforcement learning-based approach by using a multi-agent deep deterministic policy gradient algorithm. We then develop a game-theoretic solution to model the offloading and mining competition among edge devices as a potential game, and prove the existence of a pure Nash equilibrium. Simulation results demonstrate the significant system utility improvements of our proposed scheme over baseline approaches.


title: On monohedral tilings of a regular polygon

abstract: A tiling of a topological disc by topological discs is called monohedral if all tiles are congruent. Maltby (J. Combin. Theory Ser. A 66: 40-52, 1994) characterized the monohedral tilings of a square by three topological discs. Kurusa, L\'angi and V\'\i gh (Mediterr. J. Math. 17: article number 156, 2020) characterized the monohedral tilings of a circular disc by three topological discs. The aim of this note is to connect these two results by characterizing the monohedral tilings of any regular $n$-gon with at most three tiles for any $n \geq 5$.


title: Majority Vote in Social Networks: Make Random Friends or Be Stubborn to Overpower Elites

abstract: Consider a graph $G$, representing a social network. Assume that initially each node is colored either black or white, which corresponds to a positive or negative opinion regarding a consumer product or a technological innovation. In the majority model, in each round all nodes simultaneously update their color to the most frequent color among their connections. Experiments on the graph data from the real world social networks (SNs) suggest that if all nodes in an extremely small set of high-degree nodes, often referred to as the elites, agree on a color, that color becomes the dominant color at the end of the process. We propose two countermeasures that can be adopted by individual nodes relatively easily and guarantee that the elites will not have this disproportionate power to engineer the dominant output color. The first countermeasure essentially requires each node to make some new connections at random while the second one demands the nodes to be more reluctant towards changing their color (opinion). We verify their effectiveness and correctness both theoretically and experimentally. We also investigate the majority model and a variant of it when the initial coloring is random on the real world SNs and several random graph models. In particular, our results on the Erd\H{o}s-R\'{e}nyi and regular random graphs confirm or support several theoretical findings or conjectures by the prior work regarding the threshold behavior of the process. Finally, we provide theoretical and experimental evidence for the existence of a poly-logarithmic bound on the expected stabilization time of the majority model.


title: n-Qubit Operations on Sphere and Queueing Scaling Limits for Programmable Quantum Computer

abstract: We study n-qubit operation rules on (n+1)-sphere with the target to help developing a (photon or other technique) based programmable quantum computer. In the meanwhile, we derive the scaling limits (called reflecting Gaussian random fields on a (n+1)-sphere) for n-qubit quantum computer based queueing systems under two different heavy traffic regimes. The queueing systems are with multiple classes of users and batch quantum random walks over the $(n+1)$-sphere as arrival inputs. In the first regime, the qubit number $n$ is fixed and the scaling is in terms of both time and space. Under this regime, performance modeling during deriving the scaling limit in terms of balancing the arrival and service rates under first-in first-out and work conserving service policy is conducted. In the second regime, besides the time and space scaling parameters, the qubit number $n$ itself is also considered as a varying scaling parameter with the additional aim to find a suitable number of qubits for the design of a quantum computer. This regime is in contrast to the well-known Halfin-Whitt regime.

categories: quant-ph cs.ET cs.IT cs.PF math.IT math.PR


title: Real Space Quantum Cluster Formulation for the Typical Medium Theory of Anderson Localization

abstract: We develop a real space cluster extension of the typical medium theory (cluster-TMT) to study Anderson localization. By construction, the cluster-TMT approach is formally equivalent to the real space cluster extension of the dynamical mean field theory. Applying the developed method to the 3D Anderson model with a box disorder distribution, we demonstrate that cluster-TMT successfully captures the localization phenomena in all disorder regimes. As a function of the cluster size, our method obtains the correct critical disorder strength for the Anderson localization in 3D, and systematically recovers the re-entrance behavior of the mobility edge. From a general perspective, our developed methodology offers the potential to study Anderson localization at surfaces within quantum embedding theory.


title: Formulation and validation of a car-following model based on deep reinforcement learning

abstract: We propose and validate a novel car following model based on deep reinforcement learning. Our model is trained to maximize externally given reward functions for the free and car-following regimes rather than reproducing existing follower trajectories. The parameters of these reward functions such as desired speed, time gap, or accelerations resemble that of traditional models such as the Intelligent Driver Model (IDM) and allow for explicitly implementing different driving styles. Moreover, they partially lift the black-box nature of conventional neural network models. The model is trained on leading speed profiles governed by a truncated Ornstein-Uhlenbeck process reflecting a realistic leader's kinematics. This allows for arbitrary driving situations and an infinite supply of training data. For various parameterizations of the reward functions, and for a wide variety of artificial and real leader data, the model turned out to be unconditionally string stable, comfortable, and crash-free. String stability has been tested with a platoon of five followers following an artificial and a real leading trajectory. A cross-comparison with the IDM calibrated to the goodness-of-fit of the relative gaps showed a higher reward compared to the traditional model and a better goodness-of-fit.


title: Dynamic Emotions of Supporters and Opponents of Anti-racism Movement from George Floyd Protests

abstract: Social media empowers citizens to raise the voice and expressed civil outrage leads to collective action to change the society. Since social media welcomes anyone regardless of the political ideology or perspectives, social media is where the supporters and opponents of specific issue discuss. This study attempts to empirically examine a recent anti-racism movement initiated by the death of George Floyd with the lens of stance prediction and aspect-based sentiment analysis (ABSA). First, this study found the stance of the tweet and users do change over the course of the protest. Furthermore, there are more users who shifted the stance compared to those who maintained the stance. Second, both supporters and opponents expressed negative sentiment more on nine extracted aspects. This indicates that there was no significant difference of sentiment among supporters and opponents and raise a caution in predicting stance based on the sentiment. The contribution of the study is two-fold. First, ABSA was explored in the context of computational social science and second, stance prediction was first attempted at scale.


title: In the Search for the Infinite Servers Queue with Poisson Arrivals Busy Period Distribution Exponential Behavior

abstract: This paper purpose is to investigate exponential behavior conditions for the infinite servers queue with Poisson arrivals busy period length distribution. It is presented a general theoretical result that is the basis of this work. The complementary analysis rely on the infinite servers queue with Poisson arrivals busy period length distribution moments computation. In infinite servers queue with Poisson arrivals practical applications, in economic, management, and business areas, the management of the effective number of servers is essential since the physical presence of infinite servers is not viable and so it is necessary to create that condition through an adequate management of the number of servers during the busy period.


title: (Machine) Learning to Improve the Empirical Performance of Discrete Algorithms

abstract: This paper discusses a data-driven, empirically-based framework to make algorithmic decisions or recommendations without expert knowledge. We improve the performance of two algorithmic case studies: the selection of a pivot rule for the Simplex method and the selection of an all-pair shortest paths algorithm. We train machine learning methods to select the optimal algorithm for given data without human expert opinion. We use two types of techniques, neural networks and boosted decision trees. We concluded, based on our experiments, that: 1) Our selection framework recommends various pivot rules that improve overall total performance over just using a fixed default pivot rule. Over many years experts identified steepest-edge pivot rule as a favorite pivot rule. Our data analysis corroborates that the number of iterations by steepest-edge is no more than 4 percent more than the optimal selection which corroborates human expert knowledge, but this time the knowledge was obtained using machine learning. Here our recommendation system is best when using gradient boosted trees. 2) For the all-pairs shortest path problem, the models trained made a large improvement and our selection is on average .07 percent away from the optimal choice. The conclusions do not seem to be affected by the machine learning method we used. We tried to make a parallel analysis of both algorithmic problems, but it is clear that there are intrinsic differences. For example, in the all-pairs shortest path problem the graph density is a reasonable predictor, but there is no analogous single parameter for decisions in the Simplex method.


title: Computing Necessary Conditions for Near-Optimality in Capacity Expansion Planning Problems

abstract: In power systems, large-scale optimisation problems are extensively used to plan for capacity expansion at the supra-national level. However, their cost-optimal solutions are often not exploitable by decision-makers who are preferably looking for features of solutions that can accommodate their different requirements. This paper proposes a generic framework for addressing this problem. It is based on the concept of the epsilon-optimal feasible space of a given optimisation problem and the identification of necessary conditions over this space. This framework has been developed in a generic case, and an approach for solving this problem is subsequently described for a specific case where conditions are constrained sums of variables. The approach is tested on a case study about capacity expansion planning of the European electricity network to determine necessary conditions on the minimal investments in transmission, storage and generation capacity.


title: Scale invariance of electrodynamics in radio-frequency linear accelerators

abstract: This study discusses the scale transformation of electromagnetic fields and trajectories of ions in radio-frequency linear accelerators. We will obtain a condition for mechanical similarity, where Maxwell equations with source terms and the equation of motion under Lorentz force remain invariant, including space charge fields, while the trajectories of ions are accordingly scaled. The possibility of linear accelerators for extremely high current beams will be considered based on this condition.


title: Designing Counterfactual Generators using Deep Model Inversion

abstract: Explanation techniques that synthesize small, interpretable changes to a given image while producing desired changes in the model prediction have become popular for introspecting black-box models. Commonly referred to as counterfactuals, the synthesized explanations are required to contain discernible changes (for easy interpretability) while also being realistic (consistency to the data manifold). In this paper, we focus on the case where we have access only to the trained deep classifier and not the actual training data. While the problem of inverting deep models to synthesize images from the training distribution has been explored, our goal is to develop a deep inversion approach to generate counterfactual explanations for a given query image. Despite their effectiveness in conditional image synthesis, we show that existing deep inversion methods are insufficient for producing meaningful counterfactuals. We propose DISC (Deep Inversion for Synthesizing Counterfactuals) that improves upon deep inversion by utilizing (a) stronger image priors, (b) incorporating a novel manifold consistency objective and (c) adopting a progressive optimization strategy. We find that, in addition to producing visually meaningful explanations, the counterfactuals from DISC are effective at learning classifier decision boundaries and are robust to unknown test-time corruptions.


title: Simulation-based Bayesian inference for multi-fingered robotic grasping

abstract: Multi-fingered robotic grasping is an undeniable stepping stone to universal picking and dexterous manipulation. Yet, multi-fingered grippers remain challenging to control because of their rich nonsmooth contact dynamics or because of sensor noise. In this work, we aim to plan hand configurations by performing Bayesian posterior inference through the full stochastic forward simulation of the robot in its environment, hence robustly accounting for many of the uncertainties in the system. While previous methods either relied on simplified surrogates of the likelihood function or attempted to learn to directly predict maximum likelihood estimates, we bring a novel simulation-based approach for full Bayesian inference based on a deep neural network surrogate of the likelihood-to-evidence ratio. Hand configurations are found by directly optimizing through the resulting amortized and differentiable expression for the posterior. The geometry of the configuration space is accounted for by proposing a Riemannian manifold optimization procedure through the neural posterior. Simulation and physical benchmarks demonstrate the high success rate of the procedure.


title: Sequential Deep Learning Architectures for Anomaly Detection in Virtual Network Function Chains

abstract: Software-defined networking (SDN) and network function virtualization (NFV) have enabled the efficient provision of network service. However, they also raised new tasks to monitor and ensure the status of virtualized service, and anomaly detection is one of such tasks. There have been many data-driven approaches to implement anomaly detection system (ADS) for virtual network functions in service function chains (SFCs). In this paper, we aim to develop more advanced deep learning models for ADS. Previous approaches used learning algorithms such as random forest (RF), gradient boosting machine (GBM), or deep neural networks (DNNs). However, these models have not utilized sequential dependencies in the data. Furthermore, they are limited as they can only apply to the SFC setting from which they were trained. Therefore, we propose several sequential deep learning models to learn time-series patterns and sequential patterns of the virtual network functions (VNFs) in the chain with variable lengths. As a result, the suggested models improve detection performance and apply to SFCs with varying numbers of VNFs.


title: Red- and blue-detuned magneto-optical trapping with liquid crystal variable retarders

abstract: We exploit red- and blue-detuned magneto optical trapping (MOT) of $^{87}$Rb benefitting from a simplified setup based on liquid crystal variable retarders (LCVR). To maintain the trapping forces when switching from a red- to a blue-detuned MOT, the circularity of the cooling beams needs to be reversed. LCVRs allow fast polarization control and represent compact, simple and cost-efficient components which can easily be implemented in existing laser systems. This way, we achieve a blue-detuned type-II MOT for $^{87}$Rb atoms with sub-Doppler temperatures of 44.5 $\mu$K. The phase space density is increased by more than two orders of magnitude compared to the standard red-detuned type-I MOT. The setup can readily be transferred to any other system working with $^{87}$Rb.


title: 4-dimensional Chern-Simons theory and integrable field theories

abstract: These lecture notes concern the semi-holomorphic 4d Chern-Simons theory and its applications to classical integrable field theories in 2d and in particular integrable sigma-models. After introducing the main properties of the Chern-Simons theory in 3d, we will define its 4d analogue and explain how it is naturally related to the Lax formalism of integrable 2d theories. Moreover, we will explain how varying the boundary conditions imposed on this 4d theory allows to recover various occurences of integrable sigma-models through this construction, in particular illustrating this on two simple examples: the Principal Chiral Model and its Yang-Baxter deformation. These notes were written for the lectures delivered at the school "Integrability, Dualities and Deformations", that ran from 23 to 27 August 2021 in Santiago de Compostela and virtually.


title: Localizing Objects with Self-Supervised Transformers and no Labels

abstract: Localizing objects in image collections without supervision can help to avoid expensive annotation campaigns. We propose a simple approach to this problem, that leverages the activation features of a vision transformer pre-trained in a self-supervised manner. Our method, LOST, does not require any external object proposal nor any exploration of the image collection; it operates on a single image. Yet, we outperform state-of-the-art object discovery methods by up to 8 CorLoc points on PASCAL VOC 2012. We also show that training a class-agnostic detector on the discovered objects boosts results by another 7 points. Moreover, we show promising results on the unsupervised object discovery task. The code to reproduce our results can be found at https://github.com/valeoai/LOST.


title: Introduction to classical and quantum integrability

abstract: In these lecture notes we aim for a pedagogical introduction to both classical and quantum integrability. Starting from Liouville integrability and passing through Lax pair and r-matrix we discuss the construction of the conserved charges for classical integrable models taking as example the harmonic oscillator. The construction of these charges for 2D integrable field theories is also discussed using a Lax connection and the Sine-Gordon model as example. On the quantum side, the XXZ spin chain is used to explain the systematic construction of the conserved charges starting from a quantum R-matrix, solution of the quantum Yang-Baxter equation. The diagonalization of these charges is performed using the algebraic Bethe ansatz. At the end, the interpretation of the R-matrix as an S-matrix in a scattering process is also presented. These notes were written for the lectures delivered at the school "Integrability, Dualities and Deformations", that ran from 23 to 27 August 2021 in Santiago de Compostela and virtually.


title: An infinite class of Neumaier graphs and non-existence results

abstract: A Neumaier graph is a non-complete edge-regular graph containing a regular clique. A Neumaier graph that is not strongly regular is called a strictly Neumaier graph. In this work we present a new construction of strictly Neumaier graphs, and using Jacobi sums, we show that our construction produces infinitely many instances. Moreover, we prove some necessary conditions for the existence of (strictly) Neumaier graphs that allow us to show that several parameter sets are not admissible.


title: A gradient-based variable selection for binary classification in reproducing kernel Hilbert space

abstract: Variable selection is essential in high-dimensional data analysis. Although various variable selection methods have been developed, most rely on the linear model assumption. This article proposes a nonparametric variable selection method for the large-margin classifier defined by reproducing the kernel Hilbert space (RKHS). we propose a gradient-based representation of the large-margin classifier and then regularize the gradient functions by the group-lasso penalty to obtain sparse gradients that naturally lead to the variable selection. The groupwise-majorization-decent algorithm (GMD, Yang and Zou, 2015) is proposed to efficiently solve the proposed problem with a large number of parameters. We employ the strong sequential rule (Tibshirani et al., 2012) to facilitate the tuning procedure. The selection consistency of the proposed method is established by obtaining the risk bound of the estimated classifier and its gradient. Finally, we demonstrate the promising performance of the proposed method through simulations and real data illustration.


title: Convectively coupled equatorial trapped waves in stars and planets

abstract: In this paper, we have studied the convectively coupled equatorially trapped waves in rotating stars, with and without magnetic field. The equatorial trapped HD and MHD Poincar\'e, Rossby, mixed Rossby-Poincar\'e, and Kelvin waves were identified. The effects of stratification and non-traditional Coriolis force terms have been investigated. When the flow is strongly stratified, the wave frequencies of the convectively coupled model are almost the same as those of shallow water model. However, when the flow is weakly stratified, the wave frequencies are constrained by the buoyancy frequency. The non-traditional Coriolis terms affect the widths and phases of the equatorial waves. The width increases with the increasing non-traditional Coriolis parameter. Phase shift occurs when the non-traditional Coriolis parameter is included. Magnetic effect is significant when the magnetic field is strong. We have applied the model in the solar atmosphere and solar tachocline to explain the Rieger type periodicities. For the solar atmosphere, when magnetic effect is taken into account, we find that the magnetic field should be smaller than $5G$ in the solar photosphere. Otherwise, the Rieger type periodicities can be only attributed to long Rossby waves. For the solar tachocline, we find that magnetic field of the solar tachocline should be smaller than $50kG$ to observe the 160 days Rieger period. In addition, we find that the effect of the non-traditional Coriolis terms is not obvious in the solar photosphere, but its effect on the tachocline is significant.

categories: astro-ph.SR astro-ph.EP physics.ao-ph physics.flu-dyn physics.geo-ph


title: Integrable Deformations of Sigma Models

abstract: In this pedagogical review we introduce systematic approaches to deforming integrable 2-dimensional sigma models. We use the integrable principal chiral model and the conformal Wess-Zumino-Witten model as our starting points and explore their Yang-Baxter and current-current deformations. There is an intricate web of relations between these models based on underlying algebraic structures and worldsheet dualities, which is highlighted throughout. We finish with a discussion of the generalisation to other symmetric integrable models, including some original results related to Z T cosets and their deformations, and the application to string theory. This review is based on notes written for lectures delivered at the school "Integrability, Dualities and Deformations", which ran from 23 to 27 August 2021 in Santiago de Compostela and virtually.


title: Be Confident! Towards Trustworthy Graph Neural Networks via Confidence Calibration

abstract: Despite Graph Neural Networks (GNNs) have achieved remarkable accuracy, whether the results are trustworthy is still unexplored. Previous studies suggest that many modern neural networks are over-confident on the predictions, however, surprisingly, we discover that GNNs are primarily in the opposite direction, i.e., GNNs are under-confident. Therefore, the confidence calibration for GNNs is highly desired. In this paper, we propose a novel trustworthy GNN model by designing a topology-aware post-hoc calibration function. Specifically, we first verify that the confidence distribution in a graph has homophily property, and this finding inspires us to design a calibration GNN model (CaGCN) to learn the calibration function. CaGCN is able to obtain a unique transformation from logits of GNNs to the calibrated confidence for each node, meanwhile, such transformation is able to preserve the order between classes, satisfying the accuracy-preserving property. Moreover, we apply the calibration GNN to self-training framework, showing that more trustworthy pseudo labels can be obtained with the calibrated confidence and further improve the performance. Extensive experiments demonstrate the effectiveness of our proposed model in terms of both calibration and accuracy.


title: Conflicting Privacy Preference Signals in the Wild

abstract: Privacy preference signals allow users to express preferences over how their personal data is processed. These signals become important in determining privacy outcomes when they reference an enforceable legal basis, as is the case with recent signals such as the Global Privacy Control and the Transparency & Consent Framework. However, the coexistence of multiple privacy preference signals creates ambiguity as users may transmit more than one signal. This paper collects evidence about ambiguity flowing from the aforementioned two signals and the historic Do Not Track signal. We provide the first empirical evidence that ambiguous signals are sent by web users in the wild. We also show that preferences stored in the browser are reliable predictors of privacy preferences expressed in web dialogs. Finally, we provide the first evidence that popular cookie dialogs are blocked by the majority of users who adopted the Do Not Track and Global Privacy Control standards. These empirical results inform forthcoming legal debates about how to interpret privacy preference signals.


title: Multi-mission view of extragalactic black hole X-ray binaries LMC X-1 and LMC X-3: evolution of broadband spectral features

abstract: Extragalactic black hole X-ray binaries LMC X-1 and LMC X-3 are the persistent sources, usually found in soft spectral state. We present the results from spectral and timing study of these sources using MAXI, NICER, NuSTAR and AstroSat observations carried out during 2014-2020. Study of long-term MAXI lightcurve shows that the fractional variability of flux in 2-10keV is moderate ($\sim$20%) in LMCX-1 and high ($\sim$50%) in LMCX-3. The energy spectra of both sources are characterized by a disc blackbody and a Comptonization component with LMCX-1 having an additional Fe-line emission feature. NICER (0.3-10keV), NuSTAR (3.0-40keV) and AstroSat (0.5-20keV) spectral analysis collectively show that LMC X-1 remained in the soft state (disc flux contribution $f_{disc}>80$%, photon index$\Gamma\sim2.06-4.08$) throughout 2014-2020. Mass accretion rate, $\dot{M}$ of LMC X-1 calculated from bolometric luminosity (0.1-50keV) is found to be within $0.07-0.24\dot{M}_{Edd}$ (Eddington mass accretion rate). Although LMC X-3 remained in the soft state ($f_{disc}>95\%,\Gamma\sim2.3$) during most of the time, it exhibits transition into intermediate ($f_{disc}=47-73\%,\Gamma\sim2.02-2.36$) and hard state ($f_{disc}\sim26\%,\Gamma\sim1.6$). $\dot{M}$ of LMC X-3 through different spectral states varies within $0.01-0.42\dot{M}_{Edd}$. Temporal study show that the Power Density Spectra (PDS) in 0.3-10keV follow a red-noise with rms of 2% for LMC X-1 and in case of LMC X-3, it is 0.08-2.35% during the soft state, but relatively high in the intermediate(3.05-4.91%) and hard state($\sim$17.06%). From continuum-fitting method we constrain spin of LMC X-1 and LMC X-3 to be within 0.85-0.94 and 0.16-0.33 respectively and from Fe-line fitting method, spin of LMC X-1 is found to be 0.93-0.94. Finally, we discuss the implication of our findings in the context of accretion disc dynamics around the vicinity of the BHs.


title: Self-Supervised Learning for 3D Medical Image Analysis using 3D SimCLR and Monte Carlo Dropout

abstract: Self-supervised learning methods can be used to learn meaningful representations from unlabeled data that can be transferred to supervised downstream tasks to reduce the need for labeled data. In this paper, we propose a 3D self-supervised method that is based on the contrastive (SimCLR) method. Additionally, we show that employing Bayesian neural networks (with Monte-Carlo Dropout) during the inference phase can further enhance the results on the downstream tasks. We showcase our models on two medical imaging segmentation tasks: i) Brain Tumor Segmentation from 3D MRI, ii) Pancreas Tumor Segmentation from 3D CT. Our experimental results demonstrate the benefits of our proposed methods in both downstream data-efficiency and performance.


title: Large-scale Vortices in Rapidly Rotating Rayleigh-B\'enard Convection at Small Prandtl number

abstract: One prominent feature in the atmospheres of Jupiter and Saturn is the appearance of large-scale vortices. However, the sustaining mechanism of these large-scale vortices remains unclear. One possible mechanism is that these large-scale vortices are driven by rotating convection. Here we present numerical simulation results on rapidly rotating Rayleigh-B\'enard convection at a small Prandtl number $Pr=0.1$ (close to the turbulent Prandtl numbers of Jupiter and Saturn). We have identified four flow regimes in our simulation: multiple small vortices, coexisted large-scale cyclone and anticyclone, large-scale cyclone, and turbulence. The formation of large-scale vortices requires two conditions to be satisfied: the vertical Reynolds number is large ($Re_{z}\ge 400$), and the Rossby number is small ($Ro\leq 0.4$). Large-scale cyclone first appears when $Ro$ decreases to be smaller than 0.4. When $Ro$ further decreases to be smaller than 0.1, coexisted large-scale anticyclone emerges. We have studied the heat transfer in rapidly rotating convection. The result reveals that the heat transfer is more efficient in the anticyclonic region than in the cyclonic region. Besides, we find that 2D effect increases and 3D effect decreases in transporting convective flux as rotation rate increases. We find that aspect ratio has an effect on the critical Rossby number for the emergence of large-scale vortices. Our results provide helpful insights on understanding the dynamics of large-scale vortices in gas giants.

categories: physics.flu-dyn astro-ph.EP astro-ph.SR physics.ao-ph physics.geo-ph


title: Residual-based adaptivity for two-phase flow simulation in porous media using Physics-informed Neural Networks

abstract: This paper aims to provide a machine learning framework to simulate two-phase flow in porous media. The proposed algorithm is based on Physics-informed neural networks (PINN). A novel residual-based adaptive PINN is developed and compared with the residual-based adaptive refinement (RAR) method and with PINN with fixed collocation points. The proposed algorithm is expected to have great potential to be applied to different fields where adaptivity is needed. In this paper, we focus on the two-phase flow in porous media problem. We provide two numerical examples to show the effectiveness of the new algorithm. It is found that adaptivity is essential to capture moving flow fronts. We show how the results obtained through this approach are more accurate than using RAR method or PINN with fixed collocation points, while having a comparable computational cost.


title: The existence of periodic solution and asymptotic behavior of solutions for a multi-layer tumor model with a periodic provision of external nutrients

abstract: In this paper, we consider a multi-layer tumor model with a periodic provision of external nutrients. The domain occupied by tumor has a different shape (flat shape) than spherical shape which has been studied widely. The important parameters are periodic external nutrients $\Phi(t)$ and threshold concentration for proliferation $\widetilde{\sigma}$. In this paper, we give a complete classification about $\Phi(t)$ and $\widetilde{\sigma}$ according to global stability of zero equilibrium solution or global stability of the positive periodic solution. Precisely, if $\frac{1}{T} \int_{0}^{T} \Phi(t)d t\leqslant\widetilde{\sigma}$, then the zero equilibrium solution is globally stable while if $\frac{1}{T} \int_{0}^{T} \Phi(t)d t>\widetilde{\sigma}$, then there exists a unique positive T-periodic solution and it is globally stable.


title: Magnetic interlayer coupling between ferromagnetic SrRuO$_3$ layers through a SrIrO$_3$ spacer

abstract: A key element to tailor the properties of magnetic multilayers is the coupling between the individual magnetic layers. In case of skyrmion hosting multilayers, coupling of skyrmions across the magnetic layers is highly desirable. Here the magnetic interlayer coupling was studied in epitaxial all-oxide heterostructures of ferromagnetic perovskite SrRuO$_3$ layers separated by spacers of the strong spin-orbit coupling oxide SrIrO$_3$. This combination of oxide layers is being discussed as a potential candidate system to host N\'{e}el skyrmions. First order reversal curve (FORC) measurements were performed in order to distinguish between magnetic switching processes of the individual layers and to disentangle the signal of soft magnetic impurities from the samples$'$ signal. Additionally, FORC investigations enabled to determine whether the coupling between the magnetic layers is ferromagnetic or antiferromagnetic. The observed interlayer coupling strength was weak for all the heterostructures, with SrIrO$_3$ spacers between 2 monolayers and 12 monolayers thick.


title: Dirac fermion optics and directed emission from single- and bilayer graphene cavities

abstract: High-mobility graphene hosting massless charge carriers with linear dispersion provides a promising platform for electron optics phenomena. Inspired by the physics of dielectric optical micro-cavities where the photon emission characteristics can be efficiently tuned via the cavity shape, we study corresponding mechanisms for trapped Dirac fermionic resonant states in deformed micro-disk graphene billiards and directed emission from those. In such graphene devices a back-gate voltage provides an additional tunable parameter to mimic different effective refractive indices and thereby the corresponding Fresnel laws at the boundaries. Moreover, cavities based on single-layer and double-layer graphene exhibit Klein- and anti-Klein tunneling, respectively, leading to distinct differences with respect to dwell times and resulting emission profiles of the cavity states. Moreover, we find a variety of different emission characteristics depending on the position of the source where charge carriers are fed into the cavites. Combining quantum mechanical simulations with optical ray tracing and a corresponding phase-space analysis, we demonstrate strong confinement of the emitted charge carriers in the mid field of single-layer graphene systems and can relate this to a lensing effect. For bilayer graphene, trapping of the resonant states is more efficient and the emission characteristics do less depend on the source position.


title: Persistent homology and the shape of evolutionary games

abstract: For nearly three decades, spatial games have produced a wealth of insights to the study of behavior and its relation to population structure. However, as different rules and factors are added or altered, the dynamics of spatial models often become increasingly complicated to interpret. To tackle this problem, we introduce persistent homology as a rigorous framework that can be used to both define and compute higher-order features of data in a manner which is invariant to parameter choices, robust to noise, and independent of human observation. Our work demonstrates its relevance for spatial games by showing how topological features of simulation data that persist over different spatial scales reflect the stability of strategies in 2D lattice games. To do so, we analyze the persistent homology of scenarios from two games: a Prisoner's Dilemma and a SIRS epidemic model. The experimental results show how the method accurately detects features that correspond to real aspects of the game dynamics. Unlike other tools that study dynamics of spatial systems, persistent homology can tell us something meaningful about population structure while remaining neutral about the underlying structure itself. Regardless of game complexity, since strategies either succeed or fail to conform to shapes of a certain topology there is much potential for the method to provide novel insights for a wide variety of spatially extended systems in biology, social science, and physics.


title: BEdgeHealth: A Decentralized Architecture for Edge-based IoMT Networks Using Blockchain

abstract: The healthcare industry has witnessed significant transformations in e-health services by using mobile edge computing (MEC) and blockchain to facilitate healthcare operations. Many MEC-blockchain-based schemes have been proposed, but some critical technical challenges still remain, such as low quality of services (QoS), data privacy and system security vulnerabilities. In this paper, we propose a new decentralized health architecture, called BEdgeHealth that integrates MEC and blockchain for data offloading and data sharing in distributed hospital networks. First, a data offloading scheme is proposed where mobile devices can offload health data to a nearby MEC server for efficient computation with privacy awareness. Moreover, we design a data sharing scheme which enables data exchanges among healthcare users by leveraging blockchain and interplanetary file system. Particularly, a smart contract-based authentication mechanism is integrated with MEC to perform decentralized user access verification at the network edge without requiring any central authority. The real-world experiment results and evaluations demonstrate the effectiveness of the proposed BEdgeHealth architecture in terms of improved QoS with data privacy and security guarantees, compared to the existing schemes.


title: Phase-dependent controllable field generation in a ring cavity resonator

abstract: We investigate the control field phase-dependent output field transmission from a red detuned ring cavity optomechanical system. Our scheme displays a double transparency window in the presence of a strong control and a weak probe field. Additionally, we invoke an external mechanical pump to one of the movable mirrors to modulate its vibration. Complete control over the output field transmission can be achieved due to the combined effect of the amplitudes and the phases of the mechanical pump and the control field. Further, a tunable group delay of the probe pulse propagation can be obtained by the tailoring of the control field phase in the presence of a suitable mechanical drive. We further discuss the effect of control field phase on Stokes field generation via the four-wave mixing process. This scheme may find potential applications in weak signal sensing and all-optical communication purposes.


title: First Evidence of Intrinsic Alignments of Red Galaxies at $z > 1$: Cross-correlation between CFHTLenS and FastSound Samples

abstract: We report the first evidence for intrinsic alignment (IA) of red galaxies at $z>1$. We measure the gravitational shear-intrinsic ellipticity (GI) cross-correlation function at $z\sim1.3$ using galaxy positions from the FastSound spectroscopic survey and galaxy shapes from Canada-Hawaii-France telescope lensing survey data. Adopting the non-linear alignment model, we obtain a $2.4\sigma$-level detection of the IA amplitude $A^{\rm LA}=27.48_{-11.54}^{+11.53}$ (and $2.6\sigma$ with $A^{\rm LA}=29.43_{-11.49}^{+11.48}$ when weak-lensing contaminations are taken into account), which is larger than the value extrapolated from the constraints obtained at lower redshifts. Our measured IA is translated into a $\sim 20\%$ contamination to the weak lensing power spectrum for the red galaxies. This marginal detection of IA for red galaxies at $z>1$ motivates the continuing investigation of the nature of IA for weak lensing studies. Furthermore, our result provides the first step to utilize IA measurements in future high-$z$ surveys as a cosmological probe, complementary to galaxy clustering and lensing.


title: Photoexcitation of atoms by cylindrically polarized Laguerre-Gaussian beams

abstract: We analyze the photo-excitation of atoms with a single valence electron by cylindrically polarized Laguerre-Gaussian beams. Theoretical analysis is performed within the framework of first-order perturbation theory and by expanding the vector potential of the Laguerre-Gaussian beam in terms of its multipole components. For cylindrically polarized Laguerre-Gaussian beams, we show that the (magnetic) sub-components of electric-quadrupole field vary significantly in the beam cross-section with beam waist and radial distance from the beam axis. We discuss the influence of varying magnetic multipole component in the beam cross-section on the sub-level population of a localized atomic target. In addition, we calculate the total excitation rate of electric quadrupole transition ($4s \;^{2}S_{1/2} \rightarrow 3d \;^{2}D_{5/2}$) in a mesoscopic target of Ca$^{+}$ ion. These calculations shows that the total rate of excitation is sensitive to the beam waist and the distance between center of the target and the beam axis. However, the excitation by cylindrically polarized Laguerre-Gaussian beam is found more efficient in driving electric quadrupole transition in the mesoscopic atomic target than the circularly polarized beams.


title: Greedy algorithms for learning via exponential-polynomial splines

abstract: Kernel-based schemes are state-of-the-art techniques for learning by data. In this work we extend some ideas about kernel-based greedy algorithms to exponential-polynomial splines, whose main drawback consists in possible overfitting and consequent oscillations of the approximant. To partially overcome this issue, we introduce two algorithms which perform an adaptive selection of the spline interpolation points based on the minimization either of the sample residuals ($f$-greedy), or of an upper bound for the approximation error based on the spline Lebesgue function ($\lambda$-greedy). Both methods allow us to obtain an adaptive selection of the sampling points, i.e. the spline nodes. However, while the {$f$-greedy} selection is tailored to one specific target function, the $\lambda$-greedy algorithm is independent of the function values and enables us to define a priori optimal interpolation nodes.


title: High-performance Ba1-xKxFe2As2 superconducting joints for persistent current operation

abstract: Superconducting joints are one of the key technologies to make Ba1-xKxFe2As2 (Ba-122) superconducting wires or tapes for high-field applications. Herein, superconducting joints were fabricated by a simple cold-pressing method, and the joint resistance of the iron-based superconducting joint was estimated for the first time. The superconducting properties, microstructures, and elements distribution in the joint regions were investigated. At 4.2 K and 10 T, a transport critical current Ic of 105 A for the joint was obtained, and the critical current ratio (CCR= Ic-joint/Ic-tape) of the joint was 94.6%. On the other hand, the joint show very low joint resistance of 2.7x10^-13 ohm in self-field at 4.2 K. Among iron-based superconductors (IBS), this work is the first to successfully realize a superconducting joint with such high CCR and low joint resistance. This work shows great potential to apply Ba-122 in a range of practical applications, where superconducting joints are essential.


title: Self-organized swimming with odd elasticity

abstract: We theoretically investigate self-oscillating waves of an active material, which have recently been introduced as a non-symmetric part of the elastic moduli, termed odd elasticity. Using Purcell's three-link swimmer model, we reveal that an odd-elastic filament at low Reynolds number can swim in a self-organized manner and that the time-periodic dynamics are characterized by a stable limit cycle generated by elastohydrodynamic interactions. Also, we consider a noisy shape gait and derive a swimming formula for a general elastic material in the Stokes regime with its elasticity modulus being represented by a non-symmetric matrix, demonstrating that the odd elasticity produces biased net locomotion from random noise.


title: Registry dependent potential for interfaces of gold with graphitic systems

abstract: We present a semi-anisotropic interfacial potential (SAIP) designed to classically describe the interaction between gold and two-dimensional (2D) carbon allotropes such as graphene, fullerenes, or hydrocarbon molecules. The potential is able to accurately reproduce dispersion corrected density functional theory (DFT+D3) calculations performed over selected configurations: a flat graphene sheet, a benzene molecule, and a C$_{60}$ fullerene, physisorbed on the Au(111) surface. The effects of bending and of hydrogen passivation on the potential terms are discussed. The presented SAIP potential provides a noticeable improvement in the state-of-the-art description of Au-C interfaces. Also, its functional form is suitable to describe the interfacial interaction between other 2D and bulk materials.


title: Towards event aggregation for reducing the volume of logged events during IKC stages of APT attacks

abstract: Nowadays, targeted attacks like Advanced Persistent Threats (APTs) has become one of the major concern of many enterprise networks. As a common approach to counter these attacks, security staff deploy a variety of security and non-security sensors at different lines of defense (Network, Host, and Application) to track the attacker's behaviors during their kill chain. However, one of the drawbacks of this approach is the huge amount of events raised by heterogeneous security and non-security sensors which makes it difficult to analyze logged events for later processing i.e. event correlation for timely detection of APT attacks. Till now, some research papers have been published on event aggregation for reducing the volume of logged low-level events. However, most research works have been provided a method to aggregate the events of a single-type and homogeneous event source i.e. NIDS. In addition, their main focus is only on the degree to which the event volume is reduced, while the amount of security information lost during the event aggregation process is also very important. In this paper, we propose a three-phase event aggregation method to reduce the volume of logged heterogeneous events during APT attacks considering the lowest rate of loss of security information. To this aim, at first, low-level events of the sensors are clustered into some similar event groups and then, after filtering noisy event clusters, the remained clusters are summarized based on an Attribute-Oriented Induction (AOI) method in a controllable manner to reduce the unimportant or duplicated events. The method has been evaluated on the three publicly available datasets: SotM34, Bryant, and LANL. The experimental results show that the method is efficient enough in event aggregation and can reduce events volume up to 99.7\% with an acceptable level of information loss ratio (ILR).


title: Double parton distributions out of bounds in colour space

abstract: We investigate the positivity of double parton distributions with a non-trivial dependence on the parton colour. It turns out that positivity is not preserved by leading-order evolution from lower to higher scales, in contrast to the case in which parton colour is summed over. We also study the positivity properties of the distributions at small distance between the two partons, where they can be computed in terms of perturbative splitting kernels and ordinary parton densities.


title: Algebras and Banach spaces of Dirichlet series with maximal Bohr's strip

abstract: We study linear and algebraic structures in sets of Dirichlet series with maximal Bohr's strip. More precisely, we consider a set $\mathscr M$ of Dirichlet series which are uniformly continuous on the right half plane and whose strip of uniform but not absolute convergence has maximal width, i.e., $\frac{1}{2}$. Considering the uniform norm, we show that $\mathscr M$ contains an isometric copy of $\ell_1$ (except zero) and is strongly $\aleph_0$-algebrable. Also, there is a dense $G_\delta$ set such that any of its elements generates a free algebra contained in $\mathscr M\cup \{0\}$. Furthermore, we investigate $\mathscr{M}$ as a subset of the Hilbert space of Dirichlet series whose coefficients are square-summable. In this case, we prove that $\mathscr M$ contains an isometric copy of $\ell_2$ (except zero).


title: Three-Stream 3D/1D CNN for Fine-Grained Action Classification and Segmentation in Table Tennis

abstract: This paper proposes a fusion method of modalities extracted from video through a three-stream network with spatio-temporal and temporal convolutions for fine-grained action classification in sport. It is applied to TTStroke-21 dataset which consists of untrimmed videos of table tennis games. The goal is to detect and classify table tennis strokes in the videos, the first step of a bigger scheme aiming at giving feedback to the players for improving their performance. The three modalities are raw RGB data, the computed optical flow and the estimated pose of the player. The network consists of three branches with attention blocks. Features are fused at the latest stage of the network using bilinear layers. Compared to previous approaches, the use of three modalities allows faster convergence and better performances on both tasks: classification of strokes with known temporal boundaries and joint segmentation and classification. The pose is also further investigated in order to offer richer feedback to the athletes.


title: More models of walks avoiding a quadrant (extended abstract)

abstract: We continue the enumeration of plane lattice paths avoiding the negative quadrant initiated by the first author in [Bousquet-M{\'e}lou, 2016]. We solve in detail a new case, the king walks, where all 8 nearest neighbour steps are allowed. As in the two cases solved in [Bousquet-M{\'e}lou, 2016], the associated generating function is proved to differ from a simple, explicit D-finite series (related to the enumeration of walks confined to the first quadrant) by an algebraic one. The principle of the approach is the same as in [Bousquet-M{\'e}lou, 2016], but challenging theoretical and computational difficulties arise as we now handle algebraic series of larger degree. We also explain why we expect the observed algebraicity phenomenon to persist for 4 more models, for which the quadrant problem is solvable using the reflection principle.


title: Graphene bilayer and trilayer Moir\'e lattice with Rashba spin-orbit coupling

abstract: We consider twisted bilayer and trilayer graphene in the presence of Rashba spin-orbit coupling and explore the physics of Moir\'e spintronics. The electronic charge density has a sharp step right at the magic angles $\theta_m$. As a result, local spin observables (polarization and equilibrium spin currents) have sharp peaks (of width about a small fraction of 1$^\circ$) as a function of the twist angle $\theta$, and abrupt sign reversals at $\theta_m$. Thereby, the magic angle can be determined in an unprecedented accuracy. In the first chiral limit, the spin currents vanish, but the peculiar pattern of the polarization at $\theta_m$ persists. Major differences result in spintronics of twisted bilayer graphene at magic angles as compared with the spintronics of single and/or {\it un-twisted} bilayer graphene. Thus, in addition to the numerous spectacular physical phenomena already reported in twisted bilayer graphene at magic angles, new phenomena also occur in twistronic spintronics.


title: Online Aggregation of Probability Forecasts with Confidence

abstract: The paper presents numerical experiments and some theoretical developments in prediction with expert advice (PEA). One experiment deals with predicting electricity consumption depending on temperature and uses real data. As the pattern of dependence can change with season and time of the day, the domain naturally admits PEA formulation with experts having different ``areas of expertise''. We consider the case where several competing methods produce online predictions in the form of probability distribution functions. The dissimilarity between a probability forecast and an outcome is measured by a loss function (scoring rule). A popular example of scoring rule for continuous outcomes is Continuous Ranked Probability Score (CRPS). In this paper the problem of combining probabilistic forecasts is considered in the PEA framework. We show that CRPS is a mixable loss function and then the time-independent upper bound for the regret of the Vovk aggregating algorithm using CRPS as a loss function can be obtained. Also, we incorporate a ``smooth'' version of the method of specialized experts in this scheme which allows us to combine the probabilistic predictions of the specialized experts with overlapping domains of their competence.


title: A noise reduction method for force measurements in water entry experiments based on the Ensemble Empirical Mode Decomposition

abstract: In this paper a denoising strategy based on the EEMD (Ensemble Empirical Mode Decomposition) is used to reduce the background noise in non-stationary signals, which represent the forces measured in scaled model testing of the emergency water landing of aircraft, generally referred to as ditching. Ditching tests are performed at a constant horizontal speed of 12 m/s and a vertical velocity at the beginning of the impact of 0.45 m/s. The measured data are affected by a large amplitude broadband noise, which has both mechanical and electronic origin. Noise sources cannot be easily avoided or removed, since they are associated with the vibrations of the structure of the towing carriage and to the interaction between the measurement chain and the electromagnetic fields. The EEMD noise reduction method is based on the decomposition of the signal into modes and on its partial reconstruction using the residue, the signal-dominant modes and some further modes treated with a thresholding technique, which helps to retain some of the sharp features of the signal. The strategy is developed and tested first on a synthetic signal with a superimposed and known background noise. The method is then verified on the measurement of the inertial force acting on the fuselage when it is moving in air, as in this case the added mass is negligible and the denoised force should equal the product of the mass by the acceleration, both of them being known. Finally, the procedure is applied to denoise the forces measured during the actual ditching experiments. The results are superior to those obtained by other classical filtering methods, such as a moving average filter and a low-pass FIR filter, particularly due to the enhanced capabilities of the EEMD-denoising strategy here developed to preserve the sharp features of the signals and to reduce the residual low-frequency oscillations of spurious origin.


title: Learning Dynamics Models for Model Predictive Agents

abstract: Model-Based Reinforcement Learning involves learning a \textit{dynamics model} from data, and then using this model to optimise behaviour, most often with an online \textit{planner}. Much of the recent research along these lines presents a particular set of design choices, involving problem definition, model learning and planning. Given the multiple contributions, it is difficult to evaluate the effects of each. This paper sets out to disambiguate the role of different design choices for learning dynamics models, by comparing their performance to planning with a ground-truth model -- the simulator. First, we collect a rich dataset from the training sequence of a model-free agent on 5 domains of the DeepMind Control Suite. Second, we train feed-forward dynamics models in a supervised fashion, and evaluate planner performance while varying and analysing different model design choices, including ensembling, stochasticity, multi-step training and timestep size. Besides the quantitative analysis, we describe a set of qualitative findings, rules of thumb, and future research directions for planning with learned dynamics models. Videos of the results are available at https://sites.google.com/view/learning-better-models.


title: GTC/CanariCam deep mid-infrared imaging survey of northern stars within 5 pc

abstract: In this work we present the results of a direct imaging survey for brown dwarf companions around the nearest stars at the mid-infrared 10 micron range ($\lambda_{c}$=8.7$\mu$m, $\Delta\lambda$=1.1$\mu$m) using the CanariCam instrument at the 10.4 m Gran Telescopio Canarias (GTC). We imaged the 25 nearest stellar systems within 5 pc of the Sun at declinations $\delta > -25^{\circ}$ (at least half have planets from radial velocity), reaching a mean detection limit of 11.3$\pm$0.2 mag (1.5 mJy) in the Si-2 8.7$\mu$m band over a range of angular separations from 1 to 10 arcsec. This would have allowed us to uncover substellar companions at projected orbital separations between $\sim$2 and 50 au, with effective temperatures down to 600 K and masses greater than 30 $M_{Jup}$ assuming an average age of 5 Gyr and down to the deuterium-burning mass limit for objects with ages $<$1 Gyr. From the non-detection of such companions, we determined upper limits on their occurrence rate at depths and orbital separations yet unexplored by deep imaging programs. For the M dwarfs, main components of our sample, we found with a 90% confidence level that less than 20% of these low-mass stars have L and T-type brown dwarf companions with $m \gtrsim 30 M_{Jup}$ and $T_{eff} \gtrsim$ 600 K at $\sim$3.5--35 au projected orbital separations.