title: Product algebras for Galerkin discretisations of boundary integral operators and their applications

abstract: Operator products occur naturally in a range of regularized boundary integral equation formulations. However, while a Galerkin discretisation only depends on the domain space and the test (or dual) space of the operator, products require a notion of the range. In the boundary element software package Bempp we have implemented a complete operator algebra that depends on knowledge of the domain, range and test space. The aim was to develop a way of working with Galerkin operators in boundary element software that is as close to working with the strong form on paper as possible while hiding the complexities of Galerkin discretisations. In this paper, we demonstrate the implementation of this operator algebra and show, using various Laplace and Helmholtz example problems, how it significantly simplifies the definition and solution of a wide range of typical boundary integral equation problems.


title: Constraining Planetary Migration and Tidal Dissipation with Coeval Hot Jupiters

abstract: We investigate the constraints on the formation of, and tidal dissipation processes in, hot Jupiters (HJs) that can be inferred based on reliable knowledge of the age of a system or population. Particular attention is paid to the role of young systems (such as those in open clusters or star-forming regions) in such studies. For an ensemble of coeval HJ (or proto-HJ) systems, we quantify the effect of age on the distribution of orbital eccentricities with respect to orbital periods as well as the location of the observed "pile-up" feature. We expect the effects of pre-main-sequence stellar evolution to be important only if a substantial fraction of HJs approach their current orbits early in protostellar contraction (ages less than 10 Myr). Application to the HJs presently known in the cluster M67 yields constraints on the dissipation roughly consistent with those gleaned from planets in the field; for those in the Hyades and Praesepe, our results suggest a higher degree of dissipation at early times than that inferred from other populations.


title: A recurrent neural network for classification of unevenly sampled variable stars

abstract: Astronomical surveys of celestial sources produce streams of noisy time series measuring flux versus time ("light curves"). Unlike in many other physical domains, however, large (and source-specific) temporal gaps in data arise naturally due to intranight cadence choices as well as diurnal and seasonal constraints. With nightly observations of millions of variable stars and transients from upcoming surveys, efficient and accurate discovery and classification techniques on noisy, irregularly sampled data must be employed with minimal human-in-the-loop involvement. Machine learning for inference tasks on such data traditionally requires the laborious hand-coding of domain-specific numerical summaries of raw data ("features"). Here we present a novel unsupervised autoencoding recurrent neural network (RNN) that makes explicit use of sampling times and known heteroskedastic noise properties. When trained on optical variable star catalogs, this network produces supervised classification models that rival other best-in-class approaches. We find that autoencoded features learned on one time-domain survey perform nearly as well when applied to another survey. These networks can continue to learn from new unlabeled observations and may be used in other unsupervised tasks such as forecasting and anomaly detection.


title: Optimal Dynamic Sensor Subset Selection for Tracking a Time-Varying Stochastic Process

abstract: Motivated by the Internet-of-things and sensor networks for cyberphysical systems, the problem of dynamic sensor activation for the tracking of a time-varying process is examined. The tradeoff is between energy efficiency, which decreases with the number of active sensors, and fidelity, which increases with the number of active sensors. The problem of minimizing the time-averaged mean-squared error over infinite horizon is examined under the constraint of the mean number of active sensors. The proposed methods artfully combine three key ingredients: Gibbs sampling, stochastic approximation for learning, and modifications to consensus algorithms to create a high performance, energy efficient tracking mechanisms with active sensor selection. The following progression of scenarios are considered: centralized tracking of an i.i.d. process; distributed tracking of an i.i.d. process and finally distributed tracking of a Markov chain. The challenge of the i.i.d. case is that the process has a distribution parameterized by a known or unknown parameter which must be learned. The key theoretical results prove that the proposed algorithms converge to local optima for the two i.i.d process cases; numerical results suggest that global optimality is in fact achieved. The proposed distributed tracking algorithm for a Markov chain, based on Kalman-consensus filtering and stochastic approximation, is seen to offer an error performance comparable to that of a competetive centralized Kalman filter.


title: Minimal heteroclinics for a class of fourth order O.D.E. systems

abstract: We prove the existence of minimal heteroclinic orbits for a class of fourth order O.D.E. systems with variational structure. In our general set-up, the set of equilibria of these systems is a union of manifolds, and the heteroclinic orbits connect two disjoint components of this set.


title: Recent Photometric Monitoring of KIC 8462852, the Detection of a Potential Repeat of the Kepler Day 1540 Dip and a Plausible Model

abstract: This paper presents V- and g'-band observations of the F2V star KIC 8462852, which exhibited enigmatic fade patterns in Kepler mission data. We introduce a transit simulation model for interpretation of these fades, and use it to interpret an August 2017 dip as a repeat of the Kepler day 1540 dip (D1540). We suggest the August 2017 and D1540 dips may be caused by a brown dwarf and an associated ring system in a 1601-day elliptical orbit. Transiting icy moons of the proposed brown dwarf, sublimating near periapsis like comets, could provide an explanation for the significant dips observed by Kepler, as well as the recent May to October 2017 dips and the long term variation in flux detected by Simon et al. (2017). Whereas the presence of such a ring structure is attractive for its ability to explain short term fade events, we do not address how such a ring system can be created and maintained. If our speculation is correct, a brightening of about 1-2 percent should occur during October to November 2017. In addition, this scenario predicts that a set of dimming events, similar to those in 2013 (Kepler) and in 2017 (reported here), can be expected to repeat during October 2021 to January 2022 and a repeat of D1540 should occur on 27 December 2021.


title: Algebraic structures related to integrable differential equations

abstract: The survey is devoted to algebraic structures related to integrable ODEs and evolution PDEs. A description of Lax representations is given in terms of vector space decomposition of loop algebras into a direct sum of Taylor series and a complementary subalgebra. Examples of complementary subalgebras and corresponding integrable models are presented. In the framework of the bi-Hamiltonian approach compatible associative algebras related affine Dynkin diagrams are considered. A bi-Hamiltonian origin of the classical elliptic Calogero-Moser models is revealed.


title: Cycle double covers and non-separating cycles

abstract: Which $2$-regular subgraph $R$ of a cubic graph $G$ can be extended to a cycle double cover of $G$? We provide a condition which ensures that every $R$ satisfying this condition is part of a cycle double cover of $G$. As one consequence, we prove that every $2$-connected cubic graph which has a decomposition into a spanning tree and a $2$-regular subgraph $C$ consisting of $k$ circuits with $k\leq 3$, has a cycle double cover containing $C$.


title: ALMA detections of CO emission in the most luminous, heavily dust-obscured quasars at z>3

abstract: We report the results of a pilot study of CO$(4-3)$ emission line of three {\it WISE}-selected hyper-luminous, dust-obscured quasars (QSOs) with sensitive ALMA Band 3 observations. These obscured QSOs with $L_{\rm bol}>10^{14}L_\odot$ are among the most luminous objects in the universe. All three QSO hosts are clearly detected both in continuum and in CO$(4-3)$ emission line. Based on CO$(4-3)$ emission line detection, we derive the molecular gas masses ($\sim 10^{10-11}$ M$_\odot$), suggesting that these QSOs are gas-rich systems. We find that three obscured QSOs in our sample follow the similar $L'_{\rm CO}- L_{\rm FIR}$ relation as unobscured QSOs at high redshifts. We also find the complex velocity structures of CO$(4-3)$ emission line, which provide the possible evidence for gas-rich merger in W0149+2350 and possible molecular outflow in W0220+0137 and W0410$-$0913. Massive molecular outflow can blow away the obscured interstellar medium (ISM) and make obscured QSOs evolve towards the UV/optical bright, unobscured phase. Our result is consistent with the popular AGN feedback scenario involving the co-evolution between the SMBH and host galaxy.


title: Bulk diffusion in a kinetically constrained lattice gas

abstract: In the hydrodynamic regime, the evolution of a stochastic lattice gas with symmetric hopping rules is described by a diffusion equation with density-dependent diffusion coefficient encapsulating all microscopic details of the dynamics. This diffusion coefficient is, in principle, determined by a Green-Kubo formula. In practice, even when the equilibrium properties of a lattice gas are analytically known, the diffusion coefficient cannot be computed except when a lattice gas additionally satisfies the gradient condition. We develop a procedure to systematically obtain analytical approximations for the diffusion coefficient for non-gradient lattice gases with known equilibrium. The method relies on a variational formula found by Varadhan and Spohn which is a version of the Green-Kubo formula particularly suitable for diffusive lattice gases. Restricting the variational formula to finite-dimensional sub-spaces allows one to perform the minimization and gives upper bounds for the diffusion coefficient. We apply this approach to a kinetically constrained non-gradient lattice gas, viz. to the Kob-Andersen model on the square lattice.

categories: cond-mat.stat-mech math-ph math.MP math.PR nlin.CG


title: Towards a geometric variational discretization of compressible fluids: the rotating shallow water equations

abstract: This paper presents a geometric variational discretization of compressible fluid dynamics. The numerical scheme is obtained by discretizing, in a structure preserving way, the Lie group formulation of fluid dynamics on diffeomorphism groups and the associated variational principles. Our framework applies to irregular mesh discretizations in 2D and 3D. It systematically extends work previously made for incompressible fluids to the compressible case. We consider in detail the numerical scheme on 2D irregular simplicial meshes and evaluate the scheme numerically for the rotating shallow water equations. In particular, we investigate whether the scheme conserves stationary solutions, represents well the nonlinear dynamics, and approximates well the frequency relations of the continuous equations, while preserving conservation laws such as mass and total energy.


title: Canonical and Symplectic analysis of actions describing linearized gravity

abstract: By using the canonical and symplectic approaches an (nonstandard) alternative action describing linearized gravity is studied. We identify the complete set of Dirac's constraints, the counting of physical degrees of freedom is performed and the Dirac brackets are constructed. Furthermore, the symplectic analysis is developed which includes the complete set of Faddeev-Jackiw constraints and a symplectic tensor; from that symplectic matrix we show that the generalized Faddeev-Jackiw brackets and the Dirac ones coincide to each other. With all these results at hand, we prove that the number of physical degrees of freedom are eight, thus, we conclude that the theory does not describe the dynamics of linearized gravity. In addition, we also develop the symplectic analysis of standard linearized gravity and we compare the results for both standard and nonstandard theories. Finally we present some remarks and conclusions.


title: G2-structures for N=1 supersymmetric AdS4 solutions of M-theory

abstract: We study the N=1 supersymmetric solutions of D=11 supergravity obtained as a warped product of four-dimensional anti-de-Sitter space with a seven-dimensional Riemannian manifold M. Using the octonion bundle structure on M we reformulate the Killing spinor equations in terms of sections of the octonion bundle on M. The solutions then define a single complexified G2-structure on M or equivalently two real G2-structures. We then study the torsion of these G2-structures and the relationships between them.


title: Tidal Stripping and Post-Merger Relaxation of Dark Matter Halos: Causes and Consequences of Mass Loss

abstract: We study the properties of distinct dark matter halos (i.e., those that are not subhalos) that have a final virial mass $M_{\mathrm{vir}}$ at $z = 0$ less than their peak mass ($M_{\mathrm{peak}}$) in the Bolshoi-Planck cosmological simulation. We identify two primary causes of halo mass loss: relaxation after a major merger and tidal stripping by a massive neighbouring halo. Major mergers initially boost $M_{\mathrm{vir}}$ and typically cause the final halo to become more prolate and less relaxed and to have higher spin and lower NFW concentration. As the halo relaxes, high energy material from the recent merger gradually escapes beyond the virial radius, temporarily resulting in a net negative accretion rate that reduces the halo mass by $5-15\%$ on average. Halos that experience a major merger around $z = 0.4$ typically reach a minimum mass near $z = 0$. Tidal stripping mainly occurs in dense regions, and it causes halos to become less prolate and have lower spins and higher NFW concentrations. Tidally stripped halos often lose a large fraction of their peak mass ($> 20\%$) and most never recover (or even reattain a positive accretion rate). Low mass halos can be strongly affected by both post-merger mass loss and tidal stripping, while high mass halos are predominantly influenced by post-merger mass loss and show few signs of significant tidal stripping.


title: The Large Synoptic Survey Telescope as a Near-Earth Object Discovery Machine

abstract: Using the most recent prototypes, design, and as-built system information, we test and quantify the capability of the Large Synoptic Survey Telescope (LSST) to discover Potentially Hazardous Asteroids (PHAs) and Near-Earth Objects (NEOs). We empirically estimate an expected upper limit to the false detection rate in LSST image differencing, using measurements on DECam data and prototype LSST software and find it to be about $450$~deg$^{-2}$. We show that this rate is already tractable with current prototype of the LSST Moving Object Processing System (MOPS) by processing a 30-day simulation consistent with measured false detection rates. We proceed to evaluate the performance of the LSST baseline survey strategy for PHAs and NEOs using a high-fidelity simulated survey pointing history. We find that LSST alone, using its baseline survey strategy, will detect $66\%$ of the PHA and $61\%$ of the NEO population objects brighter than $H=22$, with the uncertainty in the estimate of $\pm5$ percentage points. By generating and examining variations on the baseline survey strategy, we show it is possible to further improve the discovery yields. In particular, we find that extending the LSST survey by two additional years and doubling the MOPS search window increases the completeness for PHAs to $86\%$ (including those discovered by contemporaneous surveys) without jeopardizing other LSST science goals ($77\%$ for NEOs). This equates to reducing the undiscovered population of PHAs by additional $26\%$ ($15\%$ for NEOs), relative to the baseline survey.


title: An Approach for Accelerating Incompressible Turbulent Flow Simulations Based on Simultaneous Modelling of Multiple Ensembles

abstract: The present paper deals with the problem of improving the efficiency of large scale turbulent flow simulations. The high-fidelity methods for modelling turbulent flows become available for a wider range of applications thanks to the constant growth of the supercomputers performance, however, they are still unattainable for lots of real-life problems. The key shortcoming of these methods is related to the need of simulating a long time integration interval to collect reliable statistics, while the time integration process is inherently sequential. The novel approach with modelling of multiple flow states is discussed in the paper. The suggested numerical procedure allows to parallelize the integration in time by the cost of additional computations. Multiple realizations of the same turbulent flow are performed simultaneously. This allows to use more efficient implementations of numerical methods for solving systems of linear algebraic equations with multiple right-hand sides, operating with blocks of vectors. The simple theoretical estimate for the expected simulation speedup, accounting the penalty of additional computations and the linear solver performance improvement, is presented. The two problems of modelling turbulent flows in a plain channel and in a channel with a matrix of wall-mounted cubes are used to demonstrate the correctness of the proposed estimates and efficiency of the suggested approach as a whole. The simulation speedup by a factor of 2 is shown.


title: Molecular $\Omega_b$ states

abstract: Motivated by the recent finding of five $\Omega_c$ states by the LHCb collaboration, and the successful reproduction of three of them in a recent approach searching for molecular states of meson-baryon with the quantum numbers of $\Omega_c$, we extend these ideas and make predictions for the interaction of meson-baryon in the beauty sector, searching for poles in the scattering matrix that correspond to physical states. We find several $\Omega_b$ states: two states with masses 6405~MeV and 6465~MeV for $J^P= \frac{1}{2}^-$; two more states with masses 6427~MeV and 6665~MeV for $\frac{3}{2}^-$; and three states between 6500 and 6820~MeV, degenerate with $J^P=\frac{1}{2}^-,\,\frac{3}{2}^-$, stemming from the interaction of vector-baryon in the beauty sector.


title: Symmetry approach to integrability and non-associative algebraic structures

abstract: The first part of the book is devoted to the symmetry approach to classification of scalar integrable evolution PDEs with two independent variables. In the second part systems of evolution equations with polynomial homogeneous right-hand side are considered. Such systems are associated with associative algebras and, more generally, with non-associative algebraic structures, whose structural constants are defined by the coefficients of the system under consideration. We explore which structures correspond to polynomial integrable models. Integrable systems polynomial in derivatives are described in geometric terms. The third part we consider integrable vector izotropic and anizotropic systems.


title: Using the Markov Chain Monte Carlo method to study the physical properties GeV-TeV BL Lac objects

abstract: We fit the spectral energy distributions (SEDs) of 46 GeV - TeV BL Lac objects in the frame of leptonic one-zone synchrotron self-Compton (SSC) model and investigate the physical properties of these objects. We use the Markov Chain Monte Carlo (MCMC) method to obtain the basic parameters, such as magnetic field (B), the break energy of the relativistic electron distribution ($\gamma'_{\rm{b}}$) and the electron energy spectral index. Based on the modeling results, we support the following scenarios on GeV-TeV BL Lac objects: (1) Some sources have large Doppler factors, implying other radiation mechanism should be considered. (2) Comparing with FSRQs, GeV-TeV BL Lac objects have weaker magnetic field and larger Doppler factor, which cause the ineffective cooling and shift the SEDs to higher bands. Their jet powers are around $4.0\times 10^{45}~\rm{ erg\cdot s}^{-1}$, comparing with radiation power, $5.0\times 10^{42}~\rm{ erg\cdot s}^{-1}$, indicating that only a small fraction of jet power is transformed into the emission power. (3) For some BL Lacs with large Doppler factors, their jet components could have two substructures, e.g., the fast core and the slow sheath. For most GeV-TeV BL Lacs, Kelvin-Helmholtz instabilities are suppressed by their higher magnetic fields, leading few micro-variability or intro-day variability in the optical bands. (4) Combined with a sample of FSRQs, an anti-correlation between the peak luminosity $L_{\rm {pk}}$ and the peak frequency $\nu_{\rm {pk}}$ is obtained, favoring the blazar sequence scenario. In addition, an anti-correlation between the jet power $P_{\rm {jet}}$ and the break Lorentz factor $\gamma_{\rm {b}}$ also supports the blazar sequence.


title: Multiquark States

abstract: Why do we see certain types of strongly interacting elementary particles and not others? This question was posed over 50 years ago in the context of the quark model. M. Gell-Mann and G. Zweig proposed that the known mesons were $q \bar q$ and baryons $qqq$, with quarks known at the time $u$ ("up"), $d$ ("down"), and $s$ ("strange") having charges (2/3,-1/3,-1/3). Mesons and baryons would then have integral charges. Mesons such as $qq \bar q \bar q$ and baryons such as $qqqq \bar q$ would also have integral charges. Why weren't they seen? They have now been seen, but only with additional heavy quarks and under conditions which tell us a lot about the strong interactions and how they manifest themselves. The present article describes recent progress in our understanding of such "exotic" mesons and baryons.


title: Convergence of an Explicit Iterative Leap-frog Discontinuous Galerkin Method for Time-domain Maxwell's Equations in Anisotropic Materials

abstract: We propose an explicit iterative leap-frog discontinuous Galerkin method for time-domain Maxwell's equations in anisotropic materials and derive its convergence properties. The a priori error estimates are illustrated by numerical means in some experiments. Motivated by a real application which encompasses modeling electromagnetic wave's propagation through the eye's structures, we simulate our model in a 2D domain aiming to represent a simple example of light scattering in the outer nuclear layer of the retina.


title: Coexistence of superconductivity and charge-density wave in the quasi-one-dimensional material HfTe3

abstract: We present the first experimental evidence for metallicity, superconductivity (SC) and the co-existence of charge density waves (CDW) in the quasi-one-dimensional material HfTe3. The existence of such phenomena is a typical characteristic of the transition metal chalcogenides however, without the application of hydrostatic pressure/chemical doping, it is rare for a material to exhibit the co-existence of both states. Materials such as HfTe3 can therefore provide us with a unique insight into the relationship between these multiple ordered states. By improving on the original synthesis conditions, we have successfully synthesised single phase HfTe3 and confirmed the resultant structure by performing Rietveld refinement. Using low temperature resistivity measurements, we provide the first experimental evidence of SC at ~1.4 K as well as a resistive anomaly indicative of a CDW formation at ~82 K. By the application of hydrostatic-pressure, the resistivity anomaly shifts to higher temperature. The results show that HfTe3 is a promising new material to help study the relationship between SC and CDW.


title: Univalent Wandering Domains in the Eremenko-Lyubich Class

abstract: We use the folding theorem of Bishop to construct an entire function $f$ in class $B$ and a wandering domain $U$ of $f$ such that $f$ restricted to $f^n(U)$ is univalent, for all $n \geq 0$. The components of the wandering orbit are bounded and surrounded by the postcritical set.


title: ALMA observations of molecular clouds in three group centered elliptical galaxies: NGC 5846, NGC 4636, and NGC 5044

abstract: We present new ALMA CO(2--1) observations of two well studied group-centered elliptical galaxies: NGC~4636 and NGC~5846. In addition, we include a revised analysis of Cycle 0 ALMA observations of the central galaxy in the NGC~5044 group that has been previously published. We find evidence that molecular gas, in the form of off-center orbiting clouds, is a common presence in bright group-centered galaxies (BGG). CO line widths are $\gtrsim 10$ times broader than Galactic molecular clouds, and using the reference Milky Way $X_{CO}$, the total molecular mass ranges from as low as $2.6\times 10^5 M_\odot$ in NGC~4636 to $6.1\times 10^7 M_\odot$ in NGC~5044. With these parameters the virial parameters of the molecular structures is $\gg 1$. Complementary observations of NGC~5846 and NGC~4636 using the ALMA Compact Array (ACA) do not exhibit any detection of a CO diffuse component at the sensitivity level achieved by current exposures. The origin of the detected molecular features is still uncertain, but these ALMA observations suggest that they are the end product of the hot gas cooling process and not the result of merger events. Some of the molecular clouds are associated with dust features as revealed by HST dust extinction maps suggesting that these clouds formed from dust-enhanced cooling. The global nonlinear condensation may be triggered via the chaotic turbulent field or buoyant uplift. The large virial parameter of the molecular structures and correlation with the warm ($10^3 - 10^5 K$)/hot ($\ge10^6$) phase velocity dispersion provide evidence that they are unbound giant molecular associations drifting in the turbulent field, consistently with numerical predictions of the chaotic cold accretion process. Alternatively, the observed large CO line widths may be generated by molecular gas flowing out from cloud surfaces due to heating by the local hot gas atmosphere.


title: A methodology to evaluate continuum-scale yield surfaces based on the spatial distributions of yielding at the crystal scale

abstract: Using a correlation between local yielding and a multiaxial strength-to-stiffness parameter, the continuum-scale yield surface for a polyphase, polycrystalline solid is predicted. The predicted surface explicitly accounts for microstructure through the quantification of strength-to-stiffness based on a finite element model of a crystal-scale sample. The multiaxial strength-to-stiffness is evaluated from the elastic response of the sample and the restricted slip, single-crystal yield surface. Macroscopic yielding is defined by the propogation of a yield band through the sample and is detected with the aid of a flood fill agorithm. The methodology is demonstrated with the evaluation of a plane-stress yield surface for a dual-phase superaustenitic stainless steel.


title: Comment on Gold-Nanobipyramid Enhanced Hydrogen Sensing with Plasmon Red Shifts Reaching 140 nm at 2 vol% Hydrogen Concentration

abstract: In the recently-published paper "Gold Nanobipyramid-Enhanced Hydrogen Sensing with Plasmon Red Shifts Reaching 140 nm at 2 vol% Hydrogen Concentration" by Yip et al.,[1] the authors claimed that the tuning of catalytic activities of the synthesized Au nanobipyramids(NBP)@Pd nanostructures was mainly due to the electronic effect. Here, we show that this claim is unsupported and an important reference was mis-cited in this commented publication.


title: Nested distance for stagewise-independent processes

abstract: We prove that, for two discrete-time stagewise-independent processes with a stagewise metric, the nested distance is equal to the sum of the Wasserstein distances between the marginal distributions of each stage.


title: Active Betweenness Cardinality: Algorithms and Applications

abstract: Centrality rankings such as degree, closeness, betweenness, Katz, PageRank, etc. are commonly used to identify critical nodes in a graph. These methods are based on two assumptions that restrict their wider applicability. First, they assume the exact topology of the network is available. Secondly, they do not take into account the activity over the network and only rely on its topology. However, in many applications, the network is autonomous, vast, and distributed, and it is hard to collect the exact topology. At the same time, the underlying pairwise activity between node pairs is not uniform and node criticality strongly depends on the activity on the underlying network. In this paper, we propose active betweenness cardinality, as a new measure, where the node criticalities are based on not the static structure, but the activity of the network. We show how this metric can be computed efficiently by using only local information for a given node and how we can find the most critical nodes starting from only a few nodes. We also show how this metric can be used to monitor a network and identify failed nodes.We present experimental results to show effectiveness by demonstrating how the failed nodes can be identified by measuring active betweenness cardinality of a few nodes in the system.


title: Valid Inference Corrected for Outlier Removal

abstract: Ordinary least square (OLS) estimation of a linear regression model is well-known to be highly sensitive to outliers. It is common practice to (1) identify and remove outliers by looking at the data and (2) to fit OLS and form confidence intervals and p-values on the remaining data as if this were the original data collected. This standard "detect-and-forget" approach has been shown to be problematic, and in this paper we highlight the fact that it can lead to invalid inference and show how recently developed tools in selective inference can be used to properly account for outlier detection and removal. Our inferential procedures apply to a general class of outlier removal procedures that includes several of the most commonly used approaches. We conduct simulations to corroborate the theoretical results, and we apply our method to three real data sets to illustrate how our inferential results can differ from the traditional detect-and-forget strategy. A companion R package, outference, implements these new procedures with an interface that matches the functions commonly used for inference with lm in R.


title: CTL* synthesis via LTL synthesis

abstract: We reduce synthesis for CTL* properties to synthesis for LTL. In the context of model checking this is impossible - CTL* is more expressive than LTL. Yet, in synthesis we have knowledge of the system structure and we can add new outputs. These outputs can be used to encode witnesses of the satisfaction of CTL* subformulas directly into the system. This way, we construct an LTL formula, over old and new outputs and original inputs, which is realisable if, and only if, the original CTL* formula is realisable. The CTL*-via-LTL synthesis approach preserves the problem complexity, although it might increase the minimal system size. We implemented the reduction, and evaluated the CTL*-via-LTL synthesiser on several examples.


title: Symbolic vs. Bounded Synthesis for Petri Games

abstract: Petri games are a multiplayer game model for the automatic synthesis of distributed systems. We compare two fundamentally different approaches for solving Petri games. The symbolic approach decides the existence of a winning strategy via a reduction to a two-player game over a finite graph, which in turn is solved by a fixed point iteration based on binary decision diagrams (BDDs). The bounded synthesis approach encodes the existence of a winning strategy, up to a given bound on the size of the strategy, as a quantified Boolean formula (QBF). In this paper, we report on initial experience with a prototype implementation of the bounded synthesis approach. We compare bounded synthesis to the existing implementation of the symbolic approach in the synthesis tool ADAM. We present experimental results on a collection of benchmarks, including one new benchmark family, modeling manufacturing and workflow scenarios with multiple concurrent processes.


title: Asymptotic Expansions of Fundamental Solutions in Parabolic Homogenization

abstract: For a family of second-order parabolic systems with rapidly oscillating and time-dependent periodic coefficients, we investigate the asymptotic behavior of fundamental solutions and establish sharp estimates for the remainders.


title: A Class of Control Certificates to Ensure Reach-While-Stay for Switched Systems

abstract: In this article, we consider the problem of synthesizing switching controllers for temporal properties through the composition of simple primitive reach-while-stay (RWS) properties. Reach-while-stay properties specify that the system states starting from an initial set I, must reach a goal (target) set G in finite time, while remaining inside a safe set S. Our approach synthesizes switched controllers that select between finitely many modes to satisfy the given RWS specification. To do so, we consider control certificates, which are Lyapunov-like functions that represent control strategies to achieve the desired specification. However, for RWS problems, a control Lyapunov-like function is often hard to synthesize in a simple polynomial form. Therefore, we combine control barrier and Lyapunov functions with an additional compatibility condition between them. Using this approach, the controller synthesis problem reduces to one of solving quantified nonlinear constrained problems that are handled using a combination of SMT solvers. The synthesis of controllers is demonstrated through a set of interesting numerical examples drawn from the related work, and compared with the state-of-the-art tool SCOTS. Our evaluation suggests that our approach is computationally feasible, and adds to the growing body of formal approaches to controller synthesis.


title: Notes on Fano Ratio and Portfolio Optimization

abstract: We discuss - in what is intended to be a pedagogical fashion - generalized "mean-to-risk" ratios for portfolio optimization. The Sharpe ratio is only one example of such generalized "mean-to-risk" ratios. Another example is what we term the Fano ratio (which, unlike the Sharpe ratio, is independent of the time horizon). Thus, for long-only portfolios optimizing the Fano ratio generally results in a more diversified and less skewed portfolio (compared with optimizing the Sharpe ratio). We give an explicit algorithm for such optimization. We also discuss (Fano-ratio-inspired) long-short strategies that outperform those based on optimizing the Sharpe ratio in our backtests.


title: SyGuS Techniques in the Core of an SMT Solver

abstract: We give an overview of recent techniques for implementing syntax-guided synthesis (SyGuS) algorithms in the core of Satisfiability Modulo Theories (SMT) solvers. We define several classes of synthesis conjectures and corresponding techniques that can be used when dealing with each class of conjecture.


title: Limit theorems for functionals of two independent Gaussian processes

abstract: Under certain mild conditions, some limit theorems for functionals of two independent Gaussian processes are obtained. The results apply to general Gaussian processes including fractional Brownian motion, sub-fractional Brownian motion and bi-fractional Brownian motion. A new and interesting phenomenon is that, in comparison with the results for fractional Brownian motion, extra randomness appears in the limiting distributions for Gaussian processes with nonstationary increments, say sub-fractional Brownian motion and bi-fractional Brownian. The results are obtained based on the method of moments, in which Fourier analysis, the chaining argument introduced in \cite{nx1} and a paring technique are employed.


title: A Review on Cooperative Diversity Techniques Bypassing Channel Estimation

abstract: Wireless communication technology has seen a remarkably fast evolution due to its capability to provide a quality, reliable and high-speed data transmission amongst the users. However, transmission of information in wireless channels is primarily impaired by deleterious multipath fading, which affects the quality and reliability of the system. In order to overcome the detrimental effects of fading, Multiple-Input Multiple-Output (MIMO) technology is an attractive scheme that employs multiple transceiver antennas to carry the data over the same frequency band over a variety of signal paths. This technology has shown great solutions due to its ability to provide better spectral efficiency, capacity, throughput and robustness of the data transmission. But in practice, it is impractical to install multiple antennas on small-sized devices. Hence, to overcome the limitations of MIMO gain in the future wireless networks, cooperative diversity has recently draw in attention due to its ability to circumvent the difficulties of implementing actual antenna arrays in Multiple-Input and Multiple-Output (MIMO). By exploiting the broadcast feature of the wireless medium, cooperation among multiple nearby nodes is formed for data transmission. At the receiver, the signals are either coherently or differentially detected. Coherent detection requires exact channel estimation, which is difficult to apply in a time-varying channel. Hence, when the nodes are mobile, or when the channel is inaccurately estimated, the differential detection techniques that omit channel estimation become an alternative as compared to coherent detection. This article presents a review of the differential transmission techniques for cooperative diversity networks.


title: PSIque: Next Sequence Prediction of Satellite Images using a Convolutional Sequence-to-Sequence Network

abstract: Predicting unseen weather phenomena is an important issue for disaster management. In this paper, we suggest a model for a convolutional sequence-to-sequence autoencoder for predicting undiscovered weather situations from previous satellite images. We also propose a symmetric skip connection between encoder and decoder modules to produce more comprehensive image predictions. To examine our model performance, we conducted experiments for each suggested model to predict future satellite images from historical satellite images. A specific combination of skip connection and sequence-to-sequence autoencoder was able to generate closest prediction from the ground truth image.


title: Fractional approaches for the distribution of innovation sequence of INAR(1) processes

abstract: In this paper, we present a fractional decomposition of the probability generating function of the innovation process of the first-order non-negative integer-valued autoregressive [INAR(1)] process to obtain the corresponding probability mass function. We also provide a comprehensive review of integer-valued time series models, based on the concept of thinning operators, with geometric-type marginals. In particular, we develop four fractional approaches to obtain the distribution of innovation processes of the INAR(1) model and show that the distribution of the innovations sequence has geometric-type distribution. These approaches are discussed in detail and illustrated through a few examples. Finally, using the methods presented here, we develop four new first-order non-negative integer-valued autoregressive process for autocorrelated counts with overdispersion with known marginals, and derive some properties of these models.


title: Intrinsic Analysis of the Sample Fr\'echet Mean and Sample Mean of Complex Wishart Matrices

abstract: We consider two types of averaging of complex covariance matrices, a sample mean (average) and the sample Fr\'echet mean. We analyse the performance of these quantities as estimators for the true covariance matrix via `intrinsic' versions of bias and mean square error, a methodology which takes account of geometric structure. We derive simple expressions for the intrinsic bias in both cases, and the simple average is seen to be preferable. The same is true for the asymptotic Riemannian risk, and for the Riemannian risk itself in the scalar case. Combined with a similar preference for the simple average using non-intrinsic analysis, we conclude that the simple average is preferred overall to the sample Fr\'echet mean in this context.


title: Split-Decomposition Trees with Prime Nodes: Enumeration and Random Generation of Cactus Graphs

abstract: In this paper, we build on recent results by Chauve et al. (2014) and Bahrani and Lumbroso (2017), which combined the split-decomposition, as exposed by Gioan and Paul, with analytic combinatorics, to produce new enumerative results on graphs---in particular the enumeration of several subclasses of perfect graphs (distance-hereditary, 3-leaf power, ptolemaic). Our goal was to study a simple family of graphs, of which the split-decomposition trees have prime nodes drawn from an enumerable (and manageable!) set of graphs. Cactus graphs, which we describe in more detail further down in this paper, can be thought of as trees with their edges replaced by cycles (of arbitrary lengths). Their split-decomposition trees contain prime nodes that are cycles, making them ideal to study. We derive a characterization for the split-decomposition trees of cactus graphs, produce a general template of symbolic grammars for cactus graphs, and implement random generation for these graphs, building on work by Iriza (2015).


title: Observations of a Field-Aligned Ion/Ion-Beam Instability in a Magnetized Laboratory Plasma

abstract: Collisionless coupling between super Alfv\'{e}nic ions and an ambient plasma parallel to a background magnetic field is mediated by a set of electromagnetic ion/ion-beam instabilities including the resonant right hand instability (RHI). To study this coupling and its role in parallel shock formation, a new experimental configuration at the University of California, Los Angeles utilizes high-energy and high-repetition-rate lasers to create a super-Alfv\'{e}nic field-aligned debris plasma within an ambient plasma in the Large Plasma Device (LAPD). We used a time-resolved fluorescence monochromator and an array of Langmuir probes to characterize the laser plasma velocity distribution and density. The debris ions were observed to be sufficiently super-Alfv\'{e}nic and dense to excite the RHI. Measurements with magnetic flux probes exhibited a right-hand circularly polarized frequency chirp consistent with the excitation of the RHI near the laser target. We compared measurements to 2D hybrid simulations of the experiment.


title: Limit theorems with rate of convergence under sublinear expectations

abstract: Under the sublinear expectation $\mathbb{E}[\cdot]:=\sup_{\theta\in \Theta} E_\theta[\cdot]$ for a given set of linear expectations $\{E_\theta: \theta\in \Theta\}$, we establish a new law of large numbers and a new central limit theorem with rate of convergence. We present some interesting special cases and discuss a related statistical inference problem. We also give an approximation and a representation of the $G$-normal distribution, which was used as the limit in Peng (2007)'s central limit theorem, in a probability space.


title: Algebraic surfaces with $p_g=q=1, K^2=4$ and nonhyperelliptic Albanese fibrations of genus 4

abstract: In this paper we study minimal algebraic surfaces with $p_g=q=1,K^2=4$ and nonhyperelliptic Albanese fibrations of genus 4. We construct for the first time a family of such surfaces as complete intersections of type $(2,3)$ in a $\mathbb{P}^3$-bundle over an elliptic curve. For the surfaces we construct here, the direct image of the canonical sheaf under the Albanese map is decomposable (which is a topological invariant property). Moreover we prove that, all minimal surfaces with $p_g=q=1,K^2=4$ and nonhyperelliptic Albanese fibrations of genus 4 such that the direct image of the canonical sheaf under the Albanese map is decomposable are contained in our family. As a consequence, we show that these surfaces constitute a 4-dimensional irreducible subset $\mathcal{M}$ of $\mathcal{M}_{1,1}^{4,4}$, the Gieseker moduli space of minimal surfaces with $p_g=q=1, K^2=g=4$. Moreover, the closure of $\mathcal{M}$ is an irreducible component of $\mathcal{M}_{1,1}^{4,4}$.


title: Cybersecurity Attacks in Vehicle-to-Infrastructure (V2I) Applications and their Prevention

abstract: A connected vehicle (CV) environment is composed of a diverse data collection, data communication and dissemination, and computing infrastructure systems that are vulnerable to the same cyberattacks as all traditional computing environments. Cyberattacks can jeopardize the expected safety, mobility, energy, and environmental benefits from connected vehicle applications. As cyberattacks can lead to severe traffic incidents, it has become one of the primary concerns in connected vehicle applications. In this paper, we investigate the impact of cyberattacks on the vehicle-to-infrastructure (V2I) network from a V2I application point of view. Then, we develop a novel V2I cybersecurity architecture, named CVGuard, which can detect and prevent cyberattacks on the V2I environment. In designing CVGuard, key challenges, such as scalability, resiliency and future usability were considered. A case study using a distributed denial of service (DDoS) on a V2I application, i.e., the Stop Sign Gap Assist (SSGA) application, shows that CVGuard was effective in mitigating the adverse effects created by a DDoS attack. In our case study, because of the DDoS attack, conflicts between the minor and major road vehicles occurred in an unsignalized intersection, which could have caused potential crashes. A reduction of conflicts between vehicles occurred because CVGuard was in operation. The reduction of conflicts was compared based on the number of conflicts before and after the implementation and operation of the CVGuard security platform. Analysis revealed that the strategies adopted by the CVGuard were successful in reducing the inter-vehicle conflicts by 60% where a DDoS attack compromised the SSGA application at an unsignalized intersection.


title: Online Knapsack Problem under Expected Capacity Constraint

abstract: Online knapsack problem is considered, where items arrive in a sequential fashion that have two attributes; value and weight. Each arriving item has to be accepted or rejected on its arrival irrevocably. The objective is to maximize the sum of the value of the accepted items such that the sum of their weights is below a budget/capacity. Conventionally a hard budget/capacity constraint is considered, for which variety of results are available. In modern applications, e.g., in wireless networks, data centres, cloud computing, etc., enforcing the capacity constraint in expectation is sufficient. With this motivation, we consider the knapsack problem with an expected capacity constraint. For the special case of knapsack problem, called the secretary problem, where the weight of each item is unity, we propose an algorithm whose probability of selecting any one of the optimal items is equal to $1-1/e$ and provide a matching lower bound. For the general knapsack problem, we propose an algorithm whose competitive ratio is shown to be $1/4e$ that is significantly better than the best known competitive ratio of $1/10e$ for the knapsack problem with the hard capacity constraint.


title: Promise of Commercialization: Carbon Materials for Low-Cost Perovskite Solar Cells

abstract: Perovskite solar cells (PVSCs) have attracted extensive studies due to their high power conversion efficiency (PCE) with low-cost in both raw material and processes. However, there remain obstacles that hinder the way to its commercialization. Among many drawbacks in PVSCs, we note the problems brought by the use of noble metal counter electrodes (CEs) such as gold (Au) and silver (Ag). The costly Au and Ag need high energy-consumption thermal evaporation process which can be made only with expensive evaporation equipment under vacuum. All the factors elevate the threshold of PVSCs' commercialization. Carbon material, on the other hand, is a readily available electrode candidate for the application as CE in the PVSCs. In this review, endeavors on PVSCs with low-cost carbon materials will be comprehensively discussed based on different device structures and carbon composition. We believe that the PVSCs with carbon-based CE hold the promise of commercialization of this new technology.


title: Augmented Outcome-weighted Learning for Optimal Treatment Regimes

abstract: Precision medicine is of considerable interest in clinical, academic and regulatory parties. The key to precision medicine is the optimal treatment regime. Recently, Zhou et al. (2017) developed residual weighted learning (RWL) to construct the optimal regime that directly optimize the clinical outcome. However, this method involves computationally intensive non-convex optimization, which cannot guarantee a global solution. Furthermore, this method does not possess fully semiparametrical efficiency. In this article, we propose augmented outcome-weighted learning (AOL). The method is built on a doubly robust augmented inverse probability weighted estimator, and hence constructs semiparametrically efficient regimes. Our proposed AOL is closely related to RWL. The weights are obtained from counterfactual residuals, where negative residuals are reflected to positive and accordingly their treatment assignments are switched to opposites. Convex loss functions are thus applied to guarantee a global solution and to reduce computations. We show that AOL is universally consistent, i.e., the estimated regime of AOL converges the Bayes regime when the sample size approaches infinity, without knowing any specifics of the distribution of the data. We also propose variable selection methods for linear and nonlinear regimes, respectively, to further improve performance. The performance of the proposed AOL methods is illustrated in simulation studies and in an analysis of the Nefazodone-CBASP clinical trial data.


title: Semi-classical Solutions For Fractional Schrodinger Equations With Potential Vanishing At Infinity

abstract: We study the following fractional Schr\"{o}dinger equation \begin{equation}\label{eq0.1} \varepsilon^{2s}(-\Delta)^s u + Vu = |u|^{p - 2}u,\ \ x\in\,\,\mathbb{R}^N. \end{equation} We show that if the external potential $V\in C(\mathbb{R}^N;[0,\infty))$ has a local minimum and $p\in (2 + 2s/(N - 2s), 2^*_s)$, where $2^*_s=2N/(N-2s),\,N\ge 2s$, the problem has a family of solutions concentrating at the local minimum of $V$ provided that $\liminf_{|x|\to \infty}V(x)|x|^{2s} > 0$. The proof is based on variational methods and penalized technique. {\textbf {Key words}: } fractional Schr\"{o}dinger; vanishing potential; penalized technique; variational methods.


title: Non-orthogonal Multiple Access Assisted Multi-Region Geocast

abstract: This work proposes a novel location-based multi-group multicast framework which is termed as non-orthogonal multiple access (NOMA) assisted multi-region geocast. This novel spectrum sharing framework exploits the NOMA technology to realize the simultaneous delivery of different messages to different user groups characterized by different geographical locations. The essence of the proposed framework is that the geographical information of user groups unites NOMA and multi-group multicast to enhance the spectral efficiency (SE) and the energy efficiency (EE) of wireless transmissions. Specifically, we investigate the downlink beamforming design of the proposed framework in multiple-input single-output (MISO) settings. The decoding strategy for NOMA is designed and guaranteed by users' geographical information and required quality of service (QoS). The majorization and minimization (MM) algorithm is exploited to solve the non-convex and intractable problems therein. Comprehensive numerical experiments are further provided to show that NOMA holds tremendous promise but also limitations in terms of SE and EE, compared with spatial division multiple access (SDMA) and orthogonal multiple access (OMA).


title: Ground-based Characterization of Hayabusa2 Mission Target Asteroid 162173 Ryugu: Constraining Mineralogical Composition in Preparation for Spacecraft Operations

abstract: Asteroids that are targets of spacecraft missions are interesting because they present us with an opportunity to validate ground-based spectral observations. One such object is near-Earth asteroid (NEA) (162173) Ryugu, which is the target of the Japanese Space Agency's (JAXA) Hayabusa2 sample return mission. We observed Ryugu using the 3-m NASA Infrared Telescope Facility (IRTF) on Mauna Kea, Hawaii, on July 13, 2016 to constrain the object's surface composition, meteorite analogs, and link to other asteroids in the main belt and NEA populations. We also modeled its photometric properties using archival data. Using the Lommel-Seeliger model we computed the predicted flux for Ryugu at a wide range of viewing geometries as well as albedo quantities such as geometric albedo, phase integral, and spherical Bond albedo. Our computed albedo quantities are consistent with results from Ishiguro et al. (2014). Our spectral analysis has found a near-perfect match between our spectrum of Ryugu and those of NEA (85275) 1994 LY and Mars-crossing asteroid (316720) 1998 BE7, suggesting that their surface regoliths have similar composition. We compared Ryugu's spectrum with that of main belt asteroid (302) Clarissa, the largest asteroid in the Clarissa asteroid family, suggested as a possible source of Ryugu by Campins et al. (2013). We found that the spectrum of Clarissa shows significant differences with our spectrum of Ryugu, but it is similar to the spectrum obtained by Moskovitz et al. (2013). The best possible meteorite analogs for our spectrum of Ryugu are two CM2 carbonaceous chondrites, Mighei and ALH83100.


title: Deep-Person: Learning Discriminative Deep Features for Person Re-Identification

abstract: Recently, many methods of person re-identification (Re-ID) rely on part-based feature representation to learn a discriminative pedestrian descriptor. However, the spatial context between these parts is ignored for the independent extractor to each separate part. In this paper, we propose to apply Long Short-Term Memory (LSTM) in an end-to-end way to model the pedestrian, seen as a sequence of body parts from head to foot. Integrating the contextual information strengthens the discriminative ability of local representation. We also leverage the complementary information between local and global feature. Furthermore, we integrate both identification task and ranking task in one network, where a discriminative embedding and a similarity measurement are learned concurrently. This results in a novel three-branch framework named Deep-Person, which learns highly discriminative features for person Re-ID. Experimental results demonstrate that Deep-Person outperforms the state-of-the-art methods by a large margin on three challenging datasets including Market-1501, CUHK03, and DukeMTMC-reID. Specifically, combining with a re-ranking approach, we achieve a 90.84% mAP on Market-1501 under single query setting.


title: Co-rank of weakly parafree $3$-manifold groups

abstract: Recall that a group is called large if it has a finite index subgroup which surjects onto a non-abelian free group. By work of Agol and Cooper-Long-Reid, most 3-manifold groups are large; in particular, the fundamental groups of hyperbolic 3-manifolds are large. In previous work, the first author gave examples of closed, hyperbolic 3-manifolds with arbitrarily large first homology rank but whose fundamental groups do not surject onto a non-abelian free group. We call a group very large if it surjects onto a non-abelian free group. In this paper, we consider the question of whether the groups of homology handlebodies - which are very close to being free - are very large. We show that the fundamental group of W. Thurston's tripus manifold, is not very large; it is known to be weakly parafree by Stallings' Theorem and large by the work of Cooper-Long-Reid since the tripus is a hyperbolic manifold with totally geodesic boundary. It is still unknown if a 3-manifold group that is weakly parafree of rank at least 3 must be very large. However, we more generally consider the co-rank of the fundamental group, also known as the cut number of the manifold. For each positive integer g we construct a homology handlebody Y_g of genus g whose group has co-rank equal to r(g), where r(g)=g/2 for g even and r(g)=(g+1)/2 for g odd. That is, these groups are weakly parafree of rank g and surject onto a free group of rank roughly half of g but no larger.


title: Decoding infrared imprints of quantum origins of black holes

abstract: We analyze the emission spectrum of a (fundamentally quantum) black hole in the Kerr-Newman family by assuming a discretization of black hole geometry and the holographic entropy-area relation. We demonstrate that, given the above structure of black hole entropy, a macroscopic black hole always has non-continuously separated mass states and therefore they descend down in discrete manner. We evaluate the step size of the discrete spectrum, which vanishes in the extremal limit, leading to a continuum spectrum as expected from thermal nature of black holes. This further reveals an interesting relation, in each class, between the dynamic and kinematic length scales for all black holes belonging to the Kerr-Newman family, pointing towards a possible universal character across the class, dependent only on black hole mass. Further, we have presented the computation of maximum number of emitted quanta from the black hole as well as an estimation of its lifetime. We also argue the independence of these features from the presence of additional spacetime dimensions.


title: Inner Product and Set Disjointness: Beyond Logarithmically Many Parties

abstract: A basic goal in complexity theory is to understand the communication complexity of number-on-the-forehead problems $f\colon(\{0,1\}^n)^{k}\to\{0,1\}$ with $k\gg\log n$ parties. We study the problems of inner product and set disjointness and determine their randomized communication complexity for every $k\geq\log n$, showing in both cases that $\Theta(1+\lceil\log n\rceil/\log\lceil1+k/\log n\rceil)$ bits are necessary and sufficient. In particular, these problems admit constant-cost protocols if and only if the number of parties is $k\geq n^{\epsilon}$ for some constant $\epsilon>0.$


title: An Adaptive Fuzzy-Based System to Simulate, Quantify and Compensate Color Blindness

abstract: About 8% of the male population of the world are affected by a determined type of color vision disturbance, which varies from the partial to complete reduction of the ability to distinguish certain colors. A considerable amount of color blind people are able to live all life long without knowing they have color vision disabilities and abnormalities. Nowadays the evolution of information technology and computer science, specifically image processing techniques and computer graphics, can be fundamental to aid at the development of adaptive color blindness correction tools. This paper presents a software tool based on Fuzzy Logic to evaluate the type and the degree of color blindness a person suffer from. In order to model several degrees of color blindness, herein this work we modified the classical linear transform-based simulation method by the use of fuzzy parameters. We also proposed four new methods to correct color blindness based on a fuzzy approach: Methods A and B, with and without histogram equalization. All the methods are based on combinations of linear transforms and histogram operations. In order to evaluate the results we implemented a web-based survey to get the best results according to optimize to distinguish different elements in an image. Results obtained from 40 volunteers proved that the Method B with histogram equalization got the best results for about 47% of volunteers.


title: Predicting readmission risk from doctors' notes

abstract: We develop a model using deep learning techniques and natural language processing on unstructured text from medical records to predict hospital-wide $30$-day unplanned readmission, with c-statistic $.70$. Our model is constructed to allow physicians to interpret the significant features for prediction.


title: $s$-homogeneous algebras via $s$-homogeneous triples

abstract: To study $s$-homogeneous algebras, we introduce the category of quivers with $s$-homogeneous corelations and the category of $s$-homogeneous triples. We show that both of these categories are equivalent to the category of $s$-homogeneous algebras. We prove some properties of the elements of $s$-homogeneous triples and give some consequences for $s$-Koszul algebras. Then we discuss the relations between the $s$-Koszulity and the Hilbert series of $s$-homogeneous triples. We give some application of the obtained results to $s$-homogeneous algebras with simple zero component. We describe all $s$-Koszul algebras with one relation recovering the result of Berger and all $s$-Koszul algebras with one dimensional $s$-th component. We show that if the $s$-th Veronese ring of an $s$-homogeneous algebra has two generators, then it has at least two relations. Finally, we classify all $s$-homogeneous algebras with $s$-th Veronese rings ${\bf k}\langle x,y\rangle/(xy,yx)$ and ${\bf k}\langle x,y\rangle/(x^2,y^2)$. In particular, we show that all of these algebras are not $s$-Koszul while their $s$-homogeneous duals are $s$-Koszul.


title: Cost-Effective Seed Selection in Online Social Networks

abstract: We study the min-cost seed selection problem in online social networks, where the goal is to select a set of seed nodes with the minimum total cost such that the expected number of influenced nodes in the network exceeds a predefined threshold. We propose several algorithms that outperform the previous studies both on the theoretical approximation ratios and on the experimental performance. Under the case where the nodes have heterogeneous costs, our algorithms are the first bi- criteria approximation algorithms with polynomial running time and provable logarithmic performance bounds using a general contagion model. Under the case where the users have uniform costs, our algorithms achieve logarithmic approximation ratio and provable time complexity which is smaller than that of existing algorithms in orders of magnitude. We conduct extensive experiments using real social networks. The experimental results show that, our algorithms significantly outperform the existing algorithms both on the total cost and on the running time, and also scale well to billion-scale networks.


title: Absence of Chaotic Size Dependence for Spin Glasses on Hierarchical Lattices

abstract: We investigate the question of whether chaotic size dependence occurs on hierarichical lattices and demonstrate that it is not present in these systems. Our results show that the metastate for spin glasses on hierarchical lattices is simple.


title: Laplacian Controllability of Threshold Graphs

abstract: This paper is concerned with the controllability problem of a connected threshold graph following the Laplacian dynamics. An algorithm is proposed to generate a spanning set of orthogonal Laplacian eigenvectors of the graph from a straightforward computation on its Laplacian matrix. A necessary and sufficient condition for the graph to be Laplacian controllable is then proposed. The condition suggests that the minimum number of controllers to render a connected threshold graph controllable is the maximum multiplicity of entries in the conjugate of the degree sequence determining the graph, and this minimum can be achieved by a binary control matrix. The second part of the work is the introduction of a novel class of single-input controllable graphs, which is constructed by connecting two antiregular graphs with almost the same size. This new connecting structure reduces the sum of the maximum vertex degree and the diameter by almost one half, compared to other well-known single-input controllable graphs such as the path and the antiregular graph, and has potential applications in the design of controllable graphs subject to practical edge constraints. Examples are provided to illustrate our results.


title: Closed-form expressions for effective constitutive parameters and electro/magneto-strictive tensors for bi-anisotropic metamaterials and their use in optical force density calculations

abstract: Using a multiple scattering technique, we derived closed-form expressions for effective constitutive parameters and electro/magneto-strictive tensor components for 2D bi-anisotropic metamaterials. Using the principle of virtual work, we obtained the electromagnetic stress tensor that can be used to calculate the optical force density inside such media. The analytic expressions are tested against full wave numerical simulations. Our effective medium theory is essential for providing a complete macroscopic description of the optical and opto-mechanical properties of bi-anisotropic composites.


title: Image2Mesh: A Learning Framework for Single Image 3D Reconstruction

abstract: One challenge that remains open in 3D deep learning is how to efficiently represent 3D data to feed deep networks. Recent works have relied on volumetric or point cloud representations, but such approaches suffer from a number of issues such as computational complexity, unordered data, and lack of finer geometry. This paper demonstrates that a mesh representation (i.e. vertices and faces to form polygonal surfaces) is able to capture fine-grained geometry for 3D reconstruction tasks. A mesh however is also unstructured data similar to point clouds. We address this problem by proposing a learning framework to infer the parameters of a compact mesh representation rather than learning from the mesh itself. This compact representation encodes a mesh using free-form deformation and a sparse linear combination of models allowing us to reconstruct 3D meshes from single images. In contrast to prior work, we do not rely on silhouettes and landmarks to perform 3D reconstruction. We evaluate our method on synthetic and real-world datasets with very promising results. Our framework efficiently reconstructs 3D objects in a low-dimensional way while preserving its important geometrical aspects.


title: The game of plates and olives

abstract: The game of plates and olives, introduced by Nicolaescu, begins with an empty table. At each step either an empty plate is put down, an olive is put down on a plate, an olive is removed, an empty plate is removed, or the olives on two plates that both have olives on them are combined on one of the two plates, with the other plate removed. Plates are indistinguishable from one another, as are olives, and there is an inexhaustible supply of each. The game derives from the consideration of Morse functions on the $2$-sphere. Specifically, the number of topological equivalence classes of excellent Morse functions on the $2$-sphere that have order $n$ (that is, that have $2n+2$ critical points) is the same as the number of ways of returning to an empty table for the first time after exactly $2n+2$ steps. We call this number $M_n$. Nicolaescu gave the lower bound $M_n \geq (2n-1)!! = (2/e)^{n+o(n)}n^n$ and speculated that $\log M_n \sim n\log n$. In this note we confirm this speculation, showing that $M_n \leq (4/e)^{n+o(n)}n^n$.


title: $\mathbb{F}_{q}[G]$-modules and $G$-invariant codes

abstract: If $\mathbb{F}_{q}$ is a finite field, $C$ is a vector subspace of $\mathbb{F}_{q}^{n}$ (linear code), and $G$ is a subgroup of the group of linear automorphisms of $\mathbb{F}_{q}^{n}$, $C$ is said to be $G$-invariant if $g(C)=C$ for all $g\in G$. A solution to the problem of computing all the $G$-invariant linear codes $C$ of $\mathbb{F}_{q}^{n}$ is offered. This will be referred as the invariance problem. When $n=|G|t$, we determine conditions for the existence of an isomorphism of $\mathbb{F}_{q}[G]$-modules between $\mathbb{F}_{q}^{n}$ and $\mathbb{F}_{q}[G]\times \cdots \times \mathbb{F}_{q}[G]$ ($t$-times), that preserves the Hamming weight. This reduces the invariance problem to the determination of the $\mathbb{F}_{q}[G]$-submodules of $\mathbb{F}_{q}[G]\times \cdots \times \mathbb{F}_{q}[G]$ ($t$-times). The concept of Gaussian binomial coefficient for semisimple $\mathbb{F}_{q}[G]$-modules, which is useful for counting $G$-invariant codes, is introduced. Finally, a systematic way to compute all the $G$-invariant linear codes $C\subseteq \mathbb{F}_{q}^{n}$ is provided, when $(|G|,q)=1$.


title: Invasion Percolation on Galton-Watson Trees

abstract: We consider invasion percolation on Galton-Watson trees. On almost every Galton-Watson tree, the invasion cluster almost surely contains only one infinite path. This means that for almost every Galton-Watson tree, invasion percolation induces a probability measure on infinite paths from the root. We show that under certain conditions of the progeny distribution, this measure is absolutely continuous with respect to the limit uniform measure. This confirms that invasion percolation, an efficient self-tuning algorithm, may be used to sample approximately from the limit uniform distribution. Additionally, we analyze the forward maximal weights along the backbone of the invasion cluster and prove a limit law for the process.


title: Energy-Efficient Time-Domain Vector-by-Matrix Multiplier for Neurocomputing and Beyond

abstract: We propose an extremely energy-efficient mixed-signal approach for performing vector-by-matrix multiplication in a time domain. In such implementation, multi-bit values of the input and output vector elements are represented with time-encoded digital signals, while multi-bit matrix weights are realized with current sources, e.g. transistors biased in subthreshold regime. With our approach, multipliers can be chained together to implement large-scale circuits completely in a time domain. Multiplier operation does not rely on energy-taxing static currents, which are typical for peripheral and input/output conversion circuits of the conventional mixed-signal implementations. As a case study, we have designed a multilayer perceptron, based on two layers of 10x10 four-quadrant vector-by-matrix multipliers, in 55-nm process with embedded NOR flash memory technology, which allows for compact implementation of adjustable current sources. Our analysis, based on memory cell measurements, shows that at high computing speed the drain-induced barrier lowering is a major factor limiting multiplier precision to ~6 bit. Post-layout estimates for a conservative 6-bit digital input/output NxN multiplier designed in 55 nm process, including I/O circuitry for converting between digital and time domain representations, show ~7 fJ/Op for N>200, which can be further lowered well below 1 fJ/Op for more optimal and aggressive design.


title: Direct Information Reweighted by Contact Templates: Improved RNA Contact Prediction by Combining Structural Features

abstract: It is acknowledged that co-evolutionary nucleotide-nucleotide interactions are essential for RNA structures and functions. Currently, direct coupling analysis (DCA) infers nucleotide contacts in a sequence from its homologous sequence alignment across different species. DCA and similar approaches that use sequence information alone usually yield a low accuracy, especially when the available homologous sequences are limited. Here we present a new method that incorporates a Restricted Boltzmann Machine (RBM) to augment the information on sequence co-variations with structural patterns in contact inference. We thus name our method DIRECT that stands for Direct Information REweighted by Contact Templates. Benchmark tests demonstrate that DIRECT produces a substantial enhancement of 13% in accuracy on average for contact prediction in comparison to the traditional DCA. These results suggest that DIRECT could be used for improving predictions of RNA tertiary structures and functions. The source codes and dataset of DIRECT are available at http:// http://zhao.phy.ccnu.edu.cn:8122/DIRECT/index.html.


title: Logarithmic singularities and quantum oscillations in magnetically doped topological insulators

abstract: We report magnetotransport measurements on magnetically doped (Bi,Sb)$_2$Te$_3$ films grown by molecular beam epitaxy. In Hallbar devices, logarithmic dependence on temperature and bias voltage are obseved in both the longitudinal and anomalous Hall resistance. The interplay of disorder and electron-electron interactions is found to explain quantitatively the observed logarithmic singularities and is a dominant scattering mechanism in these samples. Submicron scale devices exhibit intriguing quantum oscillations at high magnetic fields with dependence on bias voltage. The observed quantum oscillations can be attributed to bulk and surface transport.


title: Entanglement in non-local games and the hyperlinear profile of groups

abstract: We relate the amount of entanglement required to play linear-system non-local games near-optimally to the hyperlinear profile of finitely-presented groups. By calculating the hyperlinear profile of a certain group, we give an example of a finite non-local game for which the amount of entanglement required to play $\varepsilon$-optimally is at least $\Omega(1/\varepsilon^k)$, for some $k>0$. Since this function approaches infinity as $\varepsilon$ approaches zero, this provides a quantitative version of a theorem of the first author.


title: Private federated learning on vertically partitioned data via entity resolution and additively homomorphic encryption

abstract: Consider two data providers, each maintaining private records of different feature sets about common entities. They aim to learn a linear model jointly in a federated setting, namely, data is local and a shared model is trained from locally computed updates. In contrast with most work on distributed learning, in this scenario (i) data is split vertically, i.e. by features, (ii) only one data provider knows the target variable and (iii) entities are not linked across the data providers. Hence, to the challenge of private learning, we add the potentially negative consequences of mistakes in entity resolution. Our contribution is twofold. First, we describe a three-party end-to-end solution in two phases ---privacy-preserving entity resolution and federated logistic regression over messages encrypted with an additively homomorphic scheme---, secure against a honest-but-curious adversary. The system allows learning without either exposing data in the clear or sharing which entities the data providers have in common. Our implementation is as accurate as a naive non-private solution that brings all data in one place, and scales to problems with millions of entities with hundreds of features. Second, we provide what is to our knowledge the first formal analysis of the impact of entity resolution's mistakes on learning, with results on how optimal classifiers, empirical losses, margins and generalisation abilities are affected. Our results bring a clear and strong support for federated learning: under reasonable assumptions on the number and magnitude of entity resolution's mistakes, it can be extremely beneficial to carry out federated learning in the setting where each peer's data provides a significant uplift to the other.


title: AttGAN: Facial Attribute Editing by Only Changing What You Want

abstract: Facial attribute editing aims to manipulate single or multiple attributes of a face image, i.e., to generate a new face with desired attributes while preserving other details. Recently, generative adversarial net (GAN) and encoder-decoder architecture are usually incorporated to handle this task with promising results. Based on the encoder-decoder architecture, facial attribute editing is achieved by decoding the latent representation of the given face conditioned on the desired attributes. Some existing methods attempt to establish an attribute-independent latent representation for further attribute editing. However, such attribute-independent constraint on the latent representation is excessive because it restricts the capacity of the latent representation and may result in information loss, leading to over-smooth and distorted generation. Instead of imposing constraints on the latent representation, in this work we apply an attribute classification constraint to the generated image to just guarantee the correct change of desired attributes, i.e., to "change what you want". Meanwhile, the reconstruction learning is introduced to preserve attribute-excluding details, in other words, to "only change what you want". Besides, the adversarial learning is employed for visually realistic editing. These three components cooperate with each other forming an effective framework for high quality facial attribute editing, referred as AttGAN. Furthermore, our method is also directly applicable for attribute intensity control and can be naturally extended for attribute style manipulation. Experiments on CelebA dataset show that our method outperforms the state-of-the-arts on realistic attribute editing with facial details well preserved.


title: Multiplexing induced explosive synchronization in Kuramoto oscillators with inertia

abstract: Explosive synchronization (ES) of coupled oscillators on networks is shown to be originated from existence of correlation between natural frequencies of oscillators and degrees of corresponding nodes. Here, we demonstrate that ES is a generic feature of multiplex network of second-order Kuramoto oscillators and can exist in absence of a frequency-degree correlation. A monoplex network of second-order Kuramoto oscillators bearing homogeneous (heterogeneous) degree-distribution is known to display the first-order (second-order) transition to synchronization. We report that multiplexing of two such networks having homogeneous degree-distribution support the first-order transition in both the layers thereby facilitating ES. More interesting is the multiplexing of a layer bearing heterogeneous degree-distribution with another layer bearing homogeneous degree-distribution, which induces a first-order (ES) transition in the heterogeneous layer which was incapable of showing the same in the isolation. Further, we report that such induced ES transition in the heterogeneous layer of multiplex networks can be controlled by varying inter and intra-layer coupling strengths. Our findings emphasize on importance of multiplexing or impact of one layer on dynamical evolution of other layers of systems having inherent multiplex or multilevel architecture.


title: On the Approximability and Hardness of the Minimum Connected Dominating Set with Routing Cost Constraint

abstract: In the problem of minimum connected dominating set with routing cost constraint, we are given a graph $G=(V,E)$, and the goal is to find the smallest connected dominating set $D$ of $G$ such that, for any two non-adjacent vertices $u$ and $v$ in $G$, the number of internal nodes on the shortest path between $u$ and $v$ in the subgraph of $G$ induced by $D \cup \{u,v\}$ is at most $\alpha$ times that in $G$. For general graphs, the only known previous approximability result is an $O(\log n)$-approximation algorithm ($n=|V|$) for $\alpha = 1$ by Ding et al. For any constant $\alpha > 1$, we give an $O(n^{1-\frac{1}{\alpha}}(\log n)^{\frac{1}{\alpha}})$-approximation algorithm. When $\alpha \geq 5$, we give an $O(\sqrt{n}\log n)$-approximation algorithm. Finally, we prove that, when $\alpha =2$, unless $NP \subseteq DTIME(n^{poly\log n})$, for any constant $\epsilon > 0$, the problem admits no polynomial-time $2^{\log^{1-\epsilon}n}$-approximation algorithm, improving upon the $\Omega(\log n)$ bound by Du et al. (albeit under a stronger hardness assumption).


title: Formation and distribution of fragments in the spontaneous fission of 240Pu

abstract: The goal of this paper is to better understand the structure of fission fragment distributions by investigating the competition between the static structure of the collective manifold and stochastic dynamics. In particular, we study the characteristics of the tails of yield distributions, which correspond to very asymmetric fission. We use the stochastic Langevin framework to simulate the nuclear evolution after the system tunnels through the multi-dimensional potential barrier. For a representative sample of different initial configurations along the outer turning-point line, we define effective fission paths by computing a large number of Langevin trajectories. We extract the relative contribution of each such path to the fragment distribution. We then use nucleon localization functions along effective fission pathways to analyze the characteristics of prefragments at pre-scission configurations. We find that non-Newtonian Langevin trajectories, strongly impacted by the random force, produce the tails of the fission fragment distribution. The prefragments deduced from nucleon localizations are formed early and change little as the nucleus evolves towards scission. On the other hand, the system contains many nucleons that are not localized in the prefragments, even near the scission point. Such nucleons are rapidly distributed at scission to form the final fragments. Our study shows that only theoretical models of fission that account for some form of dissipative/stochastic dynamics can give an accurate description of the structure of fragment distributions. In particular, it should be nearly impossible to predict the tails of these distributions within the standard formulation of time-dependent density functional theory. At the same time, the large number of non-localized nucleons during fission suggests that adiabatic approaches are ill-suited to describe fission fragment properties.


title: Haar wavelet quasilinearization technique for doubly singular boundary value problems

abstract: The Haar wavelet based quasilinearization technique for solving a general class of singular boundary value problems is proposed. Quasilinearization technique is used to linearize nonlinear singular problem. Second rate of convergence is obtained of a sequence of linear singular problems. Numerical solution of linear singular prob- lems is obtained by Haar-wavelet method. In each iteration of quasilinearization technique, the numerical solution is updated by the Haar wavelet method. Conver- gence analysis of Haar wavelet method is discussed. The results are compared with the results obtained by the other technique and with exact solution. Eight singular problems are solved to show the applicability of the Haar wavelet quasilinearization technique.


title: Patch Correspondences for Interpreting Pixel-level CNNs

abstract: We present compositional nearest neighbors (CompNN), a simple approach to visually interpreting distributed representations learned by a convolutional neural network (CNN) for pixel-level tasks (e.g., image synthesis and segmentation). It does so by reconstructing both a CNN's input and output image by copy-pasting corresponding patches from the training set with similar feature embeddings. To do so efficiently, it makes of a patch-match-based algorithm that exploits the fact that the patch representations learned by a CNN for pixel level tasks vary smoothly. Finally, we show that CompNN can be used to establish semantic correspondences between two images and control properties of the output image by modifying the images contained in the training set. We present qualitative and quantitative experiments for semantic segmentation and image-to-image translation that demonstrate that CompNN is a good tool for interpreting the embeddings learned by pixel-level CNNs.


title: Road Extraction by Deep Residual U-Net

abstract: Road extraction from aerial images has been a hot research topic in the field of remote sensing image analysis. In this letter, a semantic segmentation neural network which combines the strengths of residual learning and U-Net is proposed for road area extraction. The network is built with residual units and has similar architecture to that of U-Net. The benefits of this model is two-fold: first, residual units ease training of deep networks. Second, the rich skip connections within the network could facilitate information propagation, allowing us to design networks with fewer parameters however better performance. We test our network on a public road dataset and compare it with U-Net and other two state of the art deep learning based road extraction methods. The proposed approach outperforms all the comparing methods, which demonstrates its superiority over recently developed state of the arts.


title: IoT based Platform as a Service for Provisioning of Concurrent Applications

abstract: The modern era has seen a speedy growth in the Internet of Things (IoT). As per statistics of 2020, twenty billion devices will be connected to the Internet. This massive increase in Internet connected devices will lead to a lot of efforts to execute critical concurrent applications such fire detection, health care based system, disaster management, high energy physics, automobiles, and medical imaging efficiently. To fasten the emergence of novel applications, this vast infrastructure requires "Platform as a Service(PaaS)" model to leverage IoT things. As a single global standard for all device types and IoT-based application domain is impracticable, we propose an IoT-based Cloud to leverage PaaS model in this paper. This model can host the concurrent application for Wireless Sensor Network (WSN). The proposed model offers the communication interface among processes by uniquely allocating network interface to a particular container.


title: Practical Synchronization Waveform for Massive Machine-Type-Communications

abstract: A general synchronization waveform type that is insensitive to frequency error for massive Internet-of-Things (mIoT) is derived in "Robust synchronization waveform design for massive low-power IoT" (IEEE Trans. on Wireless Commun., vol. 16, no. 11, pp. 7551-7559, Nov. 2017). Detailed issues such as (1) how to obtain the frequency error information from the waveform; (2) how to refine the timing estimation; and (3) how to adjust the waveform parameter to meet various spectral requirements are not addressed in the paper, and yet are of great importance for a practical application. These issues are discussed in this comment.


title: Chance-Constrained Day-Ahead Hourly Scheduling in Distribution System Operation

abstract: This paper aims to propose a two-step approach for day-ahead hourly scheduling in a distribution system operation, which contains two operation costs, the operation cost at substation level and feeder level. In the first step, the objective is to minimize the electric power purchase from the day-ahead market with the stochastic optimization. The historical data of day-ahead hourly electric power consumption is used to provide the forecast results with the forecasting error, which is presented by a chance constraint and formulated into a deterministic form by Gaussian mixture model (GMM). In the second step, the objective is to minimize the system loss. Considering the nonconvexity of the three-phase balanced AC optimal power flow problem in distribution systems, the second-order cone program (SOCP) is used to relax the problem. Then, a distributed optimization approach is built based on the alternating direction method of multiplier (ADMM). The results shows that the validity and effectiveness method.


title: Facial Dynamics Interpreter Network: What are the Important Relations between Local Dynamics for Facial Trait Estimation?

abstract: Human face analysis is an important task in computer vision. According to cognitive-psychological studies, facial dynamics could provide crucial cues for face analysis. The motion of a facial local region in facial expression is related to the motion of other facial local regions. In this paper, a novel deep learning approach, named facial dynamics interpreter network, has been proposed to interpret the important relations between local dynamics for estimating facial traits from expression sequence. The facial dynamics interpreter network is designed to be able to encode a relational importance, which is used for interpreting the relation between facial local dynamics and estimating facial traits. By comparative experiments, the effectiveness of the proposed method has been verified. The important relations between facial local dynamics are investigated by the proposed facial dynamics interpreter network in gender classification and age estimation. Moreover, experimental results show that the proposed method outperforms the state-of-the-art methods in gender classification and age estimation.


title: Unified theory for finite Markov chains

abstract: We provide a unified framework to compute the stationary distribution of any finite irreducible Markov chain or equivalently of any irreducible random walk on a finite semigroup $S$. Our methods use geometric finite semigroup theory via the Karnofsky-Rhodes and the McCammond expansions of finite semigroups with specified generators; this does not involve any linear algebra. The original Tsetlin library is obtained by applying the expansions to $P(n)$, the set of all subsets of an $n$ element set. Our set-up generalizes previous groundbreaking work involving left-regular bands (or $\mathscr{R}$-trivial bands) by Brown and Diaconis, extensions to $\mathscr{R}$-trivial semigroups by Ayyer, Steinberg, Thi\'ery and the second author, and important recent work by Chung and Graham. The Karnofsky-Rhodes expansion of the right Cayley graph of $S$ in terms of generators yields again a right Cayley graph. The McCammond expansion provides normal forms for elements in the expanded $S$. Using our previous results with Silva based on work by Berstel, Perrin, Reutenauer, we construct (infinite) semaphore codes on which we can define Markov chains. These semaphore codes can be lumped using geometric semigroup theory. Using normal forms and associated Kleene expressions, they yield formulas for the stationary distribution of the finite Markov chain of the expanded $S$ and the original $S$. Analyzing the normal forms also provides an estimate on the mixing time.


title: Load Forecasting Based Distribution System Network Reconfiguration-A Distributed Data-Driven Approach

abstract: In this paper, a short-term load forecasting approach based network reconfiguration is proposed in a parallel manner. Specifically, a support vector regression (SVR) based short-term load forecasting approach is designed to provide an accurate load prediction and benefit the network reconfiguration. Because of the nonconvexity of the three-phase balanced optimal power flow, a second-order cone program (SOCP) based approach is used to relax the optimal power flow problem. Then, the alternating direction method of multipliers (ADMM) is used to compute the optimal power flow in distributed manner. Considering the limited number of the switches and the increasing computation capability, the proposed network reconfiguration is solved in a parallel way. The numerical results demonstrate the feasible and effectiveness of the proposed approach.


title: Orthogonal and symplectic Harish-Chandra integrals and matrix product ensembles

abstract: In this paper, we highlight the role played by orthogonal and symplectic Harish-Chandra integrals in the study of real-valued matrix product ensembles. By making use of these integrals and the matrix-valued Fourier-Laplace transform, we find the explicit eigenvalue distributions for particular Hermitian anti-symmetric matrices and particular Hermitian anti-self dual matrices, involving both sums and products. As a consequence of these results, the eigenvalue probability density function of the random product structure $X_M \cdots X_1( iA) X_1^T \cdots X_M^T$, where each $X_i$ is a standard real Gaussian matrix, and $A$ is a real anti-symmetric matrix can be determined. For $M=1$ and $A$ the bidiagonal anti-symmetric matrix with 1's above the diagonal, this reclaims results of Defosseux. For general $M$, and this choice of $A$, or $A$ itself a standard Gaussian anti-symmetric matrix, the eigenvalue distribution is shown to coincide with that of the squared singular values for the product of certain complex Gaussian matrices first studied by Akemann et al. As a point of independent interest, we also include a self-contained diffusion equation derivation of the orthogonal and symplectic Harish-Chandra integrals.


title: Local Access to Huge Random Objects through Partial Sampling

abstract: Consider an algorithm performing a computation on a huge random object. Is it necessary to generate the entire object up front, or is it possible to provide query access to the object and sample it incrementally "on-the-fly"? Such an implementation should emulate the object by answering queries in a manner consistent with a random instance sampled from the true distribution. Our first set of results focus on undirected graphs with independent edge probabilities, under certain assumptions. Then, we use this to obtain the first efficient implementations for the Erdos-Renyi model and the Stochastic Block model. As in previous local-access implementations for random graphs, we support Vertex-Pair and Next-Neighbor queries. We also introduce a new Random-Neighbor query. Next, we show how to implement random Catalan objects, specifically focusing on Dyck paths (always positive random walks on the integer line). Here, we support Height queries to find the position of the walk, and First-Return queries to find the time when the walk returns to a specified height. This in turn can be used to implement Next-Neighbor queries on random rooted/binary trees, and Matching-Bracket queries on random well bracketed expressions. Finally, we define a new model that: (1) allows multiple independent simultaneous instantiations of the same implementation to be consistent with each other without communication (2) allows us to generate a richer class of random objects that do not have a succinct description. Specifically, we study uniformly random valid $q$-colorings of an input graph $G$ with max degree $\Delta$. The distribution over valid colorings is specified via a "huge" underlying graph $G$, that is far too large to be read in sub-linear time. Instead, we access $G$ through local neighborhood probes. We are able to answer queries to the color of any vertex in sub-linear time for $q > 9\Delta$.


title: Small Drone Field Experiment: Data Collection & Processing

abstract: Following an initiative formalized in April 2016 formally known as ARL West between the U.S. Army Research Laboratory (ARL) and University of Southern California's Institute for Creative Technologies (USC ICT), a field experiment was coordinated and executed in the summer of 2016 by ARL, USC ICT, and Headwall Photonics. The purpose was to image part of the USC main campus in Los Angeles, USA, using two portable COTS (commercial off the shelf) aerial drone solutions for data acquisition, for photogrammetry (3D reconstruction from images), and fusion of hyperspectral data with the recovered set of 3D point clouds representing the target area. The research aims for determining the viability of having a machine capable of segmenting the target area into key material classes (e.g., manmade structures, live vegetation, water) for use in multiple purposes, to include providing the user with a more accurate scene understanding and enabling the unsupervised automatic sampling of meaningful material classes from the target area for adaptive semi-supervised machine learning. In the latter, a target set library may be used for automatic machine training with data of local material classes, as an example, to increase the prediction chances of machines recognizing targets. The field experiment and associated data post processing approach to correct for reflectance, geo-rectify, recover the area's dense point clouds from images, register spectral with elevation properties of scene surfaces from the independently collected datasets, and generate the desired scene segmented maps are discussed. Lessons learned from the experience are also highlighted throughout the paper.


title: Backscatter Multiplicative Multiple-Access Systems: Fundamental Limits and Practical Design

abstract: In this paper, we consider a novel ambient backscatter multiple-access system, where a receiver (Rx) simultaneously detects the signals transmitted from an active transmitter (Tx) and a backscatter Tag. Specifically, the information-carrying signal sent by the Tx arrives at the Rx through two wireless channels: the direct channel from the Tx to the Rx, and the backscatter channel from the Tx to the Tag and then to the Rx. The received signal from the backscatter channel also carries the Tag's information because of the multiplicative backscatter operation at the Tag. This multiple-access system introduces a new channel model, referred to as multiplicative multiple-access channel (M-MAC). We analyze the achievable rate region of the M-MAC, and prove that its region is strictly larger than that of the conventional time-division multiple-access scheme in many cases, including, e.g., the high SNR regime and the case when the direct channel is much stronger than the backscatter channel. Hence, the multiplicative multiple-access scheme is an attractive technique to improve the throughput for ambient backscatter communication systems. Moreover, we analyze the detection error rates for coherent and noncoherent modulation schemes adopted by the Tx and the Tag, respectively, in both synchronous and asynchronous scenarios, which further bring interesting insights for practical system design.