title: Predicting maximal gaps in sets of primes

abstract: Let $q>r\ge1$ be coprime integers. Let ${\mathbb P}_c={\mathbb P}_c(q,r,{\cal H})$ be an increasing sequence of primes $p$ satisfying two conditions: (i) $p\equiv r$ (mod $q$) and (ii) $p$ starts a prime $k$-tuple with a given pattern ${\cal H}$. Let $\pi_c(x)$ be the number of primes in ${\mathbb P}_c$ not exceeding $x$. We heuristically derive formulas predicting the growth trend of the maximal gap $G_c(x)=\max_{p'\le x}(p'-p)$ between successive primes $p,p'\in{\mathbb P}_c$. Extensive computations for primes up to $10^{14}$ show that a simple trend formula $$G_c(x) \sim {x\over\pi_c(x)}\cdot(\log \pi_c(x) + O_k(1))$$ works well for maximal gaps between initial primes of $k$-tuples with $k\ge2$ (e.g., twin primes, prime triplets, etc.) in residue class $r$ (mod $q$). For $k=1$, however, a more sophisticated formula $$G_c(x) \sim {x\over\pi_c(x)}\cdot\big(\log{\pi_c^2(x)\over x}+O(\log q)\big)$$ gives a better prediction of maximal gap sizes. The latter includes the important special case of maximal gaps in the sequence of all primes ($k=1$, $q=2$, $r=1$). The distribution of appropriately rescaled maximal gaps $G_c(x)$ is close to the Gumbel extreme value distribution. Computations suggest that almost all maximal gaps satisfy a generalized strong form of Cramer's conjecture. We also conjecture that the number of maximal gaps between primes in ${\mathbb P}_c$ below $x$ is $O_k(\log x)$.


title: Automatic classification of geologic units in seismic images using partially interpreted examples

abstract: Geologic interpretation of large seismic stacked or migrated seismic images can be a time-consuming task for seismic interpreters. Neural network based semantic segmentation provides fast and automatic interpretations, provided a sufficient number of example interpretations are available. Networks that map from image-to-image emerged recently as powerful tools for automatic segmentation, but standard implementations require fully interpreted examples. Generating training labels for large images manually is time consuming. We introduce a partial loss-function and labeling strategies such that networks can learn from partially interpreted seismic images. This strategy requires only a small number of annotated pixels per seismic image. Tests on seismic images and interpretation information from the Sea of Ireland show that we obtain high-quality predicted interpretations from a small number of large seismic images. The combination of a partial-loss function, a multi-resolution network that explicitly takes small and large-scale geological features into account, and new labeling strategies make neural networks a more practical tool for automatic seismic interpretation.


title: Localizing the $E_2$ page of the Adams spectral sequence

abstract: There is only one nontrivial localization of $\pi_*S_{(p)}$ (the chromatic localization at $v_0=p$), but there are infinitely many nontrivial localizations of the Adams $E_2$ page for the sphere. The first non-nilpotent element in the $E_2$ page after $v_0$ is $b_{10}\in \mathrm{Ext}_A^{2p(p-1)-2}(\mathbb{F}_p,\mathbb{F}_p)$. We work at $p=3$ and study $b_{10}^{-1}\mathrm{Ext}_P(\mathbb{F}_3,\mathbb{F}_3)$ (where $P$ is the algebra of dual reduced powers), which agrees with the infinite summand $\mathrm{Ext}_P(\mathbb{F}_3,\mathbb{F}_3)$ of $\mathrm{Ext}_A(\mathbb{F}_3,\mathbb{F}_3)$ above a line of slope ${1\over 23}$. We compute up to the $E_9$ page of an Adams spectral sequence in the category $\mathrm{Stable}(P)$ converging to $b_{10}^{-1}\mathrm{Ext}_P(\mathbb{F}_3,\mathbb{F}_3)$, and conjecture that the spectral sequence collapses at $E_9$. We also give a complete calculation of $b_{10}^{-1}\mathrm{Ext}_P^*(\mathbb{F}_3,\mathbb{F}_3[\xi_1^3])$.


title: Semi-interactive Attention Network for Answer Understanding in Reverse-QA

abstract: Question answering (QA) is an important natural language processing (NLP) task and has received much attention in academic research and industry communities. Existing QA studies assume that questions are raised by humans and answers are generated by machines. Nevertheless, in many real applications, machines are also required to determine human needs or perceive human states. In such scenarios, machines may proactively raise questions and humans supply answers. Subsequently, machines should attempt to understand the true meaning of these answers. This new QA approach is called reverse-QA (rQA) throughout this paper. In this work, the human answer understanding problem is investigated and solved by classifying the answers into predefined answer-label categories (e.g., True, False, Uncertain). To explore the relationships between questions and answers, we use the interactive attention network (IAN) model and propose an improved structure called semi-interactive attention network (Semi-IAN). Two Chinese data sets for rQA are compiled. We evaluate several conventional text classification models for comparison, and experimental results indicate the promising performance of our proposed models.


title: Recovery of parameters of fast nonlocal heat transport in magnetic fusion plasmas: testing a model of waves with high internal reflections

abstract: Our analysis of the model [J. Phys.: Conf. Ser. 941 (2017) 012008] elaborated for interpreting the initial stage of the fast nonlocal transport events, which exhibit immediate response, in the heat diffusion time scale, of the spatial profile of electron temperature to its local perturbation, shows that the nonlocal transport by electromagnetic (EM) waves needs too high reflectivity of vacuum vessel walls to describe the experimental data. Here we try another model, which assumes high internal reflections and is compatible with the "wild cable" transport of TEM waves along magnetically-bound skeletal nanostructures. An inverse problem for recovery of the source and sink of waves, and internal reflectivity, is formulated and solved. Preliminary results of analyzing the data from tokamaks JET and TFTR, and stellarator LHD are presented.


title: List Decoding Random Euclidean Codes and Infinite Constellations

abstract: We study the list decodability of different ensembles of codes over the real alphabet under the assumption of an omniscient adversary. It is a well-known result that when the source and the adversary have power constraints $ P $ and $ N $ respectively, the list decoding capacity is equal to $ \frac{1}{2}\log\frac{P}{N} $. Random spherical codes achieve constant list sizes, and the goal of the present paper is to obtain a better understanding of the smallest achievable list size as a function of the gap to capacity. We show a reduction from arbitrary codes to spherical codes, and derive a lower bound on the list size of typical random spherical codes. We also give an upper bound on the list size achievable using nested Construction-A lattices and infinite Construction-A lattices. We then define and study a class of infinite constellations that generalize Construction-A lattices and prove upper and lower bounds for the same. Other goodness properties such as packing goodness and AWGN goodness of infinite constellations are proved along the way. Finally, we consider random lattices sampled from the Haar distribution and show that if a certain number-theoretic conjecture is true, then the list size grows as a polynomial function of the gap-to-capacity.


title: Optimization of Survey Weights under a Large Number of Conflicting Constraints

abstract: In the analysis of survey data, sampling weights are needed for consistent estimation of the population. However, the original inverse probability weights from the survey sample design are typically modified to account for non-response, to increase efficiency by incorporating auxiliary population information, and to reduce the variability in estimates due to extreme weights. It is often the case that no single set of weights can be found which successfully incorporates all of these modifications because together they induce a large number of constraints and restrictions on the feasible solution space. For example, a unique combination of categorical variables may not be present in the sample data, even if the corresponding population level information is available. Additional requirements for weights to fall within specified ranges may also lead to fewer population level adjustments being incorporated. We present a framework and accompanying computational methods to address this issue of constraint achievement or selection within a restricted space that will produce revised weights with reasonable properties. By combining concepts from generalized raking, ridge and lasso regression, benchmarking of small area estimates, augmentation of state-space equations, path algorithms, and data-cloning, this framework simultaneously selects constraints and provides diagnostics suggesting why a fully constrained solution is not possible. Combinatoric operations such as brute force evaluations of all possible combinations of constraints and restrictions are avoided. We demonstrate this framework by applying alternative methods to post-stratification for the National Survey on Drug Use and Health. We also discuss strategies for scaling up to even larger data sets. Computations were performed in R and code is available from the authors.


title: Acquisition of time-frequency localized mechanical properties of biofilms and single cells with high spatial resolution

abstract: Biofilms are a cluster of bacteria embedded in extracellular polymeric substances (EPS) that contain a complex composition of polysaccharides, proteins, and extracellular DNA (eDNA). Desirable mechanical properties of the biofilms are critical for their survival, propagation, and dispersal, and the response of mechanical properties to different treatment conditions also sheds light on biofilm control and eradication in vivo and on engineering surfaces. However, it is challenging yet important to interrogate mechanical behaviors of biofilms with a high spatial resolution because biofilms are very heterogeneous. Moreover, biofilms are viscoelastic, and their time-dependent mechanical behavior is difficult to capture. Herein, we developed a powerful technique that combines the high spatial resolution of the atomic force microscope (AFM) with a rigorous history-dependent viscoelastic analysis to deliver highly spatial-localized biofilm properties within a wide time-frequency window. By exploiting the use of static force spectroscopy in combination with an appropriate viscoelastic framework, we highlight the intensive amount of time-dependent information experimentally available that has been largely overlooked. It is shown that this technique provides a detailed nanorheological signature of the biofilms even at the single-cell level. We share the computational routines that would allow any user to perform the analysis from experimental raw data. The detailed localization of mechanical properties in space and in time-frequency domain provides insights on the understanding of biofilm stability, cohesiveness, dispersal, and control.


title: The Importance of Socio-Cultural Differences for Annotating and Detecting the Affective States of Students

abstract: The development of real-time affect detection models often depends upon obtaining annotated data for supervised learning by employing human experts to label the student data. One open question in annotating affective data for affect detection is whether the labelers (i.e., human experts) need to be socio-culturally similar to the students being labeled, as this impacts the cost feasibility of obtaining the labels. In this study, we investigate the following research questions: For affective state annotation, how does the socio-cultural background of human expert labelers, compared to the subjects, impact the degree of consensus and distribution of affective states obtained? Secondly, how do differences in labeler background impact the performance of affect detection models that are trained using these labels?


title: Analyzing a Maximum Principle for Finite Horizon State Constrained Problems via Parametric Examples. Part 1: Problems with Unilateral State Constraints

abstract: In the present paper, the maximum principle for finite horizon state constrained problems from the book by R. Vinter [\textit{Optimal Control}, Birkh\"auser, Boston, 2000; Theorem~9.3.1] is analyzed via parametric examples. The latter has origin in a recent paper by V.~Basco, P.~Cannarsa, and H.~Frankowska, and resembles the optimal growth problem in mathematical economics. The solution existence of these parametric examples is established by invoking Filippov's existence theorem for Mayer problems. Since the maximum principle is only a necessary condition for local optimal processes, a large amount of additional investigations is needed to obtain a comprehensive synthesis of finitely many processes suspected for being local minimizers. Our analysis not only helps to understand the principle in depth, but also serves as a sample of applying it to meaningful prototypes of economic optimal growth models. Problems with unilateral state constraints are studied in Part 1 of the paper. Problems with bilateral state constraints will be addressed in Part 2.


title: Connectedness of the Moduli Space of Genus 1 Planar Tropical Curves

abstract: Tropical geometry is a relatively recent field in mathematics created as a simplified model for certain problems in algebraic geometry. We introduce the definition of abstract and planar tropical curves as well as their properties, including combinatorial type and degree. We also talk about the moduli space, a geometric object that parameterizes all possible types of abstract or planar tropical curves subject to certain conditions. Our research focuses on the moduli spaces of planar tropical curves of genus one, arbitrary degree d and any number of marked, unbounded edges. We prove that these moduli spaces are connected.


title: Learning Pairwise Relationship for Multi-object Detection in Crowded Scenes

abstract: As the post-processing step for object detection, non-maximum suppression (GreedyNMS) is widely used in most of the detectors for many years. It is efficient and accurate for sparse scenes, but suffers an inevitable trade-off between precision and recall in crowded scenes. To overcome this drawback, we propose a Pairwise-NMS to cure GreedyNMS. Specifically, a pairwise-relationship network that is based on deep learning is learned to predict if two overlapping proposal boxes contain two objects or zero/one object, which can handle multiple overlapping objects effectively. Through neatly coupling with GreedyNMS without losing efficiency, consistent improvements have been achieved in heavily occluded datasets including MOT15, TUD-Crossing and PETS. In addition, Pairwise-NMS can be integrated into any learning based detectors (Both of Faster-RCNN and DPM detectors are tested in this paper), thus building a bridge between GreedyNMS and end-to-end learning detectors.


title: Integrating multi-source block-wise missing data in model selection

abstract: For multi-source data, blocks of variable information from certain sources are likely missing. Existing methods for handling missing data do not take structures of block-wise missing data into consideration. In this paper, we propose a Multiple Block-wise Imputation (MBI) approach, which incorporates imputations based on both complete and incomplete observations. Specifically, for a given missing pattern group, the imputations in MBI incorporate more samples from groups with fewer observed variables in addition to the group with complete observations. We propose to construct estimating equations based on all available information, and optimally integrate informative estimating functions to achieve efficient estimators. We show that the proposed method has estimation and model selection consistency under both fixed-dimensional and high-dimensional settings. Moreover, the proposed estimator is asymptotically more efficient than the estimator based on a single imputation from complete observations only. In addition, the proposed method is not restricted to missing completely at random. Numerical studies and ADNI data application confirm that the proposed method outperforms existing variable selection methods under various missing mechanisms.


title: 3D Human Pose Machines with Self-supervised Learning

abstract: Driven by recent computer vision and robotic applications, recovering 3D human poses has become increasingly important and attracted growing interests. In fact, completing this task is quite challenging due to the diverse appearances, viewpoints, occlusions and inherently geometric ambiguities inside monocular images. Most of the existing methods focus on designing some elaborate priors /constraints to directly regress 3D human poses based on the corresponding 2D human pose-aware features or 2D pose predictions. However, due to the insufficient 3D pose data for training and the domain gap between 2D space and 3D space, these methods have limited scalabilities for all practical scenarios (e.g., outdoor scene). Attempt to address this issue, this paper proposes a simple yet effective self-supervised correction mechanism to learn all intrinsic structures of human poses from abundant images. Specifically, the proposed mechanism involves two dual learning tasks, i.e., the 2D-to-3D pose transformation and 3D-to-2D pose projection, to serve as a bridge between 3D and 2D human poses in a type of "free" self-supervision for accurate 3D human pose estimation. The 2D-to-3D pose implies to sequentially regress intermediate 3D poses by transforming the pose representation from the 2D domain to the 3D domain under the sequence-dependent temporal context, while the 3D-to-2D pose projection contributes to refining the intermediate 3D poses by maintaining geometric consistency between the 2D projections of 3D poses and the estimated 2D poses. We further apply our self-supervised correction mechanism to develop a 3D human pose machine, which jointly integrates the 2D spatial relationship, temporal smoothness of predictions and 3D geometric knowledge. Extensive evaluations demonstrate the superior performance and efficiency of our framework over all the compared competing methods.


title: On continuous weaving fusion frames in Hilbert spaces

abstract: In this note, we first introduce the notation of weaving cfusion frames in separable Hilbert spaces. After reviewing the conditions for maintaining the weaving c-fusion frames under the bounded linear operator and also, removing vectors from these frames, we will present a necessarily and sufficient condition about c-woven and c-fusion woven. Finally, perturbation of these frames will be introduced.


title: The LFV decays of Z boson in Minimal R-symmetric Supersymmetric Standard Model

abstract: A future $Z$-factory will offer the possibility to study rare $Z$ decays $Z\rightarrow l_1l_2$, as those leading to Lepton Flavor Violation final states. In this work, by taking account of the constraints from radiative two body decays $l_2\rightarrow l_1\gamma$, we investigate the Lepton Flavor Violation decays $Z\rightarrow l_1l_2$ in the framework of Minimal R-symmetric Supersymmetric Standard Model with two benchmark points from already existing literatures. The flavor violating off-diagonal entries $\delta^{12}$, $\delta^{13}$ and $\delta^{23}$ are constrained by the current experimental bounds of $l_2\rightarrow l_1\gamma$. Considering recent experimental constraints, we also investigate Br($Z\rightarrow l_1l_2$) as a function of $M_D^W$. The numerical results show that the theoretical prediction of Br($Z\rightarrow l_1l_2$) in MRSSM are several orders of magnitude below the current experimental bounds. The Lepton Flavor Violation decays $Z\rightarrow e\tau$ and $Z\rightarrow \mu\tau$ may be promising to be observed in future experiment.


title: Operators in the Cowen-Douglas class and related topics

abstract: Linear spaces with an Euclidean metric are ubiquitous in mathematics, arising both from quadratic forms and inner products. Operators on such spaces also occur naturally. In recent years, the study of multivariate operator theory has made substantial progress. Although the study of self-adjoint operators goes back a few decades, the non self-adjoint theory has developed at a slower pace. While several approaches to this topic has been developed, the one that has been most fruitful is clearly the study of Hilbert spaces that are modules over natural function algebras like $\mathcal A({\Omega})$, where $\Omega \subseteq \mathbb C^m$ is a bounded domain, consisting of complex valued functions which are holomorphic on some open set $U$ containing $\overline{\Omega}$, the closure of $\Omega$. The book, ''Hilbert Modules over function algebra'', R. G. Douglas and V. I. Paulsen showed how to recast many of the familiar theorems of operator theory in the language of Hilbert modules. The book, ''Spectral decomposition of analytic sheaves'', J. Eschmeier and M. Putinar and the book, ''Analytic Hilbert modules'', X. Chen and K. Guo, provide an account of the achievements from the recent past. The impetus for much of what is described below comes from the interplay of operator theory with other areas of mathematics like complex geometry and representation theory of locally compact groups.


title: ALiPy: Active Learning in Python

abstract: Supervised machine learning methods usually require a large set of labeled examples for model training. However, in many real applications, there are plentiful unlabeled data but limited labeled data; and the acquisition of labels is costly. Active learning (AL) reduces the labeling cost by iteratively selecting the most valuable data to query their labels from the annotator. This article introduces a Python toobox ALiPy for active learning. ALiPy provides a module based implementation of active learning framework, which allows users to conveniently evaluate, compare and analyze the performance of active learning methods. In the toolbox, multiple options are available for each component of the learning framework, including data process, active selection, label query, results visualization, etc. In addition to the implementations of more than 20 state-of-the-art active learning algorithms, ALiPy also supports users to easily configure and implement their own approaches under different active learning settings, such as AL for multi-label data, AL with noisy annotators, AL with different costs and so on. The toolbox is well-documented and open-source on Github, and can be easily installed through PyPI.


title: Continuous k-frames and their duals

abstract: k-frames were recently introduced by Gavruta in Hilbert spaces to study atomic systems with respect to a bounded linear operator. A continuous frame is a family of vectors in Hilbert space which allows reproductions of arbitrary elements by continuous super positions. We construct a continuous k-frame, so called ck-frame along with an atomic system for this version of frames. Also we introduce a new method for obtaining the dual of a ck-frame and we prove some new results about it.


title: Byzantine Consensus under Local Broadcast Model: Tight Sufficient Condition

abstract: In this work we consider Byzantine Consensus on undirected communication graphs under the local broadcast model. In the classical point-to-point communication model the messages exchanged between two nodes $u, v$ on an edge $uv$ of $G$ are private. This allows a faulty node to send conflicting information to its different neighbours, a property called equivocation. In contrast, in the local broadcast communication model considered here, a message sent by node $u$ is received identically by all of its neighbours. This restriction to broadcast messages provides non-equivocation even for faulty nodes. In prior results [10, 11] it was shown that in the local broadcast model the communication graph must be $(\lfloor 3f/2 \rfloor +1)$-connected and have degree at least $2f$ to achieve Byzantine Consensus. In this work we show that this network condition is tight.


title: The LFV decays $B^0_{d,s}\rightarrow e\mu(e\tau,\mu\tau)$ with one neutral singlet scalar

abstract: Taking account of the constraint from radiative two body decays $l_i\rightarrow l_j\gamma$, we investigate the Lepton Flavor Violation decays $B^0_q\rightarrow\bar{l_l}l_k$ in the framework of the minimal extension of the Standard Model with one neutral singlet scalar. The couplings $C_{e\mu}$, $C_{e\tau}$ and $C_{\mu\tau}$ between the different generation leptons and scalar $S^0$ are constrained by the current bounds of $l_i\rightarrow l_j\gamma$. The numerical results show that the theoretical prediction of $B^0_q\rightarrow\bar{l_l}l_k$ strongly depend on the couplings $C_{qb}$ ($q=d\;or\; s$) between down type quarks and new scalar. The contributions from couplings $C_{uc}$, $C_{ut}$ and $C_{ct}$ between up type quark and new scalar are less dominant.


title: One loop correction to $Z\rightarrow \nu\nu$ in the Minimal R-symmetric Supersymmetric Standard Model

abstract: We analyze the one loop correction to $Z\rightarrow \nu\nu$ decay in framework of Minimal R-symmetric Supersymmetric Standard Model(MRSSM) in detail with normal and inverse neutrino mass orderings, as a function of $\tan\beta$, Dirac mass parameters $M_D^W$ and $\mu_u(\mu_d)$, slepton mass $m_l$ that parameterize the mass matrices. The numerical results indicate that the branching ratio for $Z\rightarrow \nu\nu$ decay is compatible with the experimental measurement and the SM expectation at $2\sigma$ level. For inverse neutrino mass ordering, the prediction exceeds the SM expectation at $1\sigma$ level. The prediction on $Br(Z\rightarrow \nu\nu)$ increases proportionally to $\tan\beta$ and inversely proportionally to $m_l$. For normal neutrino mass ordering, the peak value of the prediction on $Br(Z\rightarrow \nu\nu)$ exceeds the SM expectation at $1\sigma$ level.


title: Reverse Engineering the Raspberry Pi Camera V2: A study of Pixel Non-Uniformity using a Scanning Electron Microscope

abstract: In this paper we reverse engineer the Sony IMX219PQ image sensor, otherwise known as the Raspberry Pi Camera v2.0. We provide a visual reference for pixel non-uniformity by analysing variations in transistor length, microlens optic system and in the photodiode. We use these measurements to demonstrate irregularities at the microscopic level and link this to the signal variation measured as pixel non-uniformity used for unique identification of discrete image sensors.


title: ECGadv: Generating Adversarial Electrocardiogram to Misguide Arrhythmia Classification System

abstract: Deep neural networks (DNNs)-powered Electrocardiogram (ECG) diagnosis systems recently achieve promising progress to take over tedious examinations by cardiologists. However, their vulnerability to adversarial attacks still lack comprehensive investigation. The existing attacks in image domain could not be directly applicable due to the distinct properties of ECGs in visualization and dynamic properties. Thus, this paper takes a step to thoroughly explore adversarial attacks on the DNN-powered ECG diagnosis system. We analyze the properties of ECGs to design effective attacks schemes under two attacks models respectively. Our results demonstrate the blind spots of DNN-powered diagnosis systems under adversarial attacks, which calls attention to adequate countermeasures.


title: On the PIR capacity of MSR codes

abstract: We consider the problem of private information retrieval (PIR) from multiple storage nodes when the underlying database is encoded using regenerating codes, i.e., the database has the ability to recover from individual node failures. We seek optimality in terms of download rates when a user wants to retrieve a message privately (in an information-theoretic sense) from the database. We focus mostly on the more extensively explored Minimum Storage Regenerating (MSR) codes and formally state and prove the PIR capacity of such codes. We also give a simple scheme to achieve this capacity which is based on the similarly motivated construction for MDS codes. Furthermore, we give a PIR scheme for a specific class of Minimum Bandwidth Regenerating (MBR) codes that achieves the same rate as its minimum storage counterpart.


title: Universal Arduino-based experimenting system to support teaching of natural sciences

abstract: The rapid evolution of intelligent electronic devices makes information technology, computer science and electronics strongly related to the teaching of natural sciences. Today almost everybody has a smart phone that can convert light, temperature, movement, sound to numbers, therefore all these can be processed, analysed, displayed, stored, shared by software applications. The fundamental question is how education can follow this knowledge and how can education take its advantages. Components and methods of modern technology are available for education also, teachers and students can play with parts and tools which were previously used only by engineers. A good example is the very popular Arduino board which is practically an industrial microcontroller whose pins are wired to easy-to-use connectors on a printed circuit board. In this paper we show a universal system which we have developed for the Arduino platform to support experimenting and understanding of the most fundamental principles of the operation of modern devices. We show our related educational concept and discuss the most important features of the system. Open source hardware and software are available and we provide a number of video tutorials as well.


title: How accurate is an Arduino Ohmmeter?

abstract: The Arduino platform is widely used in education of physics to perform a number of different measurements. Teachers and students can build their own instruments using various sensors, the analogue-to-digital converter of the Arduino board and code to calculate and display the result. In several cases this can mean incautious reproduction of what can be found on the Internet and an in-depth understanding can be missing. Here we thoroughly analyse a frequently used resistance measurement method and show demonstration experiments as well to make it clear.


title: Kelvin modes of a skyrmion line in chiral magnets and the associated magnon transport

abstract: Magnetic skyrmions in bulk crystals are line-like topologically protected spin textures. They allow for the propagation of magnons along the skyrmion line but are localized inside the skyrmion line. Analogous to the vortex line, these propagating modes are the Kelvin modes of a skyrmion line. In crystals without an inversion center, it is known that the magnon dispersion in the ferromagnetic state is asymmetric in the wavevector. It is natural to expect that the dispersion of the Kelvin modes is also asymmetric with respect to the wavevector. We study the Kelvin modes of a skyrmion line in the ferromagnetic background. In contrast, we find that the lowest Kelvin mode is symmetric in the wavevector in the low energy region despite the inversion symmetry breaking. Other Kelvin modes below the magnon continuum are asymmetric, and most of them have a positive group velocity. Our results suggest that a skyrmion line can function as a one-way waveguide for magnons.


title: Radii problems for the generalized Mittag-Leffler functions

abstract: In this paper our aim is to investigate the radii of $\eta-$uniformly convexity, $\alpha-$convexity, $\eta-$parabolic starlikeness and strong starlikeness of order $\rho$ of the generalized Mittag-Leffler function for three different kinds of normalization by using their Hadamard factorization in such a way that the resulting functions are analytic. The basic tool of this study is the Mittag-Leffler function in series. Also we have shown that the obtained radii are the smallest positive roots of some functional equations.


title: Boundary-Aware Network for Fast and High-Accuracy Portrait Segmentation

abstract: Compared with other semantic segmentation tasks, portrait segmentation requires both higher precision and faster inference speed. However, this problem has not been well studied in previous works. In this paper, we propose a lightweight network architecture, called Boundary-Aware Network (BANet) which selectively extracts detail information in boundary area to make high-quality segmentation output with real-time( >25FPS) speed. In addition, we design a new loss function called refine loss which supervises the network with image level gradient information. Our model is able to produce finer segmentation results which has richer details than annotations.


title: Osaka Feedback Model: Isolated Disk Galaxy Simulations

abstract: We study various implementations of supernova feedback model and present the results of our `Osaka feedback model' using isolated galaxy simulations performed by the smoothed particle hydrodynamics (SPH) code {\small GADGET-3}. Our model is a modified version of Stinson et al.'s work, and we newly add the momentum kick for SN feedback rather than only thermal feedback. We incorporate the physical properties of SN remnants from the results of Chevalier and McKee \& Ostriker, such as the effective radius of SN bubble and the remnant life-time, in the form of Sedov-Taylor (ST)-like solutions with the effect of radiative cooling. Our model utilizes the local, physical parameters such as density and temperature of the ISM rather than galactic or halo properties to determine the galactic wind velocity or mass-loading factor. The Osaka model succeeds in self-regulating star formation, and naturally produces galactic outflow with variable velocities depending on the local environment and available SN energy as a function of time.An important addition to our previous work by Aoyama et al. is the implementation of the {\small CELib} chemistry library which allows us to deal with the time-dependent input of energy and metal yields for type Ia \& II supernovae (SNe) and asymptotic giant branch (AGB) stars. As initial tests of our model, we apply it to isolated galaxy simulations, and examine various galactic properties and compare with observational data including metal abundances.


title: Simple juntas for shifted families

abstract: We say that a family $\mathcal F$ of $k$-element sets is a {\it $j$-junta} if there is a set $J$ of size $j$ such that, for any $F$, its presence in $\mathcal F$ depends on its intersection with $J$ only. Approximating arbitrary families by $j$-juntas with small $j$ is a recent powerful technique in extremal set theory. The weak point of all known junta approximation results is that they work in the range $n>Ck$, where $C$ is an extremely fast growing function of the input parameters, such as the quality of approximation or the number of families we simultaneously approximate. We say that a family $\mathcal F$ is {\it shifted} if for any $F=\{x_1,\ldots, x_k\}\in \mathcal F$ and any $G =\{y_1,\ldots, y_k\}$ such that $y_i\le x_i$, we have $G\in \mathcal F$. For many extremal set theory problems, including the Erd\H os Matching Conjecture, or the Complete $t$-Intersection Theorem, it is sufficient to deal with shifted families only. In this paper, we present very general approximation by juntas results for shifted families with explicit (and essentially linear) dependency on the input parameters. The results are best possible up to some constant factors. Moreover, they give meaningful statements for almost all range of values of $n$. The proofs are shorter than the proofs of the previous approximation by juntas results and are completely self-contained. As an application of our junta approximation, we give a nearly-linear bound for the multi-family version of the Erd\H os Matching Conjecture. More precisely, we prove the following result. Let $n\ge 12sk\log(e^2s)$ and suppose that the families $\mathcal F_1,\ldots, \mathcal F_s\subset {[n]\choose k}$ do not contain $F_1\in\mathcal F_1,\ldots, F_s\in \mathcal F_s$ such that $F_i$'s are pairwise disjoint. Then $\min_{i}|\mathcal F_i|\le {n\choose k}-{n-s+1\choose k}.$


title: Logical Reasoning for Revealing the Critical Temperature through Deep Learning of Configuration Ensemble of Statistical Systems

abstract: Recently, there have been many works on the deep learning of statistical ensembles to determine the critical temperature of a possible phase transition. We analyze the detailed structure of an optimized deep learning machine and prove the basic equalities among the optimized machine parameters and the physical quantities of the statistical system. According to these equalities, we conclude that the bias parameters of the final full connection layer record the free energy of the statistical system as a function of temperature. We confirm these equalities in one- and two-dimensional Ising spin models and actually demonstrate that the deep learning machine reveals the critical temperature of the phase transition through the second difference of bias parameters, which is equivalent to the specific heat. Our results disprove the previous works claiming that the weight parameters of the full connection might play a role of the order parameter such as the spin expectation.


title: Monitoring breast cancer response to neoadjuvant chemotherapy with ultrasound signal statistics and integrated backscatter

abstract: Monitoring Neoadjuvant chemotherapy (NAC) effects is necessary to capture resistant patients and stop or change treatment. The aim of this study was to assess the tumor response at an early stage, after the first doses of the NAC, based on the variability of the backscattered ultrasound energy, and backscatter statistics. The backscatter statistics has not previously been used to monitor NAC effects. The B-mode ultrasound images and raw radio frequency data from breast tumors were obtained using an ultrasound scanner before chemotherapy and 1 week after each NAC cycle. Twenty-four malignant breast cancers, qualified for neoadjuvant treatment before surgery, were included in the study. The shape parameter of the homodyned K distribution and integrated backscatter, along with the tumor size in the longest dimension, were determined based on ultrasound data and used as markers for NAC response. Cancer tumors were assigned to responding and non-responding groups, according to histopathological evaluation, which was a reference in assessing the utility of markers. Statistical analysis was performed to rate the ability of markers to predict NAC response based on data obtained after subsequent therapeutic doses. Statistically significant differences between groups were obtained after 2, 3, 4, and 5 doses of NAC for quantitative ultrasound markers and after 5 doses for the assessment based on maximum tumor dimension. After the second and third NAC courses the marker, which was a linear combination of both quantitative ultrasound parameters, was characterized by an AUC of 0.82 and 0.91, respectively. The introduction of statistical parameters of ultrasonic backscatter to monitor the effects of chemotherapy can increase the effectiveness of monitoring and contribute to a better personalization of NAC therapy.


title: Surface impedance measurements on Nb$_{3}$Sn at high magnetic fields

abstract: Nb$_{3}$Sn is a superconductor of great relevance for perspective RF applications. We present for the first time surface impedance $Z_s$ measurements at 15 GHz and low RF field amplitude on Nb$_{3}$Sn in high magnetic fields up to 12 T, with the aim of increasing the knowledge of Nb$_{3}$Sn behavior in such conditions. $Z_s$ is a fundamental material parameter that directly gives useful information about the dissipative and reactive phenomena when the superconductor is subjected to high-frequency excitations. Therefore, we present an analysis of the measured $Z_s$ with the aim of extracting interesting data about pinning in Nb$_{3}$Sn at high frequencies. From $Z_s$ we extract the vortex motion complex resistivity to obtain the $r$-parameter and the depinning frequency $\nu_p$ in high magnetic fields. The comparison of the results with the literature shows that the measured $\nu_p$ on bulk Nb$_{3}$Sn is several times greater than that of pure Nb. This demonstrates how Nb$_{3}$Sn can be a good candidate for RF technological applications, also in high magnetic fields.


title: On the structure of locally potentially equivalent Galois representations

abstract: Suppose $\rho_1, \rho_2$ are two $\ell$-adic Galois representations of the absolute Galois group of a number field, such that the algebraic monodromy group of one of the representations is connected and the representations are locally potentially equivalent at a set of places of positive upper density. We classify such pairs of representations and show that up to twisting by some representation, it is given by a pair of representations one of which is trivial and the other abelian. Consequently, assuming that the first representation has connected algebraic monodromy group, we obtain that the representations are potentially equivalent, provided one of the following conditions hold: (a) the first representation is absolutely irreducible; (b) the ranks of the algebraic monodromy groups are equal; (c) the algebraic monodromy group of the second representation is also connected and (d) the commutant of the image of the second representation remains the same upon restriction to subgroups of finite index of the Galois group.


title: Mastering Panel 'Metrics: Causal Impact of Democracy on Growth

abstract: The relationship between democracy and economic growth is of long-standing interest. We revisit the panel data analysis of this relationship by Acemoglu, Naidu, Restrepo and Robinson (forthcoming) using state of the art econometric methods. We argue that this and lots of other panel data settings in economics are in fact high-dimensional, resulting in principal estimators -- the fixed effects (FE) and Arellano-Bond (AB) estimators -- to be biased to the degree that invalidates statistical inference. We can however remove these biases by using simple analytical and sample-splitting methods, and thereby restore valid statistical inference. We find that the debiased FE and AB estimators produce substantially higher estimates of the long-run effect of democracy on growth, providing even stronger support for the key hypothesis in Acemoglu, Naidu, Restrepo and Robinson (forthcoming). Given the ubiquitous nature of panel data, we conclude that the use of debiased panel data estimators should substantially improve the quality of empirical inference in economics.


title: Strong Evidence that the Galactic Bulge is Shining in Gamma Rays

abstract: There is growing evidence that the Galactic Center Excess identified in the $\textit{Fermi}$-LAT gamma-ray data arises from a population of faint astrophysical sources. We provide compelling supporting evidence by showing that the morphology of the excess traces the stellar over-density of the Galactic bulge. By adopting a template of the bulge stars obtained from a triaxial 3D fit to the diffuse near-infrared emission, we show that it is detected at high significance. The significance deteriorates when either the position or the orientation of the template is artificially shifted, supporting the correlation of the gamma-ray data with the Galactic bulge. In deriving these results, we have used more sophisticated templates at low-latitudes for the $\textit{Fermi}$ bubbles compared to previous work and the three-dimensional Inverse Compton (IC) maps recently released by the ${\tt GALPROP}$ team. Our results provide strong constraints on Millisecond Pulsar (MSP) formation scenarios proposed to explain the excess. We find that an $\textit{admixture formation}$ scenario, in which some of the relevant binaries are $\textit{primordial}$ and the rest are formed $\textit{dynamically}$, is preferred over a primordial-only formation scenario at $7.6\sigma$ confidence level. Our detailed morphological analysis also disfavors models of the disrupted globular clusters scenario that predict a spherically symmetric distribution of MSPs in the Galactic bulge. For the first time, we report evidence of a high energy tail in the nuclear bulge spectrum that could be the result of IC emission from electrons and positrons injected by a population of MSPs and star formation activity from the same site.

categories: astro-ph.HE astro-ph.CO astro-ph.GA hep-ph


title: Anderson localization and delocalization of massless two-dimensional Dirac electrons in random one-dimensional scalar and vector potentials

abstract: We study Anderson localization of massless Dirac electrons in two dimensions in one-dimensional random scalar and vector potentials theoretically for two different cases, in which the scalar and vector potentials are either uncorrelated or correlated. From the Dirac equation, we deduce the effective wave impedance, using which we derive the condition for total transmission and those for delocalization in our random models analytically. Based on the invariant imbedding theory, we also develop a numerical method to calculate the localization length exactly for arbitrary strengths of disorder. In addition, we derive analytical expressions for the localization length, which are extremely accurate in the weak and strong disorder limits. In the presence of both scalar and vector potentials, the conditions for total transmission and complete delocalization are generalized from the usual Klein tunneling case. We find that the incident angles at which electron waves are either completely transmitted or delocalized can be tuned to arbitrary values. When the strength of scalar potential disorder increases to infinity, the localization length also increases to infinity, both in uncorrelated and correlated cases. The detailed dependencies of the localization length on incident angle, disorder strength and energy are elucidated and the discrepancies with previous studies and some new results are discussed. All the results are explained intuitively using the concept of wave impedance.


title: On Kuznetsov-Bykovskii's formula of counting prime geodesics

abstract: We generalize a formula on the counting of prime geodesics, due to Kuznetsov-Bykovskii, used in the work of Soundararajan-Young on the prime geodesic theorem. The method works over any number field and for any congruence subgroup. We give explicit computation in the cases of principal and Hecke subgroups.


title: Fr\'echet-valued formal power series

abstract: Let $A$ be a non-projectively-pluripolar set in a Fr\'{e}chet space $E.$ We give sufficient conditions to ensure the convergence on some zero-neighbourhood in $E$ of a (sequence of) formal power series of Fr\'{e}chet-valued continuous homogeneous polynomials provided that the convergence holds at a zero-neighbourhood of each complex line $\ell_a := \C a$ for every $a \in A.$


title: A posteriori error estimates for hypersingular integral equation on spheres with spherical splines

abstract: A posteriori residual and hierarchical upper bounds for the error estimates were proved when solving the hypersingular integral equation on the unit sphere by using the Galerkin method with spherical splines. Based on these a posteriori error estimates, adaptive mesh refining procedures are used to reduce complexity and computational cost of the discrete problems. Numerical experiments illustrate our theoretical results.


title: A proof of the $C^{p^\prime}$-regularity conjecture in the plane

abstract: We establish a new oscillation estimate for solutions of nonlinear partial differential equations of elliptic, degenerate type. This new tool yields a precise control on the growth rate of solutions near their set of critical points, where ellipticity degenerates. As a consequence, we are able to prove the planar counterpart of the longstanding conjecture that solutions of the degenerate $p$-Poisson equation with a bounded source are locally of class $C^{p^\prime}=C^{1,\frac{1}{p-1}}$; this regularity is optimal.


title: Direct calculation of the strong Goos-H\"{a}nchen effect of a Gaussian light beam due to the excitation of surface plasmon polaritons in the Otto configuration

abstract: We study theoretically the influence of the surface plasmon excitation on the Goos-H\"{a}nchen lateral shift of a $p$-polarized Gaussian beam incident obliquely on a dielectric-metal bilayer in the Otto configuration. We find that the lateral shift depends sensitively on the thickness of the metal layer and the width of the incident beam, as well as on the incident angle. Near the incident angle at which surface plasmons are excited, the lateral shift changes from large negative values to large positive values as the thickness of the metal layer increases through a critical value. For wide incident beams, the maximal forward and backward lateral shifts can be as large as several hundred times of the wavelength. As the width of the incident Gaussian beam decreases, the magnitude of the lateral shift decreases rapidly, but the ratio of the width of the reflected beam to that of the incident beam, which measures the degree of the deformation of the reflected beam profile, increases. In all cases considered, we find that the reflected beam is split into two parts. We also find that the lateral shift of the transmitted beam is always positive and very weak.


title: Predicting Diffusion Reach Probabilities via Representation Learning on Social Networks

abstract: Diffusion reach probability between two nodes on a network is defined as the probability of a cascade originating from one node reaching to another node. An infinite number of cascades would enable calculation of true diffusion reach probabilities between any two nodes. However, there exists only a finite number of cascades and one usually has access only to a small portion of all available cascades. In this work, we addressed the problem of estimating diffusion reach probabilities given only a limited number of cascades and partial information about underlying network structure. Our proposed strategy employs node representation learning to generate and feed node embeddings into machine learning algorithms to create models that predict diffusion reach probabilities. We provide experimental analysis using synthetically generated cascades on two real-world social networks. Results show that proposed method is superior to using values calculated from available cascades when the portion of cascades is small.


title: On the Cauchy problem for stochastic integro-differential equations with radially O-regularly varying Levy measure

abstract: Parabolic integro-differential nondegenerate Cauchy problem is considered in the scale of L_{p} spaces of functions whose regularity is defined by a Levy measure with O-regulary varying radial profile. Existence and uniqueness of a solution is proved by deriving apriori estimates. Some probability density function estimates of the associated Levy process are used as well.


title: John--Nirenberg--Campanato Spaces

abstract: Let $p\in (1,\infty)$, $q\in[1,\infty)$, $\alpha\in [0,\infty)$ and $s$ be a non-negative integer. In this article, the authors introduce the John--Nirenberg-Campanato space $JN_{(p,q,s)_\alpha}(\mathcal{X})$, where $\mathcal{X}$ is ${\mathbb R}^n$ or any closed cube $Q_0\subsetneqq{\mathbb R}^n$, which when $\alpha=0$ and $s=0$ coincides with the $JN_p$-space introduced by F. John and L. Nirenberg in the sense of equivalent norms. The authors then give the predual space of $JN_{(p,q,s)_\alpha}(\mathcal{X})$ and a John-Nirenberg type inequality of John--Nirenberg-Campanato spaces. Moreover, the authors prove that the classical Campanato space serves as a limit space of $JN_{(p,q,s)_\alpha}(\mathcal{X})$ when $p\to \infty$.


title: Study of Lithium Rich Giants with the GALAH Spectroscopic Survey

abstract: In this article, we speculate on the possible mechanisms for Li enhancement origin in RGB stars based on a large data set of around 340,299 stars collected from the GALAH survey combined with the Gaia astrometry. Data has 51,982 low mass (M$\leq$ 2M$_{\odot}$) RGB stars with reliable atmospheric parameters. The data set shows a well populated RGB with well-defined luminosity bump and red clump with significant number of stars at each of these two key phases. We found 335 new Li-rich RGB stars with Li abundance, A(Li) $\geq$ 1.80 $\pm$ 0.14 dex, of which 20 are super Li-rich with A(Li) $\geq$ 3.20~dex. Most of them appear to be in the red clump region which, when combined with stellar evolutionary timescales on RGB, indicates that the Li enhancement origin may lie at RGB tip during He-flash rather than by external source of merging of sub-stellar objects or during luminosity bump evolution. Kinematic properties of sample stars suggest that Li-rich giants are relatively more prevalent among giants of thin disk compared to thick disk and halo.


title: Hypersurfaces with linear type singular loci

abstract: In this paper, necessary and sufficient criteria for the Jacobian ideal of a reduced hypersurface with isolated singularity to be of linear type, are presented. We prove that the gradient ideal of a reduced projective plane curve with simple singularities ($\mathrm{ADE}$) is of linear type. We show that any reduced projective quartic curve is of gradient linear type.


title: On the projections of the multifractal packing dimension for q>1

abstract: The aim of this article is to study the behaviour of the multifractal packing function $B_\mu(q)$ under projections in Euclidean space for $q>1$. We show that $B_\mu(q)$ is preserved under almost every orthogonal projection. As an application, we study the multifractal analysis of the projections of a measure. In particular, we obtain general results for the multifractal analysis of the orthogonal projections on $m$-dimensional linear subspaces of a measure $\mu$ satisfying the multifractal formalism.


title: Edge states at nematic domain walls in FeSe films

abstract: Quantum spin Hall (QSH) effect is an intriguing phenomenon arising from the helical edge states in two-dimensional topological insulators. We use molecular beam epitaxy (MBE) to prepare FeSe films with atomically sharp nematic domain boundaries, where tensile strains, nematicity suppression and topological band inversion are simultaneously achieved. Using scanning tunneling microscopy (STM), we observe edge states at the Fermi level that spatially distribute as two distinct strips in the vicinity of the domain boundaries. At the endpoint of the boundaries, a bound state at the Fermi level is further observed. The topological origin of the edge states is supported by density functional theory calculations. Our findings not only demonstrate a candidate for QSH states, but also provide a new pathway to realize topological superconductivity in a single-component film.


title: Development of the Micro Pixel Chamber with resistive electrodes

abstract: We developed a novel design of a Micro Pixel Chamber ($\mu$-PIC) with resistive electrodes for a charged-particle-tracking detector in high-rate applications. Diamond-Like Carbon (DLC) thin film is used for the cathodes. The resistivity can be controlled flexibly ($\mathrm{10^{5-7}k\Omega/sq.}$) at high uniformity. The fabrication-process was greatly improved and the resistive $\mu$-PIC could be operated at 10$\times$10 $\mathrm{cm^2}$. Resistors for the HV bias and capacitors for the AC coupling were completely removed by applying PCB and carbon-sputtering techniques, and the resistive $\mu$-PIC became a very compact detector. The performances of our new resistive $\mu$-PIC were measured in various ways. Consequently, it was possible to attain high gas gains ($\mathrm{> 10^{4}}$), high detection efficiency, and position resolution exceeding 100 $\mu$m. The spark current was suppressed, and the new resistive $\mu$-PIC was operated stably under fast-neutrons irradiation. These features offer solutions for a charged-particle-tracking detector in future high-rate applications.


title: Derivation of Correlation Dimension from Spatial Autocorrelation Functions

abstract: Spatial autocorrelation coefficients such as Moran's index proved to be an eigenvalue of the spatial correlation matrixes. An eigenvalue represents a kind of characteristic length for quantitative analysis. However, if a spatial correlation is based on self-organized evolution, complex structure, and the distributions without characteristic scale, the eigenvalue will be ineffective. In this case, the single Moran index cannot lead to reliable statistic inferences. This paper is devoted to finding advisable approach to measure spatial autocorrelation for the scale-free processes of complex systems by means of mathematical reasoning and empirical analysis. Based on relative step function as spatial contiguity function, a series of ordered spatial autocorrelation coefficients are converted into the corresponding spatial autocorrelation functions. Then the mathematical relation between spatial correlation dimension and spatial autocorrelation functions is derived by decomposition of spatial autocorrelation functions. As results, a set of useful mathematical models are constructed for spatial analysis. Using these models, we can utilize spatial correlation dimension to make simple spatial autocorrelation analysis, and use spatial autocorrelation functions to make complex spatial autocorrelation analysis for geographical phenomena. This study reveals the inherent association of fractal patterns with spatial autocorrelation processes in nature and society. The work may inspire new ideas of spatial modeling and analysis for complex geographical systems.


title: Enhancing Explainability of Neural Networks through Architecture Constraints

abstract: Prediction accuracy and model explainability are the two most important objectives when developing machine learning algorithms to solve real-world problems. The neural networks are known to possess good prediction performance, but lack of sufficient model interpretability. In this paper, we propose to enhance the explainability of neural networks through the following architecture constraints: a) sparse additive subnetworks; b) projection pursuit with orthogonality constraint; and c) smooth function approximation. It leads to an explainable neural network (xNN) with the superior balance between prediction performance and model interpretability. We derive the necessary and sufficient identifiability conditions for the proposed xNN model. The multiple parameters are simultaneously estimated by a modified mini-batch gradient descent method based on the backpropagation algorithm for calculating the derivatives and the Cayley transform for preserving the projection orthogonality. Through simulation study under six different scenarios, we compare the proposed method to several benchmarks including least absolute shrinkage and selection operator, support vector machine, random forest, extreme learning machine, and multi-layer perceptron. It is shown that the proposed xNN model keeps the flexibility of pursuing high prediction accuracy while attaining improved interpretability. Finally, a real data example is employed as a showcase application.


title: Operator splitting schemes for the two-asset Merton jump-diffusion model

abstract: This paper deals with the numerical solution of the two-dimensional time-dependent Merton partial integro-differential equation (PIDE) for the values of rainbow options under the two-asset Merton jump-diffusion model. Key features of this well-known equation are a two-dimensional nonlocal integral part and a mixed spatial derivative term. For its efficient and stable numerical solution, we study seven recent and novel operator splitting schemes of the implicit-explicit (IMEX) and the alternating direction implicit (ADI) kind. Here the integral part is always conveniently treated in an explicit fashion. The convergence behaviour and the relative performance of the seven schemes are investigated in ample numerical experiments for both European put-on-the-min and put-on-the-average options.


title: Uncertainty in marine weather routing

abstract: Weather routing methods are essential for planning routes for commercial shipping and recreational craft. This paper provides a methodology for quantifying the significance of numerical error and performance model uncertainty on the predictions returned from a weather routing algorithm. The numerical error of the routing algorithm is estimated by solving the optimum path over different discretizations of the environment. The uncertainty associated with the performance model is linearly varied in order to quantify its significance. The methodology is applied to a sailing craft routing problem: the prediction of the voyaging time for an ethnographic voyaging canoe across long distance voyages in Polynesia. We find that the average numerical error is $0.396\%$, corresponding to $1.05$ hours for an average voyage length of $266.40$ hours. An uncertainty level of $2.5 \%$ in the performance model is seen to correspond to a standard deviation of $\pm 2.41-3.08\%$ of the voyaging time. These results illustrate the significance of considering the influence of numerical error and performance uncertainty when performing a weather routing study.


title: Diophantine equations coming from binomial near-collisions

abstract: In this paper we solve the Diophantine equation $\binom{m}{l}-\binom{n}{k}=d$ (where m,n are positive integers unknowns) when (k,l)=(3,6) for various values of d and when (k,l)=(8,2) and d=1. As a byproduct of our results we will obtain that (k,l)-near collisions with difference 1 do not exist if (k,l)=(6,3), (3,6), (8,2) thus establishing a conjecture stated in the article "Binomial collisions and near collisions" published by A. Blokhuis, A. Brouwer and B. de Weger in Integers 17 (2017), A64.


title: Summarization and Visualization of Large Volumes of Broadcast Video Data

abstract: Over the past few years, there has been an astounding growth in the number of news channels as well as the amount of broadcast news video data. As a result, it is imperative that automated methods need to be developed in order to effectively summarize and store this voluminous data. Format detection of news videos plays an important role in news video analysis. Our problem involves building a robust and versatile news format detector, which identifies the different band elements in a news frame. Probabilistic progressive Hough transform has been used for the detection of band edges. The detected bands are classified as natural images, computer generated graphics (non-text) and text bands. A contrast based text detector has been used to identify the text regions from news frames. Two classifers have been trained and evaluated for the labeling of the detected bands as natural or artificial - Support Vector Machine (SVM) Classifer with RBF kernel, and Extreme Learning Machine (ELM) classifier. The classifiers have been trained on a dataset of 6000 images (3000 images of each class). The ELM classifier reports a balanced accuracy of 77.38%, while the SVM classifier outperforms it with a balanced accuracy of 96.5% using 10-fold cross-validation. The detected bands which have been fragmented due to the presence of gradients in the image have been merged using a three-tier hierarchical reasoning model. The bands were detected with a Jaccard Index of 0.8138, when compared to manually marked ground truth data. We have also presented an extensive literature review of previous work done towards news videos format detection, element band classification, and associative reasoning.


title: Fuzzy Profit Shifting: A Model for Optimal Tax-induced Transfer Pricing with Fuzzy Arm's Length Parameter

abstract: This paper proposes a model of optimal tax-induced transfer pricing with a fuzzy arm's length parameter. Fuzzy numbers provide a suitable structure for modelling the ambiguity that is intrinsic to the arm's length parameter. For the usual conditions regarding the anti-shifting mechanisms, the optimal transfer price becomes a maximising $\alpha$-cut of the fuzzy arm's length parameter. Nonetheless, we show that it is profitable for firms to choose any maximising transfer price if the probability of tax audit is sufficiently low, even if the chosen price is considered a completely non-arm's length price by tax authorities. In this case, we derive the necessary and sufficient conditions to prevent this extreme shifting strategy


title: Multiplexing More Streams in the MU-MISO Downlink by Interference Exploitation Precoding

abstract: In this paper, we study the interference exploitation precoding for the scenario where the number of streams simultaneously transmitted by the base station (BS) is larger than that of transmit antennas at the BS, and derive the optimal precoding structure by employing the pseudo inverse. We show that the optimal pre-scaling vector is equal to a linear combination of the right singular vectors that correspond to zero singular values of the coefficient matrix. By formulating the dual problem, the optimal precoding matrix can be expressed as a function of the dual variables in a closed form, and an equivalent quadratic programming (QP) formulation is further derived for computational complexity reduction. Numerical results validate our analysis and demonstrate significant performance improvements for interference exploitation precoding for the considered scenario.


title: Production Spectra of $^3$He($\pi$, $K$) Reactions with Continuum Discretized Coupled Channels

abstract: We investigate theoretically $\Lambda$ production spectra of $^3$He($\pi$, $K$) reactions at $p_\pi=$ 1.05--1.20 GeV/$c$ in the distorted-wave impulse approximation, using the continuum-discretized coupled-channel method. The production cross section of a $^3_\Lambda$H(1/2$^+$) ground state is also discussed.


title: Projection-based reduced order models for a cut finite element method in parametrized domains

abstract: This work presents a reduced order modelling technique built on a high fidelity embedded mesh finite element method. Such methods, and in particular the CutFEM method, are attractive in the generation of projection-based reduced order models thanks to their capabilities to seamlessly handle large deformations of parametrized domains. The combination of embedded methods and reduced order models allows us to obtain fast evaluation of parametrized problems, avoiding remeshing as well as the reference domain formulation, often used in the reduced order modelling for boundary fitted finite element formulations. The resulting novel methodology is presented on linear elliptic and Stokes problems, together with several test cases to assess its capability. The role of a proper extension and transport of embedded solutions to a common background is analyzed in detail.


title: Rayleigh-B\'enard convection with a melting boundary

abstract: We study the evolution of a melting front between the solid and liquid phases of a pure incompressible material where fluid motions are driven by unstable temperature gradients. In a plane layer geometry, this can be seen as classical Rayleigh-B\'enard convection where the upper solid boundary is allowed to melt due to the heat flux brought by the fluid underneath. This free-boundary problem is studied numerically in two dimensions using a phase-field approach, classically used to study the melting and solidification of alloys, which we dynamically couple with the Navier-Stokes equations in the Boussinesq approximation. The advantage of this approach is that it requires only moderate modifications of classical numerical methods. We focus on the case where the solid is initially nearly isothermal, so that the evolution of the topography is related to the inhomogeneous heat flux from thermal convection, and does not depend on the conduction problem in the solid. From a very thin stable layer of fluid, convection cells appears as the depth -- and therefore the effective Rayleigh number of the layer increases. The continuous melting of the solid leads to dynamical transitions between different convection cell sizes and topography amplitudes. The Nusselt number can be larger than its value for a planar upper boundary, due to the feedback of the topography on the flow, which can stabilize large-scale laminar convection cells.


title: Supersymmetrization of the Franke-Gorini-Kossakowski-Lindblad-Sudarshan equation

abstract: We investigate departures from the Hamiltonian dynamics such as the quantum decoherence in a supersymmetric quantum system due to its interaction with the environment. The effective description is given by the Franke-Gorini-Kossakowski-Lindblad-Sudarshan (FGKLS) equation. The prescription for supersymmetrization of this equation is elaborated and pairs of related isospectral quantum systems are built. We present some instructive examples of such a construction.


title: On the Non-uniqueness of Representations of Coxian Phase-Type Distributions

abstract: Parameter estimation in Coxian phase-type models can be challenging due to their non-unique representation leading to a multi-modal likelihood. Since each representation corresponds to a different underlying data-generating mechanism, it is of interest to identify those supported by given data (i.e., find all likelihood modes). The standard approach is to simply refit using various initial values, but this has no guarantee of working. Thus, we develop new properties specific to this class of models, and employ these to determine all the equivalent model representations. The proposed approach only requires fitting the model once, and is guaranteed to find all representations.


title: Stochastic Nicholson's blowflies delay differential equation with regime switching

abstract: In this paper, we investigate the global existence of almost surely positive solution to a stochastic Nicholson's blowflies delay differential equation with regime switching, and give the estimation of the path. The results presented in this paper extend some corresponding results in Wang et al. Stochastic Nicholson's Blowflies Delayed Differential Equations, Appl. Math. Lett. 87 (2019) 20-26.


title: Physical properties of the X-ray gas as a dynamical diagnosis for galaxy clusters

abstract: We analysed XMM-{\it Newton} EPIC data for 53 galaxy clusters. Through 2D spectral maps, we provide the most detailed and extended view of the spatial distribution of temperature (kT), pressure (P), entropy (S) and metallicity (Z) of galaxy clusters to date with the aim of correlating the dynamical state of the system to six cool-core diagnoses from the literature. With the objective of building 2D maps and resolving structures in kT, P, S and Z, we divide the data in small regions from which spectra can be extracted. Our analysis shows that when clusters are spherically symmetric the cool-cores (CC) are preserved, the systems are relaxed with little signs of perturbation, and most of the CC criteria agree. The disturbed clusters are elongated, show clear signs of interaction in the 2D maps, and most do not have a cool-core. However, 16 well studied clusters classified as CC by at least four criteria show spectral maps that appear disturbed. All of these clusters but one show clear signs of recent mergers, with a complex structure and geometry but with a cool-core that remains preserved. Thus, although very useful for CC characterization, most diagnoses are too simplistic to reproduce the overall structure and dynamics of galaxy clusters, and therefore the selection of relaxed systems according to these criteria may affect mass estimates. The complex structure of galaxy clusters can be reliably assessed through the 2D maps presented here.


title: One-view occlusion detection for stereo matching with a fully connected CRF model

abstract: In this paper, we extend the standard belief propagation (BP) sequential technique proposed in the tree-reweighted sequential method to the fully connected CRF models with the geodesic distance affinity. The proposed method has been applied to the stereo matching problem. Also a new approach to the BP marginal solution is proposed that we call one-view occlusion detection (OVOD). In contrast to the standard winner takes all (WTA) estimation, the proposed OVOD solution allows to find occluded regions in the disparity map and simultaneously improve the matching result. As a result we can perform only one energy minimization process and avoid the cost calculation for the second view and the left-right check procedure. We show that the OVOD approach considerably improves results for cost augmentation and energy minimization techniques in comparison with the standard one-view affinity space implementation. We apply our method to the Middlebury data set and reach state-of-the-art especially for median, average and mean squared error metrics.


title: Caristi-Kirk and Oettli-Th\'era Ball Spaces and applications

abstract: Based on the theory of ball spaces introduced by Kuhlmann and Kuhlmann we introduce and study Caristi-Kirk and Oettli-Th\'era ball spaces. We show that if the underlying metric space is complete, then these have a very strong property: every ball contains a singleton ball. This fact provides quick proofs for several results which are equivalent to the Caristi-Kirk Fixed Point Theorem, namely Ekeland's Variational Principles, the Oettli-Th\'era Theorem, Takahashi's Theorem and the Flower Petal Theorem.


title: Almost sure global well posedness for the BBM equation with infinite $L^{2}$ initial data

abstract: We consider the probabilistic Cauchy problem for the Benjamin-Bona-Mahony equation (BBM) on the one-dimensional torus $\mathbb{T}$ with initial data below $L^{2}(\mathbb{T})$. With respect to random initial data of strictly negative Sobolev regularity, we prove that BBM is almost surely globally well-posed. The argument employs the $I$-method to obtain an a priori bound on the growth of the `residual' part of the solution. We then discuss the stability properties of the solution map in the deterministically ill-posed regime.


title: Revisit to electrical and thermal conductivities, Lorenz number and Knudsen number in thermal QCD in a strong magnetic field

abstract: We have explored how the electrical ($\sigma_{\rm el}$) and thermal ($\kappa$) conductivities in a thermal QCD medium get affected in weak-momentum anisotropies arising either due to a strong magnetic field or due to asymptotic expansion. This study facilitates to understand the longevity of strong magnetic field through $\sigma_{el}$, Lorenz number in Wiedemann-Franz law, and the validity of equilibrium by the Knudsen number. We calculate the conductivities by solving relativistic Boltzmann transport equation in relaxation-time approximation within quasiparticle model at finite T and strong B. We have found that $\sigma_{el}$ and $\kappa$ get enhanced in a magnetic field-driven anisotropy, but $\sigma_{el}$ decreases with temperature, opposite to its faster increase in expansion-driven anisotropy. Whereas $\kappa$ increases slowly with temperature, contrary to its rapid increase in expansion-driven anisotropy. The above findings are broadly attributed to three factors: the stretching and squeezing of distribution function in anisotropies generated by the magnetic field and asymptotic expansion, respectively, the dispersion relation and resulting phase-space factor, the relaxation-time in absence and presence of strong magnetic field. So $\sigma_{\rm el}$ extracts the time-dependence of magnetic field, which decays slower than in vacuum but expansion-driven anisotropy makes the decay faster. The variation in $\kappa$ transpires that Knudsen number decreases with T but expansion-driven anisotropy reduces its value and magnetic field-driven anisotropy raises its value but to less than one, thus the system can still be in equilibrium. The ratio, $\kappa/\sigma_{el}$ in magnetic field-driven anisotropy increases linearly with temperature but with a value smaller than in expansion-driven anisotropy. Thus the Lorenz number can make the distinction between different anisotropies.


title: First results of EPOS-HQ model for open heavy flavor production in A-A collission at RHIC and LHC

abstract: The EPOS-HQ model is a framework designed for predicting heavy-flavor observables both in p-p, p-A and A-A collisions. It results from the integration of the MC$@_s$HQ and EPOS3 codes and should supersede predictions made with each of these two. In this contribution, we briefly present the new ingredients of the EPOS-HQ model as well as some preliminary results. From the comparison with experimental data, we derive some direction of improvement for the near future.


title: Deep-learning-based identification of odontogenic keratocysts in hematoxylin- and eosin-stained jaw cyst specimens

abstract: The aim of this study was to develop a digital histopathology system for identifying odontogenic keratocysts in hematoxylin- and eosin-stained tissue specimens of jaw cysts. Approximately 5000 microscopy images with 400$\times$ magnification were obtained from 199 odontogenic keratocysts, 208 dentigerous cysts, and 55 radicular cysts. A proportion of these images were used to make training patches, which were annotated as belonging to one of the following three classes: keratocysts, non-keratocysts, and stroma. The patches for the cysts contained the complete lining epithelium, with the cyst cavity being present on the upper side. The convolutional neural network (CNN) VGG16 was finetuned to this dataset. The trained CNN could recognize the basal cell palisading pattern, which is the definitive criterion for diagnosing keratocysts. Some of the remaining images were scanned and analyzed by the trained CNN, whose output was then used to train another CNN for binary classification (keratocyst or not). The area under the receiver operating characteristics curve for the entire algorithm was 0.997 for the test dataset. Thus, the proposed patch classification strategy is usable for automated keratocyst diagnosis. However, further optimization must be performed to make it suitable for practical use.


title: Operator means of probability measures

abstract: Let $\mathbb{P}$ be the complete metric space consisting of positive invertible operators on an infinite-dimensional Hilbert space with the Thompson metric. We introduce the notion of operator means of probability measures on $\mathbb{P}$, in parallel with Kubo and Ando's definition of two-variable operator means, and show that every operator mean is contractive for the $\infty$-Wasserstein distance. By means of a fixed point method we consider deformation of such operator means, and show that the deformation of any operator mean becomes again an operator mean in our sense. Based on this deformation procedure we prove a number of properties and inequalities for operator means of probability measures.


title: What comes next? Extractive summarization by next-sentence prediction

abstract: Existing approaches to automatic summarization assume that a length limit for the summary is given, and view content selection as an optimization problem to maximize informativeness and minimize redundancy within this budget. This framework ignores the fact that human-written summaries have rich internal structure which can be exploited to train a summarization system. We present NEXTSUM, a novel approach to summarization based on a model that predicts the next sentence to include in the summary using not only the source article, but also the summary produced so far. We show that such a model successfully captures summary-specific discourse moves, and leads to better content selection performance, in addition to automatically predicting how long the target summary should be. We perform experiments on the New York Times Annotated Corpus of summaries, where NEXTSUM outperforms lead and content-model summarization baselines by significant margins. We also show that the lengths of summaries produced by our system correlates with the lengths of the human-written gold standards.


title: Prototypical Metric Transfer Learning for Continuous Speech Keyword Spotting With Limited Training Data

abstract: Continuous Speech Keyword Spotting (CSKS) is the problem of spotting keywords in recorded conversations, when a small number of instances of keywords are available in training data. Unlike the more common Keyword Spotting, where an algorithm needs to detect lone keywords or short phrases like "Alexa", "Cortana", "Hi Alexa!", "Whatsup Octavia?" etc. in speech, CSKS needs to filter out embedded words from a continuous flow of speech, ie. spot "Anna" and "github" in "I know a developer named Anna who can look into this github issue." Apart from the issue of limited training data availability, CSKS is an extremely imbalanced classification problem. We address the limitations of simple keyword spotting baselines for both aforementioned challenges by using a novel combination of loss functions (Prototypical networks' loss and metric loss) and transfer learning. Our method improves F1 score by over 10%.


title: HorizonNet: Learning Room Layout with 1D Representation and Pano Stretch Data Augmentation

abstract: We present a new approach to the problem of estimating the 3D room layout from a single panoramic image. We represent room layout as three 1D vectors that encode, at each image column, the boundary positions of floor-wall and ceiling-wall, and the existence of wall-wall boundary. The proposed network, HorizonNet, trained for predicting 1D layout, outperforms previous state-of-the-art approaches. The designed post-processing procedure for recovering 3D room layouts from 1D predictions can automatically infer the room shape with low computation cost - it takes less than 20ms for a panorama image while prior works might need dozens of seconds. We also propose Pano Stretch Data Augmentation, which can diversify panorama data and be applied to other panorama-related learning tasks. Due to the limited data available for non-cuboid layout, we relabel 65 general layout from the current dataset for finetuning. Our approach shows good performance on general layouts by qualitative results and cross-validation.


title: One-loop renormalization of staple-shaped operators in continuum and lattice regularizations

abstract: In this paper we present one-loop results for the renormalization of nonlocal quark bilinear operators, containing a staple-shaped Wilson line, in both continuum and lattice regularizations. The continuum calculations were performed in dimensional regularization, and the lattice calculations for the Wilson/clover fermion action and for a variety of Symanzik-improved gauge actions. We extract the strength of the one-loop linear and logarithmic divergences (including cusp divergences), which appear in such nonlocal operators; we identify the mixing pairs which occur among some of these operators on the lattice, and we calculate the corresponding mixing coefficients. We also provide the appropriate RI'-like scheme, which disentangles this mixing nonperturbatively from lattice simulation data, as well as the one-loop expressions of the conversion factors, which turn the lattice data to the MS-bar scheme. Our results can be immediately used for improving recent nonperturbative investigations of transverse momentum-dependent distribution functions (TMDs) on the lattice. Finally, extending our perturbative study to general Wilson-line lattice operators with n cusps, we present results for their renormalization factors, including identification of mixing and determination of the corresponding mixing coefficients, based on our results for the staple operators.


title: Sharp-edged geometric obstacles in microfluidics promote deformability-based sorting of cells

abstract: Sorting cells based on their intrinsic properties is a highly desirable objective, since changes in cell deformability are often associated with various stress conditions and diseases. Deterministic lateral displacement (DLD) devices offer high precision for rigid spherical particles, while their success in sorting deformable particles remains limited due to the complexity of cell traversal in DLDs. We employ mesoscopic hydrodynamics simulations and demonstrate prominent advantages of sharp-edged DLD obstacles for probing deformability properties of red blood cells (RBCs). By consecutive sharpening of the pillar shape from circular to diamond to triangular geometry, a pronounced cell bending around an edge is achieved, serving as a deformability sensor. Bending around the edge is the primary mechanism, which governs the traversal of RBCs through such DLD device. This strategy requires an appropriate degree of cell bending by fluid stresses, which can be controlled by the flow rate, and exhibits good sensitivity to moderate changes in cell deformability. We expect that similar mechanisms should be applicable for the development of novel DLD devices that target intrinsic properties of many other cells.


title: Anatomy of the vertebral column lymphatic network in mice

abstract: Cranial lymphatic vessels (LVs) are involved in transport of fluids, macromolecules and CNS immune responses. Little information about spinal LVs is available, because these delicate structures are embedded within vertebral tissues and difficult to visualize using traditional histology. Here we reveal an extended vertebral column LV network using three-dimensional imaging of decalcified iDISCO-clarified spine segments. Spinal LVs are metameric circuits exiting along spinal nerve roots and connecting to lymph nodes and the thoracic duct. They navigate in the epidural space and the dura mater around the spinal cord, and associate with leukocytes, peripheral dorsal root and sympathetic ganglia. Spinal LVs are VEGF-C-dependent and remodel extensively after spinal cord injury. They constitute an extension to cranial circuits for meningeal fluids, but also a route for perineural fluids and a link with peripheral immune and nervous circuits. Vertebral column LVs may be potential targets to improve the maintenance and repair of 32 spinal tissues as well as gatekeepers of CNS immunity.


title: An Experience Report On Applying Software Testing Academic Results In Industry: We Need Usable Automated Test Generation

abstract: What is the impact of software engineering research on current practices in industry? In this paper, I report on my direct experience as a PhD/post-doc working in software engineering research projects, and then spending the following five years as an engineer in two different companies (the first one being the same I worked in collaboration with during my post-doc). Given a background in software engineering research, what cutting-edge techniques and tools from academia did I use in my daily work when developing and testing the systems of these companies? Regarding validation and verification (my main area of research), the answer is rather short: as far as I can tell, only FindBugs. In this paper, I report on why this was the case, and discuss all the challenging, complex open problems we face in industry and which somehow are "neglected" in the academic circles. In particular, I will first discuss what actual tools I could use in my daily work, such as JaCoCo and Selenium. Then, I will discuss the main open problems I faced, particularly related to environment simulators, unit and web testing. After that, popular topics in academia are presented, such as UML, regression and mutation testing. Their lack of impact on the type of projects I worked on in industry is then discussed. Finally, from this industrial experience, I provide my opinions about how this situation can be improved, in particular related to how academics are evaluated, and advocate for a greater involvement into open-source projects.


title: HAS-QA: Hierarchical Answer Spans Model for Open-domain Question Answering

abstract: This paper is concerned with open-domain question answering (i.e., OpenQA). Recently, some works have viewed this problem as a reading comprehension (RC) task, and directly applied successful RC models to it. However, the performances of such models are not so good as that in the RC task. In our opinion, the perspective of RC ignores three characteristics in OpenQA task: 1) many paragraphs without the answer span are included in the data collection; 2) multiple answer spans may exist within one given paragraph; 3) the end position of an answer span is dependent with the start position. In this paper, we first propose a new probabilistic formulation of OpenQA, based on a three-level hierarchical structure, i.e.,~the question level, the paragraph level and the answer span level. Then a Hierarchical Answer Spans Model (HAS-QA) is designed to capture each probability. HAS-QA has the ability to tackle the above three problems, and experiments on public OpenQA datasets show that it significantly outperforms traditional RC baselines and recent OpenQA baselines.


title: CP violating triple product asymmetries in charm decays

abstract: $CP$ violation asymmetry is expected to be very small within the Standard Model. Observation of $CP$ violation in charm sector would hint towards New Physics. Time-reversal asymmetry is sensitive to $CPV$ via the $CPT$ Theorem, and is a clean and alternative way to search for CP violation in charm decays. This is achieved by the use of triple-product correlations in 4-body charm meson decays. The results from previous searches at BaBar and LHCb Collaborations and the latest results from Belle are reported. The most recent result from Belle performed in $D^{0} \to K^{+}K^{-}\pi^{+}\pi^{-}$ decays using kinematic asymmetry variables is the first measurement of its kind.


title: Multiplayer Multi-armed Bandits for Optimal Assignment in Heterogeneous Networks

abstract: We consider an ad hoc network where multiple users access the same set of channels. The channel characteristics are unknown and could be different for each user (heterogeneous). No controller is available to coordinate channel selections by the users, and if multiple users select the same channel, they collide and none of them receive any rate (or reward). For such a completely decentralized network we develop algorithms that aim to achieve optimal network throughput. Due to lack of any direct communication between the users, we allow each user to exchange information by transmitting in a specific pattern and sense such transmissions from others. However, such transmissions and sensing for information exchange do not add to network throughput. For the wideband sensing and narrowband sensing scenarios, we first develop explore-and-commit algorithms that converge to near-optimal allocation with high probability in a small number of rounds. Building on this, we develop an algorithm that gives logarithmic regret, even when the number of users changes with time. We validate our claims through extensive experiments and show that our algorithms perform significantly better than the state-of-the-art CSM-MAB, dE3 and dE3-TS algorithms.


title: Magnetospheric Multiscale Observation of Kinetic Signatures in the Alfv\'en Vortex

abstract: Alfv\'en vortex is a multi-scale nonlinear structure which contributes to intermittency of turbulence. Despite previous explorations mostly on the spatial properties of the Alfv\'en vortex (i.e., scale, orientation, and motion), the plasma characteristics within the Alfv\'en vortex are unknown. Moreover, the connection between the plasma energization and the Alfv\'en vortex still remains unclear. Based on high resolution in-situ measurement from the Magnetospheric Multiscale (MMS) mission, we report for the first time, distinctive plasma features within an Alfv\'en vortex. This Alfv\'en vortex is identified to be two-dimensional ($k_{\bot} \gg k_{\|}$) quasi-monopole with a radius of ~10 proton gyroscales. Its magnetic fluctuations $\delta B_{\bot}$ are anti correlated with velocity fluctuations $\delta V_{\bot}$, thus the parallel current density $j_{\|}$ and flow vorticity $\omega_{\|}$ are anti-aligned. In different part of the vortex (i.e., edge, middle, center), the ion and electron temperatures are found to be quite different and they behave in the reverse trend: the ion temperature variations are correlated with $j_{\|}$, while the electron temperature variations are correlated with $\omega_{\|}$. Furthermore, the temperature anisotropies, together with the non-Maxwellian kinetic effects, exhibit strong enhancement at peaks of $|\omega_{\|}| (|j_{\|}|)$ within the vortex. Comparison between observations and numerical/theoretical results are made. In addition, the energy-conversion channels and the compressibility associated with the Alfv\'en vortex are discussed. These results may help to understand the link between coherent vortex structures and the kinetic processes, which determines how turbulence energy dissipate in the weakly-collisional space plasmas.


title: Preservation of the invariants of Lotka-Volterra equations by iterated deferred correction methods

abstract: In this paper we apply Kahan's nonstandard discretization to three dimensional Lotka-Volterra equations in bi-Hamiltonian form. The periodicity of the solutions and all polynomial and non-polynomial invariants are well preserved in long-term integration. Applying classical deferred correction method, we show that the invariants are preserved with increasing accuracy as a results of more accurate numerical solutions. Substantial speedups over the Kahan's method are achieved at each run with deferred correction method.


title: Precision Analysis of the $^{136}$Xe Two-Neutrino $\beta\beta$ Spectrum in KamLAND-Zen and Its Impact on the Quenching of Nuclear Matrix Elements

abstract: We present a precision analysis of the $^{136}$Xe two-neutrino $\beta\beta$ electron spectrum above 0.8 MeV, based on high-statistics data obtained with the KamLAND-Zen experiment. An improved formalism for the two-neutrino $\beta\beta$ rate allows us to measure the ratio of the leading and subleading $2\nu\beta\beta$ nuclear matrix elements (NMEs), $\xi^{2\nu}_{31} = -0.26^{+0.31}_{-0.25}$. Theoretical predictions from the nuclear shell model and the majority of the quasiparticle random-phase approximation (QRPA) calculations are consistent with the experimental limit. However, part of the $\xi^{2\nu}_{31}$ range allowed by the QRPA is excluded by the present measurement at the 90% confidence level. Our analysis reveals that predicted $\xi^{2\nu}_{31}$ values are sensitive to the quenching of NMEs and the competing contributions from low- and high-energy states in the intermediate nucleus. Since these aspects are also at play in neutrinoless $\beta\beta$ decay, $\xi^{2\nu}_{31}$ provides new insights towards reliable neutrinoless $\beta\beta$ NMEs.


title: Nonparametric Inverse Dynamic Models for Multimodal Interactive Robots

abstract: Direct design of a robot's rendered dynamics, such as in impedance control, is now a well-established control mode in uncertain environments. When the physical interaction port variables are not measured directly, dynamic and kinematic models are required to relate the measured variables to the interaction port variables. A typical example is serial manipulators with joint torque sensors, where the interaction occurs at the end-effector. As interactive robots perform increasingly complex tasks, they will be intermittently coupled with additional dynamic elements such as tools, grippers, or workpieces, some of which should be compensated and brought to the robot side of the interaction port, making the inverse dynamics multimodal. Furthermore, there may also be unavoidable and unmeasured external input when the desired system cannot be totally isolated. Towards semi-autonomous robots, capable of handling such applications, a multimodal Gaussian process regression approach to manipulator dynamic modelling is developed. A sampling-based approach clusters different dynamic modes from unlabelled data, also allowing the seperation of perturbed data with significant, irregular external input. The passivity of the overall approach is shown analytically, and experiments examine the performance and safety of this approach on a test actuator.


title: Semileptonic and nonleptonic decays of $D$ into tensor mesons with light-cone sum rule

abstract: Form factors of $D$ decays into $J^{PC}=2^{++}$ tensor mesons are calculated in the light-cone sum rules approach up to twist-4 distribution amplitudes of the tensor meson. The masses of the tensor mesons are comparable to that of the charm quark mass $m_c$; therefore all terms including powers of $m_T/m_c$ are kept out in the expansion of the two-particle distribution amplitude $\langle T|\bar q_{1\alpha}(x) \, q_{2\delta}(0)|0\rangle$. Branching ratios of the semileptonic $D \to T\,\mu \,\bar{\nu}_{\mu}$ decays and nonleptonic $D\to T~ P ~(P=K, \pi)$ decays are taken into consideration. A comparison is also made between our results and predictions of other methods and the existing experimental values for the nonleptonic case. The semileptonic branching ratios are typically of the order of $10^{-5}$, and the nonleptonic ones show better agreement with the experimental data in comparison to the Isgur-Scora-Grinstein-Wise predictions.


title: A Risk-Sharing Framework of Bilateral Contracts

abstract: We introduce a two-agent problem which is inspired by price asymmetry arising from funding difference. When two parties have different funding rates, the two parties deduce different fair prices for derivative contracts even under the same pricing methodology and parameters. Thus, the two parties should enter the derivative contracts with a negotiated price, and we call the negotiation a risk-sharing problem. This framework defines the negotiation as a problem that maximizes the sum of utilities of the two parties. By the derived optimal price, we provide a theoretical analysis on how the price is determined between the two parties. As well as the price, the risk-sharing framework produces an optimal amount of collateral. The derived optimal collateral can be used for contracts between financial firms and non-financial firms. However, inter-dealers markets are governed by regulations. As recommended in Basel III, it is a convention in inter-dealer contracts to pledge the full amount of a close-out price as collateral. In this case, using the optimal collateral, we interpret conditions for the full margin requirement to be indeed optimal.