title: A Parsimonious Personalized Dose Finding Model via Dimension Reduction

abstract: Learning an individualized dose rule in personalized medicine is a challenging statistical problem. Existing methods often suffer from the curse of dimensionality, especially when the decision function is estimated nonparametrically. To tackle this problem, we propose a dimension reduction framework that effectively reduces the estimation to a lower-dimensional subspace of the covariates. We exploit that the individualized dose rule can be defined in a subspace spanned by a few linear combinations of the covariates, leading to a more parsimonious model. The proposed framework does not require the inverse probability of the propensity score under observational studies due to a direct maximization of the value function. This distinguishes us from the outcome weighted learning framework, which also solves decision rules directly. Under the same framework, we further propose a pseudo-direct learning approach that focuses more on estimating the dimensionality-reduced subspace of the treatment outcome. Parameters in both approaches can be estimated efficiently using an orthogonality constrained optimization algorithm on the Stiefel manifold. Under mild regularity assumptions, the results on the asymptotic normality of the proposed estimators are established, respectively. We also derive the consistency and convergence rate for the value function under the estimated optimal dose rule. We evaluate the performance of the proposed approaches through extensive simulation studies and a warfarin pharmacogenetic dataset.


title: Attack on the Edon-K Key Encapsulation Mechanism

abstract: The key encapsulation mechanism Edon-K was proposed in response to the call for post-quantum cryptography standardization issued by the National Institute of Standards and Technologies (NIST). This scheme is inspired by the McEliece scheme but uses another family of codes defined over $\mathbb{F}_{2^{128}}$ instead of $\mathbb{F}_2$ and is not based on the Hamming metric. It allows significantly shorter public keys than the McEliece scheme. In this paper, we give a polynomial time algorithm that recovers the encapsulated secret. This attack makes the scheme insecure for the intended use. We obtain this result by observing that recovering the error in the McEliece scheme corresponding to Edon-K can be viewed as a decoding problem for the rank-metric. We show that the code used in Edon-K is in fact a super-code of a Low Rank Parity Check (LRPC) code of very small rank (1 or 2). A suitable parity-check matrix for the super-code of such low rank can be easily derived from for the public key. We then use this parity-check matrix in a decoding algorithm that was devised for LRPC codes to recover the error. Finally we explain how we decapsulate the secret once we have found the error.


title: Study of the likelihood of Alfv\'enic mode bifurcation in NSTX and predictions for ITER baseline scenarios

abstract: Rare Alfv\'enic wave transitions between fixed-frequency and chirping phases are identified in NSTX, where Alfv\'enic waves are normally observed to exhibit either chirping or avalanching responses. For those transitions, we apply a criterion [Duarte et al, Nucl. Fusion 57, 054001 (2017)] to predict the nature of fast ion redistribution in tokamaks to be in the convective or diffusive nonlinear regimes. For NSTX discharges in which the transition is not accompanied by changes in the beam deposited power or modifications in the injected radiofrequency power, it has been found that the anomalous fast ion transport is a likely mediator of the bifurcation between the fixed-frequency mode behavior and rapid chirping. For a quantitative assessment, global gyrokinetic simulations of the effects of electrostatic ion temperature gradient turbulence and trapped electron mode turbulence on chirping were pursued using the GTS code. The investigation is extended by means of predictive studies of the probable spectral behavior of Alfv\'enic eigenmodes for baseline ITER cases consisting of elmy, advanced and hybrid scenarios. It has been observed that most modes are found to be borderline between the steady and the chirping phases.


title: Ad Hoc Table Retrieval using Semantic Similarity

abstract: We introduce and address the problem of ad hoc table retrieval: answering a keyword query with a ranked list of tables. This task is not only interesting on its own account, but is also being used as a core component in many other table-based information access scenarios, such as table completion or table mining. The main novel contribution of this work is a method for performing semantic matching between queries and tables. Specifically, we (i) represent queries and tables in multiple semantic spaces (both discrete sparse and continuous dense vector representations) and (ii) introduce various similarity measures for matching those semantic representations. We consider all possible combinations of semantic representations and similarity measures and use these as features in a supervised learning model. Using a purpose-built test collection based on Wikipedia tables, we demonstrate significant and substantial improvements over a state-of-the-art baseline.


title: Effect of $\alpha$-particle irradiation on a NdFeAs(O,F) thin film

abstract: The effect of $\alpha$-particle irradiation on a NdFeAs(O,F) thin film has been investigated to determine how the introduction of defects affects basic superconducting properties, including the critical temperature $T_c$ and the upper critical field $H_{c2}$, and properties more of interest for applications, like the critical current density $J_c$ and the related pinning landscape. The irradiation-induced suppression of the film $T_c$ is significantly smaller than on a similarly damaged single crystal. Moreover $H_{c2}$ behaves differently, depending on the field orientation: for H//c the $H_{c2}$ slope monotonically increases with increasing disorder, whereas for H//ab it remains constant at low dose and it increases only when the sample is highly disordered. This suggests that a much higher damage level is necessary to drive the NdFeAs(O,F) thin film into the dirty limit. Despite the increase in the low temperature $H_{c2}$, the effects on the $J_c$(H//c) performances are moderate in the measured temperature and field ranges, with a shifting of the pinning force maximum from 4.5 T to 6 T after an irradiation of $2\times10^{15} cm^{-2}$. On the contrary, $J_c$(H//ab) is always suppressed. The analysis demonstrates that irradiation does introduce point defects acting as pinning centres proportionally to the irradiation fluence but also suppresses the effectiveness of c-axis correlated pinning present in the pristine sample. We estimate that significant performance improvements may be possible at high field or at temperatures below 10 K. The suppression of the $J_c$(H//ab) performance is not related to a decrease of the $J_c$ anisotropy as found in other superconductors. Instead it is due to the presence of point defects that decrease the efficiency of the ab-plane intrinsic pinning typical of materials with a layered structure.


title: Ultracold Rydberg Molecules

abstract: Ultracold Rydberg molecules have been extensively studied both theoretically and ex-perimentally. Here the authors review the recent realizations of various ultralong-range Rydberg molecules and macrodimers, and explore their potential for current and future applications in spectroscopy, few- and many-body interactions and quantum information processing.


title: Liquid xenon scintillation measurements and pulse shape discrimination in the LUX dark matter detector

abstract: Weakly Interacting Massive Particles (WIMPs) are a leading candidate for dark matter and are expected to produce nuclear recoil (NR) events within liquid xenon time-projection chambers. We present a measurement of the scintillation timing characteristics of liquid xenon in the LUX dark matter detector and develop a pulse shape discriminant to be used for particle identification. To accurately measure the timing characteristics, we develop a template-fitting method to reconstruct the detection times of photons. Analyzing calibration data collected during the 2013-16 LUX WIMP search, we provide a new measurement of the singlet-to-triplet scintillation ratio for electron recoils (ER) below 46~keV, and we make a first-ever measurement of the NR singlet-to-triplet ratio at recoil energies below 74~keV. We exploit the difference of the photon time spectra for NR and ER events by using a prompt fraction discrimination parameter, which is optimized using calibration data to have the least number of ER events that occur in a 50\% NR acceptance region. We then demonstrate how this discriminant can be used in conjunction with the charge-to-light discrimination to possibly improve the signal-to-noise ratio for nuclear recoils.


title: A bound for a typical differential dimension of system of linear differential equations

abstract: We prove upper and lower bounds for leading coefficient of Kolchin dimension polynomial of systems of partial linear differential equations in codimension two.


title: Inferring Information Flow in Spike-train Data Sets using a Trial-Shuffle Method

abstract: Understanding information processing in the brain requires the ability to determine the functional connectivity between the different regions of the brain. We present a method using transfer entropy to extract this flow of information between brain regions from spike-train data commonly obtained in neurological experiments. Transfer entropy is a statistical measure based in information theory that attempts to quantify the information flow from one process to another, and has been applied to find connectivity in simulated spike-train data. Due to statistical error in the estimator, inferring functional connectivity requires a method for determining significance in the transfer entropy values. We discuss the issues with numerical estimation of transfer entropy and resulting challenges in determining significance before presenting the trial-shuffle method as a viable option. The trial-shuffle method, for spike-train data that is split into multiple trials, determines significant transfer entropy values independently for each individual pair of neurons by comparing to a created baseline distribution using a rigorous statistical test. This is in contrast to either globally comparing all neuron transfer entropy values or comparing pairwise values to a single baseline value. In establishing the viability of this method by comparison to several alternative approaches in the literature, we find evidence that preserving the inter-spike-interval timing is important. We then use the trial-shuffle method to investigate information flow within a model network as we vary model parameters. This includes investigating the global flow of information within a connectivity network divided into two well-connected subnetworks, going beyond local transfer of information between pairs of neurons.


title: Tractable and Robust Modeling of Building Flexibility Using Coarse Data

abstract: Controllable building loads have the potential to increase the flexibility of power systems. A key step in developing effective and attainable load control policies is modeling the set of feasible building load profiles. In this paper, we consider buildings whose source of flexibility is their HVAC loads. We propose a data-driven method to empirically estimate a robust feasible region of the load using coarse data, that is, using only total building load and average indoor temperatures. The proposed method uses easy-to-gather coarse data and can be adapted to buildings of any type. The resulting feasible region model is robust to temperature prediction errors and is described by linear constraints. The mathematical simplicity of these constraints makes the proposed model adaptable to many power system applications, for example, economic dispatch, and optimal power flow. We validate our model using data from EnergyPlus and demonstrate its usefulness through a case study in which flexible building loads are used to balance errors in wind power forecasts.


title: The Structure and Dynamics of the Sub-parsec Scale Jet in M87 Based on 50 VLBA Observations Over 17 Years at 43 GHz

abstract: The central radio source in M87 provides the best opportunity to study jet formation because it has a large angular size for the gravitational radius of the black hole and has a bright jet that is well resolved by VLBI observations. We present intensive monitoring observations from 2007 and 2008, plus roughly annual observations that span 17 years, all made with the the Very Long Baseline Array at 43 GHz with a resolution of about 30 by 60 Rs. Our high-dynamic-range images clearly show the wide-opening-angle structure and the counter-jet. The jet and counter-jet are nearly symmetric in the inner 1.5 milli-arcseconds (mas; 0.12 pc in projection) with both being edge brightened. Both show deviations from parabolic shape in the form of an initial rapid expansion and subsequent contraction followed by further rapid expansion and, beyond the visible counter-jet, subsequent collimation. Proper motions and counter-jet/jet intensity ratios both indicate acceleration from apparent speeds of $\lesssim 0.5c$ to $\gtrsim 2c$ in the inner about 2 mas and suggest a helical flow. The jet displays a sideways shift with an approximately 8 to 10 year quasi-periodicity. The shift propagates outwards non-ballistically and significantly more slowly than the flow speed revealed by the fastest moving components. Polarization data show a systematic structure with magnetic field vectors that suggest a toroidal field close to the core.


title: CapsuleGAN: Generative Adversarial Capsule Network

abstract: We present Generative Adversarial Capsule Network (CapsuleGAN), a framework that uses capsule networks (CapsNets) instead of the standard convolutional neural networks (CNNs) as discriminators within the generative adversarial network (GAN) setting, while modeling image data. We provide guidelines for designing CapsNet discriminators and the updated GAN objective function, which incorporates the CapsNet margin loss, for training CapsuleGAN models. We show that CapsuleGAN outperforms convolutional-GAN at modeling image data distribution on MNIST and CIFAR-10 datasets, evaluated on the generative adversarial metric and at semi-supervised image classification.


title: CHANG-ES XI: Circular Polarization in the Cores of Nearby Galaxies

abstract: We detect 5 galaxies in the CHANG-ES (Continuum Halos in Nearby Galaxies -- an EVLA Survey) sample that show circular polarization (CP) at L-band in our high resolution data sets. Two of the galaxies (NGC~4388 and NGC~4845) show strong Stokes $V/I\,\equiv\,m_C\,\sim\,2$\%, two (NGC~660 and NGC~3628) have values of $m_C\sim \,0.3$\%, and NGC~3079 is a marginal detection at $m_C\sim \,0.2$\%. The two strongest $m_C$ galaxies also have the most luminous X-ray cores and the strongest internal absorption in X-rays. We have expanded on our previous Faraday conversion interpretation and analysis and provide analytical expressions for the expected $V$ signal for a general case in which the cosmic ray electron energy spectral index can take on any value. We provide examples as to how such expressions could be used to estimate magnetic field strengths and the lower energy cutoff for CR electrons. Four out of our detections are {\it resolved}, showing unique structures, including a {\it jet} in NGC~4388 and a CP `conversion disk' in NGC~4845. The conversion disk is inclined to the galactic disk but is perpendicular to a possible outflow direction. Such CP structures have never before been seen in any galaxy to our knowledge. None of the galaxy cores show linear polarization at L-band. Thus CP may provide a unique probe of physical conditions deep into radio AGNs.


title: Role of edge poloidal density asymmetry in tokamak confinement

abstract: Mass flows and radial electric field driven by edge poloidal density asymmetries can be used as a highly effective control mechanism for the edge and thus global confinement in tokamaks. The underlying physics can be demonstrated entirely within a simple magnetohydrodynamic equilibrium model, without resorting to sophisticated and usually collisionality-dependent neoclassical physics arguments. As an example, strong dependence of the low to high (LH) transition power threshold on the magnetic topology, an experimental observation still poorly understood, can be easily explained within this framework. Similar arguments indicate that the ITER fueling ports above the midplane might lead to higher input power requirements.


title: Random Relation Algebras

abstract: We develop a random model for relation algebras. We prove some preliminary results and pose questions that lay out a new direction of research.


title: Potential and spin-exchange interaction between Anderson impurities in graphene

abstract: The effective interaction between resonant magnetic Anderson impurities in graphene, mediated by conduction electrons, is studied as a function of the strength of the onsite energy level of the impurities and the amplitude of coupling to conduction electrons. The sign and character of the interaction depend on whether the impurities reside on the same or opposite sublattices. For the same (opposite) sublattice, the potential interaction is attractive (repulsive) in the weak coupling limit with $1/R^3$ dependence on the distance; the interaction reverses sign and becomes repulsive (attractive) in the strong coupling limit and displays $1/R$ behavior. The spin-exchange coupling is ferromagnetic (antiferromagnetic) at both large and small distances, but reverses sign and becomes anti-ferromagnetic (ferromagnetic) for intermediate distances. For opposite sublattices, the effective spin exchange coupling is resonantly enhanced at distances where the energy levels cross the Dirac points.


title: Controlling rare-earth ions in a nanophotonic resonator using the ac Stark shift

abstract: On-chip nanophotonic cavities will advance quantum information science and measurement because they enable efficient interaction between photons and long-lived solid-state spins, such as those associated with rare-earth ions in crystals. The enhanced photon-ion interaction creates new opportunities for all-optical control using the ac Stark shift. Toward this end, we characterize the ac Stark interaction between off-resonant optical fields and Nd$^{3+}$-ion dopants in a photonic crystal resonator fabricated from yttrium orthovanadate (YVO$_4$). Using photon echo techniques, at a detuning of 160 MHz we measure a maximum ac Stark shift of 2$\pi\times$12.3 MHz per intra-cavity photon, which is large compared to both the homogeneous linewidth ($\Gamma_h =$100 kHz) and characteristic width of isolated spectral features created through optical pumping ($\Gamma_f \approx$3 MHz). The photon-ion interaction strength in the device is sufficiently large to control the frequency and phase of the ions for quantum information processing applications. In particular, we discuss and assess the use of the cavity enhanced ac Stark shift to realize all-optical quantum memory and detection protocols. Our results establish the ac Stark shift as a powerful added control in rare-earth ion quantum technologies.


title: Matrix variate Birnbaum-Saunders distribution under elliptical models

abstract: This paper derives the elliptical matrix variate version of the well known univariate Birnbaum and Saunders distribution. A generalisation based on a matrix transformation is proposed, instead of the independent element by element representation of the Gaussian univariate version of 1969. New results on Jacobians were needed to derived the matrix variate distribution. A number of special cases are studied and some basic properties are found. Finally, an example based on real data of two populations is provided. The maximum likelihood estimates are found for a number of matrix variate generalised Birnbaum-Saunders distributions based on Kotz models. A comparison with the Gaussian kernel is also given by using a modified BIC criterion.


title: Light inflaton completing Higgs inflation

abstract: We consider the extension of the Standard Model (SM) with a light inflaton where both unitarity problem in Higgs inflation and vacuum instability problem are resolved. The linear non-minimal coupling of the inflaton to gravity leads to a significant kinetic mixing between the inflaton and the graviton such that perturbative unitarity is restored up to Planck scale. We show the correlation between unitarity scale and inflationary observables in this model and discuss how the effective Higgs inflation appears.


title: An Alternative View: When Does SGD Escape Local Minima?

abstract: Stochastic gradient descent (SGD) is widely used in machine learning. Although being commonly viewed as a fast but not accurate version of gradient descent (GD), it always finds better solutions than GD for modern neural networks. In order to understand this phenomenon, we take an alternative view that SGD is working on the convolved (thus smoothed) version of the loss function. We show that, even if the function $f$ has many bad local minima or saddle points, as long as for every point $x$, the weighted average of the gradients of its neighborhoods is one point convex with respect to the desired solution $x^*$, SGD will get close to, and then stay around $x^*$ with constant probability. More specifically, SGD will not get stuck at "sharp" local minima with small diameters, as long as the neighborhoods of these regions contain enough gradient information. The neighborhood size is controlled by step size and gradient noise. Our result identifies a set of functions that SGD provably works, which is much larger than the set of convex functions. Empirically, we observe that the loss surface of neural networks enjoys nice one point convexity properties locally, therefore our theorem helps explain why SGD works so well for neural networks.


title: Introduction to topological quantum computation with non-Abelian anyons

abstract: Topological quantum computers promise a fault tolerant means to perform quantum computation. Topological quantum computers use particles with exotic exchange statistics called non-Abelian anyons, and the simplest anyon model which allows for universal quantum computation by particle exchange or braiding alone is the Fibonacci anyon model. One classically hard problem that can be solved efficiently using quantum computation is finding the value of the Jones polynomial of knots at roots of unity. We aim to provide a pedagogical, self-contained, review of topological quantum computation with Fibonacci anyons, from the braiding statistics and matrices to the layout of such a computer and the compiling of braids to perform specific operations. Then we use a simulation of a topological quantum computer to explicitly demonstrate a quantum computation using Fibonacci anyons, evaluating the Jones polynomial of a selection of simple knots. In addition to simulating a modular circuit-style quantum algorithm, we also show how the magnitude of the Jones polynomial at specific points could be obtained exactly using Fibonacci or Ising anyons. Such an exact algorithm seems ideally suited for a proof of concept demonstration of a topological quantum computer.


title: Symmetry Enforced Self-Learning Monte Carlo Method Applied to the Holstein Model

abstract: Self-learning Monte Carlo method (SLMC), using a trained effective model to guide Monte Carlo sampling processes, is a powerful general-purpose numerical method recently introduced to speed up simulations in (quantum) many-body systems. In this work, we further improve the efficiency of SLMC by enforcing physical symmetries on the effective model. We demonstrate its effectiveness in the Holstein Hamiltonian, one of the most fundamental many-body descriptions of electron-phonon coupling. Simulations of the Holstein model are notoriously difficult due to the combination of the typical cubic scaling of fermionic Monte Carlo and the presence of extremely long autocorrelation times. Our method addresses both bottlenecks. This enables simulations on large lattices in the most difficult parameter regions, and evaluation of the critical point for the charge density wave transition at half-filling with high precision. We argue that our work opens a new research area of quantum Monte Carlo (QMC), providing a general procedure to deal with ergodicity in situations involving Hamiltonians with multiple, distinct low energy states.


title: A Treatise on Information Geometry

abstract: This is the first version of a submission made to vixra back in August 2009. It is concerned with the exploration of various ideas due to Roy Frieden of the University of Arizona in his manuscript "Physics from Fisher Information", which was originally published in December 1998 by Cambridge University Press. In addition to exploration and investigation of same, a key aim of this article is/was to further develop the ideas mentioned in the aforementioned work, as well as explore various exotic differentiable structures and their relationship to physics. Apart from the original component of this manuscript, a series of survey chapters are provided, in the interest of keeping the treatise self-contained. The first summarises the main preliminary results on the existence of non-standard structures on manifolds from the Milnor-Steenrod school. The second is a standard introduction to semi-riemannian geometry. The third introduces the language of geometric measure theory, which is important in justifying the existence of smooth solutions to variational problems with smooth structures and smooth integrands. The fourth is a short remark on PDE existence theory, which is needed for the fifth, which is essentially a typeset version of a series of lectures given by Ben Andrews and Gerhard Huisken on the Hamilton-Perelman program for proving the Geometrisation Conjecture of Bill Thurston.


title: Learning to Race through Coordinate Descent Bayesian Optimisation

abstract: In the automation of many kinds of processes, the observable outcome can often be described as the combined effect of an entire sequence of actions, or controls, applied throughout its execution. In these cases, strategies to optimise control policies for individual stages of the process might not be applicable, and instead the whole policy might have to be optimised at once. On the other hand, the cost to evaluate the policy's performance might also be high, being desirable that a solution can be found with as few interactions as possible with the real system. We consider the problem of optimising control policies to allow a robot to complete a given race track within a minimum amount of time. We assume that the robot has no prior information about the track or its own dynamical model, just an initial valid driving example. Localisation is only applied to monitor the robot and to provide an indication of its position along the track's centre axis. We propose a method for finding a policy that minimises the time per lap while keeping the vehicle on the track using a Bayesian optimisation (BO) approach over a reproducing kernel Hilbert space. We apply an algorithm to search more efficiently over high-dimensional policy-parameter spaces with BO, by iterating over each dimension individually, in a sequential coordinate descent-like scheme. Experiments demonstrate the performance of the algorithm against other methods in a simulated car racing environment.


title: Virtual Immersive Reality for Stated Preference Travel Behaviour Experiments: A Case Study of Autonomous Vehicles on Urban Roads

abstract: Stated preference experiments have been known to suffer from the lack of realism. This issue is particularly visible when the scenario doesn't have a well understood prior reference e.g. in case of the autonomous vehicles related scenarios. We present Virtual Immersive Reality Environment (VIRE) that is capable of developing highly realistic, immersive, and interactive choice scenario. We demonstrate the use of VIRE in the pedestrian preferences related to autonomous vehicles and associated infrastructure changes on urban streets of Montr\'eal. The results are compared with predominantly used approaches i.e. text-only and visual aid. We show that VIRE results in better understanding of the scenario and consistent results.


title: Semi-supervised multi-task learning for lung cancer diagnosis

abstract: Early detection of lung nodules is of great importance in lung cancer screening. Existing research recognizes the critical role played by CAD systems in early detection and diagnosis of lung nodules. However, many CAD systems, which are used as cancer detection tools, produce a lot of false positives (FP) and require a further FP reduction step. Furthermore, guidelines for early diagnosis and treatment of lung cancer are consist of different shape and volume measurements of abnormalities. Segmentation is at the heart of our understanding of nodules morphology making it a major area of interest within the field of computer aided diagnosis systems. This study set out to test the hypothesis that joint learning of false positive (FP) nodule reduction and nodule segmentation can improve the computer aided diagnosis (CAD) systems' performance on both tasks. To support this hypothesis we propose a 3D deep multi-task CNN to tackle these two problems jointly. We tested our system on LUNA16 dataset and achieved an average dice similarity coefficient (DSC) of 91% as segmentation accuracy and a score of nearly 92% for FP reduction. As a proof of our hypothesis, we showed improvements of segmentation and FP reduction tasks over two baselines. Our results support that joint training of these two tasks through a multi-task learning approach improves system performance on both. We also showed that a semi-supervised approach can be used to overcome the limitation of lack of labeled data for the 3D segmentation task.


title: CREPE: A Convolutional Representation for Pitch Estimation

abstract: The task of estimating the fundamental frequency of a monophonic sound recording, also known as pitch tracking, is fundamental to audio processing with multiple applications in speech processing and music information retrieval. To date, the best performing techniques, such as the pYIN algorithm, are based on a combination of DSP pipelines and heuristics. While such techniques perform very well on average, there remain many cases in which they fail to correctly estimate the pitch. In this paper, we propose a data-driven pitch tracking algorithm, CREPE, which is based on a deep convolutional neural network that operates directly on the time-domain waveform. We show that the proposed model produces state-of-the-art results, performing equally or better than pYIN. Furthermore, we evaluate the model's generalizability in terms of noise robustness. A pre-trained version of CREPE is made freely available as an open-source Python module for easy application.


title: Towards Realisation of Heterogeneous Earth-Observation Sensor Database Framework for the Sensor Observation Service based on PostGIS

abstract: Environmental monitoring and management systems in most cases deal with models and spatial analytics that involve the integration of in-situ and remote Geosensor observations. In-situ sensor observations and those gathered by remote sensors are usually provided by different databases and services in real-time dynamic services such as the Geo-Web Services. Thus, data have to be pulled from different databases and transferred over the network before they are fused and processed on the service middleware. This process is very massive and unnecessary communication-work load on the service middleware. Massive work load in large raster downloads from flat-file raster data sources each time a request is made and huge integration and geo-processing work load on the service middleware which could actually be better leveraged at the database This paper therefore proposes the realization of heterogeneous sensor database framework based on PostGIS for integration, geo-processing and spatial analysis of remote and in-situ sensor observations at the database level. Also discussed in this paper is how the framework can be integrated in the Sensor Observation Service (SOS) to reduce communication and massive workload on the Geospatial Web Services and as well make query request from the user end a lot more flexible. Keywords: Earth-Observation, Heterogeneous Earth-Observation Sensor Database, PostGIS , Sensor Observation Service.


title: Hubble Space Telescope Observations of Extended [O III] {\lambda}5007 Emission in Nearby QSO2s: New Constraints On AGN / Host Galaxy Interaction

abstract: We present a Hubble Space Telescope (HST) survey of extended [O III] {\lambda}5007 emission for a sample of 12 nearby (z < 0.12), luminous Type 2 quasars (QSO2s), which we use to measure the extent and kinematics of their AGN-ionized gas. We find the size of the observed [O III] regions scale with luminosity in comparison to nearby, less luminous Seyfert galaxies and radially outflowing kinematics to exist in all targets. We report an average maximum outflow radius of $\sim$600 pc, with gas continuing to be kinematically influenced by the central AGN out to an average radius of $\sim$1130 pc. These findings question the effectiveness of AGN being capable of clearing material from their host bulge in the nearby universe and suggest that disruption of gas by AGN activity may prevent star formation without requiring evacuation. Additionally, we find a dichotomy in our targets when comparing [O III] radial extent and nuclear FWHM, where QSO2s with compact [O III] morphologies typically possess broader nuclear emission-lines.


title: Building a Word Segmenter for Sanskrit Overnight

abstract: There is an abundance of digitised texts available in Sanskrit. However, the word segmentation task in such texts are challenging due to the issue of 'Sandhi'. In Sandhi, words in a sentence often fuse together to form a single chunk of text, where the word delimiter vanishes and sounds at the word boundaries undergo transformations, which is also reflected in the written text. Here, we propose an approach that uses a deep sequence to sequence (seq2seq) model that takes only the sandhied string as the input and predicts the unsandhied string. The state of the art models are linguistically involved and have external dependencies for the lexical and morphological analysis of the input. Our model can be trained "overnight" and be used for production. In spite of the knowledge lean approach, our system preforms better than the current state of the art by gaining a percentage increase of 16.79 % than the current state of the art.


title: Optimal Single Sample Tests for Structured versus Unstructured Network Data

abstract: We study the problem of testing, using only a single sample, between mean field distributions (like Curie-Weiss, Erd\H{o}s-R\'enyi) and structured Gibbs distributions (like Ising model on sparse graphs and Exponential Random Graphs). Our goal is to test without knowing the parameter values of the underlying models: only the \emph{structure} of dependencies is known. We develop a new approach that applies to both the Ising and Exponential Random Graph settings based on a general and natural statistical test. The test can distinguish the hypotheses with high probability above a certain threshold in the (inverse) temperature parameter, and is optimal in that below the threshold no test can distinguish the hypotheses. The thresholds do not correspond to the presence of long-range order in the models. By aggregating information at a global scale, our test works even at very high temperatures. The proofs are based on distributional approximation and sharp concentration of quadratic forms, when restricted to Hamming spheres. The restriction to Hamming spheres is necessary, since otherwise any scalar statistic is useless without explicit knowledge of the temperature parameter. At the same time, this restriction radically changes the behavior of the functions under consideration, resulting in a much smaller variance than in the independent setting; this makes it hard to directly apply standard methods (i.e., Stein's method) for concentration of weakly dependent variables. Instead, we carry out an additional tensorization argument using a Markov chain that respects the symmetry of the Hamming sphere.


title: Quenching of stimulated Raman scattering in the kinetic regime by external magnetic fields

abstract: We show via particle-in-cell simulations that small normalized magnetic fields ($\omega_c/\omega_p \ll 1$) applied perpendicularly to a light wave can significantly modify the evolution of backward stimulated Raman scattering (SRS) in the kinetic regime. The presence of the magnetic field increases the threshold for kinetic inflation and decreases the amount of reflectivity when SRS is driven significantly above threshold. Analysis indicates this arises because trapped electrons are accelerated as they surf across the wave, leading to the continual dissipation of the electron plasma waves over a wider range of wave amplitudes. The simulation parameters are directly relevant for SRS in inertial confinement fusion devices and indicate that approximately 30 Tesla magnetic fields might significantly reduce SRS backscatter.


title: Some determinants involving incomplete Fubini numbers

abstract: We study some properties of restricted and associated Fubini numbers. In particular, they have the natural extensions of the original Fubini numbers in the sense of determinants. We also introduce modified Bernoulli and Cauchy numbers and study characteristic properties.


title: Contrast Subgraph Mining from Coherent Cores

abstract: Graph pattern mining methods can extract informative and useful patterns from large-scale graphs and capture underlying principles through the overwhelmed information. Contrast analysis serves as a keystone in various fields and has demonstrated its effectiveness in mining valuable information. However, it has been long overlooked in graph pattern mining. Therefore, in this paper, we introduce the concept of contrast subgraph, that is, a subset of nodes that have significantly different edges or edge weights in two given graphs of the same node set. The major challenge comes from the gap between the contrast and the informativeness. Because of the widely existing noise edges in real-world graphs, the contrast may lead to subgraphs of pure noise. To avoid such meaningless subgraphs, we leverage the similarity as the cornerstone of the contrast. Specifically, we first identify a coherent core, which is a small subset of nodes with similar edge structures in the two graphs, and then induce contrast subgraphs from the coherent cores. Moreover, we design a general family of coherence and contrast metrics and derive a polynomial-time algorithm to efficiently extract contrast subgraphs. Extensive experiments verify the necessity of introducing coherent cores as well as the effectiveness and efficiency of our algorithm. Real-world applications demonstrate the tremendous potentials of contrast subgraph mining.


title: Tests about R multivariate simple linear models

abstract: Hypothesis about the parallelism of the regression lines in R multivariate simple linear models are studied in this paper. Tests on common intercept and sets of lines intersected at a fixed value, are also developed. An application in an agricultural entomology context is provided.


title: Organic Photodiodes with an Extended Responsivity using Ultrastrong Light-Matter Coupling

abstract: In organic photodiodes (OPDs) light is absorbed by excitons, which dissociate to generate photocurrent. Here, we demonstrate a novel type of OPD in which light is absorbed by polaritons, hybrid light-matter states. We demonstrate polariton OPDs operating in the ultra-strong coupling regime at visible and infrared wavelengths. These devices can be engineered to show narrow responsivity with a very weak angle-dependence. More importantly, they can be tuned to operate in a spectral range outside that of the bare exciton absorption. Remarkably, we show that the responsivity of a polariton OPD can be pushed to near infrared wavelengths, where few organic absorbers are available, with external quantum efficiencies exceeding those of a control OPD.


title: A Re-solving Heuristic with Uniformly Bounded Loss for Network Revenue Management

abstract: We consider the canonical (quantity-based) network revenue management problem, where a firm accepts or rejects incoming customer requests irrevocably in order to maximize expected revenue given limited resources. Due to the curse of dimensionality, the exact solution to this problem by dynamic programming is intractable when the number of resources is large. We study a family of re-solving heuristics that periodically re-optimize an approximation to the original problem known as the deterministic linear program (DLP), where random customer arrivals are replaced by their expectations. We find that, in general, frequently re-solving the DLP produces the same order of revenue loss as one would get without re-solving, which scales as the square root of the time horizon length and resource capacities. By re-solving the DLP at a few selected points in time and applying thresholds to the customer acceptance probabilities, we design a new re-solving heuristic whose revenue loss is uniformly bounded by a constant that is independent of the time horizon and resource capacities.


title: The least prime ideal in a given ideal class

abstract: Let $K$ be a number field with the discriminant $D_K$ and the class number $h_{K}$, which has bounded degree over $\mathbb{Q}$. By assuming GRH, we prove that every ideal class of $K$ contains a prime ideal with norm less than $h_{K}^2\log(D_K)^{2}$ and also all but $o(h_K)$ of them have a prime ideal with norm less than $h_{K}\log(D_K)^{2+\epsilon}$. For imaginary quadratic fields $K=\mathbb{Q}(\sqrt{D})$, by assuming Conjecture~\ref{piarcor} (a weak version of the pair correlation conjecure), we improve our bounds by removing a factor of $\log(D)$ from our bounds and show that these bounds are optimal.


title: HWNet v2: An Efficient Word Image Representation for Handwritten Documents

abstract: We present a framework for learning an efficient holistic representation for handwritten word images. The proposed method uses a deep convolutional neural network with traditional classification loss. The major strengths of our work lie in: (i) the efficient usage of synthetic data to pre-train a deep network, (ii) an adapted version of the ResNet-34 architecture with the region of interest pooling (referred to as HWNet v2) which learns discriminative features for variable sized word images, and (iii) a realistic augmentation of training data with multiple scales and distortions which mimics the natural process of handwriting. We further investigate the process of transfer learning to reduce the domain gap between synthetic and real domain, and also analyze the invariances learned at different layers of the network using visualization techniques proposed in the literature. Our representation leads to a state-of-the-art word spotting performance on standard handwritten datasets and historical manuscripts in different languages with minimal representation size. On the challenging IAM dataset, our method is first to report an mAP of around 0.90 for word spotting with a representation size of just 32 dimensions. Furthermore, we also present results on printed document datasets in English and Indic scripts which validates the generic nature of the proposed framework for learning word image representation.


title: High Speed SRT Divider for Intelligent Embedded System

abstract: Increasing development in embedded systems, VLSI and processor design have given rise to increased demands from the system in terms of power, speed, area, throughput etc. Most of the sophisticated embedded system applications consist of processors, which now need an arithmetic unit with the ability to execute complex division operations with maximum efficiency. Hence the speed of the arithmetic unit is critically dependent on division operation. Most of the dividers use the SRT division algorithm for division. In IoT and other embedded applications, typically radix 2 and radix 4 division algorithms are used. The proposed algorithm lies on parallel execution of various steps so as to reduce time critical path, use fuzzy logic to solve the overlap problem in quotient selection, hence reducing maximum delay and increasing the accuracy. Every logical circuit has a maximum delay on which the timing of the circuit is dependent and the path, causing the maximum delay is known as the critical path. Our approach uses the previous SRT algorithm methods to make a highly parallel pipelined design and use Mamdani model to determine a solution to the overlapping problem to reduce the overall execution time of radix 4 SRT division on 64 bits double precision floating point numbers to 281ns. The design is made using Bluespec System Verilog, synthesized and simulated using Vivado v.2016.1 and implemented on Xilinx Virtex UltraScale FPGA board.


title: Can Network Embedding of Distributional Thesaurus be Combined with Word Vectors for Better Representation?

abstract: Distributed representations of words learned from text have proved to be successful in various natural language processing tasks in recent times. While some methods represent words as vectors computed from text using predictive model (Word2vec) or dense count based model (GloVe), others attempt to represent these in a distributional thesaurus network structure where the neighborhood of a word is a set of words having adequate context overlap. Being motivated by recent surge of research in network embedding techniques (DeepWalk, LINE, node2vec etc.), we turn a distributional thesaurus network into dense word vectors and investigate the usefulness of distributional thesaurus embedding in improving overall word representation. This is the first attempt where we show that combining the proposed word representation obtained by distributional thesaurus embedding with the state-of-the-art word representations helps in improving the performance by a significant margin when evaluated against NLP tasks like word similarity and relatedness, synonym detection, analogy detection. Additionally, we show that even without using any handcrafted lexical resources we can come up with representations having comparable performance in the word similarity and relatedness tasks compared to the representations where a lexical resource has been used.


title: Approximations to galaxy star formation rate histories: properties and uses of two examples

abstract: Galaxies evolve via a complex interaction of numerous different physical processes, scales and components. In spite of this, overall trends often appear. Simplified models for galaxy histories can be used to search for and capture such emergent trends, and thus to interpret and compare results of galaxy formation models to each other and to nature. Here, two approximations are applied to galaxy integrated star formation rate histories, drawn from a semi-analytic model grafted onto a dark matter simulation. Both a lognormal functional form and principal component analysis (PCA) approximate the integrated star formation rate histories fairly well. Machine learning, based upon simplified galaxy halo histories, is somewhat successful at recovering both fits. The fits to the histories give fixed time star formation rates which have notable scatter from their true fixed time rates at final time, especially for quiescent and "green valley" galaxies, and more so for the PCA fit. For classifying galaxies into subfamilies sharing similar integrated histories, both approximations are better than using final stellar mass or specific star formation rate. Several subsamples from the simulation illustrate how these simple parameterizations can provide points of contact for comparisons between different galaxy formation samples, or more generally, models. As a side result, the halo masses of simulated galaxies with early peak star formation rate (according to the lognormal fit) lie on one of two tracks. The small fraction of galaxies with a lower halo mass at peak star formation rate appear to stall in their halo growth, even though they are central in their host halos.


title: Mass gap in the weak coupling limit of $(2+1)$ SU(2) lattice gauge theory

abstract: We develop the dual description of $2+1$ SU(2) lattice gauge theory as interacting `abelian like' electric loops by using Schwinger bosons. "Point splitting" of the lattice enables us to construct explicit Hilbert space for the gauge invariant theory which in turn makes dynamics more transparent. Using path integral representation in phase space, the interacting closed loop dynamics is analyzed in the weak coupling limit to get the mass gap.


title: Drift-Free Indoor Navigation Using Simultaneous Localization and Mapping of the Ambient Heterogeneous Magnetic Field

abstract: In the absence of external reference position information (e.g. GNSS) SLAM has proven to be an effective method for indoor navigation. The positioning drift can be reduced with regular loop-closures and global relaxation as the backend, thus achieving a good balance between exploration and exploitation. Although vision-based systems like laser scanners are typically deployed for SLAM, these sensors are heavy, energy inefficient, and expensive, making them unattractive for wearables or smartphone applications. However, the concept of SLAM can be extended to non-optical systems such as magnetometers. Instead of matching features such as walls and furniture using some variation of the ICP algorithm, the local magnetic field can be matched to provide loop-closure and global trajectory updates in a Gaussian Process (GP) SLAM framework. With a MEMS-based inertial measurement unit providing a continuous trajectory, and the matching of locally distinct magnetic field maps, experimental results in this paper show that a drift-free navigation solution in an indoor environment with millimetre-level accuracy can be achieved. The GP-SLAM approach presented can be formulated as a maximum a posteriori estimation problem and it can naturally perform loop-detection, feature-to-feature distance minimization, global trajectory optimization, and magnetic field map estimation simultaneously. Spatially continuous features (i.e. smooth magnetic field signatures) are used instead of discrete feature correspondences (e.g. point-to-point) as in conventional vision-based SLAM. These position updates from the ambient magnetic field also provide enough information for calibrating the accelerometer and gyroscope bias in-use. The only restriction for this method is the need for magnetic disturbances (which is typically not an issue indoors); however, no assumptions are required for the general motion of the sensor.


title: The Ando-Hiai inequalities for the solution of the generalized Karcher equation and related results

abstract: In this paper, we shall give a concrete relation between the generalized Karcher equation and operator means as its solution. Next, we shall show two types of the Ando-Hiai inequalities for the solution of the generalized Karcher equation. In this discussion, we also give a property of relative operator entropy.


title: Multiple Object Trajectography Using Particle Swarm Optimization Combined to Hungarian Method

abstract: The problem of simultaneous trajectography of several dynamical objects is formulated as an optimization problem. The available observations consist in a series of photographs showing undiscriminated objects. The goal is to find the object initial states so that the resulting trajectories match as well as possible the set of observations. An assignment problem is solved at each observation date by the Hungarian method, yielding a deviation cost between the simulated trajectories and the measurements. A fitness function summing the deviation costs is minimized by a particle swarm algorithm. The method is illustrated on a space orbitography application.


title: Coplanar Low-Thrust Transfer with Eclipses Using Analytical Costate Guess

abstract: Low-thrust orbital transfers are difficult to optimize by indirect methods. The main issues come from the costate guess and from the numerical propagation accuracy required by the shooting method. In the case of a coplanar minimum-time low-thrust transfer with eclipses, an analytical costate guess is proposed. The optimal control problem reduces to an unconstrained minimization problem with two unknowns. A derivative-free algorithm yields a quasi-optimal solution from scratch in a few minutes. No specific guess is necessary and the algorithm properties ensure finding the global minimum of the unconstrained problem. The method is applicable to thrust levels from large to very low and to any eclipse configuration, as exemplified on a transfer towards the geostationary orbit.


title: Numerical solution of boundary value problems for the eikonal equation in an anisotropic medium

abstract: A Dirichlet problem is considered for the eikonal equation in an anisotropic medium. The nonlinear boundary value problem (BVP) formulated in the present work is the limit of the diffusion-reaction problem with a reaction parameter tending to infinity. To solve numerically the singularly perturbed diffusion-reaction problem, monotone approximations are employed. Numerical examples are presented for a two-dimensional BVP for the eikonal equation in an anisotropic medium. The standard piecewise-linear finite-element approximation in space is used in computations.


title: Approximate Set Union Via Approximate Randomization

abstract: We develop an randomized approximation algorithm for the size of set union problem $\arrowvert A_1\cup A_2\cup...\cup A_m\arrowvert$, which given a list of sets $A_1,...,A_m$ with approximate set size $m_i$ for $A_i$ with $m_i\in \left((1-\beta_L)|A_i|, (1+\beta_R)|A_i|\right)$, and biased random generators with $Prob(x=\randomElm(A_i))\in \left[{1-\alpha_L\over |A_i|},{1+\alpha_R\over |A_i|}\right]$ for each input set $A_i$ and element $x\in A_i,$ where $i=1, 2, ..., m$. The approximation ratio for $\arrowvert A_1\cup A_2\cup...\cup A_m\arrowvert$ is in the range $[(1-\epsilon)(1-\alpha_L)(1-\beta_L), (1+\epsilon)(1+\alpha_R)(1+\beta_R)]$ for any $\epsilon\in (0,1)$, where $\alpha_L, \alpha_R, \beta_L,\beta_R\in (0,1)$. The complexity of the algorithm is measured by both time complexity, and round complexity. The algorithm is allowed to make multiple membership queries and get random elements from the input sets in one round. Our algorithm makes adaptive accesses to input sets with multiple rounds. Our algorithm gives an approximation scheme with $O(\setCount\cdot(\log \setCount)^{O(1)})$ running time and $O(\log m)$ rounds, where $m$ is the number of sets. Our algorithm can handle input sets that can generate random elements with bias, and its approximation ratio depends on the bias. Our algorithm gives a flexible tradeoff with time complexity $O\left(\setCount^{1+\xi}\right)$ and round complexity $O\left({1\over \xi}\right)$ for any $\xi\in(0,1)$.


title: Towards Principled Design of Deep Convolutional Networks: Introducing SimpNet

abstract: Major winning Convolutional Neural Networks (CNNs), such as VGGNet, ResNet, DenseNet, \etc, include tens to hundreds of millions of parameters, which impose considerable computation and memory overheads. This limits their practical usage in training and optimizing for real-world applications. On the contrary, light-weight architectures, such as SqueezeNet, are being proposed to address this issue. However, they mainly suffer from low accuracy, as they have compromised between the processing power and efficiency. These inefficiencies mostly stem from following an ad-hoc designing procedure. In this work, we discuss and propose several crucial design principles for an efficient architecture design and elaborate intuitions concerning different aspects of the design procedure. Furthermore, we introduce a new layer called {\it SAF-pooling} to improve the generalization power of the network while keeping it simple by choosing best features. Based on such principles, we propose a simple architecture called {\it SimpNet}. We empirically show that SimpNet provides a good trade-off between the computation/memory efficiency and the accuracy solely based on these primitive but crucial principles. SimpNet outperforms the deeper and more complex architectures such as VGGNet, ResNet, WideResidualNet \etc, on several well-known benchmarks, while having 2 to 25 times fewer number of parameters and operations. We obtain state-of-the-art results (in terms of a balance between the accuracy and the number of involved parameters) on standard datasets, such as CIFAR10, CIFAR100, MNIST and SVHN. The implementations are available at \href{url}{https://github.com/Coderx7/SimpNet}.


title: Observationally quantified reconnection providing a viable mechanism for active region coronal heating

abstract: The heating of the Sun's corona has been explained by several different mechanisms including wave dissipation and magnetic reconnection. While both have been shown capable of supplying the requisite power, neither has been used in a quantitative model of observations fed by measured inputs. Here we show that impulsive reconnection is capable of producing an active region corona agreeing both qualitatively and quantitatively with extreme-ultraviolet observations. We calculate the heating power proportional to the velocity difference between magnetic footpoints and the photospheric plasma, called the non-ideal velocity. The length scale of flux elements reconnected in the corona is found to be around 160 km. The differential emission measure of the model corona agrees with that derived using multi-wavelength images. Synthesized extreme-ultraviolet images resemble observations both in their loop-dominated appearance and their intensity histograms. This work provides compelling evidence that impulsive reconnection events are a viable mechanism for heating the corona.


title: Randomness of formal languages via automatic martingales

abstract: We define a notion of randomness for individual and collections of formal languages based on automatic martingales acting on sequences of words from some underlying domain. An automatic martingale bets if the incoming word belongs to the target language or not. Then randomness of both single languages and collections of languages is defined as a failure of automatic martingale to gain an unbounded capital by betting on the target language according to an incoming sequence of words, or a text. The randomness of formal languages turned out to be heavily dependent on the text. For very general classes of texts, any nonregular language happens to be random when considered individually. As for collections of languages, very general classes of texts permits nonrandomness of automatic families of languages only. On the other hand, an arbitrary computable language is be shown to be nonrandom under certain dynamic texts.


title: On the Extensions of a Discrete Valuation in a Number Field

abstract: Let $K$ be a number field defined by a monic irreducible polynomial $F(X) \in \mathbb{Z}[X]$, $p$ a fixed rational prime, and $\nu_p$ the discrete valuation associated to $p$. Assume that $\overline{F}(X)$ factors modulo $p$ into the product of powers of $r$ distinct monic irreducible polynomials. We present in this paper a condition, weaker than the known ones, which guarantees the existence of exactly $r$ valuations of $K$ extending $\nu_p$. We further specify the ramification indices and residue degrees of these extended valuations in such a way that generalizes the known estimates. Some useful remarks and computational examples are also given to highlight some improvements due to our result.


title: Sentiment Analysis on Speaker Specific Speech Data

abstract: Sentiment analysis has evolved over past few decades, most of the work in it revolved around textual sentiment analysis with text mining techniques. But audio sentiment analysis is still in a nascent stage in the research community. In this proposed research, we perform sentiment analysis on speaker discriminated speech transcripts to detect the emotions of the individual speakers involved in the conversation. We analyzed different techniques to perform speaker discrimination and sentiment analysis to find efficient algorithms to perform this task.


title: Very true pseudo-BCK algebras

abstract: In this paper we introduce the very true operators on pseudo-BCK algebras and we study their properties. We prove that the composition of two very true operators is a very true operator if and only if they commute. Moreover, given a very true bounded pseudo-BCK algebra $(A,v)$, we define the pseudo-BCK$_{vt,st}$ algebra by adding two truth-depressing hedges operators associated with $v$. We also define the very true deductive systems and the very true homomorphisms and we investigate their properties. Also, given a normal $v$-deductive system of a very true pseudo-BCK algebra $(A,v)$ we construct a very true operator on the quotient pseudo-BCK algebra $A/H$. We investigate the very true operators on some classes of pseudo-BCK algebras such as pseudo-BCK(pP) algebras, FL$_w$-algebras, pseudo-MTL algebras. Finally, we define the $Q$-Smarandache pseudo-BCK algebras and we introduce the notion of a very true operator on $Q$-Smarandache pseudo-BCK algebras.


title: Gluons, Heavy and Light Quarks in the QCD Vacuum

abstract: We are discussing the properties of the QCD vacuum which might be important especially for the understanding of hadrons with small quark core size $\sim0.3\,fm.$ We assume that at these distances the QCD vacuum can be described by the Instanton Liquid Model (ILM). At larger distances, where confinement is important, ILM should be extended to Dyons Liquid Model (DLM). The ILM has only two free parameters, average instanton size $\rho\approx0.3\,fm$ and average inter-instanton distance $R\approx1\,fm$, and can successfully describe the key features of light hadron physics. One of the important conceptual results was prediction of the momentum dependent dynamical quark mass $M\sim(packing\,\,fraction)^{1/2}\,\,\rho^{-1}\approx360\,MeV$, later confirmed numerically by evaluations in the lattice. The estimates show that gluon-instanton interaction strength is also big and is controlled by the value of dynamical gluon mass $M_{g}\approx M$. Heavy quarks interact with instantons much weaker. The heavy quark-instanton interaction strength is given by $\Delta m_{Q}\sim packing\,\,fraction\,\,\rho^{-1}\approx70\,MeV.$ Nevertheless, the direct instanton contribution to the colorless heavy-heavy quarks potential is sizable and must be taken into account. At small distances, where one-gluon exchange contribution to this potential is dominated, we have to take into account dynamical gluon mass $M_{g}$. Also, instantons are generating light-heavy quarks interactions and allow to describe the nonperturbative effects in heavy-light quarks systems.


title: Multi-Pass Streaming Algorithms for Monotone Submodular Function Maximization

abstract: We consider maximizing a monotone submodular function under a cardinality constraint or a knapsack constraint in the streaming setting. In particular, the elements arrive sequentially and at any point of time, the algorithm has access to only a small fraction of the data stored in primary memory. We propose the following streaming algorithms taking $O(\varepsilon^{-1})$ passes: ----a $(1-e^{-1}-\varepsilon)$-approximation algorithm for the cardinality-constrained problem ---- a $(0.5-\varepsilon)$-approximation algorithm for the knapsack-constrained problem. Both of our algorithms run in $O^\ast(n)$ time, using $O^\ast(K)$ space, where $n$ is the size of the ground set and $K$ is the size of the knapsack. Here the term $O^\ast$ hides a polynomial of $\log K$ and $\varepsilon^{-1}$. Our streaming algorithms can also be used as fast approximation algorithms. In particular, for the cardinality-constrained problem, our algorithm takes $O(n\varepsilon^{-1} \log (\varepsilon^{-1}\log K) )$ time, improving on the algorithm of Badanidiyuru and Vondr\'{a}k that takes $O(n \varepsilon^{-1} \log (\varepsilon^{-1} K) )$ time.


title: Generating alternating and symmetric groups with two elements of fixed order

abstract: We answer a question raised by Lanier about the possibility of generating $A_n$ and $S_n$ with two elements of order $k$, where $n \geqslant k \geqslant 3$. We show that this can always be done apart from some clear exceptions.


title: A New De-blurring Technique for License Plate Images with Robust Length Estimation

abstract: Recognizing a license plate clearly while seeing a surveillance camera snapshot is often important in cases where the troublemaker vehicle(s) have to be identified. In many real world situations, these images are blurred due to fast motion of the vehicle and cannot be recognized by the human eye. For this kind of blurring, the kernel involved can be said to be a linear uniform convolution described by its angle and length. We propose a new de-blurring technique in this paper to parametrically estimate the kernel as accurately as possible with emphasis on the length estimation process. We use a technique which employs Hough transform in estimating the kernel angle. To accurately estimate the kernel length, a novel approach using the cepstral transform is introduced. We compare the de-blurred results obtained using our scheme with those of other recently introduced blind de-blurring techniques. The comparisons corroborate that our scheme can remove a large blur from the image captured by the camera to recover vital semantic information about the license plate.


title: HyP-DESPOT: A Hybrid Parallel Algorithm for Online Planning under Uncertainty

abstract: Planning under uncertainty is critical for robust robot performance in uncertain, dynamic environments, but it incurs high computational cost. State-of-the-art online search algorithms, such as DESPOT, have vastly improved the computational efficiency of planning under uncertainty and made it a valuable tool for robotics in practice. This work takes one step further by leveraging both CPU and GPU parallelization in order to achieve near real-time online planning performance for complex tasks with large state, action, and observation spaces. Specifically, we propose Hybrid Parallel DESPOT (HyP-DESPOT), a massively parallel online planning algorithm that integrates CPU and GPU parallelism in a multi-level scheme. It performs parallel DESPOT tree search by simultaneously traversing multiple independent paths using multi-core CPUs and performs parallel Monte-Carlo simulations at the leaf nodes of the search tree using GPUs. Experimental results show that HyP-DESPOT speeds up online planning by up to several hundred times, compared with the original DESPOT algorithm, in several challenging robotic tasks in simulation.


title: Disentangling covariant Wigner functions for chiral fermions

abstract: We develop a general formalism for the quantum kinetics of chiral fermions in a background electromagnetic field based on a semiclassical expansion of covariant Wigner functions in the Planck constant $\hbar$. We demonstrate to any order of $\hbar$ that only the time-component of the Wigner function is independent while other components are explicit derivative. We further demonstrate to any order of $\hbar$ that a system of quantum kinetic equations for multiple-components of Wigner functions can be reduced to one chiral kinetic equation involving only the single-component distribution function. These are remarkable properties of the quantum kinetics of chiral fermions and will significantly simplify the description and simulation of chiral effects in heavy ion collisions and Dirac/Weyl semimetals. We present the unintegrated chiral kinetic equations in four-momenta up to $O(\hbar ^2)$ and the integrated ones in three-momenta up to $O(\hbar)$. We find that some singular terms emerge in the integration over the time component of the four-momentum, which result in a new source term contributing to the chiral anomaly, in contrast to the well-known scenario of the Berry phase term. Finally we rewrite our results in any Lorentz frame with a reference four-velocity and show how the non-trivial transformation of the distribution function in different frames emerges in a natural way.


title: Design and Implementation of the Andromeda Proof Assistant

abstract: Andromeda is an LCF-style proof assistant where the user builds derivable judgments by writing code in a meta-level programming language AML. The only trusted component of Andromeda is a minimalist nucleus (an implementation of the inference rules of an object-level type theory), which controls construction and decomposition of type-theoretic judgments. Since the nucleus does not perform complex tasks like equality checking beyond syntactic equality, this responsibility is delegated to the user, who implements one or more equality checking procedures in the meta-language. The AML interpreter requests witnesses of equality from user code using the mechanism of algebraic operations and handlers. Dynamic checks in the nucleus guarantee that no invalid object-level derivations can be constructed. %even if the AML code (or interpreter) is untrusted. To demonstrate the flexibility of this system structure, we implemented a nucleus consisting of dependent type theory with equality reflection. Equality reflection provides a very high level of expressiveness, as it allows the user to add new judgmental equalities, but it also destroys desirable meta-theoretic properties of type theory (such as decidability and strong normalization). The power of effects and handlers in AML is demonstrated by a standard library that provides default algorithms for equality checking, computation of normal forms, and implicit argument filling. Users can extend these new algorithms by providing local "hints" or by completely replacing these algorithms for particular developments. We demonstrate the resulting system by showing how to axiomatize and compute with natural numbers, by axiomatizing the untyped $\lambda$-calculus, and by implementing a simple automated system for managing a universe of types.


title: Theoretical study of the flow in a fluid damper containing high viscosity silicone oil: effects of shear-thinning and viscoelasticity

abstract: The flow inside a fluid damper where a piston reciprocates sinusoidally inside an outer casing containing high-viscosity silicone oil is simulated using a Finite Volume method, at various excitation frequencies. The oil is modelled by the Carreau-Yasuda (CY) and Phan-Thien \& Tanner (PTT) constitutive equations. Both models account for shear-thinning but only the PTT model accounts for elasticity. The CY and other generalised Newtonian models have been previously used in theoretical studies of fluid dampers, but the present study is the first to perform full two-dimensional (axisymmetric) simulations employing a viscoelastic constitutive equation. It is found that the CY and PTT predictions are similar when the excitation frequency is low, but at medium and higher frequencies the CY model fails to describe important phenomena that are predicted by the PTT model and observed in experimental studies found in the literature, such as the hysteresis of the force-displacement and force-velocity loops. Elastic effects are quantified by applying a decomposition of the damper force into elastic and viscous components, inspired from LAOS (Large Amplitude Oscillatory Shear) theory. The CY model also overestimates the damper force relative to the PTT, because it underpredicts the flow development length inside the piston-cylinder gap. It is thus concluded that (a) fluid elasticity must be accounted for and (b) theoretical approaches that rely on the assumption of one-dimensional flow in the piston-cylinder gap are of limited accuracy, even if they account for fluid viscoelasticity. The consequences of using lower-viscosity silicone oil are also briefly examined.


title: Network Centrality Analysis of Tehran Urban and Suburban Railway System

abstract: Nowadays, Tehran Urban and Suburban Railway System (TUSRS) is going to be completed by eight lines and 149 stations. This complex transportation system contains 168 links between each station pairs and 20 cross-section and Y-branch stations among all eight lines. In this study, we considered TUSRS as a complex network and undertook several analyzes based on graph theory. Examining e.g. centrality measures, we identified central stations within TUSRS. This analysis could be useful for redistributing strategy of the overcrowded stations and improving the organization of maintaining system. These findings are also promising for better designing the systems of tomorrow in other metropolitan areas in Iran.


title: Fusion of finite set distributions: Pointwise consistency and global cardinality

abstract: A recent trend in distributed multi-sensor fusion is to use random finite set filters at the sensor nodes and fuse the filtered distributions algorithmically using their exponential mixture densities (EMDs). Fusion algorithms which extend the celebrated covariance intersection and consensus based approaches are such examples. In this article, we analyse the variational principle underlying EMDs and show that the EMDs of finite set distributions do not necessarily lead to consistent fusion of cardinality distributions. Indeed, we demonstrate that these inconsistencies may occur with overwhelming probability in practice, through examples with Bernoulli, Poisson and independent identically distributed (IID) cluster processes. We prove that pointwise consistency of EMDs does not imply consistency in global cardinality and vice versa. Then, we redefine the variational problems underlying fusion and provide iterative solutions thereby establishing a framework that guarantees cardinality consistent fusion.


title: A new foundational crisis in mathematics, is it really happening?

abstract: The article reconsiders the position of the foundations of mathematics after the discovery of HoTT. Discussion that this discovery has generated in the community of mathematicians, philosophers and computer scientists might indicate a new crisis in the foundation of mathematics. By examining the mathematical facts behind HoTT and their relation with the existing foundations, we conclude that the present crisis is not one. We reiterate a pluralist vision of the foundations of mathematics. The article contains a short survey of the mathematical and historical background needed to understand the main tenets of the foundational issues.


title: Efficient GAN-Based Anomaly Detection

abstract: Generative adversarial networks (GANs) are able to model the complex highdimensional distributions of real-world data, which suggests they could be effective for anomaly detection. However, few works have explored the use of GANs for the anomaly detection task. We leverage recently developed GAN models for anomaly detection, and achieve state-of-the-art performance on image and network intrusion datasets, while being several hundred-fold faster at test time than the only published GAN-based method.


title: The Geodesic Farthest-point Voronoi Diagram in a Simple Polygon

abstract: Given a set of point sites in a simple polygon, the geodesic farthest-point Voronoi diagram partitions the polygon into cells, at most one cell per site, such that every point in a cell has the same farthest site with respect to the geodesic metric. We present an $O(n\log\log n+m\log m)$- time algorithm to compute the geodesic farthest-point Voronoi diagram of $m$ point sites in a simple $n$-gon. This improves the previously best known algorithm by Aronov et al. [Discrete Comput. Geom. 9(3):217-255, 1993]. In the case that all point sites are on the boundary of the simple polygon, we can compute the geodesic farthest-point Voronoi diagram in $O((n + m) \log \log n)$ time.


title: Generating Python Code From Object-Z Specifications

abstract: Object-Z is an object-oriented specification language which extends the Z language with classes, objects, inheritance and polymorphism that can be used to represent the specification of a complex system as collections of objects. There are a number of existing works that mapped Object-Z to C++ and Java programming languages. Since Python and Object-Z share many similarities, both are object-oriented paradigm, support set theory and predicate calculus moreover, Python is a functional programming language which is naturally closer to formal specifications, we propose a mapping from Object-Z specifications to Python code that covers some Object-Z constructs and express its specifications in Python to validate these specifications. The validations are used in the mapping covered preconditions, post-conditions, and invariants that are built using lambda function and Python's decorator. This work has found Python is an excellent language for developing libraries to map Object-Z specifications to Python.


title: A Deep Q-Learning Agent for the L-Game with Variable Batch Training

abstract: We employ the Deep Q-Learning algorithm with Experience Replay to train an agent capable of achieving a high-level of play in the L-Game while self-learning from low-dimensional states. We also employ variable batch size for training in order to mitigate the loss of the rare reward signal and significantly accelerate training. Despite the large action space due to the number of possible moves, the low-dimensional state space and the rarity of rewards, which only come at the end of a game, DQL is successful in training an agent capable of strong play without the use of any search methods or domain knowledge.


title: Post Selection Inference with Incomplete Maximum Mean Discrepancy Estimator

abstract: Measuring divergence between two distributions is essential in machine learning and statistics and has various applications including binary classification, change point detection, and two-sample test. Furthermore, in the era of big data, designing divergence measure that is interpretable and can handle high-dimensional and complex data becomes extremely important. In the paper, we propose a post selection inference (PSI) framework for divergence measure, which can select a set of statistically significant features that discriminate two distributions. Specifically, we employ an additive variant of maximum mean discrepancy (MMD) for features and introduce a general hypothesis test for PSI. A novel MMD estimator using the incomplete U-statistics, which has an asymptotically Normal distribution (under mild assumptions) and gives high detection power in PSI, is also proposed and analyzed theoretically. Through synthetic and real-world feature selection experiments, we show that the proposed framework can successfully detect statistically significant features. Last, we propose a sample selection framework for analyzing different members in the Generative Adversarial Networks (GANs) family.


title: On some geometric properties of operator spaces

abstract: In this paper we study some geometric properties like parallelism, orthogonality and semi-rotundity in the space of bounded linear operators. We completely characterize parallelism of two compact linear operators between normed linear spaces $\mathbb{X} $ and $\mathbb{Y}$, assuming $\mathbb{X}$ to be reflexive. We also characterize parallelism of two bounded linear operators between normed linear spaces $\mathbb{X} $ and $\mathbb{Y}.$ We investigate parallelism and approximate parallelism in the space of bounded linear operators defined on a Hilbert space. Using the characterization of operator parallelism, we study Birkhoff-James orthogonality in the space of compact linear operators as well as bounded linear operators. Finally, we introduce the concept of semi-rotund points (semi-rotund spaces) which generalizes the notion of exposed points (strictly convex spaces). We further study semi-rotund operators and prove that $\mathbb{B}(\mathbb{X},\mathbb{Y})$ is a semi-rotund space which is not strictly convex, if $\mathbb{X},\mathbb{Y}$ are finite-dimensional Banach spaces and $\mathbb{Y}$ is strictly convex.


title: Bogolon-mediated electron capture by impurities in hybrid Bose-Fermi systems

abstract: We investigate the processes of electron capture by a Coulomb impurity center residing in a hybrid system consisting of spatially separated two-dimensional layers of electron and Bose-condensed dipolar exciton gases coupled via the Coulomb forces. We calculate the probability of the electron capture accompanied by the emission of a single Bogoliubov excitation (bogolon), similar to regular phonon-mediated scattering in solids. Further, we study the electron capture mediated by the emission of a pair of bogolons in a single capture event and show that these processes not only should be treated in the same order of the perturbation theory, but also they give more important contribution than single bogolon-mediated capture, in contrast with regular phonon scattering.


title: The hadronic vacuum polarisation contributions to the muon $g-2$

abstract: The hadronic vacuum polarisation contributions to the anomalous magnetic moment of the muon, $a_{\mu}^{\rm had, VP}$ have been re-evaluated from the combination of $e^+e^-\rightarrow {\rm hadrons}$ cross section data. Focus has been placed on the development of a new data combination method, which fully incorporates all correlated statistical and systematic uncertainties in a bias free approach. Using these combined data have resulted in estimates of the hadronic vacuum polarisation contributions to $g-2$ of the muon of $a_{\mu}^{\rm had, \, LO \, VP} = (693.27 \pm 2.46)\times 10^{-10}$ and $a_{\mu}^{\rm had, \, NLO \, VP} = (-9.82 \pm 0.04)\times 10^{-10}$. The new estimate for the Standard Model prediction is found to be $a_{\mu}^{\rm SM} = (11\ 659 \ 182.05 \pm 3.56) \times 10^{-10}$, which is $3.7\sigma$ below the current experimental measurement. In addition, the prediction for the hadronic contribution to the QED coupling at the $Z$ boson mass has been calculated to be $\Delta\alpha_{\rm had}^{(5)}(M_Z^2)= (276.11 \pm 1.11)\times 10^{-4}$, resulting in $\alpha^{-1}(M_Z^2) = 128.946 \pm 0.015$.


title: Dark-bright soliton pairs: bifurcations and collisions

abstract: The statics, stability and dynamical properties of dark-bright soliton pairs are investigated motivated by applications in a homogeneous system of two-component repulsively interacting Bose-Einstein condensate. One of the intra-species interaction coefficients is used as the relevant parameter controlling the deviation from the integrable Manakov limit. Two different families of stationary states are identified consisting of dark-bright solitons that are either antisymmetric (out-of-phase) or asymmetric (mass imbalanced) with respect to their bright soliton. Both of the above dark-bright configurations coexist at the integrable limit of equal intra- and inter-species repulsions and are degenerate in that limit. However, they are found to bifurcate from it in a transcritical bifurcation. The latter interchanges the stability properties of the bound dark-bright pairs rendering the antisymmetric states unstable and the asymmetric ones stable past the associated critical point (and vice versa before it). Finally, on the dynamical side, it is found that large kinetic energies and thus rapid soliton collisions are essentially unaffected by the intra-species variation, while cases involving near equilibrium states or breathing dynamics are significantly modified under such a variation.

categories: nlin.PS cond-mat.quant-gas physics.atom-ph quant-ph


title: Persistence of Small Noise and Random initial conditions

abstract: The effect of small noise in a smooth dynamical system is negligible on any finite time interval. Here we study situations when it persists on intervals increasing to infinity. Such asymptotic regime occurs when the system starts from initial condition, sufficiently close to an unstable fixed point. In this case, under appropriate scaling, the trajectory converges to solution of the unperturbed system, started from a certain {\em random} initial condition. In this paper we consider the case of one dimensional diffusions on the positive half line, which often arise as scaling limits in population dynamics.


title: Finding Global Optima in Nonconvex Stochastic Semidefinite Optimization with Variance Reduction

abstract: There is a recent surge of interest in nonconvex reformulations via low-rank factorization for stochastic convex semidefinite optimization problem in the purpose of efficiency and scalability. Compared with the original convex formulations, the nonconvex ones typically involve much fewer variables, allowing them to scale to scenarios with millions of variables. However, it opens a new challenge that under what conditions the nonconvex stochastic algorithms may find the global optima effectively despite their empirical success in applications. In this paper, we provide an answer that a stochastic gradient descent method with variance reduction, can be adapted to solve the nonconvex reformulation of the original convex problem, with a \textit{global linear convergence}, i.e., converging to a global optimum exponentially fast, at a proper initial choice in the restricted strongly convex case. Experimental studies on both simulation and real-world applications on ordinal embedding are provided to show the effectiveness of the proposed algorithms.


title: Persistence and Stochastic Periodicity in the Intensity Dynamics of a Fiber Laser During the Transition to Optical Turbulence

abstract: Many natural systems display transitions among different dynamical regimes, which are difficult to identify when the data is noisy and high dimensional. A technologically relevant example is a fiber laser, which can display complex dynamical behaviors that involve nonlinear interactions of millions of cavity modes. Here we study the laminar-turbulence transition that occurs when the laser pump power is increased. By applying various data analysis tools to empirical intensity time series we characterize their persistence and demonstrate that at the transition temporal correlations can be precisely represented by a surprisingly simple model.


title: Characterization of the superconducting phase in tellurium hydride at high pressure

abstract: At present, hydrogen-based compounds constitute one of the most promising classes of materials for applications as a phonon-mediated high-temperature superconductors. Herein, the behavior of the superconducting phase in tellurium hydride (HTe) at high pressure ($p=300$ GPa) is analyzed in details, by using the isotropic Migdal-Eliashberg equations. The chosen pressure conditions are considered here as a case study which corresponds to the highest critical temperature value ($T_{c}$) in the analyzed material, as determined within recent density functional theory simulations. It is found that the Migdal-Eliashberg formalism, which constitutes a strong-coupling generalization of the Bardeen-Cooper-Schrieffer (BCS) theory, predicts that the critical temperature value ($T_{c}=52.73$ K) is higher than previous estimates of the McMillan formula. Further investigations show that the characteristic dimensionless ratios for the the thermodynamic critical field, the specific heat for the superconducting state, and the superconducting band gap exceeds the limits of the BCS theory. In this context, also the effective electron mass is not equal to the bare electron mass as provided by the BCS theory. On the basis of these findings it is predicted that the strong-coupling and retardation effects play pivotal role in the superconducting phase of HTe at 300 GPa, in agreement with similar theoretical estimates for the sibling hydrogen and hydrogen-based compounds. Hence, it is suggested that the superconducting state in HTe cannot be properly described within the mean-field picture of the BCS theory.


title: Evidence for two spin-glass transitions with magnetoelastic and magnetoelectric couplings in the multiferroic (Bi$_{1-x}$Ba$_x$)(Fe$_{1-x}$Ti$_x$)O$_3$ system

abstract: For disordered Heisenberg systems with small single ion anisotropy, two spin glass transitions below the long range ordered phase transition temperature has been predicted theoretically for compositions close to the percolation threshold. Experimental verification of these predictions is still controversial for conventional spin glasses. We show that multiferroic spin glass systems can provide a unique platform for verifying these theoretical predictions via a study of change in magnetoelastic and magnetoelectric couplings, obtained from an analysis of diffraction data, at the spin glass transition temperatures. Results of macroscopic and microscopic (x-ray and neutron scattering) measurements are presented on disordered BiFeO3, a canonical Heisenberg system with small single ion anisotropy, which reveal appearance of two spin glass phases SG1 and SG2 in coexistence with the LRO phase below the A-T and G-T lines. It is shown that the temperature dependence of the integrated intensity of the antiferromagnetic peak shows dips with respect to the Brillouin function behaviour around the SG1 and SG2 transition temperatures. The ferroelectric polarisation changes significantly at the two spin glass transition temperatures. These results, obtained using microscopic techniques, clearly demonstrate that the SG1 and SG2 transitions occur on the same magnetic sublattice and are intrinsic to the system. We also construct a phase diagram showing all the magnetic phases in BF-xBT system. While our results on the two spin glass transitions support the theoretical predictions, it also raises several open questions which need to be addressed by revisiting the existing theories of spin glass transitions by taking into account the effect of magnetoelastic and magnetoelectric couplings as well as electromagnons.


title: Quantum phase transitions and the degree of nonidentity in the system with two different species of vector bosons

abstract: We address the system with two species of vector bosons in an optical lattice. In addition to the the standard parameters characterizing such a system, we are dealing here with the "degree of atomic nonidentity", manifesting itself in the difference of tunneling amplitudes and on-site Coulomb interactions. We obtain a cascade of quantum phase transitions occurring with the increase in the degree of atomic nonidentity. In particular, we show that the phase diagram for strongly distinct atoms is qualitatively different from that for (nearly) identical atoms considered earlier. The resulting phase diagrams evolve from the images similar to the "J. Mir\'o-like paintings" to "K. Malewicz-like" ones.


title: Tetramers of two heavy and two light bosons

abstract: This article considers the bound states of two heavy and two light bosons, when a short-range force attracts the bosons of different mass, and a short-range force repel the light bosons. The existence of such four-body bound states results from the competition between these two forces. For a given strength of the attraction, the critical strength of the repulsion necessary to unbind the four particles is calculated. This study is motivated by the experimental realisation of impurity atoms immersed in an atomic Bose-Einstein condensate, and aims at determining in which regime only one boson contributes to binding two impurities.


title: Properties of Decameter IIIb-III Pairs

abstract: A large number of Type IIIb-III pairs, in which the first component is a Type IIIb burst and the second one is a Type III burst, are often recorded during decameter Type III burst storms. From the beginning of their observation, the question of whether the components of these pairs are the first and the second harmonics of radio emission or not has remained open. We discuss properties of decameter IIIb-III pairs in detail to answer this question. The components of these pairs, Type IIIb bursts and Type III bursts, have essentially different durations and polarizations. At the same time their frequency drift rates are rather close, provided that the drift rates of Type IIIb bursts are a little larger those of Type III bursts at the same frequency. Frequency ratios of the bursts at the same moment are close to two. This points at a harmonic connection of the components in IIIb-III pairs. At the same time there was a serious difficulty, namely why the first harmonic had fine frequency structure in the form of striae and the second harmonic did not have it. Recently Loi, Cairns, and Li ( Astrophys. J. 790, 67, 2014) succeeded in solving this problem. The physical aspects of observational properties of decameter IIIb-III pairs are discussed and pros and cons of harmonic character of Type IIIb bursts and Type III bursts in IIIb-III pairs are presented. We conclude that practically all properties of the IIIb-III pair components can be understood in the framework of the harmonic relation of the components of the IIIb-III pairs.


title: Multiple boundary representations of $\lambda$-harmonic functions on trees

abstract: We consider a countable tree $T$, possibly having vertices with infinite degree, and an arbitrary stochastic nearest neighbour transition operator $P$. We provide a boundary integral representation for general eigenfunctions of $P$ with eigenvalue $\lambda \in \mathbb{C}$, under the condition that the oriented edges can be equipped with complex-valued weights satisfying three natural axioms. These axioms guarantee that one can construct a $\lambda$-Poisson kernel. The boundary integral is with respect to distributions, that is, elements in the dual of the space of locally constant functions. Distributions are interpreted as finitely additive complex measures. In general, they do not extend to $\sigma$-additive measures: for this extension, a summability condition over disjoint boundary arcs is required. Whenever $\lambda$ is in the resolvent of $P$ as a self-adjoint operator on a naturally associated $\ell^2$-space and the diagonal elements of the resolvent (`Green function') do not vanish at $\lambda$, one can use the ordinary edge weights corresponding to the Green function and obtain the ordinary $\lambda$-Martin kernel. We then consider the case when $P$ is invariant under a transitive group action. In this situation, we study the phenomenon that in addition to the $\lambda$-Martin kernel, there may be further choices for the edge weights which give rise to another $\lambda$-Poisson kernel with associated integral representations. In particular, we compare the resulting distributions on the boundary. The material presented here is closely related to the contents of our `companion' paper arXiv:1802.01976


title: J0023+0307: A mega metal-poor dwarf star from SDSS/BOSS

abstract: Only a handful of stars have been identified with an iron abundance [Fe/H]<-5, and only one at [Fe/H]<-7. These stars have very large carbon-to-iron ratios, with A(C)~7.0, most likely due to fallback in core-collapse supernovae, which makes their total metallicity Z much higher than their iron abundances. The failure to find population III stars, those with no metals, has been interpreted, with support from theoretical modeling, as the result of a top-heavy initial mass function. With zero or very low metal abundance limiting radiative cooling, the formation of low-mass stars could be inhibited. Currently, the star SDSS J1029+1729 sets the potential metallicity threshold for the formation of low-mass stars at log Z/Z_0~-5. We have identified SDSS J0023+0307, a primitive star with T_eff =6188+-84 K, and log g=4.9+-0.5, an upper limit [Fe/H]< -6.6, and a carbon abundance A(C)< 6.3. In our quest to push down the metallicity threshold we find J0023+0307 to be one of the two most iron-poor stars known, and it exhibits less carbon that most of stars at [Fe/H]< -5.


title: Sedimentation Efficiency of Condensation Clouds in Substellar Atmospheres

abstract: Condensation clouds in substellar atmospheres have been widely inferred from spectra and photometric variability. Up until now, their horizontally averaged vertical distribution and mean particle size have been largely characterized using models, one of which is the eddy diffusion-sedimentation model from Ackerman & Marley (2001) that relies on a sedimentation efficiency parameter, $f_{\rm sed}$, to determine the vertical extent of clouds in the atmosphere. However, the physical processes controlling the vertical structure of clouds in substellar atmospheres are not well understood. In this work, we derive trends in $f_{\rm sed}$ across a large range of eddy diffusivities ($K_{zz}$), gravities, material properties, and cloud formation pathways by fitting cloud distributions calculated by a more detailed cloud microphysics model. We find that $f_{\rm sed}$ is dependent on $K_{zz}$, but not gravity, when $K_{zz}$ is held constant. $f_{\rm sed}$ is most sensitive to the nucleation rate of cloud particles, as determined by material properties like surface energy and molecular weight. High surface energy materials form fewer, larger cloud particles, leading to large $f_{\rm sed}$ ($>$1), and vice versa for materials with low surface energy. For cloud formation via heterogeneous nucleation, $f_{\rm sed}$ is sensitive to the condensation nuclei flux and radius, connecting cloud formation in substellar atmospheres to the objects' formation environments and other atmospheric aerosols. These insights could lead to improved cloud models that help us better understand substellar atmospheres. For example, we demonstrate that $f_{\rm sed}$ could increase with increasing cloud base depth in an atmosphere, shedding light on the nature of the brown dwarf L/T transition.