title: Bounds for modified Struve functions of the first kind and their ratios

abstract: We obtain a simple two-sided inequality for the ratio $\mathbf{L}_\nu(x)/\mathbf{L}_{\nu-1}(x)$ in terms of the ratio $I_\nu(x)/I_{\nu-1}(x)$, where $\mathbf{L}_\nu(x)$ is the modified Struve function of the first kind and $I_\nu(x)$ is the modified Bessel function of the first kind. This result allows one to use the extensive literature on bounds for $I_\nu(x)/I_{\nu-1}(x)$ to immediately deduce bounds for $\mathbf{L}_\nu(x)/\mathbf{L}_{\nu-1}(x)$. We note some consequences and obtain further bounds for $\mathbf{L}_\nu(x)/\mathbf{L}_{\nu-1}(x)$ by adapting techniques used to bound the ratio $I_\nu(x)/I_{\nu-1}(x)$. We apply these results to obtain new bounds for the condition numbers $x\mathbf{L}_\nu'(x)/\mathbf{L}_\nu(x)$, the ratio $\mathbf{L}_\nu(x)/\mathbf{L}_\nu(y)$ and the modified Struve function $\mathbf{L}_\nu(x)$ itself. Amongst other results, we obtain two-sided inequalities for $x\mathbf{L}_\nu'(x)/\mathbf{L}_\nu(x)$ and $\mathbf{L}_\nu(x)/\mathbf{L}_\nu(y)$ that are given in terms of $xI_\nu'(x)/I_\nu(x)$ and $I_\nu(x)/I_\nu(y)$, respectively, which again allows one to exploit the substantial literature on bounds for these quantities. The results obtained in this paper complement and improve existing bounds in the literature.


title: Graph-based regularization for regression problems with alignment and highly-correlated designs

abstract: Sparse models for high-dimensional linear regression and machine learning have received substantial attention over the past two decades. Model selection, or determining which features or covariates are the best explanatory variables, is critical to the interpretability of a learned model. Much of the current literature assumes that covariates are only mildly correlated. However, in many modern applications covariates are highly correlated and do not exhibit key properties (such as the restricted eigenvalue condition, restricted isometry property, or other related assumptions). This work considers a high-dimensional regression setting in which a graph governs both correlations among the covariates and the similarity among regression coefficients -- meaning there is \emph{alignment} between the covariates and regression coefficients. Using side information about the strength of correlations among features, we form a graph with edge weights corresponding to pairwise covariances. This graph is used to define a graph total variation regularizer that promotes similar weights for correlated features. This work shows how the proposed graph-based regularization yields mean-squared error guarantees for a broad range of covariance graph structures. These guarantees are optimal for many specific covariance graphs, including block and lattice graphs. Our proposed approach outperforms other methods for highly-correlated design in a variety of experiments on synthetic data and real biochemistry data.


title: A New Equation for Activity Calculation in Pulse Irradiation: Derivation, Simulation and Experimental Validation

abstract: To calculate the radioactivity of product nuclides generated in pulse irradiation, it is generally assumed that the irradiation is approximately continues in the entire irradiation period ($t_i$) and the flux of incoming irradiation particle can be obtained by averaging their intensity in each pulse period ($T$). However, this approximation fails to acknowledge the fact that the product nuclides are not created in each pulse period ($T$) evenly: they are only produced in a very short pulse width ($t_p$) and then decay in a relative long rest time ($t_r=T-t_p$). Given by the enormous number of pulses, the sum of these decays may not be negligible. To make the activity calculation in accordance with the real situation in pulse irradiation, we scrutinize the details of irradiation and decay processes in each pulse, applies the geometric series to obtain the activity superimposition of millions of pulses, and derives a novel activity equation particularly suitable for pulse irradiation. The experimental results, numerical simulations, and activity measurements from photon activation driven by a pulsed electron LINAC have confirmed the validity of this new equation. The comparison between the new and traditional equations indicates that their discrepancy could be significant under certain conditions. The limitations of the new activity equation for pulse irradiation are discussed as well.


title: Decentralized decision making for networks of uncertain systems

abstract: Distributed model predictive control (MPC) has been proven a successful method in regulating the operation of large-scale networks of constrained dynamical systems. This paper is concerned with cooperative distributed MPC in which the decision actions of the systems are usually derived by the solution of a system-wide optimization problem. However, formulating and solving such large-scale optimization problems is often a hard task which requires extensive information communication among the individual systems and fails to address privacy concerns in the network. Hence, the main challenge is to design decision policies with a prescribed structure so that the resulting system-wide optimization problem to admit a loosely coupled structure and be amendable to distributed computation algorithms. In this paper, we propose a decentralized problem synthesis scheme which only requires each system to communicate sets which bound its states evolution to neighboring systems. The proposed method alleviates concerns on privacy since this limited communication scheme does not reveal the exact characteristics of the dynamics within each system. In addition, it enables a distributed computation of the solution, making our method highly scalable. We demonstrate in a number of numerical studies, inspired by engineering and finance, the efficacy of the proposed approach which leads to solutions that closely approximate those obtained by the centralized formulation only at a fraction of the computational effort.


title: Efficient Recurrent Neural Networks using Structured Matrices in FPGAs

abstract: Recurrent Neural Networks (RNNs) are becoming increasingly important for time series-related applications which require efficient and real-time implementations. The recent pruning based work ESE suffers from degradation of performance/energy efficiency due to the irregular network structure after pruning. We propose block-circulant matrices for weight matrix representation in RNNs, thereby achieving simultaneous model compression and acceleration. We aim to implement RNNs in FPGA with highest performance and energy efficiency, with certain accuracy requirement (negligible accuracy degradation). Experimental results on actual FPGA deployments shows that the proposed framework achieves a maximum energy efficiency improvement of 35.7$\times$ compared with ESE.


title: Sur les paquets d'Arthur des groupes classiques et unitaires non quasi-d\'eploy\'es

abstract: Nous \'etendons aux groupes orthogonaux et unitaires non quasi-d\'eploy\'es sur un corps local des r\'esultats de J. Arthur et de la premi\`ere auteure \'etablis dans le cas quasi-d\'eploy\'e. En particulier, nous obtenons une classification de Langlands compl\`ete pour les repr\'esentations temp\'er\'ees dans le cas $p$-adique. Nous en d\'eduisons en utilisant l'involution d'Aubert-Schneider-Stuhler un r\'esultat de multiplicit\'e un dans les paquets unipotents, et par des m\'ethodes globales, le m\^eme r\'esultat pour les paquets unipotents dans le cas archim\'edien. We extend to non quasi-split orthogonal and unitary groups over a local field some results of J. Arthur and the first author established in the quasi-split case. In particular, we obtain a full Langlands classification for tempered representations in the $p$-adic case. Using Aubert-Schneider-Stuhler involution, we deduce from this a multiplicity one result for unipotent packets, and by global methods, the same result for unipotent packets in the archimedean case.


title: Impact of dissipation on the energy spectrum of experimental turbulence of gravity surface waves

abstract: We discuss the impact of dissipation on the development of the energy spectrum in wave turbulence of gravity surface waves with emphasis on the effect of surface contamination. We performed experiments in the Coriolis facility which is a 13-m diameter wave tank. We took care of cleaning surface contamination as well as possible considering that the surface of water exceeds 100~m$^2$. We observe that for the cleanest condition the frequency energy spectrum shows a power law decay extending up to the gravity capillary crossover (14 Hz) with a spectral exponent that is increasing with the forcing strength and decaying with surface contamination. Although slightly higher than reported previously in the literature, the exponent for the cleanest water remains significantly below the prediction from the Weak Turbulence Theory. By discussing length and time scales, we show that weak turbulence cannot be expected at frequencies above 3 Hz. We observe with a stereoscopic reconstruction technique that the increase with the forcing strength of energy spectrum beyond 3~Hz is mostly due to the formation and strenghtening of bound waves.


title: Order of tangency between manifolds

abstract: We study the order of tangency between two manifolds of same dimension and give that notion three quite different geometric interpretations. Related aspects of the order of tangency, e.g., regular separation exponents, are also discussed.


title: Cavallo's Multiplier for in situ Generation of High Voltage

abstract: A classic electrostatic induction machine, Cavallo's multiplier, is suggested for in situ production of very high voltage in cryogenic environments. The device is suitable for generating a large electrostatic field under conditions of very small load current. Operation of the Cavallo multiplier is analyzed, with quantitative description in terms of mutual capacitances between electrodes in the system. A demonstration apparatus was constructed, and measured voltages are compared to predictions based on measured capacitances in the system. The simplicity of the Cavallo multiplier makes it amenable to electrostatic analysis using finite element software, and electrode shapes can be optimized to take advantage of a high dielectric strength medium such as liquid helium. A design study is presented for a Cavallo multiplier in a large-scale, cryogenic experiment to measure the neutron electric dipole moment.


title: Exciton Control in a Room-Temperature Bulk Semiconductor with Coherent Strain Pulses

abstract: The coherent manipulation of excitons in bulk semiconductors via the lattice degrees of freedom is key to the development of acousto-optic and acousto-excitonic devices. Wide-bandgap transition metal oxides exhibit strongly bound excitons that are interesting for applications in the deep-ultraviolet, but their properties have remained elusive due to the lack of efficient generation and detection schemes in this spectral range. Here, we perform ultrafast broadband deep-ultraviolet spectroscopy on anatase TiO$_2$ single crystals at room temperature, and reveal a dramatic modulation of the exciton peak amplitude due to coherent acoustic phonons. This modulation is comparable to those of nanostructures where exciton-phonon coupling is enhanced by quantum confinement, and is accompanied by a giant exciton shift of 30-50 meV. We model these results by many-body perturbation theory and show that the deformation potential coupling within the nonlinear regime is the main mechanism for the generation and detection of the coherent acoustic phonons. Our findings pave the way to the design of exciton control schemes in the deep-ultraviolet with propagating strain pulses.


title: Edgeworth expansions for weakly dependent random variables

abstract: We discuss sufficient conditions that guarantee the existence of asymptotic expansions for the Central Limit Theorem for weakly dependent random variables including observations arising from sufficiently chaotic dynamical systems like piece-wise expanding maps, and strongly ergodic Markov chains. We primarily use spectral techniques to obtain the results.


title: Hybrid asymptotic/numerical methods for the evaluation of layer heat potentials in two dimensions

abstract: We present a hybrid asymptotic/numerical method for the accurate computation of single and double layer heat potentials in two dimensions. It has been shown in previous work that simple quadrature schemes suffer from a phenomenon called "geometrically-induced stiffness," meaning that formally high-order accurate methods require excessively small time steps before the rapid convergence rate is observed. This can be overcome by analytic integration in time, requiring the evaluation of a collection of spatial boundary integral operators with non-physical, weakly singular kernels. In our hybrid scheme, we combine a local asymptotic approximation with the evaluation of a few boundary integral operators involving only Gaussian kernels, which are easily accelerated by a new version of the fast Gauss transform. This new scheme is robust, avoids geometrically-induced stiffness, and is easy to use in the presence of moving geometries. Its extension to three dimensions is natural and straightforward, and should permit layer heat potentials to become flexible and powerful tools for modeling diffusion processes.


title: A note on the complexity of evolutionary dynamics in a classic consumer-resource model

abstract: We study how the complexity of evolutionary dynamics in the classic MacArthur consumer-resource model depends on resource uptake and utilization rates. The traditional assumption in such models is that the utilization rate of the consumer is proportional to the uptake rate. More generally, we show that if these two rates are related through a power law (which includes the traditional assumption as a special case), then the resulting evolutionary dynamics in the consumer is necessarily a simple hill-climbing process leading to an evolutionary equilibrium, regardless of the dimension of phenotype space. When utilization and uptake rates are not related by a power law, more complex evolutionary trajectories can occur, including the chaotic dynamics observed in previous studies for high-dimensional phenotype spaces. These results draw attention to the importance of distinguishing between utilization and uptake rates in consumer-resource models.


title: Comments on the double copy construction for gravitational theories

abstract: We revisit the double copy description for linearized gravity and point out various technical issues and subtleties, associated with setting up the double copy description, including the problem of matching degrees of freedom on both sides of the double copy dictionary and the related issue of the constraint between graviton and dilaton sources. We introduce and discuss possible resolutions of these issues.


title: Multi-Modal Geometric Learning for Grasping and Manipulation

abstract: This work provides an architecture that incorporates depth and tactile information to create rich and accurate 3D models useful for robotic manipulation tasks. This is accomplished through the use of a 3D convolutional neural network (CNN). Offline, the network is provided with both depth and tactile information and trained to predict the object's geometry, thus filling in regions of occlusion. At runtime, the network is provided a partial view of an object. Tactile information is acquired to augment the captured depth information. The network can then reason about the object's geometry by utilizing both the collected tactile and depth information. We demonstrate that even small amounts of additional tactile information can be incredibly helpful in reasoning about object geometry. This is particularly true when information from depth alone fails to produce an accurate geometric prediction. Our method is benchmarked against and outperforms other visual-tactile approaches to general geometric reasoning. We also provide experimental results comparing grasping success with our method.


title: Variability of Brown Dwarfs

abstract: Brown dwarfs constitute a missing link between low-mass stars and giant planets. Their atmospheres display chemical species typical of planets, and one could wonder whether they also have weather-like patterns. While brown dwarf surface features cannot be directly resolved, the photometric and spectroscopic modulations induced by these features, as they rotate in and out of view, provide a wealth of information on the evolution of their atmosphere. A review of brown dwarfs variability through the L, T and Y spectral types sequence is presented, as well as the constraints that they set on the nature of weather-like patterns on their surface.


title: Ultra-Low Latency (ULL) Networks: The IEEE TSN and IETF DetNet Standards and Related 5G ULL Research

abstract: Many network applications, e.g., industrial control, demand Ultra-Low Latency (ULL). However, traditional packet networks can only reduce the end-to-end latencies to the order of tens of milliseconds. The IEEE 802.1 Time Sensitive Networking (TSN) standard and related research studies have sought to provide link layer support for ULL networking, while the emerging IETF Deterministic Networking (DetNet) standards seek to provide the complementary network layer ULL support. This article provides an up-to-date comprehensive survey of the IEEE TSN and IETF DetNet standards and the related research studies. The survey of these standards and research studies is organized according to the main categories of flow concept, flow synchronization, flow management, flow control, and flow integrity. ULL networking mechanisms play a critical role in the emerging fifth generation (5G) network access chain from wireless devices via access, backhaul, and core networks. We survey the studies that specifically target the support of ULL in 5G networks, with the main categories of fronthaul, backhaul, and network management. Throughout, we identify the pitfalls and limitations of the existing standards and research studies. This survey can thus serve as a basis for the development of standards enhancements and future ULL research studies that address the identified pitfalls and limitations.


title: Joule overheating poisons the fractional ac Josephson effect in topological Josephson junctions

abstract: Topological Josephson junctions designed on the surface of a 3D-topological insulator (TI) harbor Majorana bound states (MBS's) among a continuum of conventional Andreev bound states. The distinct feature of these MBS's lies in the $4\pi$-periodicity of their energy-phase relation that yields a fractional ac Josephson effect and a suppression of odd Shapiro steps under $r\!f$ irradiation. Yet, recent experiments showed that a few, or only the first, odd Shapiro steps are missing, casting doubts on the interpretation. Here, we show that Josephson junctions tailored on the large bandgap 3D TI Bi$_2$Se$_3$ exhibit a fractional ac Josephson effect acting on the first Shapiro step only. With a modified resistively shunted junction model, we demonstrate that the resilience of higher order odd Shapiro steps can be accounted for by thermal poisoning driven by Joule overheating. Furthermore, we uncover a residual supercurrent at the nodes between Shapiro lobes, which provides a direct and novel signature of the current carried by the MBS. Our findings showcase the crucial role of thermal effects in topological Josephson junctions and lend support to the Majorana origin of the partial suppression of odd Shapiro steps.


title: HPL-GEM: Controlling High Pressure Laminates bulk resistivity with GEMs

abstract: We succeeded in modifying and controlling the electrical resistance of a standard High Pressure Laminate (HPL) panel through the use of a Gas Electron Multiplier (GEM) foil that has been embedded into the bulk of the HPL plate itself. Electrical characterizations were made and preliminary data show that this HPL-GEM embedded system can easily vary its bulk resistance by more than one order of magnitude. Data show that the bulk resistance change is exponential with the applied voltage to the embedded GEM.


title: Poisson geometry, monoidal Fukaya categories, and commutative Floer cohomology rings

abstract: We describe connections between concepts arising in Poisson geometry and the theory of Fukaya categories. The key concept is that of a symplectic groupoid, which is an integration of a Poisson manifold. The Fukaya category of a symplectic groupoid is monoidal, and it acts on the Fukaya categories of the symplectic leaves of the Poisson structure. Conversely, we consider a wide range of known monoidal structures on Fukaya categories and observe that they all arise from symplectic groupoids. We also use the picture developed to resolve a conundrum in Floer theory: why are some Lagrangian Floer cohomology rings commutative?


title: 8th Low-Level RF Workshop (LLRF2017)

abstract: This volume contains a subset of contributions presented at LLRF2017: the 8th Low-Level RF Workshop held in Barcelona, Spain, on October 16-19, 2017.


title: Lagrange's Theorem For Hom-Groups

abstract: Hom-groups are nonassociative generalizations of groups where the unitality and associativity are twisted by a map. We show that a Hom-group (G, {\alpha}) is a pointed idempotent quasigroup (pique). We use Cayley table of quasigroups to introduce some examples of Hom-groups. Introducing the notions of Hom-subgroups and cosets we prove Lagrange's theorem for finite Hom-groups. This states that the order of any Hom-subgroup H of a finite Hom-group G divides the order of G. We linearize Hom-groups to obtain a class of nonassociative Hopf algebras called Hom-Hopf algebras. As an application of our results, we show that the dimension of a Hom-sub-Hopf algebra of the finite dimensional Hom-group Hopf algebra KG divides the order of G. The new tools introduced in this paper could potentially have applications in theories of quasigroups, nonassociative Hopf algebras, Hom-type objects, combinatorics, and cryptography.


title: Product Characterisation towards Personalisation: Learning Attributes from Unstructured Data to Recommend Fashion Products

abstract: In this paper, we describe a solution to tackle a common set of challenges in e-commerce, which arise from the fact that new products are continually being added to the catalogue. The challenges involve properly personalising the customer experience, forecasting demand and planning the product range. We argue that the foundational piece to solve all of these problems is having consistent and detailed information about each product, information that is rarely available or consistent given the multitude of suppliers and types of products. We describe in detail the architecture and methodology implemented at ASOS, one of the world's largest fashion e-commerce retailers, to tackle this problem. We then show how this quantitative understanding of the products can be leveraged to improve recommendations in a hybrid recommender system approach.


title: PaaS Cloud: The Business Perspective

abstract: The next generation of PaaS technology accomplishes the true promise of object-oriented and 4GLs development with less effort. Now PaaS is becoming one of the core technical services for application development organizations. PaaS offers a resourceful and agile approach to develop, operate and deploy applications in a cost-effective manner. It is now turning out to be one of the preferred choices throughout the world, especially for globally distributed development environment. However it still lacks the scale of popularity and acceptance which Software-as-a-Service (SaaS) and Infrastructure-as-a-Service (IaaS) have attained. PaaS offers a promising future with novel technology architecture and evolutionary development approach. In this article, we identify the strengths, weaknesses, opportunities and threats for the PaaS industry. We then identify the various issues that will affect the different stakeholders of PaaS industry. This research will outline a set of recommendations for the PaaS practitioners to better manage this technology. For PaaS technology researchers, we also outline the number of research areas that need attention in coming future. Finally, we also included an online survey to outline PaaS technology market leaders. This will facilitate PaaS technology practitioners to have a more deep insight into market trends and technologies.


title: Pathwise approximation of Feynman path integrals using simple random walks

abstract: The aim of the presented research is to give a rigorous mathematical approach to Feynman path integrals based on strong (pathwise) approximations based on simple random walks.


title: A Feature-Driven Active Framework for Ultrasound-Based Brain Shift Compensation

abstract: A reliable Ultrasound (US)-to-US registration method to compensate for brain shift would substantially improve Image-Guided Neurological Surgery. Developing such a registration method is very challenging, due to factors such as missing correspondence in images, the complexity of brain pathology and the demand for fast computation. We propose a novel feature-driven active framework. Here, landmarks and their displacement are first estimated from a pair of US images using corresponding local image features. Subsequently, a Gaussian Process (GP) model is used to interpolate a dense deformation field from the sparse landmarks. Kernels of the GP are estimated by using variograms and a discrete grid search method. If necessary, the user can actively add new landmarks based on the image context and visualization of the uncertainty measure provided by the GP to further improve the result. We retrospectively demonstrate our registration framework as a robust and accurate brain shift compensation solution on clinical data acquired during neurosurgery.


title: On the Complexity of Testing Attainment of the Optimal Value in Nonlinear Optimization

abstract: We prove that unless P=NP, there exists no polynomial time (or even pseudo-polynomial time) algorithm that can test whether the optimal value of a nonlinear optimization problem where the objective and constraints are given by low-degree polynomials is attained. If the degrees of these polynomials are fixed, our results along with previously-known "Frank-Wolfe type" theorems imply that exactly one of two cases can occur: either the optimal value is attained on every instance, or it is strongly NP-hard to distinguish attainment from non-attainment. We also show that testing for some well-known sufficient conditions for attainment of the optimal value, such as coercivity of the objective function and closedness and boundedness of the feasible set, is strongly NP-hard. As a byproduct, our proofs imply that testing the Archimedean property of a quadratic module is strongly NP-hard, a property that is of independent interest to the convergence of the Lasserre hierarchy. Finally, we give semidefinite programming (SDP)-based sufficient conditions for attainment of the optimal value, in particular a new characterization of coercive polynomials that lends itself to an SDP hierarchy.


title: Characterization by forbidden induced subgraphs of some subclasses of chordal graphs

abstract: Chordal graphs are the graphs in which every cycle of length at least four has a chord. A set $S$ is a vertex separator for vertices $a$ and $b$ if the removal of $S$ of the graph separates $a$ and $b$ into distinct connected components. A graph $G$ is chordal if and only if every minimal vertex separator is a clique. We study subclasses of chordal graphs defined by restrictions imposed on the intersections of its minimal separator cliques. Our goal is to characterize them by forbidden induced subgraphs. Some of these classes have already been studied such as chordal graphs in which two minimal separators have no empty intersection if and only if they are equal. Those graphs are known as strictly chordal graphs and they were first introduced as block duplicate graphs by Golumbic and Peled. They were also considered in other previous works, showing that strictly chordal graphs are exactly the (gem, dart)-free graphs.


title: aKWISP: investigating short-distance interactions at sub-micron scales

abstract: The sub-micron range in the field of short distance interactions has yet to be opened to experimental investigation, and may well hold the key to understanding al least part of the dark matter puzzle. The aKWISP (advanced-KWISP) project introduces the novel Double Membrane Interaction Monitor (DMIM), a combined source-sensing device where interaction distances can be as short as 100 nm or even 10 nm, much below the 1-10 micron distance which is the lower limit encountered by current experimental efforts. aKWISP builds on the technology and the results obtained with the KWISP opto-mechanical force sensor now searching at CAST for the direct coupling to matter of solar chameleons. It will reach the ultimate quantum-limited sensitivity by exploiting an array of technologies, including operation at milli-Kelvin temperatures. Recent suggestions point at short-distance interactions studies as intriguing possibilities for the detection of axions and of new physical phenomena.


title: Big-bang nucleosynthesis and Leptogenesis in CMSSM

abstract: We have investigated the constrained minimal supersymmetric standard model with three right-handed Majorana neutrinos whether there still is a parameter region which is consistent with all existing experimental data/limits such as Leptogenesis and the dark matter abundance and we also can solve the Lithium problem. Using Casas-Ibarra parameterization, we have found that a very narrow parameter space of the complex orthogonal matrix elements where the lightest slepton can have a long lifetime, that is necessary for solving the Lithium problem. Further, under this condition, there is a parameter region that can give an explanation for the experimental observations. We have studied three cases of the right-handed neutrino mass ratio \mbox{\em (i)} $M_{2}=2 \times M_{1}$, \mbox{\em (ii)} $M_{2}=4 \times M_{1}$, \mbox{\em (iii)} $M_{2}=10 \times M_{1}$ while $M_{3}=40 \times M_{1}$ is fixed. We have obtained the mass range of the lightest right-handed neutrino mass that lies between $10^9$ GeV and $10^{11}$ GeV. The important result is that its upper limit is derived by solving the Lithium problem and the lower limit comes from Leptogenesis. Calculated low-energy observables of these parameter sets such as BR($\mu \to e \gamma$) is not yet restricted by experiments and will be verified in the near future.


title: Tidal deformability from GW170817 as a direct probe of the neutron star radius

abstract: Gravitational waves from the coalescence of two neutron stars were recently detected for the first time by the LIGO-Virgo collaboration, in event GW170817. This detection placed an upper limit on the effective tidal deformability of the two neutron stars and tightly constrained the chirp mass of the system. We report here on a new simplification that arises in the effective tidal deformability of the binary, when the chirp mass is specified. We find that, in this case, the effective tidal deformability of the binary is surprisingly independent of the component masses of the individual neutron stars, and instead depends primarily on the ratio of the chirp mass to the neutron star radius. Thus, a measurement of the effective tidal deformability can be used to directly measure the neutron star radius. We find that the upper limit on the effective tidal deformability from GW170817 implies that the radius cannot be larger than ~13km, at the 90% level, independent of the assumed masses for the component stars. The result can be applied generally, to probe the stellar radii in any neutron star-neutron star merger with a measured chirp mass. The approximate mass-independence disappears for neutron star-black hole mergers. Finally, we discuss a Bayesian inference of the equation of state that uses the measured chirp mass and tidal deformability from GW170817 combined with nuclear and astrophysical priors and discuss possible statistical biases in this inference.


title: Point process models for quasi-periodic volcanic earthquakes

abstract: Long period (LP) earthquakes are common at active volcanoes, and are ubiquitous at persistently active andesitic and dacitic subduction zone volcanoes. They provide critical information regarding the state of volcanic unrest, and their occurrence rates are key data for eruption forecasting. LPs are commonly quasi-periodic or 'anti-clustered', unlike volcano-tectonic (VT) earthquakes, so the existing Poisson point process methods used to model occurrence rates of VT earthquakes are unlikely to be optimal for LP data. We evaluate the performance of candidate formulations for LP data, based on inhomogeneous point process models with four different inter-event time distributions: exponential (IP), Gamma (IG), inverse Gaussian (IIG), and Weibull (IW). We examine how well these models explain the observed data, and the quality of retrospective forecasts of eruption time. We use a Bayesian MCMC approach to fit the models. Goodness-of-fit is assessed using Quantile-Quantile and Kolmogorov-Smirnov methods, and benchmarking against results obtained from synthetic datasets. IG and IIG models were both found to fit the data well, with the IIG model slightly outperforming the IG model. Retrospective forecasting analysis shows that the IG model performs best, with the initial preference for the IIG model controlled by catalogue incompleteness late in the sequence. The IG model fits the data significantly better than the IP model, and simulations show it produces better forecasts for highly periodic data. Simulations also show that forecast precision increases with the degree of periodicity of the earthquake process using the IG model, and so should be better for LP earthquakes than VTs. These results provide a new framework for point process modelling of volcanic earthquake time series, and verification of alternative models.


title: Join-Idle-Queue with Service Elasticity: Large-Scale Asymptotics of a Non-monotone System

abstract: We consider the model of a token-based joint auto-scaling and load balancing strategy, proposed in a recent paper by Mukherjee, Dhara, Borst, and van Leeuwaarden (SIGMETRICS '17, arXiv:1703.08373), which offers an efficient scalable implementation and yet achieves asymptotically optimal steady-state delay performance and energy consumption as the number of servers $N\to\infty$. In the above work, the asymptotic results are obtained under the assumption that the queues have fixed-size finite buffers, and therefore the fundamental question of stability of the proposed scheme with infinite buffers was left open. In this paper, we address this fundamental stability question. The system stability under the usual subcritical load assumption is not automatic. Moreover, the stability may not even hold for all $N$. The key challenge stems from the fact that the process lacks monotonicity, which has been the powerful primary tool for establishing stability in load balancing models. We develop a novel method to prove that the subcritically loaded system is stable for large enough $N$, and establish convergence of steady-state distributions to the optimal one, as $N \to \infty$. The method goes beyond the state of the art techniques -- it uses an induction-based idea and a "weak monotonicity" property of the model; this technique is of independent interest and may have broader applicability.


title: Adaptive Super-twisting Second-order Sliding Mode for Attitude Control of Quadcopter UAVs

abstract: This work addresses the modelling and control aspects for quadcopter or drone unmanned aerial vehicles (UAVs). First, the mathematical model of the drone is derived by identifying significant parameters and the negligible ones are treated as disturbances. The control design begins with the switching surface selection, then, an Adaptive Super Twisting Sliding Mode (ASTSM) Control algorithm is applied to adjust attitudes of the quadcopter under harsh conditions such as nonlinear, strong coupling, high uncertainties and disturbances. Simulation results show that the proposed controller can achieve robust operation with disturbance rejection, parametric variation adaptation as well as chattering attenuation. Comparisons with some commonly used and advanced controllers in a quadcopter model show advantages of the proposed control scheme.


title: Modeling evolution of dark matter substructure and annihilation boost

abstract: We study evolution of dark matter substructures, especially how they lose the mass and change density profile after they fall in gravitational potential of larger host halos. We develop an analytical prescription that models the subhalo mass evolution and calibrate it to results of N-body numerical simulations of various scales from very small (Earth size) to large (galaxies to clusters) halos. We then combine the results with halo accretion histories, and calculate the subhalo mass function that is physically motivated down to Earth-mass scales. Our results --- valid for arbitrary host masses and redshifts --- show reasonable agreement with those of numerical simulations at resolved scales. Our analytical model also enables self-consistent calculations of the boost factor of dark matter annhilation, which we find to increase from tens of percent at the smallest (Earth) and intermediate (dwarfs) masses to a factor of several at galaxy size, and to become as large as a factor of $\sim$10 for the largest halos (clusters) at small redshifts. Our analytical approach can accommodate substructures in the subhalos (sub-subhalos) in a consistent framework, which we find to give up to a factor of a few enhancement to the annihilation boost. Presence of the subhalos enhances the intensity of the isotropic gamma-ray background by a factor of a few, and as the result, the measurement by Fermi Large Area Telescope excludes the annihilation cross section greater than $\sim$$4\times 10^{-26}$ cm$^3$ s$^{-1}$ for dark matter masses up to $\sim$200 GeV.


title: The $HI$- and $H_{2}$-to-stellar mass correlations of late- and early-type galaxies and their consistency with the observational mass functions

abstract: We compile and carrefully homogenize local galaxy samples with available information on stellar, $\rm HI$ and/or $\rm H_{2}$ masses, and morphology. After processing the information on upper limits in the case of non gas detections, we determine the $\rm HI$- and $\rm H_{2}$-to-stellar mass relations and their $1\sigma$ scatter for both late- and early-type galaxies. The obtained relations are fitted to single or double power laws. Late-type galaxies are significantly gas richer than early-type ones, specially at high masses. The respective $\rm H_{2}$-to-$\rm HI$ mass ratios as a function of $M_{\ast}$ are discussed. Further, we constrain the full mass-dependent distribution functions of the $\rm HI$- and $\rm H_{2}$-to-stellar mass ratios. We find that they can be described by a Schechter function for late types and a (broken) Schechter + uniform function for early types. By using the observed galaxy stellar mass function and the volume-complete late-to-early-type galaxy ratio as a function of $M_{\ast}$, these empirical distribution functions are mapped into $\rm HI$ and $\rm H_{2}$ mass functions. The obtained mass functions are consistent with those inferred from large surveys. The empirical gas-to-stellar mass relations and their distributions for local late- and early-type galaxies presented here can be used to constrain models and simulations of galaxy evolution.


title: Balanced Black and White Coloring Problem on knights chessboards

abstract: Graph anticoloring problem is partial coloring problem where the main feature is the opposite rule of the graph coloring problem, i.e., if two vertices are adjacent, their assigned colors must be the same or at least one of them is uncolored. In the same way, Berge in 1972 proposed the problem of placing b black queens and w white queens on a $n \times n$ chessboard such that no two queens of different color can attack to each other, the complexity of this problem remains open. In this work we deal with the knight piece under the balance property, since this special case is the most difficult for brute force algorithms.


title: Defective and Clustered Graph Colouring

abstract: Consider the following two ways to colour the vertices of a graph where the requirement that adjacent vertices get distinct colours is relaxed. A colouring has "defect" $d$ if each monochromatic component has maximum degree at most $d$. A colouring has "clustering" $c$ if each monochromatic component has at most $c$ vertices. This paper surveys research on these types of colourings, where the first priority is to minimise the number of colours, with small defect or small clustering as a secondary goal. List colouring variants are also considered. The following graph classes are studied: outerplanar graphs, planar graphs, graphs embeddable in surfaces, graphs with given maximum degree, graphs with given maximum average degree, graphs excluding a given subgraph, graphs with linear crossing number, linklessly or knotlessly embeddable graphs, graphs with given Colin de Verdi\`ere parameter, graphs with given circumference, graphs excluding a fixed graph as an immersion, graphs with given thickness, graphs with given stack- or queue-number, graphs excluding $K_t$ as a minor, graphs excluding $K_{s,t}$ as a minor, and graphs excluding an arbitrary graph $H$ as a minor. Several open problems are discussed.


title: Detection and characterization of spin-orbit resonances in the advanced gravitational wave detectors era

abstract: In this paper, we test the performance of templates in detection and characterization of Spin-orbit resonant (SOR) binaries. We use precessing SEOBNRv3 waveforms as well as {\it four} numerical relativity (NR) waveforms to model GWs from SOR binaries and filter them through IMRPhenomD, SEOBNRv4 (non-precessing) and IMRPhenomPv2 (precessing) approximants. We find that IMRPhenomD and SEOBNRv4 recover only $\sim70\%$ of injections with fitting factor (FF) higher than 0.97 (or 90\% of injections with ${\rm FF} >0.9$).However, using the sky-maxed statistic, IMRPhenomPv2 performs magnificently better than their non-precessing counterparts with recovering $99\%$ of the injections with FFs higher than 0.97. Interestingly, injections with $\Delta \phi = 180^{\circ}$ have higher FFs ($\Delta \phi$ is the angle between the components of the black hole spins in the plane orthogonal to the orbital angular momentum) as compared to their $\Delta \phi =0^{\circ}$ and generic counterparts. This implies that we will have a slight observation bias towards $\Delta \phi=180^{\circ}$ SORs while using non-precessing templates for searches. All template approximants are able to recover most of the injected NR waveforms with FFs $>0.95$. For all the injections including NR, the error in estimating chirp mass remains below $<10\%$ with minimum error for $\Delta \phi = 180^{\circ}$ resonant binaries. The symmetric mass ratio can be estimated with errors below $15\%$. The effective spin parameter $\chi_{\rm eff}$ is measured with maximum absolute error of 0.13. The in-plane spin parameter $\chi_p$ is mostly underestimated indicating that a precessing signal will be recovered as a relatively less precessing signal. Based on our findings, we conclude that we not only need improvements in waveform models towards precession and non-quadrupole modes but also better search strategies for precessing GW signals.


title: Inverse Optimal Control from Incomplete Trajectory Observations

abstract: This article develops a methodology that enables learning an objective function of an optimal control system from incomplete trajectory observations. The objective function is assumed to be a weighted sum of features (or basis functions) with unknown weights, and the observed data is a segment of a trajectory of system states and inputs. The proposed technique introduces the concept of the recovery matrix to establish the relationship between any available segment of the trajectory and the weights of given candidate features. The rank of the recovery matrix indicates whether a subset of relevant features can be found among the candidate features and the corresponding weights can be learned from the segment data. The recovery matrix can be obtained iteratively and its rank non-decreasing property shows that additional observations may contribute to the objective learning. Based on the recovery matrix, a method for using incomplete trajectory observations to learn the weights of selected features is established, and an incremental inverse optimal control algorithm is developed by automatically finding the minimal required observation. The effectiveness of the proposed method is demonstrated on a linear quadratic regulator system and a simulated robot manipulator.


title: 180-degree phase shift of magnetoelastic waves observed by phase-resolved spin-wave tomography

abstract: We have investigated optically-excited magnetoelastic waves by phase-resolved spin-wave tomography (PSWaT). PSWaT reconstructs dispersion relation of spin waves together with their phase information by using time-resolved magneto-optical imaging for spin-wave propagation followed by an analysis based on the convolution theorem and a complex Fourier transform. In PSWaT spectra for a Bi-doped garnet film, we found a 180 degree phase shift of magnetoelastic waves at around the crossing of the dispersion relations of spin and elastic waves. The result is explained by a coupling between spin waves and elastic waves through magnetoelastic interaction. We also propose an efficient way for phase manipulation of magnetoelastic waves by rotating the orientation of magnetization less than 10 degree.


title: The classification of $n$-dimensional anticommutative algebras with $(n-3)$-dimensional annihilator

abstract: We give the classification of all $n$-dimensional anticommutative complex algebras with $(n-3)$-dimensional annihilator. Namely, we describe all central extensions of all $3$-dimensional anticommutative complex algebras.


title: On the role of secondary motions in turbulent square duct flow

abstract: We use a direct numerical simulations (DNS) database for turbulent flow in a square duct up to bulk Reynolds number $\Rey_b=40000$, to quantitatively analyze the role of secondary motions on the mean flow structure. For that purpose we derive a generalized form of the identity of Fukagata, Iwamoto and Kasagi (FIK), which allows to quantify the effect of cross-stream convection on the mean streamwise velocity, wall shear stress and bulk friction coefficient. Secondary motions are found to contribute for about $6\%$ of total friction, and to act as a self-regulating mechanism of turbulence whereby wall shear stress nonuniformities induced by corners are equalized, and universality of the wall-normal velocity profiles is established. We also carry out numerical experiments whereby the secondary motions are artificially suppressed, in which case their equalizing role is partially taken by the turbulent stresses.


title: Instability of the solitary wave solutions for the genenalized derivative Nonlinear Schr\"odinger equation in the critical frequency case

abstract: We study the stability theory of solitary wave solutions for the generalized derivative nonlinear Schr\"odinger equation $$ i\partial_{t}u+\partial_{x}^{2}u+i|u|^{2\sigma}\partial_x u=0. $$ The equation has a two-parameter family of solitary wave solutions of the form \begin{align*} \phi_{\omega,c}(x)=\varphi_{\omega,c}(x)\exp{\big\{ i\frac c2 x-\frac{i}{2\sigma+2}\int_{-\infty}^{x}\varphi^{2\sigma}_{\omega,c}(y)dy\big\}}. \end{align*} Here $ \varphi_{\omega,c}$ is some real-valued function. It was proved in \cite{LiSiSu1} that the solitary wave solutions are stable if $-2\sqrt{\omega }<c <2z_0\sqrt{\omega }$, and unstable if $2z_0\sqrt{\omega }<c <2\sqrt{\omega }$ for some $z_0\in(0,1)$. We prove the instability at the borderline case $c =2z_0\sqrt{\omega }$ for $1<\sigma<2$, improving the previous results in \cite{Fu-16-DNLS} where $3/2<\sigma<2$.


title: Very high efficiency of low cost graphite-based solar cell by improving the fill factor using optimal ion concentration in polymer electrolyte

abstract: We report the development of graphite-based solar cells using a simple method and low cost materials. Suspension of graphite powder in mineral water was simply dropped onto the surface of fluorine-doped tin oxide glass (FTO) to form a thick film. Surprisingly, using mineral waters greatly improved the efficiency of the solar cell to reach the highest efficiency of 6.97%. Due to some minerals contained, the mineral water induced the development of fibrous structure between the graphite particles which is assumed to play a role as a bridge for the photoexcited electrons to quickly move to the electrode and suppress recombination with holes. This efficiency is very attractive when considering the materials used to develop the solar cell are all low cost. Economically this may challenge the present high efficiency semiconductor-based solar cells. We achieved the high efficiency by manipulating the cell fill factor through optimizing the ion concentration in PVA.LiOH polymer electrolyte. We also propose an equation to describe the effect of LiOH concentration and efficiency and we also provide strong correlation between the cell efficiency and the polymer conductivity


title: Robust Depth Estimation from Auto Bracketed Images

abstract: As demand for advanced photographic applications on hand-held devices grows, these electronics require the capture of high quality depth. However, under low-light conditions, most devices still suffer from low imaging quality and inaccurate depth acquisition. To address the problem, we present a robust depth estimation method from a short burst shot with varied intensity (i.e., Auto Bracketing) or strong noise (i.e., High ISO). We introduce a geometric transformation between flow and depth tailored for burst images, enabling our learning-based multi-view stereo matching to be performed effectively. We then describe our depth estimation pipeline that incorporates the geometric transformation into our residual-flow network. It allows our framework to produce an accurate depth map even with a bracketed image sequence. We demonstrate that our method outperforms state-of-the-art methods for various datasets captured by a smartphone and a DSLR camera. Moreover, we show that the estimated depth is applicable for image quality enhancement and photographic editing.


title: Weakly Supervised Medical Diagnosis and Localization from Multiple Resolutions

abstract: Diagnostic imaging often requires the simultaneous identification of a multitude of findings of varied size and appearance. Beyond global indication of said findings, the prediction and display of localization information improves trust in and understanding of results when augmenting clinical workflow. Medical training data rarely includes more than global image-level labels as segmentations are time-consuming and expensive to collect. We introduce an approach to managing these practical constraints by applying a novel architecture which learns at multiple resolutions while generating saliency maps with weak supervision. Further, we parameterize the Log-Sum-Exp pooling function with a learnable lower-bounded adaptation (LSE-LBA) to build in a sharpness prior and better handle localizing abnormalities of different sizes using only image-level labels. Applying this approach to interpreting chest x-rays, we set the state of the art on 9 abnormalities in the NIH's CXR14 dataset while generating saliency maps with the highest resolution to date.


title: On irrationality of hypersurfaces In $\mathbf{P}^{n+1}$

abstract: In this note, we study various measures of irrationality for hypersurfaces in projective spaces which were recently proposed by Bastianelli, De Poi, Ein, Lazarsfeld and Ullery. In particular, we answer the question raised by Bastianelli that if $X \subset P^{n+1}$ is a very general smooth hypersurface of dimension $n$ and degree $d\geq 2n+2$, then $\text{stab.irr}(X)=\text{uni.irr}(X)=d-1$. As a corollary, we prove that $\text{irr}(X\times P^{m})=\text{irr}(X)$ for any integer $m\geq 1$.


title: Markov Chains with Maximum Return Time Entropy for Robotic Surveillance

abstract: Motivated by robotic surveillance applications, this paper studies the novel problem of maximizing the return time entropy of a Markov chain, subject to a graph topology with travel times and stationary distribution. The return time entropy is the weighted average, over all graph nodes, of the entropy of the first return times of the Markov chain; this objective function is a function series that does not admit in general a closed form. The paper features theoretical and computational contributions. First, we obtain a discrete-time delayed linear system for the return time probability distribution and establish its convergence properties. We show that the objective function is continuous over a compact set and therefore admits a global maximum; a unique globally-optimal solution is known only for complete graphs with unitary travel times. We then establish upper and lower bounds between the return time entropy and the well-known entropy rate of the Markov chain. To compute the optimal Markov chain numerically, we establish the asymptotic equality between entropy, conditional entropy and truncated entropy, and propose an iteration to compute the gradient of the truncated entropy. Finally, we apply these results to the robotic surveillance problem. Our numerical results show that, for a model of rational intruder over prototypical graph topologies and test cases, the maximum return time entropy chain performs better than several existing Markov chains.


title: On the Prony Series Representation of Stretched Exponential Relaxation

abstract: Stretched exponential relaxation is a ubiquitous feature of homogeneous glasses. The stretched exponential decay function can be derived from the diffusion-trap model, which predicts certain critical values of the fractional stretching exponent. In practical implementations of glass relaxation models, it is computationally convenient to represent the stretched exponential function as a Prony series of simple exponentials. Here, we perform a comprehensive mathematical analysis of the Prony series approximation of the stretched exponential relaxation, including optimized coefficients for certain critical values of the exponent. The fitting quality of the Prony series is analyzed as a function of the number of terms in the series. With a sufficient number of terms, the Prony series can accurately capture the time evolution of the stretched exponential function, including its "fat tail" at long times. However, it is unable to capture the divergence of the first-derivative of the stretched exponential function in the limit of zero time. We also present a frequency-domain analysis of the Prony series representation of the stretched exponential function and discuss its physical implications for the modeling of glass relaxation behavior.


title: Pressure-Induced Phase Transformation in $\beta$-Eucryptite: an X-Ray Diffraction and Density Functional Theory Study

abstract: Certain alumino-silicates display exotic properties enabled by their framework structure made of corner-sharing tetrahedral rigid units. Using \textit{in situ} diamond-anvil cell x-ray diffraction (XRD), we study the pressure-induced transformation of $\beta$ eucryptite, a prototypical alumino-silicate. $\beta$ eucryptite undergoes a phase transformation at moderate pressures, but the atomic structure of the new phase has not yet been reported. Based on density functional theory stability studies and Rietveld analysis of XRD patterns, we find that the pressure-stabilized phase belongs to the Pna2$_1$ space group. Furthermore, we discover two other possible pressure-stabilized polymorphs, P1c1 and Pca2$_1$.


title: K2 Ultracool Dwarfs Survey. III. White Light Flares are Ubiquitous in M6-L0 Dwarfs

abstract: We report the white light flare rates for 10 ultracool dwarfs (UCDs) using \textit{Kepler K2} short cadence data. Among our sample stars, two have spectral type M6, three are M7, three are M8 and two are L0. Most of our targets are old low mass stars. We identify a total of 283 flares in all of the stars in our sample, with \textit{Kepler} energies in the range log \textit{E$_{Kp}$} $\sim$(29 - 33.5) erg. Using the maximum-likelihood method of line fitting, we find that the flare frequency distribution (FFD) for each star in our sample follows a power law with slope -$\alpha$ in range -(1.3-2.0). We find that cooler objects tend to have shallower slopes. For some of our targets, the FFD follows either a broken power law, or a power law with an exponential cutoff. For the L0 dwarf 2MASS J12321827-0951502, we find a very shallow slope (-$\alpha$ $=$ -1.3) in the \textit{Kepler} energy range (0.82-130)$\times$10$^{30}$ erg: this L0 dwarf has flare rates which are comparable to the rates of high energy flares in stars of earlier spectral types. In addition, we report photometry of two superflares: one on the L0 dwarf 2MASS J12321827-0951502 and another on the M7 dwarf 2MASS J08352366+1029318. In case of 2MASS J12321827-0951502, we report a flare brightening by a factor of $\sim$144 relative to the quiescent photospheric level. Likewise, for 2MASS J08352366+1029318, we report a flare brightening by a factor of $\sim$60 relative to the quiescent photospheric level. These two superflares have bolometric (UV/optical/infrared) energies 3.6 $\times$ 10$^{33}$ erg and 8.9 $\times$ 10$^{33}$ erg respectively, while the FWHM time scales are very short, $\sim$2 minutes. We find that the M8 star TRAPPIST-1 is more active than the M8.5 dwarf: 2M03264453+1919309, but less active than another M8 dwarf (2M12215066-0843197).


title: Time dilation in relativistic quantum decay laws of moving unstable particles

abstract: The relativistic quantum decay laws of moving unstable particles are analyzed for a general class of mass distribution densities which behave as power laws near the (non-vanishing) lower bound $\mu_0$ of the mass spectrum. The survival probability $\mathcal{P}_p(t)$, the instantaneous mass $M_p(t)$ and the instantaneous decay rate $\Gamma_p(t)$ of the moving unstable particle are evaluated over short and long times for an arbitrary value $p$ of the (constant) linear momentum. The ultrarelativistic and non-relativistic limits are studied. Over long times, the survival probability $\mathcal{P}_p(t)$ is approximately related to the survival probability at rest $\mathcal{P}_0(t)$ by a scaling law. The scaling law can be interpreted as the effect of the relativistic time dilation if the asymptotic value $M_p\left(\infty\right)$ of the instantaneous mass is considered as the effective mass of the unstable particle over long times. The effective mass has magnitude $\mu_0$ at rest and moves with linear momentum $p$ or, equivalently, with constant velocity $1\Big/\sqrt{1+\mu_0^2\big/p^2}$. The instantaneous decay rate $\Gamma_p(t)$ is approximately independent of the linear momentum $p$, over long times, and, consequently, is approximately invariant by changing reference frame.


title: Inference in Probabilistic Graphical Models by Graph Neural Networks

abstract: A fundamental computation for statistical inference and accurate decision-making is to compute the marginal probabilities or most probable states of task-relevant variables. Probabilistic graphical models can efficiently represent the structure of such complex data, but performing these inferences is generally difficult. Message-passing algorithms, such as belief propagation, are a natural way to disseminate evidence amongst correlated variables while exploiting the graph structure, but these algorithms can struggle when the conditional dependency graphs contain loops. Here we use Graph Neural Networks (GNNs) to learn a message-passing algorithm that solves these inference tasks. We first show that the architecture of GNNs is well-matched to inference tasks. We then demonstrate the efficacy of this inference approach by training GNNs on a collection of graphical models and showing that they substantially outperform belief propagation on loopy graphs. Our message-passing algorithms generalize out of the training set to larger graphs and graphs with different structure.


title: Data-Driven Computational Methods: Parameter and Operator Estimations (Chapter 1)

abstract: Modern scientific computational methods are undergoing a transformative change; big data and statistical learning methods now have the potential to outperform the classical first-principles modeling paradigm. This book bridges this transition, connecting the theory of probability, stochastic processes, functional analysis, numerical analysis, and differential geometry. It describes two classes of computational methods to leverage data for modeling dynamical systems. The first is concerned with data fitting algorithms to estimate parameters in parametric models that are postulated on the basis of physical or dynamical laws. The second class is on operator estimation, which uses the data to nonparametrically approximate the operator generated by the transition function of the underlying dynamical systems. This self-contained book is suitable for graduate studies in applied mathematics, statistics, and engineering. Carefully chosen elementary examples with supplementary MATLAB codes and appendices covering the relevant prerequisite materials are provided, making it suitable for self-study.


title: Causal Inference on Discrete Data via Estimating Distance Correlations

abstract: In this paper, we deal with the problem of inferring causal directions when the data is on discrete domain. By considering the distribution of the cause $P(X)$ and the conditional distribution mapping cause to effect $P(Y|X)$ as independent random variables, we propose to infer the causal direction via comparing the distance correlation between $P(X)$ and $P(Y|X)$ with the distance correlation between $P(Y)$ and $P(X|Y)$. We infer "$X$ causes $Y$" if the dependence coefficient between $P(X)$ and $P(Y|X)$ is smaller. Experiments are performed to show the performance of the proposed method.


title: Robust Beamforming for SWIPT System with Chance Constraints

abstract: The robust beamforming problem in multiple-input single-output (MISO) downlink networks of simultaneous wireless information and power transfer (SWIPT) is studied in this paper. Adopting the time switching fashion to perform energy harvesting and information decoding respectively, we aim at maximizing the sum rate under imperfect channel state information (CSI) and the chance constraints of users' harvested energy. In view of the fact that the constraints for minimal harvested energy is not necessary to meet from time to time, this paper adopts chance constraint to model it and uses the Bernstein inequality to transform it into deterministic constraints equivalently. Recognizing the maximum sum rate problem of imperfect CSI as nonconvex problem, we transform it into finding the expectation of minimum mean square error (MMSE) equivalently in this paper, and an alternative optimization (AO) algorithm is proposed to decompose the optimization problem into two sub-problems: the transmit beamformer design and the division of switching time. The simulation results show the performance gains compared to non-robust state of the art schemes.


title: Breathing mode frequency of a strongly interacting Fermi gas across the 2D-3D dimensional crossover

abstract: We address the interplay between dimension and quantum anomaly on the breathing mode frequency of a strongly interacting Fermi gas harmonically trapped at zero temperature. Using a beyond mean-field, Gaussian pair fluctuation theory, we employ periodic boundary conditions to simulate the dimensionality of the system and impose a local density approximation, with two different schemes, to model different trapping potentials in the tightly-confined axial direction. By using a sum-rule approach, we compute the breathing mode frequency associated with a small variation of the trapping frequency along the weakly-confined transverse direction, and describe its behavior as functions of the dimensionality, from two- to three-dimensions, and of the interaction strength. We compare our predictions with previous calculations on the two-dimensional breathing mode anomaly and discuss their possible observation in ultracold Fermi gases of $^{6}$Li and $^{40}$K atoms.


title: SurvBoost: An R Package for High-Dimensional Variable Selection in the Stratified Proportional Hazards Model via Gradient Boosting

abstract: High-dimensional variable selection in the proportional hazards (PH) model has many successful applications in different areas. In practice, data may involve confounding variables that do not satisfy the PH assumption, in which case the stratified proportional hazards (SPH) model can be adopted to control the confounding effects by stratification of the confounding variable, without directly modeling the confounding effects. However, there is lack of computationally efficient statistical software for high-dimensional variable selection in the SPH model. In this work, an R package, SurvBoost, is developed to implement the gradient boosting algorithm for fitting the SPH model with high-dimensional covariate variables and other confounders. Extensive simulation studies demonstrate that in many scenarios SurvBoost can achieve a better selection accuracy and reduce computational time substantially compared to the existing R package that implements boosting algorithms without stratification. The proposed R package is also illustrated by an analysis of the gene expression data with survival outcome in The Cancer Genome Atlas (TCGA) study. In addition, a detailed hands-on tutorial for SurvBoost is provided.


title: Generative Adversarial Talking Head: Bringing Portraits to Life with a Weakly Supervised Neural Network

abstract: This paper presents Generative Adversarial Talking Head (GATH), a novel deep generative neural network that enables fully automatic facial expression synthesis of an arbitrary portrait with continuous action unit (AU) coefficients. Specifically, our model directly manipulates image pixels to make the unseen subject in the still photo express various emotions controlled by values of facial AU coefficients, while maintaining her personal characteristics, such as facial geometry, skin color and hair style, as well as the original surrounding background. In contrast to prior work, GATH is purely data-driven and it requires neither a statistical face model nor image processing tricks to enact facial deformations. Additionally, our model is trained from unpaired data, where the input image, with its auxiliary identity label taken from abundance of still photos in the wild, and the target frame are from different persons. In order to effectively learn such model, we propose a novel weakly supervised adversarial learning framework that consists of a generator, a discriminator, a classifier and an action unit estimator. Our work gives rise to template-and-target-free expression editing, where still faces can be effortlessly animated with arbitrary AU coefficients provided by the user.


title: The Ice Cap Zone: A Unique Habitable Zone for Ocean Worlds

abstract: Traditional definitions of the habitable zone assume that habitable planets contain a carbonate-silicate cycle that regulates CO2 between the atmosphere, surface, and the interior. Such theories have been used to cast doubt on the habitability of ocean worlds. However, Levi et al (2017) have recently proposed a mechanism by which CO2 is mobilized between the atmosphere and the interior of an ocean world. At high enough CO2 pressures, sea ice can become enriched in CO2 clathrates and sink after a threshold density is achieved. The presence of subpolar sea ice is of great importance for habitability in ocean worlds. It may moderate the climate and is fundamental in current theories of life formation in diluted environments. Here, we model the Levi et al. mechanism and use latitudinally-dependent non-grey energy balance and single-column radiative-convective climate models and find that this mechanism may be sustained on ocean worlds that rotate at least 3 times faster than the Earth. We calculate the circumstellar region in which this cycle may operate for G-M-stars (Teff = 2,600 to 5,800 K), extending from about 1.23 to 1.65, 0.69 to 0.954, 0.38 to 0.528 AU, 0.219 to 0.308 AU, 0.146 to 0.206 AU, and 0.0428 to 0.0617 AU for G2, K2, M0, M3, M5, and M8 stars, respectively. However, unless planets are very young and not tidally locked, our mechanism would be unlikely to apply to stars cooler than a ~M3. We predict C/O ratios for our atmospheres (about 0.5) that can be verified by the JWST mission.


title: InfyNLP at SMM4H Task 2: Stacked Ensemble of Shallow Convolutional Neural Networks for Identifying Personal Medication Intake from Twitter

abstract: This paper describes Infosys's participation in the "2nd Social Media Mining for Health Applications Shared Task at AMIA, 2017, Task 2". Mining social media messages for health and drug related information has received significant interest in pharmacovigilance research. This task targets at developing automated classification models for identifying tweets containing descriptions of personal intake of medicines. Towards this objective we train a stacked ensemble of shallow convolutional neural network (CNN) models on an annotated dataset provided by the organizers. We use random search for tuning the hyper-parameters of the CNN and submit an ensemble of best models for the prediction task. Our system secured first place among 9 teams, with a micro-averaged F-score of 0.693.


title: A conservative energy-momentum tensor in the $f(R,T)$ gravity and its implications for the phenomenology of neutron stars

abstract: The solutions for the Tolmann-Oppenheimer-Volkoff (TOV) equation bring valuable informations about the macroscopical features of compact astrophysical objects as neutron stars. They are sensitive to both the equation of state considered for nuclear matter and the background gravitational theory. In this work we construct the TOV equation for a conservative version of the $f(R,T)$ gravity. While the non-vanishing of the covariant derivative of the $f(R,T)$ energy-momentum tensor yields, in a cosmological perspective, the prediction of creation of matter throughout the universe evolution as shown by T. Harko, in the analysis of the hydrostatic equilibrium of compact astrophysical objects, this property still lacks a convincing physical explanation. The imposition of $\nabla^{\mu}T_{\mu\nu}=0$ demands a particular form for the function $h(T)$ in $f(R,T)=R+h(T)$, which is here derived. Therefore, the choice of a specific equation of state for the star matter demands a unique form of $h(T)$, manifesting a strong connection between conserved $f(R,T)$ gravity and the star matter constitution. We construct and solve the TOV equation for the general equation of state for $p=k\rho^{\Gamma}$, with $k$ being the EoS parameter, $\rho$ {\it the energy density} and $\Gamma$ is the adiabatic index. We also derive the macroscopical properties of neutron stars ($\Gamma=5/3$) within this approach.


title: Asymptotic Optimal Portfolio in Fast Mean-reverting Stochastic Environments

abstract: This paper studies the portfolio optimization problem when the investor's utility is general and the return and volatility of the risky asset are fast mean-reverting, which are important to capture the fast-time scale in the modeling of stock price volatility. Motivated by the heuristic derivation in [J.-P. Fouque, R. Sircar and T. Zariphopoulou, \emph{Mathematical Finance}, 2016], we propose a zeroth order strategy, and show its asymptotic optimality within a specific (smaller) family of admissible strategies under proper assumptions. This optimality result is achieved by establishing a first order approximation of the problem value associated to this proposed strategy using singular perturbation method, and estimating the risk-tolerance functions. The results are natural extensions of our previous work on portfolio optimization in a slowly varying stochastic environment [J.-P. Fouque and R. Hu, \emph{SIAM Journal on Control and Optimization}, 2017], and together they form a whole picture of analyzing portfolio optimization in both fast and slow environments.


title: Modeling Camera Effects to Improve Visual Learning from Synthetic Data

abstract: Recent work has focused on generating synthetic imagery to increase the size and variability of training data for learning visual tasks in urban scenes. This includes increasing the occurrence of occlusions or varying environmental and weather effects. However, few have addressed modeling variation in the sensor domain. Sensor effects can degrade real images, limiting generalizability of network performance on visual tasks trained on synthetic data and tested in real environments. This paper proposes an efficient, automatic, physically-based augmentation pipeline to vary sensor effects --chromatic aberration, blur, exposure, noise, and color cast-- for synthetic imagery. In particular, this paper illustrates that augmenting synthetic training datasets with the proposed pipeline reduces the domain gap between synthetic and real domains for the task of object detection in urban driving scenes.


title: A Robust Fault-Tolerant and Scalable Cluster-wide Deduplication for Shared-Nothing Storage Systems

abstract: Deduplication has been largely employed in distributed storage systems to improve space efficiency. Traditional deduplication research ignores the design specifications of shared-nothing distributed storage systems such as no central metadata bottleneck, scalability, and storage rebalancing. Further, deduplication introduces transactional changes, which are prone to errors in the event of a system failure, resulting in inconsistencies in data and deduplication metadata. In this paper, we propose a robust, fault-tolerant and scalable cluster-wide deduplication that can eliminate duplicate copies across the cluster. We design a distributed deduplication metadata shard which guarantees performance scalability while preserving the design constraints of shared- nothing storage systems. The placement of chunks and deduplication metadata is made cluster-wide based on the content fingerprint of chunks. To ensure transactional consistency and garbage identification, we employ a flag-based asynchronous consistency mechanism. We implement the proposed deduplication on Ceph. The evaluation shows high disk-space savings with minimal performance degradation as well as high robustness in the event of sudden server failure.


title: Poisson sigma model and semiclassical quantization of integrable systems

abstract: In this paper we outline the construction of semiclassical eigenfunctions of integrable models in terms of the semiclassical path integral for the Poisson sigma model with the target space being the phase space of the integrable system. The semiclassical path integral is defined as a formal power series with coefficients being Feynman diagrams. We also argue that in a similar way one can obtain irreducible semiclassical representations of Kontsevich's star product.


title: Attention on Attention: Architectures for Visual Question Answering (VQA)

abstract: Visual Question Answering (VQA) is an increasingly popular topic in deep learning research, requiring coordination of natural language processing and computer vision modules into a single architecture. We build upon the model which placed first in the VQA Challenge by developing thirteen new attention mechanisms and introducing a simplified classifier. We performed 300 GPU hours of extensive hyperparameter and architecture searches and were able to achieve an evaluation score of 64.78%, outperforming the existing state-of-the-art single model's validation score of 63.15%.


title: Semidefinite Outer Approximation of the Backward Reachable Set of Discrete-time Autonomous Polynomial Systems

abstract: We approximate the backward reachable set of discrete-time autonomous polynomial systems using the recently developed occupation measure approach. We formulate the problem as an infinite-dimensional linear programming (LP) problem on measures and its dual on continuous functions. Then we approximate the LP by a hierarchy of finite-dimensional semidefinite programming (SDP) programs on moments of measures and their duals on sums-of-squares polynomials. Finally we solve the SDP's and obtain a sequence of outer approximations of the backward reachable set. We demonstrate our approach on three dynamical systems. As a special case, we also show how to approximate the preimage of a compact semi-algebraic set under a polynomial map.


title: Gradient Descent with Random Initialization: Fast Global Convergence for Nonconvex Phase Retrieval

abstract: This paper considers the problem of solving systems of quadratic equations, namely, recovering an object of interest $\mathbf{x}^{\natural}\in\mathbb{R}^{n}$ from $m$ quadratic equations/samples $y_{i}=(\mathbf{a}_{i}^{\top}\mathbf{x}^{\natural})^{2}$, $1\leq i\leq m$. This problem, also dubbed as phase retrieval, spans multiple domains including physical sciences and machine learning. We investigate the efficiency of gradient descent (or Wirtinger flow) designed for the nonconvex least squares problem. We prove that under Gaussian designs, gradient descent --- when randomly initialized --- yields an $\epsilon$-accurate solution in $O\big(\log n+\log(1/\epsilon)\big)$ iterations given nearly minimal samples, thus achieving near-optimal computational and sample complexities at once. This provides the first global convergence guarantee concerning vanilla gradient descent for phase retrieval, without the need of (i) carefully-designed initialization, (ii) sample splitting, or (iii) sophisticated saddle-point escaping schemes. All of these are achieved by exploiting the statistical models in analyzing optimization algorithms, via a leave-one-out approach that enables the decoupling of certain statistical dependency between the gradient descent iterates and the data.

categories: stat.ML cs.IT cs.LG cs.NA math.IT math.OC


title: A family of Bell transformations

abstract: We introduce a family of sequence transformations, defined via partial Bell polynomials, that may be used for a systematic study of a wide variety of problems in enumerative combinatorics. This family includes some of the transformations listed in the paper by Bernstein & Sloane, now seen as transformations under the umbrella of partial Bell polynomials. Our goal is to describe these transformations from the algebraic and combinatorial points of view. We provide functional equations satisfied by the generating functions, derive inverse relations, and give a convolution formula. While the full range of applications remains unexplored, in this paper we show a glimpse of the versatility of Bell transformations by discussing the enumeration of several combinatorial configurations, including rational Dyck paths, rooted planar maps, and certain classes of permutations.


title: Unsupervised Representation Learning by Predicting Image Rotations

abstract: Over the last years, deep convolutional neural networks (ConvNets) have transformed the field of computer vision thanks to their unparalleled capacity to learn high level semantic image features. However, in order to successfully learn those features, they usually require massive amounts of manually labeled data, which is both expensive and impractical to scale. Therefore, unsupervised semantic feature learning, i.e., learning without requiring manual annotation effort, is of crucial importance in order to successfully harvest the vast amount of visual data that are available today. In our work we propose to learn image features by training ConvNets to recognize the 2d rotation that is applied to the image that it gets as input. We demonstrate both qualitatively and quantitatively that this apparently simple task actually provides a very powerful supervisory signal for semantic feature learning. We exhaustively evaluate our method in various unsupervised feature learning benchmarks and we exhibit in all of them state-of-the-art performance. Specifically, our results on those benchmarks demonstrate dramatic improvements w.r.t. prior state-of-the-art approaches in unsupervised representation learning and thus significantly close the gap with supervised feature learning. For instance, in PASCAL VOC 2007 detection task our unsupervised pre-trained AlexNet model achieves the state-of-the-art (among unsupervised methods) mAP of 54.4% that is only 2.4 points lower from the supervised case. We get similarly striking results when we transfer our unsupervised learned features on various other tasks, such as ImageNet classification, PASCAL classification, PASCAL segmentation, and CIFAR-10 classification. The code and models of our paper will be published on: https://github.com/gidariss/FeatureLearningRotNet .


title: Look Before You Leap: Bridging Model-Free and Model-Based Reinforcement Learning for Planned-Ahead Vision-and-Language Navigation

abstract: Existing research studies on vision and language grounding for robot navigation focus on improving model-free deep reinforcement learning (DRL) models in synthetic environments. However, model-free DRL models do not consider the dynamics in the real-world environments, and they often fail to generalize to new scenes. In this paper, we take a radical approach to bridge the gap between synthetic studies and real-world practices---We propose a novel, planned-ahead hybrid reinforcement learning model that combines model-free and model-based reinforcement learning to solve a real-world vision-language navigation task. Our look-ahead module tightly integrates a look-ahead policy model with an environment model that predicts the next state and the reward. Experimental results suggest that our proposed method significantly outperforms the baselines and achieves the best on the real-world Room-to-Room dataset. Moreover, our scalable method is more generalizable when transferring to unseen environments.


title: The Origins Space Telescope: Towards An Understanding of Temperate Planetary Atmospheres

abstract: The Origins Space Telescope (OST) is one of four mission concepts currently being studied by NASA in preparation for the Astrophysics 2020 Decadal Survey. With active cooling (~4 K), OST will be sensitive in mid- to far-IR wavelengths, using imaging and spectroscopy to probe the furthest reaches of our galaxies, trace the path of water through star and planet formation, and place thermochemical constraints on the atmospheres of exoplanets ranging in size from Jupiter to Earth. This contribution to the Exoplanet Science Strategy committee discusses the significant advancements that the OST Mid-Infrared Imager, Spectrometer, and Coronagraph (MISC) instrument can make in studying cool planetary atmospheres. We particularly focus on the atmospheres of transiting rocky planets in the habitable zones of mid-to-late M stars. We discuss how OST thermal infrared observations can significantly enhance our understanding of the temperature structure and molecular abundances of biologically interesting gases on these worlds, including O3, CH4, H2O, and CO2.


title: Mixed-timescale Per-group Hybrid Precoding for Multiuser Massive MIMO Systems

abstract: Considering the expensive radio frequency (RF) chain, huge training overhead and feedback burden issues in massive MIMO, in this letter, we propose a mixed-timescale per-group hybrid precoding (MPHP) scheme under an adaptive partially-connected RF precoding structure (PRPS), where the RF precoder is implemented using an adaptive connection network (ACN) and M analog phase shifters (APSs), where M is the number of antennas at the base station (BS). Exploiting the mixed-time stage channel state information (CSI) structure, the joint-design of ACN and APSs is formulated as a statistical signal-to-leakage-and-noise ratio (SSLNR) maximization problem, and a heuristic group RF precoding (GRFP) algorithm is proposed to provide a near-optimal solution. Simulation results show that the proposed design advances at better energy efficiency (EE) and lower hardware cost, CSI signaling overhead and computational complexity than the conventional hybrid precoding (HP) schemes.


title: Superconvergence Points of Integer and Fractional Derivatives of Special Hermite Interpolations and Its Applications in Solving FDEs

abstract: In this paper, we study convergence and superconvergence theory of integer and fractional derivatives of the one-point and the two-point Hermite interpolations. When considering the integer-order derivative, exponential decay of the error is proved, and superconvergence points are located, at which the convergence rates are $O(N^{-2})$ and $O(N^{-1.5})$, respectively, better than the global rate for the one-point and two-point interpolations. Here $N$ represents the degree of interpolation polynomial. It is proved that the $\alpha$-th fractional derivative of $(u-u_N)$ with $k<\alpha<k+1$, is bounded by its $(k+1)$-th derivative. Furthermore, the corresponding superconvergence points are predicted for fractional derivatives, and an eigenvalue method is proposed to calculate the superconvergence points for the Riemann-Liouville fractional derivative. In the application of the knowledge of superconvergence points to solve FDEs, we discover that a modified collocation method makes numerical solutions much more accurate than the traditional collocation method.


title: Bach Flow on Homogeneous Products

abstract: Qualitative behavior of Bach flow is established on compact four-dimensional locally homogeneous product manifolds. This is achieved by lifting to the homogeneous universal cover and, in most cases, capitalizing on the resultant group structure. The resulting system of ordinary differential equations is carefully analyzed on a case-by-case basis, with explicit solutions found in some cases. Limiting behavior of the metric and the curvature are determined in all cases. The behavior on quotients of $\mathbb{R} \times \mathbb{S}^3$ proves to be the most challenging and interesting.


title: Adaptive Sequential MCMC for Combined State and Parameter Estimation

abstract: In the case of a linear state space model, we implement an MCMC sampler with two phases. In the learning phase, a self-tuning sampler is used to learn the parameter mean and covariance structure. In the estimation phase, the parameter mean and covariance structure informs the proposed mechanism and is also used in a delayed-acceptance algorithm. Information on the resulting state of the system is given by a Gaussian mixture. In on-line mode, the algorithm is adaptive and uses a sliding window approach to accelerate sampling speed and to maintain appropriate acceptance rates. We apply the algorithm to joined state and parameter estimation in the case of irregularly sampled GPS time series data.


title: Phase amplification in optical interferometry with weak measurement

abstract: Improving the phase resolution of interferometry is crucial for high-precision measurements of various physical quantities. Systematic phase errors dominate the phase uncertainties in most realistic optical interferometers. Here we propose and experimentally demonstrate a weak measurement scheme to considerably suppress the phase uncertainties by the direct amplification of phase shift in optical interferometry. Given an initial ultra-small phase shift between orthogonal polarization states, we observe the phase amplification effect with a factor of 388. Our weak measurement scheme provides a practical approach to significantly improve the interferometric phase resolution, which is favorable for precision measurement applications.


title: Asteroseismic Investigations of the Binary System HD176465

abstract: HD 176465 is a binary system whose both components are solar-like pulsators and whose oscillation frequencies were observed by the \kepler\ mission. In this paper we have modeled the asteroseismic and spectroscopic data of the stars, and have determined their convection-zone helium abundances using the signatures left by the He{\sc ii} ionization zone on the mode frequencies. As expected we find that the components of the binary have the same age within uncertainties ($3.087 \pm 0.580$ Gyr and $3.569 \pm 0.912$ Gyr); they also have the same initial helium abundance (Y$_{\mathrm{init}}$=0.253 $\pm$ 0.006 and 0.254 $\pm$ 0.008). Their current metallicity ([Fe/H]=$-0.275 \pm 0.04$ and $-0.285 \pm 0.04$) is also the same within errors. Fits to the signature of the He{\sc ii} acoustic glitch yields current helium abundances of $Y_{\rm A} = 0.224 \pm 0.006$ and $Y_{\rm B} = 0.233 \pm 0.008$ for the two components. Analyzing the complete ensemble of models generated for this investigation we find that both the amplitude and acoustic depth of the glitch signature arising from the second helium ionization zone and the base of the convection zone (CZ) are functions of mass. We show that the acoustic depths of these glitches are positively correlated with each other. The analysis can help us to detect the internal structure and constrain the chemical compositions.


title: PyramidBox: A Context-assisted Single Shot Face Detector

abstract: Face detection has been well studied for many years and one of remaining challenges is to detect small, blurred and partially occluded faces in uncontrolled environment. This paper proposes a novel context-assisted single shot face detector, named \emph{PyramidBox} to handle the hard face detection problem. Observing the importance of the context, we improve the utilization of contextual information in the following three aspects. First, we design a novel context anchor to supervise high-level contextual feature learning by a semi-supervised method, which we call it PyramidAnchors. Second, we propose the Low-level Feature Pyramid Network to combine adequate high-level context semantic feature and Low-level facial feature together, which also allows the PyramidBox to predict faces of all scales in a single shot. Third, we introduce a context-sensitive structure to increase the capacity of prediction network to improve the final accuracy of output. In addition, we use the method of Data-anchor-sampling to augment the training samples across different scales, which increases the diversity of training data for smaller faces. By exploiting the value of context, PyramidBox achieves superior performance among the state-of-the-art over the two common face detection benchmarks, FDDB and WIDER FACE. Our code is available in PaddlePaddle: \href{https://github.com/PaddlePaddle/models/tree/develop/fluid/face_detection}{\url{https://github.com/PaddlePaddle/models/tree/develop/fluid/face_detection}}.


title: Speech Emotion Recognition Considering Local Dynamic Features

abstract: Recently, increasing attention has been directed to the study of the speech emotion recognition, in which global acoustic features of an utterance are mostly used to eliminate the content differences. However, the expression of speech emotion is a dynamic process, which is reflected through dynamic durations, energies, and some other prosodic information when one speaks. In this paper, a novel local dynamic pitch probability distribution feature, which is obtained by drawing the histogram, is proposed to improve the accuracy of speech emotion recognition. Compared with most of the previous works using global features, the proposed method takes advantage of the local dynamic information conveyed by the emotional speech. Several experiments on Berlin Database of Emotional Speech are conducted to verify the effectiveness of the proposed method. The experimental results demonstrate that the local dynamic information obtained with the proposed method is more effective for speech emotion recognition than the traditional global features.


title: Assessing Shape Bias Property of Convolutional Neural Networks

abstract: It is known that humans display "shape bias" when classifying new items, i.e., they prefer to categorize objects based on their shape rather than color. Convolutional Neural Networks (CNNs) are also designed to take into account the spatial structure of image data. In fact, experiments on image datasets, consisting of triples of a probe image, a shape-match and a color-match, have shown that one-shot learning models display shape bias as well. In this paper, we examine the shape bias property of CNNs. In order to conduct large scale experiments, we propose using the model accuracy on images with reversed brightness as a metric to evaluate the shape bias property. Such images, called negative images, contain objects that have the same shape as original images, but with different colors. Through extensive systematic experiments, we investigate the role of different factors, such as training data, model architecture, initialization and regularization techniques, on the shape bias property of CNNs. We show that it is possible to design different CNNs that achieve similar accuracy on original images, but perform significantly different on negative images, suggesting that CNNs do not intrinsically display shape bias. We then show that CNNs are able to learn and generalize the structures, when the model is properly initialized or data is properly augmented, and if batch normalization is used.


title: New measurements of excitation functions of 186W(p,x) nuclear reactions up to 65 MeV. Production of 178W/178mTa generator

abstract: New experimental excitation functions for proton induced reactions on natW are presented in the 32- 65 MeV energy range. The cross sections for natW(p,xn)186,184,183,182m,182g,181Re, naW(p,x)178W, natW(p,x)183,182,180m,177,176,175Ta, 175Hf, 177Lu were measured via an activation method by using a stacked-foil irradiation technique and high resolution gamma-ray spectrometry. The results were compared with predicted values obtained with the nuclear reaction model code TALYS (results taken from TENDL-2014 and TENDL-2015 on-line library). Production routes of the medically relevant radionuclides 186Re, 178W 178Ta and 181W are discussed.


title: A Distributed Stochastic Gradient Tracking Method

abstract: In this paper, we study the problem of distributed multi-agent optimization over a network, where each agent possesses a local cost function that is smooth and strongly convex. The global objective is to find a common solution that minimizes the average of all cost functions. Assuming agents only have access to unbiased estimates of the gradients of their local cost functions, we consider a distributed stochastic gradient tracking method. We show that, in expectation, the iterates generated by each agent are attracted to a neighborhood of the optimal solution, where they accumulate exponentially fast (under a constant step size choice). More importantly, the limiting (expected) error bounds on the distance of the iterates from the optimal solution decrease with the network size, which is a comparable performance to a centralized stochastic gradient algorithm. Numerical examples further demonstrate the effectiveness of the method.


title: Fast Semantic Segmentation on Video Using Block Motion-Based Feature Interpolation

abstract: Convolutional networks optimized for accuracy on challenging, dense prediction tasks are prohibitively slow to run on each frame in a video. The spatial similarity of nearby video frames, however, suggests opportunity to reuse computation. Existing work has explored basic feature reuse and feature warping based on optical flow, but has encountered limits to the speedup attainable with these techniques. In this paper, we present a new, two part approach to accelerating inference on video. First, we propose a fast feature propagation technique that utilizes the block motion vectors present in compressed video (e.g. H.264 codecs) to cheaply propagate features from frame to frame. Second, we develop a novel feature estimation scheme, termed feature interpolation, that fuses features propagated from enclosing keyframes to render accurate feature estimates, even at sparse keyframe frequencies. We evaluate our system on the Cityscapes and CamVid datasets, comparing to both a frame-by-frame baseline and related work. We find that we are able to substantially accelerate segmentation on video, achieving near real-time frame rates (20.1 frames per second) on large images (960 x 720 pixels), while maintaining competitive accuracy. This represents an improvement of almost 6x over the single-frame baseline and 2.5x over the fastest prior work.


title: A Supplementary Condition for the Convergence of the Control Policy during Adaptive Dynamic Programming

abstract: Reinforcement learning based adaptive/approximate dynamic programming (ADP) is a powerful technique to determine an approximate optimal controller for a dynamical system. These methods bypass the need to analytically solve the nonlinear Hamilton-Jacobi-Bellman equation, whose solution is often to difficult to determine but is needed to determine the optimal control policy. ADP methods usually employ a policy iteration algorithm that evaluates and improves a value function at every step to find the optimal control policy. Previous works in ADP have been lacking a stronger condition that ensures the convergence of the policy iteration algorithm. This paper provides a sufficient but not necessary condition that guarantees the convergence of an ADP algorithm. This condition may provide a more solid theoretical framework for ADP-based control algorithm design for nonlinear dynamical systems.


title: Passivity and Evolutionary Game Dynamics

abstract: This paper investigates an energy conservation and dissipation -- passivity -- aspect of dynamic models in evolutionary game theory. We define a notion of passivity using the state-space representation of the models, and we devise systematic methods to examine passivity and to identify properties of passive dynamic models. Based on the methods, we describe how passivity is connected to stability in population games and illustrate stability of passive dynamic models using numerical simulations.