title: Neural Architecture Search For Fault Diagnosis

abstract: Data-driven methods have made great progress in fault diagnosis, especially deep learning method. Deep learning is suitable for processing big data, and has a strong feature extraction ability to realize end-to-end fault diagnosis systems. However, designing neural network architecture requires rich professional knowledge and debugging experience, and a lot of experiments are needed to screen models and hyperparameters, increasing the difficulty of developing deep learning models. Frortunately, neural architecture search (NAS) is developing rapidly, and is becoming one of the next directions for deep learning. In this paper, we proposed a NAS method for fault diagnosis using reinforcement learning. A recurrent neural network is used as an agent to generate network architecture. The accuracy of the generated network on the validation dataset is fed back to the agent as a reward, and the parameters of the agent are updated through the strategy gradient algorithm. We use PHM 2009 Data Challenge gearbox dataset to prove the effectiveness of proposed method, and obtain state-of-the-art results compared with other artificial designed network structures. To author's best knowledge, it's the first time that NAS has been applied in fault diagnosis.


title: Generalized list colouring of graphs

abstract: This paper disproves a conjecture of Wang, Wu, Yan and Xie, and answers in negative a question in Dvorak, Pekarek and Sereni. In return, we pose five open problems.


title: Fluctuation Theorem as a special case of Girsanov Theorem

abstract: Stochastic thermodynamics is an important development in the direction of finding general thermodynamic principles for non-equilibrium systems. We believe stochastic thermodynamics has the potential to benefit from the measure-theoretic framework of stochastic differential equations. Towards this, in this work, we show that Fluctuation Theorem (FT) is a special case of the Girsanov theorem, which is an important result in the theory of stochastic differential equations. We report that by employing Girsanov transformation of measures between the forward and the reversed dynamics of a general class of Langevin dynamic systems, we arrive at the Integral Fluctuation Relation. Following the same approach, we derive the FT also for the overdamped case. Our derivation is applicable to both transient and steady state conditions and can also incorporate diffusion coefficients varying as a function of state and time. We expect that the proposed method will be an easy route towards deriving the FT irrespective of the complexity and non-linearity of the system.


title: Action-Manipulation Attacks Against Stochastic Bandits: Attacks and Defense

abstract: Due to the broad range of applications of stochastic multi-armed bandit model, understanding the effects of adversarial attacks and designing bandit algorithms robust to attacks are essential for the safe applications of this model. In this paper, we introduce a new class of attack named action-manipulation attack. In this attack, an adversary can change the action signal selected by the user. We show that without knowledge of mean rewards of arms, our proposed attack can manipulate Upper Confidence Bound (UCB) algorithm, a widely used bandit algorithm, into pulling a target arm very frequently by spending only logarithmic cost. To defend against this class of attacks, we introduce a novel algorithm that is robust to action-manipulation attacks when an upper bound for the total attack cost is given. We prove that our algorithm has a pseudo-regret upper bounded by $\mathcal{O}(\max\{\log T,A\})$, where $T$ is the total number of rounds and $A$ is the upper bound of the total attack cost.


title: Physical Parameters of Late-type Contact Binaries in the Northern Catalina Sky Survey

abstract: We present the physical parameters of 2335 late-type contact binary (CB) systems extracted from the Catalina Sky Survey (CSS). Our sample was selected from the CSS Data Release 1 by strictly limiting the prevailing temperature uncertainties and light-curve fitting residuals, allowing us to almost eliminate any possible contaminants. We developed an automatic Wilson--Devinney-type code to derive the relative properties of CBs based on their light-curve morphology. By adopting the distances derived from CB (orbital) period--luminosity relations (PLRs), combined with the well-defined mass--luminosity relation for the systems' primary stars and assuming solar metallicity, we calculated the objects' masses, radii, and luminosities. Our sample of fully eclipsing CBs contains 1530 W-, 710 A-, and 95 B-type CBs. A comparison with literature data and with the results from different surveys confirms the accuracy and coherence of our measurements. The period distributions of the various CB subtypes are different, hinting at a possible evolutionary sequence. W-type CBs are clearly located in a strip in the total mass versus mass ratio plane, while A-type CBs may exhibit a slightly different dependence. There are no significant differences among the PLRs of A- and W-type CBs, but the PLR zero points are affected by their mass ratios and fill-out factors. Determination of zero-point differences for different types of CBs may help us improve the accuracy of the resulting PLRs. We demonstrate that automated approaches to deriving CB properties could be a powerful tool for application to the much larger CB samples expected to result from future surveys.


title: Finiteness in Polygonal Billiards on Hyperbolic Plane

abstract: \textsc{J. Hadamard} studied the geometric properties of geodesic flows on surfaces of negative curvature, thus initiating "Symbolic Dynamics". In this article, we follow the same geometric approach to study the geodesic trajectories of billiards in "rational polygons" on the hyperbolic plane. We particularly show that the billiard dynamics resulting thus are just 'Subshifts of Finite Type' or their dense subsets. We further show that 'Subshifts of Finite Type' play a central role in subshift dynamics and while discussing the topological structure of the space of all subshifts, we demonstrate that they approximate any shift dynamics.


title: A Two-Stage Decomposition Approach for AC Optimal Power Flow

abstract: The alternating current optimal power flow (AC-OPF) problem is critical to power system operations and planning, but it is generally hard to solve due to its nonconvex and large-scale nature. This paper proposes a scalable decomposition approach in which the power network is decomposed into a master network and a number of subnetworks, where each network has its own AC-OPF subproblem. This formulates a two-stage optimization problem and requires only a small amount of communication between the master and subnetworks. The key contribution is a smoothing technique that renders the response of a subnetwork differentiable with respect to the input from the master problem, utilizing properties of the barrier problem formulation that naturally arises when subproblems are solved by a primal-dual interior-point algorithm. Consequently, existing efficient nonlinear programming solvers can be used for both the master problem and the subproblems. The advantage of this framework is that speedup can be obtained by processing the subnetworks in parallel, and it has convergence guarantees under reasonable assumptions. The formulation is readily extended to instances with stochastic subnetwork loads. Numerical results show favorable performance and illustrate the scalability of the algorithm which is able to solve instances with more than 11 million buses.


title: Fast and linear-time string matching algorithms based on the distances of $q$-gram occurrences

abstract: Given a text $T$ of length $n$ and a pattern $P$ of length $m$, the string matching problem is a task to find all occurrences of $P$ in $T$. In this study, we propose an algorithm that solves this problem in $O((n + m)q)$ time considering the distance between two adjacent occurrences of the same $q$-gram contained in $P$. We also propose a theoretical improvement of it which runs in $O(n + m)$ time, though it is not necessarily faster in practice. We compare the execution times of our and existing algorithms on various kinds of real and artificial datasets such as an English text, a genome sequence and a Fibonacci string. The experimental results show that our algorithm is as fast as the state-of-the-art algorithms in many cases, particularly when a pattern frequently appears in a text.


title: On-line non-overlapping camera calibration net

abstract: We propose an easy-to-use non-overlapping camera calibration method. First, successive images are fed to a PoseNet-based network to obtain ego-motion of cameras between frames. Next, the pose between cameras are estimated. Instead of using a batch method, we propose an on-line method of the inter-camera pose estimation. Furthermore, we implement the entire procedure on a computation graph. Experiments with simulations and the KITTI dataset show the proposed method to be effective in simulation.


title: Tailoring the Spectral Absorption Coefficient of a Blended Plasmonic Nanofluid Using a Customized Genetic Algorithm

abstract: Recently, plasmonic nanofluids (i.e., a suspension of plasmonic nanoparticles in a base fluid) have been widely employed in direct-absorption solar collectors because the localized surface plasmon supported by plasmonic nanoparticles can greatly improve the direct solar thermal conversion performance. Considering that the surface plasmon resonance frequency of metallic nanoparticles, such as gold, silver, and aluminum, is usually located in the ultraviolet to visible range, the absorption coefficient of a plasmonic nanofluid must be spectrally tuned for full utilization of the solar radiation in a broad spectrum. In the present study, a modern design process in the form of a genetic algorithm (GA) is applied to the tailoring of the spectral absorption coefficient of a plasmonic nanofluid. To do this, the major components of a conventional GA, such as the gene description, fitness function for the evaluation, crossover, and mutation function, are modified to be suitable for the inverse problem of tailoring the spectral absorption coefficient of a plasmonic nanofluid. By applying the customized GA, we obtained an optimal combination for a blended nanofluid with the desired spectral distribution of the absorption coefficient, specifically a uniform distribution, solar-spectrum-like distribution, and a step-function-like distribution. The resulting absorption coefficient of the designed plasmonic nanofluid is in good agreement with the prescribed spectral distribution within about 10\% to 20\% of error when six types of nanoparticles are blended. Finally, we also investigate how the inhomogeneous broadening effect caused by the fabrication uncertainty of the nanoparticles changes their optimal combination.


title: Detection and Analysis of Drive-by Downloads and Malicious Websites

abstract: A drive by download is a download that occurs without users action or knowledge. It usually triggers an exploit of vulnerability in a browser to downloads an unknown file. The malicious program in the downloaded file installs itself on the victims machine. Moreover, the downloaded file can be camouflaged as an installer that would further install malicious software. Drive by downloads is a very good example of the exponential increase in malicious activity over the Internet and how it affects the daily use of the web. In this paper, we try to address the problem caused by drive by downloads from different standpoints. We provide in depth understanding of the difficulties in dealing with drive by downloads and suggest appropriate solutions. We propose machine learning and feature selection solutions to remedy the the drive-by download problem. Experimental results reported 98.2% precision, 98.2% F-Measure and 97.2% ROC area.


title: Complete Finsler spaces of constant negative Ricci curvature

abstract: Here, using the projectively invariant pseudo-distance and Schwarzian derivative, it is shown that every connected complete Finsler space of the constant negative Ricci scalar is reversible. In particular, every complete Randers metric of constant negative Ricci (or flag) curvature is Riemannian.


title: The Benefits of Probability-Proportional-to-Size Sampling in Cluster-Randomized Experiments

abstract: In a cluster-randomized experiment, treatment is assigned to clusters of individual units of interest--households, classrooms, villages, etc.--instead of the units themselves. The number of clusters sampled and the number of units sampled within each cluster is typically restricted by a budget constraint. Previous analysis of cluster randomized experiments under the Neyman-Rubin potential outcomes model of response have assumed a simple random sample of clusters. Estimators of the population average treatment effect (PATE) under this assumption are often either biased or not invariant to location shifts of potential outcomes. We demonstrate that, by sampling clusters with probability proportional to the number of units within a cluster, the Horvitz-Thompson estimator (HT) is invariant to location shifts and unbiasedly estimates PATE. We derive standard errors of HT and discuss how to estimate these standard errors. We also show that results hold for stratified random samples when samples are drawn proportionally to cluster size within each stratum. We demonstrate the efficacy of this sampling scheme using a simulation based on data from an experiment measuring the efficacy of the National Solidarity Programme in Afghanistan.


title: Boundary Curvature Effect on the Wrinkling of Thin Suspended Films

abstract: In this letter, we demonstrate a relation between the boundary curvature $\kappa$ and the wrinkle wavelength $\lambda$ of a thin suspended film under boundary confinement. Experiments are done with nanocrystalline diamond films of thickness $t \approx 184$~nm grown on glass substrates. By removing portions of the substrate after growth, suspended films with circular boundaries of radius $R$ ranging from approximately 30 to 811 $\mu$m are made. Due to residual stresses, the portions of film attached to the substrate are of compressive prestrain $\epsilon_0 \approx 11 \times 10^{-4}$ and the suspended portions of film are azimuthally wrinkled at their boundary. We find that $\lambda$ monotonically decreases with $\kappa$ and present a model predicting that $\lambda \propto t^{1/2}(\epsilon_0 + \Delta R \kappa)^{-1/4}$, where $\Delta R$ denotes a penetration depth over which strain relaxes at a boundary. This relation is in agreement with our experiments and may be adapted to other systems such as plant leaves. Also, we establish a novel method for measuring residual compressive strain in thin films.


title: Masses of double neutron star mergers

abstract: Aims. The mass discrepancy between the observed population of double neutron star binaries by radio pulsar observations and gravitational-wave observation requires an explanation. Methods. Binary population synthesis calculations are performed, and their results are compared with the radio and the gravitational-wave observations simultaneously. Results. Simulations of binary evolution are used to link different observations of double neutron star binaries with each other. The progenitor of GW190425 is investigated in more detail. A distribution of masses and merger times of the possible progenitors is presented. Conclusions. A mass discrepancy between the radio pulsars in the Milky Way with another neutron star companion and the inferred masses from gravitational-wave observations of those kind of merging systems is naturally found in binary evolution.


title: Indirect Adversarial Attacks via Poisoning Neighbors for Graph Convolutional Networks

abstract: Graph convolutional neural networks, which learn aggregations over neighbor nodes, have achieved great performance in node classification tasks. However, recent studies reported that such graph convolutional node classifier can be deceived by adversarial perturbations on graphs. Abusing graph convolutions, a node's classification result can be influenced by poisoning its neighbors. Given an attributed graph and a node classifier, how can we evaluate robustness against such indirect adversarial attacks? Can we generate strong adversarial perturbations which are effective on not only one-hop neighbors, but more far from the target? In this paper, we demonstrate that the node classifier can be deceived with high-confidence by poisoning just a single node even two-hops or more far from the target. Towards achieving the attack, we propose a new approach which searches smaller perturbations on just a single node far from the target. In our experiments, our proposed method shows 99% attack success rate within two-hops from the target in two datasets. We also demonstrate that m-layer graph convolutional neural networks have chance to be deceived by our indirect attack within m-hop neighbors. The proposed attack can be used as a benchmark in future defense attempts to develop graph convolutional neural networks with having adversary robustness.


title: CNN-based approach for glaucoma diagnosis using transfer learning and LBP-based data augmentation

abstract: Glaucoma causes an irreversible damage to retinal nerve fibers which results in vision loss, if undetected in early stage. Therefore, diagnosis of glaucoma in its early stage may prevent further vision loss. In this paper, we propose a convolutional neural network (CNN) based approach for automated glaucoma diagnosis by employing retinal fundus images. This approach employs transfer learning technique and local binary pattern (LBP) based data augmentation. In the proposed approach, we employ Alexnet as a pre-trained CNN model which is used for transfer learning. Initially, the proposed approach divides the fundus image dataset into training and testing data. Further, the color fundus images in training and testing data are separated into red (R), green (G), and blue (B) channels. Additionally, the LBP-based data augmentation is performed on training data. Specifically, we compute LPBs for each of the channel. Finally, the augmented training data is used to train the CNN model via transfer learning. In testing stage, the R, G, and B channels of test image are fed to the trained CNN model which generates 3 decisions. We employ a decision level fusion technique to combine the decisions obtained from the trained CNN model. The experimental evaluation of the proposed approach on the public RIM-ONE fundus image database, achieves state-of-the-art performance for glaucoma diagnosis.


title: Communication-Efficient Distributed SVD via Local Power Iterations

abstract: We study distributed computing of the truncated singular value decomposition problem. We develop an algorithm that we call \texttt{LocalPower} for improving communication efficiency. Specifically, we uniformly partition the dataset among $m$ nodes and alternate between multiple (precisely $p$) local power iterations and one global aggregation. In the aggregation, we propose to weight each local eigenvector matrix with orthogonal Procrustes transformation (OPT). As a practical surrogate of OPT, sign-fixing, which uses a diagonal matrix with $\pm 1$ entries as weights, has better computation complexity and stability in experiments. We theoretically show that under certain assumptions \texttt{LocalPower} lowers the required number of communications by a factor of $p$ to reach a constant accuracy. We also show that the strategy of periodically decaying $p$ helps obtain high-precision solutions. We conduct experiments to demonstrate the effectiveness of \texttt{LocalPower}.


title: The Schwarzian derivative and conformal transformation on Finsler manifolds

abstract: Thurston, in 1986, discovered that the Schwarzian derivative has mysterious properties similar to the curvature on a manifold. After his work, there are several approaches to develop this notion on Riemannian manifolds. Here, a tensor field is identified in the study of global conformal diffeomorphisms on Finsler manifolds as a natural generalization of the Schwarzian derivative. Then, a natural definition of a Mobius mapping on Finsler manifolds is given and its properties are studied. In particular, it is shown that Mobius mappings are mappings that preserve circles and vice versa. Therefore, if a forward geodesically complete Finsler manifold admits a Mobius mapping, then the indicatrix is conformally diffeomorphic to the Euclidean sphere $ S^{n-1}$ in $ \mathbb{R}^n $. In addition, if a forward geodesically complete absolutely homogeneous Finsler manifold of scalar flag curvature admits a non-trivial change of Mobius mapping, then it is a Riemannian manifold of constant sectional curvature.


title: Ion Exchange in Silicate Glasses: Physics of Ion Concentration, Residual Stress, and Refractive Index Profiles

abstract: A systematic review of main physical effects generated by ion exchange in silicate glasses is presented. Ion concentration distributions, residual stress profiles, and refractive index effects are discussed with particular attention on the physical and mathematical underpinnings of the ion exchange process. The study has the purpose of presenting a scientific foundation to enable future developments in the field. In this respect, the objective of this article is more educational than to present new research results. Appendixes are included to consider the detail of some specific topics without disrupting the continuity of the overall discussion. Despite the review approach of this study, some original topics are included, such as the concentration distribution with variable boundary conditions and residual stress profile with anomalies due to different relaxation mechanisms, viz., either isochoric and non-isochoric and stress driven or free energy driven. The time scale of the different relaxation mechanisms results, for some specific glass chemical compositions, in the appearance of a subsurface compression maximum progressively moving apart from glass surface and eventually turning in a reversal from compression into tensile state upon prolonged ion exchange processes. Finally, a broad discussion on optical effects induced by ion exchange is presented, paying particular attention to the possibility of experimental determination of residual stress profile.


title: Gaseous dynamical friction under radiative feedback: do intermediate-mass black holes speed up or down?

abstract: Coalescence of intermediate-mass black holes (IMBHs) as a result of the migration toward galactic centers via dynamical friction may contribute to the formation of supermassive BHs. Here we reinvestigate the gaseous dynamical friction, which was claimed to be inefficient with radiative feedback from BHs in literature, by performing 3D radiation-hydrodynamics simulations that solve the flow structure in the vicinity of BHs. We consider a $10^4~M_\odot$ BH moving at the velocity $V_{\rm flow}$ through the homogeneous medium with metallicity $Z$ in the range of $0-0.1~Z_\odot$ and density $n_{\infty}$. We show that, if $n_{\infty} \lesssim 10^{6}~{\rm cm^{-3}}$ and $V_{\rm flow} \lesssim 60~{\rm km~s^{-1}}$, the BH is accelerated forward because of the gravitational pull from a dense shell ahead of an ionized bubble around the BH, regardless of the value of $Z$. If $n_{\infty} \gtrsim 10^{6}~{\rm cm^{-3}}$, however, our simulation shows the opposite result. The ionized bubble and associating shell temporarily appear, but immediately go downstream with significant ram pressure of the flow. They eventually converge into a massive downstream wake, which gravitationally drags the BH backward. The BH decelerates over the timescale of $\sim 0.01$~Myr, much shorter than the dynamical timescale in galactic disks. Our results suggest that IMBHs that encounter the dense clouds rapidly migrate toward galactic centers, where they possibly coalescence with others.


title: Task-space Synergies for Reaching using Upper-limb Prostheses

abstract: Synergistic prostheses enable the coordinated movement of the human-prosthetic arm, as required by activities of daily living. This is achieved by coupling the motion of the prosthesis to the human command, such as the residual limb movement in motion-based interfaces. Previous studies demonstrated that developing human-prosthetic synergies in joint-space must consider individual motor behaviour and the intended task to be performed, requiring personalisation and task calibration. In this work, an alternative synergy-based strategy, utilising a synergistic relationship expressed in task-space, is proposed. This task-space synergy has the potential to replace the need for personalisation and task calibration with a model-based approach requiring knowledge of the individual user's arm kinematics, the anticipated hand motion during the task and voluntary information from the prosthetic user. The proposed method is compared with surface electromyography-based and joint-space synergy-based prosthetic interfaces in a study of motor behaviour and task performance on able-bodied subjects using a VR-based transhumeral prosthesis. Experimental results showed that for a set of forward reaching tasks the proposed task-space synergy achieves comparable performance to joint-space synergies without the need to rely on time-consuming calibration processes or human motor learning. Case study results with an amputee subject motivate the further development of the proposed task-space synergy method.


title: A proposal to fight tornadoes with multiple connected balloons

abstract: A tornado is an extreme weather condition that can cause enormous damage to human society. In this paper, we propose a low cost, environmentally friendly method to fight against tornadoes using clusters of connected balloons. Deployed over a sufficiently wide area, these large balloons may be able to reduce the wind speed of the tornado, block and disrupt the convective flow of air, and destroy the tornado.


title: Probing blue long persistent Eu2+/Nd3+ activated (Ca$_{0.9}$Si$_{0.1}$)Al$_2$O$_6$ emerging phosphor based novel invisible ink for knocking out fake security papers and combat counterfeiting

abstract: In the current work, a novel class of (Ca0.9Si0.1)Al2O6:Eu2+,Nd3+ long persistent phosphor was synthesized using high-temperature (1000-1500oC) solid-state reaction method. Further, the phosphor material was formulated into a phosphorescent ink using epoxy based colloidal suspension. Interestingly the phosphorescent ink absorbs ambient energies such as white light from D65 and tungsten lamps, sunlight, artificial lights etc., stores up the energy and emit as efficient blue colored photoluminescence (PL) apparently visible for many hours for dark adapted human eye. However, the optimum PL emission centred at ~450 nm has been registered at ~365 nm excitation that is attributed to transitions from 4f65d1 to 4f7 energy levels of Eu2+ ions. The longer persistence of PL aroused due to trapping and de-trapping of holes at Nd3+ sites situated proximate to valance band. In the optimised stoichiometric composition of the phosphor, the persistence lasted for >4 hours in extremely dark conditions. Furthermore, the phenomenon of long persistence enables to create/design unique hidden markings on goods and certificates using quick response (QR) code patterns. Systematic studies have been performed on various templates kept under UV (~365 nm) excitation to identify fake currencies, barcodes and combat counterfeiting effectively.


title: Seasonal and Trend Forecasting of Tourist Arrivals: An Adaptive Multiscale Ensemble Learning Approach

abstract: The accurate seasonal and trend forecasting of tourist arrivals is a very challenging task. In the view of the importance of seasonal and trend forecasting of tourist arrivals, and limited research work paid attention to these previously. In this study, a new adaptive multiscale ensemble (AME) learning approach incorporating variational mode decomposition (VMD) and least square support vector regression (LSSVR) is developed for short-, medium-, and long-term seasonal and trend forecasting of tourist arrivals. In the formulation of our developed AME learning approach, the original tourist arrivals series are first decomposed into the trend, seasonal and remainders volatility components. Then, the ARIMA is used to forecast the trend component, the SARIMA is used to forecast seasonal component with a 12-month cycle, while the LSSVR is used to forecast remainder volatility components. Finally, the forecasting results of the three components are aggregated to generate an ensemble forecasting of tourist arrivals by the LSSVR based nonlinear ensemble approach. Furthermore, a direct strategy is used to implement multi-step-ahead forecasting. Taking two accuracy measures and the Diebold-Mariano test, the empirical results demonstrate that our proposed AME learning approach can achieve higher level and directional forecasting accuracy compared with other benchmarks used in this study, indicating that our proposed approach is a promising model for forecasting tourist arrivals with high seasonality and volatility.


title: Properties of Prime Products and Robin's Inequality

abstract: There are many formulations of problems that have been proven to be equivalent to the Riemann Hypothesis in modern mathematics. In this paper we look at the formulation of an inequality derived by Robin in 1984 that proves the Riemann Hypothesis true if and only if Robin's Inequality is true for all $n\geq 5040$. In this paper we look at explicit products of prime numbers and show that if a given $n=p_1^{k_1}p_2^{k_2}\ldots p_m^{k_m}$ satisfies Robin's Inequality, then $n=P_1^{k_1}P_2^{k_2}\ldots P_m^{k_m}$ with $p_j \leq P_j$ also satisfies the inequality. We also then offer two conjectures that if proven individually could imply some interesting properties of the Riemann Hypothesis and Robin's Inequality.


title: Dust Destruction by Charging: A Possible Origin of Grey Extinction Curves of Active Galactic Nuclei

abstract: Observed extinction curves of active galactic nuclei (AGNs) are significantly different from those observed in the Milky Way. The observations require preferential removal of small grains at the AGN environment; however, the physics for this remains unclear. In this paper, we propose that dust destruction by charging, or Coulomb explosion, may be responsible for AGN extinction curves. Harsh AGN radiation makes a dust grain highly charged through photoelectric emission, and grain fission via the Coulomb explosion occurs when the electrostatic tensile stress of a charge grain exceeds its tensile strength. We show that the Coulomb explosion can preferentially remove both small silicate and graphite grains and successfully reproduce both flat extinction curves and the absence of 2175\AA~bump.


abstract: Recent efforts in Dialogue State Tracking (DST) for task-oriented dialogues have progressed toward open-vocabulary or generation-based approaches where the models can generate slot value candidates from the dialogue history itself. These approaches have shown good performance gain, especially in complicated dialogue domains with dynamic slot values. However, they fall short in two aspects: (1) they do not allow models to explicitly learn signals across domains and slots to detect potential dependencies among (domain, slot) pairs; and (2) existing models follow auto-regressive approaches which incur high time cost when the dialogue evolves over multiple domains and multiple turns. In this paper, we propose a novel framework of Non-Autoregressive Dialog State Tracking (NADST) which can factor in potential dependencies among domains and slots to optimize the models towards better prediction of dialogue states as a complete set rather than separate slots. In particular, the non-autoregressive nature of our method not only enables decoding in parallel to significantly reduce the latency of DST for real-time dialogue response generation, but also detect dependencies among slots at token level in addition to slot and domain level. Our empirical results show that our model achieves the state-of-the-art joint accuracy across all domains on the MultiWOZ 2.1 corpus, and the latency of our model is an order of magnitude lower than the previous state of the art as the dialogue history extends over time.


title: Influence Function based Data Poisoning Attacks to Top-N Recommender Systems

abstract: Recommender system is an essential component of web services to engage users. Popular recommender systems model user preferences and item properties using a large amount of crowdsourced user-item interaction data, e.g., rating scores; then top-$N$ items that match the best with a user's preference are recommended to the user. In this work, we show that an attacker can launch a data poisoning attack to a recommender system to make recommendations as the attacker desires via injecting fake users with carefully crafted user-item interaction data. Specifically, an attacker can trick a recommender system to recommend a target item to as many normal users as possible. We focus on matrix factorization based recommender systems because they have been widely deployed in industry. Given the number of fake users the attacker can inject, we formulate the crafting of rating scores for the fake users as an optimization problem. However, this optimization problem is challenging to solve as it is a non-convex integer programming problem. To address the challenge, we develop several techniques to approximately solve the optimization problem. For instance, we leverage influence function to select a subset of normal users who are influential to the recommendations and solve our formulated optimization problem based on these influential users. Our results show that our attacks are effective and outperform existing methods.


title: On homeomorphism groups and the set-open topology

abstract: In this paper we focus on the set-open topologies on the group $\mathcal{H}(X)$ of all self-homeomorphisms of a topological space $X$ which yield continuity of both the group operations, product and inverse function. As a consequence, we make the more general case of Dijkstra's theorem. In this case a homogeneously encircling family $\mathcal{B}$ consists of regular open sets and the closure of every set from $\mathcal{B}$ is contained in the finite union of connected sets from $\mathcal{B}$. Also we proved that the zero-cozero topology of $\mathcal{H}(X)$ is the relativisation to $\mathcal{H}(X)$ of the compact-open topology of $\mathcal{H}(\beta X)$ for any Tychonoff space $X$ and every homogeneous zero-dimensional space $X$ can be represented as the quotient space of a topological group with respect to a closed subgroup.


title: Toward Low-Cost and Stable Blockchain Networks

abstract: Envisioned to be the future of secured distributed systems, blockchain networks have received increasing attention from both the industry and academia in recent years. However, blockchain mining processes demand high hardware costs and consume a vast amount of energy (studies have shown that the amount of energy consumed in Bitcoin mining is almost the same as the electricity used in Ireland). To address the high mining cost problem of blockchain networks, in this paper, we propose a blockchain mining resources allocation algorithm to reduce the mining cost in PoW-based (proof-of-work-based) blockchain networks. We first propose an analytical queueing model for general blockchain networks. In our queueing model, transactions arrive randomly to the queue and are served in a batch manner with unknown service rate probability distribution and agnostic to any priority mechanism. Then, we leverage the Lyapunov optimization techniques to propose a dynamic mining resources allocation algorithm (DMRA), which is parameterized by a tuning parameter $K>0$. We show that our algorithm achieves an $[O(1/K), O(K)]$ cost-optimality-gap-vs-delay tradeoff. Our simulation results also demonstrate the effectiveness of DMRA in reducing mining costs.


title: Variability and the size-luminosity relation of the intermediate mass AGN in NGC 4395

abstract: We present the variability study of the lowest-luminosity Seyfert 1 galaxy NGC 4395 based on the photometric monitoring campaigns in 2017 and 2018. Using 22 ground-based and space telescopes, we monitored NGC 4395 with a $\sim$5 minute cadence during a period of 10 days and obtained light curves in the UV, V, J, H, and K/Ks bands as well as the H$\alpha$ narrow-band. The RMS variability is $\sim$0.13 mag on \emph{Swift}-UVM2 and V filter light curves, decreasing down to $\sim$0.01 mag on K filter. After correcting for continuum contribution to the H$\alpha$ narrow-band, we measured the time lag of the H$\alpha$ emission line with respect to the V-band continuum as ${55}^{+27}_{-31}$ to ${122}^{+33}_{-67}$ min. in 2017 and ${49}^{+15}_{-14}$ to ${83}^{+13}_{-14}$ min. in 2018, depending on the assumption on the continuum variability amplitude in the H$\alpha$ narrow-band. We obtained no reliable measurements for the continuum-to-continuum lag between UV and V bands and among near-IR bands, due to the large flux uncertainty of UV observations and the limited time baseline. We determined the AGN monochromatic luminosity at 5100\AA\ $\lambda L_\lambda = \left(5.75\pm0.40\right)\times 10^{39}\,\mathrm{erg\,s^{-1}}$, after subtracting the contribution of the nuclear star cluster. While the optical luminosity of NGC 4395 is two orders of magnitude lower than that of other reverberation-mapped AGNs, NGC 4395 follows the size-luminosity relation, albeit with an offset of 0.48 dex ($\geq$2.5$\sigma$) from the previous best-fit relation of Bentz et al. (2013).


title: Probabilistic principal component and linear discriminant analysis of interaction of 1064 nm CW laser with ICP plasma: Anti-Stokes cooling due to analogue black hole

abstract: Pattern recognition and machine learning techniques are known to play an emerging role for studying underlying physics of light-matter interaction. While on the subject, photon condensation and Anti-Stokes cooling remain compelling phenomena of laser-underdense plasma interaction. In this work, the interaction of a 1064 nm continuum-wave laser with inductively-coupled plasma (ICP) of Mercury(Hg) has been studied by probabilistic pattern recognition over observed time resolved spectral data.3D vector fields obtained by the probabilistic pattern recognition over spectral database proficiently illuminates the analogue black hole (ABH) and its sonic horizon in which homoclinic orbit Whistler waves are exited. Vector spectra show the collective behavior of condensation and anti-Stokes cooling is due to the electronic transitions of Hg1 trapped by the phonon sink of ABH of laser. Modeling of subsonic phonon vector spectrum obtained by the probabilistic linear discriminant analysis(PLDA) estimates the 0.70 nanoKelvin of temperature of the sink region of ABH.


title: An Efficient Transfer Learning Framework for Multiagent Reinforcement Learning

abstract: Transfer Learning has shown great potential to enhance single-agent Reinforcement Learning (RL) efficiency. Similarly, Multiagent RL (MARL) can also be accelerated if agents can share knowledge with each other. However, it remains a problem of how an agent should learn from other agents. In this paper, we propose a novel Multiagent Policy Transfer Framework (MAPTF) to improve MARL efficiency. MAPTF learns which agent's policy is the best to reuse for each agent and when to terminate it by modeling multiagent policy transfer as the option learning problem. Furthermore, in practice, the option module can only collect all agent's local experiences for update due to the partial observability of the environment. While in this setting, each agent's experience may be inconsistent with each other, which may cause the inaccuracy and oscillation of the option-value's estimation. Therefore, we propose a novel option learning algorithm, the successor representation option learning to solve it by decoupling the environment dynamics from rewards and learning the option-value under each agent's preference. MAPTF can be easily combined with existing deep RL and MARL approaches, and experimental results show it significantly boosts the performance of existing methods in both discrete and continuous state spaces.


title: Discovery of a pre-cataclysmic binary with unusual chromaticity of the eclipsed white dwarf by the GPX Survey

abstract: We report the discovery of a relatively bright eclipsing binary system, which consists of a white dwarf and a main sequence K7 star with clear signs of chromospheric and spot activity. The light curve of this system shows $\sim0.2$mag ellipsoidal variability with a period of 0.297549d and a short total eclipse of the white dwarf. Based on our analysis of the spectral and photometric data, we estimated the parameters of the system. The K7V star is tidally deformed but does not fill its Roche lobe (the filling factor is about 0.86). The orbital inclination is $i=73^\circ.1\pm 0^\circ.2$, the mass ratio is $q=M_2/M_1\approx 0.88$. The parameters of the K7V star are $M_2\approx 0.64$M$_{\odot}$, $R_2=0.645\pm 0.012$R$_{\odot}$, $T_2\approx 4070$K. The parameters of the white dwarf are $M_1\approx 0.72$M$_{\odot}$, $R_1=0.013\pm 0.003$R$_{\odot}$, $T_1=8700\pm 1100$K. Photometric observations in different bands revealed that the maximum depth of the eclipse is in the \textit{SDSS r} filter, which is unusual for a system of a white dwarf and a late main sequence star. We suspect that this system is a product of the evolution of a common envelope binary star, and that the white dwarf accretes the stellar wind from the secondary star (the so-called low-accretion rate polar, LARP).


title: A Fixed point view: A Model-Based Clustering Framework

abstract: With the inflation of the data, clustering analysis, as a branch of unsupervised learning, lacks unified understanding and application of its mathematical law. Based on the view of fixed point, this paper restates the model-based clustering and proposes a unified clustering framework. In order to find fixed points as cluster centers, the framework iteratively constructs the contraction map, which strongly reveals the convergence mechanism and interconnections among algorithms. By specifying a contraction map, Gaussian mixture model (GMM) can be mapped to the framework as an application. We hope the fixed point framework will help the design of future clustering algorithms.


title: Theoretical study of the spectroscopic properties of mendelevium ($Z=101$)

abstract: Using recently developed version of the configuration interaction method for atoms with open shells we calculate electron structure and spectroscopic properties of the mendelevium atom (Md, $Z=101$). These include energy levels, first and second ionisation potentials, electron affinity, hyperfine structure and electric dipole transition amplitudes between ground and low lying states of opposite parity. The accuracy of the calculations is controlled by performing similar calculations for lighter analog of mendelevium, thulium atom and comparing the results with experiment and other calculations. The calculations for Md are to address the lack of experimental data and help in planing and interpreting the measurements.


title: Emotion Recognition Through Observer's Physiological Signals

abstract: Emotion recognition based on physiological signals is a hot topic and has a wide range of applications, like safe driving, health care and creating a secure society. This paper introduces a physiological dataset PAFEW, which is obtained using movie clips from the Acted Facial Expressions in the Wild (AFEW) dataset as stimuli. To establish a baseline, we use the electrodermal activity (EDA) signals in this dataset and extract 6 features from each signal series corresponding to each movie clip to recognize 7 emotions, i.e., Anger, Disgust, Fear, Happy, Surprise, Sad and Neutral. Overall, 24 observers participated in our collection of the training set, including 19 observers who participated in only one session watching 80 videos from 7 classes and 5 observers who participated multiple times and watched all the videos. All videos were presented in an order balanced fashion. Leave-one-observer-out was employed in this classification task. We report the classification accuracy of our baseline, a three-layer network, on this initial training set while training with signals from all participants, only single participants and only multiple participants. We also investigate the recognition accuracy of grouping the dataset by arousal or valence, which achieves 68.66% and 72.72% separately. Finally, we provide a two-step network. The first step is to classify the features into high/low arousal or positive/negative valence by a network. Then the arousal/valence middle output of the first step is concatenated with feature sets as input of the second step for emotion recognition. We found that adding arousal or valence information can help to improve the classification accuracy. In addition, the information of positive/negative valence boosts the classification accuracy to a higher degree on this dataset.


title: A Case for Humans-in-the-Loop: Decisions in the Presence of Erroneous Algorithmic Scores

abstract: The increased use of algorithmic predictions in sensitive domains has been accompanied by both enthusiasm and concern. To understand the opportunities and risks of these technologies, it is key to study how experts alter their decisions when using such tools. In this paper, we study the adoption of an algorithmic tool used to assist child maltreatment hotline screening decisions. We focus on the question: Are humans capable of identifying cases in which the machine is wrong, and of overriding those recommendations? We first show that humans do alter their behavior when the tool is deployed. Then, we show that humans are less likely to adhere to the machine's recommendation when the score displayed is an incorrect estimate of risk, even when overriding the recommendation requires supervisory approval. These results highlight the risks of full automation and the importance of designing decision pipelines that provide humans with autonomy.


title: Architecture of three-planet systems predicted from the observed protoplanetary disk of HL Tau

abstract: A number of protoplanetary disks observed with ALMA potentially provide direct examples of initial conditions for planetary systems. In particular, the HL Tau disk has been intensively studied, and its rings/gaps are conventionally interpreted to be a result of unseen massive planets embedded in the gaps. Based on this interpretation, we carried out N-body simulations to investigate orbital evolution of planets within the protoplanetary disk and after the disk dispersal. Before the disk dispersal, our N-body simulations include both migration and mass-growth of the planet coupled with evolution of the disk. By varying the disk parameters, we produce a variety of widely-separated planetary systems consisting of three super-Jupiters at the end of disk dispersal. We found the outer planet is more massive than the inner one, and the migration of the innermost planet is inefficient due to the accretion of outer planet(s). We also showed how the final configuration and the final planetary mass depend on disk parameters. The migration is found to be convergent and no planet-pair has a period ratio less than 2. After the disk dispersal, we switch to pure gravitational N-body simulations and integrate the orbits up to 10 Gyr. Most simulated systems remain stable for at least 10 Gyr. We discuss implications of our result in terms of the observed widely-separated planetary systems HR 8799 and PDS 70.


title: Efficient Deep Reinforcement Learning via Adaptive Policy Transfer

abstract: Transfer Learning (TL) has shown great potential to accelerate Reinforcement Learning (RL) by leveraging prior knowledge from past learned policies of relevant tasks. Existing transfer approaches either explicitly computes the similarity between tasks or select appropriate source policies to provide guided explorations for the target task. However, how to directly optimize the target policy by alternatively utilizing knowledge from appropriate source policies without explicitly measuring the similarity is currently missing. In this paper, we propose a novel Policy Transfer Framework (PTF) to accelerate RL by taking advantage of this idea. Our framework learns when and which source policy is the best to reuse for the target policy and when to terminate it by modeling multi-policy transfer as the option learning problem. PTF can be easily combined with existing deep RL approaches. Experimental results show it significantly accelerates the learning process and surpasses state-of-the-art policy transfer methods in terms of learning efficiency and final performance in both discrete and continuous action spaces.


title: Statistical Inversion Using Sparsity and Total Variation Prior And Monte Carlo Sampling Method For Diffuse Optical Tomography

abstract: In this paper, we formulate the reconstruction problem in diffuse optical tomography (DOT) in a statistical setting for determining the optical parameters, scattering and absorption, from boundary photon density measurements. A special kind of adaptive Metropolis algorithm for the reconstruction procedure using sparsity and total variation prior is presented. Finally, a simulation study of this technique with different regularization functions and its comparison to the deterministic Iteratively Regularized Gauss Newton method shows the effectiveness and stability of the method.


title: Feasibility of Video-based Sub-meter Localization on Resource-constrained Platforms

abstract: While the satellite-based Global Positioning System (GPS) is adequate for some outdoor applications, many other applications are held back by its multi-meter positioning errors and poor indoor coverage. In this paper, we study the feasibility of real-time video-based localization on resource-constrained platforms. Before commencing a localization task, a video-based localization system downloads an offline model of a restricted target environment, such as a set of city streets, or an indoor shopping mall. The system is then able to localize the user within the model, using only video as input. To enable such a system to run on resource-constrained embedded systems or smartphones, we (a) propose techniques for efficiently building a 3D model of a surveyed path, through frame selection and efficient feature matching, (b) substantially reduce model size by multiple compression techniques, without sacrificing localization accuracy, (c) propose efficient and concurrent techniques for feature extraction and matching to enable online localization, (d) propose a method with interleaved feature matching and optical flow based tracking to reduce the feature extraction and matching time in online localization. Based on an extensive set of both indoor and outdoor videos, manually annotated with location ground truth, we demonstrate that sub-meter accuracy, at real-time rates, is achievable on smart-phone type platforms, despite challenging video conditions.


title: Cellular UAV-to-Device Communications: Trajectory Design and Mode Selection by Multi-agent Deep Reinforcement Learning

abstract: In the current unmanned aircraft systems (UASs) for sensing services, unmanned aerial vehicles (UAVs) transmit their sensory data to terrestrial mobile devices over the unlicensed spectrum. However, the interference from surrounding terminals is uncontrollable due to the opportunistic channel access. In this paper, we consider a cellular Internet of UAVs to guarantee the Quality-of-Service (QoS), where the sensory data can be transmitted to the mobile devices either by UAV-to-Device (U2D) communications over cellular networks, or directly through the base station (BS). Since UAVs' sensing and transmission may influence their trajectories, we study the trajectory design problem for UAVs in consideration of their sensing and transmission. This is a Markov decision problem (MDP) with a large state-action space, and thus, we utilize multi-agent deep reinforcement learning (DRL) to approximate the state-action space, and then propose a multi-UAV trajectory design algorithm to solve this problem. Simulation results show that our proposed algorithm can achieve a higher total utility than policy gradient algorithm and single-agent algorithm.


title: Enlarging Discriminative Power by Adding an Extra Class in Unsupervised Domain Adaptation

abstract: In this paper, we study the problem of unsupervised domain adaptation that aims at obtaining a prediction model for the target domain using labeled data from the source domain and unlabeled data from the target domain. There exists an array of recent research based on the idea of extracting features that are not only invariant for both domains but also provide high discriminative power for the target domain. In this paper, we propose an idea of empowering the discriminativeness: Adding a new, artificial class and training the model on the data together with the GAN-generated samples of the new class. The trained model based on the new class samples is capable of extracting the features that are more discriminative by repositioning data of current classes in the target domain and therefore drawing the decision boundaries more effectively. Our idea is highly generic so that it is compatible with many existing methods such as DANN, VADA, and DIRT-T. We conduct various experiments for the standard data commonly used for the evaluation of unsupervised domain adaptations and demonstrate that our algorithm achieves the SOTA performance for many scenarios.


title: How to choose master integrals

abstract: The standard procedure when evaluating integrals of a given family of Feynman integrals, corresponding to some Feynman graph, is to construct an algorithm which provides the possibility to write any particular integral as a linear combination of so-called master integrals. To do this, public (AIR, FIRE, REDUZE, LiteRed, KIRA) and private codes based on solving integration by parts relations are used. However, the choice of the master integrals provided by these codes is not always optimal. We present an algorithm to improve a given basis of the master integrals, as well as its computer implementation; see also a competitive variant [1].


title: Meta Segmentation Network for Ultra-Resolution Medical Images

abstract: Despite recent progress on semantic segmentation, there still exist huge challenges in medical ultra-resolution image segmentation. The methods based on multi-branch structure can make a good balance between computational burdens and segmentation accuracy. However, the fusion structure in these methods require to be designed elaborately to achieve desirable result, which leads to model redundancy. In this paper, we propose Meta Segmentation Network (MSN) to solve this challenging problem. With the help of meta-learning, the fusion module of MSN is quite simple but effective. MSN can fast generate the weights of fusion layers through a simple meta-learner, requiring only a few training samples and epochs to converge. In addition, to avoid learning all branches from scratch, we further introduce a particular weight sharing mechanism to realize a fast knowledge adaptation and share the weights among multiple branches, resulting in the performance improvement and significant parameters reduction. The experimental results on two challenging ultra-resolution medical datasets BACH and ISIC show that MSN achieves the best performance compared with the state-of-the-art methods.


title: Relaxed Gauss-Newton methods with applications to electrical impedance tomography

abstract: As second-order methods, Gauss--Newton-type methods can be more effective than first-order methods for the solution of nonsmooth optimization problems with expensive-to-evaluate smooth components. Such methods, however, often do not converge. Motivated by nonlinear inverse problems with nonsmooth regularization, we propose a new Gauss--Newton-type method with inexact relaxed steps. We prove that the method converges to a set of disjoint critical points given that the linearisation of the forward operator for the inverse problem is sufficiently precise. We extensively evaluate the performance of the method on electrical impedance tomography (EIT).


title: Sharp Weak Bounds for p-adic Hardy operators on p-adic Linear Spaces

abstract: In this article, we establish sharp weak endpoint estimates for $p$-adic fractional Hardy operator. In addition, sharp weak bounds for p-adic Hardy operator on p-adic central Morrey space are also obtained.


title: Tree-structured Attention with Hierarchical Accumulation

abstract: Incorporating hierarchical structures like constituency trees has been shown to be effective for various natural language processing (NLP) tasks. However, it is evident that state-of-the-art (SOTA) sequence-based models like the Transformer struggle to encode such structures inherently. On the other hand, dedicated models like the Tree-LSTM, while explicitly modeling hierarchical structures, do not perform as efficiently as the Transformer. In this paper, we attempt to bridge this gap with "Hierarchical Accumulation" to encode parse tree structures into self-attention at constant time complexity. Our approach outperforms SOTA methods in four IWSLT translation tasks and the WMT'14 English-German translation task. It also yields improvements over Transformer and Tree-LSTM on three text classification tasks. We further demonstrate that using hierarchical priors can compensate for data shortage, and that our model prefers phrase-level attentions over token-level attentions.


title: Nonequilibrium effects in the Casimir force between two similar metallic plates kept at different temperatures

abstract: We study the Casimir pressure between two similar plates of finite thickness kept at different temperatures in the case when the dielectric permittivity of the plates depends on temperature. It is suggested to consider the dielectric permittivity at two different temperatures as the permittivities of two dissimilar bodies, thus allowing to apply the theory of Casimir forces out of thermal equilibrium developed earlier in the literature. Following this approach, we show that, in addition to the equilibrium contribution to the nonequilibrium Casimir pressure, a proper nonequilibrium contribution arises for temperature-dependent dielectric permittivities. Furthermore, the equilibrium contribution in this case does not equal the mean of the equilibrium Casimir pressures at the temperatures of the plates. As an application, the total nonequilibrium Casimir pressure between two gold plates and between two titanium plates is calculated as a function of the plate thickness and their separation, using the Drude and the plasma model. For plate separations ranging from 0.5 to 2 micrometers, the relative difference between the theoretical predictions for these two models reaches 39%. The proper nonequilibrium term may be as large as 4% of the magnitude of the total nonequilibrium pressure.


title: Towards the merger of Hawking radiating black holes

abstract: We discuss the merger process of binary black holes with Hawking radiation taken into account. Besides the redshifted radiation to infinity, binary black holes can exchange radiation between themselves, which is first redshifted and then blueshifted when it propagates from one hole to the other. The exchange rate should be large when the temperature-divergent horizons are penetrating each other to form a single horizon with unique temperature. This will cause non-negligible mass and angular momentum transfer between the black holes during the merging process of the horizons. We further argue in the large mass ratio limit that the light hole whose local evaporation is enhanced by the competing redshift-blueshift effects will probably evaporate or decay completely before reaching the the horizon of the heavy one. We also discuss the possibility of testing Hawking radiation and even exploring the information loss puzzle in gravitational wave observations.


title: The uniqueness of covers for widely generalized line graphs

abstract: As a natural generalization of line graphs, Hoffman line graphs were defined by Woo and Neumaier. Especially, Hoffman line graphs are closely related to the smallest eigenvalue of graphs, and the uniqueness of strict covers of a Hoffman line graph plays a key role in such a study. In this paper, we prove a theorem for the uniqueness of strict covers under a condition which can be checked in finite time. Our result gives a generalization and a short proof for the main part of [Ars Math.~Contemp. \textbf{1} (2008) 81--98].


title: Formation of low-spinning 100 Msun black holes

abstract: It is speculated that a merger of two massive stellar-origin BHs in a dense stellar environment may lead to the formation of a massive BH in the pair-instability mass gap (50-135 Msun). Such a merger-formed BH is expected to typically have a high spin (a=0.7). If such a massive BH acquires another BH it may lead to another merger detectable by LIGO/Virgo in gravitational waves. Acquiring a companion may be hindered by gravitational-wave kick/recoil, which accompanies the first merger and may quickly remove the massive BH from its parent globular or nuclear cluster. We test whether it is possible for a massive merger-formed BH in the pair-instability gap to be retained in its parent cluster and have low spin. Such a BH would be indistinguishable from a primordial BH. We employed results from numerical relativity calculations of black hole mergers to explore the range of gravitational-wave recoil velocities for various combinations of merging BH masses and spins. We compared merger-formed massive BH speeds with typical escape velocities from globular and nuclear clusters. We show that a globular cluster is highly unlikely to form and retain a 100 Msun BH if the spin of the BH is low (a<0.3) as such BHs acquire high recoil speeds (>200 km/s) that exceed typical escape speeds from globular clusters (50 km/s). However, a very low-spinning (a=0.1) and massive (100 Msun) BH could be formed and retained in a galactic nuclear star cluster. Even though such massive merger-formed BHs with such low spins acquire high speeds during formation (400 km/s), they may avoid ejection since massive nuclear clusters have high escape velocities (300-500 km/s). A future detection of a massive BH in the pair-instability mass gap with low spin would therefore not be proof of the existence of primordial BHs, which are sometimes claimed to have low spins and arbitrarily high masses.


title: Selected strong decays of pentaquark State $P_c(4312)$ in a chiral constituent quark model

abstract: The newly confirmed pentaquark state $P_c(4312)$ has been treated as a weakly bound $(\Sigma_c\bar{D})$ state by a well-established chiral constituent quark model and by a dynamical calculation on quark degrees of freedom where the quark exchange effect is accounted for. The obtained mass $4308$ MeV agrees with data. In this work, the selected strong decays of the $P_c(4312)$ state are studied with the obtained wave function. It is shown that the width of the $\Lambda_c\bar{D}^*$ decay is overwhelmed and the branching ratios of the $p\,\eta_c$ and $p\,J/\psi$ decays are both less than 1 percentage.


title: Effective Hamiltonian model for helically constrained quantum systems within adiabatic perturbation theory: application to the Chirality-Induced Spin Selectivity (CISS) Effect

abstract: The chirality-induced spin selectivity (CISS) effect has been confirmed experimentally for a large class of organic molecules. Adequately modeling the effect remains a challenging task, with both phenomenological models and first-principle simulations yielding inconclusive results. Building upon a previously presented model by K. Michaeli and R. Naaman (J. Phys. Chem C 123, 17043 (2019)) we systematically investigate an effective 1-dimensional model derived as the limit of a 3-dimensional quantum system with strong confinement and including spin-orbit coupling. Having a simple analytic structure, such models can be considered a minimal setup for the description of spin-dependent effects. We use adiabatic perturbation theory to provide a mathematically sound approximation procedure applicable to a large class of spin-dependent continuum models. We take advantage of the models simplicity by analyzing its structure to gain a better understanding how the occurrence and magnitude of spin polarization effects relate to the model's parameters and geometry.


title: Progressive Identification of True Labels for Partial-Label Learning

abstract: Partial-label learning (PLL) is a typical weakly supervised learning problem, where each training instance is equipped with a set of candidate labels among which only one is the true label. Most existing methods elaborately designed learning objectives as constrained optimizations that must be solved in specific manners, making their computational complexity a bottleneck for scaling up to big data. The goal of this paper is to propose a novel framework of PLL with flexibility on the model and optimization algorithm. More specifically, we propose a novel estimator of the classification risk, theoretically analyze the classifier-consistency, and establish an estimation error bound. Then we propose a progressive identification algorithm for approximately minimizing the proposed risk estimator, where the update of the model and identification of true labels are conducted in a seamless manner. The resulting algorithm is model-independent and loss-independent, and compatible with stochastic optimization. Thorough experiments demonstrate it sets the new state of the art.


title: Majority Logic Decoding for Certain Schubert Codes Using Lines in Schubert Varieties

abstract: In this article, we consider Schubert codes, linear codes associated to Schubert varieties, and discuss minimum weight codewords for dual Schubert codes. The notion of lines in Schubert varieties is looked closely at, and it has been proved that the supports of the minimum weight codewords of the dual Schubert codes lie on lines and any three points on a line in Schubert variety correspond to the support of some minimum weight parity check for the Schubert code. We use these lines in Schubert varieties to construct orthogonal parity checks for certain Schubert codes and use them for majority logic decoding. In some special cases, we can correct approximately up to $\lfloor (d-1)/2\rfloor$ many errors where $d$ is the minimum distance of the code.


title: Bilinear spherical maximal functions of product type

abstract: In this paper we introduce and study a bilinear spherical maximal function of product type in the spirit of bilinear Calder\'{o}n-Zygmund theory. This operator is different from the bilinear spherical maximal function considered by Geba et al. We deal with lacunary and full versions of this operator, and we prove weighted estimates with respect to bilinear weights. Our approach involves sparse forms following ideas by Lacey, but we also use other techniques to handle the particular triplet $(2,2,1)$.


title: The Geometry of Sign Gradient Descent

abstract: Sign-based optimization methods have become popular in machine learning due to their favorable communication cost in distributed optimization and their surprisingly good performance in neural network training. Furthermore, they are closely connected to so-called adaptive gradient methods like Adam. Recent works on signSGD have used a non-standard "separable smoothness" assumption, whereas some older works study sign gradient descent as steepest descent with respect to the $\ell_\infty$-norm. In this work, we unify these existing results by showing a close connection between separable smoothness and $\ell_\infty$-smoothness and argue that the latter is the weaker and more natural assumption. We then proceed to study the smoothness constant with respect to the $\ell_\infty$-norm and thereby isolate geometric properties of the objective function which affect the performance of sign-based methods. In short, we find sign-based methods to be preferable over gradient descent if (i) the Hessian is to some degree concentrated on its diagonal, and (ii) its maximal eigenvalue is much larger than the average eigenvalue. Both properties are common in deep networks.


title: Positive entropy using Hecke operators at a single place

abstract: We prove the following statement: Let $X=\text{SL}_n(\mathbb{Z})\backslash \text{SL}_n(\mathbb{R})$, and consider the standard action of the diagonal group $A<\text{SL}_n(\mathbb{R})$ on it. Let $\mu$ be an $A$-invariant probability measure on $X$, which is a limit $$ \mu=\lambda\lim_i|\phi_i|^2dx, $$ where $\phi_i$ are normalized eigenfunctions of the Hecke algebra at some fixed place $p$, and $\lambda>0$ is some positive constant. Then any regular element $a\in A$ acts on $\mu$ with positive entropy on almost every ergodic component. We also prove a similar result for lattices coming from division algebras over $\mathbb{Q}$, and derive a quantum unique ergodicity result for the associated locally symmetric spaces. This generalizes a result of Brooks and Lindenstrauss.


title: Exploiting characteristics in stationary action problems

abstract: Connections between the principle of stationary action and optimal control, and between established notions of minimax and viscosity solutions, are combined to describe trajectories of energy conserving systems as solutions of corresponding Cauchy problems defined with respect to attendant systems of characteristics equations.


title: Nanoscale magnetic resonance imaging of proteins in a single cell

abstract: Magnetic resonance imaging (MRI) is a non-invasive and label-free technique widely used in medical diagnosis and life science research, and its success has benefited greatly from continuing efforts on enhancing contrast and resolution. Here we reported nanoscale MRI in a single cell using an atomic-size quantum sensor. With nitrogen-vacancy center in diamond, the intracellular protein ferritin has been imaged with a spatial resolution of ~ 10 nanometers, and ferritin-containing organelles were co-localized by correlative MRI and electron microscopy. Comparing to the current micrometer resolution in current state-of-art conventional MRI, our approach represents a 100-fold enhancement, and paves the way for MRI of intracellular proteins.


title: Simulataneous Control of Wave Systems

abstract: In this paper, we study the simultaneous controllability of wave systems in an open domain of R d , d $\in$ N *. We obtain a partial controllability result on a co-finite dimensional space for wave equations coupled by a single control function. We use microlocal defect measures and the unique continuation property of eigenfunctions to prove that an appropriate observability inequality holds for wave equations with space varying and different speeds coupled by a single control function. For the unique continuation property of eigenfunctions, we construct a counterexample to show that in some metrics, the unique continuation property does not hold. Moreover, we study different conditions to ensure the unique continuation property. We also extend our result to the case of constant coefficients and possibly multiple control functions. In this context, we prove the controllability property is equivalent to an appropriate Kalman rank condition.


title: Translating Between Wavelet Tree and Wavelet Matrix Construction

abstract: The wavelet tree (Grossi et al. [SODA, 2003]) and wavelet matrix (Claude et al. [Inf. Syst., 2015]) are compact data structures with many applications such as text indexing or computational geometry. By continuing the recent research of Fischer et al. [ALENEX, 2018], we explore the similarities and differences of these heavily related data structures with focus on their construction. We develop a data structure to modify construction algorithms for either the wavelet tree or matrix to construct instead the other. This modification is efficient, in that it does not worsen the asymptotic time and space requirements of any known wavelet tree or wavelet matrix construction algorithm.


title: Primality tests, linear recurrent sequences and the Pell equation

abstract: We study new primality tests based on linear recurrent sequences of degree two exploiting a matricial approach. The classical Lucas test arises as a particular case and we see how it can be easily improved. Moreover, this approach shows clearly how the Lucas pseudoprimes are connected to the Pell equation and the Brahamagupta product. We also introduce a new specific primality test, which we will call generalized Pell test. We perform some numerical computations on the new primality tests and, for the generalized Pell test, we do not any pseudoprime up to $10^{10}$.


title: Geometrical Optical Illusion via Sub-Riemannian Geodesics in the Roto-Translation Group

abstract: We present a neuro-mathematical model for geometrical optical illusions (GOIs), a class of illusory phenomena that consists in a mismatch of geometrical properties of the visual stimulus and its associated percept. They take place in the visual areas V1/V2 whose functional architecture have been modelled in previous works by Citti and Sarti as a Lie group equipped with a sub-Riemannian (SR) metric. Here we extend their model proposing that the metric responsible for the cortical connectivity is modulated by the modelled neuro-physiological response of simple cells to the visual stimulus, hence providing a more biologically plausible model that takes into account a presence of visual stimulus. Illusory contours in our model are described as geodesics in the new metric. The model is confirmed by numerical simulations, where we compute the geodesics via SR-Fast Marching.


title: Distributed Algorithms that Solve Boolean Equations with Local and Differential Privacies

abstract: In this paper, we propose distributed algorithms that solve a system of Boolean equations over a network, where each node in the network possesses only one Boolean equation from the system. The Boolean equation assigned at any particular node is a {\em private} equation known to this node only, and the nodes aim to compute the exact set of solutions to the system without exchanging their local equations. We show that each private Boolean equation can be locally lifted to a linear algebraic equation under a basis of Boolean vectors, leading to a network linear equation that is distributedly solvable using existing distributed linear equation algorithms as a subroutine. A number of exact or approximate solutions to the induced linear equation are then computed at each node from different initial values. The solutions to the original Boolean equations are eventually computed locally via a Boolean vector search algorithm. We prove that given solvable Boolean equations, when the initial values of the nodes for the distributed linear equation solving step are i.i.d selected according to a uniform distribution in a high-dimensional cube, our algorithms return the exact solution set of the Boolean equations at each node with high probability. Furthermore, we present an algorithm for distributed verification of the satisfiability of Boolean equations, and prove its correctness. Finally, we show that by utilizing linear equation solvers with differential privacy to replace the in-network computing routines, the overall distributed Boolean equation algorithms can be made differentially private. Under the standard Laplace mechanism, we prove an explicit level of noises that can be injected in the linear equation steps for ensuring a prescribed level of differential privacy.


title: Comments on "Spatio-Temporal Gaussian Process Models for Extended and Group Object Tracking with Irregular Shapes"

abstract: In the study "Spatio-temporal Gaussian process models for extended and group object tracking with irregular shapes" (IEEE Trans. Veh. Tech., vol. 68, no. 3, pp. 2137-2151, Mar. 2019), the extended object tracking problem was tackled by an approach based on spatio-temporal Gaussian processes (STGP). The performance of the proposed STGP-based trackers was comparatively evaluated through simulations and real data together with another state-of-the-art method (referred to as GP-EKF) proposed in "Extended target tracking using Gaussian processes" (IEEE Trans. Signal Process., vol. 63, no. 16, pp. 4165-4178, Aug. 2015). Unfortunately, we recognized that there are major errors in the implementation of the experiments presented in the STGP paper, which led to incorrect performance evaluation results. In this correspondence, our aim is to share the correct results of these experiments and to respond to some claims regarding GP-EKF, which we believe, would contribute to a better understanding of the methods.


title: Light-matter interaction in van der Waals heterostructures

abstract: Even if individual two-dimensional materials own various interesting and unexpected properties, the stacking of such layers leads to van der Waals solids which unite the characteristics of two dimensions with novel features originating from the interlayer interactions. In this topical review, we cover fabrication and characterization of van der Waals heterosructures with a focus on heterobilayers made of monolayers of semiconducting transition metal dichalcogenides. Experimental and theoretical techniques to investigate those heterobilayers are introduced. Most recent findings focusing on different transition metal dichalcogenides heterostructures are presented and possible optical transitions between different valleys, appearance of moire patterns and signatures of moire excitons are discussed. The fascinating and fast growing research on van der Waals hetero-bilayers provide promising insights required for their application as emerging quantum-nano materials.


title: Generating function for Naturalized Series: The case of Ordered Motzkin Words

abstract: We continue to consider the ordered lexicographic sequence, which is constructed according to the formal characteristics of a series of natural numbers. For analysis, we selected balanced parentheses with zeros, Motzkin words. As you know, generating functions allow you to work with combinatorial objects by analytical methods. Motzkin words are enumerated by Motzkin numbers, for the generation of which there is a corresponding generating function. In our case, restrictions are imposed on Motzkin words, for example, there are no leading zeros in bracket sets. The purpose of this article is to obtain the generating function of such modified Motzkin words.


title: The convex invertible cone structure of positive real odd rational matrix functions

abstract: Positive real odd matrix functions, often referred to as positive real lossless matrix functions, play an important role in many applications in multi-port electrical systems. In this paper we present closer analogues to some of the known results for the scalar, one-port, case in the multi-port setting. Specifically, we determine necessary and sufficient conditions for the well studied partial fraction formula to represent functions in the class of positive real odd matrix functions, and explicit minimal state space realization formulas for the inverse (admittance) of a function in this class, which itself is also a positive real odd matrix function. Doing so, enables us to provide a partial analogue of the poles-zero interlacing behavior from the scalar case.


title: White Dwarfs as Advanced Physics Laboratories. The Axion case

abstract: The shape of the luminosity function of white dwarfs (WDLF) is sensitive to the characteristic cooling time and, therefore, it can be used to test the existence of additional sources or sinks of energy such as those predicted by alternative physical theories. However, because of the degeneracy between the physical properties of white dwarfs and the properties of the Galaxy, the star formation history (SFH) and the IMF, it is almost always possible to explain any anomaly as an artifact introduced by the star formation rate. To circumvent this problem there are at least two possibilities, the analysis of the WDLF in populations with different stories, like disc and halo, and the search of effects not correlated with the SFH. These procedures are illustrated with the case of axions.


title: Complete classification of trapping coins for quantum walks on the 2D square lattice

abstract: One of the unique features of discrete-time quantum walks is called trapping, meaning the inability of the quantum walker to completely escape from its initial position, albeit the system is translationally invariant. The effect is dependent on the dimension and the explicit form of the local coin. A four state discrete-time quantum walk on a square lattice is defined by its unitary coin operator, acting on the four dimensional coin Hilbert space. The well known example of the Grover coin leads to a partial trapping, i.e., there exists some escaping initial state for which the probability of staying at the initial position vanishes. On the other hand, some other coins are known to exhibit strong trapping, where such escaping state does not exist. We present a systematic study of coins leading to trapping, explicitly construct all such coins for discrete-time quantum walks on the 2D square lattice, and classify them according to the structure of the operator and the manifestation of the trapping effect. We distinguish three types of trapping coins exhibiting distinct dynamical properties, as exemplified by the existence or non-existence of the escaping state and the area covered by the spreading wave-packet.


title: Dissecting Neural ODEs

abstract: Continuous deep learning architectures have recently re-emerged as Neural Ordinary Differential Equations (Neural ODEs). This infinite-depth approach theoretically bridges the gap between deep learning and dynamical systems, offering a novel perspective. However, deciphering the inner working of these models is still an open challenge, as most applications apply them as generic black-box modules. In this work we "open the box", further developing the continuous-depth formulation with the aim of clarifying the influence of several design choices on the underlying dynamics.


title: The Stellar Variability Noise Floor for Transiting Exoplanet Photometry with PLATO

abstract: One of the main science motivations for the ESA PLAnetary Transit and Oscillations (PLATO) mission is to measure exoplanet transit radii with 3% precision. In addition to flares and starspots, stellar oscillations and granulation will enforce fundamental noise floors for transiting exoplanet radius measurements. We simulate light curves of Earth-sized exoplanets transiting continuum intensity images of the Sun taken by the HMI instrument aboard SDO to investigate the uncertainties introduced on the exoplanet radius measurements by stellar granulation and oscillations. After modeling the solar variability with a Gaussian process, we find that the amplitude of solar oscillations and granulation is of order 100 ppm -- similar to the depth of an Earth transit -- and introduces a fractional uncertainty on the depth of transit of 0.73% assuming four transits are observed over the mission duration. However, when we translate the depth measurement into a radius measurement of the planet, we find a much larger radius uncertainty of 3.6%. This is due to a degeneracy between the transit radius ratio, the limb-darkening, and the impact parameter caused by the inability to constrain the transit impact parameter in the presence of stellar variability. We find that surface brightness inhomogeneity due to photospheric granulation contributes a lower limit of only 2 ppm to the photometry in-transit. The radius uncertainty due to granulation and oscillations, combined with the degeneracy with the transit impact parameter, accounts for a significant fraction of the error budget of the PLATO mission, before detector or observational noise is introduced to the light curve. If it is possible to constrain the impact parameter or to obtain follow-up observations at longer wavelengths where limb-darkening is less significant, this may enable higher precision radius measurements.


title: Continuity from neutron matter to color-superconducting quark matter with $^3 P_2$ superfluidity

abstract: I clarify how the concept of quark-hadron continuity, which was previously considered in the context of the asymptotic color-flavor locked phase with ideal three-flavor symmetry, is applied to two-flavor matter. Our observation is that neutron star matter can continuously be connected to two-flavor color-superconducting (2SC) phase with an additional condensate of $d$-quarks in the $^3P_2$ channel. I discuss here two aspects of this novel phase. First, I introduce the notion of continuity based on the patterns of symmetry breaking and the corresponding order parameters, and then explain qualitatively the physics mechanism of $d$-quark $^3P_2$ pairing in analogy to nuclear physics. Our finding serves as the theoretical underpinning for the phenomenological construction of the equation of state in the neutron star environment.


title: A note on uniquely 10-colorable graphs

abstract: We prove for k at most 10, that every graph of chromatic number k with a unique k-coloring admits a clique minor of order k.


title: Extended calculations of energy levels and transition rates for singly ionized lanthanide elements I: Pr - Gd

abstract: Lanthanide elements play important roles as an opacity source in the ejected material from neutron star mergers. Accurate and complete atomic data are necessary to evaluate the opacities and to analyze the observed data. In this paper, we perform extended, ab-initio atomic calculations from Pr II (Z=59) to Gd II (Z=64). By using multiconfiguration Dirac-Hartree-Fock and relativistic configuration-interaction methods, implemented in the general-purpose relativistic atomic structure package GRASP2K, we calculate the energy levels and transition data of electric dipole transitions. These computations are based on strategies (with small variations) of Nd II published by Gaigalas et al. (2019). Accuracy of data is evaluated by comparing computed energy levels with the NIST database or other works. For the energy levels, we obtain the average relative accuracy of 8%, 12%, 6%, 8%, and 7% for Pr II, Pm II, Sm II, Eu II, and Gd II ions, respectively as compared with the NIST data. Accuracy of energy transfer to the wavelength as 3%, 14% and 11% for Pr II, Eu II and Gd II. Our computed E1 type transition probabilities are in good agreement with experimental values presented by other authors especially for strong transitions.


title: Pressure effect on the order-disorder transformation in L1$_0$ FeNi

abstract: The ordered phase of the FeNi system is known to have promising magnetic properties as a rare-earth-free permanent magnet. Understanding the parameter space that controls the order-disorder transformation is important to find growth conditions that stabilize the $L1_0$ phase. Magnetic properties and chemical order-disorder transformation of FeNi are investigated as a function of lattice expansion using first-principles theory. Thermodynamic and magnetic calculations are performed using long-range order parameters, which is used to find the ordering temperature. Negative pressure promotes ordering and thus synthetic routes involving an increase in volume is expected to expand the stability field of the $L1_0$ phase.


title: The role of disc torques in forming resonant planetary systems

abstract: The most accurate method for modelling planetary migration and hence the formation of resonant systems is using hydrodynamical simulations. Usually, the force (torque) acting on a planet is calculated using the forces from the gas disc and the star, while the gas accelerations are computed using the pressure gradient, the star, and the planet's gravity, ignoring its own gravity. For the non-migrating the neglect of the disc gravity results in a consistent torque calculation while for the migrating case it is inconsistent. We aim to study how much this inconsistent torque calculation can affect the final configuration of a two-planet system. Our focus will be on low-mass planets because most of the multi-planetary systems, discovered by the Kepler survey, have masses around 10 Earth masses. Performing hydrodynamical simulations of planet-disc interaction, we measure the torques on non-migrating and migrating planets for various disc masses as well as density and temperature slopes with and without considering the disc self-gravity. Using this data, we find a relation that quantifies the inconsistency, use it in an N-body code, and perform an extended parameter study modelling the migration of a planetary system with different planet mass ratios and disc surface densities, in order to investigate the impact of the torque inconsistency on the architecture of the planetary system. Not considering disc self-gravity produces an artificially larger torque on the migrating planet that can result in tighter planetary systems. The deviation of this torque from the correct value is larger in discs with steeper surface density profiles. In hydrodynamical modelling of multi-planetary systems, it is crucial to account for the torque correction, otherwise the results favour more packed systems.


title: EXTraS discovery of an X-ray superflare from an L dwarf

abstract: We present the first detection of an X-ray flare from an ultracool dwarf of spectral class L. The event was identified in the EXTraS database of XMM-Newton variable sources, and its optical counterpart, J0331-27, was found through a cross-match with the Dark Energy Survey Year 3 release. Next to an earlier four-photon detection of Kelu-1, J0331-27 is only the second L dwarf detected in X-rays, and much more distant than other ultracool dwarfs with X-ray detections (photometric distance of 240 pc). From an optical spectrum with the VIMOS instrument at the VLT, we determine the spectral type of J0331-27 to be L1. The X-ray flare has an energy of E_X,F ~ 2x10^33 erg, placing it in the regime of superflares. No quiescent emission is detected, and from 2.5 Msec of XMM data we derive an upper limit of L_X,qui < 10^27 erg/s. The flare peak luminosity L_X,peak = 6.3x10^29 erg/s, flare duration tau_decay ~ 2400 s, and plasma temperature (~16 MK) are similar to values observed in X-ray flares of M dwarfs. This shows that strong magnetic reconnection events and the ensuing plasma heating are still present even in objects with photospheres as cool as ~2100 K. However, the absence of any other flares above the detection threshold of E_X,F ~2.5x10^32 erg in a total of ~2.5 Ms of X-ray data yields a flare energy number distribution inconsistent with the canonical power law dN/dE ~ E^-2, suggesting that magnetic energy release in J0331-27 -- and possibly in all L dwarfs -- takes place predominantly in the form of giant flares.


title: History-dependent percolation on multiplex networks

abstract: The structure of interconnected systems and its impact on the system dynamics is a much-studied cross-disciplinary topic. Although various critical phenomena have been found in different models, the study on the connections between different percolation transitions is still lacking. Here we propose a unified framework to study the origins of the discontinuous transitions of the percolation process on interacting networks. The model evolves in generations with the result of the present percolation depending on the previous state and thus is history-dependent. Both theoretical analysis and Monte Carlo simulations reveal that the nature of the transition remains the same at finite generations but exhibits an abrupt change for the infinite generation. We use brain functional correlation and morphological similarity data to show that our model also provides a general method to explore the network structure and can contribute to many practical applications, such as detecting the abnormal structures of human brain networks.


title: Dynamic Coded Caching in Wireless Networks

abstract: We consider distributed and dynamic caching of coded content at small base stations (SBSs) in an area served by a macro base station (MBS). Specifically, content is encoded using a maximum distance separable code and cached according to a time-to-live (TTL) cache eviction policy, which allows coded packets to be removed from the caches at periodic times. Mobile users requesting a particular content download coded packets from SBSs within communication range. If additional packets are required to decode the file, these are downloaded from the MBS. We formulate an optimization problem that is efficiently solved numerically, providing TTL caching policies minimizing the overall network load. We demonstrate that distributed coded caching using TTL caching policies can offer significant reductions in terms of network load when request arrivals are bursty. We show how the distributed coded caching problem utilizing TTL caching policies can be analyzed as a specific single cache, convex optimization problem. Our problem encompasses static caching and the single cache as special cases. We prove that, interestingly, static caching is optimal under a Poisson request process, and that for a single cache the optimization problem has a surprisingly simple solution.


title: Lattice dynamics and spin-phonon interaction in strained NiO films

abstract: NiO thin films with various strains were grown on SrTiO3 (STO) and MgO substrates using a pulsed laser deposition technique. The films were characterized using an x-ray diffractometer, atomic force microscopy, and infrared reflectance spectroscopy. The films grown on STO (001) substrate show a compressive in-plane strain which increases as the film thickness is reduced, resulting in an increase of the NiO phonon frequency. On the other hand, a tensile strain was detected in the NiO film grown on MgO (001) substrate which induces a softening of the phonon frequency. Overall, the variation of in-plane strain from -0.36% to +0.48% yields the decrease of the phonon frequency from 409.6 cm-1 to 377.5 cm-1 which occurs due to the ~1% change of the inter-atomic distances. The magnetic exchange -driven phonon splitting Delta(W) in three different sample, with relaxed (i.e. zero) strain, 0.36% compressive and 0.48% tensile strain was measured as a function of temperature. The Delta(W) increases on cooling in NiO relaxed film as in the previously published work on a bulk crystal. The splitting increases on cooling also in 0.48% tensile strained film, but Delta(W) is systematically 3-4 cm-1 smaller than in relaxed film. Since the phonon splitting is proportional to the non-dominant magnetic exchange interaction J1, the reduction of phonon splitting in tensile-strained film was explained by a diminishing J1 with lattice expansion. Increase of Delta(W) on cooling can be also explained by rising of J1 with reduced temperature.


title: Realtime Index-Free Single Source SimRank Processing on Web-Scale Graphs

abstract: Given a graph G and a node u in G, a single source SimRank query evaluates the similarity between u and every node v in G. Existing approaches to single source SimRank computation incur either long query response time, or expensive pre-computation, which needs to be performed again whenever the graph G changes. Consequently, to our knowledge none of them is ideal for scenarios in which (i) query processing must be done in realtime, and (ii) the underlying graph G is massive, with frequent updates. Motivated by this, we propose SimPush, a novel algorithm that answers single source SimRank queries without any pre-computation, and at the same time achieves significantly higher query processing speed than even the fastest known index-based solutions. Further, SimPush provides rigorous result quality guarantees, and its high performance does not rely on any strong assumption of the underlying graph. Specifically, compared to existing methods, SimPush employs a radically different algorithmic design that focuses on (i) identifying a small number of nodes relevant to the query, and subsequently (ii) computing statistics and performing residue push from these nodes only. We prove the correctness of SimPush, analyze its time complexity, and compare its asymptotic performance with that of existing methods. Meanwhile, we evaluate the practical performance of SimPush through extensive experiments on 8 real datasets. The results demonstrate that SimPush consistently outperforms all existing solutions, often by over an order of magnitude. In particular, on a commodity machine, SimPush answers a single source SimRank query on a web graph containing over 133 million nodes and 5.4 billion edges in under 62 milliseconds, with 0.00035 empirical error, while the fastest index-based competitor needs 1.18 seconds.


title: Zero Order Stochastic Weakly Convex Composite Optimization

abstract: In this paper we consider stochastic weakly convex composite problems, however without the existence of a stochastic subgradient oracle. We present a derivative free algorithm that uses a two point approximation for computing a gradient estimate of the smoothed function. We prove convergence at a similar rate as state of the art methods, however with a larger constant, and report some numerical results showing the effectiveness of the approach.


title: Performance Analysis of Single-Cell Adaptive Data Rate-Enabled LoRaWAN

abstract: LoRaWAN enables massive connectivity for Internet-of-Things applications. Many published works employ stochastic geometry to derive outage models of LoRaWAN over fading channels assuming fixed transmit power and distance-based spreading factor (SF) allocation. However, in practice, LoRaWAN employs the Adaptive Data Rate (ADR) mechanism, which dynamically adjusts SF and transmit power of nodes based on channel state. The community addressed the performance of ADR using simulations, but analytical models have not been introduced. In this letter, we seek to close this gap. We build over an analytical LoRaWAN model to consider the performance of steady-state ADR-enabled LoRaWAN. We derive outage expressions and an optimization procedure to maximize the number of users under reliability constraints. Results show that power allocation reduces interference and improves network capacity while reducing average power.