title: Learning Cost Functions for Optimal Transport

abstract: Inverse optimal transport (OT) refers to the problem of learning the cost function for OT from observed transport plan or its samples. In this paper, we derive an unconstrained convex optimization formulation of the inverse OT problem, which can be further augmented by any customizable regularization. We provide a comprehensive characterization of the properties of inverse OT, including uniqueness of solutions. We also develop two numerical algorithms, one is a fast matrix scaling method based on the Sinkhorn-Knopp algorithm for discrete OT, and the other one is a learning based algorithm that parameterizes the cost function as a deep neural network for continuous OT. The novel framework proposed in the work avoids repeatedly solving a forward OT in each iteration which has been a thorny computational bottleneck for the bi-level optimization in existing inverse OT approaches. Numerical results demonstrate promising efficiency and accuracy advantages of the proposed algorithms over existing state-of-the-art methods.


title: Incomparable actions of free groups

abstract: Suppose that $X$ is a Polish space, $E$ is a countable Borel equivalence relation on $X$, and $\mu$ is an $E$-invariant Borel probability measure on $X$. We consider the circumstances under which for every countable non-abelian free group $\Gamma$, there is a Borel sequence $(\cdot_r)_{r \in \mathbb{R}}$ of free actions of $\Gamma$ on $X$, generating subequivalence relations $E_r$ of $E$ with respect to which $\mu$ is ergodic, with the further property that $(E_r)_{r \in \mathbb{R}}$ is an increasing sequence of relations which are pairwise incomparable under $\mu$-reducibility. In particular, we show that if $E$ satisfies a natural separability condition, then this is the case as long as there exists a free Borel action of a countable non-abelian free group on $X$, generating a subequivalence relation of $E$ with respect to which $\mu$ is ergodic.


title: Another determinantal inequality involving partial traces

abstract: Let $A$ be a positive semidefinite $m\times m$ block matrix with each block $n$-square, then the following determinantal inequality for partial traces holds \[ (\mathrm{tr} A)^{mn} - \det(\mathrm{tr}_2 A)^n \ge \bigl| \det A - \det(\mathrm{tr}_1 A)^m \bigr|, \] where $\mathrm{tr}_1$ and $\mathrm{tr}_2$ stand for the first and second partial trace, respectively. This result improves a recent result of Lin [14].


title: Measurable perfect matchings for acyclic locally countable Borel graphs

abstract: We characterize the structural impediments to the existence of Borel perfect matchings for acyclic locally countable Borel graphs admitting a Borel selection of finitely many ends from their connected components. In particular, this yields the existence of Borel matchings for such graphs of degree at least three. As a corollary, it follows that acyclic locally countable Borel graphs of degree at least three generating $\mu$-hyperfinite equivalence relations admit $\mu$-measurable matchings. We establish the analogous result for Baire measurable matchings in the locally finite case, and provide a counterexample in the locally countable case.


title: On-chip spectrometer using stratified waveguides filters

abstract: We present an ultra-compact single-shot spectrometer on silicon platform with broad operation bandwidth and high resolution. It consists of 32 stratified waveguide filters (SWFs) with diverse transmission spectra for sampling the unknown spectrum of the input signal and a specially designed ultra-compact structure for splitting the incident signal into 32 filters with low imbalance. Each SWF has a footprint less than 1um x 30um, while the 1x32 splitter and 32 filters in total occupy an area of about 35um x 260um, which to the best of our knowledge, is the smallest footprint spectrometer realized on silicon photonic platform. Experimental characteristics of the fabricated spectrometer demonstrate a broad operating bandwidth of 180nm centered at 1550nm and narrowband peaks with 0.45nm Full-Width-Half-Maximum (FWHM) can be clearly resolved. This concept can also be implemented using other material platforms for operation in optical spectral bands of interest for various applications.


title: Measure reducibility of countable Borel equivalence relations

abstract: We show that every basis for the countable Borel equivalence relations strictly above $\mathbb{E}_0$ under measure reducibility is uncountable, thereby ruling out natural generalizations of the Glimm-Effros dichotomy. We also push many known results concerning the abstract structure of the measure reducibility hierarchy to its base, using arguments substantially simpler than those previously employed.


title: A new hybrid approach for crude oil price forecasting: Evidence from multi-scale data

abstract: Faced with the growing research towards crude oil price fluctuations influential factors following the accelerated development of Internet technology, accessible data such as Google search volume index are increasingly quantified and incorporated into forecasting approaches. In this paper, we apply multi-scale data that including both GSVI data and traditional economic data related to crude oil price as independent variables and propose a new hybrid approach for monthly crude oil price forecasting. This hybrid approach, based on divide and conquer strategy, consists of K-means method, kernel principal component analysis and kernel extreme learning machine , where K-means method is adopted to divide input data into certain clusters, KPCA is applied to reduce dimension, and KELM is employed for final crude oil price forecasting. The empirical result can be analyzed from data and method levels. At the data level, GSVI data perform better than economic data in level forecasting accuracy but with opposite performance in directional forecasting accuracy because of Herd Behavior, while hybrid data combined their advantages and obtain best forecasting performance in both level and directional accuracy. At the method level, the approaches with K-means perform better than those without K-means, which demonstrates that divide and conquer strategy can effectively improve the forecasting performance.


title: Noncommutative Furstenberg boundary

abstract: We introduce and study the notions of boundary actions and of the Furstenberg boundary of a discrete quantum group. As for classical groups, properties of boundary actions turn out to encode significant properties of the operator algebras associated with the discrete quantum group in question; for example we prove that if the action on the Furstenberg boundary is faithful, the quantum group C*-algebra admits at most one KMS-state for the scaling automorphism group. To obtain these results we develop a version of Hamana's theory of injective envelopes for quantum group actions, and establish several facts on relative amenability for quantum subgroups. We then show that the Gromov boundary actions of free orthogonal quantum groups, as studied by Vaes and Vergnioux, are also boundary actions in our sense; we obtain this by proving that these actions admit unique stationary states. Moreover, we prove these actions are faithful, hence conclude a new unique KMS-state property in the general case, and a new proof of unique trace property when restricted to the unimodular case. We prove equivalence of simplicity of the crossed products of all boundary actions of a given discrete quantum group, and use it to obtain a new simplicity result for the crossed product of the Gromov boundary actions of free orthogonal quantum groups.


title: An Efficient MPC Algorithm For Switched Systems with Minimum Dwell Time Constraints

abstract: This paper presents an efficient suboptimal model predictive control (MPC) algorithm for nonlinear switched systems subject to minimum dwell time constraints (MTC). While MTC are required for most physical systems due to stability, power and mechanical restrictions, MPC optimization problems with MTC are challenging to solve. To efficiently solve such problems, the on-line MPC optimization problem is decomposed into a sequence of simpler problems, which include two nonlinear programs (NLP) and a rounding step, as typically done in mixed-integer optimal control (MIOC). Unlike the classical approach that embeds MTC in a mixed-integer linear program (MILP) with combinatorial constraints in the rounding step, our proposal is to embed the MTC in one of the NLPs using move blocking. Such a formulation can speedup on-line computations by employing recent move blocking algorithms for NLP problems and by using a simple sum-up-rounding (SUR) method for the rounding step. An explicit upper bound of the integer approximation error for the rounding step is given. In addition, a combined shrinking and receding horizon strategy is developed to satisfy closed-loop MTC. Recursive feasibility is proven using a $l$-step control invariant ($l$-CI) set, where $l$ is the minimum dwell time step length. An algorithm to compute $l$-CI sets for switched linear systems off-line is also presented. Numerical studies show significant speed-up and comparable control performance of the proposed MPC algorithm against the classical approach, though at the cost of sub-optimal solutions.


title: Minimal mass blow-up solutions to rough nonlinear Schroedinger equations

abstract: We study the focusing mass-critical rough nonlinear Schroedinger equations, where the stochastic integration is taken in the sense of controlled rough path. We obtain the global well-posedness if the mass of initial data is below that of the ground state. Moreover, the existence of minimal mass blow-up solutions is also obtained in both dimensions one and two. In particular, these yield that the mass of ground state is exactly the threshold of global well-posedness and blow-up of solutions in the stochastic focusing mass-critical case. Similar results are also obtained for a class of nonlinear Schroedinger equations with lower order perturbations.


title: New Lower Bounds for Tverberg Partitions with Tolerance in the Plane

abstract: Let $P$ be a set $n$ points in a $d$-dimensional space. Tverberg's theorem says that, if $n$ is at least $(k-1)(d+1)+1$, then $P$ can be partitioned into $k$ sets whose convex hulls intersect. Partitions with this property are called {\em Tverberg partitions}. A partition has tolerance $t$ if the partition remains a Tverberg partition after removal of any set of $t$ points from $P$. Tolerant Tverberg partitions exist in any dimension provided that $n$ is sufficiently large. Let $N(d,k,t)$ be the smallest value of $n$ such that tolerant Tverberg partitions exist for any set of $n$ points in $\mathbb{R}^d$. Only few exact values of $N(d,k,t)$ are known. In this paper we establish a new tight bound for $N(2,2,2)$. We also prove many new lower bounds on $N(2,k,t)$ for $k\ge 2$ and $t\ge 1$.


title: Multi-Branch Learning for Weakly-Labeled Sound Event Detection

abstract: There are two sub-tasks implied in the weakly-supervised SED: audio tagging and event boundary detection. Current methods which combine multi-task learning with SED requires annotations both for these two sub-tasks. Since there are only annotations for audio tagging available in weakly-supervised SED, we design multiple branches with different learning purposes instead of pursuing multiple tasks. Similar to multiple tasks, multiple different learning purposes can also prevent the common feature which the multiple branches share from overfitting to any one of the learning purposes. We design these multiple different learning purposes based on combinations of different MIL strategies and different pooling methods. Experiments on the DCASE 2018 Task 4 dataset and the URBAN-SED dataset both show that our method achieves competitive performance.


title: Theory of multiple quantum coherence signals in dilute thermal gases

abstract: Manifestations of dipole-dipole interactions in dilute thermal gases are difficult to sense because of strong inhomogeneous broadening. Recent experiments reported signatures of such interactions in fluorescence detection-based measurements of multiple quantum coherence (MQC) signals, with many characteristic features hitherto unexplained. We develop an original open quantum systems theory of MQC in dilute thermal gases, which allows us to resolve this conundrum. Our theory accounts for the vector character of the atomic dipoles as well as for driving laser pulses of arbitrary strength, includes the far-field coupling between the dipoles, which prevails in dilute ensembles, and effectively incorporates atomic motion via a disorder average. We show that collective decay processes -- which were ignored in previous treatments employing the electrostatic form of dipolar interactions -- play a key role in the emergence of MQC signals.


title: Active Lighting Recurrence by Parallel Lighting Analogy for Fine-Grained Change Detection

abstract: This paper studies a new problem, namely active lighting recurrence (ALR) that physically relocalizes a light source to reproduce the lighting condition from single reference image for a same scene, which may suffer from fine-grained changes during twice observations. ALR is of great importance for fine-grained visual inspection and change detection, because some phenomena or minute changes can only be clearly observed under particular lighting conditions. Therefore, effective ALR should be able to online navigate a light source toward the target pose, which is challenging due to the complexity and diversity of real-world lighting and imaging processes. To this end, we propose to use the simple parallel lighting as an analogy model and based on Lambertian law to compose an instant navigation ball for this purpose. We theoretically prove the feasibility, i.e., equivalence and convergence, of this ALR approach for realistic near point light source and small near surface light source. Besides, we also theoretically prove the invariance of our ALR approach to the ambiguity of normal and lighting decomposition. The effectiveness and superiority of the proposed approach have been verified by both extensive quantitative experiments and challenging real-world tasks on fine-grained change detection of cultural heritages. We also validate the generality of our approach to non-Lambertian scenes.


title: Book-Ahead & Supply Management for Ridesourcing Platforms

abstract: Ridesourcing platforms recently introduced the ``schedule a ride'' service where passengers may reserve (book-ahead) a ride in advance of their trip. Reservations give platforms precise information that describes the start time and location of anticipated future trips; in turn, platforms can use this information to adjust the availability and spatial distribution of the driver supply. In this article, we propose a framework for modeling/analyzing reservations in time-varying stochastic ridesourcing systems. We consider that the driver supply is distributed over a network of geographic regions and that book-ahead rides have reach time priority over non-reserved rides. First, we propose a state-dependent admission control policy that assigns drivers to passengers; this policy ensures that the reach time service requirement would be attained for book-ahead rides. Second, given the admission control policy and reservations information in each region, we predict the ``target" number of drivers that is required (in the future) to probabilistically guarantee the reach time service requirement for stochastic non-reserved rides. Third, we propose a reactive dispatching/rebalancing mechanism that determines the adjustments to the driver supply that are needed to maintain the targets across regions. For a specific reach time quality of service, simulation results using data from Lyft rides in Manhattan exhibit how the number of idle drivers decreases with the fraction of book-ahead rides. We also observe that the non-stationary demand (ride request) rate varies significantly across time; this rapid variation further illustrates that time-dependent models are needed for operational analysis of ridesourcing systems.


title: Relative big polynomial rings

abstract: Let $K$ be the field of Laurent series with complex coefficients, let $\mathcal{R}$ be the inverse limit of the standard-graded polynomial rings $K[x_1, \ldots, x_n]$, and let $\mathcal{R}^{\flat}$ be the subring of $\mathcal{R}$ consisting of elements with bounded denominators. In previous joint work with Erman and Sam, we showed that $\mathcal{R}$ and $\mathcal{R}^{\flat}$ (and many similarly defined rings) are abstractly polynomial rings, and used this to give new proofs of Stillman's conjecture. In this paper, we prove the complementary result that $\mathcal{R}$ is a polynomial algebra over $\mathcal{R}^{\flat}$.


title: String stable integral control design for vehicle platoons with disturbances

abstract: This paper presents a control design with integral action for vehicle platoons with disturbance that ensures string stability of the closed loop and disturbance rejection. The addition of integral action and a coordinate change allows to develop sufficient smoothness conditions on the closed loop system to ensure that the closed loop system using the proposed controller is string stable in the presence of time-varying disturbances and able able to reject constant disturbances. In addition, bounds for the tracking error of the platoon configuration are also given. Further, a case study is considered together with a suitable controller structure, which satisfies the required smoothness conditions. Simulation results illustrate the performance of the closed loop.


title: Global well-posedness to three-dimensional full compressible magnetohydrodynamic equations with vacuum

abstract: This paper studies the Cauchy problem for three-dimensional viscous, compressible, and heat conducting magnetohydrodynamic equations with vacuum as far field density. We prove the global existence and uniqueness of strong solutions provided that the quantity $\|\rho_0\|_{L^\infty}+\|b_0\|_{L^3}$ is suitably small and the viscosity coefficients satisfy $3\mu>\lambda$. Here, the initial velocity and initial temperature could be large. The assumption on the initial density do not exclude that the initial density may vanish in a subset of $\mathbb{R}^3$ and that it can be of a nontrivially compact support. Our result is an extension of the works of Fan and Yu \cite{FY09} and Li et al. \cite{LXZ13}, where the local strong solutions in three dimensions and the global strong solutions for isentropic case were obtained, respectively. The analysis is based on some new mathematical techniques and some new useful energy estimates. This paper can be viewed as the first result concerning the global existence of strong solutions with vacuum at infinity in some classes of large data in higher dimension.


title: Communication-Efficient Edge AI: Algorithms and Systems

abstract: Artificial intelligence (AI) has achieved remarkable breakthroughs in a wide range of fields, ranging from speech processing, image classification to drug discovery. This is driven by the explosive growth of data, advances in machine learning (especially deep learning), and easy access to vastly powerful computing resources. Particularly, the wide scale deployment of edge devices (e.g., IoT devices) generates an unprecedented scale of data, which provides the opportunity to derive accurate models and develop various intelligent applications at the network edge. However, such enormous data cannot all be sent from end devices to the cloud for processing, due to the varying channel quality, traffic congestion and/or privacy concerns. By pushing inference and training processes of AI models to edge nodes, edge AI has emerged as a promising alternative. AI at the edge requires close cooperation among edge devices, such as smart phones and smart vehicles, and edge servers at the wireless access points and base stations, which however result in heavy communication overheads. In this paper, we present a comprehensive survey of the recent developments in various techniques for overcoming these communication challenges. Specifically, we first identify key communication challenges in edge AI systems. We then introduce communication-efficient techniques, from both algorithmic and system perspectives for training and inference tasks at the network edge. Potential future research directions are also highlighted.


title: Revealing Rejuvenated Disorder States towards Crystallization in a Supercooled Metallic Glass-Forming Liquid

abstract: We report a metadynamics simulation study of crystallization in a deep undercooled metallic glass-forming liquid by developing appropriate collective variables. Through a combined analysis of free energy surface (FES) and atomic-level behaviors, a picture of an abnormal-endothermic crystallization process is revealed: rejuvenated disorder states with less local fivefold-symmetry and fast dynamics form firstly by changing the local chemical order around Cu atoms, which then act as the precursor for the nucleation of well-ordered crystallites. This process reflects a complex energy landscape with well-separated glassy and crystal basins, giving rise to the direct evidence of intrinsic frustration against crystallization in deep undercooled metallic glass forming liquids. Moreover, the rejuvenated disorder states with distinct physical behaviors offer great opportunities to tailor the performances of metallic glass by controlling the thermal history of a metallic melt.


title: Nonmyopic Gaussian Process Optimization with Macro-Actions

abstract: This paper presents a multi-staged approach to nonmyopic adaptive Gaussian process optimization (GPO) for Bayesian optimization (BO) of unknown, highly complex objective functions that, in contrast to existing nonmyopic adaptive BO algorithms, exploits the notion of macro-actions for scaling up to a further lookahead to match up to a larger available budget. To achieve this, we generalize GP upper confidence bound to a new acquisition function defined w.r.t. a nonmyopic adaptive macro-action policy, which is intractable to be optimized exactly due to an uncountable set of candidate outputs. The contribution of our work here is thus to derive a nonmyopic adaptive epsilon-Bayes-optimal macro-action GPO (epsilon-Macro-GPO) policy. To perform nonmyopic adaptive BO in real time, we then propose an asymptotically optimal anytime variant of our epsilon-Macro-GPO policy with a performance guarantee. We empirically evaluate the performance of our epsilon-Macro-GPO policy and its anytime variant in BO with synthetic and real-world datasets.


title: Vehicle Tracking in Wireless Sensor Networks via Deep Reinforcement Learning

abstract: Vehicle tracking has become one of the key applications of wireless sensor networks (WSNs) in the fields of rescue, surveillance, traffic monitoring, etc. However, the increased tracking accuracy requires more energy consumption. In this letter, a decentralized vehicle tracking strategy is conceived for improving both tracking accuracy and energy saving, which is based on adjusting the intersection area between the fixed sensing area and the dynamic activation area. Then, two deep reinforcement learning (DRL) aided solutions are proposed relying on the dynamic selection of the activation area radius. Finally, simulation results show the superiority of our DRL aided design.


title: Aggregation dynamics of methacrylate binary and ternary biomimetic polymers in solution

abstract: Using detailed atomistic simulations,we explore the conformational landscape of aggregates formed by biomimetic antimicrobial(AM) binary methacrylate copolymers,with hydrophobic and charged functional groups and the role of inclusion of polar functional groups on such aggregate morphologies.The effect of sequence of the constituent functional groups on aggregate conformation is also studied by considering random and block sequences along the polymer backbone.Our results suggest that block binary copolymers form large spherical aggregates with effective shielding of hydrophobic groups by charged groups.In contrast, random binary copolymers tend to form more bundle-like structures with exposed hydrophobic groups.The strong aggregation of binary polymers is driven primarily by attractive interactions between hydrophobic groups.However,replacing some of the hydrophobic groups with overall charge neutral polar groups weakens the aggregate considerably, leading to increased conformational fluctuations and formation of loose-packed,open aggregates,particularly in case of random ternary polymers.Interaction energy calculations strongly suggest that the role of inclusion of polar groups is two-fold: (1)to reduce possible strong local concentration of hydrophobic groups and smear the hydrophobicity along the polymer backbone to increase the solubility of the polymers (2)to compensate the loss of attractive hydrophobic interactions by forming attractive electrostatic interactions with charged groups and contribute to aggregation formation,albeit weak.Given that most of the naturally occurring AM peptides have contributions from all the three functional groups,this study elucidates the functionally tuneable role of inclusion of polar groups in the way AM agents interact with each other in solution phase,which can eventually dictate their partitioning behavior into bacterial and mammalian membranes.


title: Incorporating Effective Global Information via Adaptive Gate Attention for Text Classification

abstract: The dominant text classification studies focus on training classifiers using textual instances only or introducing external knowledge (e.g., hand-craft features and domain expert knowledge). In contrast, some corpus-level statistical features, like word frequency and distribution, are not well exploited. Our work shows that such simple statistical information can enhance classification performance both efficiently and significantly compared with several baseline models. In this paper, we propose a classifier with gate mechanism named Adaptive Gate Attention model with Global Information (AGA+GI), in which the adaptive gate mechanism incorporates global statistical features into latent semantic features and the attention layer captures dependency relationship within the sentence. To alleviate the overfitting issue, we propose a novel Leaky Dropout mechanism to improve generalization ability and performance stability. Our experiments show that the proposed method can achieve better accuracy than CNN-based and RNN-based approaches without global information on several benchmarks.


title: Temporal Sparse Adversarial Attack on Sequence-based Gait Recognition

abstract: Gait recognition is widely used in social security applications due to its advantages in long-distance human identification. Recently, sequence-based methods have achieved high accuracy by learning abundant temporal and spatial information. However, their robustness under adversarial attacks has not been clearly explored. In this paper, we demonstrate that the state-of-the-art gait recognition model is vulnerable to such attacks. To this end, we propose a novel temporal sparse adversarial attack method. Different from previous additive noise models which add perturbations on original samples, we employ a generative adversarial network based architecture to semantically generate adversarial high-quality gait silhouettes or video frames. Moreover, by sparsely substituting or inserting a few adversarial gait silhouettes, the proposed method ensures its imperceptibility and achieves a high attack success rate. The experimental results show that if only one-fortieth of the frames are attacked, the accuracy of the target model drops dramatically.


title: Hadronic decays of the $J/\psi$ meson

abstract: In this work we propose mainly phenomenological and theoretical models that allow, together with the use of the available experimental data, the calculation of form factors (FFs) and decay amplitudes. We focus our attention firstly on the $J/\psi$ decay into a pair of pions, showing that this process does not proceed only electromagnetically as believed so far, due to the presence of a non-negligible mixed strong-EM contribution to the total branching ratio. Moreover we consider the decay of the $J/\psi$ meson into a pair of baryon-antibaryon, where we separate, for the first time, the single strong, electromagnetic and strong-EM contributions to the total BR and the relative Feynman amplitudes, obtaining also the relative phase between the strong and the electromagnetic ones.


title: Guided Constrained Policy Optimization for Dynamic Quadrupedal Robot Locomotion

abstract: Deep reinforcement learning (RL) uses model-free techniques to optimize task-specific control policies. Despite having emerged as a promising approach for complex problems, RL is still hard to use reliably for real-world applications. Apart from challenges such as precise reward function tuning, inaccurate sensing and actuation, and non-deterministic response, existing RL methods do not guarantee behavior within required safety constraints that are crucial for real robot scenarios. In this regard, we introduce guided constrained policy optimization (GCPO), an RL framework based upon our implementation of constrained proximal policy optimization (CPPO) for tracking base velocity commands while following the defined constraints. We also introduce schemes which encourage state recovery into constrained regions in case of constraint violations. We present experimental results of our training method and test it on the real ANYmal quadruped robot. We compare our approach against the unconstrained RL method and show that guided constrained RL offers faster convergence close to the desired optimum resulting in an optimal, yet physically feasible, robotic control behavior without the need for precise reward function tuning.


title: Kernel interpolation with continuous volume sampling

abstract: A fundamental task in kernel methods is to pick nodes and weights, so as to approximate a given function from an RKHS by the weighted sum of kernel translates located at the nodes. This is the crux of kernel density estimation, kernel quadrature, or interpolation from discrete samples. Furthermore, RKHSs offer a convenient mathematical and computational framework. We introduce and analyse continuous volume sampling (VS), the continuous counterpart -- for choosing node locations -- of a discrete distribution introduced in (Deshpande & Vempala, 2006). Our contribution is theoretical: we prove almost optimal bounds for interpolation and quadrature under VS. While similar bounds already exist for some specific RKHSs using ad-hoc node constructions, VS offers bounds that apply to any Mercer kernel and depend on the spectrum of the associated integration operator. We emphasize that, unlike previous randomized approaches that rely on regularized leverage scores or determinantal point processes, evaluating the pdf of VS only requires pointwise evaluations of the kernel. VS is thus naturally amenable to MCMC samplers.


title: APRIL : a novel Algorithm for Particle Reconstruction at ILC

abstract: The current developments for future electron-positron colliders are driven by the Particle Flow concept. In these developments, high granularity calorimeters play a central role. This presentation will focus on a new Particle Flow Algorithm (PFA) developed for high granularity calorimeters, and especially for the Semi-Digital Hadronic CALorimeter (SDHCAL) option of the International Large Detector (ILD) project. The first PFA for ILD was PandoraPFA. This new PFA (APRIL) is based on the PandoraPFA Software Development Kit, but implements a different clustering inspired from the ARBOR PFA approach. This proceeding will describe briefly the APRIL algorithm and discuss its performance against that of PandoraPFA.


title: On the mean value property of fractional harmonic functions

abstract: As well known, harmonic functions satisfy the mean value property, namely the average of the function over a ball is equal to its value at the center. This fact naturally raises the question on whether this is a characterizing feature of balls, namely whether a set for which all harmonic functions satisfy the mean value property is necessarily a ball. This question was investigated by several authors, and was finally elegantly, completely and positively settled by \"Ulk\"u Kuran, with an artful use of elementary techniques. This classical problem has been recently fleshed out by Giovanni Cupini, Nicola Fusco, Ermanno Lanconelli and Xiao Zhong who proved a quantitative stability result for the mean value formula, showing that a suitable "mean value gap" (measuring the normalized difference between the average of harmonic functions on a given set and their pointwise value) is bounded from below by the Lebesgue measure of the "gap" between the set and the ball (and, consequently, by the Fraenkel asymmetry of the set). That is, if a domain "almost" satisfies the mean value property, then it must be necessarily close to a ball. Here we investigate the nonlocal counterparts of these results. In particular we will prove a classification result and a stability result, establishing that: if fractional harmonic functions enjoy a suitable exterior average property for a given domain, then the domain is necessarily a ball, a suitable "nonlocal mean value gap" is bounded from below by an appropriate measure of the difference between the set and the ball. Differently from the classical case, some of our arguments rely on purely nonlocal properties, with no classical counterpart, such as the fact that "all functions are locally fractional harmonic up to a small error".


title: On Boolean gates in fungal colony

abstract: A fungal colony maintains its integrity via flow of cytoplasm along mycelium network. This flow, together with possible coordination of mycelium tips propagation, is controlled by calcium waves and associated waves of electrical potential changes. We propose that these excitation waves can be employed to implement a computation in the mycelium networks. We use FitzHugh-Nagumo model to imitate propagation of excitation in a single colony of Aspergillus niger. Boolean values are encoded by spikes of extracellular potential. We represent binary inputs by electrical impulses on a pair of selected electrodes and we record responses of the colony from sixteen electrodes. We derive sets of two-inputs-on-output logical gates implementable the fungal colony and analyse distributions of the gates.


title: Towards field-programmable photonic gate arrays

abstract: We review some of the basic principles, fundamentals, technologies, architectures and recent advances leading to thefor the implementation of Field Programmable Photonic Field Arrays (FPPGAs).


title: Concurrent Kleene Algebra with Observations: from Hypotheses to Completeness

abstract: Concurrent Kleene Algebra (CKA) extends basic Kleene algebra with a parallel composition operator, which enables reasoning about concurrent programs. However, CKA fundamentally misses tests, which are needed to model standard programming constructs such as conditionals and $\mathsf{while}$-loops. It turns out that integrating tests in CKA is subtle, due to their interaction with parallelism. In this paper we provide a solution in the form of Concurrent Kleene Algebra with Observations (CKAO). Our main contribution is a completeness theorem for CKAO. Our result resorts on a more general study of CKA "with hypotheses", of which CKAO turns out to be an instance: this analysis is of independent interest, as it can be applied to extensions of CKA other than CKAO.


title: Incoherent $\rho$ meson photoproduction in ultraperipheral nuclear collisions at the LHC

abstract: Using the Gribov-Glauber model for photon-nucleus scattering and a generalization of the vector meson dominance model for the hadronic structure of the photon, we make predictions for the cross section of incoherent $\rho$ photoproduction in Pb-Pb ultraperipheral collisions (UPCs) in the Large Hadron Collider kinematics. We find that the effect of the inelastic nuclear shadowing is significant and leads to an additional 25% suppression of the incoherent cross section. Comparing our predictions to those of the STARlight Monte Carlo framework, we observe very significant differences.


title: Time-Varying Internal Models and Regulation of Unknown Harmonic Signals for Regular Linear Systems

abstract: We introduce general results on well-posedness and output regulation of regular linear systems with nonautonomous controllers. We present a generalization of the internal model principle for time-dependent controllers with asymptotically converging parameters. This general result is utilised in controller design for output tracking and disturbance rejection of harmonic signals with unknown frequencies. Our controller can be flexibly combined with different frequency estimation methods. The results are illustrated in rejection of unknown harmonic disturbances for a one-dimensional boundary controlled heat equation.


title: Investigating Typed Syntactic Dependencies for Targeted Sentiment Classification Using Graph Attention Neural Network

abstract: Targeted sentiment classification predicts the sentiment polarity on given target mentions in input texts. Dominant methods employ neural networks for encoding the input sentence and extracting relations between target mentions and their contexts. Recently, graph neural network has been investigated for integrating dependency syntax for the task, achieving the state-of-the-art results. However, existing methods do not consider dependency label information, which can be intuitively useful. To solve the problem, we investigate a novel relational graph attention network that integrates typed syntactic dependency information. Results on standard benchmarks show that our method can effectively leverage label information for improving targeted sentiment classification performances. Our final model significantly outperforms state-of-the-art syntax-based approaches.


title: Dark soliton of polariton condensates under nonresonant PT-symmetric pumping

abstract: A quantum system in complex potentials obeying parity-time (PT ) symmetry could exhibit all real spectra, starting out in non-Hermitian quantum mechanics. The key physics behind a PT-symmetric system consists of the balanced gain and loss of the complex potential. In this work, we plan to include the nonequilibrium nature, i.e. the intrinsic kinds of gain and loss of a system, to a PT-symmetric many-body quantum system, with the emphasis on the combined effects of non-Hermitian due to nonequilibrium nature and PT symmetry in determining the properties of a system. In this end, we investigate the static and dynamical properties of a dark soliton of a polariton Bose-Einstein condensate under the PT-symmetric non-resonant pumping by solving the driven-dissipative GrossPitaevskii equation both analytically and numerically. We derive the equation of motion for the center of mass of the dark solitons center analytically with the help of the Hamiltonian approach. The resulting equation captures how the combination of the open-dissipative character and PT-symmetry affects the properties of dark soliton, i.e. the soliton relaxes by blending with the background at a finite time. Further numerical solutions are in excellent agreement with the analytical results.


title: Multiple orthogonal geodesic chords and a proof of Seifert's conjecture on brake orbits

abstract: Using nonsmooth critical point theory, we prove the existence of at least N orthogonal geodesic chords in a class of Riemannian N-disk with strongly concave boundary. This yields a proof of a celebrated conjecture by H.Seifert on the number of brake orbits in a potential well of a natural Lagrangian/Hamiltonian system.


title: Proof-of-Concept of Uncompressed 4K Video Transmission from Drone through mmWave

abstract: Drones are attracting increasing attention in varieties of research fields because of their flexibility and are expected to be applied to a wide range of potential applications, among which the super-high-resolution video surveillance system using drones especially gains the authors research attention. Surveillance systems using cameras with fixed locations always suffer the blind spots due to the blockage or inappropriate deployments. Instead, by using the drones equipped with cameras, the surveillance performance can be drastically improved due to their high mobilities. The video quality is also a key factor of the surveillance performance. In face recognition, one of the most important surveillance applications, the uncompressed video can greatly improve the detection accuracy, but it is difficult to transmit uncompressed video in real time due to the huge data sizes. To address the issue, we propose to use the ultra-high speed mmWave communication for the video transmission from drones. Moreover, due to the limited battery energy and computing power in drones, we introduce the edge computing and propose to offload all the computation from the drones to the ground station. In addition, a proof-of-concept prototype hardware of the proposed uncompressed 4K video transmission system from drones through mmWave is developed, and the experiments results are consistent with the system design expectations.


title: Fair and Decentralized Exchange of Digital Goods

abstract: We construct a privacy-preserving, distributed and decentralized marketplace where parties can exchange data for tokens. In this market, buyers and sellers make transactions in a blockchain and interact with a third party, called notary, who has the ability to vouch for the authenticity and integrity of the data. We introduce a protocol for the data-token exchange where neither party gains more information than what it is paying for, and the exchange is fair: either both parties gets the other's item or neither does. No third party involvement is required after setup, and no dispute resolution is needed.


title: A Positive and Energy Stable Numerical Scheme for the Poisson-Nernst-Planck-Cahn-Hilliard Equations with Steric Interactions

abstract: We consider numerical methods for the Poisson-Nernst-Planck-Cahn-Hilliard (PNPCH) equations with steric interactions. We propose a novel energy stable numerical scheme that respects mass conservation and positivity at the discrete level. Existence and uniqueness of the solution to the proposed nonlinear scheme are established by showing that the solution is a unique minimizer of a convex functional over a closed, convex domain. The positivity of numerical solutions is further theoretically justified by the singularity of the entropy terms, which prevents the minimizer from approaching zero concentrations. A further numerical analysis proves discrete free-energy dissipation. Extensive numerical tests are performed to validate that the numerical scheme is first-order accurate in time and second-order accurate in space, and is capable of preserving the desired properties, such as mass conservation, positivity, and free energy dissipation, at the discrete level. Moreover, the PNPCH equations and the proposed scheme are applied to study charge dynamics and self-assembled nanopatterns in highly concentrated electrolytes that are widely used in electrochemical energy devices. Numerical results demonstrate that the PNPCH equations and our numerical scheme are able to capture nanostructures, such as lamellar patterns and labyrinthine patterns in electric double layers and the bulk, and multiple time relaxation with multiple time scales. In addition, we numerically characterize the interplay between cross steric interactions of short range and the concentration gradient regularization, and their impact on the development of nanostructures in the equilibrium state.

categories: math.NA cs.NA physics.chem-ph physics.comp-ph


title: The Soft State of the Black Hole Transient Source MAXI J1820+070: Emission from the Edge of the Plunge Region?

abstract: The Galactic black hole X-ray binary MAXI J1820+070 had a bright outburst in 2018 when it became the second brightest X-ray source in the Sky. It was too bright for X-ray CCD instruments such as XMM-Newton and Chandra, but was well observed by photon-counting instruments such as NICER and NuSTAR. We report here on the discovery of an excess emission component during the soft state. It is best modelled with a blackbody spectrum in addition to the regular disk emission, modelled either as diskbb or kerrbb. Its temperature varies from about 0.9 to 1.1 keV which is about 30 to 80 per cent higher than the inner disc temperature of diskbb. Its flux varies between 4 and 12 percent of the disc flux. Simulations of magnetised accretion discs have predicted the possibility of excess emission associated with a non-zero torque at the Innermost Stable Circular Orbit (ISCO) about the black hole, which from other NuSTAR studies lies at about 5 gravitational radii or about 60 km (for a black hole mass is 8 M). In this case the emitting region at the ISCO has a width varying between 1.3 and 4.6 km and would encompass the start of the plunge region where matter begins to fall freely into the black hole.


title: Communication-Efficient Decentralized Learning with Sparsification and Adaptive Peer Selection

abstract: Distributed learning techniques such as federated learning have enabled multiple workers to train machine learning models together to reduce the overall training time. However, current distributed training algorithms (centralized or decentralized) suffer from the communication bottleneck on multiple low-bandwidth workers (also on the server under the centralized architecture). Although decentralized algorithms generally have lower communication complexity than the centralized counterpart, they still suffer from the communication bottleneck for workers with low network bandwidth. To deal with the communication problem while being able to preserve the convergence performance, we introduce a novel decentralized training algorithm with the following key features: 1) It does not require a parameter server to maintain the model during training, which avoids the communication pressure on any single peer. 2) Each worker only needs to communicate with a single peer at each communication round with a highly compressed model, which can significantly reduce the communication traffic on the worker. We theoretically prove that our sparsification algorithm still preserves convergence properties. 3) Each worker dynamically selects its peer at different communication rounds to better utilize the bandwidth resources. We conduct experiments with convolutional neural networks on 32 workers to verify the effectiveness of our proposed algorithm compared to seven existing methods. Experimental results show that our algorithm significantly reduces the communication traffic and generally select relatively high bandwidth peers.


title: Interpretable Crowd Flow Prediction with Spatial-Temporal Self-Attention

abstract: Crowd flow prediction has been increasingly investigated in intelligent urban computing field as a fundamental component of urban management system. The most challenging part of predicting crowd flow is to measure the complicated spatial-temporal dependencies. A prevalent solution employed in current methods is to divide and conquer the spatial and temporal information by various architectures (e.g., CNN/GCN, LSTM). However, this strategy has two disadvantages: (1) the sophisticated dependencies are also divided and therefore partially isolated; (2) the spatial-temporal features are transformed into latent representations when passing through different architectures, making it hard to interpret the predicted crowd flow. To address these issues, we propose a Spatial-Temporal Self-Attention Network (STSAN) with an ST encoding gate that calculates the entire spatial-temporal representation with positional and time encodings and therefore avoids dividing the dependencies. Furthermore, we develop a Multi-aspect attention mechanism that applies scaled dot-product attention over spatial-temporal information and measures the attention weights that explicitly indicate the dependencies. Experimental results on traffic and mobile data demonstrate that the proposed method reduces inflow and outflow RMSE by 16% and 8% on the Taxi-NYC dataset compared to the SOTA baselines.


title: A New Family of Boundary-Domain Integral Equations for the Dirichlet Problem of the Diffusion Equation in Inhomogeneous Media with $H^{-1}(\Omega)$ Source Term on Lipschitz Domains

abstract: The interior Dirichlet boundary value problem for the diffusion equation in non-homogeneous media is reduced to a system of Boundary-Domain Integral Equations (BDIEs) employing the parametrix obtained in (Fresneda-Portillo, 2019) different from (Chkadua et. al 2009). We further extend the results obtained in (Fresneda-Portillo, 2019) for the mixed problem in a smooth domain with $L^{2}(\Omega)$ right hand side to Lipschitz domains and source term $f$ in the Sobolev space $H^{-1}(\Omega)$, where neither the classical nor the canonical co-normal derivatives are well defined. Equivalence between the system of BDIEs and the original BVP is proved along with their solvability and solution uniqueness in appropriate Sobolev spaces.


title: A New Unified Deep Learning Approach with Decomposition-Reconstruction-Ensemble Framework for Time Series Forecasting

abstract: A new variational mode decomposition (VMD) based deep learning approach is proposed in this paper for time series forecasting problem. Firstly, VMD is adopted to decompose the original time series into several sub-signals. Then, a convolutional neural network (CNN) is applied to learn the reconstruction patterns on the decomposed sub-signals to obtain several reconstructed sub-signals. Finally, a long short term memory (LSTM) network is employed to forecast the time series with the decomposed sub-signals and the reconstructed sub-signals as inputs. The proposed VMD-CNN-LSTM approach is originated from the decomposition-reconstruction-ensemble framework, and innovated by embedding the reconstruction, single forecasting, and ensemble steps in a unified deep learning approach. To verify the forecasting performance of the proposed approach, four typical time series datasets are introduced for empirical analysis. The empirical results demonstrate that the proposed approach outperforms consistently the benchmark approaches in terms of forecasting accuracy, and also indicate that the reconstructed sub-signals obtained by CNN is of importance for further improving the forecasting performance.


title: Molecular dynamics and continuum analyses of the electrokinetic zeta potential in nanostructured slit channels

abstract: This work presents a theoretical and numerical study of electrokinetic flow and the zeta potential for the case of slit channels with nanoscale roughness of dimensions comparable to the Debye length, by employing molecular dynamics simulations and continuum-level analyses. A simple analytical model for considering the effects of the surface roughness is proposed by employing matched asymptotic solutions for the charge density and fluid flow field, and matching conditions that satisfy electroneutrality and the Onsager reciprocal relation between the electroosmotic flow rate and streaming current. The proposed analytical model quantitatively accounts for results from molecular dynamics simulations that consider the presence of ion solvation shells and surface hydration layers. Our analysis indicates that a simultaneous knowledge of the electroosmotic and pressure-driven flow rate or streaming current can be instrumental to determine unambiguously the zeta potential and the characteristic surface roughness height.


title: Azimuthal modulation of electromagnetically-induced transparency by using asymmetrical Laguerre-Gaussian beams

abstract: Recently, the generation and detection of structured light field have drawn a great deal of attention, due to their importance in high-capacity optical data storage and quantum technology. In this letter, we explore the azimuthal modulation of optical transparency in a four-level double-V type quantum system near a plasmonic nanostructure. A Laguerre-Gaussian beam and the interaction of the system with free-space vacuum modes have been employed to create the phase-dependent absorption of a non-vortex probe field. First, we demonstrate how to identify the azimuthal index associated with the conventional LG beam via measuring the probe absorption so that the phase information of such a beam gets encoded on the spatially-dependent absorption profile with angularly-distributed lobes. Also, a spatially-varying optical transparency can be formed, due to the periodic variation of the absorption spectrum. Then, asymmetrical Laguerre-Gaussian beams are used to extend the selective spatial transparency mechanism to asymmetric spatially-structured windows, allowing for optical manipulation of spatial modes at arbitrary positions. Moreover, we investigate the influence of the asymmetric parameter on the features of the spatial inhomogeneities and show how the beams enable us to imprint the phase information of the orbital angular momentum at a desired position.


title: Monodromy of projections of hypersurfaces

abstract: Let $X$ be an irreducible, reduced complex projective hypersurface of degree $d$. A point $P$ not contained in $X$ is called uniform if the monodromy group of the projection of $X$ from $P$ is isomorphic to the symmetric group $S_d$. We prove that the locus of non--uniform points is finite when $X$ is smooth or a general projection of a smooth variety. In general, it is contained in a finite union of linear spaces of codimension at least $2$, except possibly for a special class of hypersurfaces with singular locus linear in codimension $1$. Moreover, we generalise a result of Fukasawa and Takahashi on the finiteness of Galois points.


title: FMore: An Incentive Scheme of Multi-dimensional Auction for Federated Learning in MEC

abstract: Promising federated learning coupled with Mobile Edge Computing (MEC) is considered as one of the most promising solutions to the AI-driven service provision. Plenty of studies focus on federated learning from the performance and security aspects, but they neglect the incentive mechanism. In MEC, edge nodes would not like to voluntarily participate in learning, and they differ in the provision of multi-dimensional resources, both of which might deteriorate the performance of federated learning. Also, lightweight schemes appeal to edge nodes in MEC. These features require the incentive mechanism to be well designed for MEC. In this paper, we present an incentive mechanism FMore with multi-dimensional procurement auction of K winners. Our proposal FMore not only is lightweight and incentive compatible, but also encourages more high-quality edge nodes with low cost to participate in learning and eventually improve the performance of federated learning. We also present theoretical results of Nash equilibrium strategy to edge nodes and employ the expected utility theory to provide guidance to the aggregator. Both extensive simulations and real-world experiments demonstrate that the proposed scheme can effectively reduce the training rounds and drastically improve the model accuracy for challenging AI tasks.


title: Online Statistics Teaching and Learning

abstract: For statistics courses at all levels, teaching and learning online poses challenges in different aspects. Particular online challenges include how to effectively and interactively conduct exploratory data analyses, how to incorporate statistical programming, how to include individual or team projects, and how to present mathematical derivations efficiently and effectively. This article draws from the authors' experience with seven different online statistics courses to address some of the aforementioned challenges. One course is an online exploratory data analysis course taught at Bowling Green State University. A second course is an upper level Bayesian statistics course taught at Vassar College and shared among 10 liberal arts colleges through a hybrid model. We alo describes a five-course MOOC specialization on Coursera, offered by Duke University.


title: Geometrical and energy scaling in the pulsed laser deposition plasma during epitaxial growth of FeSe thin films

abstract: Pulsed laser deposition (PLD) is a versatile technique for growing epitaxial heterostructures of a wide variety of novel materials combinations. Achieving low-defect layers with atomically sharp interfaces by PLD requires careful management of crystal growth conditions. Control over the flux of depositing particles, their kinetic energy, and the substrate temperature, is generally sufficient to obtain high-quality single crystal epitaxy in vacuum. In this article, we show that measurements of plasma parameters such as the electron temperature, the electron density, and the Mach number of the plasma fluid expansion may provide additional insights into the tuning of the crystal growth and processing environment of PLD. We report Langmuir probe measurements during growth of FeSe thin films on (100)-oriented SrTiO$_3$. We discuss two distinct plasma regimes that are accessible by KrF laser ablation of FeSe when pulse energies in the 20-630 mJ range are used. The two regimes can be characterized by the relative number of photons to absorbing centers in the laser-plasma interaction volume. Thin films grown under the conditions created by these two distinct plasma regimes are analyzed by x-ray diffraction and x-ray reflectivity and their epitaxial configurations correlated to the PLD conditions associated with these plasmas.


title: Radius-Dependent Homogeneous Strain in Uncoalesced GaN Nanowires

abstract: We investigate the strain state of ensembles of thin and nearly coalescence-free self-assembled GaN nanowires prepared by plasma-assisted molecular beam epitaxy on Ti/Al$_{2}$O$_{3}(0001)$ substrates. The shifts of Bragg peaks in high-resolution X-ray diffraction profiles reveal the presence of a homogeneous tensile strain in the out-of-plane direction. This strain is inversely proportional to the average nanowire radius and results from the surface stress acting on the nanowire sidewalls. The superposition of strain from nanowires with different radii in the same ensemble results in a broadening of the Bragg peaks that mimics an inhomogeneous strain on a macroscopic scale. The nanowire ensembles show a small blueshift of the bound-exciton transitions in photoluminescence spectra, reflecting the existence of a compensating in-plane compressive strain, as further supported by grazing incidence x-ray diffraction measurements carried out at a synchrotron. By combining X-ray diffraction and photoluminescence spectroscopy, the surface stress components $f_{x}$ and $f_{z}$ of the air-exposed GaN$\{1\bar100\}$ planes that constitute the nanowire sidewalls are determined experimentally to be 2.25 and $-0.7$~N/m, respectively.


title: Automatic Data Augmentation via Deep Reinforcement Learning for Effective Kidney Tumor Segmentation

abstract: Conventional data augmentation realized by performing simple pre-processing operations (\eg, rotation, crop, \etc) has been validated for its advantage in enhancing the performance for medical image segmentation. However, the data generated by these conventional augmentation methods are random and sometimes harmful to the subsequent segmentation. In this paper, we developed a novel automatic learning-based data augmentation method for medical image segmentation which models the augmentation task as a trial-and-error procedure using deep reinforcement learning (DRL). In our method, we innovatively combine the data augmentation module and the subsequent segmentation module in an end-to-end training manner with a consistent loss. Specifically, the best sequential combination of different basic operations is automatically learned by directly maximizing the performance improvement (\ie, Dice ratio) on the available validation set. We extensively evaluated our method on CT kidney tumor segmentation which validated the promising results of our method.


title: Blowing-up solutions of a time-space fractional semi-linear equation with a structural damping and a nonlocal in time nonlinearity

abstract: In this paper, we investigate the semilinear equation with a time-space fractional structural damping and a nonlocal in time nonlinearity \begin{equation*} {\mathbf{D}}_{0|t}^{1+\alpha_1}u + (-\Delta)^\sigma u+(-\Delta )^\delta\mathbf{D}_{0|t}^{\alpha _2} u = I_{0|t}^{1-\gamma }|u|^{p}, \qquad (t,x)\in (0,\infty) \times \mathbb{R}^N, \end{equation*} where $p>1$, $\alpha_i, \gamma\ in (0,1)$, $\delta, \sigma \in (0,1)$, ${\mathbf{D}}_{0|t}^{\alpha_i}$ is the Caputo fractional derivative and $I_{0|t}^{1-\gamma }$ is the Riemann-Liouville fractional integral of order $1-\gamma$. We prove the non-existence of global solutions if \begin{equation*} 1<p\leqslant \frac{2(2+\alpha_1-\gamma)}{(\frac{\alpha_1+1}{\sigma} N+2\gamma-2\alpha_1-2)_+ }+1, \end{equation*} for any space dimension $N\geqslant 1.$ Then, we extend the result to the system \begin{align*} &{\mathbf{D}}_{0|t}^{1+\alpha_1}u + (-\Delta)^{\sigma_1} u + (-\Delta )^{\delta_1}{\mathbf{D}}_{0|t}^{\alpha_2} u = I_{0|t}^{1-\gamma _{1}}|v|^{p},\qquad (t,x)\in (0,\infty) \times \mathbb{R}^N, \\ &{\mathbf{D}}_{0|t}^{1+\beta _{1}}v+(-\Delta)^{\sigma_2} v + (-\Delta )^{\delta_2}{\mathbf{D}}_{0|t}^{\beta_2}v = I_{0|t}^{1-\gamma_2}|u|^{q},\qquad (t,x)\in (0,\infty )\times \mathbb{R}^{N}, \end{align*} where $p,q>1$, $0<\delta_i$, $\sigma_i<1$ and $\gamma_2\in (0,1)$. Also, we present the necessary conditions for the existence of local or global solutions.


title: The candy wrapper problem -- a temporal multiscale approach for pde/pde systems

abstract: We describe a temporal multiscale approach for the simulation of long-term processes with short-term influences involving partial differential equations. The specific problem under consideration is a growth process in blood vessels. The \emph{Candy Wrapper Process} describes a restenosis in a vessel that has previously be widened by inserting a stent. The development of a new stenosis takes place on a long time horizon (months) while the acting forces are mainly given by the pulsating blood flow. We describe a coupled pde model and a finite element simulation that is used as basis for our multiscale approach, which is based on averaging the long scale equation and approximating the fast scale impact by localized periodic-in-time problems. Numerical test cases in prototypical 3d configurations demonstrate the power of the approach.


title: Structural Combinatorial of Network Information System of Systems based on Evolutionary Optimization Method

abstract: The network information system is a military information network system with evolution characteristics. Evolution is a process of replacement between disorder and order, chaos and equilibrium. Given that the concept of evolution originates from biological systems, in this article, the evolution of network information architecture is analyzed by genetic algorithms, and the network information architecture is represented by chromosomes. Besides, the genetic algorithm is also applied to find the optimal chromosome in the architecture space. The evolutionary simulation is used to predict the optimal scheme of the network information architecture and provide a reference for system construction.


title: Compression with wildcards: All spanning trees with prescribed vertex-degrees

abstract: By processing all minimal cutsets of a graph G, and by using novel wildcards, all spanning trees of G can be compactly encoded. Surprisingly, a 1986 algorithm of Winter seems to achieve (Conjecture 2) exactly the same compression, although in totally unrelated ways! Although Winter's algorithm is faster, our method adapts to relevant variations of the task, e.g. to the enumeration of only those spanning trees that obey prescribed degree conditions


title: Robust Multimodal Brain Tumor Segmentation via Feature Disentanglement and Gated Fusion

abstract: Accurate medical image segmentation commonly requires effective learning of the complementary information from multimodal data. However, in clinical practice, we often encounter the problem of missing imaging modalities. We tackle this challenge and propose a novel multimodal segmentation framework which is robust to the absence of imaging modalities. Our network uses feature disentanglement to decompose the input modalities into the modality-specific appearance code, which uniquely sticks to each modality, and the modality-invariant content code, which absorbs multimodal information for the segmentation task. With enhanced modality-invariance, the disentangled content code from each modality is fused into a shared representation which gains robustness to missing data. The fusion is achieved via a learning-based strategy to gate the contribution of different modalities at different locations. We validate our method on the important yet challenging multimodal brain tumor segmentation task with the BRATS challenge dataset. With competitive performance to the state-of-the-art approaches for full modality, our method achieves outstanding robustness under various missing modality(ies) situations, significantly exceeding the state-of-the-art method by over 16% in average for Dice on whole tumor segmentation.


title: Knot diagrams on a punctured sphere as a model of string figures

abstract: A string figure is topologically a trivial knot lying on an imaginary plane orthogonal to the fingers with some crossings. The fingers prevent cancellation of these crossings. As a mathematical model of string figure we consider a knot diagram on the $xy$-plane in $xyz$-space missing some straight lines parallel to the $z$-axis. These straight lines correspond to fingers. We study minimal number of crossings of these knot diagrams under Reidemeister moves missing these lines.


title: Actively Mapping Industrial Structures with Information Gain-Based Planning on a Quadruped Robot

abstract: In this paper, we develop an online active mapping system to enable a quadruped robot to autonomously survey large physical structures. We describe the perception, planning and control modules needed to scan and reconstruct an object of interest, without requiring a prior model. The system builds a voxel representation of the object, and iteratively determines the Next-Best-View (NBV) to extend the representation, according to both the reconstruction itself and to avoid collisions with the environment. By computing the expected information gain of a set of candidate scan locations sampled on the as-sensed terrain map, as well as the cost of reaching these candidates, the robot decides the NBV for further exploration. The robot plans an optimal path towards the NBV, avoiding obstacles and un-traversable terrain. Experimental results on both simulated and real-world environments show the capability and efficiency of our system. Finally we present a full system demonstration on the real robot, the ANYbotics ANYmal, autonomously reconstructing a building facade and an industrial structure.


title: Robotic modeling of snake traversing large, smooth obstacles reveals stability benefits of body compliance

abstract: Snakes can move through almost any terrain. Although their locomotion on flat surfaces using planar gaits is inherently stable, when snakes deform their body out of plane to traverse complex terrain, maintaining stability becomes a challenge. On trees and desert dunes, snakes grip branches or brace against depressed sand for stability. However, how they stably surmount obstacles like boulders too large and smooth to gain such anchor points is less understood. Similarly, snake robots are challenged to stably traverse large, smooth obstacles for search and rescue and building inspection. Our recent study discovered that snakes combine body lateral undulation and cantilevering to stably traverse large steps. Here, we developed a snake robot with this gait and snake-like anisotropic friction and used it as a physical model to understand stability principles. The robot traversed steps as high as a third of its body length rapidly and stably. However, on higher steps, it was more likely to fail due to more frequent rolling and flipping over, which was absent in the snake with a compliant body. Adding body compliance reduced the robot roll instability by statistically improving surface contact, without reducing speed. Besides advancing understanding of snake locomotion, our robot achieved high traversal speed surpassing most previous snake robots and approaching snakes, while maintaining high traversal probability.


title: Workshop Report: Brightest Light Initiative (March 27-29 2019, OSA Headquarters, Washington, D.C.)

abstract: This Brightest Light Initiative (BLI) Workshop Report captures the important research ideas and recommendations for enabling that work developed by over 100 leading scientists at a community-initiated workshop held March 27-29, 2019 in Washington, DC. Workshop attendees developed an understanding of key opportunities, as well as gaps in current technologies and capabilities, for science enabled by the highest-intensity lasers.

categories: physics.plasm-ph physics.acc-ph physics.atom-ph physics.optics quant-ph


title: Connections between statistical practice in elementary particle physics and the severity concept as discussed in Mayo's Statistical Inference as Severe Testing

abstract: For many years, philosopher-of-statistics Deborah Mayo has been advocating the concept of severe testing as a key part of hypothesis testing. Her recent book, Statistical Inference as Severe Testing, is a comprehensive exposition of her arguments in the context of a historical study of many threads of statistical inference, both frequentist and Bayesian. Her foundational point of view is called error statistics, emphasizing frequentist evaluation of the errors called Type I and Type II in the Neyman-Pearson theory of frequentist hypothesis testing. Since the field of elementary particle physics (also known as high energy physics) has strong traditions in frequentist inference, one might expect that something like the severity concept was independently developed in the field. Indeed, I find that, at least operationally (numerically), we high-energy physicists have long interpreted data in ways that map directly onto severity. Whether or not we subscribe to Mayo's philosophical interpretations of severity is a more complicated story that I do not address here.


title: The bursty star formation history of the Fornax dwarf spheroidal galaxy revealed with the HST

abstract: We present a new derivation of the star formation history (SFH) of the dSph galaxy Fornax in two central regions, characterised by unprecedented precision and age resolution. It reveals that star formation has proceeded in sharp bursts separated by periods of low-level or quiescent activity. The SFH was derived through colour-magnitude diagram (CMD) fitting of two extremely deep Hubble Space Telescope CMDs, sampling the centre and one core radius. The attained age resolution allowed us to single out a major star formation episode at early times, a second strong burst $4.6\pm0.4$ Gyr ago and recent intermittent episodes $\sim2-0.2$ Gyr ago. Detailed testing with mock stellar populations was used to estimate the duration of the main bursts and study the occurrence of low-level star formation between them. The SFHs in both regions show common features, with activity at the same epochs and similar age-metallicity relationship. However, clear indications of a spatial gradient were also found, with mean age increasing with radius and star formation episodes being more prolonged in the centre. While some galaxy evolution models predict bursty SFHs in dwarf galaxies and thus a secular origin of the observed SFH cannot be excluded in Fornax, other evidence points to possible mergers or interactions as the cause of its bursty SFH. In particular, we calculated the Fornax orbit relative to the closest dwarfs and the Milky Way and observed a correspondence between the main intermediate-age and young events and peri-passages of Fornax around the Milky Way, possibly indicating tidally-induced star formation.


title: $L^p$ eigenfunction bounds for fractional Schr\"odinger operators on manifolds

abstract: This paper is dedicated to $L^p$ bounds on eigenfunctions of a Sch\"odinger-type operator $(-\Delta_g)^{\alpha/2} +V$ on closed Riemannian manifolds for critically singular potentials $V$. The operator $(-\Delta_g)^{\alpha/2}$ is defined spectrally in terms of the eigenfunctions of $-\Delta_g$. We obtain also quasimodes and spectral clusters estimates. As an application, we derive Strichartz estimates for the fractional wave equation $(\partial_t^2+(-\Delta_g)^{\alpha/2}+V)u=0$. The wave kernel techniques recently developed by Bourgain-Shao-Sogge-Yao and Shao-Yao play a key role in this paper. We construct a new reproducing operator with several local operators and some good error terms. Moreover, we shall prove that these local operators satisfy certain variable coefficient versions of the "uniform Sobolev estimates" by Kenig-Ruiz-Sogge. These enable us to handle the critically singular potentials $V$ and prove the quasimode estimates.


title: Bayesian Computing in the Undergraduate Statistics Curriculum

abstract: Bayesian statistics has gained great momentum since the computational developments of the 1990s. Gradually, advances in Bayesian methodology and software have made Bayesian techniques much more accessible to applied statisticians and, in turn, have potentially transformed Bayesian education at the undergraduate level. This article provides an overview on the various options for implementing Bayesian computational methods motivated to achieve particular learning outcomes. The advantages and disadvantages of each computational method are described based on the authors' experience in using these methods in the classroom. The goal is to present guidance on the choice of computation for the instructors who are introducing Bayesian methods in their undergraduate statistics curriculum.


title: Ill-posedness for the Maxwell-Dirac system below charge in space dimension three and lower

abstract: The Maxwell-Dirac system describes the interaction of an electron with its self-induced electromagnetic field. In space dimension $d=3$ the system is charge-critical, that is, $L^2$-critical for the spinor with respect to scaling, and local well-posedness is known almost down to the critical regularity. In the charge-subcritical dimensions $d=1,2$, global well-posedness is known in the charge class. Here we prove that these results are sharp (or almost sharp, if $d=3$), by demonstrating ill-posedness below the charge regularity. In fact, for $d \le 3$ we exhibit a spinor datum belonging to $H^s(\mathbb R^d)$ for $s<0$, and to $L^p(\mathbb R^d)$ for $1 \le p < 2$, but not to $L^2(\mathbb R^d)$, which does not admit any local solution that can be approximated by smooth solutions in a reasonable sense.


title: Safe Screening for the Generalized Conditional Gradient Method

abstract: The conditional gradient method (CGM) has been widely used for fast sparse approximation, having a low per iteration computational cost for structured sparse regularizers. We explore the sparsity acquiring properties of a generalized CGM (gCGM), where the constraint is replaced by a penalty function based on a gauge penalty; this can be done without significantly increasing the per-iteration computation, and applies to general notions of sparsity. Without assuming bounded iterates, we show $O(1/t)$ convergence of the function values and gap of gCGM. We couple this with a safe screening rule, and show that at a rate $O(1/(t\delta^2))$, the screened support matches the support at the solution, where $\delta \geq 0$ measures how close the problem is to being degenerate. In our experiments, we show that the gCGM for these modified penalties have similar feature selection properties as common penalties, but with potentially more stability over the choice of hyperparameter.


title: Joint Transmission and Computing Scheduling for Status Update with Mobile Edge Computing

abstract: Age of Information (AoI), defined as the time elapsed since the generation of the latest received update, is a promising performance metric to measure data freshness for real-time status monitoring. In many applications, status information needs to be extracted through computing, which can be processed at an edge server enabled by mobile edge computing (MEC). In this paper, we aim to minimize the average AoI within a given deadline by jointly scheduling the transmissions and computations of a series of update packets with deterministic transmission and computing times. The main analytical results are summarized as follows. Firstly, the minimum deadline to guarantee the successful transmission and computing of all packets is given. Secondly, a \emph{no-wait computing} policy which intuitively attains the minimum AoI is introduced, and the feasibility condition of the policy is derived. Finally, a closed-form optimal scheduling policy is obtained on the condition that the deadline exceeds a certain threshold. The behavior of the optimal transmission and computing policy is illustrated by numerical results with different values of the deadline, which validates the analytical results.


title: Linearly dependent and concise subsets of a Segre variety depending on k factors

abstract: We study linearly dependent subsets with prescribed cardinality, $s$, of a multiprojective space. If the set $S$ is a circuit, we give an upper bound on the number of factors of the minimal multiprojective space containing $S$, while if $S$ has higher dependency this may be not true without strong assumptions. We describe the dependent subsets $S$ with $\#S=6$.


title: General theory of interpolation error estimates on anisotropic meshes

abstract: We propose a general theory of estimating interpolation error for smooth functions in two and three dimensions. In our theory, the error of interpolation is bound in terms of the diameter of a simplex and a geometric parameter. In the two-dimensional case, our geometric parameter is equivalent to the circumradius of a triangle. In the three-dimensional case, our geometric parameter also represents the flatness of a tetrahedron. Through the introduction of the geometric parameter, the error estimates newly obtained can be applied to cases that violate the maximum-angle condition.


title: Checking Phylogenetic Decisiveness in Theory and in Practice

abstract: Suppose we have a set $X$ consisting of $n$ taxa and we are given information from $k$ loci from which to construct a phylogeny for $X$. Each locus offers information for only a fraction of the taxa. The question is whether this data suffices to construct a reliable phylogeny. The decisiveness problem expresses this question combinatorially. Although a precise characterization of decisiveness is known, the complexity of the problem is open. Here we relate decisiveness to a hypergraph coloring problem. We use this idea to (1) obtain lower bounds on the amount of coverage needed to achieve decisiveness, (2) devise an exact algorithm for decisiveness, (3) develop problem reduction rules, and use them to obtain efficient algorithms for inputs with few loci, and (4) devise an integer linear programming formulation of the decisiveness problem, which allows us to analyze data sets that arise in practice.


title: Constructing fast approximate eigenspaces with application to the fast graph Fourier transforms

abstract: We investigate numerically efficient approximations of eigenspaces associated to symmetric and general matrices. The eigenspaces are factored into a fixed number of fundamental components that can be efficiently manipulated (we consider extended orthogonal Givens or scaling and shear transformations). The number of these components controls the trade-off between approximation accuracy and the computational complexity of projecting on the eigenspaces. We write minimization problems for the single fundamental components and provide closed-form solutions. Then we propose algorithms that iterative update all these components until convergence. We show results on random matrices and an application on the approximation of graph Fourier transforms for directed and undirected graphs.


title: Stochastic production planning with regime switching

abstract: This paper considers a stochastic production planning problem with regime switching. There are two regimes corresponding to different economic cycles. A factory is planning its production so as to minimize production costs. We analyze this problem through the value function approach. The optimal production is characterized through the solution of an elliptic system of partial differential equations which is shown to have a solution.


title: Testing the Agreement of Trees with Internal Labels

abstract: The input to the agreement problem is a collection $P = \{T_1, T_2, \dots , T_k\}$ of phylogenetic trees, called input trees, over partially overlapping sets of taxa. The question is whether there exists a tree $T$, called an agreement tree, whose taxon set is the union of the taxon sets of the input trees, such that for each $i \in \{1, 2, \dots , k\}$, the restriction of $T$ to the taxon set of $T_i$ is isomorphic to $T_i$. We give a $O(n k (\sum_{i \in [k]} d_i + \log^2(nk)))$ algorithm for a generalization of the agreement problem in which the input trees may have internal labels, where $n$ is the total number of distinct taxa in $P$, $k$ is the number of trees in $P$, and $d_i$ is the maximum number of children of a node in $T_i$.


title: Operator inference for non-intrusive model reduction of systems with non-polynomial nonlinear terms

abstract: This work presents a non-intrusive model reduction method to learn low-dimensional models of dynamical systems with non-polynomial nonlinear terms that are spatially local and that are given in analytic form. In contrast to state-of-the-art model reduction methods that are intrusive and thus require full knowledge of the governing equations and the operators of a full model of the discretized dynamical system, the proposed approach requires only the non-polynomial terms in analytic form and learns the rest of the dynamics from snapshots computed with a potentially black-box full-model solver. The proposed method learns operators for the linear and polynomially nonlinear dynamics via a least-squares problem, where the given non-polynomial terms are incorporated in the right-hand side. The least-squares problem is linear and thus can be solved efficiently in practice. The proposed method is demonstrated on three problems governed by partial differential equations, namely the diffusion-reaction Chafee-Infante model, a tubular reactor model for reactive flows, and a batch-chromatography model that describes a chemical separation process. The numerical results provide evidence that the proposed approach learns reduced models that achieve comparable accuracy as models constructed with state-of-the-art intrusive model reduction methods that require full knowledge of the governing equations.


title: Realization of Lie algebras of high dimension via pseudo-bosonic operators

abstract: The present paper is the third contribution of a series of works, where we investigate pseudo--bosonic operators and their connections with finite dimensional Lie algebras. We show that all finite dimensional nilpotent Lie algebras (over the complex field) can be realized by central extensions of Lie algebras of pseudo-bosonic operators. This result is interesting, because it provides new examples of dynamical systems for nilpotent Lie algebras of any dimension. One could ask whether these operators are intrinsic with the notion of nilpotence or not, but this is false. In fact we exibit both a simple Lie algebra and a solvable nonnilpotent Lie algebra, which can be realized in terms of pseudo-bosonic operators.


title: The Higman operations and embeddings of recursive groups

abstract: In the context of Higman embeddings of recursive groups into finitely presented groups we suggest an algorithm which uses Higman operations to explicitly constructs the specific recursively enumerable sets of integer sequences arising during the embeddings. This makes the constructive Higman embedding a doable task for certain wide classes of groups. Specific auxiliary operations are introduced to make the work with Higman operations a simpler and more intuitive procedure. Also, an automated mechanism of constructive embeddings of countable groups into 2-generator groups is mentioned.


title: Thermodynamics in Rotating Anti-de Sitter Black Holes with Massive Scalar Field in Three Dimensions

abstract: We investigate the tendency in variations of the dual CFT$_2$ when a rotating AdS$_3$ black hole changes by fluxes transferred by the scattering of a massive scalar field according to the AdS/CFT correspondence. The conserved quantities of the black hole are definitely constrained by the extremal condition. Moreover, the laws of thermodynamics provide a direction for the changes in the conserved quantities, so the black hole cannot be extremal under the scattering. This is naturally preferred. By the relationship between the rotating AdS$_3$ black hole and dual CFT$_2$, we find that such changes in the black hole constrain variations in eigenstates of the dual CFT$_2$. Furthermore, the tendency in the variations is closely related to the laws of thermodynamics.


title: Microreversibility and quantum transport in Aharonov-Bohm rings

abstract: The consequences of microreversibility for the linear and nonlinear transport properties of systems subjected to external magnetic fields are systematically investigated in Aharonov-Bohm rings connected to two, three, and four terminals. Within the independent electron approximation, the cumulant generating function, which fully specifies the statistics of the nonequilibrium currents, is expressed in terms of the scattering matrix of these circuits. The time-reversal symmetry relations up to the third responses of the currents and the fourth cumulants are analytically investigated and numerically tested as a function of the magnetic flux. The validity of such relations is thus firmly confirmed in this class of open quantum systems.


title: Impurity temperature screening in stellarators close to quasisymmetry

abstract: Impurity temperature screening is a favorable neoclassical phenomenon involving an outward radial flux of impurity ions from the core of fusion devices. Quasisymmetric magnetic fields are intrinsically ambipolar and give rise to temperature screening for low enough $\eta^{-1}\equiv d\ln n/d\ln T$. In contrast, neoclassical fluxes in generic stellarators will depend on the radial electric field, which is predicted to be inward for ion-root plasmas, potentially leading to impurity accumulation. Here we examine the impurity particle flux in a number of approximately quasisymmetric stellarator configurations and parameter regimes while varying the amount of symmetry-breaking in the magnetic field. The amount of symmetry-breaking (from perfect quasisymmetry) required to lose temperature screening is dependent on the particular configuration and flux surface. Neoclassical fluxes have been obtained using the SFINCS drift-kinetic equation solver with the electrostatic potential $\Phi=\Phi(\psi)$, where $\psi$ is a flux surface label. Results indicate that achieving temperature screening is possible, but unlikely, at reactor-relevant conditions in the core. Thus, the small departures from symmetry in nominally quasisymmetric stellarators are large enough to significantly alter the neoclassical impurity transport. A further look at these fluxes when compared to a gyro-Bohm turbulence estimate suggests that neoclassical fluxes are small in optimized configurations compared to turbulent fluxes. Therefore, although neoclassical impurity accumulation is expected in most situations, the strength of the turbulence may render this irrelevant.


title: Performance of Oscillating Plasma Thrusters

abstract: Traditional expressions and definitions describing performance of plasma thrusters, including the thrust, specific impulse, and the thruster efficiency, assume a steady state plasma flow with a constant flow velocity. However, it is very common for these thrusters that the plasma exhibits unstable behavior resulting in time-variations of the thrust and the exhaust velocity. For example, in Hall thrusters, the ionization instability leads to strong oscillations of the discharge current (so-called breathing oscillations), plasma density, ion energy, and as a result the ion flow. In this paper, we revisit the formulation of the thrust and the thrust efficiency to account for time variations of the ion parameters including the phase shift between the ion energy and the ion flow. For sinusoidal oscillations it was found thrust can potentially change more than 20%. It is shown that by modulating ion energy at specific amplitudes, thrust can be maximized in such regimes. Finally, an expression for the thruster efficiency of the modulating thruster is derived to show a mechanism for inefficiencies in such thrusters.


title: Numerical Analysis of a High-Order Scheme for Nonlinear Fractional Differential Equations with Uniform Accuracy

abstract: We introduce a high-order numerical scheme for fractional ordinary differential equations with the Caputo derivative. The method is developed by dividing the domain into a number of subintervals, and applying the quadratic interpolation on each subinterval. The method is shown to be unconditionally stable, and for general nonlinear equations, the uniform sharp numerical order $3-\nu$ can be rigorously proven for sufficiently smooth solutions at all time steps. The proof provides a general guide for proving the sharp order for higher-order schemes in the nonlinear case. Some numerical examples are given to validate our theoretical results.


title: Multiplicative Decomposition of Heterogeneity in Mixtures of Continuous Distributions

abstract: A system's heterogeneity (\textit{diversity}) is the effective size of its event space, and can be quantified using the R\'enyi family of indices (also known as Hill numbers in ecology or Hannah-Kay indices in economics), which are indexed by an elasticity parameter $q \geq 0$. Under these indices, the heterogeneity of a composite system (the $\gamma$-heterogeneity) is decomposable into heterogeneity arising from variation \textit{within} and \textit{between} component subsystems (the $\alpha$- and $\beta$-heterogeneity, respectively). Since the average heterogeneity of a component subsystem should not be greater than that of the pooled system, we require that $\gamma \geq \alpha$. There exists a multiplicative decomposition for R\'enyi heterogeneity of composite systems with discrete event spaces, but less attention has been paid to decomposition in the continuous setting. We therefore describe multiplicative decomposition of the R\'enyi heterogeneity for continuous mixture distributions under parametric and non-parametric pooling assumptions. Under non-parametric pooling, the $\gamma$-heterogeneity must often be estimated numerically, but the multiplicative decomposition holds such that $\gamma \geq \alpha$ for $q > 0$. Conversely, under parametric pooling, $\gamma$-heterogeneity can be computed efficiently in closed-form, but the $\gamma \geq \alpha$ condition holds reliably only at $q=1$. Our findings will further contribute to heterogeneity measurement in continuous systems.


title: Partially Observed Dynamic Tensor Response Regression

abstract: In modern data science, dynamic tensor data is prevailing in numerous applications. An important task is to characterize the relationship between such dynamic tensor and external covariates. However, the tensor data is often only partially observed, rendering many existing methods inapplicable. In this article, we develop a regression model with partially observed dynamic tensor as the response and external covariates as the predictor. We introduce the low-rank, sparsity and fusion structures on the regression coefficient tensor, and consider a loss function projected over the observed entries. We develop an efficient non-convex alternating updating algorithm, and derive the finite-sample error bound of the actual estimator from each step of our optimization algorithm. Unobserved entries in tensor response have imposed serious challenges. As a result, our proposal differs considerably in terms of estimation algorithm, regularity conditions, as well as theoretical properties, compared to the existing tensor completion or tensor response regression solutions. We illustrate the efficacy of our proposed method using simulations, and two real applications, a neuroimaging dementia study and a digital advertising study.


title: Model-assisted estimation through random forests in finite population sampling

abstract: In surveys, the interest lies in estimating finite population parameters such as population totals and means. In most surveys, some auxiliary information is available at the estimation stage. This information may be incorporated in the estimation procedures to increase their precision. In this article, we use random forests to estimate the functional relationship between the survey variable and the auxiliary variables. In recent years, random forests have become attractive as National Statistical Offices have now access to a variety of data sources, potentially exhibiting a large number of observations on a large number of variables. We establish the theoretical properties of model-assisted procedures based on random forests and derive corresponding variance estimators. A model-calibration procedure for handling multiple survey variables is also discussed. The results of a simulation study suggest that the proposed point and estimation procedures perform well in term of bias, efficiency, and coverage of normal-based confidence intervals, in a wide variety of settings. Finally, we apply the proposed methods using data on radio audiences collected by M\'ediam\'etrie, a French audience company.


title: Amortised Learning by Wake-Sleep

abstract: Models that employ latent variables to capture structure in observed data lie at the heart of many current unsupervised learning algorithms, but exact maximum-likelihood learning for powerful and flexible latent-variable models is almost always intractable. Thus, state-of-the-art approaches either abandon the maximum-likelihood framework entirely, or else rely on a variety of variational approximations to the posterior distribution over the latents. Here, we propose an alternative approach that we call amortised learning. Rather than computing an approximation to the posterior over latents, we use a wake-sleep Monte-Carlo strategy to learn a function that directly estimates the maximum-likelihood parameter updates. Amortised learning is possible whenever samples of latents and observations can be simulated from the generative model, treating the model as a "black box". We demonstrate its effectiveness on a wide range of complex models, including those with latents that are discrete or supported on non-Euclidean spaces.


title: High-resolution 4D acquisition of freely swimming human sperm cells without staining

abstract: We present a new acquisition method that enables high-resolution, fine-detail full reconstruction of the three-dimensional movement and structure of individual human sperm cells swimming freely. We achieve both retrieval of the three-dimensional refractive-index profile of the sperm head, revealing its fine internal organelles and time-varying orientation, and the detailed four-dimensional localization of the thin, highly-dynamic flagellum of the sperm cell. Live human sperm cells were acquired during free swim using a high-speed off-axis holographic system that does not require any moving elements or cell staining. The reconstruction is based solely on the natural movement of the sperm cell and a novel set of algorithms, enabling the detailed four-dimensional recovery. Using this refractive-index imaging approach, we believe we have detected an area in the cell that is attributed to the centriole. This method has great potential for both biological assays and clinical use of intact sperm cells.


title: Generalized theory of pseudomodes for exact descriptions of non-Markovian quantum processes

abstract: We develop an exact framework to describe the non-Markovian dynamics of an open quantum system interacting with an environment modeled by a generalized spectral density function. The approach relies on mapping the initial system onto an auxiliary configuration, comprising the original open system coupled to a small number of discrete modes, which in turn are each coupled to an independent Markovian reservoir. Based on the connection between the discrete modes and the poles of the spectral density function, we show how expanding the system using the discrete modes allows for the full inclusion non-Markovian effects within an enlarged open system whose dynamics is governed by an exact Lindblad master equation. Initially we apply this result to obtain a generalization of the pseudomode method [B. M. Garraway, Phys. Rev. A 55, 2290 (1997)] in cases where the spectral density function has a Lorentzian structure. For many other types of spectral density function, we extend our proof to show that an open system dynamics may be modeled physically using discrete modes which admit a non-Hermitian coupling to the system, and for such cases determine the equivalent master equation to no longer be of Lindblad form. For applications involving two discrete modes, we demonstrate how to convert between pathological and Lindblad forms of the master equation using the techniques of the pseudomode method.