title: Estimation of lactate threshold with machine learning techniques in recreational runners

abstract: Lactate threshold is considered an essential parameter when assessing performance of elite and recreational runners and prescribing training intensities in endurance sports. However, the measurement of blood lactate concentration requires expensive equipment and the extraction of blood samples, which are inconvenient for frequent monitoring. Furthermore, most recreational runners do not have access to routine assessment of their physical fitness by the aforementioned equipment so they are not able to calculate the lactate threshold without resorting to an expensive and specialized centre. Therefore, the main objective of this study is to create an intelligent system capable of estimating the lactate threshold of recreational athletes participating in endurance running sports. The solution here proposed is based on a machine learning system which models the lactate evolution using recurrent neural networks and includes the proposal of standardization of the temporal axis as well as a modification of the stratified sampling method. The results show that the proposed system accurately estimates the lactate threshold of 89.52% of the athletes and its correlation with the experimentally measured lactate threshold is very high (R=0,89). Moreover, its behaviour with the test dataset is as good as with the training set, meaning that the generalization power of the model is high. Therefore, in this study a machine learning based system is proposed as alternative to the traditional invasive lactate threshold measurement tests for recreational runners.


title: Optimal Bipartite Network Clustering

abstract: We study bipartite community detection in networks, or more generally the network biclustering problem. We present a fast two-stage procedure based on spectral initialization followed by the application of a pseudo-likelihood classifier twice. Under mild regularity conditions, we establish the weak consistency of the procedure (i.e., the convergence of the misclassification rate to zero) under a general bipartite stochastic block model. We show that the procedure is optimal in the sense that it achieves the optimal convergence rate that is achievable by a biclustering oracle, adaptively over the whole class, up to constants. This is further formalized by deriving a minimax lower bound over a class of biclustering problems. The optimal rate we obtain sharpens some of the existing results and generalizes others to a wide regime of average degree growth, from sparse networks with average degrees growing arbitrarily slowly to fairly dense networks with average degrees of order $\sqrt{n}$. As a special case, we recover the known exact recovery threshold in the $\log n$ regime of sparsity. To obtain the consistency result, as part of the provable version of the algorithm, we introduce a sub-block partitioning scheme that is also computationally attractive, allowing for distributed implementation of the algorithm without sacrificing optimality. The provable algorithm is derived from a general class of pseudo-likelihood biclustering algorithms that employ simple EM type updates. We show the effectiveness of this general class by numerical simulations.


title: You Watch, You Give, and You Engage: A Study of Live Streaming Practices in China

abstract: Despite gaining traction in North America, live streaming has not reached the popularity it has in China, where livestreaming has a tremendous impact on the social behaviors of users. To better understand this socio-technological phenomenon, we conducted a mixed methods study of live streaming practices in China. We present the results of an online survey of 527 live streaming users, focusing on their broadcasting or viewing practices and the experiences they find most engaging. We also interviewed 14 active users to explore their motivations and experiences. Our data revealed the different categories of content that was broadcasted and how varying aspects of this content engaged viewers. We also gained insight into the role reward systems and fan group-chat play in engaging users, while also finding evidence that both viewers and streamers desire deeper channels and mechanisms for interaction in addition to the commenting, gifting, and fan groups that are available today.


title: $W^{1,p}$ regularity of solutions to Kolmogorov equation and associated Feller semigroup

abstract: In $\mathbb R^d$, $d \geq 3$, consider the divergence and the non-divergence form operators \begin{equation} \tag{$i$} - \nabla \cdot a \cdot \nabla + b \cdot \nabla, \end{equation} \begin{equation} \tag{$ii$} - a \cdot \nabla^2 + b \cdot \nabla, \end{equation} where $a=I+c \mathsf{f} \otimes \mathsf{f}$, the vector fields $\nabla_i \mathsf{f}$ ($i=1,2,\dots,d$) and $b$ are form-bounded (this includes the sub-critical class $[L^d + L^\infty]^d$ as well as vector fields having critical-order singularities). We characterize quantitative dependence on $c$ and the values of the form-bounds of the $L^q \rightarrow W^{1,qd/(d-2)}$ regularity of the resolvents of the operator realizations of ($i$), ($ii$) in $L^q$, $q \geq 2 \vee ( d-2)$ as (minus) generators of positivity preserving $L^\infty$ contraction $C_0$ semigroups. The latter allows to run an iteration procedure $L^p \rightarrow L^\infty$ that yields associated with ($i$), ($ii$) $L^q$-strong Feller semigroups.


title: Multistage stochastic programs with a random number of stages: dynamic programming equations, solution methods, and application to portfolio selection

abstract: We introduce the class of multistage stochastic optimization problems with a random number of stages. For such problems, we show how to write dynamic programming equations and detail the Stochastic Dual Dynamic Programming algorithm to solve these equations. Finally, we consider a portfolio selection problem over an optimization period of random duration. For several instances of this problem, we show the gain obtained using a policy that takes the random duration of the number of stages into account over a policy built taking a fixed number of stages (namely the maximal possible number of stages).


title: Reputation is required for cooperation to emerge in dynamic networks

abstract: Melamed, Harrell, and Simpson have recently reported on an experiment which appears to show that cooperation can arise in a dynamic network without reputational knowledge, i.e., purely via dynamics [1]. We believe that their experimental design is actually not testing this, in so far as players do know the last action of their current partners before making a choice on their own next action and subsequently deciding which link to cut. Had the authors given no information at all, the result would be a decline in cooperation as shown in [2].


title: A Role for Turbulence in Circumgalactic Precipitation

abstract: Abundant observational evidence indicates that the cooling time t_cool of the hot ambient medium pervading a massive galaxy does not drop much below 10 times the freefall time t_ff at any radius. Theoretical models have accounted for this finding by hypothesizing that cold clouds start to condense out of the ambient medium when t_cool/t_ff < 10 and fuel a strong black-hole feedback response that reheats the ambient gas, but those models have not yet been able to provide a simple explanation for the origin of the critical t_cool/t_ff ratio. This paper explores a heuristic model for condensation that links the critical ratio to turbulent driving of gravity-wave oscillations. In the linear regime, internal gravity waves are thermally unstable in a thermally balanced medium. Buoyancy oscillations in a balanced medium with t_cool/t_ff therefore grow until they saturate without condensing at an amplitude that depends on t_cool/t_ff. However, in a medium with 10 < t_cool/t_ff < 20, turbulence with a velocity dispersion roughly half the galaxy's stellar velocity dispersion can drive those oscillations into condensation. Intriguingly, this is indeed the gas-phase velocity dispersion observed among galaxy-cluster cores that contain multiphase gas. It is therefore possible that both the critical t_cool/t_ff ratio for condensation of ambient gas and the level of turbulence in that gas are determined by coupling between condensation, feedback, and turbulence. Such a system can converge to a well-regulated equilibrium state, as long as the fraction of feedback energy that goes into turbulence is significantly less than the fraction that goes more directly into heat.


title: Anderson Localization for Radial Tree Graphs With Random Branching Numbers

abstract: We prove Anderson localization for the discrete Laplace operator on radial tree graphs with random branching numbers. Our method relies on the representation of the Laplace operator as the direct sum of half-line Jacobi matrices whose entries are non-degenerate, independent, identically distributed random variables with singular distributions.


title: Optimality of multi-refraction dividend strategies in the dual model

abstract: We consider the multi-refraction strategies in two equivalent versions of the optimal dividend problem in the dual (spectrally positive L\'evy) model. The first problem is a variant of the bail-out case where both dividend payments and capital injections must be absolutely continuous with respect to the Lebesgue measure. The second is an extension of Avanzi et al. [4] where a strategy is a combination of two absolutely continuous dividend payments with different upper bounds and different transaction costs. In both problems, it is shown to be optimal to refract the process at two thresholds, with the optimally controlled process being the multi-refracted L\'evy process recently studied by Czarna et al. [9]. The optimal strategy and the value function are succinctly written in terms of a version of the scale function. Numerical results are also given.


title: Generic asymptotics of resonance counting function for Schr\"odinger point interactions

abstract: We study the leading coefficient in the asymptotical formula $ N (R) = \frac{W}{\pi} R + O (1) $, $ R \to \infty $, for the resonance counting function $ N (R)$ of Schr\"odinger Hamiltonians with point interactions. For such Hamiltonians, the Weyl-type and non-Weyl-type asymptotics of $N (R)$ was introduced recently in a paper by J. Lipovsk\'y and V. Lotoreichik (2017). In this note, we prove that the Weyl-type asymptotics is generic.


title: False discovery rate control for multiple testing based on p-values with c\`adl\`ag distribution functions

abstract: For multiple testing based on p-values with c\`{a}dl\`{a}g distribution functions, we propose an FDR procedure "BH+" with proven conservativeness. BH+ is at least as powerful as the BH procedure when they are applied to super-uniform p-values. Further, when applied to mid p-values, BH+ is more powerful than it is applied to conventional p-values. An easily verifiable necessary and sufficient condition for this is provided. BH+ is perhaps the first conservative FDR procedure applicable to mid p-values. BH+ is applied to multiple testing based on discrete p-values in a methylation study, an HIV study and a clinical safety study, where it makes considerably more discoveries than the BH procedure.


title: Matroids and Codes with the Rank Metric

abstract: We study the relationship between a q-analogue of matroids and linear codes with the rank metric in the vector space of matrices with entries in a finite field. We prove a Greene type identity for the rank generating function of these matroidal structures and the rank weight enumerator of these linear codes. As an application, we give a combinatorial proof of a MacWilliams type identity for Delsarte rank-metric codes.


title: Stability Engineering of Halide Perovskite via Machine Learning

abstract: Perovskite stability is of the core importance and difficulty in current research and application of perovskite solar cells. Nevertheless, over the past century, the formability and stability of perovskite still relied on simplified factor based on human knowledge, such as the commonly used tolerance factor t. Combining machine learning (ML) with first-principles density functional calculations, we proposed a strategy to firstly calculate the decomposition energies, considered to be closely related to thermodynamic stability, of 354 kinds halide perovskites, establish the machine learning relationship between decomposition energy and compositional ionic radius and investigate the stabilities of 14,190 halide double perovskites. The ML-predicted results enable us to rediscover a series of stable rare earth metal halide perovskites (up to ~1000 kinds), indicating the generalization of this model and further provide elemental and concentration suggestion for improving the stability of mixed perovskite.


title: Chiral kinetic approach to the chiral magnetic effect in isobaric collisions

abstract: Based on the chiral kinetic approach using quarks and antiquarks from a multiphase transport model as initial conditions, we study the chiral magnetic effect, i.e., the magnetic field induced separation of charged particles in the transverse plane, in non-central isobaric collisions of Zr$+$Zr and Ru$+$Ru, which have the same atomic number but different proton numbers. For the observable $\gamma^{OS}-\gamma^{SS}$ related to the difference between the correlations of particles of opposite charges and of same charges, we find a difference between the two collision systems if the magnetic field has a long lifetime of 0.6 fm$/c$ and the observable is evaluated using the initial reaction plane. This signal of the chiral magnetic effect becomes smaller and comparable to the background contributions from elliptic flow if the event plane determined from particle emission angles is used. For the other observable given by the $R(\Delta S)$ correlator related to the distribution of average charge separation in a collision, the signal due to the chiral magnetic effect is found to depend less on whether the reaction or event plane is used in the analysis, and their difference between the two isobaric collision systems is thus a more robust observable.


title: Optimal Boundary Kernels and Weightings for Local Polynomial Regression

abstract: Kernel smoothers are considered near the boundary of the interval. Kernels which minimize the expected mean square error are derived. These kernels are equivalent to using a linear weighting function in the local polynomial regression. It is shown that any kernel estimator that satisfies the moment conditions up to order $m$ is equivalent to a local polynomial regression of order $m$ with some non-negative weight function if and only if the kernel has at most $m$ sign changes. A fast algorithm is proposed for computing the kernel estimate in the boundary region for an arbitrary placement of data points.

categories: stat.ME eess.SP math.ST physics.data-an stat.TH


title: Decoy-state quantum key distribution with a leaky source

abstract: In recent years, there has been a great effort to prove the security of quantum key distribution (QKD) with a minimum number of assumptions. Besides its intrinsic theoretical interest, this would allow for larger tolerance against device imperfections in the actual implementations. However, even in this device-independent scenario, one assumption seems unavoidable, that is, the presence of a protected space devoid of any unwanted information leakage in which the legitimate parties can privately generate, process and store their classical data. In this paper we relax this unrealistic and hardly feasible assumption and introduce a general formalism to tackle the information leakage problem in most of existing QKD systems. More specifically, we prove the security of optical QKD systems using phase and intensity modulators in their transmitters, which leak the setting information in an arbitrary manner. We apply our security proof to cases of practical interest and show key rates similar to those obtained in a perfectly shielded environment. Our work constitutes a fundamental step forward in guaranteeing implementation security of quantum communication systems.


title: Robustness to incorrect system models in stochastic control

abstract: In stochastic control applications, typically only an ideal model (controlled transition kernel) is assumed and the control design is based on the given model, raising the problem of performance loss due to the mismatch between the assumed model and the actual model. Toward this end, we study continuity properties of discrete-time stochastic control problems with respect to system models (i.e., controlled transition kernels) and robustness of optimal control policies designed for incorrect models applied to the true system. We study both fully observed and partially observed setups under an infinite horizon discounted expected cost criterion. We show that continuity and robustness cannot be established under weak and setwise convergences of transition kernels in general, but that the expected induced cost is robust under total variation. By imposing further assumptions on the measurement models and on the kernel itself (such as continuous convergence), we show that the optimal cost can be made continuous under weak convergence of transition kernels as well. Using these continuity properties, we establish convergence results and error bounds due to mismatch that occurs by the application of a control policy which is designed for an incorrectly estimated system model to a true model, thus establishing positive and negative results on robustness.Compared to the existing literature, we obtain strictly refined robustness results that are applicable even when the incorrect models can be investigated under weak convergence and setwise convergence criteria (with respect to a true model), in addition to the total variation criteria. These entail positive implications on empirical learning in (data-driven) stochastic control since often system models are learned through empirical training data where typically weak convergence criterion applies but stronger convergence criteria do not.


title: A Family of Second-Order Energy-Stable Schemes for Cahn-Hilliard Type Equations

abstract: We focus on the numerical approximation of the Cahn-Hilliard type equations, and present a family of second-order unconditionally energy-stable schemes. By reformulating the equation into an equivalent system employing a scalar auxiliary variable, we approximate the system at the time step $(n+\theta)$ ($n$ denoting the time step index and $\theta$ is a real-valued parameter), and devise a family of corresponding approximations that are second-order accurate and unconditionally energy stable. This family of approximations contains the often-used Crank-Nicolson scheme and the second-order backward differentiation formula as particular cases. We further develop an efficient solution algorithm for the resultant discrete system of equations to overcome the difficulty caused by the unknown scalar auxiliary variable. The final algorithm requires only the solution of four de-coupled individual Helmholtz type equations within each time step, which involve only constant and time-independent coefficient matrices that can be pre-computed. A number of numerical examples are presented to demonstrate the performance of the family of schemes developed herein. We note that this family of second-order approximations can be readily applied to devise energy-stable schemes for other types of gradient flows when combined with the auxiliary variable approaches.


title: Identifying and Estimating Principal Causal Effects in Multi-site Trials

abstract: Randomized trials are often conducted with separate randomizations across multiple sites such as schools, voting districts, or hospitals. These sites can differ in important ways, including the site's implementation, local conditions, and the composition of individuals. An important question in practice is whether---and under what assumptions---researchers can leverage this cross-site variation to learn more about the intervention. We address these questions in the principal stratification framework, which describes causal effects for subgroups defined by post-treatment quantities. We show that researchers can estimate certain principal causal effects via the multi-site design if they are willing to impose the strong assumption that the site-specific effects are uncorrelated with the site-specific distribution of stratum membership. We motivate this approach with a multi-site trial of the Early College High School Initiative, a unique secondary education program with the goal of increasing high school graduation rates and college enrollment. Our analyses corroborate previous studies suggesting that the initiative had positive effects for students who would have otherwise attended a low-quality high school, although power is limited.


title: Zero-Shot Object Detection: Learning to Simultaneously Recognize and Localize Novel Concepts

abstract: Current Zero-Shot Learning (ZSL) approaches are restricted to recognition of a single dominant unseen object category in a test image. We hypothesize that this setting is ill-suited for real-world applications where unseen objects appear only as a part of a complex scene, warranting both the `recognition' and `localization' of an unseen category. To address this limitation, we introduce a new \emph{`Zero-Shot Detection'} (ZSD) problem setting, which aims at simultaneously recognizing and locating object instances belonging to novel categories without any training examples. We also propose a new experimental protocol for ZSD based on the highly challenging ILSVRC dataset, adhering to practical issues, e.g., the rarity of unseen objects. To the best of our knowledge, this is the first end-to-end deep network for ZSD that jointly models the interplay between visual and semantic domain information. To overcome the noise in the automatically derived semantic descriptions, we utilize the concept of meta-classes to design an original loss function that achieves synergy between max-margin class separation and semantic space clustering. Furthermore, we present a baseline approach extended from recognition to detection setting. Our extensive experiments show significant performance boost over the baseline on the imperative yet difficult ZSD problem.


title: Sufficient Conditions for a Linear Estimator to be a Local Polynomial Regression

abstract: It is shown that any linear estimator that satisfies the moment conditions up to order $p$ is equivalent to a local polynomial regression of order $p$ with some non-negative weight function if and only if the kernel has at most $p$ sign changes. If the data points are placed symmetrically about the estimation point, a linear weighting function is equivalent to the standard quadratic weighting function.


title: Deep Multiple Instance Learning for Zero-shot Image Tagging

abstract: In-line with the success of deep learning on traditional recognition problem, several end-to-end deep models for zero-shot recognition have been proposed in the literature. These models are successful to predict a single unseen label given an input image, but does not scale to cases where multiple unseen objects are present. In this paper, we model this problem within the framework of Multiple Instance Learning (MIL). To the best of our knowledge, we propose the first end-to-end trainable deep MIL framework for the multi-label zero-shot tagging problem. Due to its novel design, the proposed framework has several interesting features: (1) Unlike previous deep MIL models, it does not use any off-line procedure (e.g., Selective Search or EdgeBoxes) for bag generation. (2) During test time, it can process any number of unseen labels given their semantic embedding vectors. (3) Using only seen labels per image as weak annotation, it can produce a bounding box for each predicted labels. We experiment with the NUS-WIDE dataset and achieve superior performance across conventional, zero-shot and generalized zero-shot tagging tasks.


title: Revealing weak spin-orbit coupling effects on charge carriers in a $\pi$-conjugated polymer

abstract: We measure electrically detected magnetic resonance (EDMR) on organic light-emitting diodes (OLEDs) made of the polymer poly[2-methoxy-5-(2-ethylhexyloxy)-1,4-phenylenevinylene] (MEH-PPV) at room temperature and high magnetic fields, where spectral broadening of the resonance due to spin-orbit coupling (SOC) exceeds that due to the local hyperfine fields. Density-functional-theory calculations on an open-shell model of the material reveal g-tensors of charge-carrier spins in the lowest unoccupied (electron) and highest occupied (hole) molecular orbitals. These tensors are used for simulations of magnetic resonance line-shapes. Besides providing the first quantification and direct observation of SOC effects on charge-carrier states in these weakly SO-coupled hydrocarbons, this procedure demonstrates that spin-related phenomena in these materials are fundamentally monomolecular in nature.


title: The world of research has gone berserk: modeling the consequences of requiring "greater statistical stringency" for scientific publication

abstract: In response to growing concern about the reliability and reproducibility of published science, researchers have proposed adopting measures of greater statistical stringency, including suggestions to require larger sample sizes and to lower the highly criticized p<0.05 significance threshold. While pros and cons are vigorously debated, there has been little to no modeling of how adopting these measures might affect what type of science is published. In this paper, we develop a novel optimality model that, given current incentives to publish, predicts a researcher's most rational use of resources in terms of the number of studies to undertake, the statistical power to devote to each study, and the desirable pre-study odds to pursue. We then develop a methodology that allows one to estimate the reliability of published research by considering a distribution of preferred research strategies. Using this approach, we investigate the merits of adopting measures of `greater statistical stringency' with the goal of informing the ongoing debate.


title: Deep Predictive Coding Neural Network for RF Anomaly Detection in Wireless Networks

abstract: Intrusion detection has become one of the most critical tasks in a wireless network to prevent service outages that can take long to fix. The sheer variety of anomalous events necessitates adopting cognitive anomaly detection methods instead of the traditional signature-based detection techniques. This paper proposes an anomaly detection methodology for wireless systems that is based on monitoring and analyzing radio frequency (RF) spectrum activities. Our detection technique leverages an existing solution for the video prediction problem, and uses it on image sequences generated from monitoring the wireless spectrum. The deep predictive coding network is trained with images corresponding to the normal behavior of the system, and whenever there is an anomaly, its detection is triggered by the deviation between the actual and predicted behavior. For our analysis, we use the images generated from the time-frequency spectrograms and spectral correlation functions of the received RF signal. We test our technique on a dataset which contains anomalies such as jamming, chirping of transmitters, spectrum hijacking, and node failure, and evaluate its performance using standard classifier metrics: detection ratio, and false alarm rate. Simulation results demonstrate that the proposed methodology effectively detects many unforeseen anomalous events in real time. We discuss the applications, which encompass industrial IoT, autonomous vehicle control and mission-critical communications services.


title: Characterization and Optical Properties of Erbium doped As2S3 Films Prepared by Multi-layer Magnetron Sputtering

abstract: As2S3 film doped with erbium is prepared using multi-layer magnetron sputtering. The optical properties were measured by reflectance spectroscopy, and its chemical composition is examined by x-ray photoelectron, Rutherford backscattering, and Raman spectroscopy. The results show that the refractive index and absorption coefficient follow closely to a sputtered As2S3 film, and there are no detectable Er-S clusters and photo-induced As2O3 in the film. Rutherford backscattering spectroscopy shows that the film is homogeneous, and revealed the concentration level of erbium, and the stoichiometry of the film. The deposition method was used to fabricate an integrated Erdoped As2S3 Mach-Zehnder Interferometer and the presence of active erbium ions in the waveguide is evident from the green luminescence it emitted when it was pumped by 1488 nm diode laser. This method is attractive because the doping process can produce an Er:As2S3 film that is close to the ideal stoichiometry of As2S3 with lower risk of photo-decomposed As2O3 crystals developing on the surface when the as-deposited film is exposed to the environment.


title: Between homogeneous and inhomogeneous Navier-Stokes systems: the issue of stability

abstract: We construct large velocity vector solutions to the three dimensional inhomogeneous Navier-Stokes system. The result is proved via the stability of two dimensional solutions with constant density, under the assumption that initial density is point-wisely close to a constant. Key elements of our approach are estimates in the maximal regularity regime and the Lagrangian coordinates. Considerations are done in the whole $\R^3$.


title: The crucial role of ground-based, Doppler measurements for the future of exoplanet science

abstract: We outline the important role that ground-based, Doppler monitoring of exoplanetary systems will play in advancing our theories of planet formation and dynamical evolution. A census of planetary systems requires a well designed survey to be executed over the course of a decade or longer. A coordinated survey to monitor several thousand targets each at ~1000 epochs (~3-5 million new observations) will require roughly 40 dedicated spectrographs. We advocate for improvements in data management, data sharing, analysis techniques, and software testing, as well as possible changes to the funding structures for exoplanet science.


title: Constant-Time Predictive Distributions for Gaussian Processes

abstract: One of the most compelling features of Gaussian process (GP) regression is its ability to provide well-calibrated posterior distributions. Recent advances in inducing point methods have sped up GP marginal likelihood and posterior mean computations, leaving posterior covariance estimation and sampling as the remaining computational bottlenecks. In this paper we address these shortcomings by using the Lanczos algorithm to rapidly approximate the predictive covariance matrix. Our approach, which we refer to as LOVE (LanczOs Variance Estimates), substantially improves time and space complexity. In our experiments, LOVE computes covariances up to 2,000 times faster and draws samples 18,000 times faster than existing methods, all without sacrificing accuracy.


title: Distributed Cache Enabled V2X Networks: Proposals, Research Trends and Challenging Issues

abstract: Nowadays, the internet of vehicles (IoV) has been evolved into the stage of vehicle to everything (V2X). However, the majority of existing work focuses on the motor-vehicles. In contrast, the sharing bicycle system is vastly and rapidly deployed as a feasible internet of things (IoT) application scene for the last mile problem (e.g., from station to home/office). Moreover, the internet access of current V2X is relied on the back-haul to roadside unit (RSU) connections. In this paper, other than prior work, we propose a versatile V2X system with a distributed framework and heterogeneous caching method. All the vehicles and devices on-the-road (motor-vehicle, non-motor-vehicle, pedestrian, etc.) are comprehensively included in the proposed networks. We further introduce a heterogeneous cache method for effective wireless transmission while utilizing the massive connected devices. The potential research trends on achieving high speed transmission, deep programming dedicated network slicing are highlighted as well as the big data, machine learning (ML), fog computing based image recognition and reconstruction, to provide some insights for future studies. Finally, the challenging issues, i.e., urban street canyon path loss and channel models, ultra reliable communication and low latency requirements, are discussed.


title: Characterizing the velocity of a wandering black hole and properties of the surrounding medium using convolutional neural networks

abstract: We present a method for estimating the velocity of a wandering black hole and the equation of state for the gas around, based on a catalog of numerical simulations. The method uses machine learning methods based on convolutional neural networks applied to the classification of images resulting from numerical simulations. Specifically we focus on the supersonic velocity regime and choose the direction of the black hole to be parallel to its spin. We build a catalog of 900 simulations by numerically solving Euler's equations onto the fixed space-time background of a black hole, for two parameters: the adiabatic index $\Gamma$ with values in the range [1.1, 5/3], and the asymptotic relative velocity of the black hole with respect to the surroundings $v_{\infty}$, with values within $[0.2, 0.8]c$. For each simulation we produce a 2D image of the gas density once the process of accretion has approached a stationary regime. The results obtained show that the implemented Convolutional Neural Networks are capable to classify correctly the adiabatic index $87.78\%$ of the time within an uncertainty of $\pm 0.0284$ while the prediction of the velocity is correct $96.67\%$ of the times within an uncertainty of $\pm 0.03c$. We expect that this combination of a massive number of numerical simulations and machine learning methods will help analyze more complicated scenarios related to future high resolution observations of black holes, like those from the Event Horizon Telescope.


title: Inverse Design of Compact Multimode Cavity Couplers

abstract: Efficient coupling between on-chip sources and cavities plays a key role in silicon photonics. However, despite the importance of this basic functionality, there are few systematic design tools to simultaneously control coupling between multiple modes in a compact resonator and a single waveguide. Here, we propose a large-scale adjoint optimization approach to produce wavelength-scale waveguide--cavity couplers operating over tunable and broad frequency bands. We numerically demonstrate couplers discovered by this method that can achieve critical, or nearly critical, coupling between multi-ring cavities and a single waveguide at up to six widely separated wavelengths spanning the $560$--$1500$~nm range of interest for on-chip nonlinear optical devices.


title: Heuristics for vehicle routing problems: Sequence or set optimization?

abstract: We investigate a structural decomposition for the capacitated vehicle routing problem (CVRP) based on vehicle-to-customer "assignment" and visits "sequencing" decision variables. We show that an heuristic search focused on assignment decisions with a systematic optimal choice of sequences (using Concorde TSP solver) during each move evaluation is promising but requires a prohibitive computational effort. We therefore introduce an intermediate search space, based on the dynamic programming procedure of Balas & Simonetti, which finds a good compromise between intensification and computational efficiency. A variety of speed-up techniques are proposed for a fast exploration: neighborhood reductions, dynamic move filters, memory structures, and concatenation techniques. Finally, a tunneling strategy is designed to reshape the search space as the algorithm progresses. The combination of these techniques within a classical local search, as well as in the unified hybrid genetic search (UHGS) leads to significant improvements of solution accuracy. New best solutions are found for surprisingly small instances with as few as 256 customers. These solutions had not been attained up to now with classic neighborhoods. Overall, this research permits to better evaluate the respective impact of sequence and assignment optimization, proposes new ways of combining the optimization of these two decision sets, and opens promising research perspectives for the CVRP and its variants.


title: Simulation of angular resolved RABBITT measurements in noble gas atoms

abstract: We simulate angular resolved RABBITT (Reconstruction of Attosecond Beating By Interference of Two-photon Transitions) measurements on valence shells of noble gas atoms (Ne, Ar, Kr, and Xe). Our non-perturbative numerical simulation is based on solution of the time-dependent Schr\"odinger equation for a target atom driven by an ionizing XUV and dressing IR fields. From these simulations we extract the angular dependent magnitude and phase of the RABBITT oscillations and deduce the corresponding angular anisotropy {\beta} parameter and Wigner time delay ${\tau}_W$ for the single XUV photon absorption which initiates the RABBITT process. Said {\beta} and ${\tau}_W$ parameters are compared with calculations in the random phase approximation with exchange (RPAE) which includes inter-shell correlation. This comparison is used to test various effective potentials employed in the one-electron TDSE. In lighter atoms (Ne and Ar), several effective potentials are found to provide accurate simulation of RABBITT measurements for a wide range of photon energies up to 100 eV above the valence shell threshold. In heavier atoms (Kr and Xe), the onset of strong correlation with the d-shell restricts the validity of the single active electron approximation to several tens of eV above the valence shell threshold.


title: A Meaning-based Statistical English Math Word Problem Solver

abstract: We introduce MeSys, a meaning-based approach, for solving English math word problems (MWPs) via understanding and reasoning in this paper. It first analyzes the text, transforms both body and question parts into their corresponding logic forms, and then performs inference on them. The associated context of each quantity is represented with proposed role-tags (e.g., nsubj, verb, etc.), which provides the flexibility for annotating an extracted math quantity with its associated context information (i.e., the physical meaning of this quantity). Statistical models are proposed to select the operator and operands. A noisy dataset is designed to assess if a solver solves MWPs mainly via understanding or mechanical pattern matching. Experimental results show that our approach outperforms existing systems on both benchmark datasets and the noisy dataset, which demonstrates that the proposed approach understands the meaning of each quantity in the text more.


title: The compression body graph has infinite diameter

abstract: We show that the compression body graph has infinite diameter.


title: Study of Total Electron Content-TEC and electron density profile during geomagnetic storms

abstract: Total Electron Content (TEC) and electron density are the basic parameters, which determine the major properties of the Ionosphere. Detail study of the ionospheric TEC and electron density variations has been carried out during geomagnetic storms, with longitude and latitude, for four different locations: (24{\deg}W-14{\deg}W, 25{\deg}S-10{\deg}S); (53{\deg}W- 46{\deg}W, 04{\deg}N-14{\deg}N); (161{\deg}E-165{\deg}E, 42{\deg}S-34{\deg}S), and (135{\deg}W- 120{\deg}W, 39{\deg}S-35{\deg}S) using the COSMIC satellite data. In order to find the background conditions of the ionosphere, the solar wind parameters such as north-south component of inter planetary magnetic field (Bz), plasma velocity (Vsw), AE, Dst and Kp indices, have also been correlated with the TEC and electron density. The results illustrates that the observed TEC and electron density profile significantly vary with longitudes and latitudes as well. This study illustrates that the values of TEC and the vertical electron density profile are influenced by the solar wind parameters associated with solar activities. The peak value of electron density and TEC increase as the geomagnetic storms becomes stronger. Similarly, the electron density profile vary with altitudes which peaks around the altitude range of about 180-280 km, depending on the strength of geomagnetic storms. The results clearly show that the peak electron density shifted to higher altitude (from about 180 km to 300 km) as the geomagnetic disturbances becomes stronger.


title: Dynamic-structured Semantic Propagation Network

abstract: Semantic concept hierarchy is still under-explored for semantic segmentation due to the inefficiency and complicated optimization of incorporating structural inference into dense prediction. This lack of modeling semantic correlations also makes prior works must tune highly-specified models for each task due to the label discrepancy across datasets. It severely limits the generalization capability of segmentation models for open set concept vocabulary and annotation utilization. In this paper, we propose a Dynamic-Structured Semantic Propagation Network (DSSPN) that builds a semantic neuron graph by explicitly incorporating the semantic concept hierarchy into network construction. Each neuron represents the instantiated module for recognizing a specific type of entity such as a super-class (e.g. food) or a specific concept (e.g. pizza). During training, DSSPN performs the dynamic-structured neuron computation graph by only activating a sub-graph of neurons for each image in a principled way. A dense semantic-enhanced neural block is proposed to propagate the learned knowledge of all ancestor neurons into each fine-grained child neuron for feature evolving. Another merit of such semantic explainable structure is the ability of learning a unified model concurrently on diverse datasets by selectively activating different neuron sub-graphs for each annotation at each step. Extensive experiments on four public semantic segmentation datasets (i.e. ADE20K, COCO-Stuff, Cityscape and Mapillary) demonstrate the superiority of our DSSPN over state-of-the-art segmentation models. Moreoever, we demonstrate a universal segmentation model that is jointly trained on diverse datasets can surpass the performance of the common fine-tuning scheme for exploiting multiple domain knowledge.


title: Memory Slices: A Modular Building Block for Scalable, Intelligent Memory Systems

abstract: While reduction in feature size makes computation cheaper in terms of latency, area, and power consumption, performance of emerging data-intensive applications is determined by data movement. These trends have introduced the concept of scalability as reaching a desirable performance per unit cost by using as few number of units as possible. Many proposals have moved compute closer to the memory. However, these efforts ignored maintaining a balance between bandwidth and compute rate of an architecture, with those of applications, which is a key principle in designing scalable large systems. This paper proposes the use of memory slices, a modular building block for scalable memory systems integrated with compute, in which performance scales with memory size (and volume of data). The slice architecture utilizes a programmable memory interface feeding a systolic compute engine with high reuse rate. The modularity feature of slice-based systems is exploited with a partitioning and data mapping strategy across allocated memory slices where training performance scales with the data size. These features enable shifting the most pressure to cheap compute units rather than expensive memory accesses or transfers via interconnection network. An application of the memory slices to a scale-out memory system is accelerating the training of recurrent, convolutional, and hybrid neural networks (RNNs and RNNs+CNN) that are forming cloud workloads. The results of our cycle-level simulations show that memory slices exhibits a superlinear speedup when the number of slices increases. Furthermore, memory slices improve power efficiency to 747 GFLOPs/J for training LSTMs. While our current evaluation uses memory slices with 3D packaging, a major value is that slices can also be constructed with a variety of packaging options, for example with DDR-based memory units.


title: A Preliminary Analysis of the Shangri-La Bolide on 2017 Oct 4

abstract: At 12:07 UT (8:07 pm China Standard Time) on 2017 Oct 4, a bright bolide was widely observed in the Shangri-La region in the Province of Yunnan, China (Figure 1). The event was well observed by the general public as it took place on the night of the Mid Autumn Festival which associates with moon gazing. Sonic booms and ground shaking were reported in an area about a thousand square kilometers wide to the northwest of the Shangri-La City. Data from the U.S. government sensor suggested that the impact energy of the event is approximately 0.54 kt TNT equivalent, with the terminus of the bolide positioned at $28.1^\circ$ N, $99.4^\circ$ E. This is the largest observed bolide event overland since the bolide event took place in Mauritania on 2016 Jun 27 (1.2 kt). In this Research Note we present a preliminary analysis of this event.


title: Modelling sparsity, heterogeneity, reciprocity and community structure in temporal interaction data

abstract: We propose a novel class of network models for temporal dyadic interaction data. Our goal is to capture a number of important features often observed in social interactions: sparsity, degree heterogeneity, community structure and reciprocity. We propose a family of models based on self-exciting Hawkes point processes in which events depend on the history of the process. The key component is the conditional intensity function of the Hawkes Process, which captures the fact that interactions may arise as a response to past interactions (reciprocity), or due to shared interests between individuals (community structure). In order to capture the sparsity and degree heterogeneity, the base (non time dependent) part of the intensity function builds on compound random measures following Todeschini et al. (2016). We conduct experiments on a variety of real-world temporal interaction data and show that the proposed model outperforms many competing approaches for link prediction, and leads to interpretable parameters.


title: Impacts of Dirty Data: and Experimental Evaluation

abstract: Data quality issues have attracted widespread attention due to the negative impacts of dirty data on data mining and machine learning results. The relationship between data quality and the accuracy of results could be applied on the selection of the appropriate algorithm with the consideration of data quality and the determination of the data share to clean. However, rare research has focused on exploring such relationship. Motivated by this, this paper conducts an experimental comparison for the effects of missing, inconsistent and conflicting data on classification and clustering algorithms. Based on the experimental findings, we provide guidelines for algorithm selection and data cleaning.


title: Computing Special $L$-Values of Certain Modular Forms with Complex Multiplication

abstract: In this expository paper, we illustrate two explicit methods which lead to special $L$-values of certain modular forms admitting complex multiplication (CM), motivated in part by properties of $L$-functions obtained from Calabi-Yau manifolds defined over $\mathbb Q$.


title: Lyapunov Functions for First-Order Methods: Tight Automated Convergence Guarantees

abstract: We present a novel way of generating Lyapunov functions for proving linear convergence rates of first-order optimization methods. Our approach provably obtains the fastest linear convergence rate that can be verified by a quadratic Lyapunov function (with given states), and only relies on solving a small-sized semidefinite program. Our approach combines the advantages of performance estimation problems (PEP, due to Drori & Teboulle (2014)) and integral quadratic constraints (IQC, due to Lessard et al. (2016)), and relies on convex interpolation (due to Taylor et al. (2017c;b)).


title: Reconfiguring spanning and induced subgraphs

abstract: Subgraph reconfiguration is a family of problems focusing on the reachability of the solution space in which feasible solutions are subgraphs, represented either as sets of vertices or sets of edges, satisfying a prescribed graph structure property. Although there has been previous work that can be categorized as subgraph reconfiguration, most of the related results appear under the name of the property under consideration; for example, independent set, clique, and matching. In this paper, we systematically clarify the complexity status of subgraph reconfiguration with respect to graph structure properties.


title: Formal Analysis of Non-functional Properties for a Cooperative Automotive System

abstract: Modeling and analysis of nonfunctional requirements is crucial in automotive systems. EAST-ADL is an architectural language dedicated to safety-critical automotive system design. We have previously modified EAST-ADL to include energy constraints and transformed energy-aware timed (ET) behaviors modeled in SIMULINK/STATEFLOW into UPPAAL models amenable to formal verification. Previous work is extended in this paper by including support for SIMULINK DESIGN VERIFIER (SDV), i.e., the ET constraints are translated into proof objective models that can be verified using SDV. Furthermore, probabilistic extension of EAST-ADL constraints is defined and the semantics of the extended constraints is translated into verifiable UPPAAL models with stochastic semantics for formal verification. A set of mapping rules are proposed to facilitate the guarantee of translation. Verification & Validation are performed on the extended timing and energy constraints using SDV and UPPAAL-SMC. Our approach is demonstrated on a cooperative automotive system case study.


title: Renewable Energy Integration in Distribution System -- Synchrophasor Sensor based Big Data Analysis, Visualization, and System Operation

abstract: Due to the large volume of heterogeneous data provided by both the customer and the grid side, a big data visualization platform is built to discover the hidden useful knowledge for smart grid (SG) operation, control and situation awareness. An open source cluster computing framework based on Apache Spark is considering to discover the hidden knowledge of the bag-data. And the data is transmitted with a Open System Interconnection (OSI) model to the data visualization platform in a high-speed communication architecture. Google Earth and global geographic information system (GIS) are used to design the visualization platform and realize the result with the test bench. The alternating direction method of multipliers (ADMM) is used to compute the optimal power flow in distributed manner. The proposed network reconfiguration is solved in parallel manner with the limited switches and the strong computing capability. Further more, a multi-timescale operation approach is modeled with a three-phase distributed system, which consists the hourly scheduling at substation level and the minutes power flow operation at feeder level. The system cost with renewable generation is minimized at the substation level. The given error distribution model of renewable generation is simulated with a chance constraint, and the derived deterministic form is modeled with Gaussian mixture model (GMM) with genetic algorithm-based expectationmaximization (GAEM). The system cost is further reduced with the OPF in real-time (RT) scheduling. The semidefinite programming (SDP) is used to relax the nonconvexity of the three-phase unbalanced distribution system into a convex problem, which make sure to achieve the global optimal result. In the parallel manner, the ADMM is realizing getting the results in a short time.


title: Real-time Detection, Tracking, and Classification of Moving and Stationary Objects using Multiple Fisheye Images

abstract: The ability to detect pedestrians and other moving objects is crucial for an autonomous vehicle. This must be done in real-time with minimum system overhead. This paper discusses the implementation of a surround view system to identify moving as well as static objects that are close to the ego vehicle. The algorithm works on 4 views captured by fisheye cameras which are merged into a single frame. The moving object detection and tracking solution uses minimal system overhead to isolate regions of interest (ROIs) containing moving objects. These ROIs are then analyzed using a deep neural network (DNN) to categorize the moving object. With deployment and testing on a real car in urban environments, we have demonstrated the practical feasibility of the solution. The video demos of our algorithm have been uploaded to Youtube: https://youtu.be/vpoCfC724iA, https://youtu.be/2X4aqH2bMBs


title: Sampling Conditions for Conforming Voronoi Meshing by the VoroCrust Algorithm

abstract: We study the problem of decomposing a volume bounded by a smooth surface into a collection of Voronoi cells. Unlike the dual problem of conforming Delaunay meshing, a principled solution to this problem for generic smooth surfaces remained elusive. VoroCrust leverages ideas from $\alpha$-shapes and the power crust algorithm to produce unweighted Voronoi cells conforming to the surface, yielding the first provably-correct algorithm for this problem. Given an $\epsilon$-sample on the bounding surface, with a weak $\sigma$-sparsity condition, we work with the balls of radius $\delta$ times the local feature size centered at each sample. The corners of this union of balls are the Voronoi sites, on both sides of the surface. The facets common to cells on opposite sides reconstruct the surface. For appropriate values of $\epsilon$, $\sigma$ and $\delta$, we prove that the surface reconstruction is isotopic to the bounding surface. With the surface protected, the enclosed volume can be further decomposed into an isotopic volume mesh of fat Voronoi cells by generating a bounded number of sites in its interior. Compared to state-of-the-art methods based on clipping, VoroCrust cells are full Voronoi cells, with convexity and fatness guarantees. Compared to the power crust algorithm, VoroCrust cells are not filtered, are unweighted, and offer greater flexibility in meshing the enclosed volume by either structured grids or random samples.


title: Electron hole instability in linearly sub-critical plasmas

abstract: Electron holes (EH) are highly stable non-linear structures met omnipresently in driven collision-less hot plasmas. A mechanism destabilizing small perturbations into holes is essential for an often witnessed but less understood sub-critically driven intermittent plasma turbulence. In this paper we show how a tiny, eddy-like, non-topological seed fluctuation can trigger an unstable evolution deep in the linearly damped region, a process being controlled by the trapping non-linearity and hence being beyond the realm of the Landau scenario. After a (transient) transition phase modes of the privileged spectrum of cnoidal EH are excited which in the present case consist of a solitary electron hole (SEH), two counter-propagating "Langmuir" modes (plasma oscillation), and an ion acoustic mode. A quantitative explanation involves employing non-linear eigen-modes, yielding a non-linear dispersion relation with a forbidden regime and the negative energy character of the SEH, properties being inherent in Schamel's model of undamped Vlasov-Poisson structures identified here as lowest order trapped particle equilibria. An important role in the final adaption of nonlinear plasma eigenmodes is played by a deterministic response of trapped electrons which facilitates transfer of energy from electron thermal energy to an ion acoustic non-uniformity, accelerating the SEH and positioning it into the right place assigned by the theory.


title: K-Theory of Hilbert Schemes as a Formal Quantum Field Theory

abstract: We define a notion of formal quantum field theory and associate a formal quantum field theory to K-theoretical intersection theories on Hilbert schemes of points on algebraic surfaces. This enables us to find an effective way to compute K-theoretical intersection theories on Hilbert schemes via a connection to Macdonald polynomials and vertex operators.


title: Two-Qubit Stabilizer Circuits with Recovery I: Existence

abstract: Understanding how a stabilizer circuit responds to different input is important to formulating an effective strategy for resource management. In this paper, we further investigate the many ways of using stabilizer operations to generate a single qubit output from a two-qubit state. In particular, by restricting the input to certain product states, we discover probabilistic operations capable of transforming stabilizer circuit outputs back into stabilizer circuit inputs. These secondary operations are ideally suited for recovery purposes and require only one extra resource ancilla to succeed. As a result of reusing qubits in this manner, we present an alternative to the original state preparation process that can lower the overall costs of executing any two-qubit stabilizer procedure involving non-stabilizer resources.


title: Load Balancing for 5G Ultra-Dense Networks using Device-to-Device Communications

abstract: Load balancing is an effective approach to address the spatial-temporal fluctuation problem of mobile data traffic for cellular networks. The existing schemes that focus on channel borrowing from neighboring cells cannot be directly applied to future 5G wireless networks, because the neighboring cells will reuse the same spectrum band in 5G systems. In this paper, we consider an orthogonal frequency division multiple access~(OFDMA) ultra-dense small cell network, where Device-to-Device~(D2D) communication is advocated to facilitate load balancing without extra spectrum. Specifically, the data traffic can be effectively offloaded from a congested small cell to other underutilized small cells by D2D communications. The problem is naturally formulated as a joint resource allocation and D2D routing problem that maximizes the system sum-rate. To efficiently solve the problem, we decouple the problem into a resource allocation subproblem and a D2D routing subproblem. The two subproblems are solved iteratively as a monotonic optimization problem and a complementary geometric programming problem, respectively. Simulation results show that the data sum-rate in the neighboring small cells increases 20\% on average by offloading the data traffic in the congested small cell to the neighboring small cell base stations~(SBSs).


title: On isomorphisms between weighted $L^p$-algebras

abstract: Let $G$ be a locally compact group and $\omega$ be a continuous weight on $G$. In this paper, we first characterize bicontinuous biseparating algebra isomorphisms between weighted $L^p$-algebras. As a result we extend previous results of Edwards, Strichartz and Parrot on algebra isomorphisms between $L^p$-algebras to the setting of weighted $L^p$-algebras. We then study the automorphisms of certain weighted $L^p$-algebras on integers, applying known results on composition operators on classical function spaces.


title: A Kernel Theory of Modern Data Augmentation

abstract: Data augmentation, a technique in which a training set is expanded with class-preserving transformations, is ubiquitous in modern machine learning pipelines. In this paper, we seek to establish a theoretical framework for understanding data augmentation. We approach this from two directions: First, we provide a general model of augmentation as a Markov process, and show that kernels appear naturally with respect to this model, even when we do not employ kernel classification. Next, we analyze more directly the effect of augmentation on kernel classifiers, showing that data augmentation can be approximated by first-order feature averaging and second-order variance regularization components. These frameworks both serve to illustrate the ways in which data augmentation affects the downstream learning model, and the resulting analyses provide novel connections between prior work in invariant kernels, tangent propagation, and robust optimization. Finally, we provide several proof-of-concept applications showing that our theory can be useful for accelerating machine learning workflows, such as reducing the amount of computation needed to train using augmented data, and predicting the utility of a transformation prior to training.


title: High-Resolution Observations of Flares in an Arch Filament System

abstract: We study five sequential solar flares (SOL2015-08-07) occurring in Active Region 12396 observed with the Goode Solar Telescope (GST) at the BBSO, complemented by IRIS and SDO observations. The main flaring region is an arch filament system (AFS) consisting of multiple bundles of dark filament threads enclosed by semi-circular flare ribbons. We study the magnetic configuration and evolution of the active region by constructing coronal magnetic field models based on SDO/HMI magnetograms using two independent methods, i.e., the nonlinear force-free field (NLFFF) extrapolation and the flux rope insertion method. The models consist of multiple flux ropes with mixed signs of helicity, i.e., positive (negative) in the northern (southern) region, which is consistent with the GST observations of multiple filament bundles. The footprints of quasi-separatrix layers (QSLs) derived from the extrapolated NLFFF compare favorably with the observed flare ribbons. An interesting double-ribbon fine structure located at the east border of the AFS is consistent with the fine structure of the QSL's footprint. Moreover, magnetic field lines traced along the semi-circular footprint of a dome-like QSL surrounding the AFS are connected to the regions of significant helicity and Poynting flux injection. The maps of magnetic twist show that positive twist became dominant as time progressed, which is consistent with the injection of positive helicity before the flares. We hence conclude that these circular shaped flares are caused by 3D magnetic reconnection at the QSLs associated with the AFS possessing mixed signs of helicity.


title: Weak units, universal cells, and coherence via universality for bicategories

abstract: Poly-bicategories generalise planar polycategories in the same way as bicategories generalise monoidal categories. In a poly-bicategory, the existence of enough 2-cells satisfying certain universal properties (representability) induces coherent algebraic structure on the 2-graph of single-input, single-output 2-cells. A special case of this theory was used by Hermida to produce a proof of strictification for bicategories. No full strictification is possible for higher-dimensional categories, seemingly due to problems with 2-cells that have degenerate boundaries; it was conjectured by C. Simpson that semi-strictification excluding units may be possible. We study poly-bicategories where 2-cells with degenerate boundaries are barred, and show that we can recover the structure of a bicategory through a different construction of weak units. We prove that the existence of these units is equivalent to the existence of 1-cells satisfying lower-dimensional universal properties, and study the relation between preservation of units and universal cells. Then, we introduce merge-bicategories, a variant of poly-bicategories with more composition operations, which admits a natural monoidal closed structure, giving access to higher morphisms. We derive equivalences between morphisms, transformations, and modifications of representable merge-bicategories and the corresponding notions for bicategories. Finally, we prove a semi-strictification theorem for representable merge-bicategories with a choice of composites and units.


title: A Globally Asymptotically Stable Polynomial Vector Field with Rational Coefficients and no Local Polynomial Lyapunov Function

abstract: We give an explicit example of a two-dimensional polynomial vector field of degree seven that has rational coefficients, is globally asymptotically stable, but does not admit an analytic Lyapunov function even locally.


title: High resolution observations of sympathetic eruptions by NVST

abstract: We investigate two sympathetic filament eruptions observed by the New Vacuum Solar Telescope (NVST) on 2015 October 15. The full picture of the eruptions is obtained from the corresponding SDO/AIA observations. The two filaments start from active region NOAA 12434 in the north and end in one large quiescent filament channel in the south. The left filament erupts firstly, followed by the right filament eruption about 10 minutes later. Clear twist structure and rotating motion are observed in both filaments during the eruption. Both eruptions are failed, since the filaments firstly rise up, then flow towards the south and merge into the southern large quiescent filament. We also observe repeating activations of mini filaments below the right filament after its eruption. Using magnetic field models constructed based on SDO/HMI magnetograms by flux rope insertion method, we find that the left filament eruption is likely to be triggered by kink instability, while weakening of overlying magnetic fields due to magnetic reconnection at an X-point between the two filament systems might play an important role in the onset of the right filament eruption.


title: Distributed Caching for Complex Querying of Raw Arrays

abstract: As applications continue to generate multi-dimensional data at exponentially increasing rates, fast analytics to extract meaningful results is becoming extremely important. The database community has developed array databases that alleviate this problem through a series of techniques. In-situ mechanisms provide direct access to raw data in the original format---without loading and partitioning. Parallel processing scales to the largest datasets. In-memory caching reduces latency when the same data are accessed across a workload of queries. However, we are not aware of any work on distributed caching of multi-dimensional raw arrays. In this paper, we introduce a distributed framework for cost-based caching of multi-dimensional arrays in native format. Given a set of files that contain portions of an array and an online query workload, the framework computes an effective caching plan in two stages. First, the plan identifies the cells to be cached locally from each of the input files by continuously refining an evolving R-tree index. In the second stage, an optimal assignment of cells to nodes that collocates dependent cells in order to minimize the overall data transfer is determined. We design cache eviction and placement heuristic algorithms that consider the historical query workload. A thorough experimental evaluation over two real datasets in three file formats confirms the superiority -- by as much as two orders of magnitude -- of the proposed framework over existing techniques in terms of cache overhead and workload execution time.


title: A falling Corona model for the anomalous behavior of the broad emission lines in NGC 5548

abstract: NGC 5548 has been intensively monitored by the AGN Space Telescope and Optical Reverberation Mapping collaboration. Approximately after half of the light curves, the correlation between the broad emission lines and the lag-corrected ultraviolet continua becomes weak. This anomalous behavior is accompanied by an increase of soft X-ray emission. We propose a simple model to understand this anomalous behavior, i.e., the corona might fall down, thereby increasing the covering fraction of the inner disk. Therefore, X-ray and extreme ultraviolet emission suffer from spectral variations. The ultraviolet continua variations are driven by both X-ray and extreme ultraviolet variations. Consequently, the spectral variability induced by the falling corona would dilute the correlation between the broad emission lines and the ultraviolet continua. Our model can explain many additional observational facts, including the dependence of the anomalous behavior on velocity and ionization energy. We also show that the time lag and correlation between the X-ray and the ultraviolet variations change as NGC 5548 displays the anomalous behavior. The time lag is dramatically longer than the expectation from disk reprocessing if the anomalous behavior is properly excluded. During the anomalous state, the time lag approaches the light-travel timescale of disk reprocessing albeit with a much weaker correlation. We speculate that the time lag in the normal state is caused by reprocessing of the broad line region gas. As NGC 5548 enters the abnormal state, the contribution of the broad line region gas is smaller; the time lag reflects disk reprocessing. We also discuss alternative scenarios.


title: Salient Objects in Clutter: Bringing Salient Object Detection to the Foreground

abstract: We provide a comprehensive evaluation of salient object detection (SOD) models. Our analysis identifies a serious design bias of existing SOD datasets which assumes that each image contains at least one clearly outstanding salient object in low clutter. The design bias has led to a saturated high performance for state-of-the-art SOD models when evaluated on existing datasets. The models, however, still perform far from being satisfactory when applied to real-world daily scenes. Based on our analyses, we first identify 7 crucial aspects that a comprehensive and balanced dataset should fulfill. Then, we propose a new high quality dataset and update the previous saliency benchmark. Specifically, our SOC (Salient Objects in Clutter) dataset, includes images with salient and non-salient objects from daily object categories. Beyond object category annotations, each salient image is accompanied by attributes that reflect common challenges in real-world scenes. Finally, we report attribute-based performance assessment on our dataset.


title: A Dataset and Architecture for Visual Reasoning with a Working Memory

abstract: A vexing problem in artificial intelligence is reasoning about events that occur in complex, changing visual stimuli such as in video analysis or game play. Inspired by a rich tradition of visual reasoning and memory in cognitive psychology and neuroscience, we developed an artificial, configurable visual question and answer dataset (COG) to parallel experiments in humans and animals. COG is much simpler than the general problem of video analysis, yet it addresses many of the problems relating to visual and logical reasoning and memory -- problems that remain challenging for modern deep learning architectures. We additionally propose a deep learning architecture that performs competitively on other diagnostic VQA datasets (i.e. CLEVR) as well as easy settings of the COG dataset. However, several settings of COG result in datasets that are progressively more challenging to learn. After training, the network can zero-shot generalize to many new tasks. Preliminary analyses of the network architectures trained on COG demonstrate that the network accomplishes the task in a manner interpretable to humans.


title: Holomorphic sectional curvature, nefness and Miyaoka-Yau type inequality

abstract: On a compact K\"ahler manifold, we introduce a notion of almost nonpositivity for the holomorphic sectional curvature, which by definition is weaker than the existence of a K\"ahler metric with semi-negative holomorphic sectional curvature. We prove that a compact K\"ahler manifold of almost nonpositive holomorphic sectional curvature has a nef canonical line bundle, contains no rational curves and satisfies some Miyaoka-Yau type inequalities. In the course of the discussions, we attach a real value to any fixed K\"ahler class which, up to a constant factor depending only on the dimension of manifold, turns out to be an upper bound for the nef threshold.


title: Turbulence, Accretion Braking Torques and Efficient Jets Without Magnetocentrifugal Acceleration: Core Concepts

abstract: I discuss three mutually-supportive notions or assumptions regarding jets and accretion. The first is magnetocentrifugal acceleration (MCA), the overwhelmingly favored mechanism for the production of jets in most steady accreting systems. The second is the zero-torque inner boundary condition. The third is that effective viscous dissipation is like real dissipation, leading directly to heating. All three assumptions fit nicely together in a manner that is simple, persuasive, and mutually-consistent. All, I argue, are incorrect. For concreteness I focus on protostars. Magnetohydrodynamic (MHD) turbulence in accretion is not a sink of energy, but a reservoir, capable of doing mechanical work directly and therefore efficiently, rather than solely through ohmic ("viscous") heating. Advection of turbulence energy reduces the effective radiative efficiency, and may help solve the missing boundary-layer emission problem. The angular momentum problem, whereby accretion spins up a protostar to breakup, is resolved by allowing direct viscous coupling to the protostar, permitting substantially greater energy to be deposited into the accretion flow than otherwise possible. This goes not into heat, but into a turbulent, tangled, buoyant toroidal magnetic field. I argue that there is neither an angular momentum problem nor an efficiency problem that MCA is needed to solve. Moreover, the turbulent magnetic field has ample ability not just to collimate but to accelerate gas, first radially inwards through tension forces and then vertically through pressure forces, without any MCA mechanism. I suggest then that jets, particularly the most powerful and well-collimated protostellar jets, are not magnetocentrifugally driven.


title: On the Iwasawa asymptotic class number formula for $\mathbb{Z}_p^r\rtimes\mathbb{Z}_p$-extensions

abstract: Let $p$ be an odd prime and $F_{\infty,\infty}$ a $p$-adic Lie extension of a number field $F$ with Galois group isomorphic to $\mathbb{Z}_p^r\rtimes\mathbb{Z}_p$, $r\geq 1$. Under certain assumptions, we prove an asymptotic formula for the growth of $p$-exponents of the class groups in the said $p$-adic Lie extension. This generalizes a previous result of Lei, where he establishes such a formula in the case $r=1$. An important and new ingredient towards extending Lei's result rests on an asymptotic formula for a finitely generated (not necessarily torsion) $\mathbb{Z}_p[[\mathbb{Z}_p^r]]$-module which we will also establish in this paper. We then continue studying the growth of $p$-exponents of the class groups under more restrictive assumptions and show that there is an asymptotic formula in our noncommutative $p$-adic Lie extension analogous to a refined formula of Monsky (which is for the commutative extension) in a special case.


title: Expected Time to Extinction of SIS Epidemic Model Using Quasy Stationary Distribution

abstract: We study that the breakdown of epidemic depends on some parameters, that is expressed in epidemic reproduction ratio number. It is noted that when $R_0 $ exceeds 1, the stochastic model have two different results. But, eventually the extinction will be reached even though the major epidemic occurs. The question is how long this process will reach extinction. In this paper, we will focus on the Markovian process of SIS model when major epidemic occurs. Using the approximation of quasi--stationary distribution, the expected mean time of extinction only occurs when the process is one step away from being extinct. Combining the theorm from Ethier and Kurtz, we use CLT to find the approximation of this quasi distribution and successfully determine the asymptotic mean time to extinction of SIS model without demography.


title: Evolution of superconductivity in ultrathin NbS2

abstract: We report a systematic study of thickness-dependent superconductivity and carrier transport properties in exfoliated layered 2H-NbS2. Hall-effect measurements reveal 2H-NbS2 in its normal state to be a p-type metal with hole mobility of 1-3 cm2/Vs. The superconducting transition temperature is found to decrease with thickness. We find that the suppression of superconductivity is due to disorder resulting from the incorporation of atmospheric oxygen and a reduced hole density. Cross-section transmission electron microscope (TEM) imaging reveals a chemical change of NbS2 in ambient conditions, resulting in the formation of amorphous oxide layers sandwiching crystalline layered NbS2. Though few-nm-thick 2H-NbS2 completely converts to amorphous oxide in ambient conditions, PMMA encapsulation prevents further chemical change and preserves superconductivity.


title: Skyrmions in models with pions and rho mesons

abstract: A problem with the standard Skyrme model is that Skyrmion binding energies are around $15\%$, being much larger than the order $1\%$ binding energies of the nuclei that they aim to describe. Here we consider theories that extend the standard Skyrme model of pions by including rho mesons, via dimensional deconstruction of Yang-Mills theory with an extra dimension. We report the first results of parallel numerical computations of multi-Skyrmions in theories of this type, including a model that reduces Skyrmion energies below those of the standard Skyrme model whilst retaining exactly the same Faddeev-Bogomolny energy bound. We compute all Skyrmions with baryons numbers up to 12 and find that the inclusion of rho mesons reduces binding energies to less than $4\%,$ and therefore moves Skyrmion theory closer to experimental data. Furthermore, we find that this dramatic reduction in binding energies is obtained without changing the qualitative features of the Skyrmions, such as their symmetries.


title: Kinematic dynamo action of a precession driven flow based on the results of water experiments and hydrodynamic simulations

abstract: The project DRESDYN (DREsden Sodium facility for DYNamo and thermohydraulic studies) conducted at Helmholtz-Zentrum Dresden-Rossendorf (HZDR) provides a new platform for a variety of liquid sodium experiments devoted to problems of geo- and astrophysical magnetohydrodynamics. The most ambitious experiment within this project is a precession driven dynamo experiment that currently is under construction. It consists of a cylinder filled with liquid sodium that simultaneously rotates around two axes. The experiment is motivated by the idea of a precession-driven flow as a complementary energy source for the geodynamo or the ancient lunar dynamo. In the present study we address numerical and experimental examinations in order to identify parameter regions where the onset of magnetic field excitation will be most probable. Both approaches show that in the strongly nonlinear regime the flow is essentially composed of the directly forced primary Kelvin mode and higher modes in terms of standing inertial waves that arise from nonlinear self-interactions. A peculiarity is the resonance-like emergence of an axisymmetric mode that represents a double roll structure in the meridional plane, which, however, only occurs in a very limited range of the precession ratio. This axisymmetric mode turns out to be beneficial for dynamo action, and kinematic simulations of the magnetic field evolution induced by the time-averaged flow exhibit magnetic field excitation at critical magnetic Reynolds numbers around ${\rm{Rm}}^{\rm{c}}\approx 430$, which is well within the range of the planned liquid sodium experiment.


title: Coexisting localized and itinerant gapless excitations in a quantum spin liquid candidate 1T-TaS$_2$

abstract: To reveal the nature of elementary excitations in a quantum spin liquid (QSL), we measured low temperature thermal conductivity and specific heat of 1T-TaS$_2$, a QSL candidate material with frustrated triangular lattice of spin-1/2. The nonzero temperature linear specific heat coefficient $\gamma$ and the finite residual linear term of the thermal conductivity in the zero temperature limit $\kappa_0/T=\kappa/T(T\rightarrow 0)$ are clearly resolved. This demonstrates the presence of highly mobile gapless excitations, which is consistent with fractionalized spinon excitations that form a Fermi surface. Remarkably, an external magnetic field strongly suppresses $\gamma$, whereas it enhances $\kappa_0/T$. This unusual contrasting behavior in the field dependence of specific heat and thermal conductivity can be accounted for by the presence of two types of gapless excitations with itinerant and localized characters, as recently predicted theoretically (I. Kimchi et al., arXiv:1803.00013 (2018)). This unique feature of 1T-TaS$_2$ provides new insights into the influence of quenched disorder on the QSL.


title: Tractability properties of the weighted star discrepancy of the Halton sequence

abstract: We study the weighted star discrepancy of the Halton sequence. In particular, we show that the Halton sequence achieves strong polynomial tractability for the weighted star discrepancy for product weights $(\gamma_j)_{j \ge 1}$ under the mildest condition on the weight sequence known so far for explicitly constructive sequences. The condition requires $\sup_{d \ge 1} \max_{\emptyset \not= \mathfrak{u} \subseteq [d]} \prod_{j \in \mathfrak{u}} (j \gamma_j) < \infty$. The same result holds for Niederreiter sequences and for other types of digital sequences. Our results are true also for the weighted unanchored discrepancy.


title: Parameterized Low-Rank Binary Matrix Approximation

abstract: We provide a number of algorithmic results for the following family of problems: For a given binary m\times n matrix A and integer k, decide whether there is a "simple" binary matrix B which differs from A in at most k entries. For an integer r, the "simplicity" of B is characterized as follows. - Binary r-Means: Matrix B has at most r different columns. This problem is known to be NP-complete already for r=2. We show that the problem is solvable in time 2^{O(k\log k)}\cdot(nm)^{O(1)} and thus is fixed-parameter tractable parameterized by k. We prove that the problem admits a polynomial kernel when parameterized by r and k but it has no polynomial kernel when parameterized by k only unless NP\subseteq coNP/poly. We also complement these result by showing that when being parameterized by r and k, the problem admits an algorithm of running time 2^{O(r\cdot \sqrt{k\log{(k+r)}})}(nm)^{O(1)}, which is subexponential in k for r\in O(k^{1/2 -\epsilon}) for any \epsilon>0. - Low GF(2)-Rank Approximation: Matrix B is of GF(2)-rank at most r. This problem is known to be NP-complete already for r=1. It also known to be W[1]-hard when parameterized by k. Interestingly, when parameterized by r and k, the problem is not only fixed-parameter tractable, but it is solvable in time 2^{O(r^{ 3/2}\cdot \sqrt{k\log{k}})}(nm)^{O(1)}, which is subexponential in k. - Low Boolean-Rank Approximation: Matrix B is of Boolean rank at most r. The problem is known to be NP-complete for k=0 as well as for r=1. We show that it is solvable in subexponential in k time 2^{O(r2^r\cdot \sqrt{k\log k})}(nm)^{O(1)}.


title: Model-based Verification and Validation of an Autonomous Vehicle System

abstract: The software development for Cyber-Physical Systems (CPS), e.g., autonomous vehicles, requires both functional and non-functional quality assurance to guarantee that the CPS operates safely and effectively. EAST-ADL is a domain specific architectural language dedicated to safety-critical automotive embedded system design. We have previously modified EAST-ADL to include energy constraints and transformed energy-aware real-time (ERT) behaviors modeled in EAST-ADL/STATEFLOW into UPPAAL models amenable to formal verification. Previous work is extended in this paper by including support for SIMULINK and an integration of Simulink/Stateflow within a same tool-chain. Simulink/Stateflow models are transformed, based on extended ERT constraints in EAST-ADL, into verifiable UPPAAL models with stochastic semantics and integrate the translation with formal statistical analysis techniques: Probabilistic extension of EAST-ADL constraints is defined as a semantics denotation. A set of mapping rules is proposed to facilitate the guarantee of translation. Formal analysis on both functional- and non-functional properties is performed using SIMULINK DESIGN VERIFIER/UPPAAL-SMC. The analysis techniques are validated and demonstrated on the autonomous traffic sign recognition vehicle case study.


title: The growth of interest for astronomical X-ray polarimetry

abstract: Astronomical X-ray polarimetry was first explored in the end of the 60's by pioneering rocket instruments. The craze arising from the first discoveries on stellar and supernova remnant X-ray polarization led to the addition of X-ray polarimeters on-board of early satellites. Unfortunately, the inadequacy of the diffraction and scattering technologies required to measure polarization with respect to the constraints driven by X-ray mirrors and detectors, coupled to long integration times, slowed down the field for almost 40 years. Thanks to the development of new, highly sensitive, compact X-ray polarimeters in the beginning of the 2000's, the possibility to observe astronomical X-ray polarization is rising again and scientists are now ready to explore the high energy sky thanks to modern X-ray polarimeters. In the forthcoming years, several X-ray missions (both rockets, balloons and satellites) will open a new observational windows. A wind of renewal blows over the area of X-ray polarimetry and this paper presents for the first time a quantitative assessment, all based on scientific literature, of the growth of interest for astronomical X-ray polarimetry.


title: Efficient Decoding Schemes for Noisy Non-Adaptive Group Testing when Noise Depends on Number of Items in Test

abstract: The goal of non-adaptive group testing is to identify at most $d$ defective items from $N$ items, in which a test of a subset of $N$ items is positive if it contains at least one defective item, and negative otherwise. However, in many cases, especially in biological screening, the outcome is unreliable due to biochemical interaction; i.e., \textit{noise.} Consequently, a positive result can change to a negative one (false negative) and vice versa (false positive). In this work, we first consider the dilution effect in which \textit{the degree of noise depends on the number of items in the test}. Two efficient schemes are presented for identifying the defective items in time linearly to the number of tests needed. Experimental results validate our theoretical analysis. Specifically, setting the error precision of 0.001 and $d\leq16$, our proposed algorithms always identify all defective items in less than 7 seconds for $N=2^{33}\approx 9$ billion.


title: An Irreducible Polynomial Functional Basis of Two-dimensional Eshelby Tensors

abstract: Representation theorems for both isotropic and anisotropic functions are of prime importance in both theoretical and applied mechanics. The Eshelby inclusion problem is very fundamental, and is of particular importance in the design of advanced functional composite materials. In this paper, we discuss about two-dimensional Eshelby tensors (denoted as $M^{(2)}$). Eshelby tensors satisfy the minor index symmetry $M_{ijkl}^{(2)}=M_{jikl}^{(2)}=M_{ijlk}^{(2)}$ and have wide applications in many fields of mechanics. In view of the representation of two-dimensional irreducible tensors in complex field, we obtain a minimal integrity basis of ten isotropic invariants of $M^{(2)}$. Remarkably, note that an integrity basis is always a functional basis, we further confirm that the minimal integrity basis is also an irreducible function basis of isotropic invariants of $M^{(2)}$.


title: Lipschitz Constrained GANs via Boundedness and Continuity

abstract: One of the challenges in the study of Generative Adversarial Networks (GANs) is the difficulty of its performance control. Lipschitz constraint is essential in guaranteeing training stability for GANs. Although heuristic methods such as weight clipping, gradient penalty and spectral normalization have been proposed to enforce Lipschitz constraint, it is still difficult to achieve a solution that is both practically effective and theoretically provably satisfying a Lipschitz constraint. In this paper, we introduce the boundedness and continuity ($BC$) conditions to enforce the Lipschitz constraint on the discriminator functions of GANs. We prove theoretically that GANs with discriminators meeting the BC conditions satisfy the Lipschitz constraint. We present a practically very effective implementation of a GAN based on a convolutional neural network (CNN) by forcing the CNN to satisfy the $BC$ conditions (BC-GAN). We show that as compared to recent techniques including gradient penalty and spectral normalization, BC-GANs not only have better performances but also lower computational complexity.


title: Business Cycles in Economics

abstract: The business cycles are generated by the oscillating macro-/micro-/nano- economic output variables in the economy of the scale and the scope in the amplitude/frequency/phase/time domains in the economics. The accurate forward looking assumptions on the business cycles oscillation dynamics can optimize the financial capital investing and/or borrowing by the economic agents in the capital markets. The book's main objective is to study the business cycles in the economy of the scale and the scope, formulating the Ledenyov unified business cycles theory in the Ledenyov classic and quantum econodynamics.


title: Modeling the effects of telephone nursing on healthcare utilization

abstract: Background: Telephone nursing is the first line of contact for many care-seekers and aims at optimizing the performance of the healthcare system by supporting and guiding patients to the correct level of care and reduce the amount of unscheduled visits. Good statistical models that describe the effects of telephone nursing are important in order to study its impact on healthcare resources and evaluate changes in telephone nursing procedures. Objective: To develop a valid model that captures the complex relationships between the nurse's recommendations, the patients' intended actions and the patients' health seeking behavior. Using the model to estimate the effects of telephone nursing on patient behavior, healthcare utilization, and infer potential cost savings. Methods: Bayesian ordinal regression modeling of data from randomly selected patients that received telephone nursing. Inference is based on Markov Chain Monte Carlo methods, model selection using the Watanabe-Akaike Information Criteria, and model validation using posterior predictive checks on standard discrepancy measures. Results and Conclusions: We present a robust Bayesian ordinal regression model that predicts 76% of the patients' healthcare utilization after telephone nursing and we found no evidence of model deficiencies. The model reveals a risk reducing behavior and the effect of the telephone nursing recommendation is 7 times higher than the effect of the patient's intended action prior to consultation if the recommendation is the highest level of care. But the effect of the nurse's recommendation is lower, or even non-existing, if the recommendation is self-care. Telephone nursing was found to have a constricting effect on healthcare utilization, however, the compliance to nurse's recommendation is closely tied to perceptions of risk, emphasizing the importance to address caller's needs of reassurance.


title: Plastic flow and localization in an amorphous material: experimental interpretation of the fluidity

abstract: We present a thorough study of the plastic response of a granular material progressively loaded. We study experimentally the evolution of the plastic field from a homogeneous one to an heterogeneous one and its fluctuations in term of incremental strain. We show that the plastic field can be decomposed in two components evolving on two decoupled strain increment scales. We argue that the slowly varying part of the field can be identified to the so-called fluidity field introduced recently to interpret the rheological behaviour of amorphous materials. This fluidity field progressively concentrates along a macroscopic direction corresponding to the Mohr-Coulomb angle.


title: Vulnerability of Deep Learning

abstract: The Renormalisation Group (RG) provides a framework in which it is possible to assess whether a deep-learning network is sensitive to small changes in the input data and hence prone to error, or susceptible to adversarial attack. Distinct classification outputs are associated with different RG fixed points and sensitivity to small changes in the input data is due to the presence of relevant operators at a fixed point. A numerical scheme, based on Monte Carlo RG ideas, is proposed for identifying the existence of relevant operators and the corresponding directions of greatest sensitivity in the input data. Thus, a trained deep-learning network may be tested for its robustness and, if it is vulnerable to attack, dangerous perturbations of the input data identified.

categories: stat.ML cond-mat.dis-nn cond-mat.stat-mech cs.AI cs.LG


title: Magnetic shocks and substructures excited by torsional Alfv\'en wave interactions in merging expanding flux tubes

abstract: Vortex motions are frequently observed on the solar photosphere. These motions may play a key role in the transport of energy and momentum from the lower atmosphere into the upper solar atmosphere, contributing to coronal heating. The lower solar atmosphere also consists of complex networks of flux tubes that expand and merge throughout the chromosphere and upper atmosphere. We perform numerical simulations to investigate the behaviour of vortex driven waves propagating in a pair of such flux tubes in a non-force-free equilibrium with a realistically modelled solar atmosphere. The two flux tubes are independently perturbed at their footpoints by counter-rotating vortex motions. When the flux tubes merge, the vortex motions interact both linearly and nonlinearly. The linear interactions generate many small-scale transient magnetic substructures due to the magnetic stress imposed by the vortex motions. Thus, an initially monolithic tube is separated into a complex multi-threaded tube due to the photospheric vortex motions. The wave interactions also drive a superposition that increases in amplitude until it exceeds the local Mach number and produces shocks that propagate upwards with speeds of approximately $ 50$ km s$^{-1}$. The shocks act as conduits transporting momentum and energy upwards, and heating the local plasma by more than an order of magnitude, with peak temperature approximately $60,000$ K. Therefore, we present a new mechanism for the generation of magnetic waveguides from the lower solar atmosphere to the solar corona. This wave guide appears as the result of interacting perturbations in neighbouring flux tubes. Thus, the interactions of photospheric vortex motions is a potentially significant mechanism for energy transfer from the lower to upper solar atmosphere.


title: Multiple structure and symmetry types in narrow temperature and magnetic field ranges in two-dimensional Cr2Ge2Te6 crystal

abstract: Multiple structure and symmetry types and their transformations are discovered in quasi-two-dimensional (quasi-2D) Cr2Ge2Te6 crystal in surprisingly very narrow temperature range of 2 K and magnetic field range of 0.07 T using a homebuilt magnetic force microscope (MFM). A series of basic domain patterns are extracted from the MFM images. Some of them seem unique to 2-D materials as they are not observed in 3-D materials, such as self-fitting disks, and fine ladder structure within Y-connected walls. Based on these findings, a phase map is drawn for the magnetic phase structures. The symmetry of the patterns are discussed. The results are not only important in developing new theories but also highly desirable in applications.


title: A constant-ratio approximation algorithm for a class of hub-and-spoke network design problems and metric labeling problems: star metric case

abstract: Transportation networks frequently employ hub-and-spoke network architectures to route flows between many origin and destination pairs. Hub facilities work as switching points for flows in large networks. In this study, we deal with a problem, called the single allocation hub-and-spoke network design problem. In the problem, the goal is to allocate each non-hub node to exactly one of given hub nodes so as to minimize the total transportation cost. The problem is essentially equivalent to another combinatorial optimization problem, called the metric labeling problem. The metric labeling problem was first introduced by Kleinberg and Tardos in 2002, motivated by application to segmentation problems in computer vision and related areas. In this study, we deal with the case where the set of hubs forms a star, which arises especially in telecommunication networks. We propose a polynomial-time randomized approximation algorithm for the problem, whose approximation ratio is less than 5.281. Our algorithms solve a linear relaxation problem and apply dependent rounding procedures.


title: Multi-gradient fluids

abstract: An internal energy function of the mass density, the volumetric entropy and their gradients at n-order generates the representation of multi-gradient fluids. Thanks to Hamilton's principle, we obtain a thermodynamical form of the equation of motion which generalizes the case of perfect compressible fluids. First integrals of flows are extended cases of perfect compressible fluids. The equation of motion and the equation of energy are written for dissipative cases, and are compatible with the second law of thermodynamics.