title: Alma Twenty-six Arcmin^2 Survey Of Goods-s At One-millimeter (asagao): Source Catalog And Number Counts

abstract: We present the survey design, data reduction, construction of images, and source catalog of the Atacama Large Millimeter/submillimeter Array (ALMA) twenty-six arcmin^2 survey of GOODS-S at one-millimeter (ASAGAO). ASAGAO is a deep (1sigma ~ 61 uJy/beam for a 250 klambda-tapered map with a synthesized beam size of 0.51" x 0.45") and wide area (26 arcmin^2) survey on a contiguous field at 1.2 mm. By combining with ALMA archival data in the GOODS-South field, we obtained a deeper map in the same region (1sigma ~ 30 uJy/beam for a deep region with a 250 klambda-taper, and a synthesized beam size of 0.59" x 0.53"), providing the largest sample of sources (25 sources at >=5.0sigma, 45 sources at >=4.5sigma) among ALMA blank-field surveys to date. The number counts shows that 52(+11 -8)% of the extragalactic background light at 1.2 mm is resolved into discrete sources at S1.2m > 135 uJy. We create infrared (IR) luminosity functions (LFs) in the redshift range of z = 1-3 from the ASAGAO sources with KS-band counterparts, and constrain the faintest luminosity of the LF at 2.0 < z < 3.0. The LFs are consistent with previous results based on other ALMA and SCUBA-2 observations, which suggest a positive luminosity evolution and negative density evolution with increasing redshift. We find that obscured star-formation of sources with IR luminosities of log(L(IR)/Lsun)} ~> 11.8 account for ~~60%-90% of the z ~ 2 cosmic star-formation rate density.


title: Shared Multi-Task Imitation Learning for Indoor Self-Navigation

abstract: Deep imitation learning enables robots to learn from expert demonstrations to perform tasks such as lane following or obstacle avoidance. However, in the traditional imitation learning framework, one model only learns one task, and thus it lacks of the capability to support a robot to perform various different navigation tasks with one model in indoor environments. This paper proposes a new framework, Shared Multi-headed Imitation Learning(SMIL), that allows a robot to perform multiple tasks with one model without switching among different models. We model each task as a sub-policy and design a multi-headed policy to learn the shared information among related tasks by summing up activations from all sub-policies. Compared to single or non-shared multi-headed policies, this framework is able to leverage correlated information among tasks to increase performance.We have implemented this framework using a robot based on NVIDIA TX2 and performed extensive experiments in indoor environments with different baseline solutions. The results demonstrate that SMIL has doubled the performance over nonshared multi-headed policy.


title: Searching for Biosignatures in Exoplanetary Impact Ejecta

abstract: With the number of confirmed rocky exoplanets increasing steadily, their characterisation and the search for exoplanetary biospheres is becoming an increasingly urgent issue in astrobiology. We aim to investigate the possibility of characterising an exoplanet (in terms of habitability, geology, presence of life etc.) by studying material ejected from the surface during an impact event. For given parameters characterising the impact event, we estimate the escaping mass and assess its subsequent collisional evolution in a circumstellar orbit, assuming a Sun-like host star. We calculate the fractional luminosity of the dust as a function of time after the impact event and study its detectability with current and future instrumentation. We consider the possibility to constrain the dust composition, giving information on the geology or the presence of a biosphere. As examples, we investigate whether calcite, silica or ejected microorganisms could be detected. For a 20 km diameter impactor, we find that the dust mass escaping the exoplanet is roughly comparable to the zodiacal dust. The collisional evolution is best modelled by considering two independent dust populations, a spalled population consisting of non-melted ejecta evolving on timescales of millions of years, and dust recondensed from melt or vapour evolving on much shorter timescales. While the presence of dust can potentially be inferred with current telescopes, studying its composition requires advanced instrumentation not yet available. The direct detection of biological matter turns out to be extremely challenging. Despite considerable difficulties (small dust masses, noise such as exozodiacal dust etc.), studying dusty material ejected from an exoplanetary surface might become an interesting complement to atmospheric studies in the future.


title: Fine-Grained Representation Learning and Recognition by Exploiting Hierarchical Semantic Embedding

abstract: Object categories inherently form a hierarchy with different levels of concept abstraction, especially for fine-grained categories. For example, birds (Aves) can be categorized according to a four-level hierarchy of order, family, genus, and species. This hierarchy encodes rich correlations among various categories across different levels, which can effectively regularize the semantic space and thus make prediction less ambiguous. However, previous studies of fine-grained image recognition primarily focus on categories of one certain level and usually overlook this correlation information. In this work, we investigate simultaneously predicting categories of different levels in the hierarchy and integrating this structured correlation information into the deep neural network by developing a novel Hierarchical Semantic Embedding (HSE) framework. Specifically, the HSE framework sequentially predicts the category score vector of each level in the hierarchy, from highest to lowest. At each level, it incorporates the predicted score vector of the higher level as prior knowledge to learn finer-grained feature representation. During training, the predicted score vector of the higher level is also employed to regularize label prediction by using it as soft targets of corresponding sub-categories. To evaluate the proposed framework, we organize the 200 bird species of the Caltech-UCSD birds dataset with the four-level category hierarchy and construct a large-scale butterfly dataset that also covers four level categories. Extensive experiments on these two and the newly-released VegFru datasets demonstrate the superiority of our HSE framework over the baseline methods and existing competitors.


title: Thermodynamics and Phase Transitions of Nonlinear Electrodynamics Black Holes in an Extended Phase Space

abstract: We investigate the thermodynamic behavior of nonlinear electrodynamics (NLED) black holes in an extended phase space, which includes the cosmological constant $\Lambda=-3/l^{2}$ and dimensionful couplings $a$ in NLED as thermodynamic variables. For a generic NLED black hole with the charge $Q$, we find that the Smarr relation is satisfied in the extended phase space, and the state of equation can be written as $Tl=\tilde{T}\left( r_{+}/l,Q/l,al^{-c} \right) $, where $\left[ a\right] =L^{c}$, and $T$ and $r_{+}$ are the temperature and horizon radius of the black hole, respectively. For some values of $Q/l$ and $al^{-c}$, the phase structure of the black hole is uniquely determined. Focusing on Born-Infeld and iBorn-Infeld AdS black holes, we obtain the corresponding phase diagrams in the $a/l^{2}$-$Q/l$ plane, which provides a new viewpoint towards the black holes' phase structure and critical behavior. For Born-Infeld black holes, the critical line and the region, where a reentrant phase transition occurs, in the $a/l^{2}$-$Q/l$ plane are both finite and terminate at $\left\{ \tilde{a}_{c}\text{, }\tilde{Q}_{c}\right\} \simeq\left\{ 0.069\text{, }0.37\right\} $. However for iBorn-Infeld black holes, the critical line and the reentrant phase transition region in the $a/l^{2}$-$Q/l$ plane are semi-infinite and extend to $Q/l=\infty$. We also examine thermal and electrical stabilities of Born-Infeld and iBorn-Infeld AdS black holes.


title: Alternating Iterative Secure Structure between Beamforming and Power Allocation for UAV-aided Directional Modulation Networks

abstract: In unmanned aerial vehicle (UAV) networks, directional modulation (DM) is adopted to improve the secrecy rate (SR) performance. Alice, a ground base station, behaves as a control center, and Bob is a UAV of flying along a linear flight trajectory who optimizes its SR performance by dynamically adjusting its beamforming vectors and power allocation (PA) strategy. Per fixed time interval during the Bob's flight process, the transmit beamforming vectors for useful messages and AN projection are given by the rule of maximizing signal-to-leakage-and-noise ratio (Max-SLNR) and maximizing AN-and-leakage-to-noise ratio (ANLNR), and the optimal PA strategy is based on maximizing SR (Max-SR). More importantly, an alternating iterative structure (AIS) between beamforming and PA is proposed to further improve the SR performance. Simulation results show that the proposed AIS converges rapidly, and can achieve substantial SR gains over Max-SLNR plus Max-ANLNR with fixed PA such as PA factor $\beta=$ 0.5, and 0.9. In particular, in the case of small-scale antenna array, the SR performance gain achieved by the proposed AIS is more attractive. Additionally, as the number of antennas tends to be large-scale, the average SR performance of the proposed AIS approaches an SR ceil.


title: Internal Kinematics of UM 461 and CTS 1020

abstract: We have used integral field spectroscopy to study the internal kinematics of the H\,{\sc ii} galaxies CTS 1020 and UM 461. We based our analysis on the velocity and velocity dispersion maps, and intensity-velocity dispersion ($I-\sigma$) and velocity-velocity dispersion ($V_r-\sigma$) diagrams. We found that the motion in both star-forming knots of UM 461 has different patterns, suggesting a weak kinematical connection between the knots. The overall kinematics of the galaxy is probably affected by stellar feedback. CTS 1020 has an ordered motion with a gradient compatible with a disc rotating at $\sim50$ km s$^{-1}$, though the velocity field is disturbed. In both galaxies the highest and lowest $\sigma$ values are distributed in the outer parts and are associated with the diffuse gas that permeates the galaxies. UM 461 has a ring-like structure with small regions of increasing $\sigma$ in the eastern knot, which resemble what we could expect in a collect and collapse scenario of star formation. We found that UM 461 seems to be more susceptible to stellar feedback, whereas in CTS 1020 the gravitational potential dominates.


title: Concentration-adjustable micromixer using droplet injection into a microchannel

abstract: A novel micromixing technique that exploit a thrust of droplets into the mixing interface is developed. The technique enhances the mixing by injecting immiscible droplets in a mixing channel and the methodology enables a control of the mixing level simply by changing the droplet injection frequency. We experimentally characterize the mixing performance with various droplet injection frequencies, channel geometries, and diffusion coefficients. Consequently, it is revealed that the mixing level increases with the injection frequency, the droplet-diameter-to-channel-width ratio, and the diffusion coefficient. Moreover, the mixing level is found to be a linear function of the droplet volume fraction in the mixing section. The results suggest that the developed technique can produce a large amount of sample solution whose concentration is arbitrary and precisely controllable with a simple and stable operation.


title: On the approximability of the stable matching problem with ties of size two

abstract: The stable matching problem is one of the central problems of algorithmic game theory. If participants are allowed to have ties, the problem of finding a stable matching of maximum cardinality is an NP-hard problem, even when the ties are of size two. Moreover, in this setting it is UGC-hard to provide an approximation for the maximum cardinality stable matching problem with a constant factor smaller than 4/3. In this paper, we give a tight analysis of an approximation algorithm given by Huang and Kavitha for the maximum cardinality stable matching problem with ties of size two, demonstrating an improved 4/3-approximation factor.


title: A Record Linkage Model Incorporating Relational Data

abstract: In this paper we introduce a novel Bayesian approach for linking multiple social networks in order to discover the same real world person having different accounts across networks. In particular, we develop a latent model that allow us to jointly characterize the network and linkage structures relying in both relational and profile data. In contrast to other existing approaches in the machine learning literature, our Bayesian implementation naturally provides uncertainty quantification via posterior probabilities for the linkage structure itself or any function of it. Our findings clearly suggest that our methodology can produce accurate point estimates of the linkage structure even in the absence of profile information, and also, in an identity resolution setting, our results confirm that including relational data into the matching process improves the linkage accuracy. We illustrate our methodology using real data from popular social networks such as Twitter, Facebook, and YouTube.


title: Multicast Triangular Semilattice Network

abstract: We investigate the structure of the code graph of a multicast network that has a characteristic shape of an inverted equilateral triangle. We provide a criterion that determines the validity of a receiver placement within the code graph, present invariance properties of the determinants corresponding to receiver placements under symmetries, and provide a complete study of these networks' receivers and required field sizes up to a network of 4 sources. We also improve on various definitions related to code graphs.


title: Ridge Rerandomization: An Experimental Design Strategy in the Presence of Collinearity

abstract: Randomization ensures that observed and unobserved covariates are balanced, on average. However, randomizing units to treatment and control often leads to covariate imbalances in realization, and such imbalances can inflate the variance of estimators of the treatment effect. One solution to this problem is rerandomization---an experimental design strategy that randomizes units until some balance criterion is fulfilled---which yields more precise estimators of the treatment effect if covariates are correlated with the outcome. Most rerandomization schemes in the literature utilize the Mahalanobis distance, which may not be preferable when covariates are correlated or vary in importance. As an alternative, we introduce an experimental design strategy called ridge rerandomization, which utilizes a modified Mahalanobis distance that addresses collinearities among covariates and automatically places a hierarchy of importance on the covariates according to their eigenstructure. This modified Mahalanobis distance has connections to principal components and the Euclidean distance, and---to our knowledge---has remained unexplored. We establish several theoretical properties of this modified Mahalanobis distance and our ridge rerandomization scheme. These results guarantee that ridge rerandomization is preferable over randomization and suggest when ridge rerandomization is preferable over standard rerandomization schemes. We also provide simulation evidence that suggests that ridge rerandomization is particularly preferable over typical rerandomization schemes in high-dimensional or high-collinearity settings.


title: Binary Black Hole Mergers from Globular Clusters: the Impact of Globular Cluster Properties

abstract: The dense environment of globular clusters (GCs) can facilitate the formation of binary black holes (BBHs), some of which can merge with gravitational waves (GW) within the age of the Universe. We have performed a survey of Monte-Carlo simulations following the dynamical evolution of GCs with different masses, sizes and binary fractions and explored the impact of the host GC properties on the formation of BBH mergers. We find that the number of BBH mergers from GCs is determined by the GC's initial mass, size and primordial binary fraction. We identify two groups of BBH mergers: a primordial group whose formation does not depend on cluster's dynamics and a dynamical group whose formation is driven by the cluster's dynamical evolution. We show how the BBH origin affects the BBH mergers' main properties such as the chirp mass and merging time distributions. We provide analytic expressions for the dependence of the number of BBH mergers from individual GCs on the main cluster's structural properties and the time evolution of the merger rates of these BBHs. These expressions provide an essential ingredient for a general framework allowing to estimate the merger rate density. Using the relations found in our study, we find a local merger rate density of 0.18$-$1.8 ${\rm Gpc}^{-3}{\rm yr}^{-1}$ for primordial BBH mergers and 0.6$-$18 ${\rm Gpc}^{-3}{\rm yr}^{-1}$ for dynamical BBH mergers, depending on the GC mass and size distributions, initial binary fraction and the number density of GCs in the Universe.


title: Relaxation algorithms for matrix completion, with applications to seismic travel-time data interpolation

abstract: Travel time tomography is used to infer the underlying three-dimensional wavespeed structure of the Earth by fitting seismic travel time data collected at surface stations. Data interpolation and denoising techniques are important pre-processing steps that use prior knowledge about the data, including parsimony in the frequency and wavelet domains, low-rank structure of matricizations, and local smoothness. We show how local smoothness structure can be combined with low rank constraints using level-set optimization formulations, and develop a new relaxation algorithm that can efficiently solve these joint problems. In the seismology setting, we use the approach to interpolate missing stations and de-noise observed stations. The new approach is competitive with alternative algorithms, and offers new functionality to interpolate observed data using both smoothness and low rank structure in the presence of data fitting constraints.


title: In situ measurement of MWA primary beam variation using ORBCOMM

abstract: We provide the first in situ measurements of antenna element (tile) beam shapes of the Murchison Widefield Array (MWA), a low radio-frequency interferometer and an SKA precursor. Most current MWA processing pipelines use an assumed beam shape, errors in which can cause absolute and relative flux density errors, as well as polarisation 'leakage'. This makes understanding the primary beam of paramount importance, especially for sensitive experiments such as a measurement of the 21 cm line from the epoch of reionisation (EoR). The calibration requirements for measuring the EoR 21 cm line are so extreme that tile to tile beam variations may affect our ability to make a detection. Measuring the primary beam shape from visibilities alone is challenging, as multiple instrumental, atmospheric, and astrophysical factors contribute to uncertainties in the data. Building on the methods of Neben et al. (2015), we tap directly into the receiving elements of the MWA before any digitisation or correlation of the signal. Using ORBCOMM satellite passes we are able to produce all-sky maps for 4 separate tiles in the XX polarisation. We find good agreement with the cutting-edge 'fully' embedded element (FEE) model of Sokolowski et al. (2017), and observe that the MWA beamformers consistently recreate beam shapes to within ~1dB in the reliable areas of our beam maps. We also clearly observe the effects of a missing dipole from a tile in one of our beam maps, and show that the FEE model is able to reproduce this modified beam shape. We end by motivating and outlining additional onsite experiments to further constrain the primary beam behaviour.


title: Feasibility of 5G mm-wave communication for connected autonomous vehicles

abstract: The internet-of-things (IoT) environment holds different intelligent components networked together and will enable seamless data communication between the connected components. Connected autonomous vehicles or CAVs are major components of the IoT, and the smooth, reliable, and safe operation of CAVs demands a reliable wireless communication system, which can ensure high connectivity, high throughput and low communication latency. The 5G millimeter-wave or mm-wave communication network offers such benefits, which can be the enabler of CAVs, especially for dense congested areas. In this research, we have evaluated the 5G mm-wave and Dedicated Short Range Communication (DSRC) for different CAV applications in a Network Simulator-3 (ns-3). For the CAV applications, we have evaluated the end-to-end latency, packet loss, and data rate (for both the CAV receiver and transmitter) of the 5G mm-wave. We have found that 5G mm-wave can support CAV safety applications by ensuring lower latency compared to the required minimum latency of 200 milliseconds for the forward collision warning application. For mobility applications, we have found that 5G mm-wave can support multiple CAVs with a high data receive rate, which is enough for real-time high definition video streaming for in-vehicle infotainment, with mean packet delay of 13 milliseconds. The findings from this study establish that 5G mm-wave can be the enabler of future CAVs in congested areas. Using the evaluation framework developed in this paper, public agencies can evaluate 5G mm-wave to support CAVs in congested areas, such as heavy pedestrian areas like downtown, commercial zones, under their jurisdiction.


title: System identification calorimetry

abstract: We report a model-based method for quantifying heat flow and storage in thermal systems using data from multiple thermal sensors. This approach avoids stringent requirements on the system geometry and sensor positions and enables calorimetry to be performed under a broader range of circumstances than is accessible with existing calorimeters, such as when non-linear heat transfer occurs, when spatially separated heat sources are active, or when multiple thermal masses participate. Using experimental data from a model thermal system, this paper provides a tutorial on the construction of non-linear lumped element heat transfer models and the use of system identification to estimate the parameters of these models from calibration data. The calibrated models are then used to estimate unknown energy inputs to the thermal system from sensor data. Our best model enabled the measurement of the total input energy with 0.02% accuracy; the instantaneous input power could be measured with a root-mean-square error of 10% of the average input power.


title: A machine-learning solver for modified diffusion equations

abstract: A feed-forward neural network has a remarkable property which allows the network itself to be a universal approximator for any functions.Here we present a universal, machine-learning based solver for multi-variable partial differential equations. The algorithm approximates the target functions by neural networks and adjusts the network parameters to approach the desirable solutions.The idea can be easily adopted for dealing with multi-variable, coupled integrodifferential equations, such as those in the self-consistent field theory for predicting polymer microphase-separated structures.


title: On the level of modular curves that give rise to isolated $j$-invariants

abstract: We say a closed point $x$ on a curve $C$ is sporadic if $C$ has only finitely many closed points of degree at most $\operatorname{deg}(x)$ and that $x$ is isolated if it is not in a family of effective degree $d$ divisors parametrized by $\mathbb{P}^1$ or a positive rank abelian variety (see Section 4 for more precise definitions and a proof that sporadic points are isolated). Motivated by well-known classification problems concerning rational torsion of elliptic curves, we study sporadic and isolated points on the modular curves $X_1(N)$. In particular, we show that any non-cuspidal non-CM sporadic, respectively isolated, point $x \in X_1(N)$ maps down to a sporadic, respectively isolated, point on a modular curve $X_1(d)$, where $d$ is bounded by a constant depending only on $j(x)$. Conditionally, we show that $d$ is bounded by a constant depending only on the degree of $\mathbb{Q}(j(x))$, so in particular there are only finitely many $j$-invariants of bounded degree that give rise to sporadic or isolated points.


title: Low Rank Regularization: A Review

abstract: Low rank regularization, in essence, involves introducing a low rank or approximately low rank assumption for matrix we aim to learn, which has achieved great success in many fields including machine learning, data mining and computer version. Over the last decade, much progress has been made in theories and practical applications. Nevertheless, the intersection between them is very slight. In order to construct a bridge between practical applications and theoretical research, in this paper we provide a comprehensive survey for low rank regularization. We first review several traditional machine learning models using low rank regularization, and then show their (or their variants) applications in solving practical issues, such as non-rigid structure from motion and image denoising. Subsequently, we summarize the regularizers and optimization methods that achieve great success in traditional machine learning tasks but are rarely seen in solving practical issues. Finally, we provide a discussion and comparison for some representative regularizers including convex and non-convex relaxations. Extensive experimental results demonstrate that non-convex regularizers can provide a large advantage over the nuclear norm, the regularizer widely used in solving practical issues.


title: Hydra games for recursively Mahlo operations

abstract: A hydra game is proposed, and the fact that every hydra eventually die out is shown to be equivalent (over a weak arithmetic) to the 1-consistency of set theory KPM for recursively Mahlo universes.


title: Adaptive Sampling for Convex Regression

abstract: In this paper, we introduce the first principled adaptive-sampling procedure for learning a convex function in the $L_\infty$ norm, a problem that arises often in the behavioral and social sciences. We present a function-specific measure of complexity and use it to prove that, for each convex function $f_{\star}$, our algorithm nearly attains the information-theoretically optimal, function-specific error rate. We also corroborate our theoretical contributions with numerical experiments, finding that our method substantially outperforms passive, uniform sampling for favorable synthetic and data-derived functions in low-noise settings with large sampling budgets. Our results also suggest an idealized "oracle strategy", which we use to gauge the potential advance of any adaptive-sampling strategy over passive sampling, for any given convex function.


title: Darboux evaluations for hypergeometric functions with the projective monodromy PSL(2,F7)

abstract: Algebraic hypergeometric functions can be compactly expressed as radical functions on pull-back curves where the monodromy group is simpler, say, a finite cyclic group. These so-called Darboux evaluations were already considered for algebraic 2F1-functions. This article presents Darboux evaluations for the classical case of 3F2-functions with the projective monodromy group PSL(2,F7). As an application, appealing modular evaluations of the same 3F2-functions are derived.


title: Discrete Structural Planning for Neural Machine Translation

abstract: Structural planning is important for producing long sentences, which is a missing part in current language generation models. In this work, we add a planning phase in neural machine translation to control the coarse structure of output sentences. The model first generates some planner codes, then predicts real output words conditioned on them. The codes are learned to capture the coarse structure of the target sentence. In order to obtain the codes, we design an end-to-end neural network with a discretization bottleneck, which predicts the simplified part-of-speech tags of target sentences. Experiments show that the translation performance are generally improved by planning ahead. We also find that translations with different structures can be obtained by manipulating the planner codes.


title: A Density Functional Tight Binding Layer for Deep Learning of Chemical Hamiltonians

abstract: Current neural networks for predictions of molecular properties use quantum chemistry only as a source of training data. This paper explores models that use quantum chemistry as an integral part of the prediction process. This is done by implementing self-consistent-charge Density-Functional-Tight-Binding (DFTB) theory as a layer for use in deep learning models. The DFTB layer takes, as input, Hamiltonian matrix elements generated from earlier layers and produces, as output, electronic properties from self-consistent field solutions of the corresponding DFTB Hamiltonian. Backpropagation enables efficient training of the model to target electronic properties. Two types of input to the DFTB layer are explored, splines and feed-forward neural networks. Because overfitting can cause models trained on smaller molecules to perform poorly on larger molecules, regularizations are applied that penalize non-monotonic behavior and deviation of the Hamiltonian matrix elements from those of the published DFTB model used to initialize the model. The approach is evaluated on 15,700 hydrocarbons by comparing the root mean square error in energy and dipole moment, on test molecules with 8 heavy atoms, to the error from the initial DFTB model. When trained on molecules with up to 7 heavy atoms, the spline model reduces the test error in energy by 60% and in dipole moments by 42%. The neural network model performs somewhat better, with error reductions of 67% and 59% respectively. Training on molecules with up to 4 heavy atoms reduces performance, with both the spline and neural net models reducing the test error in energy by about 53% and in dipole by about 25%.


title: Weight Learning in a Probabilistic Extension of Answer Set Programs

abstract: LPMLN is a probabilistic extension of answer set programs with the weight scheme derived from that of Markov Logic. Previous work has shown how inference in LPMLN can be achieved. In this paper, we present the concept of weight learning in LPMLN and learning algorithms for LPMLN derived from those for Markov Logic. We also present a prototype implementation that uses answer set solvers for learning as well as some example domains that illustrate distinct features of LPMLN learning. Learning in LPMLN is in accordance with the stable model semantics, thereby it learns parameters for probabilistic extensions of knowledge-rich domains where answer set programming has shown to be useful but limited to the deterministic case, such as reachability analysis and reasoning about actions in dynamic domains. We also apply the method to learn the parameters for probabilistic abductive reasoning about actions.


title: General conversion method for constrained systems

abstract: We reformulate in a systematic way the conversional approach in its most general and compact form. We present a new definition of generalized Dirac bracket directly in terms of the super-observables commuting with the basic BFV-BRST charge.


title: Simulating topological tensor networks with Majorana qubits

abstract: The realization of topological quantum phases of matter remains a key challenge to condensed matter physics and quantum information science. In this work, we demonstrate that progress in this direction can be made by combining concepts of tensor network theory with Majorana device technology. Considering the topological double semion string-net phase as an example, we exploit the fact that the representation of topological phases by tensor networks can be significantly simpler than their description by lattice Hamiltonians. The building blocks defining the tensor network are tailored to realization via simple units of capacitively coupled Majorana bound states. In the case under consideration, this yields a remarkably simple blueprint of a synthetic double semion string-net, and one may be optimistic that the required device technology will be available soon. Our results indicate that the implementation of tensor network structures via mesoscopic quantum devices opens up a powerful novel avenue toward the realization and quantum simulation of synthetic topological quantum matter.


title: Hybrid Powerline/Wireless Diversity for Smart Grid Communications: Design Challenges and Real-time Implementation

abstract: The demand for energy is growing at an unprecedented pace that is much higher than the energy generation capacity growth rate using both conventional and green technologies.In particular, the electric power sector is consistently rated among the most dynamic growth markets over all other energy markets. Distributed (decentralized) energy generation based on renewable energy sources is an efficient and reliable solution to serve such huge energy demand growth [1]. However, to manage dynamic and complex distributed systems, a massive amount of data has to be measured, collected, exchanged and processed in real time. Smart grids manage an intelligent energy delivery network enabled two-way communications between data concentrators operated by utility companies and smart meters installed at the end users. In particular, dynamic power-grid loading and peak load management are the two main driving forces for bidirectional communications over the grid. Narrowband power line communications (NB-PLC) and wireless communications in the unlicensed frequency band (sub-1 GHz or 2.4 GHz) are the two main communications systems adopted to support the growing smart grid applications. Moreover, since NB-PLC and unlicensed wireless links experience channel and interference with markedly different statistics, transmitting the same information signal concurrently over both links significantly enhances the smart grid communications reliability. In this article, we compare various diversity combining schemes for simultaneous power line and wireless transmissions. Furthermore, we developed a real-time testbed for the hybrid PLC/wireless system to demonstrate the performance enhancement achieved by PLC/wireless diversity combining over a single link performance.


title: Black hole demographics from TDE modeling

abstract: The occurrence rate of tidal disruption events (TDEs) by survey missions depend on the black hole mass function of the galaxies, properties of the stellar cusp and mass of the central black hole. Using a power law density profile with Kroupa mass function, we solve the steady state Fokker-Planck to calculate the theoretical capture rate of stars by the black hole. Using a steady accretion model, the Schechter black hole mass function (BHMF) and the cosmological parameters, we calculate the detection rate of TDEs for various surveys which is then fit with the observed TDE rates to extract the Schechter parameters. The rate tension between observation ($\sim 10^{-5}~{\rm yr^{-1}}$) and theory ($\sim 10^{-4}~{\rm yr^{-1}}$) for individual galaxies is explained by the statistical average over the BHMF.


title: An Extension of Waring's Problem in $\mathbb{Z}_{p^k}$

abstract: In this paper, we will investigate the solvability of the equation $x_1^k + x_2^k + \ldots + x_s^k = n$, $n\in \mathbb{Z}_{p^k}$, $x_1,...,x_s\in \mathcal{A}$, $\mathcal{A}\subseteq \mathbb{Z}_{p^k}$. We will give a upper bound of the number of solutions for $s=2$.


title: Revised Exoplanet Radii and Habitability Using Gaia Data Release 2

abstract: Accurate stellar properties are crucial for determining exoplanet characteristics. Gaia DR2 presents revised distances, luminosities, and radii for 1.6 billion stars. Here, we report the calculation of revised radii and densities for 320 exoplanets using this data and present updated calculations of the incident flux received by 690 known exoplanets. This allows the likelihood that those planets orbit in the habitable zone of their host stars to be reassessed. As a result of this analysis, three planets can be added to the catalogue of potentially habitable worlds: HIP~67537~b, HD~148156~b, and HD~106720~b. In addition, the changed parameterisation of BD~+49~898 means that its planet, BD~+49~898~b, now receives an incident flux that places it outside the optimistic habitable zone region, as defined by \citep{Kopparapu2013,Kopparapu2014}. We find that use of the new \textit{Gaia} data results in a mean increase in calculated exoplanet radius of 3.76\%. Previously, CoRoT-3 b had been reported as having the highest density of all known exoplanets. Here, we use updated information to revise the calculated density of CoRoT-3~b from 26.4$g\:cm^{-3}$ to 17.3$\pm2.9g\:cm^{-3}$. We also report the densest exoplanet in our dataset, KELT-1~b, with a density of 23.7$\pm4.0g\:cm^{-3}$. Overall, our results highlight the importance of ensuring the the parameterisation of known exoplanets be revisited whenever significant improvements are made to the precision of the stellar parameters upon which they are based.


title: Stable almost complex structures on certain $10$-manifolds

abstract: Let $M$ be a $10$-dimensional closed oriented smooth manifold. Set $$\mathcal{D}_{M} := \{ x \in H^{2}(M; \Z/2) \mid x^{2} + w_{2}(M) x \in \rho_{2} ( TH^{4}(M;\Z) ) \}.$$ Suppose that $H_{1}(M;\Z)=0$ and $\mathcal{D}_{M} \subset \rho_{2}( H^{2}(M; \Z) )$. Then the necessary and sufficient conditions for $M$ to admit a stable almost complex structure are determined in terms of the characteristic classes and cohomology ring of $M$.


title: Microscopic optical potentials for calcium isotopes

abstract: We construct nucleonic microscopic optical potentials by combining the Green's function approach with the coupled-cluster method for $\rm{^{40}Ca}$ and $\rm{^{48}Ca}$. For the computation of the ground-state of $\rm{^{40}Ca}$ and $\rm{^{48}Ca}$, we use the coupled-cluster method in the singles-and-doubles approximation, while for the A = $\pm 1$ nuclei we use particle-attached/removed equation-of-motion method truncated at two-particle-one-hole and one-particle-two-hole excitations, respectively. Our calculations are based on the chiral nucleon-nucleon and three-nucleon interaction $\rm{NNLO_{sat}}$, which reproduces the charge radii of $^{40}$Ca and $^{48}$Ca, and the chiral nucleon-nucleon interaction $\rm{NNLO_{opt}}$. In all cases considered here, we observe that the overall form of the neutron scattering cross section is reproduced for both interactions, but the imaginary part of the potential, which reflects the loss of flux in the elastic channel, is negligible. The latter points to neglected many-body correlations that would appear beyond the coupled-cluster truncation level considered in this work. We show that, by artificially increasing the parameter $\eta$ in the Green's function, practical results can be further improved.


title: Evolution of the $M_{\bullet} - \sigma$ relation

abstract: Black holes at the centers of the galaxies grow mainly by the processes of accretion, mergers, and consumption of stars. In the case of gas accretion with cooling sources, the flow is momentum driven, after which the black hole reaches a saturated mass, and subsequently, it grows only by consumption of stars. In addition, we include the effect of mergers on the growth of black hole spin and mass and study its evolution as a function of redshift in a $\Lambda$CDM cosmology using an initial seed mass and spin distribution functions that we have derived. For the stellar ingestion, we have assumed a power-law density profile for the galaxy in our framework of a new relativistic loss cone theory that includes the effect of the black hole spin. We predict the impact of the evolution on the $M_{\bullet} - \sigma$ relation and compare it with available observations.


title: Learning Linear Transformations for Fast Arbitrary Style Transfer

abstract: Given a random pair of images, an arbitrary style transfer method extracts the feel from the reference image to synthesize an output based on the look of the other content image. Recent arbitrary style transfer methods transfer second order statistics from reference image onto content image via a multiplication between content image features and a transformation matrix, which is computed from features with a pre-determined algorithm. These algorithms either require computationally expensive operations, or fail to model the feature covariance and produce artifacts in synthesized images. Generalized from these methods, in this work, we derive the form of transformation matrix theoretically and present an arbitrary style transfer approach that learns the transformation matrix with a feed-forward network. Our algorithm is highly efficient yet allows a flexible combination of multi-level styles while preserving content affinity during style transfer process. We demonstrate the effectiveness of our approach on four tasks: artistic style transfer, video and photo-realistic style transfer as well as domain adaptation, including comparisons with the state-of-the-art methods.


title: Text-to-Image-to-Text Translation using Cycle Consistent Adversarial Networks

abstract: Text-to-Image translation has been an active area of research in the recent past. The ability for a network to learn the meaning of a sentence and generate an accurate image that depicts the sentence shows ability of the model to think more like humans. Popular methods on text to image translation make use of Generative Adversarial Networks (GANs) to generate high quality images based on text input, but the generated images don't always reflect the meaning of the sentence given to the model as input. We address this issue by using a captioning network to caption on generated images and exploit the distance between ground truth captions and generated captions to improve the network further. We show extensive comparisons between our method and existing methods.


title: Constructions of maximally recoverable local reconstruction codes via function fields

abstract: Local Reconstruction Codes (LRCs) allow for recovery from a small number of erasures in a local manner based on just a few other codeword symbols. A maximally recoverable (MR) LRC offers the best possible blend of such local and global fault tolerance, guaranteeing recovery from all erasure patterns which are information-theoretically correctable given the presence of local recovery groups. In an $(n,r,h,a)$-LRC, the $n$ codeword symbols are partitioned into $r$ disjoint groups each of which include $a$ local parity checks capable of locally correcting $a$ erasures. MR LRCs have received much attention recently, with many explicit constructions covering different regimes of parameters. Unfortunately, all known constructions require a large field size that exponential in $h$ or $a$, and it is of interest to obtain MR LRCs of minimal possible field size. In this work, we develop an approach based on function fields to construct MR LRCs. Our method recovers, and in most parameter regimes improves, the field size of previous approaches. For instance, for the case of small $r \ll \epsilon \log n$ and large $h \ge \Omega(n^{1-\epsilon})$, we improve the field size from roughly $n^h$ to $n^{\epsilon h}$. For the case of $a=1$ (one local parity check), we improve the field size quadratically from $r^{h(h+1)}$ to $r^{h \lfloor (h+1)/2 \rfloor}$ for some range of $r$. The improvements are modest, but more importantly are obtained in a unified manner via a promising new idea.


title: A Ramsey-type theorem for the matching number regarding connected graphs

abstract: A major line of research is discovering Ramsey-type theorems, which are results of the following form: given a graph parameter $\rho$, every graph $G$ with sufficiently large $\rho(G)$ contains a `well-structured' induced subgraph $H$ with large $\rho(H)$. The classical Ramsey's theorem deals with the case when the graph parameter under consideration is the number of vertices; there is also a Ramsey-type theorem regarding connected graphs. Given a graph $G$, the matching number and the induced matching number of $G$ is the maximum size of a matching and an induced matching, respectively, of $G$. In this paper, we formulate Ramsey-type theorems for the matching number and the induced matching number regarding connected graphs. Along the way, we obtain a Ramsey-type theorem for the independence number regarding connected graphs as well.


title: Eruptions from quiet Sun coronal bright points. I. Observations

abstract: Observations of the full lifetime of CBPs in data taken with the AIA on board SDO in four passbands, He II 304 A, Fe IX/X 171 A, Fe XII 193 A, and Fe XVIII 94 A are investigated for the occurrence of plasma ejections, micro-flaring, mini-filament eruptions and mini coronal mass ejections (mini-CMEs). First and foremost, our study shows that the majority (76%) of quiet Sun CBPs (31 out of 42 CBPs) produce at least one eruption during their lifetime. From 21 eruptions in 11 CBPs, 18 occur in average ~17 hrs after the CBP formation for an average lifetime of the CBPs in AIA 193 A of ~21 hrs. This time delay in the eruption occurrence coincides in each BP with the convergence and cancellation phase of the CBP bipole evolution during which the CBPs become smaller until they fully disappear. The remaining three happen 4 - 6 hrs after the CBP formation. In sixteen out of 21 eruptions the magnetic convergence and cancellation involve the CBP main bipoles, while in three eruptions one of the BP magnetic fragments and a pre-existing fragment of opposite polarity converge and cancel. In one BP with two eruptions cancellation was not observed. The CBP eruptions involve in most cases the expulsion of chromospheric material either as elongated filamentary structure (mini-filament, MF) or a volume of cool material (cool plasma cloud, CPC), together with the CBP or higher overlying hot loops. Coronal waves were identified during three eruptions. A micro-flaring is observed beneath all erupting MFs/CPCs. It remains uncertain whether the destabilised MF causes the micro-flaring or the destabilisation and eruption of the MF is triggered by reconnection beneath the filament. In most eruptions, the cool erupting plasma obscures partially or fully the micro-flare until the erupting material moves away from the CBP. From 21 eruptions 11 are found to produce mini-CMEs.


title: The Discrete-Dual Minimal-Residual Method (DDMRes) for Weak Advection-Reaction Problems in Banach spaces

abstract: We propose and analyse a minimal-residual method in discrete dual norms for approximating the solution of the advection-reaction equation in a weak Banach-space setting. The weak formulation allows for the direct approximation of solutions in the Lebesgue $L^p$-space, $1<p<\infty$. The greater generality of this weak setting is natural when dealing with rough data and highly irregular solutions, and when enhanced qualitative features of the approximations are needed. We first present a rigorous analysis of the well-posedness of the underlying continuous weak formulation, under natural assumptions on the advection-reaction coefficients. The main contribution is the study of several discrete subspace pairs guaranteeing the discrete stability of the method and quasi-optimality in $L^p$, and providing numerical illustrations of these findings, including the elimination of Gibbs phenomena, computation of optimal test spaces, and application to 2-D advection.


title: Weak coupling theory of topological Hall effect

abstract: Topological Hall effect (THE) caused by a noncoplanar spin texture characterized by a scalar spin chirality is often described by the Berry phase, or the associated effective magnetic field. This picture is appropriate when the coupling, $M$, of conduction electrons to the spin texture is strong (strong-coupling regime) and the adiabatic condition is satisfied. However, in the weak-coupling regime, where the coupling $M$ is smaller than the electrons' scattering rate, the adiabatic condition is not satisfied and the Berry phase picture does not hold. In such a regime, the relation of the effective magnetic field to the spin texture can be "nonlocal", in contrast to the "local" relation in the strong coupling case. Focusing on the case of continuous but general spin texture, we investigate the THE in various characteristic regions in the weak-coupling regime, namely, (1) diffusive and local, (2) diffusive and nonlocal, and (3) ballistic. In the presence of spin relaxation, there arise two more regions in the "weakest-coupling" regime: (1') diffusive and local, and (2') diffusive and nonlocal. We derived the analytic expression of the topological Hall conductivity (THC) for each region, and found that the condition for the locality of the effective field is governed by transverse spin diffusion of electrons. In region 1, where the spin relaxation is negligible and the effective field is local, the THC is found to be proportional to $M$, instead of $M^3$ of the weakest-coupling regime. In the diffusive, nonlocal regions (2 and 2'), the effective field is given by a spin chirality formed by "effective spins" that the electrons see during their diffusive motion. Applying the results to a skyrmion lattice, we found the THC decreses as the skyrmion density is increased in region 2', reflecting the nonlocality of the effective field, and shows a maximum at the boundary to the "local" region.


title: Absolute and arbitrary orientation of single molecule shapes

abstract: DNA origami is a modular platform for the combination of molecular and colloidal components to create optical, electronic, and biological devices. Integration of such nanoscale devices with microfabricated connectors and circuits is challenging: large numbers of freely diffusing devices must be fixed at desired locations with desired alignment. We present a DNA origami molecule whose energy landscape on lithographic binding sites has a unique maximum. This property enables device alignment within 3.2$^{\circ}$ on SiO$_2$. Orientation is absolute (all degrees of freedom are specified) and arbitrary (every molecule's orientation is independently specified). The use of orientation to optimize device performance is shown by aligning fluorescent emission dipoles within microfabricated optical cavities. Large-scale integration is demonstrated via an array of 3,456 DNA origami with 12 distinct orientations, which indicates the polarization of excitation light.

categories: physics.app-ph cond-mat.mes-hall cond-mat.soft physics.atm-clus physics.bio-ph


title: MT-VAE: Learning Motion Transformations to Generate Multimodal Human Dynamics

abstract: Long-term human motion can be represented as a series of motion modes---motion sequences that capture short-term temporal dynamics---with transitions between them. We leverage this structure and present a novel Motion Transformation Variational Auto-Encoders (MT-VAE) for learning motion sequence generation. Our model jointly learns a feature embedding for motion modes (that the motion sequence can be reconstructed from) and a feature transformation that represents the transition of one motion mode to the next motion mode. Our model is able to generate multiple diverse and plausible motion sequences in the future from the same input. We apply our approach to both facial and full body motion, and demonstrate applications like analogy-based motion transfer and video synthesis.


title: Trapping low-mass planets at the inner edge of the protostellar disc

abstract: The formation of multiple close-in low-mass exoplanets is still a mystery. The challenge is to build a system wherein the outermost planet is beyond 0.2 AU from the star. Here we investigate how the prescription for type I planet migration affects the ability to trap multiple planets in a resonant chain near the inner edge of the protostellar disc. A sharp edge modelled as a hyperbolic tangent function coupled with supersonic corrections to the classical type I migration torques results in the innermost planets being pushed inside the cavity through resonant interaction with farther planets because migration is starward at slightly supersonic eccentricities. Planets below a few Earth masses are generally trapped in a resonant chain with the outermost planet near the disc edge, but long-term stability is not guaranteed. For more massive planets the migration is so fast that the eccentricity of the innermost resonant pair is excited to highly supersonic levels due to decreased damping on the innermost planet as it is pushed inside the cavity; collisions frequently occur and the system consists one or two intermediate-mass planets residing closer to the star than the disc's inner edge. We found a neat pileup of resonant planets outside the disc edge only if the corotation torque does not rapidly diminish at high eccentricity. We call for detailed studies on planet migration near the disc's inner edge, which is still uncertain, and for an improved understanding of eccentricity damping and disc torques in the supersonic regime.


title: Machine Learning for Heterogeneous Ultra-Dense Networks with Graphical Representations

abstract: Heterogeneous ultra-dense network (H-UDN) is envisioned as a promising solution to sustain the explosive mobile traffic demand through network densification. By placing access points, processors, and storage units as close as possible to mobile users, H-UDNs bring forth a number of advantages, including high spectral efficiency, high energy efficiency, and low latency. Nonetheless, the high density and diversity of network entities in H-UDNs introduce formidable design challenges in collaborative signal processing and resource management. This article illustrates the great potential of machine learning techniques in solving these challenges. In particular, we show how to utilize graphical representations of H-UDNs to design efficient machine learning algorithms.


title: Tuning the topological insulator states of artificial graphene

abstract: We develop a robust, non-perturbative approach to study the band structure of artificial graphene. Artificial graphene, as considered here, is generated by imposing a superlattice structure on top of a two dimensional hole gas in a semiconductor heterostructure, where the hole gas naturally possesses large spin-orbit coupling. Via tuning of the system parameters we demonstrate how best to exploit the spin-orbit coupling to generate time reversal symmetry-protected topological insulator phases. Our major conclusion is the identification of a second set of topological Dirac bands in the band structure (with spin Chern number $C=3$), which were not reliably obtainable in previous perturbative approaches to artificial graphene. Importantly, the second Dirac bands host more desirable features than the previously studied first set of Dirac bands (with $C=1$). Moreover, we find that upon tuning of the system parameters, we can drive the system to the highly desirable regime of the topological flat band. We discuss the possibilities this opens up for exotic, strongly correlated phases.


title: Multiscale core-periphery structure in a global liner shipping network

abstract: Maritime transport accounts for a majority of trades in volume, of which 70% in value is carried by container ships that transit regular routes on fixed schedules in the ocean. In the present paper, we analyse a data set of global liner shipping as a network of ports. In particular, we construct the network of the ports as the one-mode projection of a bipartite network composed of ports and ship routes. Like other transportation networks, global liner shipping networks may have core-periphery structure, where a core and a periphery are groups of densely and sparsely interconnected nodes, respectively. Core-periphery structure may have practical implications for understanding the robustness, efficiency and uneven development of international transportation systems. We develop an algorithm to detect core-periphery pairs in a network, which allows one to find core and peripheral nodes on different scales and uses a configuration model that accounts for the fact that the network is obtained by the one-mode projection of a bipartite network. We also found that most ports are core (as opposed to peripheral) ports and that ports in some countries in Europe, America and Asia belong to a global core-periphery pair across different scales, whereas ports in other countries do not.


title: SciSports: Learning football kinematics through two-dimensional tracking data

abstract: SciSports is a Dutch startup company specializing in football analytics. This paper describes a joint research effort with SciSports, during the Study Group Mathematics with Industry 2018 at Eindhoven, the Netherlands. The main challenge that we addressed was to automatically process empirical football players' trajectories, in order to extract useful information from them. The data provided to us was two-dimensional positional data during entire matches. We developed methods based on Newtonian mechanics and the Kalman filter, Generative Adversarial Nets and Variational Autoencoders. In addition, we trained a discriminator network to recognize and discern different movement patterns of players. The Kalman-filter approach yields an interpretable model, in which a small number of player-dependent parameters can be fit; in theory this could be used to distinguish among players. The Generative-Adversarial-Nets approach appears promising in theory, and some initial tests showed an improvement with respect to the baseline, but the limits in time and computational power meant that we could not fully explore it. We also trained a Discriminator network to distinguish between two players based on their trajectories; after training, the network managed to distinguish between some pairs of players, but not between others. After training, the Variational Autoencoders generated trajectories that are difficult to distinguish, visually, from the data. These experiments provide an indication that deep generative models can learn the underlying structure and statistics of football players' trajectories. This can serve as a starting point for determining player qualities based on such trajectory data.


title: Moving Object Segmentation in Jittery Videos by Stabilizing Trajectories Modeled in Kendall's Shape Space

abstract: Moving Object Segmentation is a challenging task for jittery/wobbly videos. For jittery videos, the non-smooth camera motion makes discrimination between foreground objects and background layers hard to solve. While most recent works for moving video object segmentation fail in this scenario, our method generates an accurate segmentation of a single moving object. The proposed method performs a sparse segmentation, where frame-wise labels are assigned only to trajectory coordinates, followed by the pixel-wise labeling of frames. The sparse segmentation involving stabilization and clustering of trajectories in a 3-stage iterative process. At the 1st stage, the trajectories are clustered using pairwise Procrustes distance as a cue for creating an affinity matrix. The 2nd stage performs a block-wise Procrustes analysis of the trajectories and estimates Frechet means (in Kendall's shape space) of the clusters. The Frechet means represent the average trajectories of the motion clusters. An optimization function has been formulated to stabilize the Frechet means, yielding stabilized trajectories at the 3rd stage. The accuracy of the motion clusters are iteratively refined, producing distinct groups of stabilized trajectories. Next, the labels obtained from the sparse segmentation are propagated for pixel-wise labeling of the frames, using a GraphCut based energy formulation. Use of Procrustes analysis and energy minimization in Kendall's shape space for moving object segmentation in jittery videos, is the novelty of this work. Second contribution comes from experiments performed on a dataset formed of 20 real-world natural jittery videos, with manually annotated ground truth. Experiments are done with controlled levels of artificial jitter on videos of SegTrack2 dataset. Qualitative and quantitative results indicate the superiority of the proposed method.


title: The elementary construction of formal anafunctors

abstract: This article gives an elementary and formal 2-categorical construction of a bicategory of right fractions analogous to anafunctors, starting from a 2-category equipped with a family of covering maps that are fully faithful and co-fully faithful.


title: Spectral bounds for non-uniform hypergraphs using weighted clique expansion

abstract: Hypergraphs are an invaluable tool to understand many hidden patterns in large data sets. Among many ways to represent hypergraph, one useful representation is that of weighted clique expansion. In this paper, we consider this representation for non-uniform hypergraphs. We generalize the spectral results for uniform hypergraphs to non-uniform hypergraphs and show that they extend in a natural way. We provide a bound on the largest eigenvalue with respect to the average degree of neighbours of a vertex in a graph. We also prove an inequality on the boundary of a vertex set in terms of the largest and second smallest eigenvalue and use it to obtain bounds on some connectivity parameters of the hypergraph.


title: Three-particle electron-hole complexes in two-dimensional electron systems

abstract: Three-particle complexes consisting of two holes in the completely filled zero electron Landau level and an excited electron in the unoccupied first Landau level are investigated in a quantum Hall insulator. The distinctive features of these three-particle complexes are an electron-hole mass symmetry and the small energy gap of the quantum Hall insulator itself. Theoretical calculations of the trion energy spectrum in a quantizing magnetic field predict that, besides the ground state, trions feature a hierarchy of excited bound states. In agreement with the theoretical simulations, we observe new photoluminescence lines related to the excited trion states. A relatively small energy gap allows the binding of three-particle complexes with magnetoplasma oscillations and formation of plasmarons. The plasmaron properties are investigated experimentally.


title: Branched projective structures on a Riemann surface and logarithmic connections

abstract: We study the set ${\mathcal P}_S$ consisting of all branched holomorphic projective structures on a compact Riemann surface $X$ of genus $g \geq 1$ and with a fixed branching divisor $S:= \sum_{i=1}^d n_i\cdot x_i$, where $x_i \in X$. Under the hypothesis that $n_i=1$, for all $i$, with $d$ a positive even integer such that $d \neq 2g-2$, we show that ${\mathcal P}_S$ coincides with a subset of the set of all logarithmic connections with singular locus $S$, satisfying certain geometric conditions, on the rank two holomorphic jet bundle $J^1(Q)$, where $Q$ is a fixed holomorphic line bundle on $X$ such that $Q^{\otimes 2}= TX\otimes {\mathcal O}_X(S)$. The space of all logarithmic connections of the above type is an affine space over the vector space $H^0(X, K^{\otimes 2}_X \otimes {\mathcal O}_X(S))$ of dimension $3g-3+d$. We conclude that ${\mathcal P}_S$ is a subset of this affine space that has codimenison $d$ at a generic point.


title: Laser writing of bright colours on near-percolation plasmonic reflector arrays

abstract: Colouration by surface nanostructuring has attracted a great deal of attention by the virtue of making use of environment-friendly recyclable materials and generating non-bleaching colours [1-8]. Recently, it was found possible to delegate the task of colour printing to laser post-processing that modifies carefully designed and fabricated nanostructures [9,10]. Here we take the next crucial step in the development of structural colour printing by dispensing with preformed nanostructures and using instead near-percolation metal films atop dielectric-metal sandwiches, i.e., near-percolation plasmonic reflector arrays. Scanning rapidly (~ 20 {\mu}m/s) across 4-nm-thin island-like gold films supported by 30-nm-thin silica layers atop 100-nm-thick gold layers with a strongly focused Ti-sapphire laser beam, while adjusting the average laser power from 1 to 10 mW, we produce bright colours varying from green to red by laser-heating-induced merging and reshaping of gold islands. Selection of strongly heated islands and their reshaping, both originating from the excitation of plasmonic resonances, are strongly influenced by the polarization direction of laser illumination, so that the colours produced are well pronounced only when viewed with the same polarization. Conversely, the laser colour writing with circular polarizations results in bright polarization-independent colour images. The fabrication procedure for near-percolation reflector arrays is exceedingly simple and scalable to mass production, while the laser-induced modification occurs inherently with the subwavelength resolution. This unique combination of remarkable features makes the approach developed for laser colour writing readily amenable for practical implementation and use in diverse applications ranging from nanoscale patterning for security marking to large-scale colour printing for decoration.


title: A Method for Quickly Bounding the Optimal Objective Value of an OPF Problem using a Semidefinite Relaxation and a Local Solution

abstract: Optimal power flow (OPF) is an important problem in the operation of electric power systems. Due to the OPF problem's non-convexity, there may exist multiple local optima. Certifiably obtaining the global solution is important for certain applications of OPF problems. Many global optimization techniques compute an optimality gap that compares the achievable objective value corresponding to the feasible point from a local solution algorithm with the objective value bound from a convex relaxation technique. Rather than the traditional practice of completely separating the local solution and convex relaxation computations, this paper proposes a method that exploits information from a local solution to speed the computation of an objective value bound using a semidefinite programming (SDP) relaxation. The improvement in computational tractability comes with the trade-off of reduced tightness for the resulting objective value bound. Numerical experiments illustrate this trade-off, with the proposed method being faster but weaker than the SDP relaxation and slower but tighter than second-order cone programming (SOCP) and quadratic convex (QC) relaxations for many large test cases.


title: Explicit construction of optimal locally recoverable codes of distance 5 and 6 via binary constant weight codes

abstract: It was shown in \cite{GXY18} that the length $n$ of a $q$-ary linear locally recoverable code with distance $d\ge 5$ is upper bounded by $O(dq^3)$. Thus, it is a challenging problem to construct $q$-ary locally recoverable codes with distance $d\ge 5$ and length approaching the upper bound. The paper \cite{GXY18} also gave an algorithmic construction of $q$-ary locally recoverable codes with locality $r$ and length $n=\Omega_r(q^2)$ for $d=5$ and $6$, where $\Omega_r$ means that the implicit constant depends on locality $r$. In the present paper, we present an explicit construction of $q$-ary locally recoverable codes of distance $d= 5$ and $6$ via binary constant weight codes. It turns out that (i) our construction is simpler and more explicit; and (ii) lengths of our codes are larger than those given in \cite{GXY18}.


title: Time-Multiplexed Parsing in Marking-based Network Telemetry

abstract: Network telemetry is a key capability for managing the health and efficiency of a large-scale network. Alternate Marking Performance Measurement (AM-PM) is a recently introduced approach that accurately measures the packet loss and delay in a network using a small overhead of one or two bits per data packet. This paper introduces a novel time-multiplexed parsing approach that enables a practical and accurate implementation of AM-PM in network devices, while requiring just a single bit per packet. Experimental results are presented, based on a hardware implementation, and a software P4-based implementation.


title: Deep Retinex Decomposition for Low-Light Enhancement

abstract: Retinex model is an effective tool for low-light image enhancement. It assumes that observed images can be decomposed into the reflectance and illumination. Most existing Retinex-based methods have carefully designed hand-crafted constraints and parameters for this highly ill-posed decomposition, which may be limited by model capacity when applied in various scenes. In this paper, we collect a LOw-Light dataset (LOL) containing low/normal-light image pairs and propose a deep Retinex-Net learned on this dataset, including a Decom-Net for decomposition and an Enhance-Net for illumination adjustment. In the training process for Decom-Net, there is no ground truth of decomposed reflectance and illumination. The network is learned with only key constraints including the consistent reflectance shared by paired low/normal-light images, and the smoothness of illumination. Based on the decomposition, subsequent lightness enhancement is conducted on illumination by an enhancement network called Enhance-Net, and for joint denoising there is a denoising operation on reflectance. The Retinex-Net is end-to-end trainable, so that the learned decomposition is by nature good for lightness adjustment. Extensive experiments demonstrate that our method not only achieves visually pleasing quality for low-light enhancement but also provides a good representation of image decomposition.


title: Commutation matrices and Commutation tensors

abstract: The commutation matrix was first introduced in statistics as a transposition matrix by Murnaghan in 1938. In this paper, we first investigate the commutation matrix which is employed to transform a matrix into its transpose. We then extend the concept of the commutation matrix to commutation tensor and use the commutation tensor to achieve the unification of the two formulae of the linear preserver of the matrix rank, a classical result of Marcus in 1971.


title: Algebraic Properties of Quasilinear Two-Dimensional Lattices connected with integrability

abstract: In the article a classification method for nonlinear integrable equations with three independent variables is discussed based on the notion of the integrable reductions. We call the equation integrable if it admits a large class of reductions being Darboux integrable systems of hyperbolic type equations with two independent variables. The most natural and convenient object to be studied within the frame of this scheme is the class of two dimensional lattices generalizing the well-known Toda lattice. In the present article we deal with the quasilinear lattices of the form $u_{n,xy}=\alpha(u_{n+1} ,u_n,u_{n-1} )u_{n,x}u_{n,y} + \beta(u_{n+1},u_n,u_{n-1})u_{n,x}+\gamma(u_{n+1} ,u_n,u_{n-1} )u_{n,y}+\delta(u_{n+1} ,u_n,u_{n-1})$. We specify the coefficients of the lattice assuming that there exist cutting off conditions which reduce the lattice to a Darboux integrable hyperbolic type system of the arbitrarily high order. Under some extra assumption of nondegeneracy we described the class of the lattices integrable in the sense indicated above. There are new examples in the obtained list of chains.


title: Selective-area chemical beam epitaxy of in-plane InAs one-dimensional channels grown on InP(001), InP(111)B, and InP(110) surfaces

abstract: We report on the selective-area chemical beam epitaxial growth of InAs in-plane, one-dimensional (1-D) channels using patterned SiO$_{2}$-coated InP(001), InP(111)B, and InP(110) substrates to establish a scalable platform for topological superconductor networks. Top-view scanning electron micrographs show excellent surface selectivity and dependence of major facet planes on the substrate orientations and ridge directions, and the ratios of the surface energies of the major facet planes were estimated. Detailed structural properties and defects in the InAs nanowires (NWs) were characterized by transmission electron microscopic analysis of cross-sections perpendicular to the NW ridge direction and along the NW ridge direction. Electrical transport properties of the InAs NWs were investigated using Hall bars, a field effect mobility device, a quantum dot, and an Aharonov-Bohm loop device, which reflect the strong spin-orbit interaction and phase-coherent transport characteristic in the selectively grown InAs systems. This study demonstrates that selective-area chemical beam epitaxy is a scalable approach to realize semiconductor 1-D channel networks with the excellent surface selectivity and this material system is suitable for quantum transport studies.


title: In Search of Recent Disruption of (3200) Phaethon: Model Implication and Hubble Space Telescope Search

abstract: Near-Earth asteroid (3200) Phaethon is notable for its association to a strong annual meteor shower, the Geminids, indicative of one or more episodes of mass ejection in the past. The mechanism of Phaethon's past activity is not yet understood. Here we present a Hubble Space Telescope (HST) search of meter-sized fragments in the vicinity of Phaethon, carried out during Phaethon's historic approach to the Earth in mid-December of 2017. Numerical simulation conducted to guide HST's pointing also show that the dynamical evolution of Phaethon-originated particles is quick, as ejected materials take no longer than $\sim250$ yr to spread to the entire orbit of Phaethon. Our search was completed down to 4-meter-class limit (assuming Phaethon-like albedo) and was expected to detect 0.035% particles ejected by Phaethon in the last several decades. The negative result of our search capped the total mass loss of Phaethon over the past few dozen orbits to be $10^{12}$ kg at $3\sigma$ level, taking the best estimates of size power-law from meteor observations and spacecraft data. Our result also implies a millimeter-sized dust flux of $<10^{-12} \mathrm{m^{-2} s^{-1}}$ within 0.1 au of Phaethon, suggesting that any Phaethon-bound mission is unlikely to encounter dense dust clouds.


title: Performance of a GridPix detector based on the Timepix3 chip

abstract: A GridPix readout for a TPC based on the Timepix3 chip is developed for future applications at a linear collider. The GridPix detector consists of a gaseous drift volume read out by a single Timepix3 chip with an integrated amplification grid. Its performance is studied in a test beam with 2.5 GeV electrons. The GridPix detector detects single ionization electrons with high efficiency. The Timepix3 chip allowed for high sample rates and time walk corrections. Diffusion is found to be the dominating error on the track position measurement both in the pixel plane and in the drift direction, and systematic distortions in the pixel plane are below 10 $\mu$m. Using a truncated sum, an energy loss (dE/dx) resolution of 4.1% is found for an effective track length of 1 m.


title: AKARI mid-infrared slitless spectroscopic survey of star-forming galaxies at $z\lesssim 0.5$

abstract: Deep MIR surveys have revealed numerous strongly star-forming galaxies at redshift z~<2. Their MIR fluxes are produced by a combination of continuum and PAH emission features. The PAH features can dominate the total MIR flux, but are difficult to measure without spectroscopy. We aim to study star-forming galaxies by using a blind spectroscopic survey at MIR wavelengths to understand evolution of their star formation rate (SFR) and specific SFR up to z~=0.5, by paying particular attention to their PAH properties. We conducted a low-resolution (R~=50) slitless spectroscopic survey at 5-13um of 9um flux-selected sources (>0.3 mJy) around the North Ecliptic Pole with the Infrared Camera (IRC) onboard AKARI. After removing 11 AGN candidates by using the IRC photometry, we identified 48 PAH galaxies with PAH 6.2, 7.7, and 8.6um features at z<0.5. The rest-frame optical-MIR spectral energy distributions (SEDs) based on CFHT and AKARI/IRC imaging covering 0.37-18um were produced, and analysed in conjunction with the PAH spectroscopy. We defined the PAH enhancement by using the luminosity ratio of the 7.7um PAH feature over the 3.5um stellar component of the SEDs. The rest-frame SEDs of all PAH galaxies have a universal shape with stellar and 7.7um bumps, except that the PAH enhancement significantly varies as a function of the PAH luminosities. We identified a PAH-enhanced population at z~>0.35, whose SEDs and luminosities are typical of luminous infrared galaxies. They show particularly larger PAH enhancement at high luminosity, implying that they are vigorous star-forming galaxies with elevated specific SFR. Our composite starburst model that combines a very young and optically very thick starburst with a very old population can successfully reproduce most of their SED characteristics, although we could not confirm this optically thick component from our spectral analysis.


title: Geometry and symmetry in quantum Boltzmann machine

abstract: Quantum Boltzmann machine extends the classical Boltzmann machine learning to the quantum regime, which makes its power to simulate the quantum states beyond the classical probability distributions. We develop the BFGS algorithm to study the corresponding optimization problem in quantum Boltzmann machine, especially focus on the target states being a family of states with parameters. As an typical example, we study the target states being the real symmetric two-qubit pure states, and we find two obvious features shown in the numerical results on the minimal quantum relative entropy: First, the minimal quantum relative entropy in the first and the third quadrants is zero; Second, the minimal quantum relative entropy is symmetric with the axes $y=x$ and $y=-x$ even with one qubit hidden layer. Then we theoretically prove these two features from the geometric viewpoint and the symmetry analysis. Our studies show that the traditional physical tools can be used to help us to understand some interesting results from quantum Boltzmann machine learning.


title: A Hamiltonian for Massless and Zero-Energy States

abstract: We present a non-hermitian Hamiltonian which can be employed to explain a condensed matter system with effectively massless and zero energy states. We analyze the 2D tunneling problem and derive the transmission and reflection coefficients for the massless and zero energy states. We find that the transmission coefficient of the massless and zero-energy particles in this case is consistent with non-chiral tunneling of quasiparticles which is in contrast to the Klein tunneling known for electrons described by the Dirac equation. Our analysis predicts that the massless electron can be reflected as a hole-like state whereas this transition does not occur for the zero-energy state. Experimental observations are needed to test whether the presented Hamiltonian can be realized in such a condensed matter system.


title: Metric for Rotating object in Infrared Corrected Nonlocal Gravity Model

abstract: Here, we derive the metric for the spacetime around rotating object for the gravity action having nonlocal correction of $R\Box^{-2} R $ to the Einstein-Hilbert action. Starting with the generic stationary, axisymmetric metric, we solve the equations of motion in linearized gravity limit for the modified action including energy-momentum tensor of the rotating mass. We also derive the rotating metric from the static metric using the Demanski-Janis-Newmann algorithm. Finally, we obtain the constraint on the value of $M$ by calculating the frame dragging effect in our theory and comparing it with that of General Relativity and Gravity Probe B results, where $M$ is the mass scale of the theory.


title: The Peculiarities of the Cosmological Models Based on Non-Linear Classical and Phantom Fields with Minimal Interaction. I. The Cosmological Model Based on Scalar Singlet

abstract: A detailed comparative qualitative analysis and numerical simulation of evolution of the cosmological models based on classical and phantom scalar fields with self-action was performed. The phase portraits of the dynamic systems of classical and phantom fields and their projections to the Poincare sphere were constructed. It was shown that the phase trajectories of the corresponding dynamic systems can be split by bifurcation trajectories into 2,4 or 6 different dynamic streams corresponding to different pairwise symmetric histories of the Universe depending on the parameters of the scalar field's model. The phase space of such systems becomes multiply connected, the ranges of negative total effective energy unavailable for motion, getting appear there. In the case when attracting centers are situated inside these ranges, the phase trajectories of the classical scalar field in the infinite future tend to limit cycles, winding onto the boundaries of these ranges. The phase trajectories of the scalar field, in turn, get away from the boundaries of ranges with null effective energy and in the infinite future are wound onto one of the symmetrical focuses (centers). Thus, the situations when the Universe, in case of classical scalar field, begins its history with the inflation and ends it up in the Euclidean future, or, in the case of phantom scalar field, in opposite, has the Euclidian start and proceeds to inflation mode after anomalous burst of the acceleration, both become possible. The potentials of scalar fields on the surface of null curvature are distinct from zero and thereby define a certain vacuum state.


title: Learning A Shared Transform Model for Skull to Digital Face Image Matching

abstract: Human skull identification is an arduous task, traditionally requiring the expertise of forensic artists and anthropologists. This paper is an effort to automate the process of matching skull images to digital face images, thereby establishing an identity of the skeletal remains. In order to achieve this, a novel Shared Transform Model is proposed for learning discriminative representations. The model learns robust features while reducing the intra-class variations between skulls and digital face images. Such a model can assist law enforcement agencies by speeding up the process of skull identification, and reducing the manual load. Experimental evaluation performed on two pre-defined protocols of the publicly available IdentifyMe dataset demonstrates the efficacy of the proposed model.


title: Small Sample Learning in Big Data Era

abstract: As a promising area in artificial intelligence, a new learning paradigm, called Small Sample Learning (SSL), has been attracting prominent research attention in the recent years. In this paper, we aim to present a survey to comprehensively introduce the current techniques proposed on this topic. Specifically, current SSL techniques can be mainly divided into two categories. The first category of SSL approaches can be called "concept learning", which emphasizes learning new concepts from only few related observations. The purpose is mainly to simulate human learning behaviors like recognition, generation, imagination, synthesis and analysis. The second category is called "experience learning", which usually co-exists with the large sample learning manner of conventional machine learning. This category mainly focuses on learning with insufficient samples, and can also be called small data learning in some literatures. More extensive surveys on both categories of SSL techniques are introduced and some neuroscience evidences are provided to clarify the rationality of the entire SSL regime, and the relationship with human learning process. Some discussions on the main challenges and possible future research directions along this line are also presented.


title: Effect of filter shape on excess noise performance in continuous variable quantum key distribution with Gaussian modulation

abstract: An attractive implementation for quantum cryptography is the continuous variable variation, as it relies on standard telecommunication components. Modulating the quantum signal using a Gaussian format is attractive since it has been proven to be secure. This work investigates the effect of the roll-off of a root raised cosine pulse shaping and matched filter on the excess noise performance of a Gaussian modulated quantum key distribution system in a simulated back to back configuration. Contrary to intuition, it is found that the roll-off parameter does not significantly impact the performance of the system.


title: Threshold and jet radius joint resummation for single-inclusive jet production

abstract: The QCD predictions for single-inclusive jet production are computed with joint resummation of threshold and jet radius logarithms. The results are compared to those based on fixed order perturbation theory up to next-to-next-to-leading order and to data by the CMS collaboration measured in proton-proton collisions at the LHC at $\sqrt{S}=8$ TeV. The joint resummation results are in remarkable agreement with the data.


title: Ubiquitous ram pressure stripping in the Coma cluster of galaxies

abstract: We report the detection of H$\alpha$ trails behind three new intermediate-mass irregular galaxies in the NW outskirts of the nearby cluster of galaxies Abell 1656 (Coma). Hints that these galaxies possess an extended component were found in earlier, deeper H$\alpha$ observations carried out with the Subaru telescope. However the lack of a simultaneous $r$-band exposure, together with the presence of strong stellar ghosts in the Subaru images, prevented us from quantifying the detections. We therefore devoted one full night of H$\alpha$ observation to each of the three galaxies using the San Pedro Martir 2.1m telescope. One-sided tails of H$\alpha$ emission of 10-20 kpc projected size were detected, suggesting an ongoing ram pressure stripping event. We added these 3 new sources of extended ionized gas (EIG) added to the 12 found by Yagi et al. (2010), NGC 4848 (Fossati et al. 2012), and NGC 4921 whose ram pressure stripping is certified by HI asymmetry. This brings the number sources with H$\alpha$ trails to 17 gaseous tails out of 27 (63 \%) late-type galaxies (LTG) galaxies members of the Coma cluster with direct evidence of ram pressure stripping. The 27 LTG galaxies, among these the 17 with extended H$\alpha$ tails, have kinematic properties that are different from the rest of the early-type galaxy (ETG) population of the c ore of the Coma cluster, as they deviate in the phase-space diagram $\Delta$V/$\sigma$ versus $r/R_{200}$.


title: Automatic Airway Segmentation in chest CT using Convolutional Neural Networks

abstract: Segmentation of the airway tree from chest computed tomography (CT) images is critical for quantitative assessment of airway diseases including bronchiectasis and chronic obstructive pulmonary disease (COPD). However, obtaining an accurate segmentation of airways from CT scans is difficult due to the high complexity of airway structures. Recently, deep convolutional neural networks (CNNs) have become the state-of-the-art for many segmentation tasks, and in particular the so-called Unet architecture for biomedical images. However, its application to the segmentation of airways still remains a challenging task. This work presents a simple but robust approach based on a 3D Unet to perform segmentation of airways from chest CTs. The method is trained on a dataset composed of 12 CTs, and tested on another 6 CTs. We evaluate the influence of different loss functions and data augmentation techniques, and reach an average dice coefficient of 0.8 between the ground-truth and our automated segmentations.


title: Two constructions of optimal pairs of linear codes for resisting side channel and fault injection attacks

abstract: Direct sum masking (DSM) has been proposed as a counter-measure against side-channel attacks (SCA) and fault injection attacks (FIA), which are nowadays important domains of cryptanalysis. DSM needs two linear codes whose sum is direct and equals a whole space $\Bbb F_q^n$. The minimum distance of the former code and the dual distance of the latter should be as large as possible, given their length and dimensions. But the implementation needs in practice to work with words obtained by appending, to each codeword $y$ of the latter code, the source word from which $y$ is the encoding. Let $\mathcal C_1$ be an $[n, k]$ linear code over the finite field $\Bbb F_q$ with generator matrix $G$ and let $\mathcal C_2$ be the linear code over the finite field $\Bbb F_q$ with generator matrix $[G, I_k]$. It is then highly desired to construct optimal pairs of linear codes satisfying that $d(\mathcal C_2^\perp)= d(\mathcal C_1^\perp)$. In this paper, we employ the primitive irreducible cyclic codes to derive two constructions of optimal pairs of linear codes for resisting SCA and FIA, where the security parameters are determined explicitly. To the best of our knowledge, it is the first time that primitive irreducible cyclic codes are used to construct (optimal) pairs of codes. As a byproduct, we obtain the weight enumerators of the codes $\mathcal C_1, \mathcal C_2, \mathcal C_1^\perp$, and $\mathcal C_2^\perp$ in our both constructions.


title: A note on eigenvalue bounds for Schr\"odinger operators

abstract: We obtain a new bound on the location of eigenvalues for a non-self-adjoint Schr\"odinger operator with complex-valued potentials by obtaining a weighted $L^2$ estimate for the resolvent of the Laplacian.


title: Bringing Together Dynamic Geometry Software and the Graphics Processing Unit

abstract: We equip dynamic geometry software (DGS) with a user-friendly method that enables massively parallel calculations on the graphics processing unit (GPU). This interplay of DGS and GPU opens up various applications in education and mathematical research. The GPU-aided discovery of mathematical properties, interactive visualizations of algebraic surfaces (raycasting), the mathematical deformation of images and footage in real-time, and computationally demanding numerical simulations of PDEs are examples from the long and versatile list of new domains that our approach makes accessible within a DGS. We ease the development of complex (mathematical) visualizations and provide a rapid-prototyping scheme for general-purpose computations (GPGPU). The possibility to program both CPU and GPU with the use of only one high-level (scripting) programming language is a crucial aspect of our concept. We embed shader programming seamlessly within a high-level (scripting) programming environment. The aforementioned requires the symbolic process of the transcompilation of a high-level programming language into shader programming language for GPU and, in this article, we address the challenge of the automatic translation of a high-level programming language to a shader language of the GPU. To maintain platform independence and the possibility to use our technology on modern devices, we focus on a realization through WebGL.


title: NFFT meets Krylov methods: Fast matrix-vector products for the graph Laplacian of fully connected networks

abstract: The graph Laplacian is a standard tool in data science, machine learning, and image processing. The corresponding matrix inherits the complex structure of the underlying network and is in certain applications densely populated. This makes computations, in particular matrix-vector products, with the graph Laplacian a hard task. A typical application is the computation of a number of its eigenvalues and eigenvectors. Standard methods become infeasible as the number of nodes in the graph is too large. We propose the use of the fast summation based on the nonequispaced fast Fourier transform (NFFT) to perform the dense matrix-vector product with the graph Laplacian fast without ever forming the whole matrix. The enormous flexibility of the NFFT algorithm allows us to embed the accelerated multiplication into Lanczos-based eigenvalues routines or iterative linear system solvers and even consider other than the standard Gaussian kernels. We illustrate the feasibility of our approach on a number of test problems from image segmentation to semi-supervised learning based on graph-based PDEs. In particular, we compare our approach with the Nystr\"om method. Moreover, we present and test an enhanced, hybrid version of the Nystr\"om method, which internally uses the NFFT.


title: ACE of Spades in the IoT Security Game: A Flexible IPsec Security Profile for Access Control

abstract: The Authentication and Authorization for Constrained Environments (ACE) framework provides fine-grained access control in the Internet of Things, where devices are resource-constrained and with limited connectivity. The ACE framework defines separate profiles to specify how exactly entities interact and what security and communication protocols to use. This paper presents the novel ACE IPsec profile, which specifies how a client establishes a secure IPsec channel with a resource server, contextually using the ACE framework to enforce authorized access to remote resources. The profile makes it possible to establish IPsec Security Associations, either through their direct provisioning or through the standard IKEv2 protocol. We provide the first Open Source implementation of the ACE IPsec profile for the Contiki OS and test it on the resource-constrained Zolertia Firefly platform. Our experimental performance evaluation confirms that the IPsec profile and its operating modes are affordable and deployable also on constrained IoT platforms.


title: The joint distribution of pin-point plant cover data: a reparametrized Dirichlet -- multinomial distribution

abstract: A reparametrized Dirichlet-multinomial distribution is introduced, and the covariance matrix, as well as, the algorithm for calculating the PDF for n species are provided. The distribution is suited for modelling the joint distribution of pin-point cover data of spatially aggregated plant species, and the parametrization ensures that the degree of spatially aggregation is modelled by a parameter and the mean cover of the species is modelled by the other parameters. This last property is convenient for using the distribution in Bayesian hierarchical models where the mean cover of the different species typically is modelled as latent variables.


title: The Pauli sum rules imply BSM physics

abstract: Some 67 years ago (1951) Wolfgang Pauli mooted the three sum rules: \[ \sum_n (-1)^{2S_n} g_n = 0; \qquad \sum_n (-1)^{2S_n} g_n \; m_n^2 =0; \qquad \sum_n (-1)^{2S_n} g_n \; m_n^4=0. \] These three sum rules are intimately related to both the Lorentz invariance and the finiteness of the zero-point stress-energy tensor. Further afield, these three constraints are also intimately related to the existence of finite QFTs ultimately based on fermi--bose cancellations. (Supersymmetry is neither necessary nor sufficient for the existence of these finite QFTs; though softly but explicitly broken supersymmetry or mis-aligned supersymmetry can be used as a book-keeping device to keep the calculations manageable.) In the current article I shall instead take these three Pauli sum rules as given, assume their exact non-perturbative validity, contrast them with the observed standard model particle physics spectrum, and use them to extract as much model-independent information as possible regarding beyond standard model (BSM) physics.


title: On the possible mechanism of radio emission of polars

abstract: We suggest a mechanism for generation of radio emission from polars. It is based on the cyclotron radiation of thermal electrons at the background of fluctuating magnetic field. The source of fluctuations is Alfv\'{e}n wave turbulence. Expressions for the radiation spectrum and degree of polarization are obtained. By the example of the polar AM~Her, the fluxes of radio emission from the accretion stream are calculated. Within the framework of proposed emission model, with realistic plasma characteristics, it is possible to reproduce observed radiation fluxes in the VLA bands.


title: Optimal additive Schwarz preconditioning for adaptive 2D IGA boundary element methods

abstract: We define and analyze (local) multilevel diagonal preconditioners for isogeometric boundary elements on locally refined meshes in two dimensions. Hypersingular and weakly-singular integral equations are considered. We prove that the condition number of the preconditioned systems of linear equations is independent of the mesh-size and the refinement level. Therefore, the computational complexity, when using appropriate iterative solvers, is optimal. Our analysis is carried out for closed and open boundaries and numerical examples confirm our theoretical results.


title: On the Brun spectral sequence for topological Hochschild homology

abstract: We generalize a spectral sequence of Brun for the computation of topological Hochschild homology. The generalized version computes the $E$-homology of $THH(A;B)$, where $E$ is a ring spectrum, $A$ is a commutative $S$-algebra and $B$ is a connective commutative $A$-algebra. The input of the spectral sequence are the topological Hochschild homology groups of $B$ with coefficients in the $E$-homology groups of $B \wedge_A B$. The mod $p$ and $v_1$ topological Hochschild homology of connective complex $K$-theory has been computed by Ausoni and later again by Rognes, Sagave and Schlichtkrull. We present an alternative, short computation using the generalized Brun spectral sequence.