title: Could Solar Radiation Pressure Explain 'Oumuamua's Peculiar Acceleration?

abstract: `Oumuamua (1I/2017 U1) is the first object of interstellar origin observed in the Solar System. Recently, \citet{Micheli2018} reported that `Oumuamua showed deviations from a Keplerian orbit at a high statistical significance. The observed trajectory is best explained by an excess radial acceleration $\Delta a \propto r^{-2}$, where $r$ is the distance of `Oumuamua from the Sun. Such an acceleration is naturally expected for comets, driven by the evaporating material. However, recent observational and theoretical studies imply that `Oumuamua is not an active comet. We explore the possibility that the excess acceleration results from Solar radiation pressure. The required mass-to-area ratio is $(m/A)\approx 0.1$ g cm$^{-2}$. For a thin sheet this requires a thickness of $\approx 0.3-0.9$ mm. We find that although extremely thin, such an object would survive an interstellar travel over Galactic distances of $\sim 5$ kpc, withstanding collisions with gas and dust-grains as well as stresses from rotation and tidal forces. We discuss the possible origins of such an object. Our general results apply to any light probes designed for interstellar travel.


title: Empirical Evaluation of Contextual Policy Search with a Comparison-based Surrogate Model and Active Covariance Matrix Adaptation

abstract: Contextual policy search (CPS) is a class of multi-task reinforcement learning algorithms that is particularly useful for robotic applications. A recent state-of-the-art method is Contextual Covariance Matrix Adaptation Evolution Strategies (C-CMA-ES). It is based on the standard black-box optimization algorithm CMA-ES. There are two useful extensions of CMA-ES that we will transfer to C-CMA-ES and evaluate empirically: ACM-ES, which uses a comparison-based surrogate model, and aCMA-ES, which uses an active update of the covariance matrix. We will show that improvements with these methods can be impressive in terms of sample-efficiency, although this is not relevant any more for the robotic domain.


title: Black hole entropy and quantum mechanics

abstract: We give a brief overview of black hole entropy, covering a few main developments since Bekenstein's original proposal


title: Exoplanet Transits with Next-Generation Radio Telescopes

abstract: Nearly everything we know about extrasolar planets to date comes from optical astronomy. While exoplanetary aurorae are predicted to be bright at low radio frequencies (< 1 GHz), we consider the effect of an exoplanet transit on radio emission from the host star. As radio emission from solar-like stars is concentrated in active regions, a planet occulting a starspot can cause a disproportionately deep transit which should be detectable with major radio arrays currently under development, such as the Square Kilometre Array (SKA). We calculate the radiometric sensitivity of the SKA stages and components, finding that SKA2-Mid can expect to detect transits around the very nearest solar-like stars and many cool dwarfs. The shape of this radiometric light curve will be affected by scintillation and lensing from the planet's magnetosphere and thereby encode magnetospheric parameters. Furthermore, these transits will also probe the distribution of stellar activity across a star's surface, and will help scrub out contamination from stellar activity on exoplanet transmission spectra and radial velocity spectra. This radio window on exoplanets and their host stars is therefore a valuable complement to existing optical tools.


title: Charging of a Single InAs QD with Electrically-Injected Holes using a Lateral Electric Field

abstract: InAs/GaAs quantum dots (QDs) and quantum dot molecules (QDMs) are self-assembled semiconductor nanostructures that can trap a single electron or hole with well-defined spin projections. QDs and QDMs have excellent optical properties and have long been of interest for incorporation into quantum optoelectronic devices ranging from single photon sources to multi-bit quantum computers. The properties of single QDs, or carriers confined within those QDs, can be tuned by external electric fields, which provides an important tool for the development of scalable and tunable devices. Deterministic charging of a QD with a single electron or hole is an important tool for quantum devices and is well-established under the application of growth-direction electric fields in a diode structure. Here, we report a new charging mechanism for a single QD in a 3-electrode device that does not contain a vertical diode and can be used to control electric field profiles in two dimensions. We fabricate the device with E-beam lithography and characterize photoluminescence from single QDs under different bias configurations. Using a combination of experimental data and COMSOL band structure calculations, we explain how the charging originates in the electrical-injection of holes induced by lateral electric fields. We discuss the applications for this device and the potential for full 2-D electric field control of a single QD and QDM.


title: Conditioning and backward errors of eigenvalues of homogeneous matrix polynomials under M\"{o}bius transformations

abstract: M\"{o}bius transformations have been used in numerical algorithms for computing eigenvalues and invariant subspaces of structured generalized and polynomial eigenvalue problems (PEPs). These transformations convert problems with certain structures arising in applications into problems with other structures and whose eigenvalues and invariant subspaces are easily related to the ones of the original problem. Thus, an algorithm that is efficient and stable for some particular structure can be used for solving efficiently another type of structured problem via an adequate M\"{o}bius transformation. A key question in this context is whether these transformations may change significantly the conditioning of the problem and the backward errors of the computed solutions, since, in that case, their use may lead to unreliable results. We present the first general study on the effect of M\"{o}bius transformations on the eigenvalue condition numbers and backward errors of approximate eigenpairs of PEPs. By using the homogeneous formulation of PEPs, we are able to obtain two clear and simple results. First, we show that, if the matrix inducing the M\"{o}bius transformation is well conditioned, then such transformation approximately preserves the eigenvalue condition numbers and backward errors when they are defined with respect to perturbations of the matrix polynomial which are small relative to the norm of the polynomial. However, if the perturbations in each coefficient of the matrix polynomial are small relative to the norm of that coefficient, then the corresponding eigenvalue condition numbers and backward errors are preserved approximately by the M\"{o}bius transformations induced by well-conditioned matrices only if a penalty factor, depending on those coefficients, is moderate. It is important to note that these simple results are no longer true if a non-homogeneous formulation is used.


title: Balanced ideals and domains of discontinuity of Anosov representations

abstract: We consider the action of Anosov subgroups of a semi-simple Lie group on the associated flag manifolds. A systematic approach to construct cocompact domains of discontinuity for this action was given by Kapovich, Leeb and Porti in arXiv:1306.3837. For $\Delta$-Anosov representations, we prove that every cocompact domain of discontinuity arises from this construction, up to a few exceptions in low rank. Then we compute which flag manifolds admit these domains and, in some cases, the number of domains. We also find a new compactification for locally symmetric spaces arising from maximal representations into $\mathrm{Sp}(4n+2, \mathbb{R})$.


title: Parsing Coordination for Spoken Language Understanding

abstract: Typical spoken language understanding systems provide narrow semantic parses using a domain-specific ontology. The parses contain intents and slots that are directly consumed by downstream domain applications. In this work we discuss expanding such systems to handle compound entities and intents by introducing a domain-agnostic shallow parser that handles linguistic coordination. We show that our model for parsing coordination learns domain-independent and slot-independent features and is able to segment conjunct boundaries of many different phrasal categories. We also show that using adversarial training can be effective for improving generalization across different slot types for coordination parsing.


title: Automatic Identification and Ranking of Emergency Aids in Social Media Macro Community

abstract: Online social microblogging platforms including Twitter are increasingly used for aiding relief operations during disaster events. During most of the calamities that can be natural disasters or even armed attacks, non-governmental organizations look for critical information about resources to support effected people. Despite the recent advancement of natural language processing with deep neural networks, retrieval and ranking of short text becomes a challenging task because a lot of conversational and sympathy content merged with the critical information. In this paper, we address the problem of categorical information retrieval and ranking of most relevance information while considering the presence of short-text and multilingual languages that arise during such events. Our proposed model is based on the formation of embedding vector with the help of textual and statistical preprocessing, and finally, entire training 2,100,000 vectors were normalized using feed-forward neural network for need and availability tweets. Another important contribution of this paper lies in novel weighted Ranking Key algorithm based on top five general terms to rank the classified tweets in most relevance with classification. Lastly, we test our model on Nepal Earthquake dataset (contains short text and multilingual language tweets) and achieved 6.81% of mean average precision on 5,250,000 unlabeled embedding vectors of disaster relief tweets.


title: Information Bottleneck Methods for Distributed Learning

abstract: We study a distributed learning problem in which Alice sends a compressed distillation of a set of training data to Bob, who uses the distilled version to best solve an associated learning problem. We formalize this as a rate-distortion problem in which the training set is the source and Bob's cross-entropy loss is the distortion measure. We consider this problem for unsupervised learning for batch and sequential data. In the batch data, this problem is equivalent to the information bottleneck (IB), and we show that reduced-complexity versions of standard IB methods solve the associated rate-distortion problem. For the streaming data, we present a new algorithm, which may be of independent interest, that solves the rate-distortion problem for Gaussian sources. Furthermore, to improve the results of the iterative algorithm for sequential data we introduce a two-pass version of this algorithm. Finally, we show the dependency of the rate on the number of samples $k$ required for Gaussian sources to ensure cross-entropy loss that scales optimally with the growth of the training set.


title: Non-isobaric Thermal Instability

abstract: Multiphase media have very complex structure and evolution. Accurate numerical simulations are necessary to make advances in our understanding of this rich physics. Because simulations can capture both the linear and nonlinear evolution of perturbations with a relatively wide range of sizes, it is important to thoroughly understand the stability of condensation and acoustic modes between the two extreme wavelength limits of isobaric and isochoric instability as identified by Field (1965). Partially motivated by a recent suggestion that large non-isobaric clouds can `shatter' into tiny cloudlets, we revisit the linear theory to survey all possible regimes of thermal instability. We uncover seven regimes in total, one of which allows three unstable condensation modes. Using the code Athena++, we determine the numerical requirements to properly evolve small amplitude perturbations of the entropy mode into the nonlinear regime. Our 1D numerical simulations demonstrate that for a typical AGN cooling function, the nonlinear evolution of a single eigenmode in an isobarically unstable plasma involves increasingly larger amplitude oscillations in cloud size, temperature and density as the wavelength increases. Such oscillations are the hallmark behavior of non-isobaric multiphase gas dynamics and may be observable as correlations between changes in brightness and the associated periodic redshifts and blueshifts in systems that can be spatially resolved. Intriguingly, we discuss regimes and derive characteristic cloud sizes for which the saturation process giving rise to these oscillations can be so energetic that the cloud may indeed break apart. However, we dub this process `splattering' instead of `shattering', as it is a different fragmentation mechanism triggered when the cloud suddenly `lands' on the stable cold branch of the equilibrium curve.


title: The maximum luminosity of fast radio bursts

abstract: Under the assumption that fast radio bursts (FRBs) are from coherent curvature emission powered by the dissipation of magnetic energy in the magnetosphere of neutron stars, we predict a maximum isotropic equivalent luminosity of (~2x10^{47} erg/s) min({\rho}^2, B{\rho}^{4/3}{\nu}^{-2/3}), where {\rho} is the curvature radius of the magnetic field lines near the source region in units of 10 km, B is the local magnetic field strength in units of 10^{16} G, and {\nu} is the FRB wave frequency in units of GHz. This is because the electric field responsible for accelerating the radiating particles becomes close to the quantum critical field strength and is then quickly shielded by Schwinger pairs within a nano-second. Future observations should be able to confirm this cut-off luminosity and hence provide a unique probe of the source location and magnetic field strength. We discuss how the maximum luminosity affects the observed flux distributions for repeating bursts from the same object and for the entire FRB population.


title: Pulsating low-mass white dwarfs in the frame of new evolutionary sequences VI. Thin H-envelope sequences and asteroseismology of ELMV stars revisited

abstract: Some low-mass white-dwarf (LMWD) stars with H atmospheres show long-period g-mode pulsations, comprising the class of pulsating WDs called extremely low-mass variable (ELMV) stars. It is generally believed that these stars have thick H envelopes. However, from stellar evolution considerations, the existence of LMWDs with thin H envelopes is also possible. We present a thorough asteroseismological analysis of ELMV stars based on a complete set of fully evolutionary models that represents low-mass He-core WD stars harboring a range of H envelope thicknesses. Although there are currently nine ELMVs, here we only focus on those that exhibit more than three periods and whose periods do not show significant uncertainties. We considered g-mode adiabatic pulsation periods for low-mass He-core WD models with $M_*$ in the range [0.1554-0.4352]$M_{\odot}$, $T_{\rm eff}$ in the range [6000-10000]K, and H envelope thicknesses in the range -5.8<log($M_{\rm H}/M_*$)<-1.7. We explore the effects of employing different H-envelope thicknesses on the adiabatic pulsation properties of low-mass He-core WD models, and perform period-to-period fits to ELMVs to search for a representative model. We found that the mode-trapping effects of g modes depend sensitively on $M_{\rm H}$, with the trapping cycle and trapping amplitude larger for thinner H envelopes. Also, the asymptotic period spacing is longer for thinner H envelopes. Finally, we found asteroseismological models (when possible) characterized by canonical (thick) and by thin H envelope, with $T_{\rm eff}$ and $M_*$ in agreement with the spectroscopic determinations. The fact that we have found asteroseismological solutions with H envelopes thinner than canonical gives a clue of the possible scenario of formation of these stars. Indeed, in the light of our results, some of these stars could have been formed by binary evolution through unstable mass loss.


title: Machine learning determination of dynamical parameters: The Ising model case

abstract: We train a set of Restricted Boltzmann Machines (RBMs) on one- and two-dimensional Ising spin configurations at various values of temperature, generated using Monte Carlo simulations. We validate the training procedure by monitoring several estimators, including measurements of the log-likelihood, with the corresponding partition functions estimated using annealed importance sampling. The effects of various choices of hyper-parameters on training the RBM are discussed in detail, with a generic prescription provided. Finally, we present a closed form expression for extracting the values of couplings, for every $n$-point interaction between the visible nodes of an RBM, in a binary system such as the Ising model. We aim at using this study as the foundation for further investigations of less well-known systems.


title: Riemann problems and dispersive shocks in self-focusing media

abstract: The dynamical behavior resulting from an initial discontinuity in focusing media is investigated using a combination of numerical simulations and Whitham modulation theory for the focusing nonlinear Schrodinger equation. Initial conditions with a jump in either or both the amplitude and the local wavenumber are considered. It is shown analytically and numerically that the space-time plane divides into expanding domains in which the solution is described by a slow modulation of genus-zero, genus-one or genus-two solutions, their precise arrangement depending on the specifics of the initial datum


title: Stability-certified reinforcement learning: A control-theoretic perspective

abstract: We investigate the important problem of certifying stability of reinforcement learning policies when interconnected with nonlinear dynamical systems. We show that by regulating the input-output gradients of policies, strong guarantees of robust stability can be obtained based on a proposed semidefinite programming feasibility problem. The method is able to certify a large set of stabilizing controllers by exploiting problem-specific structures; furthermore, we analyze and establish its (non)conservatism. Empirical evaluations on two decentralized control tasks, namely multi-flight formation and power system frequency regulation, demonstrate that the reinforcement learning agents can have high performance within the stability-certified parameter space, and also exhibit stable learning behaviors in the long run.


title: Critique of Critical Shear Crack Theory (CSCT) for fib Model Code Articles on Shear Strength and Size Effect of RC Beams

abstract: The size effect of Muttoni et al.'s CSCT is shown to be very close and asymptotically identical to the energetic size effect law (SEL) extensively verified experimentally and theoretically (due to appear in the 2019 ACI Code for both beam shear and punching). However, the CSCT derivation and calculation procedure obfuscates the mechanics of failure. It is shown to rest on six scientifically untenable hypotheses, which makes CSCT untrustworthy outside the testing range. It would have to be taught to students as a dogma. These conclusions are supported by experimentally calibrated finite element simulations of crack path and width, stress distributions and localizations during failure, and strain energy release. The simulations also show the CSCT to be incompatible with the 'strut-and-tie' model modernized to include the size effect in the compression strut. Finally, further deficiencies are pointed out for the Modified Compression Field Theory (MCFT), currently embedded in the Model Code.


title: Efficient Distributed Hessian Free Algorithm for Large-scale Empirical Risk Minimization via Accumulating Sample Strategy

abstract: In this paper, we propose a Distributed Accumulated Newton Conjugate gradiEnt (DANCE) method in which sample size is gradually increasing to quickly obtain a solution whose empirical loss is under satisfactory statistical accuracy. Our proposed method is multistage in which the solution of a stage serves as a warm start for the next stage which contains more samples (including the samples in the previous stage). The proposed multistage algorithm reduces the number of passes over data to achieve the statistical accuracy of the full training set. Moreover, our algorithm in nature is easy to be distributed and shares the strong scaling property indicating that acceleration is always expected by using more computing nodes. Various iteration complexity results regarding descent direction computation, communication efficiency and stopping criteria are analyzed under convex setting. Our numerical results illustrate that the proposed method outperforms other comparable methods for solving learning problems including neural networks.


title: Complete IRAC mapping of the CFHTLS-DEEP, MUSYC AND NMBS-II FIELDS

abstract: The IRAC mapping of the NMBS-II fields program is an imaging survey at 3.6 and 4.5$\mu$m with the Spitzer Infrared Array Camera (IRAC). The observations cover three Canada-France-Hawaii Telescope Legacy Survey Deep (CFHTLS-D) fields, including one also imaged by AEGIS, and two MUSYC fields. These are then combined with archival data from all previous programs into deep mosaics. The resulting imaging covers a combined area of about 3 $deg^2$, with at least $\sim$2 hr integration time for each field. In this work, we present our data reduction techniques and document the resulting coverage maps at 3.6 and 4.5$\mu$m. All of the images are W-registered to the reference image, which is either the z-band stack image of the 25\% best seeing images from the CFHTLS-D for CFHTLS-D1, CFHTLS-D3, and CFHTLS-D4, or the K-band images obtained at the Blanco 4-m telescope at CTIO for MUSYC1030 and MUSYC1255. We make all images and coverage maps described herein publicly available via the Spitzer Science Center.


title: Neural Network-Based Approach to Phase Space Integration

abstract: Monte Carlo methods are widely used in particle physics to integrate and sample probability distributions (differential cross sections or decay rates) on multi-dimensional phase spaces. We present a Neural Network (NN) algorithm optimized to perform this task. The algorithm has been applied to several examples of direct relevance for particle physics, including situations with non-trivial features such as sharp resonances and soft/collinear enhancements. Excellent performance has been demonstrated in all examples, with the properly trained NN achieving unweighting efficiencies of between 30% and 75%. In contrast to traditional Monte Carlo algorithms such as VEGAS, the NN-based approach does not require that the phase space coordinates be aligned with resonant or other features in the cross section.


title: Computational modeling of decay and hybridization in superconducting circuits

abstract: We present a framework for modeling superconducting circuits that integrates classical microwave analysis with circuit quantization. Our framework enables the calculation of the lossy eigenmodes of superconducting circuits, and we demonstrate the method by analyzing several circuits relevant to multiplexed, Purcell filtered transmon readout architectures. We show that the transmon relaxation times obtained by our method agree with the established approximation $T_1 \approx C/\text{Re}[Y(i\omega_q)]$ away from environmental resonances and do not vanish on resonance. We also show that the hybridization of the modes in the readout circuit is highly sensitive to the bandwidth of the Purcell filter.


title: VanQver: The Variational and Adiabatically Navigated Quantum Eigensolver

abstract: The accelerated progress in manufacturing noisy intermediate-scale quantum (NISQ) computing hardware has opened the possibility of exploring its application in transforming approaches to solving computationally challenging problems. The important limitations common among all NISQ computing technologies are the absence of error correction and the short coherence time, which limit the computational power of these systems. Shortening the required time of a single run of a quantum algorithm is essential for reducing environment-induced errors and for the efficiency of the computation. We have investigated the ability of a variational version of adiabatic quantum computation (AQC) to generate an accurate state more efficiently compared to existing adiabatic methods. The standard AQC method uses a time-dependent Hamiltonian, connecting the initial Hamiltonian with the final Hamiltonian. In the current approach, a navigator Hamiltonian is introduced which has a non-zero amplitude only in the middle of the annealing process. Both the initial and navigator Hamiltonians are determined using variational methods. A hermitian cluster operator, inspired by coupled-cluster theory and truncated to single and double excitations/de-excitations, is used as a navigator Hamiltonian. A comparative study of our variational algorithm (VanQver) with that of standard AQC, starting with a Hartree--Fock Hamiltonian, is presented. The results indicate that the introduction of the navigator Hamiltonian significantly improves the annealing time required to achieve chemical accuracy by two to three orders of magnitude. The efficiency of the method is demonstrated in the ground-state energy estimation of molecular systems, namely, H$_2$, P4, and LiH.


title: A new lattice invariant for lattices in totally disconnected locally compact groups

abstract: We introduce and explore a natural rank for totally disconnected locally compact groups called the bounded conjugacy rank. This rank is shown to be a lattice invariant for lattices in sigma compact totally disconnected locally compact groups; that is to say, for a given sigma compact totally disconnected locally compact group, some lattice has bounded conjugacy rank n if and only if every lattice has bounded conjugacy rank n. Several examples are then presented.


title: Investigation of re-entrant relaxor behaviour in lead cobalt niobate ceramic

abstract: The temperature dependent dielectric properties revealed re-entrant relaxor behaviour (Tm ~130 K and 210 K for 1 kHz) below a high temperature diffused phase transition, Tc ~270 K in lead cobalt niobate (PCN). Multiple positive/negative magnetodielectric effect and deviation from straight line at ~130 K is observed in temperature dependence of inverse susceptibility, which depicts origin of frustration. Microstructure examination depicts closely packed grains with grain size ~8-10 microm and XRD pattern revealed single phase pseudo cubic crystal structure having Pm3m symmetry with lattice constant ~4.0496(2) {\AA}. Rietveld Refinement on XRD data yields larger value of thermal parameters, implying Pb and O are disordered along <111> and <110> directions respectively. Observation of A1g (780 cm-1) mode in Raman spectroscopy and F-spot in SAED pattern along <110> unit axis in TEM suggests presence of nano scale 1:1 Co and Nb non-stoichiometric chemical ordering (CORs), akin to lead magnesium niobate (PMN). K-edge XANES spectra reveals the presence of cobalt in two oxidation states (Co2+ and Co3+); whereas, niobium exists in Nb3+ state. Therefore, these local-average structural properties suggest chemical, structural and spatial heterogeneities. Such multiple heterogeneities are believed to play a crucial role in producing re-entrant relaxor behaviour.


title: Learning and Interpreting Multi-Multi-Instance Learning Networks

abstract: We introduce an extension of the multi-instance learning problem where examples are organized as nested bags of instances (e.g., a document could be represented as a bag of sentences, which in turn are bags of words). This framework can be useful in various scenarios, such as text and image classification, but also supervised learning over graphs. As a further advantage, multi-multi instance learning enables a particular way of interpreting predictions and the decision function. Our approach is based on a special neural network layer, called bag-layer, whose units aggregate bags of inputs of arbitrary size. We prove theoretically that the associated class of functions contains all Boolean functions over sets of sets of instances and we provide empirical evidence that functions of this kind can be actually learned on semi-synthetic datasets. We finally present experiments on text classification, on citation graphs, and social graph data, which show that our model obtains competitive results with respect to accuracy when compared to other approaches such as convolutional networks on graphs, while at the same time it supports a general approach to interpret the learnt model, as well as explain individual predictions.


title: Noise Sensitivity of Local Descriptors vs ConvNets: An application to Facial Recognition

abstract: The Local Binary Patterns (LBP) is a local descriptor proposed by Ojala et al to discriminate texture due to its discriminative power. However, the LBP is sensitive to noise and illumination changes. Consequently, several extensions to the LBP such as Median Binary Pattern (MBP) and methods such as Local Directional Pattern (LDP) have been proposed to address its drawbacks. Though studies by Zhou et al, suggest that the LDP exhibits poor performance in presence of random noise. Recently, convolution neural networks (ConvNets) were introduced which are increasingly becoming popular for feature extraction due to their discriminative power. This study aimed at evaluating the sensitivity of ResNet50, a ConvNet pre-trained model and local descriptors (LBP and LDP) to noise using the Extended Yale B face dataset with 5 different levels of noise added to the dataset. In our findings, it was observed that despite adding different levels of noise to the dataset, ResNet50 proved to be more robust than the local descriptors (LBP and LDP).


title: The issue with the initial state in quantum mechanics

abstract: In the conventional formulation of quantum mechanics, the initial description is given only for the physical system under study. It factors out the state for the experimenter. We argue that such description is incomplete and can lead to statements which can in theory be meaningless. We propose that within a complete description, the initial state must include the state of the experimenter. With such formulation quantum mechanics provides joint probabilities for conjointly observed events, rather than a probability conditional on some initial state for the system under study. This feature is desirable, as with quantum mechanics, statements on what happened in the past may have no meaning in the present.


title: Generalized Persistence Diagrams for Persistence Modules over Posets

abstract: When a category $\mathcal{C}$ satisfies certain conditions, we define the notion of rank invariant for arbitrary poset-indexed functors $F:\mathbf{P} \rightarrow \mathcal{C}$ from a category theory perspective. This generalizes the standard notion of rank invariant as well as Patel's recent extension. Specifically, the barcode of any interval decomposable persistence modules $F:\mathbf{P} \rightarrow \mathbf{vec}$ of finite dimensional vector spaces can be extracted from the rank invariant by the principle of inclusion-exclusion. Generalizing this idea allows freedom of choosing the indexing poset $\mathbf{P}$ of $F: \mathbf{P} \rightarrow \mathcal{C}$ in defining Patel's generalized persistence diagram of $F$. Of particular importance is the fact that the generalized persistence diagram of $F$ is defined regardless of whether $F$ is interval decomposable or not. By specializing our idea to zigzag persistence modules, we also show that the barcode of a Reeb graph can be obtained in a purely set-theoretic setting without passing to the category of vector spaces. This leads to a promotion of Patel's semicontinuity theorem about type $\mathcal{A}$ persistence diagram to Lipschitz continuity theorem for the category of sets.


title: An Acceleration Scheme to The Local Directional Pattern

abstract: This study seeks to improve the running time of the Local Directional Pattern (LDP) during feature extraction using a newly proposed acceleration scheme to LDP. LDP is considered to be computationally expensive. To confirm this, the running time of the LDP to gray level co-occurrence matrix (GLCM) were it was established that the running time for LDP was two orders of magnitude higher than that of the GLCM. In this study, the performance of the newly proposed acceleration scheme was evaluated against LDP and Local Binary patter (LBP) using images from the publicly available extended Cohn-Kanade (CK+) dataset. Based on our findings, the proposed acceleration scheme significantly improves the running time of the LDP by almost 3 times during feature extraction


title: Beyond Complex Langevin Equations: a Progress Report

abstract: After a short review of one of proposals to avoid complex stochastic processes in Complex Langevin studies, the recent progress in the former is reported. In particular, the new developments allow now to construct positive and normalizable representations for gaussian quantum mechanical, as well as field theoretical, path integrals directly in the Minkowski time. A relation to the idea of thimbles is also discussed.


title: Spectrogram-channels u-net: a source separation model viewing each channel as the spectrogram of each source

abstract: Sound source separation has attracted attention from Music Information Retrieval(MIR) researchers, since it is related to many MIR tasks such as automatic lyric transcription, singer identification, and voice conversion. In this paper, we propose an intuitive spectrogram-based model for source separation by adapting U-Net. We call it Spectrogram-Channels U-Net, which means each channel of the output corresponds to the spectrogram of separated source itself. The proposed model can be used for not only singing voice separation but also multi-instrument separation by changing only the number of output channels. In addition, we propose a loss function that balances volumes between different sources. Finally, we yield performance that is state-of-the-art on both separation tasks.


title: Whetstone: A Method for Training Deep Artificial Neural Networks for Binary Communication

abstract: This paper presents a new technique for training networks for low-precision communication. Targeting minimal communication between nodes not only enables the use of emerging spiking neuromorphic platforms, but may additionally streamline processing conventionally. Low-power and embedded neuromorphic processors potentially offer dramatic performance-per-Watt improvements over traditional von Neumann processors, however programming these brain-inspired platforms generally requires platform-specific expertise which limits their applicability. To date, the majority of artificial neural networks have not operated using discrete spike-like communication. We present a method for training deep spiking neural networks using an iterative modification of the backpropagation optimization algorithm. This method, which we call Whetstone, effectively and reliably configures a network for a spiking hardware target with little, if any, loss in performance. Whetstone networks use single time step binary communication and do not require a rate code or other spike-based coding scheme, thus producing networks comparable in timing and size to conventional ANNs, albeit with binarized communication. We demonstrate Whetstone on a number of image classification networks, describing how the sharpening process interacts with different training optimizers and changes the distribution of activity within the network. We further note that Whetstone is compatible with several non-classification neural network applications, such as autoencoders and semantic segmentation. Whetstone is widely extendable and currently implemented using custom activation functions within the Keras wrapper to the popular TensorFlow machine learning framework.


title: Common lines modeling for reference free ab-initio reconstruction in cryo-EM

abstract: We consider the problem of estimating an unbiased and reference-free ab-inito model for non-symmetric molecules from images generated by single-particle cryo-electron microscopy. The proposed algorithm finds the globally optimal assignment of orientations that simultaneously respects all common lines between all images. The contribution of each common line to the estimated orientations is weighted according to a statistical model for common lines' detection errors. The key property of the proposed algorithm is that it finds the global optimum for the orientations given the common lines. In particular, any local optima in the common lines energy landscape do not affect the proposed algorithm. As a result, it is applicable to thousands of images at once, very robust to noise, completely reference free, and not biased towards any initial model. A byproduct of the algorithm is a set of measures that allow to asses the reliability of the obtained ab-initio model. We demonstrate the algorithm using class averages from two experimental data sets, resulting in ab-initio models with resolutions of 20A or better, even from class averages consisting of as few as three raw images per class.


title: Auger-induced charge migration

abstract: Novel perspectives of controlling molecular systems have recently arisen from the possibility of generating attosecond pulses in the ultraviolet regime and tailoring electron dynamics in its natural time scale. The cornerstone mechanism is the so-called charge migration, he production of a coherent charge transfer with subfemtosecond oscillations across a molecule. Typically, charge migration is induced by the ionization of valence molecular orbitals. However, recent technological developments allow the generation of attosecond pulses in the x-ray regime. In this case, the absorption of photons creates core-hole states. In light elements, core-hole states mainly decay by Auger processes that, driven by electron correlations, involve valence orbitals. We theoretically demonstrate in a fluoroacetylene molecule a double-hole charge migration triggered by attosecond core-electron photoionization, followed by Auger electron relaxations. This opens a new route for inducing with x rays charge transfer processes in the subfemtosecond time scale.


title: Euler's Function on Products of Primes in Progressions

abstract: We study generalizations of some results of Jean-Louis Nicolas regarding the relation between small values of Euler's function $\varphi(n)$ and the Riemann Hypothesis. Among other things, we prove that for $1\leq q\leq 10$ and for $q=12, 14$, the generalized Riemann Hypothesis for the Dedekind zeta function of the cyclotomic field $\mathbb{Q}(e^{2\pi i/q})$ is true if and only if for all integers $k\geq 1$ we have \[\frac{\bar{N}_k}{\varphi(\bar{N}_k)(\log(\varphi(q)\log{\bar{N}_k}))^{\frac{1}{\varphi(q)}}} > \frac{1}{C(q,1)}.\] Here $\bar{N}_k$ is the product of the first $k$ primes in the arithmetic progression $p\equiv 1~({\rm mod}~{q})$ and $C(q, 1)$ is the constant appearing in the asymptotic formula \[\prod_{\substack{p \leq x \\ p \equiv 1~({\rm mod}~{q})}} \left(1 - \frac{1}{p}\right) \sim \frac{C(q, 1)}{(\log{x})^\frac{1}{\varphi(q)}},\] as $x\rightarrow\infty$. We also prove that, for $q\leq 400,000$ and integers $a$ coprime to $q$, the analogous inequality \[\frac{\bar{N}_k}{\varphi(\bar{N}_k)(\log(\varphi(q)\log{\bar{N}_k}))^{\frac{1}{\varphi(q)}}} > \frac{1}{C(q,a)}\] holds for infinitely many values of $k$. If in addition $a$ is a not a square modulo $q$, then there are infinitely many $k$ for which this inequality holds and also infinitely many $k$ for which this inequality fails.


title: Semantic Mapping with Simultaneous Object Detection and Localization

abstract: We present a filtering-based method for semantic mapping to simultaneously detect objects and localize their 6 degree-of-freedom pose. For our method, called Contextual Temporal Mapping (or CT-Map), we represent the semantic map as a belief over object classes and poses across an observed scene. Inference for the semantic mapping problem is then modeled in the form of a Conditional Random Field (CRF). CT-Map is a CRF that considers two forms of relationship potentials to account for contextual relations between objects and temporal consistency of object poses, as well as a measurement potential on observations. A particle filtering algorithm is then proposed to perform inference in the CT-Map model. We demonstrate the efficacy of the CT-Map method with a Michigan Progress Fetch robot equipped with a RGB-D sensor. Our results demonstrate that the particle filtering based inference of CT-Map provides improved object detection and pose estimation with respect to baseline methods that treat observations as independent samples of a scene.


title: Algebraic tests of general Gaussian latent tree models

abstract: We consider general Gaussian latent tree models in which the observed variables are not restricted to be leaves of the tree. Extending related recent work, we give a full semi-algebraic description of the set of covariance matrices of any such model. In other words, we find polynomial constraints that characterize when a matrix is the covariance matrix of a distribution in a given latent tree model. However, leveraging these constraints to test a given such model is often complicated by the number of constraints being large and by singularities of individual polynomials, which may invalidate standard approximations to relevant probability distributions. Illustrating with the star tree, we propose a new testing methodology that circumvents singularity issues by trading off some statistical estimation efficiency and handles cases with many constraints through recent advances on Gaussian approximation for maxima of sums of high-dimensional random vectors. Our test avoids the need to maximize the possibly multimodal likelihood function of such models and is applicable to models with larger number of variables. These points are illustrated in numerical experiments.


title: Synthesizing Symmetric Lenses

abstract: Lenses are programs that can be run both "front to back" and "back to front," allowing updates to either their source or their target data to be transferred in both directions. Lenses have been extensively studied, extended, and applied. Recent work has demonstrated how techniques from type-directed program synthesis can be used to efficiently synthesize a simple class of lenses---bijective lenses over string data---given a pair of types (regular expressions) and examples. We extend this synthesis algorithm to a broader class of lenses, called simple symmetric lenses, including all bijective lenses, all of the popular category of "asymmetric" lenses, and a subset of the "symmetric lenses" proposed by Hofmann et al. Intuitively, simple symmetric lenses allow some information to be present on one side but not the other and vice versa. They are of independent theoretical interest, being the largest class of symmetric lenses that do not use persistent internal state. Synthesizing simple symmetric lenses is more challenging than synthesizing bijective lenses: Since some of the information on each side can be "disconnected" from the other side, there will typically be many lenses that agree with a given example. To guide the search process, we use stochastic regular expressions and information theory to estimate the amount of information propagated by a candidate lens, preferring lenses that propagate more information, as well as user annotations marking parts of the source and target formats as either irrelevant or essential. We describe an implementation of simple symmetric lenses and our synthesis procedure as extensions to the Boomerang language. We evaluate its performance on 48 benchmark examples drawn from Flash Fill, Augeas, and the bidirectional programming literature. Our implementation can synthesize each of these lenses in under 30 seconds.


title: Learn to Fly: A Distributed Mechanism for Joint 3D Placement and Users Association in UAVs-assisted Networks

abstract: In this paper, we study the joint 3D placement of unmanned aerial vehicles (UAVs) and UAVs-users association under bandwidth limitation and quality of service constraint. In particular, in order to allow to UAVs to dynamically improve their 3D locations in a distributed fashion while maximizing the network's sum-rate, we break the underlying optimization into 3 subproblems where we separately solve the 2D UAVs positioning, the altitude optimization, and the UAVs-users association. First, given fixed 3D positions of UAVs, we propose a fully distributed matching based association that alleviates the bottlenecks of bandwidth allocation and guarantees the required quality of service. Next, to address the 2D positions of UAVs, we adopt a modified version of K-means algorithm, with a distributed implementation, where UAVs dynamically change their 2D positions in order to reach the barycenter of the cluster that is composed of the served ground users. In order to optimize the UAVs altitudes, we study a naturally defined game-theoretic version of the problem and show that under fixed UAVs 2D coordinates, a predefined association scheme, and limited-interference, the UAVs altitudes game is a non-cooperative potential game where the players (UAVs) can maximize the limited-interference sum-rate by only optimizing a local utility function. Therefore, we adopt the best response dynamics to reach a Nash equilibrium of the game which is also a local optimum of the social welfare function. Our simulation results show that, using the proposed approach, the network's sum rate of the studied scenario is improved as compared with the trivial case where the classical version of K-means is adopted and users are assigned, at each iteration, to the closest UAV.


title: On Some Configurations of Oppositely Charged Trapped Vortices in the Plane

abstract: Our aim in the present work is to identify all the possible standing wave configurations involving few vortices of different charges in an atomic Bose-Einstein condensate (BEC). In this effort, we deploy the use of a computational algebra approach in order to identify stationary multi-vortex states with up to 6 vortices. The use of invariants and symmetries enables deducing a set of equations in elementary symmetric polynomials, which can then be fully solved via computational algebra packages within Maple. We retrieve a number of previously identified configurations, including collinear ones and polygonal (e.g. quadrupolar and hexagonal) ones. However, importantly, we also retrieve a configuration with 4 positive charges and 2 negative ones which is unprecedented, to the best of our knowledge, in BEC studies. We corroborate these predictions via numerical computations in the fully two-dimensional PDE system of the Gross-Pitaevskii type which characterizes the BEC at the mean-field level.


title: Automatic differentiation in ML: Where we are and where we should be going

abstract: We review the current state of automatic differentiation (AD) for array programming in machine learning (ML), including the different approaches such as operator overloading (OO) and source transformation (ST) used for AD, graph-based intermediate representations for programs, and source languages. Based on these insights, we introduce a new graph-based intermediate representation (IR) which specifically aims to efficiently support fully-general AD for array programming. Unlike existing dataflow programming representations in ML frameworks, our IR naturally supports function calls, higher-order functions and recursion, making ML models easier to implement. The ability to represent closures allows us to perform AD using ST without a tape, making the resulting derivative (adjoint) program amenable to ahead-of-time optimization using tools from functional language compilers, and enabling higher-order derivatives. Lastly, we introduce a proof of concept compiler toolchain called Myia which uses a subset of Python as a front end.


title: On invariant $1$-dimensional representations of a finite $W$-algebra

abstract: Let $\mathfrak{g}$ be a simple Lie algebra over $\mathbb{C}$ and $G$ be the corresponding simply connected algebraic group. Consider a nilpotent element $e\in \mathfrak{g}$, the corresponding element $\chi=(e, \bullet)$ in $\mathfrak{g}^*$, and the coadjoint orbit $\mathbb{O}=G\chi$. We are interested in the set $\mathfrak{J}\mathfrak{d}^1(\mathcal{W})$ of codimension $1$ ideals $J\subset \mathcal{W}$ in a finite $W$-algebra $\mathcal{W}=U(\mathfrak{g}, e)$. We have a natural action of the component group $\Gamma=Z_G(\chi)/Z_G^\circ(\chi)$ on $\mathfrak{J}\mathfrak{d}^1(\mathcal{W})$. Denote the set of $\Gamma$-stable points of $\mathfrak{J}\mathfrak{d}^1(\mathcal{W})$ by $\mathfrak{J}\mathfrak{d}^{1}(\mathcal{W})^{\Gamma}$. For a classical $\mathfrak{g}$ Premet and Topley proved that $\mathfrak{J}\mathfrak{d}^{1}(\mathcal{W})^{\Gamma}$ is isomorphic to an affine space. In this paper we will give an easier and shorter proof of this fact.


title: A Theoretical Framework of Approximation Error Analysis of Evolutionary Algorithms

abstract: In the empirical study of evolutionary algorithms, the solution quality is evaluated by either the fitness value or approximation error. The latter measures the fitness difference between an approximation solution and the optimal solution. Since the approximation error analysis is more convenient than the direct estimation of the fitness value, this paper focuses on approximation error analysis. However, it is straightforward to extend all related results from the approximation error to the fitness value. Although the evaluation of solution quality plays an essential role in practice, few rigorous analyses have been conducted on this topic. This paper aims at establishing a novel theoretical framework of approximation error analysis of evolutionary algorithms for discrete optimization. This framework is divided into two parts. The first part is about exact expressions of the approximation error. Two methods, Jordan form and Schur's triangularization, are presented to obtain an exact expression. The second part is about upper bounds on approximation error. Two methods, convergence rate and auxiliary matrix iteration, are proposed to estimate the upper bound. The applicability of this framework is demonstrated through several examples.


title: Third boundary of the Shastry-Sutherland Model by Numerical Diagonalization

abstract: The Shastry-Sutherland model --- the $S=1/2$ Heisenberg antiferromagnet on the square lattice accompanied by orthogonal dimerized interactions --- is studied by the numerical-diagonalization method. Large-scale calculations provide results for larger clusters that have not been reported yet. The present study successfully captures the phase boundary between the dimer and plaquette-singlet phases and clarifies that the spin gap increases once when the interaction forming the square lattice is increased from the boundary. Our calculations strongly suggest that in addition to the edge of the dimer phase given by $J_{2}/J_{1}\sim 0.675$ and the edge of the N$\acute{\rm e}$el-ordered phase given by $J_{2}/J_{1}\sim 0.76$, there exists a third boundary ratio $J_{2}/J_{1}\sim 0.70$ that divides the intermediate region into two parts, where $J_{1}$ and $J_{2}$ denote dimer and square-lattice interactions, respectively.


title: A simple controller for the transition maneuver of a tail-sitter drone

abstract: This paper presents a controller for the transition maneuver of a tail-sitter drone. The tail-sitter model considers aerodynamic terms whereas the proposed controller considers the time-scale separation between drone attitude and position dynamics. The controller design is based on Lyapunov approach and linear saturation functions. Simulations experiments demonstrate the effectiveness of the derived theoretical results.


title: The diverse galaxy counts in the environment of high-redshift massive black holes in Horizon-AGN

abstract: High-redshift quasars are believed to reside in highly biased regions of the Universe, where black hole (BH) growth is sustained by an enhanced number of mergers and by being at the intersection of filaments bringing fresh gas. This assumption should be supported by an enhancement of the number counts of galaxies in the field of view of quasars. While the current observations of quasar environments do not lead to a consensus on a possible excess of galaxies, the future missions JWST, WFIRST, and Euclid will provide new insights on quasar environments, and will substantially increase the number of study-cases. We are in a crucial period, where we need to both understand the current observations and predict how upcoming missions will improve our understanding of BH environments. Using the large-scale simulation Horizon-AGN, we find that statistically the most massive BHs reside in environments with the largest number counts of galaxies. However, we find a large variance in galaxy number counts, and some massive BHs do not show enhanced counts in their neighborhood. Interestingly, some massive BHs have a very close galaxy companion but no further enhancement of the galaxy number counts at larger scales, in agreement with recent observations. We find that AGN feedback in the surrounding galaxies is able to decrease their luminosity and stellar mass, and therefore to make them un-observable when using restrictive galaxy selection criteria. Radiation from the quasars can spread over large distances, which could affect the formation history of surrounding galaxies, but a careful analysis of these processes requires radiative transfer simulations.


title: Automatic Graphics Program Generation using Attention-Based Hierarchical Decoder

abstract: Recent progress on deep learning has made it possible to automatically transform the screenshot of Graphic User Interface (GUI) into code by using the encoder-decoder framework. While the commonly adopted image encoder (e.g., CNN network), might be capable of extracting image features to the desired level, interpreting these abstract image features into hundreds of tokens of code puts a particular challenge on the decoding power of the RNN-based code generator. Considering the code used for describing GUI is usually hierarchically structured, we propose a new attention-based hierarchical code generation model, which can describe GUI images in a finer level of details, while also being able to generate hierarchically structured code in consistency with the hierarchical layout of the graphic elements in the GUI. Our model follows the encoder-decoder framework, all the components of which can be trained jointly in an end-to-end manner. The experimental results show that our method outperforms other current state-of-the-art methods on both a publicly available GUI-code dataset as well as a dataset established by our own.


title: Hyperplane arrangements and mixed Hodge numbers of the Milnor fiber

abstract: For each complex central essential hyperplane arrangement $\mathcal{A}$, let $F_{\mathcal{A}}$ denote its Milnor fiber. We use Tevelev's theory of tropical compactifications to study invariants related to the mixed Hodge structure on the cohomology of $F_{\mathcal{A}}$. We prove that the map taking each arrangement $\mathcal{A}$ to the Hodge-Deligne polynomial of $F_{\mathcal{A}}$ is locally constant on the realization space of any loop-free matroid. When $\mathcal{A}$ consists of distinct hyperplanes, we also give a combinatorial description for the homotopy type of the boundary complex of any simple normal crossing compactification of $F_{\mathcal{A}}$. As a direct consequence, we obtain a combinatorial formula for the top weight cohomology of $F_{\mathcal{A}}$, recovering a result of Dimca and Lehrer.


title: Odderon, HEGS model and LHC data

abstract: We show that the impact of the maximal odderon amplitude at t=0 and centre-of-momentum energy 13 TeV is small. We obtain a value of rho(t=0) at 13 TeV of the order of 0.12. The real part of the odderon amplitude grows like log(s/s0) at high energies, and is calculated from the analytic properties of the amplitude. In the framework of the HEGS model, taking the same intercept for the odderon and the pomeron leads to a good fit of the new LHC data at 13 TeV. We also show that the main effect of the odderon can be seen in the region of the diffraction minimum of the differential elastic cross section. The form and energy dependence of the odderon amplitude determined in the HEGS model reproduce the characteristics of the diffraction minimum at 7, 8 and 13 TeV.


title: Earthquake catalog-based machine learning identification of laboratory fault states and the effects of magnitude of completeness

abstract: Machine learning regression can predict macroscopic fault properties such as shear stress, friction, and time to failure using continuous records of fault zone acoustic emissions. Here we show that a similar approach is successful using event catalogs derived from the continuous data. Our methods are applicable to catalogs of arbitrary scale and magnitude of completeness. We investigate how machine learning regression from an event catalog of laboratory earthquakes performs as a function of the catalog magnitude of completeness. We find that strong model performance requires a sufficiently low magnitude of completeness, and below this magnitude of completeness model performance saturates.


title: Sub-Hermitian Geometry and the Quantitative Newlander-Nirenberg Theorem

abstract: Given a finite collection of $C^1$ complex vector fields on a $C^2$ manifold $M$ such that they and their complex conjugates span the complexified tangent space at every point, the classical Newlander-Nirenberg theorem gives conditions on the vector fields so that there is a complex structure on $M$ with respect to which the vector fields are $T^{0,1}$. In this paper, we give intrinsic, diffeomorphic invariant, necessary and sufficient conditions on the vector fields so that they have a desired level of regularity with respect to this complex structure (i.e., smooth, real analytic, or have Zygmund regularity of some finite order). By addressing this in a quantitative way we obtain a holomorphic analog of the quantitative theory of sub-Riemannian geometry initiated by Nagel, Stein, and Wainger. We call this sub-Hermitian geometry. Moreover, we proceed more generally and obtain similar results for manifolds which have an associated formally integrable elliptic structure. This allows us to introduce a setting which generalizes both the real and complex theories.


title: Human-Robot Trust Integrated Task Allocation and Symbolic Motion planning for Heterogeneous Multi-robot Systems

abstract: This paper presents a human-robot trust integrated task allocation and motion planning framework for multi-robot systems (MRS) in performing a set of tasks concurrently. A set of task specifications in parallel are conjuncted with MRS to synthesize a task allocation automaton. Each transition of the task allocation automaton is associated with the total trust value of human in corresponding robots. Here, the human-robot trust model is constructed with a dynamic Bayesian network (DBN) by considering individual robot performance, safety coefficient, human cognitive workload and overall evaluation of task allocation. Hence, a task allocation path with maximum encoded human-robot trust can be searched based on the current trust value of each robot in the task allocation automaton. Symbolic motion planning (SMP) is implemented for each robot after they obtain the sequence of actions. The task allocation path can be intermittently updated with this DBN based trust model. The overall strategy is demonstrated by a simulation with 5 robots and 3 parallel subtask automata.


title: Revisiting CFR+ and Alternating Updates

abstract: The CFR+ algorithm for solving imperfect information games is a variant of the popular CFR algorithm, with faster empirical performance on a range of problems. It was introduced with a theoretical upper bound on solution error, but subsequent work showed an error in one step of the proof. We provide updated proofs to recover the original bound.


title: The Origin of Double-peaked Narrow Lines in Active Galactic Nuclei. IV. Association with Galaxy Mergers

abstract: Double-peaked narrow emission lines in active galactic nucleus (AGN) spectra can be produced by AGN outflows, rotation, or dual AGNs, which are AGN pairs in ongoing galaxy mergers. Consequently, double-peaked narrow AGN emission lines are useful tracers of the coevolution of galaxies and their supermassive black holes, as driven by AGN feedback and AGN fueling. We investigate this concept further with follow-up optical longslit observations of a sample of 95 Sloan Digital Sky Survey (SDSS) galaxies that have double-peaked narrow AGN emission lines. Based on a kinematic analysis of the longslit spectra, we confirm previous work that finds that the majority of double-peaked narrow AGN emission lines are associated with outflows. We also find that eight of the galaxies have companion galaxies with line-of-sight velocity separations < 500 km/s and physical separations <30 kpc. Since we find evidence of AGNs in both galaxies, all eight of these systems are compelling dual AGN candidates. Galaxies with double-peaked narrow AGN emission lines occur in such galaxy mergers at least twice as often as typical active galaxies. Finally, we conclude that at least 3% of SDSS galaxies with double-peaked narrow AGN emission lines are found in galaxy mergers where both galaxies are resolved in SDSS imaging.


title: Quantifying Learning Guarantees for Convex but Inconsistent Surrogates

abstract: We study consistency properties of machine learning methods based on minimizing convex surrogates. We extend the recent framework of Osokin et al. (2017) for the quantitative analysis of consistency properties to the case of inconsistent surrogates. Our key technical contribution consists in a new lower bound on the calibration function for the quadratic surrogate, which is non-trivial (not always zero) for inconsistent cases. The new bound allows to quantify the level of inconsistency of the setting and shows how learning with inconsistent surrogates can have guarantees on sample complexity and optimization difficulty. We apply our theory to two concrete cases: multi-class classification with the tree-structured loss and ranking with the mean average precision loss. The results show the approximation-computation trade-offs caused by inconsistent surrogates and their potential benefits.


title: Efficiently Combining Human Demonstrations and Interventions for Safe Training of Autonomous Systems in Real-Time

abstract: This paper investigates how to utilize different forms of human interaction to safely train autonomous systems in real-time by learning from both human demonstrations and interventions. We implement two components of the Cycle-of-Learning for Autonomous Systems, which is our framework for combining multiple modalities of human interaction. The current effort employs human demonstrations to teach a desired behavior via imitation learning, then leverages intervention data to correct for undesired behaviors produced by the imitation learner to teach novel tasks to an autonomous agent safely, after only minutes of training. We demonstrate this method in an autonomous perching task using a quadrotor with continuous roll, pitch, yaw, and throttle commands and imagery captured from a downward-facing camera in a high-fidelity simulated environment. Our method improves task completion performance for the same amount of human interaction when compared to learning from demonstrations alone, while also requiring on average 32% less data to achieve that performance. This provides evidence that combining multiple modes of human interaction can increase both the training speed and overall performance of policies for autonomous systems.


title: Mobile Sensor Data Anonymization

abstract: Motion sensors such as accelerometers and gyroscopes measure the instant acceleration and rotation of a device, in three dimensions. Raw data streams from motion sensors embedded in portable and wearable devices may reveal private information about users without their awareness. For example, motion data might disclose the weight or gender of a user, or enable their re-identification. To address this problem, we propose an on-device transformation of sensor data to be shared for specific applications, such as monitoring selected daily activities, without revealing information that enables user identification. We formulate the anonymization problem using an information-theoretic approach and propose a new multi-objective loss function for training deep autoencoders. This loss function helps minimizing user-identity information as well as data distortion to preserve the application-specific utility. The training process regulates the encoder to disregard user-identifiable patterns and tunes the decoder to shape the output independently of users in the training set. The trained autoencoder can be deployed on a mobile or wearable device to anonymize sensor data even for users who are not included in the training dataset. Data from 24 users transformed by the proposed anonymizing autoencoder lead to a promising trade-off between utility and privacy, with an accuracy for activity recognition above 92% and an accuracy for user identification below 7%.


title: Unsupervised Multi-Target Domain Adaptation: An Information Theoretic Approach

abstract: Unsupervised domain adaptation (uDA) models focus on pairwise adaptation settings where there is a single, labeled, source and a single target domain. However, in many real-world settings one seeks to adapt to multiple, but somewhat similar, target domains. Applying pairwise adaptation approaches to this setting may be suboptimal, as they fail to leverage shared information among multiple domains. In this work we propose an information theoretic approach for domain adaptation in the novel context of multiple target domains with unlabeled instances and one source domain with labeled instances. Our model aims to find a shared latent space common to all domains, while simultaneously accounting for the remaining private, domain-specific factors. Disentanglement of shared and private information is accomplished using a unified information-theoretic approach, which also serves to establish a stronger link between the latent representations and the observed data. The resulting model, accompanied by an efficient optimization algorithm, allows simultaneous adaptation from a single source to multiple target domains. We test our approach on three challenging publicly-available datasets, showing that it outperforms several popular domain adaptation methods.


title: On the Identifiability of the Influence Model for Stochastic Spatiotemporal Spread Processes

abstract: The influence model is a discrete-time stochastic model that succinctly captures the interactions of a network of Markov chains. The model produces a reduced-order representation of the stochastic network, and can be used to describe and tractably analyze probabilistic spatiotemporal spread dynamics, and hence has found broad usage in network applications such as social networks, traffic management, and failure cascades in power systems. This paper provides sufficient and necessary conditions for the identifiability of the influence model, and also develops estimators for the model structure through exploiting the model's special properties. In addition, we analyze conditions for the identifiability of the partially observed influence model (POIM), for which not all of the sites can be measured.


title: Birkhoff normal form and long time existence for periodic gravity water waves

abstract: We consider the gravity water waves system with a periodic one-dimensional interface in infinite depth, and prove a rigorous reduction of these equations to Birkhoff normal form up to degree four. This proves a conjecture of Zakharov-Dyachenko [55] suggested by the formal Birkhoff integrability of the water waves Hamiltonian truncated at order four. As a consequence, we also obtain a long-time stability result: periodic perturbations of a flat interface that are of size $\varepsilon$ in a sufficiently smooth Sobolev space lead to solutions that remain regular and small up to times of order $\varepsilon^{-3}$. This time scale is expected to be optimal. Main difficulties in the proof are the presence of non-trivial resonant four-waves interactions, the so-called Benjamin-Feir resonances, the small divisors arising from near-resonances and the quasilinear nature of the equations. Some of the main ingredients that we use are: (1) a reduction procedure to constant coefficient operators up to smoothing remainders, that, together with the verification of key algebraic cancellations of the system, implies the integrability of the equations at non-negative orders; (2) a Poincar\'e-Birkhoff normal form of the smoothing remainders that deals with near-resonances; (3) an a priori algebraic identification argument of the above Poincar\'e- Birkhoff normal form equations with the formal Hamiltonian computations of [55, 19, 27, 17], that allows us to handle the Benjamin-Feir resonances.


title: Deep Convolutional Neural Network Applied to Quality Assessment for Video Tracking

abstract: Surveillance videos often suffer from blur and exposure distortions that occur during acquisition and storage, which can adversely influence following automatic image analysis results on video-analytic tasks. The purpose of this paper is to deploy an algorithm that can automatically assess the presence of exposure distortion in videos. In this work we to design and build one architecture for deep learning applied to recognition of distortions in a video. The goal is to know if the video present exposure distortions. Such an algorithm could be used to enhance or restoration image or to create an object tracker distortion-aware.


title: Estimators for Multivariate Information Measures in General Probability Spaces

abstract: Information theoretic quantities play an important role in various settings in machine learning, including causality testing, structure inference in graphical models, time-series problems, feature selection as well as in providing privacy guarantees. A key quantity of interest is the mutual information and generalizations thereof, including conditional mutual information, multivariate mutual information, total correlation and directed information. While the aforementioned information quantities are well defined in arbitrary probability spaces, existing estimators add or subtract entropies (we term them $\Sigma H$ methods). These methods work only in purely discrete space or purely continuous case since entropy (or differential entropy) is well defined only in that regime. In this paper, we define a general graph divergence measure ($\mathbb{GDM}$), as a measure of incompatibility between the observed distribution and a given graphical model structure. This generalizes the aforementioned information measures and we construct a novel estimator via a coupling trick that directly estimates these multivariate information measures using the Radon-Nikodym derivative. These estimators are proven to be consistent in a general setting which includes several cases where the existing estimators fail, thus providing the only known estimators for the following settings: (1) the data has some discrete and some continuous-valued components (2) some (or all) of the components themselves are discrete-continuous mixtures (3) the data is real-valued but does not have a joint density on the entire space, rather is supported on a low-dimensional manifold. We show that our proposed estimators significantly outperform known estimators on synthetic and real datasets.


title: Motivic zeta functions of hyperplane arrangements

abstract: For each central essential hyperplane arrangement $\mathcal{A}$ over an algebraically closed field, let $Z_\mathcal{A}^{\hat\mu}(T)$ denote the Denef-Loeser motivic zeta function of $\mathcal{A}$. We prove a formula expressing $Z_\mathcal{A}^{\hat\mu}(T)$ in terms of the Milnor fibers of related hyperplane arrangements. We use this formula to show that the map taking each complex arrangement $\mathcal{A}$ to the Hodge-Deligne specialization of $Z_{\mathcal{A}}^{\hat\mu}(T)$ is locally constant on the realization space of any loop-free matroid. We also prove a combinatorial formula expressing the motivic Igusa zeta function of $\mathcal{A}$ in terms of the characteristic polynomials of related arrangements.


title: Measure and Dimension of Sums and Products

abstract: We investigate the Lebesgue measure, Hausdorff dimension, and Fourier dimension of sets of the form $RY + Z, $ where $R \subseteq (0,\infty)$ and $Y, Z \subseteq \mathbb{R}^d$. We prove a theorem on the Lebesgue measure and Hausdorff dimension of $RY+Z$; The theorem is a generalized variant of some theorems of Wolff and Oberlin in which $Y$ is the unit sphere, but its proof is much simpler. We also prove a deeper existence theorem: For each $\alpha \in [0,1]$ and for each non-empty compact set $R \subseteq (0,\infty)$, there exists a compact set $Y \subseteq [1,2]$ such that $\dim_F(Y) = \dim_H(Y) = \overline{\dim_M}(Y) = \alpha$ and $\dim_F(RY) \geq \min\{ 1, \dim_F(R) + \dim_F(Y)\}$. This theorem verifies a weak form of a more general conjecture, and it can be used to produce new Salem sets from old ones.


title: Uniqueness and regularity for the Navier-Stokes-Cahn-Hilliard system

abstract: The motion of two contiguous incompressible and viscous fluids is described within the diffuse interface theory by the so-called Model H. The system consists of the Navier-Stokes equations, which are coupled with the Cahn-Hilliard equation associated to the Ginzburg-Landau free energy with physically relevant logarithmic potential. This model is studied in bounded smooth domain in R^d, d=2 and d=3, and is supplemented with a no-slip condition for the velocity, homogeneous Neumann boundary conditions for the order parameter and the chemical potential, and suitable initial conditions. We study uniqueness and regularity of weak and strong solutions. In a two-dimensional domain, we show the uniqueness of weak solutions and the existence and uniqueness of global strong solutions originating from an initial velocity u_0 in V, namely u_0 in H_0^1 such that div u_0=0. In addition, we prove further regularity properties and the validity of the instantaneous separation property. In a three-dimensional domain, we show the existence and uniqueness of local strong solutions with initial velocity u_0 in V.


title: Coherent systems of probability measures on graphs for representations of free Frobenius towers

abstract: First formally defined by Borodin and Olshanski, a coherent system on a graded graph is a sequence of probability measures which respect the action of certain down/up transition functions between graded components. In one common example of such a construction, each measure is the Plancherel measure for the symmetric group $S_{n}$ and the down transition function is induced from the inclusions $S_{n} \hookrightarrow S_{n+1}$. In this paper we generalize the above framework to the case where $\{A_n\}_{n \geq 0}$ is any free Frobenius tower and $A_n$ is no longer assumed to be semisimple. In particular, we describe two coherent systems on graded graphs defined by the representation theory of $\{A_n\}_{n \geq 0}$ and connect one of these systems to a family of central elements of $\{A_n\}_{n \geq 0}$. When the algebras $\{A_n\}_{n \geq 0}$ are not semisimple, the resulting coherent systems reflect the duality between simple $A_n$-modules and indecomposable projective $A_n$-modules.


title: Efficient and Trustworthy Social Navigation Via Explicit and Implicit Robot-Human Communication

abstract: In this paper, we present a planning framework that uses a combination of implicit (robot motion) and explicit (visual/audio/haptic feedback) communication during mobile robot navigation. First, we developed a model that approximates both continuous movements and discrete behavior modes in human navigation, considering the effects of implicit and explicit communication on human decision making. The model approximates the human as an optimal agent, with a reward function obtained through inverse reinforcement learning. Second, a planner uses this model to generate communicative actions that maximize the robot's transparency and efficiency. We implemented the planner on a mobile robot, using a wearable haptic device for explicit communication. In a user study of an indoor human-robot pair of orthogonal crossing situation, the robot was able to actively communicate its intent to users in order to avoid collisions and facilitate efficient trajectories. Results showed that the planner generated plans that were easier to understand, reduced users' effort, and increased users' trust of the robot, compared to simply performing collision avoidance. The key contribution of this work is the integration and analysis of explicit communication (together with implicit communication) for social navigation.


title: The Duration of Optimal Stopping Problems

abstract: Optimal stopping problems give rise to random distributions describing how many applicants the decision-maker will sample or interview before choosing one, a quantity sometimes referred to as the search time or process duration. This research note surveys several variants of optimal stopping problems, extends earlier results in various directions, and shows how many interviews are expected to be conducted in various settings. The focus is on problems that require a decision-maker to choose a candidate from a pool of sequential applicants with no recall, in the vein of previously studied Cayley-Moser, Secretary and Sultan's Dowry problems.


title: MCA-based Rule Mining Enables Interpretable Inference in Clinical Psychiatry

abstract: Development of interpretable machine learning models for clinical healthcare applications has the potential of changing the way we understand, treat, and ultimately cure, diseases and disorders in many areas of medicine. These models can serve not only as sources of predictions and estimates, but also as discovery tools for clinicians and researchers to reveal new knowledge from the data. High dimensionality of patient information (e.g., phenotype, genotype, and medical history), lack of objective measurements, and the heterogeneity in patient populations often create significant challenges in developing interpretable machine learning models for clinical psychiatry in practice. In this paper we take a step towards the development of such interpretable models. First, by developing a novel categorical rule mining method based on Multivariate Correspondence Analysis (MCA) capable of handling datasets with large numbers of features, and second, by applying this method to build transdiagnostic Bayesian Rule List models to screen for psychiatric disorders using the Consortium for Neuropsychiatric Phenomics dataset. We show that our method is not only at least 100 times faster than state-of-the-art rule mining techniques for datasets with 50 features, but also provides interpretability and comparable prediction accuracy across several benchmark datasets.


title: Blue-enhanced Supercontinuum Generation in a Fluorine-doped Graded-index Multimode Fiber

abstract: Blue-enhanced, broadband supercontinuum generation in a 50 um core fluorine-doped graded-index multimode fiber is demonstrated by pumping with a 1064 nm picosecond source. Multi-octave supercontinuum spectrum extending from ~400 nm to 2,400 nm is achieved in a relatively short fiber. The measured spectrum and multimode nonlinear properties are compared to those of a germanium-doped graded-index multimode fiber. It is shown that engineering the glass composition and thus its chromatic dispersion and absorption characteristics, provides a novel and effective approach to further controlling the nonlinear dynamics in highly multimoded parabolic-index fibers.


title: Topological metal bands with double-triple-point fermions in optical lattices

abstract: Novel fermionic quasiparticles with integer pseudospins in some energy bands, such as pseudospin-1 triple-point fermions, recently attract increasing interest since they are beyond the conventional spin-$1/2$ Dirac and Weyl counterparts. In this paper, we propose a class of pseudospin-1 fermioic excitations emerging in topological metal bands, dubbed double-triple-point (DTP) fermions. We first present a general three-band continuum model with $C_4$ symmetry in three dimensions, which has three types of threefold degenerate points in the bands classified by their topological charges $C=\pm4,\pm2,0$, respectively. They are dubbed DTPs as spin-1 generalization of double-Weyl points. We then construct two-dimensional and three-dimensional tight-binding lattice models of topological metal bands with exotic DTP fermions near the DTPs. In two dimensions, the band gaps close at a trivial DTP with zero Berry phase, which occurs at the transition between the normal and topological insulator phases. In three dimensions, the topological properties of three different DTP fermions in lattice systems are further investigated, and the effects of breaking $C_4$ symmetry are also studied, which generally leads to splitting each quadratic DTP into two linear triple points and gives topological phase diagrams. Using ultracold fermionic atoms in optical lattices, the proposed models can be realized and the topological properties of the DTP fermions can be detected.


title: Mass Measurements of Neutron-Rich Gallium Isotopes Refine Production of Nuclei of the First r-Process Abundance Peak in Neutron Star Merger Calculations

abstract: We report mass measurements of neutron-rich Ga isotopes $^{80-85}$Ga with TRIUMF's Ion Trap for Atomic and Nuclear science (TITAN). The measurements determine the masses of $^{80-83}$Ga in good agreement with previous measurements. The masses of $^{84}$Ga and $^{85}$Ga were measured for the first time. Uncertainties between $25-48$ keV were reached. The new mass values reduce the nuclear uncertainties associated with the production of A $\approx$ 84 isotopes by the \emph{r}-process for astrophysical conditions that might be consistent with a binary neutron star (BNS) merger producing a blue kilonova. Our nucleosynthesis simulations confirm that BNS merger may contribute to the first abundance peak under moderate neutron-rich conditions with electron fractions $Y_e=0.35-0.38$.


title: Monitoring the shape of weather, soundscapes, and dynamical systems: a new statistic for dimension-driven data analysis on large data sets

abstract: Dimensionality-reduction methods are a fundamental tool in the analysis of large data sets. These algorithms work on the assumption that the "intrinsic dimension" of the data is generally much smaller than the ambient dimension in which it is collected. Alongside their usual purpose of mapping data into a smaller dimension with minimal information loss, dimensionality-reduction techniques implicitly or explicitly provide information about the dimension of the data set. In this paper, we propose a new statistic that we call the $\kappa$-profile for analysis of large data sets. The $\kappa$-profile arises from a dimensionality-reduction optimization problem: namely that of finding a projection into $k$-dimensions that optimally preserves the secants between points in the data set. From this optimal projection we extract $\kappa,$ the norm of the shortest projected secant from among the set of all normalized secants. This $\kappa$ can be computed for any $k$; thus the tuple of $\kappa$ values (indexed by dimension) becomes a $\kappa$-profile. Algorithms such as the Secant-Avoidance Projection algorithm and the Hierarchical Secant-Avoidance Projection algorithm, provide a computationally feasible means of estimating the $\kappa$-profile for large data sets, and thus a method of understanding and monitoring their behavior. As we demonstrate in this paper, the $\kappa$-profile serves as a useful statistic in several representative settings: weather data, soundscape data, and dynamical systems data.


title: Three Lectures on Complexity and Black Holes

abstract: Given at PiTP 2018 summer program entitled "From Qubits to Spacetime." The first lecture describes the meaning of quantum complexity, the analogy between entropy and complexity, and the second law of complexity. Lecture two reviews the connection between the second law of complexity and the interior of black holes. I discuss how firewalls are related to periods of non-increasing complexity which typically only occur after an exponentially long time. The final lecture is about the thermodynamics of complexity, and "uncomplexity" as a resource for doing computational work. I explain the remarkable power of "one clean qubit," in both computational terms and in space-time terms. The lectures can also be found online at \url{https://static.ias.edu/pitp/2018/node/1796.html} .


title: New test vector for Waldspurger's period integral, relative trace formula, and hybrid subconvexity bounds

abstract: In this paper we give quantitative local test vectors for Waldspurger's period integral (i.e., a toric period on $\text{GL}_2$) in new cases with joint ramifications. The construction involves minimal vectors, rather than newforms and their variants. This paper gives a uniform treatment for the matrix algebra and division algebra cases under mild assumptions, and establishes an explicit relation between the size of the local integral and the finite conductor $C(\pi\times\pi_{\chi^{-1}})$. As an application, we combine the test vector results with the relative trace formula, and prove a hybrid type subconvexity bound which can be as strong as the Weyl bound in proper range.


title: Radiative-capture cross sections for the $^{139}$La($n,\gamma$) reactionusing thermal neutrons and structural properties of $^{140}$La

abstract: A set of prompt partial $\gamma$-ray production cross sections from thermal neutron-capture were measured for the $^{139}$La($n,\gamma$) reaction using a guided beam of subthermal (thermal and cold) neutrons incident upon a $^{\rm nat}$La$_{2}$O$_{3}$ target at the Prompt Gamma Activation Analysis facility of the Budapest Research Reactor. Absolute $^{140}$La cross sections were determined relative to the well-known comparator $^{35}$Cl($n,\gamma$) cross sections from the irradiation of a stoichiometric $^{\rm nat}$LaCl$_{3}$ sample. The total cross section for radiative thermal neutron-capture on $^{139}$La from the sum of experimentally measured cross sections observed to directly feed the $^{140}$La ground state was determined to be $\sigma_{0}^{\rm expt} = 8.58(50)$ b. To assess completeness of the decay scheme and as a consistency check, the measured cross sections for transitions feeding the ground state from levels below a critical energy of $E_{c} = 285$ keV were combined with a modeled contribution accounting for ground-state feeding from the quasicontinuum to arrive at a total cross section of $\sigma_{0} = 9.36(74)$ b. In addition, a neutron-separation energy of $S_{n} = 5161.005(21)$ keV was determined from a least-squares fit of the measured primary $\gamma$-ray energies to the low-lying levels of the $^{140}$La decay scheme. Furthermore, several nuclear structure improvements are proposed for the decay scheme. The measured cross-section and separation-energy results are comparable to earlier measurements of these quantities.


title: Characterizing vortex tangle properties in steady-state He II counterflow using particle tracking velocimetry

abstract: Historically, there is little faith in particle tracking velocimetry (PTV) as a tool to make quantitative measurements of thermal counterflow in He II, since tracer particle motion is complicated by influences from the normal fluid, superfluid, and quantized vortex lines, or a combination thereof. Recently, we introduced a scheme for differentiating particles trapped on vortices (G1) from particles entrained by the normal fluid (G2). In this paper, we apply this scheme to demonstrate the utility of PTV for quantitative measurements of vortex dynamics in He II counterflow. We estimate $\ell$, the mean vortex line spacing, using G2 velocity data, and $c_2$, a parameter related to the mean curvature radius of vortices and energy dissipation in quantum turbulence, using G1 velocity data. We find that both estimations show good agreement with existing measurements that were obtained using traditional experimental methods. This is of particular consequence since these parameters likely vary in space, and PTV offers the advantage of spatial resolution. We also show a direct link between power-law tails in transverse particle velocity probability density functions (PDFs) and reconnection of vortex lines on which G1 particles are trapped.


title: Crystal orientation-dependent polarization state of high-order harmonics

abstract: We analyze the crystal orientation-dependent polarization state of extreme ultraviolet (XUV) high-order harmonics from bulk magnesium oxide crystals subjected to intense linearly polarized laser fields. We find that only along high-symmetry directions in crystals high-order harmonics follow the polarization direction of the laser field. In general, the polarization direction of high-order harmonics deviates from that of the laser field, and the deviation amplitude depends on the crystal orientation, harmonic order and the strength of the laser field. We use a real-space electron trajectory model to understand the crystal orientation-dependent polarization state of XUV harmonics. The polarization analysis allows us to track the motion of strong-field-driven electron in conduction bands in two dimensions. These results have implications in all-optical probing of atomic-scale structure in real-space, electronic band-structure in momentum space, and in the possibility of generating attosecond pulses with time-dependent polarization in a compact setup.


title: Distributed Stochastic Approximation for Solving Network Optimization Problems Under Random Quantization

abstract: We study distributed optimization problems over a network when the communication between the nodes is constrained, and so information that is exchanged between the nodes must be quantized. This imperfect communication poses a fundamental challenge, and this imperfect communication, if not properly accounted for, prevents the convergence of these algorithms. Our first contribution in this paper is to propose a modified consensus-based gradient method for solving such problems using random (dithered) quantization. This algorithm can be interpreted as a distributed variant of a well-known two-time-scale stochastic algorithm. We then study the convergence and derive upper bounds on the rates of convergence of the proposed method as a function of the bandwidths available between the nodes and the underlying network topology, for both convex and strongly convex objective functions. Our results complement for existing literature where such convergence and explicit formulas of the convergence rates are missing. Finally, we provide numerical simulations to compare the convergence properties of the distributed gradient methods with and without quantization for solving the well-known regression problems over networks, for both quadratic and absolute loss functions.


title: Inner amenable groupoids and central sequences

abstract: We introduce inner amenability for discrete p.m.p. groupoids and investigate its basic properties, examples, and the connection with central sequences in the full group of the groupoid or central sequences in the von Neumann algebra associated with the groupoid. Among other things, we show that every free ergodic p.m.p. compact action of an inner amenable group gives rise to an inner amenable orbit equivalence relation. We also obtain an analogous result for compact extensions of equivalence relations which either are stable or have a non-trivial central sequence in their full group.


title: Bounds on the Dimension of Ext for Finite Groups of Lie Type

abstract: Let $G$ be a finite group of Lie type defined in characteristic $p$, and let $k$ be an algebraically closed field of characteristic $r>0$. We will assume that $r \neq p$ (so, we are in the non-defining characteristic case). Let $V$ be a finite-dimensional irreducible left $kG$-module. In 2011, Guralnick and Tiep found bounds on the dimension of $H^1(G,V)$ in non-defining characteristic, which are independent of $V$. The aim of this paper is to generalize the work of Gurlanick and Tiep. We assume that $G$ is split and use methods of modular Harish-Chandra theory to find bounds on the dimension of $\mathrm{Ext^1}$ between certain irreducible $kG$-modules. We then use Dipper and Du's algorithms to illustrate our bounds in a series of examples.


title: Analysis of KNN Information Estimators for Smooth Distributions

abstract: KSG mutual information estimator, which is based on the distances of each sample to its k-th nearest neighbor, is widely used to estimate mutual information between two continuous random variables. Existing work has analyzed the convergence rate of this estimator for random variables whose densities are bounded away from zero in its support. In practice, however, KSG estimator also performs well for a much broader class of distributions, including not only those with bounded support and densities bounded away from zero, but also those with bounded support but densities approaching zero, and those with unbounded support. In this paper, we analyze the convergence rate of the error of KSG estimator for smooth distributions, whose support of density can be both bounded and unbounded. As KSG mutual information estimator can be viewed as an adaptive recombination of KL entropy estimators, in our analysis, we also provide convergence analysis of KL entropy estimator for a broad class of distributions.

categories: stat.ML cs.IT cs.LG math.IT math.ST stat.TH


title: Decoding Schemes for Foliated Sparse Quantum Error Correcting Codes

abstract: Foliated quantum codes are a resource for fault-tolerant measurement-based quantum error correction for quantum repeaters and for quantum computation. They represent a general approach to integrating a range of possible quantum error correcting codes into larger fault-tolerant networks. Here we present an efficient heuristic decoding scheme for foliated quantum codes, based on message passing between primal and dual code 'sheets'. We test this decoder on two different families of sparse quantum error correcting code: turbo codes and bicycle codes, and show reasonably high numerical performance thresholds. We also present a construction schedule for building such code states.


title: Short-segment heart sound classification using an ensemble of deep convolutional neural networks

abstract: This paper proposes a framework based on deep convolutional neural networks (CNNs) for automatic heart sound classification using short-segments of individual heart beats. We design a 1D-CNN that directly learns features from raw heart-sound signals, and a 2D-CNN that takes inputs of two- dimensional time-frequency feature maps based on Mel-frequency cepstral coefficients (MFCC). We further develop a time-frequency CNN ensemble (TF-ECNN) combining the 1D-CNN and 2D-CNN based on score-level fusion of the class probabilities. On the large PhysioNet CinC challenge 2016 database, the proposed CNN models outperformed traditional classifiers based on support vector machine and hidden Markov models with various hand-crafted time- and frequency-domain features. Best classification scores with 89.22% accuracy and 89.94% sensitivity were achieved by the ECNN, and 91.55% specificity and 88.82% modified accuracy by the 2D-CNN alone on the test set.


title: Stochastic electron heating in the laser and quasi-static electric and magnetic fields

abstract: The dynamics of relativistic electrons in the intense laser radiation and quasi-static electromagnetic fields both along and across to the laser propagating direction are studied in the 3/2 dimensional Hamiltonian framework. It is shown that the unperturbed oscillations of the relativistic electron in these electric fields could exhibit a long tail of harmonics which makes an onset of stochastic electron motion be a primary candidate for electron heating. The Poincar\'e mappings describing the electron motions in the laser and electric fields only are derived from which the criterions for instability are obtained. It follows that for both transverse and longitudinal electric fields, there exist upper limits of the stochastic electron energy depending on the laser intensity and electric field strength. Specifically, these maximum stochastic energies are enhanced by a strong laser intensity but weak electric field. Such stochastic heating would be reduced by the superluminal phase velocity in both cases. The impacts of the magnetic fields on the electron dynamics are different for these two cases and discussed qualitatively. These analytic results are confirmed by the numerical simulations of solving the 3/2D Hamiltonian equations directly.


title: Sampling of Planar Curves: Theory and Fast Algorithms

abstract: We introduce a continuous domain framework for the recovery of a planar curve from a few samples. We model the curve as the zero level set of a trigonometric polynomial. We show that the exponential feature maps of the points on the curve lie on a low-dimensional subspace. We show that the null-space vector of the feature matrix can be used to uniquely identify the curve, given a sufficient number of samples. The worst-case theoretical guarantees show that the number of samples required for unique recovery depends on the bandwidth of the underlying trigonometric polynomial, which is a measure of the complexity of the curve. We introduce an iterative algorithm that relies on the low-rank property of the feature maps to recover the curves when the samples are noisy or when the true bandwidth of the curve is unknown. We also demonstrate the preliminary utility of the proposed curve representation in the context of image segmentation.