title: Geometrical contributions to the exchange constants: Free electrons with spin-orbit interaction

abstract: Using thermal quantum field theory we derive an expression for the exchange constant that resembles Fukuyama's formula for the orbital magnetic susceptibility (OMS). Guided by this formal analogy between the exchange constant and OMS we identify a contribution to the exchange constant that arises from the geometrical properties of the band structure in mixed phase space. We compute the exchange constants for free electrons and show that the geometrical contribution is generally important. Our formalism allows us to study the exchange constants in the presence of spin-orbit interaction (SOI). Thereby, we find sizable differences between the exchange constants of helical and cycloidal spin spirals. Furthermore, we discuss how to calculate the exchange constants based on a gauge-field approach in the case of the Rashba model with an additional exchange splitting and show that the exchange constants obtained from this gauge-field approach are in perfect agreement with those obtained from the quantum field theoretical method.


title: Curvature dependence of relativistic epicyclic frequencies in static, axially symmetric spacetimes

abstract: The sum of squared epicyclic frequencies of nearly circular motion ($\omega_r^2+\omega_\theta^2$) in axially symmetric configurations of Newtonian gravity is known to depend both on the matter density and on the angular velocity profile of circular orbits. It was recently found that this sum goes to zero at the photon orbits of Schwarzschild and Kerr spacetimes. However, these are the only relativistic configurations for which such result exists in the literature. Here, we extend the above formalism in order to describe the analogous relation for geodesic motion in arbitrary static, axially symmetric, asymptotically flat solutions of general relativity. The sum of squared epicyclic frequencies is found to vanish at photon radii of vacuum solutions. In the presence of matter, we obtain that $\omega_r^2+\omega_\theta^2>0$ for perturbed timelike circular geodesics on the equatorial plane if the strong energy condition holds for the matter-energy fluid of spacetime; in vacuum, the allowed region for timelike circular geodesic motion is characterized by the inequality above. The results presented here may be of use to shed light on general issues concerning the stability of circular orbits once they approach photon radii, mainly the ones corresponding to stable photon motion.


title: Terahertz and higher-order Brunel harmonics: from tunnel to multiphoton ionization regime in tailored fields

abstract: Brunel radiation appears as a result of a two-step process of photo-ionization and subsequent acceleration of electron, without the need of electron recollision. We show that for generation of Brunel harmonics at all frequencies the subcycle ionization dynamics is of critical importance. Namely, such harmonics disappear at low pump intensities when the ionization dynamics depends only on the slow envelope (so called multiphoton ionization regime) and not on the instantaneous field. Nevertheless, if the pump pulse contains incommensurate frequencies, Brunel mechanism does generate new frequencies even in the multiphoton ionization regime.


title: Density-functional calculations of multivalency-driven formation of Te-based monolayer materials with superior electronic and optical properties

abstract: Contemporary science is witnessing a rapid expansion of the two-dimensional (2D) materials family, each member possessing intriguing emergent properties of fundamental and practical importance. Using the particle-swarm optimization method in combination with first-principles density functional theory calculations, here wepredict a new category of 2D monolayers named tellurene, composed of the metalloid element Te, with stable 1T-MoS2-like ( {\alpha}-Te), and metastable tetragonal (\b{eta}-Te) and 2H-MoS2-like ({\gamma}-Te) structures. The underlying formation mechanism of such tri-layer arrangements is uniquely rooted in the multivalent nature of Te, with the central-layer Te behaving more metal-like (e.g., Mo), and the two outer layers more semiconductor-like (e.g.,S). In particular, the {\alpha}-Te phase can be spontaneously obtained from the magic thicknesses truncated along the [001] direction of the trigonal structure of bulk Te. Furthermore, both the {\alpha}- and \b{eta}-Te phases possess electron and hole mobilities much higher than MoS2, as well as salient optical absorption properties. These findings effectively extend the realm of 2D materials to group-VI monolayers, and provide a new and generic formation mechanism for designing 2D materials.


title: Implications for Post-Processing Nucleosynthesis of Core-Collapse Supernova Models with Lagrangian Particles

abstract: We investigate core-collapse supernova (CCSN) nucleosynthesis with self-consistent, axisymmetric (2D) simulations performed using the radiation-hydrodynamics code Chimera. Computational costs have traditionally constrained the evolution of the nuclear composition within multidimensional CCSN models to, at best, a 14-species $\alpha$-network capable of tracking only $(\alpha,\gamma)$ reactions from $^{4}$He to $^{60}$Zn. Such a simplified network limits the ability to accurately evolve detailed composition and neutronization or calculate the nuclear energy generation rate. Lagrangian tracer particles are commonly used to extend the nuclear network evolution by incorporating more realistic networks in post-processing nucleosynthesis calculations. However, limitations such as poor spatial resolution of the tracer particles, inconsistent thermodynamic evolution, including misestimation of expansion timescales, and uncertain determination of the multidimensional mass-cut at the end of the simulation impose uncertainties inherent to this approach. We present a detailed analysis of the impact of such uncertainties for four self-consistent axisymmetric CCSN models initiated from stellar metallicity, non-rotating progenitors of 12 $M_\odot$, 15 $M_\odot$, 20 $M_\odot$, and 25 $M_\odot$ and evolved with the smaller $\alpha$-network to more than 1 s after the launch of an explosion.


title: 1.5 bit-per-stage 8-bit Pipelined CMOS A/D Converter for Neuromophic Vision Processor

abstract: Neuromorphic vision processor is an electronic implementation of vision algorithm processor on semiconductor. To image the world, a low-power CMOS image sensor array is required in the vision processor. The image sensor array is typically formed through photo diodes and analog to digital converter (ADC). To achieve low power acquisition, a low-power mid-resolution ADC is necessary. In this paper, a 1.8V, 8-bit, 166MS/s pipelined ADC was proposed in a 0.18 um CMOS technology. The ADC used operational amplifier sharing architecture to reduce power consumption and achieved maximum DNL of 0.24 LSB, maximum INL of 0.35 LSB, at a power consumption of 38.9mW. When input frequency is 10.4MHz, it achieved an SNDR 45.9dB, SFDR 50dB, and an ENOB of 7.33 bit.


title: Deep Reinforcement Learning for Robotic Manipulation-The state of the art

abstract: The focus of this work is to enumerate the various approaches and algorithms that center around application of reinforcement learning in robotic ma- ]]nipulation tasks. Earlier methods utilized specialized policy representations and human demonstrations to constrict the policy. Such methods worked well with continuous state and policy space of robots but failed to come up with generalized policies. Subsequently, high dimensional non-linear function approximators like neural networks have been used to learn policies from scratch. Several novel and recent approaches have also embedded control policy with efficient perceptual representation using deep learning. This has led to the emergence of a new branch of dynamic robot control system called deep r inforcement learning(DRL). This work embodies a survey of the most recent algorithms, architectures and their implementations in simulations and real world robotic platforms. The gamut of DRL architectures are partitioned into two different branches namely, discrete action space algorithms(DAS) and continuous action space algorithms(CAS). Further, the CAS algorithms are divided into stochastic continuous action space(SCAS) and deterministic continuous action space(DCAS) algorithms. Along with elucidating an organ- isation of the DRL algorithms this work also manifests some of the state of the art applications of these approaches in robotic manipulation tasks.


title: Robotic Haptic Proxies for Collaborative Virtual Reality

abstract: We propose a new approach for interaction in Virtual Reality (VR) using mobile robots as proxies for haptic feedback. This approach allows VR users to have the experience of sharing and manipulating tangible physical objects with remote collaborators. Because participants do not directly observe the robotic proxies, the mapping between them and the virtual objects is not required to be direct. In this paper, we describe our implementation, various scenarios for interaction, and a preliminary user study.


title: Pulsating low-mass white dwarfs in the frame of new evolutionary sequences: IV. The secular rate of period change

abstract: We present a theoretical assessment of the expected temporal rates of change of periods ($\dot{\Pi}$) for low-mass ($M_{\star}/M_{\sun} \lesssim 0.45$) and extremely low-mass (ELM, $M_{\star}/M_{\sun} \lesssim 0.18-0.20$) white-dwarf stars, based on fully evolutionary low-mass He-core white dwarf and pre-white dwarf models. Our analysis is based on a large set of adiabatic periods of radial and nonradial pulsation modes computed on a suite of low-mass He-core white dwarf and pre-white dwarf models with masses ranging from $0.1554$ to $0.4352 M_{\sun}$. We compute the secular rates of period change of radial ($\ell= 0$) and nonradial ($\ell= 1, 2$) $g$ and $p$ modes for stellar models representative of ELMV and pre-ELMV stars, as well as for stellar objects that are evolving just before the occurrence of CNO flashes at the early cooling branches. We found that the theoretically expected magnitude of $\dot{\Pi}$ of $g$ modes for pre-ELMVs are by far larger than for ELMVs. In turn, $\dot{\Pi}$ of $g$ modes for models evolving before the occurrence of CNO flashes are larger than the maximum values of the rates of period change predicted for pre-ELMV stars. Regarding $p$ and radial modes, we found that the larger absolute values of $\dot{\Pi}$ correspond to pre-ELMV models. We conclude that any eventual measurement of a rate of period change for a given pulsating low-mass pre-white dwarf or white dwarf star could shed light about its evolutionary status. Also, in view of the systematic difficulties in the spectroscopic classification of stars of the ELM Survey, an eventual measurement of $\dot{\Pi}$ could help to confirm that a given pulsating star is an authentic low-mass white dwarf and not a star from another stellar population.


title: A Chandra X-ray census of the interacting binaries in old open clusters - Collinder 261

abstract: We present the first X-ray study of Collinder 261 (Cr261), which at an age of 7 Gyr is one of the oldest open clusters known in the Galaxy. Our observation with the Chandra X-ray Observatory is aimed at uncovering the close interacting binaries in Cr261, and reaches a limiting X-ray luminosity of Lx ~ 4e29 erg/s (0.3-7 keV) for stars in the cluster. We detect 107 sources within the cluster half-mass radius r_h, and we estimate that among the sources with Lx >~ 1e30 erg/s, ~26 are associated with the cluster. We identify a mix of active binaries and candidate active binaries, candidate cataclysmic variables, and stars that have "straggled" from the main locus of Cr261 in the colour-magnitude diagram. Based on a deep optical source catalogue of the field, we estimate that Cr261 has an approximate mass of 6500 M_sun, roughly the same as the old open cluster NGC6791. The X-ray emissivity of Cr261 is similar to that of other old open clusters, supporting the trend that they are more luminous in X-rays per unit mass than old populations of higher (globular clusters) and lower (the local neighbourhood) stellar density. This implies that the dynamical destruction of binaries in the densest environments is not solely responsible for the observed differences in X-ray emissivity.


title: High Energy Gamma Rays from Nebulae Associated with Extragalactic Microquasars and Ultra-Luminous X-ray Sources

abstract: In the extragalactic sky, microquasars and ultra-luminous X-ray sources (ULXs) are known as energetic compact objects locating at off-nucleus positions in galaxies. Some of these objects are associated with expanding bubbles with a velocity of 80-250 ${\rm km~s^{-1}}$. We investigate the shock acceleration of particles in those expanding nebulae. The nebulae having fast expansion velocity $\gtrsim120~{\rm km~s^{-1}}$ are able to accelerate cosmic rays up to $\sim100$ TeV. If 10% of the shock kinetic energy goes into particle acceleration, powerful nebulae such as the microquasar S26 in NGC 7793 would emit gamma rays up to several tens TeV with a photon index of $\sim2$. These nebulae will be good targets for future Cherenkov Telescope Array observations given its sensitivity and angular resolution. They would also contribute to $\sim7$% of the unresolved cosmic gamma-ray background radiation at $\ge0.1~{\rm GeV}$. In contrast, particle acceleration in slowly expanding nebulae $\lesssim120~{\rm km~s^{-1}}$ would be less efficient due to ion-neutral collisions and result in softer spectra at $\gtrsim10$ GeV.


title: A Covert Queueing Channel in Round Robin Schedulers

abstract: We study a covert queueing channel (CQC) between two users sharing a round robin scheduler. Such a covert channel can arise when users share a resource such as a computer processor or a router arbitrated by a round robin policy. We present an information-theoretic framework to model and derive the maximum reliable data transmission rate, i.e., the capacity of this channel, for both noiseless and noisy setups. Our results show that seemingly isolated users can communicate with a high rate over the covert channel and demonstrate the possibility of significant information leakage and privacy threats brought by CQCs in round robin schedulers. Moreover, we propose practical finite-length code constructions, which achieve the capacity limit.


title: Enabling large-scale viscoelastic calculations via neural network acceleration

abstract: One of the most significant challenges involved in efforts to understand the effects of repeated earthquake cycle activity are the computational costs of large-scale viscoelastic earthquake cycle models. Computationally intensive viscoelastic codes must be evaluated thousands of times and locations, and as a result, studies tend to adopt a few fixed rheological structures and model geometries, and examine the predicted time-dependent deformation over short (<10 yr) time periods at a given depth after a large earthquake. Training a deep neural network to learn a computationally efficient representation of viscoelastic solutions, at any time, location, and for a large range of rheological structures, allows these calculations to be done quickly and reliably, with high spatial and temporal resolution. We demonstrate that this machine learning approach accelerates viscoelastic calculations by more than 50,000%. This magnitude of acceleration will enable the modeling of geometrically complex faults over thousands of earthquake cycles across wider ranges of model parameters and at larger spatial and temporal scales than have been previously possible.


title: Spectropolarimetry of the 2012 outburst of SN 2009ip: a bi-polar explosion in a dense, disk-like CSM

abstract: We present a sequence of eight spectropolarimetric observations monitoring the geometric evolution of the late phase of the major 2012 outburst of SN 2009ip. These were acquired with the FORS2 polarimeter mounted on ESO VLT. The continuum was polarised at 0.3-0.8 per cent throughout the observations, showing that the photosphere deviated substantially from spherical symmetry by 10-15 per cent. Significant line polarisation is detected for both hydrogen and helium at high velocities. The similarity in the polarised signal between these elements indicates that they form in the same location in the ejecta. The line polarisation (p$\sim$1-1.5 per cent) at low velocities revealed the presence of a highly-aspherical hydrogen and helium rich circumstellar medium (CSM). Monte Carlo simulations of the observed polarimetry were performed in an effort to constrain the shape of the CSM. The simulations imply that the polarimetry can be understood within the framework of a disk-like CSM inclined by 14$\pm$2 degrees out of the line of sight, obscuring the photosphere only at certain epochs. The varying temporal evolution of polarisation at high and low velocities indicated that the fast-moving ejecta expanded with a preferred direction orthogonal to that of the CSM.


title: SenseGen: A Deep Learning Architecture for Synthetic Sensor Data Generation

abstract: Our ability to synthesize sensory data that preserves specific statistical properties of the real data has had tremendous implications on data privacy and big data analytics. The synthetic data can be used as a substitute for selective real data segments,that are sensitive to the user, thus protecting privacy and resulting in improved analytics.However, increasingly adversarial roles taken by data recipients such as mobile apps, or other cloud-based analytics services, mandate that the synthetic data, in addition to preserving statistical properties, should also be difficult to distinguish from the real data. Typically, visual inspection has been used as a test to distinguish between datasets. But more recently, sophisticated classifier models (discriminators), corresponding to a set of events, have also been employed to distinguish between synthesized and real data. The model operates on both datasets and the respective event outputs are compared for consistency. In this paper, we take a step towards generating sensory data that can pass a deep learning based discriminator model test, and make two specific contributions: first, we present a deep learning based architecture for synthesizing sensory data. This architecture comprises of a generator model, which is a stack of multiple Long-Short-Term-Memory (LSTM) networks and a Mixture Density Network. second, we use another LSTM network based discriminator model for distinguishing between the true and the synthesized data. Using a dataset of accelerometer traces, collected using smartphones of users doing their daily activities, we show that the deep learning based discriminator model can only distinguish between the real and synthesized traces with an accuracy in the neighborhood of 50%.


title: A southern-sky total intensity source catalogue at 2.3 GHz from S-band Polarisation All-Sky Survey data

abstract: The S-band Polarisation All-Sky Survey (S-PASS) has observed the entire southern sky using the 64-metre Parkes radio telescope at 2.3GHz with an effective bandwidth of 184MHz. The surveyed sky area covers all declinations $\delta\leq 0^\circ$. To analyse compact sources the survey data have been re-processed to produce a set of 107 Stokes $I$ maps with 10.75arcmin resolution and the large scale emission contribution filtered out. In this paper we use these Stokes $I$ images to create a total intensity southern-sky extragalactic source catalogue at 2.3GHz. The source catalogue contains 23,389 sources and covers a sky area of 16,600deg$^2$, excluding the Galactic plane for latitudes $|b|<10^\circ$. Approximately 8% of catalogued sources are resolved. S-PASS source positions are typically accurate to within 35arcsec. At a flux density of 225mJy the S-PASS source catalogue is more than 95% complete, and $\sim$94% of S-PASS sources brighter than 500mJy beam$^{-1}$ have a counterpart at lower frequencies.


title: Integrating Reviews into Personalized Ranking for Cold Start Recommendation

abstract: Item recommendation task predicts a personalized ranking over a set of items for each individual user. One paradigm is the rating-based methods that concentrate on explicit feedbacks and hence face the difficulties in collecting them. Meanwhile, the ranking-based methods are presented with rated items and then rank the rated above the unrated. This paradigm takes advantage of widely available implicit feedback. It, however, usually ignores a kind of important information: item reviews. Item reviews not only justify the preferences of users, but also help alleviate the cold-start problem that fails the collaborative filtering. In this paper, we propose two novel and simple models to integrate item reviews into Bayesian personalized ranking. In each model, we make use of text features extracted from item reviews using word embeddings. On top of text features we uncover the review dimensions that explain the variation in users' feedback and these review factors represent a prior preference of users. Experiments on six real-world data sets show the benefits of leveraging item reviews on ranking prediction. We also conduct analyses to understand the proposed models.


title: On a possibility to extract a signal from heavy gluino cascade decays via isolated muon detection at the LHC energies

abstract: 800 GeV gluino pair production is simulated using ISASUSY 1.1 MC-program. The cascade decay of the gluino into two lightest neutralinos, $\widetilde{g} \rightarrow \widetilde{\chi}_2^0 X$, $\widetilde{\chi}_2^0 \rightarrow \mu^+ \mu^- \widetilde{\chi}_1^0$, is studied for a convenient set of SUSY parameters: $\mu$ = 200 GeV, tan$\beta$ = 2, average squark masses of 1600 GeV, charged Higgs mass of 500 GeV, and top quark mass of 140 GeV. With such parameters the mass difference between $\widetilde{\chi}_2^0$ and $\widetilde{\chi}_1^0$ is less than the mass of Z. As a background, $t\bar{t}$, Z + jets, WW, WZ and ZZ events were generated using ISAJET 6.50 MC-program. The granularities and energy resolutions of the calorimeters in the pseudorapidity range of $\mid\eta\mid$ < 5 are taken into account. If missing transverse momentum is greater than 200 GeV, the number of jets in the calorimeters is 4 or more, the invariant mass of two isolated muons is less than 80 GeV, the background in the muon transverse momentum distribution is negligible.


title: An integrated approach to grey relational analysis, analytic hierarchy process and data envelopment analysis

abstract: Purpose: This paper aims to propose an integration of the analytic hierarchy process (AHP) and data envelopment analysis (DEA) methods in a multiattribute grey relational analysis (GRA) methodology in which the attribute weights are completely unknown and the attribute values take the form of fuzzy numbers. Design/methodology/approach: This research has been organized to proceed along the following steps: computing the grey relational coefficients for alternatives with respect to each attribute using a fuzzy GRA methodology. Grey relational coefficients provide the required (output) data for additive DEA models; computing the priority weights of attributes using the AHP method to impose weight bounds on attribute weights in additive DEA models; computing grey relational grades using a pair of additive DEA models to assess the performance of each alternative from the optimistic and pessimistic perspectives; and combining the optimistic and pessimistic grey relational grades using a compromise grade to assess the overall performance of each alternative. Findings: The proposed approach provides a more reasonable and encompassing measure of performance, based on which the overall ranking position of alternatives is obtained. An illustrated example of a nuclear waste dump site selection is used to highlight the usefulness of the proposed approach. Originality/value: This research is a step forward to overcome the current shortcomings in the weighting schemes of attributes in a fuzzy multiattribute GRA methodology. Keywords: Data envelopment analysis, Grey relational analysis, Analytic hierarchy process, Multiple attribute decision making, Fuzzy numbers, Weighting.


title: Covert Communication with Finite Blocklength in AWGN Channels

abstract: Covert communication is to achieve a reliable transmission from a transmitter to a receiver while guaranteeing an arbitrarily small probability of this transmission being detected by a warden. In this work, we study the covert communication in AWGN channels with finite blocklength, in which the number of channel uses is finite. Specifically, we analytically prove that the entire block (all available channel uses) should be utilized to maximize the effective throughput of the transmission subject to a predetermined covert requirement. This is a nontrivial result because more channel uses results in more observations at the warden for detecting the transmission. We also determine the maximum allowable transmit power per channel use, which is shown to decrease as the blocklength increases. Despite the decrease in the maximum allowable transmit power per channel use, the maximum allowable total power over the entire block is proved to increase with the blocklength, which leads to the fact that the effective throughput increases with the blocklength.


title: Reshetnyak rigidity for Riemannian manifolds

abstract: We prove two rigidity theorems for maps between Riemannian manifolds. First, we prove that a Lipschitz map $f:M\to N$ between two oriented Riemannian manifolds, whose differential is almost everywhere an orientation-preserving isometry, is an isometric immersion. This theorem was previously proved using regularity theory for conformal maps; we give a new, simple proof, by generalizing the Piola identity for the cofactor operator. Second, we prove that if there exists a sequence of mapping $f_n:M\to N$, whose differentials converge in $L^p$ to the set of orientation-preserving isometries, then there exists a subsequence converging to an isometric immersion. These results are generalizations of celebrated rigidity theorems by Liouville (1850) and Reshetnyak (1967) from Euclidean to Riemannian settings. Finally, we describe applications of these theorems to non-Euclidean elasticity and to convergence notions of manifolds.


title: Stable and Controllable Neural Texture Synthesis and Style Transfer Using Histogram Losses

abstract: Recently, methods have been proposed that perform texture synthesis and style transfer by using convolutional neural networks (e.g. Gatys et al. [2015,2016]). These methods are exciting because they can in some cases create results with state-of-the-art quality. However, in this paper, we show these methods also have limitations in texture quality, stability, requisite parameter tuning, and lack of user controls. This paper presents a multiscale synthesis pipeline based on convolutional neural networks that ameliorates these issues. We first give a mathematical explanation of the source of instabilities in many previous approaches. We then improve these instabilities by using histogram losses to synthesize textures that better statistically match the exemplar. We also show how to integrate localized style losses in our multiscale framework. These losses can improve the quality of large features, improve the separation of content and style, and offer artistic controls such as paint by numbers. We demonstrate that our approach offers improved quality, convergence in fewer iterations, and more stability over the optimization.


title: Reclassifying symbiotic stars with 2MASS and WISE: An atlas of spectral energy distribution

abstract: We present a new updated catalogue of Galactic and extragalactic symbiotic stars (SySts). Since the last catalogue of SySts (Belczynski et al. 2000), the number of known SySts has significantly increased. Our new catalogue contains 316 known and 82 candidates SySts. Of the confirmed Systs 252 are located in our Galaxy and 64 in nearby galaxies. This reflects an increase of $\sim$50% in the population of Galactic SySts and $\sim$400% in the population of extragalactic SySts. The spectral energy distribution (SED) of 334 (known and candidates) SySts have been constructed using the 2MASS and WISE data. These SEDs are used to provide a robust reclassification in scheme of S- (74%), D- (15%) and D'-types (2.5%). The SEDs of S- and D-type peak between 0.8 and 1.6$\mu$m and between 1.6 and 4$\mu$m, respectively, whereas those of D'-type exhibit a plateau profile. Moreover, we provide the first compilation of SySts that exhibit the OVI Raman-scattered line at 6830\AA. Our analysis shows that 55% of the Galactic SySts exhibit that line in their spectrum, whereas this percentage is different from galaxy to galaxy.


title: Chentsov's theorem for exponential families

abstract: Chentsov's theorem characterizes the Fisher information metric on statistical models as essentially the only Riemannian metric that is invariant under sufficient statistics. This implies that each statistical model is naturally equipped with a geometry, so Chentsov's theorem explains why many statistical properties can be described in geometric terms. However, despite being one of the foundational theorems of statistics, Chentsov's theorem has only been proved previously in very restricted settings or under relatively strong regularity and invariance assumptions. We therefore prove a version of this theorem for the important case of exponential families. In particular, we characterise the Fisher information metric as the only Riemannian metric (up to rescaling) on an exponential family and its derived families that is invariant under independent and identically distributed extensions and canonical sufficient statistics. Our approach is based on the central limit theorem, so it gives a unified proof for both discrete and continuous exponential families, and it is less technical than previous approaches.

categories: math.ST cs.IT math.DG math.IT math.PR stat.TH


title: On the Role of a Market Maker in Networked Cournot Competition

abstract: We study Cournot competition among firms in a networked marketplace that is centrally managed by a market maker. In particular, we study a situation in which a market maker facilitates trade between geographically separate markets via a constrained transport network. Our focus is on understanding the consequences of the design of the market maker and on providing tools for optimal design. To that end we provide a characterization of the equilibrium outcomes of the game between the firms and the market maker. Our results highlight that the equilibrium structure is impacted dramatically by the market maker's objective - depending on the objective there may be a unique equilibrium, multiple equilibria, or no equilibria. Further, the game may be a potential game (as in the case of classical Cournot competition) or not. Beyond characterizing the equilibria of the game, we provide an approach for designing the market maker in order to optimize a design objective (e.g., social welfare) at the equilibrium of the game. Additionally, we use our results to explore the value of transport (trade) and the efficiency of the market maker (as compared to a single, aggregate market).


title: Computing a tree having a small vertex cover

abstract: We consider a new Steiner tree problem, called vertex-cover-weighted Steiner tree problem. This problem defines the weight of a Steiner tree as the minimum weight of vertex covers in the tree, and seeks a minimum-weight Steiner tree in a given vertex-weighted undirected graph. Since it is included by the Steiner tree activation problem, the problem admits an O(log n)-approximation algorithm in general graphs with n vertices. This approximation factor is tight up to a constant because it is NP-hard to achieve an o(log n)-approximation for the vertex-cover-weighted Steiner tree problem on general graphs even if the given vertex weights are uniform and a spanning tree is required instead of a Steiner tree. In this paper, we present constant-factor approximation algorithms for the problem with unit disk graphs and with graphs excluding a fixed minor. For the latter graph class, our algorithm can be also applied for the Steiner tree activation problem.


title: OVI 6830\AA\ Imaging Polarimetry of Symbiotic Stars

abstract: I present here the first results from an ongoing pilot project with the 1.6 m telescope at the OPD, Brasil, aimed at the detection of the OVI $\lambda$6830 line via linear polarization in symbiotic stars. The main goal is to demonstrate that OVI imaging polarimetry is an efficient technique for discovering new symbiotic stars. The OVI $\lambda$6830 line is found in 5 out of 9 known symbiotic stars, in which the OVI line has already been spectroscopically confirmed, with at least 3-$\sigma$ detection. Three new symbiotic star candidates have also been found.


title: Nested Hilbert schemes on surfaces: Virtual fundamental class

abstract: We construct virtual fundamental classes on nested Hilbert schemes of points and curves in complex nonsingular projective surfaces. These classes recover the virtual classes of Seiberg-Witten theory as well as the (reduced) stable theory, and play a crucial role in the reduced Donaldson-Thomas theory of local-surface-threefolds that we study in [GSY17b] (arXiv:1807.05697). We show that certain integrals against the virtual fundamental classes of punctual nested Hilbert schemes are expressed as integrals over the products of the Hilbert scheme of points. We are able to find explicit formulas for some of these integrals by relating them to Carlsson-Okounkov's vertex operator formulas.


title: On likely solutions of the stable matching problem with unequal numbers of men and women

abstract: Following up a recent work by Ashlagi, Kanoria and Leshno, we study a stable matching problem with unequal numbers of men and women, and independent uniform preferences. The asymptotic formulas for the expected number of stable matchings, and for the probabilities of one point--concentration for the range of husbands' total ranks and for the range of wives' total ranks are obtained.


title: Irreducibility of quantum representations of mapping class groups with boundary

abstract: We prove that the Witten--Reshetikhin--Turaev $\mathrm{SU}(2)$ quantum representations of mapping class groups are always irreducible in the case of surfaces equipped with colored banded points, provided that at least one banded point is colored by one. We thus generalize a well--known result due to J. Roberts.


title: Localized Donaldson-Thomas theory of surfaces

abstract: Let $S$ be a projective simply connected complex surface and $\mathcal{L}$ be a line bundle on $S$. We study the moduli space of stable compactly supported 2-dimensional sheaves on the total spaces of $\mathcal{L}$. The moduli space admits a $\mathbb{C}^*$-action induced by scaling the fibers of $\mathcal{L}$. We identify certain components of the fixed locus of the moduli space with the moduli space of torsion free sheaves and the nested Hilbert schemes on $S$. We define the localized Donaldson-Thomas invariants of $\mathcal{L}$ by virtual localization in the case that $\mathcal{L}$ twisted by the anti-canonical bundle of $S$ admits a nonzero global section. When $p_g(S)>0$, in combination with Mochizuki's formulas, we are able to express the localized DT invariants in terms of the invariants of the nested Hilbert schemes defined by the authors in [GSY17a], the Seiberg-Witten invariants of $S$, and the integrals over the products of Hilbert schemes of points on $S$. When $\mathcal{L}$ is the canonical bundle of $S$, the Vafa-Witten invariants defined recently by Tanaka-Thomas, can be extracted from these localized DT invariants. VW invariants are expected to have modular properties as predicted by S-duality.


title: The emergence of torsion in the continuum limit of distributed edge-dislocations - erratum

abstract: This note refers to our previous paper "The emergence of torsion in the continuum limit of distributed edge-dislocations". It identifies and fixes an error in the notion of convergence of Weitzenb\"ock manifolds defined in the paper, and in the proof of the well-definiteness of this notion of convergence.


title: Mott transition and magnetism of the triangular-lattice Hubbard model with next-nearest-neighbor hopping

abstract: The variational cluster approximation is used to study the isotropic triangular-lattice Hubbard model at half filling, taking into account the nearest-neighbor ($t_1$) and next-nearest-neighbor ($t_2$) hopping parameters for magnetic frustrations. We determine the ground-state phase diagram of the model. In the strong correlation regime, the 120$^\circ$ N\'eel and stripe ordered phases appear, and a nonmagnetic insulating phase emerges in between. In the intermediate correlation regime, the nonmagnetic insulating phase expands to a wider parameter region, which goes into a paramagnetic metallic phase in the weak correlation regime. The critical phase boundary of the Mott metal-insulator transition is discussed in terms of the van Hove singularity evident in the calculated density of states and single-particle spectral function.


title: High-Order Localized Dissipation Weighted Compact Nonlinear Scheme for Shock- and Interface-Capturing in Compressible Flows

abstract: Simulations of single- and multi-species compressible flows with shock waves and discontinuities are conducted using a weighted compact nonlinear scheme (WCNS) with a newly developed sixth order localized dissipative interpolation. In smooth regions, the scheme applies the central nonlinear interpolation with minimum dissipation to well resolve fluctuating flow features while in regions containing discontinuities and high wavenumber features, the scheme suppresses spurious numerical oscillations by hybridizing the central interpolation with the more dissipative upwind-biased nonlinear interpolation. In capturing material interfaces between species of different densities, a quasi-conservative five equation model that can conserve mass of each species is used to prevent pressure oscillations across the interfaces. Compared to upwind-biased interpolations with classical nonlinear weights and improved weights, and the interpolation with adaptive central-upwind weights for scale-separation, it is shown that WCNS with the proposed localized dissipative interpolation has better performance to simultaneously capture discontinuities and resolve smooth features.


title: On-off Analog Beamforming for Massive MIMO

abstract: This paper investigates a new analog beamforming architecture for massive multiple-input multiple-output (MIMO) systems, where each of the multiple transmit antennas is switched to be on or off to form a beam according to the channel state information at transmitters. This on-off analog beamforming (OABF) scheme has the advantages in terms of both hardware complexities and algorithmic complexities. OABF can completely remove the high-cost, power-consuming and bulky analog phase-shifters that are extensively employed by traditional analog beamforming schemes; it only requires the deployment of low-cost analog switches that are easy to implement. Moreover, we show that the beams formed by such simple antenna on-off switch operations can achieve rather good performances with low-complexity beamforming algorithms. Specifically, we first propose two optimal signal-to-noise ratio maximization algorithms to determine the on-off state of each switch under the per-antenna power constraint and the total power constraint, respectively. After that, we theoretically prove that OABF can achieve the full diversity gain and the full array gain with complexities up to a polynomial order. Numerical results are consistent with our theoretical analysis. We believe that the simple structure of OABF makes massive MIMO much easier to implement in real systems.


title: Laser annealing heals radiation damage in avalanche photodiodes

abstract: Avalanche photodiodes (APDs) are a practical option for space-based quantum communications requiring single-photon detection. However, radiation damage to APDs significantly increases their dark count rates and reduces their useful lifetimes in orbit. We show that high-power laser annealing of irradiated APDs of three different models (Excelitas C30902SH, Excelitas SLiK, and Laser Components SAP500S2) heals the radiation damage and substantially restores low dark count rates. Of nine samples, the maximum dark count rate reduction factor varies between 5.3 and 758 when operating at minus 80 degrees Celsius. The illumination power to reach these reduction factors ranges from 0.8 to 1.6 W. Other photon detection characteristics, such as photon detection efficiency, timing jitter, and afterpulsing probability, remain mostly unaffected. These results herald a promising method to extend the lifetime of a quantum satellite equipped with APDs.


title: Ultra-thin Acoustic Metasurface-Based Schroeder Diffuser

abstract: Schroeder diffuser is a classical design, proposed over 40 years ago, for artificially creating optimal and predictable sound diffuse reflection. It has been widely adopted in architectural acoustics and it has also shown substantial potential in noise control, ultrasound imaging, microparticle manipulation, among others. The conventional Schroeder diffuser, however, has a considerable thickness on the order of one wavelength, severely impeding its applications for low frequency sound. In this paper, a new class of ultra-thin and planar Schroeder diffusers are proposed based on the concept of acoustic metasurface. Both numerical and experimental results demonstrate satisfactory sound diffuse reflection produced from the metasurface-based Schroeder diffuser despite it being one order of magnitude thinner than the conventional one. The proposed design not only offer promising building blocks with great potential to profoundly impact architectural acoustics and related fields, but also constitutes a major step towards real-world applications of acoustic metasurfaces.


title: Controlling charge transport mechanisms in molecular junctions: Distilling thermally-induced hopping from coherent-resonant conduction

abstract: The electrical conductance of molecular junctions may strongly depend on the temperature, and weakly on molecular length, under two distinct mechanisms: phase-coherent resonant conduction, with charges proceeding via delocalized molecular orbitals, and incoherent thermally-assisted multi-step hopping. While in the case of coherent conduction the temperature dependence arises from the broadening of the Fermi distribution in the metal electrodes, in the latter case it corresponds to electron-vibration interaction effects on the junction. With the objective to distill the thermally-activated hopping component, thus expose intrinsic electron-vibration interaction phenomena on the junction, we suggest the design of molecular junctions with "spacers", extended anchoring groups that act to filter out phase-coherent resonant electrons. Specifically, we study the electrical conductance of fixed-gap and variable-gap junctions that include a tunneling block, with spacers at the boundaries. Using numerical simulations and analytical considerations, we demonstrate that in our design, resonant conduction is suppressed. As a result, the electrical conductance is dominated by two (rather than three) mechanisms: superexchange (deep tunneling), and multi-step thermally-induced hopping. We further exemplify our analysis on DNA junctions with an A:T block serving as a tunneling barrier. Here, we show that the electrical conductance is insensitive to the number of G:C base-pairs at the boundaries. This indicates that the tunneling-to-hopping crossover revealed in such sequences truly corresponds to the properties of the A:T barrier.


title: Convex Constrained Semialgebraic Volume Optimization: Application in Systems and Control

abstract: In this paper, we generalize the chance optimization problems and introduce constrained volume optimization where enables us to obtain convex formulation for challenging problems in systems and control. We show that many different problems can be cast as a particular cases of this framework. In constrained volume optimization, we aim at maximizing the volume of a semialgebraic set under some semialgebraic constraints. Building on the theory of measures and moments, a sequence of semidefinite programs are provided, whose sequence of optimal values is shown to converge to the optimal value of the original problem. We show that different problems in the area of systems and control that are known to be nonconvex can be reformulated as special cases of this framework. Particularly, in this work, we address the problems of probabilistic control of uncertain systems as well as inner approximation of region of attraction and invariant sets of polynomial systems. Numerical examples are presented to illustrate the computational performance of the proposed approach.


title: A New Torsion Pendulum for Gravitational Reference Sensor Technology Development

abstract: We report on the design and sensitivity of a new torsion pendulum for measuring the performance of ultra-precise inertial sensors and for the development of associated technologies for space-based gravitational wave observatories and geodesy missions. The apparatus comprises a 1 m-long, 50 um-diameter, tungsten fiber that supports an inertial member inside a vacuum system. The inertial member is an aluminum crossbar with four hollow cubic test masses at each end. This structure converts the rotation of the torsion pendulum into translation of the test masses. Two test masses are enclosed in capacitive sensors which provide readout and actuation. These test masses are electrically insulated from the rest of the cross-bar and their electrical charge is controlled by photoemission using fiber-coupled ultraviolet light emitting diodes. The capacitive readout measures the test mass displacement with a broadband sensitivity of 30 nm / sqrt(Hz), and is complemented by a laser interferometer with a sensitivity of about 0.5 nm / sqrt(Hz). The performance of the pendulum, as determined by the measured residual torque noise and expressed in terms of equivalent force acting on a single test mass, is roughly 200 fN / sqrt(Hz) around 2 mHz, which is about a factor of 20 above the thermal noise limit of the fiber.


title: Electronic structure study of screw dislocation core energetics in Aluminum and core energetics informed forces in a dislocation aggregate

abstract: We use a real-space formulation of orbital-free DFT to study the core energetics and core structure of an isolated screw dislocation in Aluminum. Using a direct energetics based approach, we estimate the core size of a perfect screw dislocation to be $\approx$ 7 $|{\bf b}|$, which is considerably larger than previous estimates of $1-3~|{\bf b}|$ based on displacement fields. The perfect screw upon structural relaxation dissociates into two Shockley partials with partial separation distances of 8.2 \AA~and 6.6 \AA~ measured from the screw and edge component differential displacement plots, respectively. Similar to a previous electronic structure study on edge dislocation, we find that the core energy of the relaxed screw dislocation is strongly dependent on macroscopic deformations. Next, we use the edge and screw dislocation core energetics data with physically reasonable assumptions to develop a continuum energetics model for an aggregate of dislocations. Further, we use this model in a discrete dislocation network, and from the variations of the core energy with respect to the nodal positions of the network, we obtain the nodal core force which can directly be incorporated into discrete dislocation dynamics frameworks. We analyze and classify the nodal core force into three different contributions based on their decay behavior. Two of these contributions to the core force, both arising from the core energy dependence on macroscopic deformations, are not accounted for in currently used discrete dislocation dynamics models which assume the core energy to be a constant. Using case studies involving simple dislocation structures, we demonstrate that the contribution to the core force from the core energy dependence on macroscopic deformations can be significant in comparison to the elastic Peach-Koehler force even up to distances of $10-15$ nm between dislocation structures.


title: Topological recursion of Eynard-Orantin and the Harmonic Oscillator

abstract: We apply the Chekhov-Eynard-Orantin topological recursion to the curve corresponding to the quantum harmonic oscillator and demonstrate that the result is equivalent to the WKB wave function. We also show that using the multi-differentials obtained by the topological recursion from the harmonic oscillator curve, one generates naturally the so-called Poincar\'e polynomials associated with the orbifolds of the metric ribbon graphs.


title: Twisting perturbed parafermions

abstract: The near-collinear expansion of scattering amplitudes in maximally supersymmetric Yang-Mills theory at strong coupling is governed by the dynamics of stings propagating on the five sphere. The pentagon transitions in the operator product expansion which systematize the series get reformulated in terms of matrix elements of branch-point twist operators in the two-dimensional O(6) nonlinear sigma model. The facts that the latter is an asymptotically free field theory and that there exists no local realization of twist fields prevents one from explicit calculation of their scaling dimensions and operator product expansion coefficients. This complication is bypassed making use of the equivalence of the sigma model to the infinite-level limit of WZNW models perturbed by current-current interactions, such that one can use conformal symmetry and conformal perturbation theory for systematic calculations. Presently, to set up the formalism, we consider the O(3) sigma model which is reformulated as perturbed parafermions.


title: Accelerated Evaluation of Automated Vehicles Using Piecewise Mixture Models

abstract: The process to certify highly Automated Vehicles has not yet been defined by any country in the world. Currently, companies test Automated Vehicles on public roads, which is time-consuming and inefficient. We proposed the Accelerated Evaluation concept, which uses a modified statistics of the surrounding vehicles and the Importance Sampling theory to reduce the evaluation time by several orders of magnitude, while ensuring the evaluation results are statistically accurate. In this paper, we further improve the accelerated evaluation concept by using Piecewise Mixture Distribution models, instead of Single Parametric Distribution models. We developed and applied this idea to forward collision control system reacting to vehicles making cut-in lane changes. The behavior of the cut-in vehicles was modeled based on more than 403,581 lane changes collected by the University of Michigan Safety Pilot Model Deployment Program. Simulation results confirm that the accuracy and efficiency of the Piecewise Mixture Distribution method outperformed single parametric distribution methods in accuracy and efficiency, and accelerated the evaluation process by almost four orders of magnitude.


title: Prototypal Analysis and Prototypal Regression

abstract: Prototypal analysis is introduced to overcome two shortcomings of archetypal analysis: its sensitivity to outliers and its non-locality, which reduces its applicability as a learning tool. Same as archetypal analysis, prototypal analysis finds prototypes through convex combination of the data points and approximates the data through convex combination of the archetypes, but it adds a penalty for using prototypes distant from the data points for their reconstruction. Prototypal analysis can be extended---via kernel embedding---to probability distributions, since the convexity of the prototypes makes them interpretable as mixtures. Finally, prototypal regression is developed, a robust supervised procedure which allows the use of distributions as either features or labels.


title: On the critical branching random walk III: the critical dimension

abstract: In this paper, we study the critical branching random walk in the critical dimension, $Z^4$. We provide the asymptotics of the probability of visiting a fixed finite subset and the range of the critical branching random walk conditioned on the total number of offsprings. We also prove that conditioned on visiting, the first visiting point converges in distribution.


title: Feature Selection based on PCA and PSO for Multimodal Medical Image Fusion using DTCWT

abstract: Multimodal medical image fusion helps to increase efficiency in medical diagnosis. This paper presents multimodal medical image fusion by selecting relevant features using Principle Component Analysis (PCA) and Particle Swarm Optimization techniques (PSO). DTCWT is used for decomposition of the images into low and high frequency coefficients. Fusion rules such as combination of minimum, maximum and simple averaging are applied to approximate and detailed coefficients. The fused image is reconstructed by inverse DTCWT. Performance metrics are evaluated and it shows that DTCWT-PCA performs better than DTCWT-PSO in terms of Structural Similarity Index Measure (SSIM) and Cross Correlation (CC). Computation time and feature vector size is reduced in DTCWT-PCA compared to DTCWT-PSO for feature selection which proves robustness and storage capacity.


title: High baryon densities achievable in the fragmentation regions at RHIC and LHC

abstract: We use the McLerran-Venugopalan model of the glasma energy-momentum tensor to compute the rapidity loss and excitation energy of the colliding nuclei in the fragmentation regions followed by a space-time picture to obtain their energy and baryon densities. At the top RHIC energy we find baryon densities up to 3 baryons/fm$^3$, which is 20 times that of atomic nuclei. Assuming the formation of quark-gluon plasma, we find initial temperatures of 200 to 300 MeV and baryon chemical potentials of order 1 GeV. Assuming a roughly adiabatic expansion it would imply trajectories in the $T-\mu$ plane which would straddle a possible critical point.


title: A parallel approach to bi-objective integer programming

abstract: To obtain a better understanding of the trade-offs between various objectives, Bi-Objective Integer Programming (BOIP) algorithms calculate the set of all non-dominated vectors and present these as the solution to a BOIP problem. Historically, these algorithms have been compared in terms of the number of single-objective IPs solved and total CPU time taken to produce the solution to a problem. This is equitable, as researchers can often have access to widely differing amounts of computing power. However, the real world has recently seen a large uptake of multi-core processors in computers, laptops, tablets and even mobile phones. With this in mind, we look at how to best utilise parallel processing to improve the elapsed time of optimisation algorithms. We present two methods of parallelising the recursive algorithm presented by Ozlen, Burton and MacRae. Both new methods utilise two threads and improve running times. One of the new methods, the Meeting algorithm, halves running time to achieve near-perfect parallelisation. The results are compared with the efficiency of parallelisation within the commercial IP solver IBM ILOG CPLEX, and the new methods are both shown to perform better.


title: Sparse Optimization for Robust and Efficient Loop Closing

abstract: It is essential for a robot to be able to detect revisits or loop closures for long-term visual navigation.A key insight explored in this work is that the loop-closing event inherently occurs sparsely, that is, the image currently being taken matches with only a small subset (if any) of previous images. Based on this observation, we formulate the problem of loop-closure detection as a sparse, convex $\ell_1$-minimization problem. By leveraging fast convex optimization techniques, we are able to efficiently find loop closures, thus enabling real-time robot navigation. This novel formulation requires no offline dictionary learning, as required by most existing approaches, and thus allows online incremental operation. Our approach ensures a unique hypothesis by choosing only a single globally optimal match when making a loop-closure decision. Furthermore, the proposed formulation enjoys a flexible representation with no restriction imposed on how images should be represented, while requiring only that the representations are "close" to each other when the corresponding images are visually similar. The proposed algorithm is validated extensively using real-world datasets.


title: Hardy-Littlewood inequalities on compact quantum groups of Kac type

abstract: The Hardy-Littlewood inequality on $\mathbb{T}$ compares the $L^p$-norm of a function with a weighted $\ell^p$-norm of its Fourier coefficients. The approach has recently been studied for compact homogeneous spaces and we study a natural analogue in the framework of compact quantum groups. Especially, in the case of the reduced group $C^*$-algebras and free quantum groups, we establish explicit $L^p-\ell^p$ inequalities through inherent information of underlying quantum group, such as growth rate and rapid decay property. Moreover, we show sharpness of the inequalities in a large class, including $C(G)$ with compact Lie group, $C_r^*(G)$ with polynomially growing discrete group and free quantum groups $O_N^+$, $S_N^+$.


title: One-step Estimation of Networked Population Size with Anonymity Using Respondent-Driven Capture-Recapture and Hashing

abstract: Estimates of population size for hidden and hard-to-reach individuals are of particular interest to health officials when health problems are concentrated in such populations. Efforts to derive these estimates are often frustrated by a range of factors including social stigma or an association with illegal activities that ordinarily preclude conventional survey strategies. This paper builds on and extends prior work that proposed a method to meet these challenges. Here we describe a rigorous formalization of a one-step, network-based population estimation procedure that can be employed under conditions of anonymity. The estimation procedure is designed to be implemented alongside currently accepted strategies for research with hidden populations. Simulation experiments are described that test the efficacy of the method across a range of implementation conditions and hidden population sizes. The results of these experiments show that reliable population estimates can be derived for hidden, networked population as large as 12,500 and perhaps larger for one family of random graphs. As such, the method shows potential for cost-effective implementation health and disease surveillance officials concerned with hidden populations. Limitations and future work are discussed in the concluding section.


title: Validating and describing linked data portals using shapes

abstract: Linked data portals need to be able to advertise and describe the structure of their content. A sufficiently expressive and intuitive schema language will allow portals to communicate these structures. Validation tools will aid in the publication and maintenance of linked data and increase their quality. Two schema language proposals have recently emerged for describing the structures of RDF graphs: Shape Expressions (ShEx) and Shapes Constraint Language (SHACL). In this paper we describe how these formalisms can be used in the development of a linked data portal to describe and validate its contents. As a use case, we specify a data model inspired by the WebIndex data model, a medium size linked data portal, using both ShEx and SHACL, and we propose a benchmark that can generate compliant test data structures of any size. We then perform some preliminary experiments showing performance of one validation engine based on ShEx.


title: Generic Cospark of a Matrix Can Be Computed in Polynomial Time

abstract: The cospark of a matrix is the cardinality of the sparsest vector in the column space of the matrix. Computing the cospark of a matrix is well known to be an NP hard problem. Given the sparsity pattern (i.e., the locations of the non-zero entries) of a matrix, if the non-zero entries are drawn from independently distributed continuous probability distributions, we prove that the cospark of the matrix equals, with probability one, to a particular number termed the generic cospark of the matrix. The generic cospark also equals to the maximum cospark of matrices consistent with the given sparsity pattern. We prove that the generic cospark of a matrix can be computed in polynomial time, and offer an algorithm that achieves this.


title: Nonstandard second-order formulation of the LWR model

abstract: We present a second-order formulation of the LWR model based on Phillips' model (Phillips, 1979); but the model is nonstandard with a hyperreal infinitesimal relaxation time. Since the original Phillips model is unstable with three different definitions of stability in both Eulerian and Lagrangian coordinates, we cannot use traditional methods to prove the equivalence between the second-order model, which can be considered the zero-relaxation limit of Phillips' model, and the LWR model, which is the equilibrium counterpart of Phillips' model. Instead, we resort to a nonstandard method based on the equivalence relationship between second-order continuum and car-following models established in (Jin, 2016) and prove that the nonstandard model and the LWR model are equivalent, since they have the same anisotropic car-following model and stability property. We further derive conditions for the nonstandard model to be forward-traveling and collision-free, prove that the collision-free condition is consistent with but more general than the CFL condition (Courant et al., 1928), and demonstrate that only anisotropic and symplectic Euler discretization methods lead to physically meaningful solutions. We numerically solve the lead-vehicle problem and show that the nonstandard second-order model has the same shock and rarefaction wave solutions as the LWR model for both Greenshields and triangular fundamental diagrams; for a non-concave fundamental diagram we show that the collision-free condition, but not the CFL condition, yields physically meaningful results. Finally we present a correction method to eliminate negative speeds and collisions in general second-order models, and verify the method with a numerical example.


title: Lower Bounds on the Number of Writing Operations by ILIFC with Inversion Cells

abstract: Index-less Indexed Flash Code (ILIFC) is a coding scheme for flash memories, in which one bit of a data sequence is stored in a slice consisting of several cells but the index of the bit is stored implicitly. Although several modified ILIFC schemes have been proposed, in this research we consider an ILIFC with inversion cells(I-ILIFC). The I-ILIFC reduces the total number of cell level changes at each writing request. Computer simulation is used to show that the I-ILIFC improves the average performance of the ILIFC in many cases. This paper presents our derivation of the lower bounds on the number of writing operations by I-ILIFC and shows that the worst-case performance of the I-ILIFC is better than that of the ILIFC if the code length is sufficiently large. Additionally, we consider the tight lower bounds thereon. The results show that the threshold of the code length that determines whether the I-ILIFC improves the worst-case performance of the ILIFC is smaller than that in the first lower bounds.


title: Transfinite Version of Welter's Game

abstract: We study the transfinite version of Welter's Game, a combinatorial game played on a belt divided into squares numbered with general ordinal. In particular, we give a straight-forward solution for the transfinite version, based on those of the transfinite version of Nim and the original version of Welter's Game.


title: Factorizations and Hardy-Rellich-Type Inequalities

abstract: The principal aim of this note is to illustrate how factorizations of singular, even-order partial differential operators yield an elementary approach to classical inequalities of Hardy-Rellich-type. More precisly, introducing the two-parameter $n$-dimensional homogeneous scalar differential expressions $T_{\alpha,\beta} := - \Delta + \alpha |x|^{-2} x \cdot \nabla + \beta |x|^{-2}$, $\alpha, \beta \in \mathbb{R}$, $x \in \mathbb{R}^n \backslash \{0\}$, $n \in \mathbb{N}$, $n \geq 2$, and its formal adjoint, denoted by $T_{\alpha,\beta}^+$, we show that nonnegativity of $T_{\alpha,\beta}^+ T_{\alpha,\beta}$ on $C_0^{\infty}(\mathbb{R}^n \backslash \{0\})$ implies the fundamental inequality, \begin{align} \int_{\mathbb{R}^n} [(\Delta f)(x)]^2 \, d^n x &\geq [(n - 4) \alpha - 2 \beta] \int_{\mathbb{R}^n} |x|^{-2} |(\nabla f)(x)|^2 \, d^n x \notag \\ & \quad - \alpha (\alpha - 4) \int_{\mathbb{R}^n} |x|^{-4} |x \cdot (\nabla f)(x)|^2 \, d^n x \notag \\ & \quad + \beta [(n - 4) (\alpha - 2) - \beta] \int_{\mathbb{R}^n} |x|^{-4} |f(x)|^2 \, d^n x, \notag \end{align} for $f \in C^{\infty}_0(\mathbb{R}^n \backslash \{0\})$. A particular choice of values for $\alpha$ and $\beta$ yields known Hardy-Rellich-type inequalities, including the classical Rellich inequality and an inequality due to Schmincke. By locality, these inequalities extend to the situation where $\mathbb{R}^n$ is replaced by an arbitrary open set $\Omega \subseteq \mathbb{R}^n$ for functions $f \in C^{\infty}_0(\Omega \backslash \{0\})$. Perhaps more importantly, we will indicate that our method, in addition to being elementary, is quite flexible when it comes to a variety of generalized situations involving the inclusion of remainder terms and higher-order situations.


title: Wigner functions for gauge equivalence classes of unitary irreducible representations of noncommutative quantum mechanics

abstract: While Wigner functions forming phase space representation of quantum states is a well-known fact, their construction for noncommutative quantum mechanics (NCQM) remains relatively lesser known, in particular with respect to gauge dependencies. This paper deals with the construction of Wigner functions of NCQM for a system of 2-degrees of freedom using 2-parameter families of gauge equivalence classes of unitary irreducible representations (UIRs) of the Lie group $\g$ which has been identified as the kinematical symmetry group of NCQM in an earlier paper. This general construction of Wigner functions for NCQM, in turn, yields the special cases of Landau and symmetric gauges of NCQM.


title: Co-segmentation for Space-Time Co-located Collections

abstract: We present a co-segmentation technique for space-time co-located image collections. These prevalent collections capture various dynamic events, usually by multiple photographers, and may contain multiple co-occurring objects which are not necessarily part of the intended foreground object, resulting in ambiguities for traditional co-segmentation techniques. Thus, to disambiguate what the common foreground object is, we introduce a weakly-supervised technique, where we assume only a small seed, given in the form of a single segmented image. We take a distributed approach, where local belief models are propagated and reinforced with similar images. Our technique progressively expands the foreground and background belief models across the entire collection. The technique exploits the power of the entire set of image without building a global model, and thus successfully overcomes large variability in appearance of the common foreground object. We demonstrate that our method outperforms previous co-segmentation techniques on challenging space-time co-located collections, including dense benchmark datasets which were adapted for our novel problem setting.


title: Lateral transport properties of thermally excited magnons in yttrium iron garnet films

abstract: Spin information carried by magnons is attractive for computing technology and the development of magnon-based computing circuits is of great interest. However, magnon transport in insulators has been challenging, different from the clear physical picture for spin transport in conductors. Here we investigate the lateral transport properties of thermally excited magnons in yttrium iron garnet (YIG), a model magnetic insulator. Polarity reversals of detected spins in non-local geometry devices have been experimentally observed and are strongly dependent on temperature, YIG film thickness, and injector-detector separation distance. A competing two-channel transport model for thermally excited magnons is proposed, which is qualitatively consistent with the spin signal behavior. In addition to the fundamental significance for thermal magnon transport, our work furthers the development of magnonics by creating an easily accessible magnon source with controllable transport


title: NLIE for the Sausage model

abstract: The sausage model, first proposed by Fateev, Onofri, and Zamolodchikov, is a deformation of the O(3) sigma model preserving integrability. The target space is deformed from the sphere to "sausage" shape by a deformation parameter \nu. This model is defined by a factorizable S-matrix which is obtained by deforming that of the O(3) sigma model by a parameter \lambda. Clues for the deformed sigma model are provided by various UV and IR information through the thermodynamic Bethe ansatz (TBA) analysis based on the S-matrix. Application of TBA to the sausage model is, however, limited to the case of 1/\lambda integer where the coupled integral equations can be truncated to a finite number. In this paper, we propose a finite set of nonliear integral equations (NLIEs), which are applicable to generic value of \lambda. Our derivation is based on T-Q relations extracted from the truncated TBA equations. For consistency check, we compute next-leading order corrections of the vacuum energy and extract the S-matrix information in the IR limit. We also solved the NLIE both analytically and numerically in the UV limit to get the effective central charge and compared with that of the zero-mode dynamics to obtain exact relation between \nu and \lambda. This paper is a tribute to the memory of Prof. Petr Kulish.


title: The variational-relaxation algorithm for finding quantum bound states

abstract: I describe a simple algorithm for numerically finding the ground state and low-lying excited states of a quantum system. The algorithm is an adaptation of the relaxation method for solving Poisson's equation, and is fundamentally based on the variational principle. It is especially useful for two-dimensional systems with nonseparable potentials, for which simpler techniques are inapplicable yet the computation time is minimal.


title: Interior Eigensolver for Sparse Hermitian Definite Matrices Based on Zolotarev's Functions

abstract: This paper proposes an efficient method for computing selected generalized eigenpairs of a sparse Hermitian definite matrix pencil $(A,B)$. Based on Zolotarev's best rational function approximations of the signum function and conformal mapping techniques, we construct the best rational function approximation of a rectangular function supported on an arbitrary interval via function compositions with partial fraction representations. This new best rational function approximation can be applied to construct spectrum filters of $(A,B)$ with a smaller number of poles than a direct construction without function compositions. Combining fast direct solvers and the shift-invariant generalized minimal residual method, a hybrid fast algorithm is proposed to apply spectral filters efficiently. Compared to the state-of-the-art algorithm FEAST, the proposed rational function approximation is more efficient when sparse matrix factorizations are required to solve multi-shift linear systems in the eigensolver, since the smaller number of matrix factorizations is needed in our method. The efficiency and stability of the proposed method are demonstrated by numerical examples from computational chemistry.


title: Deep Reinforcement Learning for Visual Object Tracking in Videos

abstract: In this paper we introduce a fully end-to-end approach for visual tracking in videos that learns to predict the bounding box locations of a target object at every frame. An important insight is that the tracking problem can be considered as a sequential decision-making process and historical semantics encode highly relevant information for future decisions. Based on this intuition, we formulate our model as a recurrent convolutional neural network agent that interacts with a video overtime, and our model can be trained with reinforcement learning (RL) algorithms to learn good tracking policies that pay attention to continuous, inter-frame correlation and maximize tracking performance in the long run. The proposed tracking algorithm achieves state-of-the-art performance in an existing tracking benchmark and operates at frame-rates faster than real-time. To the best of our knowledge, our tracker is the first neural-network tracker that combines convolutional and recurrent networks with RL algorithms.


title: Arithmetic degrees for dynamical systems over function fields of characteristic zero

abstract: We study arithmetic degree of a dominant rational self-map on a smooth projective variety over a function field of characteristic zero. We see that the notion of arithmetic degree and some related problems over function fields are interpreted into geometric ones. We give another proof of the theorem that the arithmetic degree at any point is smaller than or equal to the dynamical degree. We give a sufficient condition for an arithmetic degree to coincide with the dynamical degree, and prove that any self-map has so many points whose arithmetic degrees are equal to the dynamical degree. We study dominant rational self-maps on projective spaces in detail.


title: Dirac CP phase in the neutrino mixing matrix and the Froggatt-Nielsen mechanism with ${\rm \bf det} \bf [M_\nu]=0$

abstract: We discuss the Dirac CP violating phase $\delta_{CP}$ in the Froggatt-Nielsen model for a neutrino mass matrix $M_\nu$ imposing a condition ${\rm det} [M_\nu]=0$. This additional condition restricts the CP violating phase $\delta_{CP}$ drastically. We find that the phase $\delta_{CP}$ is predicted in the region of $\pm (0.4- 2.9)$ radian, which is consistent with the recent T2K and NO$\nu$A data. There is a remarkable correlation between $\delta_{CP}$ and $\sin^2\theta_{23}$. The phase $\delta_{CP}$ converges on $\sim \pm \pi/2$ if $\sin^2\theta_{23}$ is larger than $0.5$. Thus, accurate measurements of $\sin^2\theta_{23}$ are important for a test of our model. The effective mass $m_{ee}$ for the neutrinoless double beta decay is predicted in the rage of $3.3-4.0$ meV.


title: Deep Submodular Functions

abstract: We start with an overview of a class of submodular functions called SCMMs (sums of concave composed with non-negative modular functions plus a final arbitrary modular). We then define a new class of submodular functions we call {\em deep submodular functions} or DSFs. We show that DSFs are a flexible parametric family of submodular functions that share many of the properties and advantages of deep neural networks (DNNs). DSFs can be motivated by considering a hierarchy of descriptive concepts over ground elements and where one wishes to allow submodular interaction throughout this hierarchy. Results in this paper show that DSFs constitute a strictly larger class of submodular functions than SCMMs. We show that, for any integer $k>0$, there are $k$-layer DSFs that cannot be represented by a $k'$-layer DSF for any $k'<k$. This implies that, like DNNs, there is a utility to depth, but unlike DNNs, the family of DSFs strictly increase with depth. Despite this, we show (using a "backpropagation" like method) that DSFs, even with arbitrarily large $k$, do not comprise all submodular functions. In offering the above results, we also define the notion of an antitone superdifferential of a concave function and show how this relates to submodular functions (in general), DSFs (in particular), negative second-order partial derivatives, continuous submodularity, and concave extensions. To further motivate our analysis, we provide various special case results from matroid theory, comparing DSFs with forms of matroid rank, in particular the laminar matroid. Lastly, we discuss strategies to learn DSFs, and define the classes of deep supermodular functions, deep difference of submodular functions, and deep multivariate submodular functions, and discuss where these can be useful in applications.


title: Harmonic Eisenstein Series of Weight One

abstract: In this short note, we will construct a harmonic Eisenstein series of weight one, whose image under the $\xi$-operator is a weight one Eisenstein series studied by Hecke.


title: Compositeness Of Quasibound States from Effective Field Theory

abstract: The internal structure of exotic hadron candidates are studied from the viewpoint of the compositeness of near-threshold states. We focus on the Weinberg's weak-binding relation between the experimental observables and the compositeness of the state. First we extend the relation for quasibound states within the framework of the effective field theory. Next, considering the correction term, we develop a systematic method to estimate the uncertainty of compositeness in the weak-binding relation. Finally, applying the extended relation and the error evaluation method, we conclude that the structure of Lambda(1405) is dominated by the Kbar N component.


title: Dimensionality-strain phase diagram of strontium iridates

abstract: The competition between spin-orbit coupling, bandwidth ($W$) and electron-electron interaction ($U$) makes iridates highly susceptible to small external perturbations, which can trigger the onset of novel types of electronic and magnetic states. Here we employ {\em first principles} calculations based on density functional theory and on the constrained random phase approximation to study how dimensionality and strain affect the strength of $U$ and $W$ in (SrIrO$_3$)$_m$/(SrTiO$_3$) superlattices. The result is a phase diagram explaining two different types of controllable magnetic and electronic transitions, spin-flop and insulator-to-metal, connected with the disruption of the $J_{eff}=1/2$ state which cannnot be understood within a simplified local picture.


title: Convex Hulls for the Unit Commitment Polytope

abstract: In this paper, we consider the polyhedral structure of the unit commitment polytope. In particular, we provide the convex hull results for the problem under the following different settings: 1) the convex hulls for the integrated minimum-up/-down time and ramping polytope under the general $T$ time period setting in which the ramping rate equals to the gap between the generation upper and lower bounds and equals to half of the gap between the generation upper and lower bounds, respectively, 2) the convex hull for the integrated minimum-up/-down time and ramping-up polytope for the problem under the general $T$ time period setting, and 3) the convex hull for the integrated minimum-up/-down time and ramping-down polytope for the problem under the general $T$ time period setting.


title: Thermoelectric properties of $\beta$-As, Sb and Bi monolayers

abstract: Monolayer semiconductors of group-VA elements (As, Sb, Bi) with graphenelike buckled structure offer a potential to achieve nanoscale electronic, optoelectronic and thermoelectric devices. Motivated by recently-fabricated Sb monolayer, we systematically investigate the thermoelectric properties of $\beta$-As, Sb and Bi monolayers by combining the first-principles calculations and semiclassical Boltzmann transport theory. The generalized gradient approximation (GGA) plus spin-orbit coupling (SOC) is adopted for the electron part, and GGA is employed for the phonon part. It is found that SOC has important influences on their electronic structures, especially for Bi monolayer, which can induce observable SOC effects on electronic transport coefficients. More specifically, SOC not only has detrimental influences on electronic transport coefficients, but also produces enhanced effects. The calculated lattice thermal conductivity decreases gradually from As to Bi monolayer, and the corresponding room-temperature sheet thermal conductance is 161.10 $\mathrm{W K^{-1}}$, 46.62 $\mathrm{W K^{-1}}$ and 16.02 $\mathrm{W K^{-1}}$, which can be converted into common lattice thermal conductivity by dividing by the thickness of 2D material. The sheet thermal conductance of Bi monolayer is lower than one of other 2D materials, such as semiconducting transition-metal dichalcogenide monolayers and orthorhombic group IV-VI monolayers. A series of scattering time is employed to estimate the thermoelectric figure of merit $ZT$. It is found that the n-type doping has more excellent thermoelectric properties than p-type doping for As and Bi monolayer, while the comparative $ZT$ between n- and p-type doping is observed in Bi monolayer. These results can stimulate further experimental works to open the new field for thermoelectric devices based on monolayer of group-VA elements.


title: Design of racetrack ring resonator based dielectric laser accelerators

abstract: In this paper, we propose a novel design of dielectric laser-driven accelerator (DLA) utilizing evanescent electric field of racetrack ring resonator structures. Driven by laser light with the correctly designed optical phase window, sustained acceleration of electrons with controlled deflection is shown. Based on this design, we calculate an acceleration from 30 keV to 148.312 keV in 104.655 {\mu}m using a cascaded 11-stage racetrack ring resonators. This new idea poses a solution for on-chip integration of many DLA stages, while maintains high average accelerating gradients, providing a potential practical realization for "accelerator on a chip".


title: Variable selection for clustering with Gaussian mixture models: state of the art

abstract: The mixture models have become widely used in clustering, given its probabilistic framework in which its based, however, for modern databases that are characterized by their large size, these models behave disappointingly in setting out the model, making essential the selection of relevant variables for this type of clustering. After recalling the basics of clustering based on a model, this article will examine the variable selection methods for model-based clustering, as well as presenting opportunities for improvement of these methods.


title: Sparse phase retrieval of one-dimensional signals by Prony's method

abstract: In this paper, we show that sparse signals f representable as a linear combination of a finite number N of spikes at arbitrary real locations or as a finite linear combination of B-splines of order m with arbitrary real knots can be almost surely recovered from O(N^2) Fourier intensity measurements up to trivial ambiguities. The constructive proof consists of two steps, where in the first step the Prony method is applied to recover all parameters of the autocorrelation function and in the second step the parameters of f are derived. Moreover, we present an algorithm to evaluate f from its Fourier intensities and illustrate it at different numerical examples.


title: Singularities of the area preserving curve shortening flow with a Neumann free boundary condition

abstract: We consider the area preserving curve shortening flow with Neumann free boundary conditions outside of a convex domain or at a straight line. We give a criterion on initial curves that guarantees the appearance of a singularity in finite time. We prove that the singularity is of type II. Furthermore, if these initial curves are convex, then an appropriate rescaling at the finite maximal time of existence yields a grim reaper or half a grim reaper as limit flow. We construct examples of initial curves satisfying the mentioned criterion.


title: On Supernilpotent Algebras

abstract: We establish a characterization of supernilpotent Mal'cev algebras which generalizes the affine structure of abelian Mal'cev algebras and the recent characterization of 3-supernilpotent Mal'cev algebras. We then show that for varieties in which the two-generated free algebra is finite: (1) neutrality of the higher commutators is equivalent to congruence meet-semidistributivity, and (2) the class of varieties which interpret a Mal'cev term in every supernilpotent algebra is equivalent to the existence of a weak difference term. We then establish properties of the higher commutator in the aforementioned second class of varieties.


title: Shifted tableau switchings and shifted Littlewood-Richardson coefficients

abstract: We provide two shifted analogues of the tableau switching process due to Benkart, Sottile, and Stroomer, the shifted tableau switching process and the modified shifted tableau switching process. They are performed by applying a sequence of specially contrived elementary transformations called {\em switches} and turn out to have some spectacular properties. For instance, the maps induced from these algorithms are involutive and behave very nicely with respect to shifted Young tableaux whose reading words satisfy the lattice property. As an application, we give combinatorial interpretations of Schur $P$- and $Q$-function identities. We also demonstrate the relationship between the shifted tableau switching process and the shifted $J$-operation due to Worley.


title: A Hybrid Approach for Secured Optimal Power Flow and Voltage Stability with TCSC Placement

abstract: This paper proposes a hybrid technique for secured optimal power flow coupled with enhancing voltage stability with FACTS device installation. The hybrid approach of Improved Gravitational Search algorithm (IGSA) and Firefly algorithm (FA) performance is analyzed by optimally placing TCSC controller. The algorithm is implemented in MATLAB working platform and the power flow security and voltage stability is evaluated with IEEE 30 bus transmission systems. The optimal results generated are compared with those available in literature and the superior performance of algorithm is depicted as minimum generation cost, reduced real power losses along with sustaining voltage stability.


title: Characterizations of input-to-state stability for infinite-dimensional systems

abstract: We prove characterizations of input-to-state stability (ISS) for a large class of infinite-dimensional control systems, including some classes of evolution equations over Banach spaces, time-delay systems, ordinary differential equations (ODE), switched systems. These characterizations generalize well-known criteria of ISS, proved by Sontag and Wang for ODE systems. For the special case of differential equations in Banach spaces we prove even broader criteria for ISS and apply these results to show that (under some mild restrictions) the existence of a non-coercive ISS Lyapunov functions implies ISS. We introduce the new notion of strong ISS which is equivalent to ISS in the ODE case, but which is strictly weaker than ISS in the infinite-dimensional setting and prove several criteria for the sISS property. At the same time, we show by means of counterexamples, that many characterizations, which are valid in the ODE case, are not true for general infinite-dimensional systems.


title: A tale of clusters: No resolvable periodicity in the terrestrial impact cratering record

abstract: Rampino & Caldeira (2015) carry out a circular spectral analysis (CSA) of the terrestrial impact cratering record over the past 260 million years (Ma), and suggest a ~26 Ma periodicity of impact events. For some of the impacts in that analysis, new accurate and high-precision ("robust"; 2SE<2%) 40Ar-39Ar ages have recently been published, resulting in significant age shifts. In a CSA of the updated impact age list, the periodicity is strongly reduced. In a CSA of a list containing only impacts with robust ages, we find no significant periodicity for the last 500 Ma. We show that if we relax the assumption of a fully periodic impact record, assuming it to be a mix of a periodic and a random component instead, we should have found a periodic component if it contributes more than ~80% of the impacts in the last 260 Ma. The difference between our CSA and the one by Rampino & Caldeira (2015) originates in a subset of "clustered" impacts (i.e., with overlapping ages). The ~26 Ma periodicity seemingly carried by these clusters alone is strongly significant if tested against a random distribution of ages, but this significance disappears if it is tested against a distribution containing (randomly-spaced) clusters. The presence of a few impact age clusters (e.g., from asteroid break-up events) in an otherwise random impact record can thus give rise to false periodicity peaks in a CSA. There is currently no evidence for periodicity in the impact record.


title: CommAI: Evaluating the first steps towards a useful general AI

abstract: With machine learning successfully applied to new daunting problems almost every day, general AI starts looking like an attainable goal. However, most current research focuses instead on important but narrow applications, such as image classification or machine translation. We believe this to be largely due to the lack of objective ways to measure progress towards broad machine intelligence. In order to fill this gap, we propose here a set of concrete desiderata for general AI, together with a platform to test machines on how well they satisfy such desiderata, while keeping all further complexities to a minimum.


title: Microlocal analysis and beyond

abstract: We shall explain how the idea of microlocal analysis of the seventies has been reformulated in the framework of sheaf theory in the eighties and then applied to various branches of mathematics, such as linear partial differential equations or symplectic topology.


title: Maximum and minimum of support functions

abstract: For given continuous functions $\gamma_{{}_{i}}: S^{n}\to \mathbb{R}_{+}$ (where $i=1, 2$), the functions $\gamma_{{}_{max}}$ and $\gamma_{{}_{min}}$ can be defined as natural way. In this paper, we show that the Wulff shape associated to $\gamma_{{}_{max}}$ is the convex hull of the union of Wulff shapes associated to $\gamma_{{}_1}$ and $\gamma_{{}_2}$ , if $\gamma_{{}_1}$ and $\gamma_{{}_2}$ are convex integrands. And, the Wulff shape associated to $\gamma_{{}_{min}}$ is the intersection of Wulff shapes associated to $\gamma_{{}_1}$ and $\gamma_{{}_2}$. Moreover, relationships between their dual Wulff shapes are given.


title: Apparent quasar disc sizes in the "bird's nest" paradigm

abstract: Quasar microlensing effects make it possible to measure the accretion disc sizes around distant supermassive black holes that are still well beyond the spatial resolution of contemporary instrumentation. The sizes measured with this technique appear inconsistent with the standard accretion disc model. Not only are the measured accretion disc sizes larger, but their dependence on wavelength is in most cases completely different from the predictions of the standard model. We suggest that these discrepancies may arise not from non-standard accretion disc structure or systematic errors, as it was proposed before, but rather from scattering and reprocession of the radiation of the disc. In particular, the matter falling from the gaseous torus and presumably feeding the accretion disc may at certain distances become ionized and produce an extended halo that is free from colour gradients. A simple analytical model is proposed assuming that a geometrically thick translucent inflow acts as a scattering mirror changing the apparent spatial properties of the disc. This inflow may be also identified with the broad line region or its inner parts. Such a model is able to explain the basic properties of the apparent disc sizes, primarily their large values and their shallow dependence on wavelength. The only condition required is to scatter significant portion of the luminosity of the disc. This can easily be fulfilled if the scattering inflow has large geometrical thickness and clumpy structure.


title: Transverse Anderson localization of surface plasmon polaritons

abstract: We investigate the effect of disorder on the propagation of surface plasmon polaritons in arrays of evanescently coupled dielectric loaded surface plasmon polariton waveguides. Diagonal disorder is implemented by randomly varying the heights of the waveguides. Real-space as well as Fourier-space images of the surface plasmon polariton intensity distribution in the waveguide arrays are recorded by leakage radiation microscopy. With these techniques we experimentally demonstrate the transverse localization of surface plasmon polaritons with increasing disorder.


title: Simulating low dimensional QCD with Lefschetz thimbles

abstract: Non-perturbative lattice QCD calculations at non vanishing baryon number density are hampered by the QCD sign problem. The path integral, that in lattice QCD is calculated numerically, becomes highly oscillating. One possible solution is the Lefschetz thimble approach. It requires a deformation of the original integration domain into a manifold embedded in complex space. For properly chosen integration manifolds ("thimbles") the sign problem is drastically alleviated. For some bosonic and fermionic models this approach has been shown to work. Here we apply the thimble disretization to (0+1)-dimensional QCD with standard staggerd quarks and disscuss issues that may arrise in higher dimensions.


title: Gustafson-Rakha-Type Elliptic Hypergeometric Series

abstract: We prove a multivariable elliptic extension of Jackson's summation formula conjectured by Spiridonov. The trigonometric limit case of this result is due to Gustafson and Rakha. As applications, we obtain two further multivariable elliptic Jackson summations and two multivariable elliptic Bailey transformations. The latter four results are all new even in the trigonometric case.


title: Isospin-0 $\pi\pi$ scattering from twisted mass lattice QCD

abstract: We present results for the isospin-0 $\pi\pi$ s-wave scattering length calculated in twisted mass lattice QCD. We use three $N_f = 2$ ensembles with unitary pion mass at its physical value, 240~MeV and 330~MeV respectively. We also use a large set of $N_f = 2 + 1 +1$ ensembles with unitary pion masses varying in the range of 230~MeV - 510~MeV at three different values of the lattice spacing. A mixed action approach with the Osterwalder-Seiler action in the valence sector is adopted to circumvent the complications arising from isospin symmetry breaking of the twisted mass quark action. Due to the relatively large lattice artefacts in the $N_f = 2 + 1 +1$ ensembles, we do not present the scattering lengths for these ensembles. Instead, taking the advantage of the many different pion masses of these ensembles, we qualitatively discuss the pion mass dependence of the scattering properties of this channel based on the results from the $N_f = 2 + 1 +1$ ensembles. The scattering length is computed for the $N_f = 2$ ensembles and the chiral extrapolation is performed. At the physical pion mass, our result $M_\pi a^\mathrm{I=0}_0 = 0.198(9)(6)$ agrees reasonably well with various experimental measurements and theoretical predictions.


title: On a Fractional Oscillator Equation with Natural Boundary Conditions

abstract: We prove existence of solutions for a nonlinear fractional oscillator equation with both left Riemann-Liouville and right Caputo fractional derivatives subject to natural boundary conditions. The proof is based on a transformation of the problem into an equivalent lower order fractional boundary value problem followed by the use of an upper and lower solutions method. To succeed with such approach, we first prove a result on the monotonicity of the right Caputo derivative.


title: Twin-beam real-time position estimation of micro-objects in 3D

abstract: Various optical methods for measuring positions of micro-objects in 3D have been reported in the literature. Nevertheless, majority of them are not suitable for real-time operation, which is needed, for example, for feedback position control. In this paper, we present a method for real-time estimation of the position of micro-objects in 3D; the method is based on twin-beam illumination and it requires only a very simple hardware setup whose essential part is a standard image sensor without any lens. Performance of the proposed method is tested during a micro-manipulation task in which the estimated position served as a feedback for the controller. The experiments show that the estimate is accurate to within ~3 um in the lateral position and ~7 um in the axial distance with the refresh rate of 10 Hz. Although the experiments are done using spherical objects, the presented method could be modified to handle non-spherical objects as well.


title: Enhancing Photon Correlations through Plasmonic Strong Coupling

abstract: There is an increasing scientific and technological interest on the design and implementation of nanoscale sources of quantum light. Here, we investigate the quantum statistics of the light scattered from a plasmonic nanocavity coupled to a mesoscopic ensemble of emitters under low coherent pumping. We present an analytical description of the intensity correlations taking place in these systems, and unveil the fingerprint of plasmon-exciton-polaritons in them. Our findings reveal that plasmonic cavities are able to retain and enhance excitonic nonlinearities even when the number of emitters is large. This makes plasmonic strong coupling a promising route for generating nonclassical light beyond the single emitter level.