title: On the Self-Similarity of Natural Stochastic Textures

abstract: Self-similarity is the essence of fractal images and, as such, characterizes natural stochastic textures. This paper is concerned with the property of self-similarity in the statistical sense in the case of fully-textured images that contain both stochastic texture and structural (mostly deterministic) information. We firstly decompose a textured image into two layers corresponding to its texture and structure, and show that the layer representing the stochastic texture is characterized by random phase of uniform distribution, unlike the phase of the structured information which is coherent. The uniform distribution of the the random phase is verified by using a suitable hypothesis testing framework. We proceed by proposing two approaches to assessment of self-similarity. The first is based on patch-wise calculation of the mutual information, while the second measures the mutual information that exists across scales. Quantifying the extent of self-similarity by means of mutual information is of paramount importance in the analysis of natural stochastic textures that are encountered in medical imaging, geology, agriculture and in computer vision algorithms that are designed for application on fully-textures images.


title: Single crystals of the layered dichalcogenides MoS2 and WS2 grown by liquid phase transport

abstract: The growth of single crystals of MoS2 and WS2 by materials transport through a liquid salt flux made from a low melting mixture of NaCl and CsCl is described. The crystals are characterized by single crystal X-ray diffraction, which reveals that they are 2H-MoS2 with a small percentage (about 3%) of 3R intergrowths. The 2H-WS2 crystals display less than 1% of 3R intergrowths. Photoluminescence spectra of exfoliated monolayers suggest that MoS2 grown by this method has superior crystallinity in comparison to commercially-available MoS2 crystals.


title: Revisiting Thin Silicon for Photovoltaics: A Technoeconomic Perspective

abstract: Adopting thin Si wafers for PV reduces capital expenditure (capex) and manufacturing cost, and accelerates the growth of PV manufacturing. There are two key questions about thin Si today: (a) how much can we still benefit economically from thinning wafers? (b) what are the technological challenges to transition to thin wafers? In this work, we re-evaluate the benefits and challenges of thin Si for current and future PV modules using a comprehensive techno-economic framework that couples device simulation, bottom-up cost modeling, and a cash-flow growth model. When adopting an advanced technology concept that features sufficiently good surface passivation, similarly high efficiencies are achievable for 50-um wafers as for 160-um ones. We then quantify the economic benefits for thin Si wafers in terms of poly-Si-to-module manufacturing capex, module cost, and levelized cost of electricity (LCOE) for utility PV systems. Particularly, LCOE favors thinner wafers for all investigated device architectures, and can potentially be reduced by more than 5% from the value of 160-um wafers. With further improvements in module efficiency, an advanced thin-wafer device concept with 50-um wafers could reduce manufacturing capex by 48%, module cost by 28%, and LCOE by 24%. Furthermore, we apply a sustainable growth model to investigate PV deployment scenarios in 2030. It is found that the state-of-the-art industry concept could not achieve the climate targets even with very aggressive financial scenarios, therefore the capex reduction benefit of thin wafers is needed to facilitate more rapid PV growth. Lastly, we discuss the remaining technological challenges and areas for innovation to enable high-yield manufacturing of high-efficiency PV modules with thin Si wafers.


title: 3-Lie bialgebras and 3-pre-Lie algebras induced by involutive derivations

abstract: In this paper, we study the structure of 3-Lie algebras with involutive derivations. We prove that if $A$ is an $m$-dimensional 3-Lie algebra with an involutive derivation $D$, then there exists a compatible 3-pre-Lie algebra $(A, \{ , , , \}_D)$ such that $A$ is the sub-adjacent 3-Lie algebra, and there is a local cocycle $3$-Lie bialgebraic structure on the $2m$-dimensional semi-direct product 3-Lie algebra $A\ltimes_{ad^*} A^*$, which is associated to the adjoint representation $(A, ad)$. By means of involutive derivations, the skew-symmetric solution of the 3-Lie classical Yang-Baxter equation in the 3-Lie algebra $A\ltimes_{ad^*}A^*$, a class of 3-pre-Lie algebras, and eight and ten dimensional local cocycle 3-Lie bialgebras are constructed.


title: An intriguing hyperelliptic Shimura curve quotient of genus 16

abstract: Let $F$ be the maximal totally real subfield of $\mathbf{Q}(\zeta_{32})$, the cyclotomic field of $32$nd roots of unity. Let $D$ be the quaternion algebra over $F$ ramified exactly at the unique prime above $2$ and 7 of the real places of $F$. Let $\mathcal{O}$ be a maximal order in $D$, and $X_0^D(1)$ the Shimura curve attached to $\mathcal{O}$. Let $C = X_0^D(1)/\langle w_D \rangle$, where $w_D$ is the unique Atkin-Lehner involution on $X_0^D(1)$. We show that the curve $C$ has several striking features. First, it is a hyperelliptic curve of genus $16$, whose hyperelliptic involution is exceptional. Second, there are $34$ Weierstrass points on $C$, and exactly half of these points are CM points; they are defined over the Hilbert class field of the unique CM extension $E/F$ of class number $17$ contained in $\mathbf{Q}(\zeta_{64})$, the cyclotomic field of $64$th roots of unity. Third, the normal closure of the field of $2$-torsion of the Jacobian of $C$ is the Harbater field $N$, the unique Galois number field $N/\mathbf{Q}$ unramified outside $2$ and $\infty$, with Galois group $\mathrm{Gal}(N/\mathbf{Q})\simeq F_{17} = \mathbf{Z}/17\mathbf{Z} \rtimes (\mathbf{Z}/17\mathbf{Z})^\times$. In fact, the Jacobian $\mathrm{Jac}(X_0^D(1))$ has the remarkable property that each of its simple factors has a $2$-torsion field whose normal closure is the field $N$. Finally, and perhaps the most striking fact about $C$, is that it is also hyperelliptic over $\mathbf{Q}$.


title: Heegaard Floer homology and cosmetic surgeries in $S^3$

abstract: If a knot $K$ in $S^3$ admits a pair of truly cosmetic surgeries, we show that the surgery slopes are either $\pm 2$ or $\pm 1/q$ for some value of $q$ that is explicitly determined by the knot Floer homology of $K$. Moreover, in the former case the genus of $K$ must be two, and in the latter case there is bound relating $q$ to the genus and the Heegaard Floer thickness of $K$. As a consequence, we show that the cosmetic crossing conjecture holds for alternating knots (or more generally, Heegaard Floer thin knots) with genus not equal to two. We also show that the conjecture holds for any knot $K$ for which each prime summand of $K$ has at most 16 crossings; our techniques rule out cosmetic surgeries in this setting except for slopes $\pm 1$ and $\pm 2$ on a small number of knots, and these remaining examples can be checked by comparing hyperbolic invariants. These results make use of the surgery formula for Heegaard Floer homology, which has already proved to be a powerful tool for obstructing cosmetic surgeries; we get stronger obstructions than previously known by considering the full graded theory. We make use of a new graphical interpretation of knot Floer homology and the surgery formula in terms of immersed curves, which makes the grading information we need easier to access.


title: Gateway Placement and Selection Solutions in WMNs: A Survey

abstract: Due to the high demand of Internet access by users, and the tremendous success of wireless technologies, Wireless Mesh Networks (WMNs) have become a promising solution. IGW Placement and Selection (GPS) are significantly investigated problems to achieve QoS requirements, network performance, and reduce deployment cost in WMNs. Best effort is made to classify different works in the literature based on network characteristics. At first, one of the most principal capabilities of WMNs, which is taking advantage of using multi-radio routers in a multi-channel network, is studied. In this article, GPS protocols considering a definition of three types of WMN are investigated based on channel-radio association including Single Radio Single Channel (SRSC), Single Radio Multi- Channel (SRMC), and Multi-Radio Multi-Channel (MRMC) WMNs. Furthermore, a classification regarding static and dynamic channel allocation policies is derived. In addition, the reported works from the viewpoint of network solutions are classified. The first perspective is: centralized, distributed or hybrid architectures. Following this classification, the studies are categorized regarding optimization techniques, which are operation research-based solutions, heuristic algorithms, and meta-heuristic-based algorithms.


title: Usability Squared: Principles for doing good systems research in robotics

abstract: Despite recent major advances in robotics research, massive injections of capital into robotics startups, and significant market appetite for robotic solutions, large-scale real-world deployments of robotic systems remain relatively scarce outside of heavy industry and (recently) warehouse logistics. In this paper, we posit that this scarcity comes from the difficulty of building even merely functional, first-pass robotic applications without a dizzying breadth and depth of expertise, in contrast to the relative ease with which non-experts in cloud computing can build complex distributed applications that function reasonably well. We trace this difficulty in application building to the paucity of good systems research in robotics, and lay out a path toward enabling application building by centering usability in systems research in two different ways: privileging the usability of the abstractions defined in systems research, and ensuring that the research itself is usable by application developers in the context of evaluating it for its applicability to their target domain by following principles of realism, empiricism, and exhaustive explication. In addition, we make some suggestions for community-level changes, incentives, and initiatives to create a better environment for systems work in robotics.


title: Global Convergence of Least Squares EM for Demixing Two Log-Concave Densities

abstract: This work studies the location estimation problem for a mixture of two rotation invariant log-concave densities. We demonstrate that Least Squares EM, a variant of the EM algorithm, converges to the true location parameter from a randomly initialized point. We establish the explicit convergence rates and sample complexity bounds, revealing their dependence on the signal-to-noise ratio and the tail property of the log-concave distribution. Moreover, we show that this global convergence property is robust under model mis-specification. Our analysis generalizes previous techniques for proving the convergence results for Gaussian mixtures. In particular, we make use of an angle-decreasing property for establishing global convergence of Least Squares EM beyond Gaussian settings, as $\ell_2$ distance contraction no longer holds globally for general log-concave mixtures.


title: Sparsity Invariance for Convex Design of Distributed Controllers

abstract: We address the problem of designing optimal linear time-invariant (LTI) sparse controllers for LTI systems, which corresponds to minimizing a norm of the closed-loop system subject to sparsity constraints on the controller structure. This problem is NP-hard in general and motivates the development of tractable approximations. We characterize a class of convex restrictions based on a new notion of Sparsity Invariance (SI). The underlying idea of SI is to design sparsity patterns for transfer matrices Y(s) and X(s) such that any corresponding controller K(s)=Y(s)X(s)^-1 exhibits the desired sparsity pattern. For sparsity constraints, the approach of SI goes beyond the notion of Quadratic Invariance (QI): 1) the SI approach always yields a convex restriction; 2) the solution via the SI approach is guaranteed to be globally optimal when QI holds and performs at least as well as considering a nearest QI subset. Moreover, the notion of SI naturally applies to designing structured static controllers, while QI is not utilizable. Numerical examples show that even for non-QI cases, SI can recover solutions that are 1) globally optimal and 2) strictly more performing than previous methods.


title: Brain Controllability: not a slam dunk yet

abstract: In our recent article (Tu et al., Warnings and caveats in brain controllability, arXiv:1705.08261) we provided quantitative evidence to show that there are warnings and caveats in the way Gu and collaborators (Gu et al. Controllability of structural brain networks. Nature communications 6 (2015): 8414) define brain controllability. The comment by Pasqualetti et al. (Pasqualetti et al. RE: Warnings and Caveats in Brain Controllability. NeuroImage 297 (2019), 586-588) confirms the need to go beyond the methodology and approach presented in Gu et al. original work. In fact, they recognize that the source of confusion is due to the fact that assessing controllability via numerical analysis typically leads to ill-conditioned problems, and thus often generates results that are difficult to interpret. This is indeed the first warning we discussed: our work was not meant to prove that brain networks are not controllable from one node, rather we wished to highlight that the one node controllability framework and all consequent results were not properly justified based on the methodology presented in Gu et al. We used in our work the same method of Gu et al. not because we believe it is the best methodology, but because we extensively investigated it with the aim of replicating, testing and extending their results. And the warning and caveats we have proposed are the results of this investigation.


title: Fast Magnetic Reconnection: Secondary Tearing Instability and Role of the Hall Term

abstract: Magnetic reconnection provides the primary source for explosive energy release, plasma heating and particle acceleration in many astrophysical environments. The last years witnessed a revival of interest in the MHD tearing instability as a driver for efficient reconnection. It has been established that, provided the current sheet aspect ratio becomes small enough ($a/L \sim S^{-1/3}$ for a given Lundquist number $S\gg 1$), reconnection occurs on ideal Alfv\'en timescales and becomes independent on $S$. Here we investigate, by means of two-dimensional simulations, the \emph{ideal} tearing instability in the Hall-MHD regime, which is appropriate when the width of the resistive layer $\delta$ becomes comparable to the ion inertial length $d_i$. Moreover, we study in detail the spontaneous development and reconnection of secondary current sheets, which for high $S$ naturally adjust to the ideal aspect ratio and hence their evolution proceeds very rapidly. For moderate low $S$, the aspect ratio tends to the Sweet-Parker scaling ($a/L \sim S^{-1/2}$), in order to fulfill the condition $\delta \ll a$ necessary for the onset of a tearing instability. When the Hall term is included, the reconnection rate of this secondary nonlinear phase is enhanced and, depending on the ratio $d_i/\delta$, can be twice with respect to the pure MHD case, and up to ten times larger than the linear phase. Therefore, the evolution of the tearing instability in thin current sheets in the Hall-MHD regime naturally leads to an explosive disruption of the reconnecting site and to energy release on super-Alfv\'enic timescales, as required to explain astrophysical observations.


title: The Origin of Parity Changes in the Solar Cycle

abstract: Although sunspots have been systematically observed on the Sun's surface over the last four centuries, their magnetic properties have been revealed and documented only since the early 1900s. Sunspots typically appear in pairs of opposite magnetic polarities which have a systematic orientation. This polarity orientation is opposite across the equator -- a trend that has persisted over the last century over which magnetic field observations exist. Taken together with the configuration of the global poloidal field of the Sun -- that governs the heliospheric open flux and cosmic ray flux at Earth -- this phenomenon is consistent with the dipolar parity state of an underlying magnetohydrodynamic dynamo. Although hemispheric asymmetry in the emergence of sunspots is observed in the Sun, a parity shift has never been observed. We simulate hemispheric asymmetry through the introduction of random fluctuations in a computational dynamo model of the solar cycle and demonstrate that changes in parity are indeed possible in long-term simulations covering thousands of years. Quadrupolar modes are found to exist over a significant fraction of the simulated time. In particular, we find that a parity shift in the underlying nature of the sunspot cycle is more likely to occur when sunspot activity dominates in any one hemisphere for a time which is significantly longer than the cycle period. We establish causal pathways connecting hemispheric asymmetry and cross-equatorial phase-shifts to parity flips in the underlying dynamo mechanism. Our findings indicate that the solar cycle may have resided in quadrupolar parity states in the distant past, and provides a possible pathway for predicting parity flips in the future.


title: Characterizing the Exact Behaviors of Temporal Difference Learning Algorithms Using Markov Jump Linear System Theory

abstract: In this paper, we provide a unified analysis of temporal difference learning algorithms with linear function approximators by exploiting their connections to Markov jump linear systems (MJLS). We tailor the MJLS theory developed in the control community to characterize the exact behaviors of the first and second order moments of a large family of temporal difference learning algorithms. For both the IID and Markov noise cases, we show that the evolution of some augmented versions of the mean and covariance matrix of the TD estimation error exactly follows the trajectory of a deterministic linear time-invariant (LTI) dynamical system. Applying the well-known LTI system theory, we obtain closed-form expressions for the mean and covariance matrix of the TD estimation error at any time step. We provide a tight matrix spectral radius condition to guarantee the convergence of the covariance matrix of the TD estimation error, and perform a perturbation analysis to characterize the dependence of the TD behaviors on learning rate. For the IID case, we provide an exact formula characterizing how the mean and covariance matrix of the TD estimation error converge to the steady state values. For the Markov case, we use our formulas to explain how the behaviors of TD learning algorithms are affected by learning rate and the underlying Markov chain. For both cases, upper and lower bounds for the mean square TD error are provided. The mean square TD error is shown to converge linearly to an exact limit.


title: Meta-learning Pseudo-differential Operators with Deep Neural Networks

abstract: This paper introduces a meta-learning approach for parameterized pseudo-differential operators with deep neural networks. With the help of the nonstandard wavelet form, the pseudo-differential operators can be approximated in a compressed form with a collection of vectors. The nonlinear map from the parameter to this collection of vectors and the wavelet transform are learned together from a small number of matrix-vector multiplications of the pseudo-differential operator. Numerical results for Green's functions of elliptic partial differential equations and the radiative transfer equations demonstrate the efficiency and accuracy of the proposed approach.


title: Hedetniemi's conjecture is asymptotically false

abstract: Extending a recent breakthrough of Shitov, we prove that the chromatic number of the tensor product of two graphs can be a constant factor smaller than the minimum chromatic number of the two graphs. More precisely, we prove that there exists an absolute constant $\delta>0$ such that for all $c$ sufficiently large, there exist graphs $G$ and $H$ with chromatic number at least $(1+\delta)c$ for which $\chi(G \times H) \le c$.


title: Interpolated Adversarial Training: Achieving Robust Neural Networks without Sacrificing Too Much Accuracy

abstract: Adversarial robustness has become a central goal in deep learning, both in the theory and the practice. However, successful methods to improve the adversarial robustness (such as adversarial training) greatly hurt generalization performance on the unperturbed data. This could have a major impact on how the adversarial robustness affects real world systems (i.e. many may opt to forego robustness if it can improve accuracy on the unperturbed data). We propose Interpolated Adversarial Training, which employs recently proposed interpolation based training methods in the framework of adversarial training. On CIFAR-10, adversarial training increases the standard test error (when there is no adversary) from 4.43% to 12.32%, whereas with our Interpolated adversarial training we retain the adversarial robustness while achieving a standard test error of only 6.45%. With our technique, the relative increase in the standard error for the robust model is reduced from 178.1% to just 45.5%. Moreover, we provide mathematical analysis of Interpolated Adversarial Training to confirm its efficiencies and demonstrate its advantages in terms of robustness and generalization.


title: A Low-rank Solver for the Stochastic Unsteady Navier-Stokes Problem

abstract: We study a low-rank iterative solver for the unsteady Navier-Stokes equations for incompressible flows with a stochastic viscosity. The equations are discretized using the stochastic Galerkin method, and we consider an all-at-once formulation where the algebraic systems at all the time steps are collected and solved simultaneously. The problem is linearized with Picard's method. To efficiently solve the linear systems at each step, we use low-rank tensor representations within the Krylov subspace method, which leads to significant reductions in storage requirements and computational costs. Combined with effective mean-based preconditioners and the idea of inexact solve, we show that only a small number of linear iterations are needed at each Picard step. The proposed algorithm is tested with a model of flow in a two-dimensional symmetric step domain with different settings to demonstrate the computational efficiency.


title: Recovering the parameters underlying the Lorenz-96 chaotic dynamics

abstract: Climate projections suffer from uncertain equilibrium climate sensitivity. The reason behind this uncertainty is the resolution of global climate models, which is too coarse to resolve key processes such as clouds and convection. These processes are approximated using heuristics in a process called parameterization. The selection of these parameters can be subjective, leading to significant uncertainties in the way clouds are represented in global climate models. Here, we explore three deep network algorithms to infer these parameters in an objective and data-driven way. We compare the performance of a fully-connected network, a one-dimensional and, a two-dimensional convolutional networks to recover the underlying parameters of the Lorenz-96 model, a non-linear dynamical system that has similar behavior to the climate system.


title: Magneto-hydrodynamical origin of eclipsing time variations in post-common-envelope binaries for solar mass secondaries

abstract: Eclipsing time variations have been observed for a wide range of binary systems, including post-common-envelope binaries. A frequently proposed explanation, apart from the possibility of having a third body, is the effect of magnetic activity, which may alter the internal structure of the secondary star, particularly its quadrupole moment, and thereby cause quasi-periodic oscillations. Here we present two compressible non-ideal magneto-hydrodynamical (MHD) simulations of the magnetic dynamo in a solar mass star, one of them with three times the solar rotation rate ("slow rotator"), the other one with twenty times the solar rotation rate ("rapid rotator"), to account for the high rotational velocities in close binary systems. For the slow rotator, we find that both the magnetic field and the stellar quadrupole moment change in a quasi-periodic manner, leading to O-C (observed - corrected times of the eclipse) variations of ~0.025 s. For the rapid rotator, the behavior of the magnetic field as well as the quadrupole moment changes become considerably more complex, due to the less coherent dynamo solution. The resulting O-C variations are of the order 0.13 s. The observed system V471~Tau shows two modes of eclipsing time variations, with amplitudes of 151 s and 20 s, respectively. However, the current simulations may not capture all relevant effects due to the neglect of the centrifugal force and self-gravity. Considering the model limitations and that the rotation of V471 Tau is still a factor of 2.5 faster than our rapid rotator, it may be conceivable to reach the observed magnitudes.


title: SEntNet: Source-aware Recurrent Entity Network for Dialogue Response Selection

abstract: Dialogue response selection is an important part of Task-oriented Dialogue Systems (TDSs); it aims to predict an appropriate response given a dialogue context. Obtaining key information from a complex, long dialogue context is challenging, especially when different sources of information are available, e.g., the user's utterances, the system's responses, and results retrieved from a knowledge base (KB). Previous work ignores the type of information source and merges sources for response selection. However, accounting for the source type may lead to remarkable differences in the quality of response selection. We propose the Source-aware Recurrent Entity Network (SEntNet), which is aware of different information sources for the response selection process. SEntNet achieves this by employing source-specific memories to exploit differences in the usage of words and syntactic structure from different information sources (user, system, and KB). Experimental results show that SEntNet obtains 91.0% accuracy on the Dialog bAbI dataset, outperforming prior work by 4.7%. On the DSTC2 dataset, SEntNet obtains an accuracy of 41.2%, beating source unaware recurrent entity networks by 2.4%.


title: Providentia -- A Large-Scale Sensor System for the Assistance of Autonomous Vehicles and Its Evaluation

abstract: The environmental perception of an autonomous vehicle is limited by its physical sensor ranges and algorithmic performance, as well as by occlusions that degrade its understanding of an ongoing traffic situation. This not only poses a significant threat to safety and limits driving speeds, but it can also lead to inconvenient maneuvers. Intelligent Infrastructure Systems can help to alleviate these problems. An Intelligent Infrastructure System can fill in the gaps in a vehicle's perception and extend its field of view by providing additional detailed information about its surroundings, in the form of a digital model of the current traffic situation, i.e. a digital twin. However, detailed descriptions of such systems and working prototypes demonstrating their feasibility are scarce. In this paper, we propose a hardware and software architecture that enables such a reliable Intelligent Infrastructure System to be built. We have implemented this system in the real world and demonstrate its ability to create an accurate digital twin of an extended highway stretch, thus enhancing an autonomous vehicle's perception beyond the limits of its on-board sensors. Furthermore, we evaluate the accuracy and reliability of the digital twin by using aerial images and earth observation methods for generating ground truth data.


title: Abundance of primordial black holes with local non-Gaussianity in peak theory

abstract: We discuss the effect of local type non-Gaussianity on the abundance of primordial black holes (PBH) based on the peak theory. We provide the PBH formation criterion based on the so-called compaction function and use the peak theory statistics associated with the curvature perturbation with the local type non-Gaussianity. Providing a method to estimate the PBH abundance, we demonstrate the effects of non-Gaussianity. It is explicitly shown that the value of non-linear parameter $|f_{\rm NL}| \sim 1$ induces a similar effect to a few factors of difference in the amplitude of the power spectrum.


title: Excitation of Kerr quasinormal modes in extreme--mass-ratio inspirals

abstract: If a small compact object orbits a black hole, it is known that it can excite the black hole's quasinormal modes (QNMs), leading to high-frequency oscillations ("wiggles") in the radiated field at $\mathcal{J}^+$, and in the radiation-reaction self-force acting on the object after its orbit passes through periapsis. Here we survey the phenomenology of these wiggles across a range of black hole spins and equatorial orbits. In both the scalar-field and gravitational cases we find that wiggles are a generic feature across a wide range of parameter space, and that they are observable in field perturbations at fixed spatial positions, in the self-force, and in radiated fields at $\mathcal{J}^+$. For a given charge or mass of the small body, the QNM excitations have the highest amplitudes for systems with a highly spinning central black hole, a prograde orbit with high eccentricity, and an orbital periapsis close to the light ring. The QNM amplitudes depend smoothly on the orbital parameters, with only very small amplitude changes when the orbit's (discrete) frequency spectrum is tuned to match QNM frequencies. The association of wiggles with QNM excitations suggest that they represent a situation where the \emph{nonlocal} nature of the self-force is particularly apparent, with the wiggles arising as result of QNM excitation by the compact object near periapsis, and then encountered later in the orbit. Astrophysically, the effects of wiggles at $\mathcal{J}^+$ might allow direct observation of Kerr QNMs in extreme-mass-ratio inspiral (EMRI) binary black hole systems, potentially enabling new tests of general relativity.


title: Floors are Flat: Leveraging Semantics for Real-Time Surface Normal Prediction

abstract: We propose 4 insights that help to significantly improve the performance of deep learning models that predict surface normals and semantic labels from a single RGB image. These insights are: (1) denoise the "ground truth" surface normals in the training set to ensure consistency with the semantic labels; (2) concurrently train on a mix of real and synthetic data, instead of pretraining on synthetic and finetuning on real; (3) jointly predict normals and semantics using a shared model, but only backpropagate errors on pixels that have valid training labels; (4) slim down the model and use grayscale instead of color inputs. Despite the simplicity of these steps, we demonstrate consistently improved results on several datasets, using a model that runs at 12 fps on a standard mobile phone.


title: Voltage matching, \'etendue and ratchet steps in advanced concept solar cells

abstract: Many advanced solar cell concepts propose surpassing the Shockley Queisser (SQ) limit by introducing multiple quasi-Fermi level separations that are arranged in series and/or in parallel. Exceeding the SQ limit with any parallel arrangement involves intermediate states that deliver additional charge carriers at, ideally, the same electro-chemical potential as the other elements in the parallel network. This can be thought of as voltage matching individual parallel components and in intermediate band materials is intricately linked to solar concentration and \'etendue mismatch between absorption and emission. Generally, to achieve voltage matching under sub-optimal conditions, an additional degree of freedom in the absorption thresholds of the material through a carrier relaxation or ratchet step is required. We explain why the ideal ratchet step decreases with solar concentration and how it depends on radiative efficiency and emission \'etendue of the individual transitions. For solar cell concepts that use Auger type carrier-carrier interactions or molecular triplet states for energetic up- or down-conversion, ideal bandgap combinations and achievable efficiencies also depend on interaction rates. We show that Auger assisted solar cells suffer more strongly from finite interaction rates than carrier multiplication devices.


title: Back-Projection based Fidelity Term for Ill-Posed Linear Inverse Problems

abstract: Ill-posed linear inverse problems appear in many image processing applications, such as deblurring, super-resolution and compressed sensing. Many restoration strategies involve minimizing a cost function, which is composed of fidelity and prior terms, balanced by a regularization parameter. While a vast amount of research has been focused on different prior models, the fidelity term is almost always chosen to be the least squares (LS) objective, that encourages fitting the linearly transformed optimization variable to the observations. In this paper, we examine a different fidelity term, which has been implicitly used by the recently proposed iterative denoising and backward projections (IDBP) framework. This term encourages agreement between the projection of the optimization variable onto the row space of the linear operator and the pseudo-inverse of the linear operator ("back-projection") applied on the observations. We analytically examine the difference between the two fidelity terms for Tikhonov regularization and identify cases (such as a badly conditioned linear operator) where the new term has an advantage over the standard LS one. Moreover, we demonstrate empirically that the behavior of the two induced cost functions for sophisticated convex and non-convex priors, such as total-variation, BM3D, and deep generative models, correlates with the obtained theoretical analysis.


title: Empirical Predictions for the Period Distribution of Planets to be Discovered by the Transiting Exoplanet Survey Satellite

abstract: Launched on April 2018, NASA's Transiting Exoplanet Survey Satellite (TESS) has been performing a wide-field survey for exoplanets orbiting stars with a goal of producing a rich database for follow-on studies. Here we present estimates of the detected exoplanet orbital periods in the 2-minute cadence mode during the TESS mission. For a two-transit detection criterion, the expected mean value of the most frequently detected orbital period is 5.01 days with the most frequently detected range of 2.12 to 11.82 days in the region with observation of 27 days. Near the poles where the observational duration is 351 days, the expected mean orbital period is 10.93 days with the most frequently detected range being from 3.35 to 35.65 days. For one-transit, the most frequently detected orbital period is 8.17 days in the region with observation of 27 days and 11.25 days near the poles. For the entire TESS mission containing several sectors, we estimate that the mean value of orbital period is 8.47 days for two-transit and 10.09 days for one-transit, respectively. If TESS yields a planet population substantially different from what's predicted here, the underlying planet occurrence rates are likely different between the stellar sample probed by TESS and that by Kepler.


title: ASAC: Active Sensing using Actor-Critic models

abstract: Deciding what and when to observe is critical when making observations is costly. In a medical setting where observations can be made sequentially, making these observations (or not) should be an active choice. We refer to this as the active sensing problem. In this paper, we propose a novel deep learning framework, which we call ASAC (Active Sensing using Actor-Critic models) to address this problem. ASAC consists of two networks: a selector network and a predictor network. The selector network uses previously selected observations to determine what should be observed in the future. The predictor network uses the observations selected by the selector network to predict a label, providing feedback to the selector network (well-selected variables should be predictive of the label). The goal of the selector network is then to select variables that balance the cost of observing the selected variables with their predictive power; we wish to preserve the conditional label distribution. During training, we use the actor-critic models to allow the loss of the selector to be "back-propagated" through the sampling process. The selector network "acts" by selecting future observations to make. The predictor network acts as a "critic" by feeding predictive errors for the selected variables back to the selector network. In our experiments, we show that ASAC significantly outperforms state-of-the-arts in two real-world medical datasets.


title: Impact of Stellar Superflares on Planetary Habitability

abstract: High-energy radiation caused by exoplanetary space weather events from planet-hosting stars can play a crucial role in conditions promoting or destroying habitability in addition to the conventional factors. In this paper, we present the first quantitative impact evaluation system of stellar flares on the habitability factors with an emphasis on the impact of Stellar Proton Events. We derive the maximum flare energy from stellar starspot sizes and examine the impacts of flare associated ionizing radiation on CO$_2$, H$_2$, N$_2$+O$_2$ --rich atmospheres of a number of well-characterized terrestrial type exoplanets. Our simulations based on the Particle and Heavy Ion Transport code System [PHITS] suggest that the estimated ground level dose for each planet in the case of terrestrial-level atmospheric pressure (1 bar) for each exoplanet does not exceed the critical dose for complex (multi-cellular) life to persist, even for the planetary surface of Proxima Centauri b, Ross-128 b and TRAPPIST-1 e. However, when we take into account the effects of the possible maximum flares from those host stars, the estimated dose reaches fatal levels at the terrestrial lowest atmospheric depth on TRAPPIST-1 e and Ross-128 b. Large fluxes of coronal XUV radiation from active stars induces high atmospheric escape rates from close-in exoplanets suggesting that the atmospheric depth can be substantially smaller than that on the Earth. In a scenario with the atmospheric thickness of 1/10 of Earth's, the radiation dose from close-in planets including Proxima Centauri b and TRAPPIST-1 e reach near fatal dose levels with annual frequency of flare occurrence from their hoststars.


title: Panoptic Image Annotation with a Collaborative Assistant

abstract: This paper aims to reduce the time to annotate images for panoptic segmentation, which requires annotating segmentation masks and class labels for all object instances and stuff regions. We formulate our approach as a collaborative process between an annotator and an automated assistant who take turns to jointly annotate an image using a predefined pool of segments. Actions performed by the annotator serve as a strong contextual signal. The assistant intelligently reacts to this signal by annotating other parts of the image on its own, which reduces the amount of work required by the annotator. We perform thorough experiments on the COCO panoptic dataset, both in simulation and with human annotators. These demonstrate that our approach is significantly faster than the recent machine-assisted interface of [4], and 2.4x to 5x faster than manual polygon drawing. Finally, we show on ADE20k that our method can be used to efficiently annotate new datasets, bootstrapping from a very small amount of annotated data.


title: Fully compensated Kondo effect for a two-channel spin S=1 impurity

abstract: We study the low-temperature properties of the generalized Anderson impurity model in which two localized configurations, one with two doublets and the other with a triplet, are mixed by two degenerate conduction channels. By using the numerical renormalization group and the non-crossing approximation, we analyze the impurity entropy, its spectral density, and the equilibrium conductance for several values of the model parameters. Marked differences with respect to the conventional one-channel spin $s=1/2$ Anderson model, that can be traced as hallmarks of an impurity spin $S=1$, are found in the Kondo temperature, the width and position of the charge transfer peak, as well as the temperature dependence of the equilibrium conductance. Furthermore, we analyze the rich effects of a single-ion magnetic anisotropy $D$ on the Kondo behavior. In particular, as shown before, for large enough positive $D$ the system behaves as a "non-Landau" Fermi liquid that cannot be adiabatically connected to a non-interacting system turning off the interactions. For negative $D$ the Kondo effect is strongly suppressed. The model studied is suitable for a comprehensive analysis for recent investigations of a single Ni impurity embedded into an Au chain.


title: On infinite iterations of the functor of idempotent probability measures

abstract: In this paper we establish that the functor of idempotent probability measures acting in the category of compacta and their continous mappings is perfect metrisable


title: Spatially Resolved Interstellar Medium and Highly-Excited Dense Molecular Gas in the Most Luminous Quasar at z=6.327

abstract: Among more than 200 quasars known at $z\gtrsim6$, only one object, J0100+2802 (z=6.327), was found hosting a $>10^{10}M_\odot$ super-massive black hole (SMBH). In order to investigate the host galaxy properties of J0100+2802, we performed multi-band ALMA observations, aiming at mapping the dust continuum, [CII] and CO(6-5) emission lines with sub-kiloparsec scale resolution, as well as detecting high-J CO lines in CO(11-10), CO(10-9), and CO(7-6). The galaxy size is measured to be $R_{\rm major}=3.6\pm0.2$ kpc from the high resolution continuum observations. No ordered motion on kilo-parsec scales was found in both [CII] and CO(6-5) emissions. The velocity dispersion is measured to be 161$\pm$7 km s$^{-1}$, $\sim$3 times smaller than that estimated from the local M-$\sigma$ relation. In addition, we found that the CO emission is more concentrate (a factor of 1.8$\pm$0.4) than the [CII] emission. Together with CO(2-1) detected by VLA, we measured the CO Spectral Line Energy Distribution (SLED), which is best fit by a two-components model, including a cool component at $\sim24$ K with a density of $n_{\rm (H_2)}$=10$^{4.5}$ cm$^{-3}$, and a warm component at $\sim224$ K with a density of $n_{\rm (H_2)}$=10$^{3.6}$ cm$^{-3}$. We also fitted the dust continuum with a graybody model, which indicates that it has either a high dust emissivity $\beta\gtrsim2$ or a hot dust temperature $T_{\rm dust}\gtrsim60$ K, or a combination of both factors. The highly excited CO emission and hot dust temperature suggest that the powerful AGN in J0100+2802 could contribute to the gas and dust heating although future observations are needed to confirm this.


title: Siegel disks of the tangent family

abstract: We study Siegel disks in the dynamics of functions from the tangent family. In particular, we prove that a forward invariant Siegel disk is unbounded if and only if it contains at least one asymptotic value on the boundary. Our argument is elementary and function-theoretic. Moreover, by using quasiconformal surgery we also construct functions in the above family with bounded Siegel disks.


title: Sticky Brownian Motion and its Numerical Solution

abstract: Sticky Brownian motion is the simplest example of a diffusion process that can spend finite time both in the interior of a domain and on its boundary. It arises in various applications such as in biology, materials science, and finance. This article spotlights the unusual behavior of sticky Brownian motions from the perspective of applied mathematics, and provides tools to efficiently simulate them. We show that a sticky Brownian motion arises naturally for a particle diffusing on $\mathbb{R}_+$ with a strong, short-ranged potential energy near the origin. This is a limit that accurately models mesoscale particles, those with diameters $\approx 100$nm-$10\mu$m, which form the building blocks for many common materials. We introduce a simple and intuitive sticky random walk to simulate sticky Brownian motion, that also gives insight into its unusual properties. In parameter regimes of practical interest, we show this sticky random walk is two to five orders of magnitude faster than alternative methods to simulate a sticky Brownian motion. We outline possible steps to extend this method towards simulating multi-dimensional sticky diffusions.

categories: math.NA cond-mat.stat-mech cs.NA math.PR


title: Three-Dimensional Fourier Scattering Transform and Classification of Hyperspectral Images

abstract: Recent developments in machine learning and signal processing have resulted in many new techniques that are able to effectively capture the intrinsic yet complex properties of hyperspectral imagery. Tasks ranging from anomaly detection to classification can now be solved by taking advantage of very efficient algorithms which have their roots in representation theory and in computational approximation. Time-frequency methods are one example of such techniques. They provide means to analyze and extract the spectral content from data. On the other hand, hierarchical methods such as neural networks incorporate spatial information across scales and model multiple levels of dependencies between spectral features. Both of these approaches have recently been proven to provide significant advances in the spectral-spatial classification of hyperspectral imagery. The 3D Fourier scattering transform, which is introduced in this paper, is an amalgamation of time-frequency representations with neural network architectures. It leverages the benefits provided by the Short-Time Fourier Transform with the numerical efficiency of deep learning network structures. We test the proposed method on several standard hyperspectral datasets, and we present results that indicate that the 3D Fourier scattering transform is highly effective at representing spectral content when compared with other state-of-the-art spectral-spatial classification methods.


title: Neural Theorem Provers Do Not Learn Rules Without Exploration

abstract: Neural symbolic processing aims to combine the generalization of logical learning approaches and the performance of neural networks. The Neural Theorem Proving (NTP) model by Rocktaschel et al (2017) learns embeddings for concepts and performs logical unification. While NTP is promising and effective in predicting facts accurately, we have little knowledge how well it can extract true relationship among data. To this end, we create synthetic logical datasets with injected relationships, which can be generated on-the-fly, to test neural-based relation learning algorithms including NTP. We show that it has difficulty recovering relationships in all but the simplest settings. Critical analysis and diagnostic experiments suggest that the optimization algorithm suffers from poor local minima due to its greedy winner-takes-all strategy in identifying the most informative structure (proof path) to pursue. We alter the NTP algorithm to increase exploration, which sharply improves performance. We argue and demonstate that it is insightful to benchmark with synthetic data with ground-truth relationships, for both evaluating models and revealing algorithmic issues.


title: SU(2) and SU(1,1) Y-Maps in Loop Quantum Gravity

abstract: In this paper we first provide the proof of $SU(2)$ Y-Map convergence. Then, by using $SU(1,1)$ LQG simplicity constraints we define $SU(1,1)$ Y-Map from infinitely differentiable with a compact support functions on $SU(1,1)$ to the functions (not necessarily square integrable) on $SL(2,C)$, and prove its convergence as well.


title: Curl-free magnetic fields for stellarator optimization

abstract: This paper describes a new and efficient method of defining an annular region of a curl-free magnetic field with specific physics and coil properties that can be used in stellarator design. Three statements define the importance: (1) Codes can follow an optimized curl-free initial state to a final full-pressure equilibrium. The large size of the optimization space of stellarators, approximately fifty externally-produced distributions of magnetic field, makes success in finding a global optimum largely determined by the starting point. (2) The design of a stellarator is actually improved when the central region of the plasma has rapid transport with the confinement provided by a surrounding annulus of magnetic surfaces with low transport. (3) The stellarator is unique among all fusion concepts, inertial as well as magnetic, in not using the plasma itself to provide an essential part of its confinement concept. This permits reliable computational design, which opens a path to faster, cheaper, and more certain achievement of fusion energy.


title: Evaluation of frictional suspension system performance in three-piece bogie by the multi-body-system approach

abstract: The nonlinear damping characteristics of friction wedges in the secondary suspension of a freight bogie are investigated considering unidirectional contact and non-smooth frictional forces. In the proposed MultiBody System model, each wedge has six Degrees of Freedom with corresponding inertial properties. The geometry of wedge, as well as, wedge angle, toe-in condition, and clearances between wedge/bolster and wedge/side frame are also considered in the modeling. The methodology for the identification of contact parameters is presented to achieve a smooth response and efficient numerical solutions. The response of the system with the proposed method is compared with the quasi-static methodology of wedge simulation and limitations of alternative methods are highlighted. The model is then applied to evaluate the effects of different design parameters of the wedge system on vertical and lateral response of the system.


title: Neutron diffraction strain tomography: 2D axisymmetric sample geometry

abstract: We consider a method for neutron residual strain/stress analysis combining features of transmission strain tomography (e.g. through Braggs edge imaging) and neutron diffraction techniques. Analogous to transmission strain tomography, our approach observes integrals of strain along the path of a beam, however, rather than transmitted beams, we observe a diffracted beam. The practical feasibility of neutron diffraction tomography is demonstrated in case of a 2D axisymmetric sample where the reconstruction of the stress field is determined from the measurement of a strain profile. The reconstruction is compared to a set of measurements made using a traditional neutron stress scanning technique in order to validate the new method. Similar strain accuracy was achieved by the new technique within a reasonable measurement time. The practical and theoretical challenges regarding the technique are also discussed.


title: General equilibrium in a heterogeneous-agent incomplete-market economy with many consumption goods and a risk-free bond

abstract: We study a pure-exchange incomplete-market economy with heterogeneous agents. In each period, the agents choose how much to save (i.e., invest in a risk-free bond), how much to consume, and which bundle of goods to consume while their endowments are fluctuating. We focus on a competitive stationary equilibrium (CSE) in which the wealth distribution is invariant, the agents maximize their expected discounted utility, and both the prices of consumption goods and the interest rate are market-clearing. Our main contribution is to extend some general equilibrium results to an incomplete-market Bewley-type economy with many consumption goods. Under mild conditions on the agents' preferences, we show that the aggregate demand for goods depends only on their relative prices and that the aggregate demand for savings is homogeneous of degree in prices, and we prove the existence of a CSE. When the agents' preferences can be represented by a CES (constant elasticity of substitution) utility function with an elasticity of substitution that is higher than or equal to one, we prove that the CSE is unique. Under the same preferences, we show that a higher inequality of endowments does not change the equilibrium prices of goods, and decreases the equilibrium interest rate. Our results shed light on the impact of market incompleteness on the properties of general equilibrium models.


title: Using historical utility outage data to compute overall transmission grid resilience

abstract: Given increasing risk from climate-induced natural hazards, there is growing interest in the development of methods that can quantitatively measure resilience in power systems. This work quantifies resilience in electric power transmission networks in a new and comprehensive way that can represent the multiple processes of resilience. A novel aspect of this approach is the use of empirical data to develop the probability distributions that drive the model. This paper demonstrates the approach by measuring the impact of one potential improvement to a power system. Specifically, we measure the impact of additional distributed generation on power system resilience.


title: A gray-box approach for curriculum learning

abstract: Curriculum learning is often employed in deep reinforcement learning to let the agent progress more quickly towards better behaviors. Numerical methods for curriculum learning in the literature provides only initial heuristic solutions, with little to no guarantee on their quality. We define a new gray-box function that, including a suitable scheduling problem, can be effectively used to reformulate the curriculum learning problem. We propose different efficient numerical methods to address this gray-box reformulation. Preliminary numerical results on a benchmark task in the curriculum learning literature show the viability of the proposed approach.


title: A Temporal Sequence Learning for Action Recognition and Prediction

abstract: In this work\footnote {This work was supported in part by the National Science Foundation under grant IIS-1212948.}, we present a method to represent a video with a sequence of words, and learn the temporal sequencing of such words as the key information for predicting and recognizing human actions. We leverage core concepts from the Natural Language Processing (NLP) literature used in sentence classification to solve the problems of action prediction and action recognition. Each frame is converted into a word that is represented as a vector using the Bag of Visual Words (BoW) encoding method. The words are then combined into a sentence to represent the video, as a sentence. The sequence of words in different actions are learned with a simple but effective Temporal Convolutional Neural Network (T-CNN) that captures the temporal sequencing of information in a video sentence. We demonstrate that a key characteristic of the proposed method is its low-latency, i.e. its ability to predict an action accurately with a partial sequence (sentence). Experiments on two datasets, \textit{UCF101} and \textit{HMDB51} show that the method on average reaches 95\% of its accuracy within half the video frames. Results, also demonstrate that our method achieves compatible state-of-the-art performance in action recognition (i.e. at the completion of the sentence) in addition to action prediction.


title: Wiener index and Harary index on pancyclic graphs

abstract: Wiener index and Harary index are two classic and well-known topological indices for the characterization of molecular graphs. Recently, Yu et al. \cite{YYSX} established some sufficient conditions for a graph to be pancyclic in terms of the edge number, the spectral radius and the signless Laplacian spectral radius of the graph. In this paper, we give some sufficient conditions for a graph being pancyclic in terms of the Wiener index, the Harary index, the distance spectral radius and the Harary spectral radius of a graph.


title: Dual-band in situ molecular spectroscopy using single-sized Al-disk perfect absorbers

abstract: In this manuscript, we have demonstrated that the antenna-enhanced infrared vibrational spectroscopy can be a powerful platform for high-sensitivity detection of ultrathin layer of polymer molecules and in situ monitoring of their crosslinking reaction as well as their kinetics of network formation. Our precisely designed dual-band perfect absorber (DPA) made of singles-sized resonator array shows strongly-enhanced nearfield intensity and polarization independence, providing a powerful antenna platform for the bond-selective vibrational sensing.


title: A new approach to forecast service parts demand by integrating user preferences into multi-objective optimization

abstract: Service supply chain management is to prepare spare parts for failed products under warranty. Their goal is to reach agreed service level at the minimum cost. We convert this business problem into a preference based multi-objective optimization problem, where two quality criteria must be simultaneously optimized. One criterion is accuracy of demand forecast and the other is service level. Here we propose a general framework supporting solving preference-based multi-objective optimization problems (MOPs) by multi-gradient descent algorithm (MGDA), which is well suited for training deep neural network. The proposed framework treats agreed service level as a constrained criterion that must be met and generate a Pareto-optimal solution with highest forecasting accuracy. The neural networks used here are two Encoder-Decoder LSTM modes: one is used for pre-training phase to learn distributed representation of former generations' service parts consumption data, and the other is used for supervised learning phase to generate forecast quantities of current generations' service parts. Evaluated under the service parts consumption data in Lenovo Group Ltd, the proposed method clearly outperform baseline methods.


title: Gr\"obner bases and the Cohen-Macaulay property of Li's double determinantal varieties

abstract: We consider double determinantal varieties, a special case of Nakajima quiver varieties. Li conjectured that double determinantal varieties are normal, irreducible, Cohen-Macaulay varieties whose defining ideals have a Gr\"obner basis given by their natural generators. We use liaison theory to prove this conjecture in a manner that generalizes results for mixed ladder determinantal varieties. We also give a formula for the dimension of a double determinantal variety.


title: Stacked Capsule Autoencoders

abstract: Objects are composed of a set of geometrically organized parts. We introduce an unsupervised capsule autoencoder (SCAE), which explicitly uses geometric relationships between parts to reason about objects. Since these relationships do not depend on the viewpoint, our model is robust to viewpoint changes. SCAE consists of two stages. In the first stage, the model predicts presences and poses of part templates directly from the image and tries to reconstruct the image by appropriately arranging the templates. In the second stage, SCAE predicts parameters of a few object capsules, which are then used to reconstruct part poses. Inference in this model is amortized and performed by off-the-shelf neural encoders, unlike in previous capsule networks. We find that object capsule presences are highly informative of the object class, which leads to state-of-the-art results for unsupervised classification on SVHN (55%) and MNIST (98.7%). The code is available at https://github.com/google-research/google-research/tree/master/stacked_capsule_autoencoders


title: A Fusion Adversarial Underwater Image Enhancement Network with a Public Test Dataset

abstract: Underwater image enhancement algorithms have attracted much attention in underwater vision task. However, these algorithms are mainly evaluated on different data sets and different metrics. In this paper, we set up an effective and pubic underwater test dataset named U45 including the color casts, low contrast and haze-like effects of underwater degradation and propose a fusion adversarial network for enhancing underwater images. Meanwhile, the well-designed the adversarial loss including Lgt loss and Lfe loss is presented to focus on image features of ground truth, and image features of the image enhanced by fusion enhance method, respectively. The proposed network corrects color casts effectively and owns faster testing time with fewer parameters. Experiment results on U45 dataset demonstrate that the proposed method achieves better or comparable performance than the other state-of-the-art methods in terms of qualitative and quantitative evaluations. Moreover, an ablation study demonstrates the contributions of each component, and the application test further shows the effectiveness of the enhanced images.


title: Robust Weyl points in a 1D superlattice with transverse spin-orbit coupling

abstract: Weyl points, synthetic magnetic monopoles in the 3D momentum space, are the key features of topological Weyl semimetals. The observation of Weyl points in ultracold atomic gases usually relies on the realization of high-dimensional spin-orbit coupling (SOC) for two pseudospin states (% \textit{i.e.,} spin-1/2), which requires complex laser configurations and precise control of laser parameters, thus has not been realized in experiment. Here we propose that robust Wely points can be realized using 1D triple-well superlattices (spin-1/three-band systems) with 2D transverse SOC achieved by Raman-assisted tunnelings. The presence of the third band is responsible to the robustness of the Weyl points against system parameters (e.g., Raman laser polarization, phase, incident angle, etc.). Different from a spin-1/2 system, the non-trivial topology of Weyl points in such spin-1 system is characterized by both the spin vector and tensor textures, which can be probed using momentum-resolved Rabi spectroscopy. Our proposal provides a simple yet powerful platform for exploring Weyl physics and related high-dimensional topological phenomena using high pseudospin ultracold atoms.


title: A Survey of Optimization Methods from a Machine Learning Perspective

abstract: Machine learning develops rapidly, which has made many theoretical breakthroughs and is widely applied in various fields. Optimization, as an important part of machine learning, has attracted much attention of researchers. With the exponential growth of data amount and the increase of model complexity, optimization methods in machine learning face more and more challenges. A lot of work on solving optimization problems or improving optimization methods in machine learning has been proposed successively. The systematic retrospect and summary of the optimization methods from the perspective of machine learning are of great significance, which can offer guidance for both developments of optimization and machine learning research. In this paper, we first describe the optimization problems in machine learning. Then, we introduce the principles and progresses of commonly used optimization methods. Next, we summarize the applications and developments of optimization methods in some popular machine learning fields. Finally, we explore and give some challenges and open problems for the optimization in machine learning.


title: Spatio-Temporal Fusion Networks for Action Recognition

abstract: The video based CNN works have focused on effective ways to fuse appearance and motion networks, but they typically lack utilizing temporal information over video frames. In this work, we present a novel spatio-temporal fusion network (STFN) that integrates temporal dynamics of appearance and motion information from entire videos. The captured temporal dynamic information is then aggregated for a better video level representation and learned via end-to-end training. The spatio-temporal fusion network consists of two set of Residual Inception blocks that extract temporal dynamics and a fusion connection for appearance and motion features. The benefits of STFN are: (a) it captures local and global temporal dynamics of complementary data to learn video-wide information; and (b) it is applicable to any network for video classification to boost performance. We explore a variety of design choices for STFN and verify how the network performance is varied with the ablation studies. We perform experiments on two challenging human activity datasets, UCF101 and HMDB51, and achieve the state-of-the-art results with the best network.


title: Topological Phase Transition in an Archetypal f-electron Correlated System: Ce

abstract: A typical f-electron Kondo lattice system Ce exhibits the well-known isostructural transition, the so-called gamma-alpha transition, accompanied by an enormous volume collapse. Most interestingly, we have discovered that a topological-phase transition also takes place in elemental Ce, concurrently with the gamma-alpha transition. Based on the dynamical mean-field theory approach combined with density functional theory, we have unravelled that the non-trivial topology in alpha-Ce is driven by the f-d band inversion, which arises from the formation of coherent 4f band around the Fermi level. We captured the formation of the 4f quasi-particle band that is responsible for the Lifshitz transition and the non-trivial Z2 topology establishment across the phase boundary. This discovery provides a concept of 'topology switch' for topological Kondo systems. The 'on' and 'off' switching knob in Ce is versatile in a sense that it is controlled by available pressure (around 1 GPa) at room temperature.


title: Maximal Cohen-Macaulay modules over a noncommutative 2-dimensional singularity

abstract: We study properties of graded maximal Cohen-Macaulay modules over an N-graded locally finite, Auslander Gorenstein, and Cohen-Macaulay algebra of dimension two. As a consequence, we extend a part of McKay correspondence in dimension two to a more general setting.


title: Sudakov suppression of the Balitsky-Kovchegov kernel

abstract: To sum high energy leading logarithms in a consistent way, one has to impose the strong ordering in both projectile rapidity and dense target rapidity simultaneously, which results in a kinematically improved Balitsky-Kovchegov(BK) equation. We find that beyond this strong ordering region, the important sub-leading double logarithms arise at high order due to the incomplete cancellation between real corrections and virtual corrections in a t-channel calculation. Based on this observation, we further argue that these double logarithms are the Sudakov type ones, and thus can be resummed into an exponential leading to a Sudakov suppressed BK equation.


title: Homogeneous Network Embedding for Massive Graphs via Reweighted Personalized PageRank

abstract: Given an input graph G and a node v in G, homogeneous network embedding (HNE) maps the graph structure in the vicinity of v to a compact, fixed-dimensional feature vector. This paper focuses on HNE for massive graphs, e.g., with billions of edges. On this scale, most existing approaches fail, as they incur either prohibitively high costs, or severely compromised result utility. Our proposed solution, called Node-Reweighted PageRank (NRP), is based on a classic idea of deriving embedding vectors from pairwise personalized PageRank (PPR) values. Our contributions are twofold: first, we design a simple and efficient baseline HNE method based on PPR that is capable of handling billion-edge graphs on commodity hardware; second and more importantly, we identify an inherent drawback of vanilla PPR, and address it in our main proposal NRP. Specifically, PPR was designed for a very different purpose, i.e., ranking nodes in G based on their relative importance from a source node's perspective. In contrast, HNE aims to build node embeddings considering the whole graph. Consequently, node embeddings derived directly from PPR are of suboptimal utility. The proposed NRP approach overcomes the above deficiency through an effective and efficient node reweighting algorithm, which augments PPR values with node degree information, and iteratively adjusts embedding vectors accordingly. Overall, NRP takes O(mlogn) time and O(m) space to compute all node embeddings for a graph with m edges and n nodes. Our extensive experiments that compare NRP against 18 existing solutions over 7 real graphs demonstrate that NRP achieves higher result utility than all the solutions for link prediction, graph reconstruction and node classification, while being up to orders of magnitude faster. In particular, on a billion-edge Twitter graph, NRP terminates within 4 hours, using a single CPU core.


title: Atiyah-Segal Derived Completions for Equivariant Algebraic G-Theory and K-Theory

abstract: In the mid 1980s, while working on establishing completion theorems for equivariant Algebraic K- Theory similar to the well-known Atiyah-Segal completion theorem for equivariant topological K-theory, the late Robert Thomason found the strong finiteness conditions that are required in such theorems to be too restrictive. Then he made a conjecture on the existence of a completion theorem in the sense of Atiyah and Segal for equivariant Algebraic G-theory, for actions of linear algebraic groups on schemes that holds without any of the strong finiteness conditions that are required in such theorems proven by him, and also appearing in the original Atiyah-Segal theorem. The main goal of the present paper is to provide a proof of this conjecture in as broad a context as possible, making use of the technique of derived completion, and to consider several of the applications. Our solution is broad enough to allow actions by all linear algebraic groups, irrespective of whether they are connected or not, and acting on any quasi-projective scheme of finite type over a field, irrespective of whether they are regular or projective. This allows us therefore to consider the Equivariant Algebraic G-Theory of large classes of varieties like all Toric varieties (for the action of a torus) and all Spherical varieties (for the action of a reductive group). Restricting to actions by split tori, we are also able to consider actions on Algebraic Spaces. These enable us to obtain a wide range of applications, some of which are briefly sketched and which we plan to explore in detail in the future. A comparison of our results with previously known results, none of which made use of derived completions, shows that without the use of derived completions, one can only obtain results which are indeed very restrictive.


title: An Optimal Test for the Additive Model with Discrete or Categorical Predictors

abstract: In multivariate nonparametric regression the additive models are very useful when a suitable parametric model is difficult to find. The backfitting algorithm is a powerful tool to estimate the additive components. However, due to complexity of the estimators, the asymptotic $p$-value of the associated test is difficult to calculate without a Monte Carlo simulation. Moreover, the conventional tests assume that the predictor variables are strictly continuous. In this paper, a new test is introduced for the additive components with discrete or categorical predictors, where the model may contain continuous covariates. This method is also applied to the semiparametric regression to test the goodness-of-fit of the model. These tests are asymptotically optimal in terms of the rate of convergence, as they can detect a specific class of contiguous alternatives at a rate of $n^{-1/2}$. An extensive simulation study is presented to support the theoretical results derived in this paper. Finally, the method is applied to a real data to model the diamond price based on its quality attributes and physical measurements.


title: A multigrid-reduction-in-time solver with a new two-level convergence for unsteady fractional Laplacian problems

abstract: The multigrid-reduction-in-time (MGRIT) technique has proven to be successful in achieving higher run-time speedup by exploiting parallelism in time. The goal of this article is to develop and analyze a MGRIT algorithm, using FCF-relaxation with time-dependent time-grid propagators, to seek the finite element approximations of unsteady fractional Laplacian problems. The multigrid with line smoother proposed in [L. Chen, R. H. Nochetto, E. Ot{\'a}rola, A. J. Salgado, Math. Comp. 85 (2016) 2583--2607] is chosen to be the spatial solver. Motivated by [B. S. Southworth, SIAM J. Matrix Anal. Appl. 40 (2019) 564--608], we provide a new temporal eigenvalue approximation property and then deduce a generalized two-level convergence theory which removes the previous unitary diagonalization assumption on the fine and coarse time-grid propagators required in [X. Q. Yue, S. Shu, X. W. Xu, W. P. Bu, K. J. Pan, Comput. Math. Appl. 78 (2019) 3471--3484]. Numerical computations are included to confirm the theoretical predictions and demonstrate the sharpness of the derived convergence upper bound.


title: Understanding Natural Language Instructions for Fetching Daily Objects Using GAN-Based Multimodal Target-Source Classification

abstract: In this paper, we address multimodal language understanding for unconstrained fetching instruction in domestic service robots context. A typical fetching instruction such as "Bring me the yellow toy from the white shelf" requires to infer the user intention, that is what object (target) to fetch and from where (source). To solve the task, we propose a Multimodal Target-source Classifier Model (MTCM), which predicts the region-wise likelihood of target and source candidates in the scene. Unlike other methods, MTCM can handle regionwise classification based on linguistic and visual features. We evaluated our approach that outperformed the state-of-the-art method on a standard data set. In addition, we extended MTCM with Generative Adversarial Nets (MTCM-GAN), and enabled simultaneous data augmentation and classification.


title: Baryon Physics and Tight Coupling Approximation in Boltzmann Codes

abstract: We provide two derivations of the baryonic equations that can be straightforwardly implemented in existing Einstein--Boltzmann solvers. One of the derivations begins with an action principle, while the other exploits the conservation of the stress-energy tensor. While our result is manifestly covariant and satisfies the Bianchi identities, we point out that this is not the case for the implementation of the seminal work by Ma and Bertschinger and in the existing Boltzmann codes. We also study the tight coupling approximation up to the second order without choosing any gauge using the covariant full baryon equations. We implement the improved baryon equations in a Boltzmann code and investigate the change in the estimate of cosmological parameters by performing an MCMC analysis. With the covariantly correct baryon equations of motion, we find 1% deviation for the best fit values of the cosmological parameters that should be taken into account. While in this paper, we study the Lambda-CDM model only, our baryon equations can be easily implemented in other models and various modified gravity theories.


title: Sample-Efficient Neural Architecture Search by Learning Action Space

abstract: Neural Architecture Search (NAS) has emerged as a promising technique for automatic neural network design. However, existing MCTS based NAS approaches often utilize manually designed action space, which is not directly related to the performance metric to be optimized (e.g., accuracy), leading to sample-inefficient explorations of architectures. To improve the sample efficiency, this paper proposes Latent Action Neural Architecture Search (LaNAS), which learns actions to recursively partition the search space into good or bad regions that contain networks with similar performance metrics. During the search phase, as different action sequences lead to regions with different performance, the search efficiency can be significantly improved by biasing towards the good regions. On three NAS tasks, empirical results demonstrate that LaNAS is at least an order more sample efficient than baseline methods including evolutionary algorithms, Bayesian optimizations, and random search. When applied in practice, both one-shot and regular LaNAS consistently outperform existing results. Particularly, LaNAS achieves 99.0% accuracy on CIFAR-10 and 80.8% top1 accuracy at 600 MFLOPS on ImageNet in only 800 samples, significantly outperforming AmoebaNet with 33x fewer samples. Our code is publicly available at https://github.com/facebookresearch/LaMCTS.


title: Low-temperature thermal conductivity and magnetic transitions of the kagom\'{e}-staircase compound Ni$_3$V$_2$O$_8$

abstract: Kagom\'{e}-staircase compound Ni$_3$V$_2$O$_8$ is an attractive multiferroic material exhibiting rich phase diagrams. However, the magnetic properties and magnetic transitions have been studied only above 1.3 K. In this work, we study the thermal conductivity ($\kappa$) of Ni$_3$V$_2$O$_8$ single crystals at low temperatures down to 0.3 K and in magnetic fields up to 14 T. In zero field, the magnetic transitions from the low-temperature incommensurate (LTI) phase to the commensurate phase (C) and then to a second commensurate phase (C') yield anomalies in $\kappa(T)$ curves at $T\rm_{LC}$ = 3.7 K and $T\rm_{CC'}$ = 2.0 K, respectively, which indicates a significant phonon scattering by the critical spin fluctuations. When the field is applied along the $a$ axis, the field dependence of $\kappa$ displays four anomalies associated with different magnetic transitions and reveals an undetected magnetic state at subKelvin temperatures. In addition, the $\kappa(B)$ curves are found to depend not only on the history but also on the magnitude of applying field. When the field is applied along the $b$ axis, a high-field phase locating above the LTI and high-temperature incommensurate (HTI) phases is revealed.


title: NLH: A Blind Pixel-level Non-local Method for Real-world Image Denoising

abstract: Non-local self similarity (NSS) is a powerful prior of natural images for image denoising. Most of existing denoising methods employ similar patches, which is a patch-level NSS prior. In this paper, we take one step forward by introducing a pixel-level NSS prior, i.e., searching similar pixels across a non-local region. This is motivated by the fact that finding closely similar pixels is more feasible than similar patches in natural images, which can be used to enhance image denoising performance. With the introduced pixel-level NSS prior, we propose an accurate noise level estimation method, and then develop a blind image denoising method based on the lifting Haar transform and Wiener filtering techniques. Experiments on benchmark datasets demonstrate that, the proposed method achieves much better performance than previous non-deep methods, and is still competitive with existing state-of-the-art deep learning based methods on real-world image denoising. The code is publicly available at https://github.com/njusthyk1972/NLH.


title: Minimax Density Estimation on Sobolev Spaces With Dominating Mixed Smoothness

abstract: We study minimax density estimation on the product space $\mathbb{R}^{d_1}\times\mathbb{R}^{d_2}$. We consider $L^p$-risk for probability density functions defined over regularity spaces that allow for different level of smoothness in each of the variables. Precisely, we study probabilities on Sobolev spaces with dominating mixed-smoothness. We provide the rate of convergence that is optimal even for the classical Sobolev spaces.


title: Multi-type Resource Allocation with Partial Preferences

abstract: We propose multi-type probabilistic serial (MPS) and multi-type random priority (MRP) as extensions of the well known PS and RP mechanisms to the multi-type resource allocation problem (MTRA) with partial preferences. In our setting, there are multiple types of divisible items, and a group of agents who have partial order preferences over bundles consisting of one item of each type. We show that for the unrestricted domain of partial order preferences, no mechanism satisfies both sd-efficiency and sd-envy-freeness. Notwithstanding this impossibility result, our main message is positive: When agents' preferences are represented by acyclic CP-nets, MPS satisfies sd-efficiency, sd-envy-freeness, ordinal fairness, and upper invariance, while MRP satisfies ex-post-efficiency, sd-strategy-proofness, and upper invariance, recovering the properties of PS and RP.


title: Why pushing a bell does not produce a sound

abstract: The difference between "beating" and "pushing" results in the perception that a push just makes the object move as a whole, while a beat produces also a sound. Through a detailed analysis of the physics underlying such everyday experiences, we identify the \emph{strength}, the \emph{duration} and the \emph{softness} of the applied contact force, as the main (measurable) characteristics that mark such difference. The strength determines the final velocity $\Delta v$ achieved by the body. The duration $2\tau$ compares to the time $\tau_\ell$ the sound takes to cross the body. The softness $\gamma$ (a positive exponent) results from the shape in time of the contact force. Those three elements enter the formula for the intensity of the sound produced. The relevant role of the softness is stressed and specific values are calculated for a thin metallic bar, chosen as the simplest possible model system.


title: Accelerating field enhancement due to ion motion in plasma wakefield accelerators

abstract: Ion motion in plasma wakefield accelerators can cause temporal increase of the longitudinal electric field shortly before the wave breaks. The increase is caused by re-distribution of the wave energy in transverse direction and may be important for correct interpretation of experimental results and acceleration of high-quality beams.


title: Presupernova neutrino signals as potential probes of neutrino mass hierarchy

abstract: We assess the potential of using presupernova neutrino signals at the Jiangmen Underground Neutrino Observatory (JUNO) to probe the yet-unknown neutrino mass hierarchy. Using models for stars of 12, 15, 20, and 25 solar masses, we find that if the electron antineutrino signals from such a star can be predicted precisely and the star is within ~440-880 pc, the number of events of electron antineutrino captures on protons detected within one day of its explosion allows to determine the hierarchy at the > ~95% confidence level. For determination at this level using such signals from Betelgeuse, which is at a distance of ~222 pc, the uncertainty in the predicted number of signals needs to be < ~14-30%. In view of more realistic uncertainties, we discuss and advocate a model-independent determination using both electron neutrino and antineutrino signals from Betelgeuse. This method is feasible if the cosmogenic background for neutrino-electron scattering events can be reduced by a factor of ~2.5-10 from the current estimate. Such reduction might be achieved by using coincidence of the background events, the exploration of which for JUNO is highly desirable.


title: Algebraization of Mochizuki's anabelian variation of ring structures, perfectoid geometry and formal groups

abstract: Let $M$ be a multiplicative monoid with identity. Then I show that there is a universal one dimensional formal group law equipped with an action of $M$. If $M$ is $p$-perfect (i.e. $m\mapsto m^p$ is an isomorphism for some prime number $p$) then the universal $M$-formal group law comes equipped with a natural Frobenius endomorphism. There are a number of concrete applications of this result. If $K$ is a $p$-adic field and $\mathcal{O}=\mathcal{O}_K$ is the multiplicative monoid of the ring of integers of $K$, then there is a universal formal group (over a suitable (non-zero) ring) which is equipped with an action of the multiplicative monoid $\mathcal{O}$. Lubin-Tate formal groups arise from this universal monoid formal group law. This has applications to Mochizuki's anabelian ideas: if two p-adic fields have isomorphic absolute Galois groups then they have isomorphic multiplicative monoids $\mathcal{O}$ (but possibly non-isomorphic ring structures). The existence of the universal monoid formal group law for the monoid $\mathcal{O}$ implies that the additive structures of a ring can be interpolated into a universal algebraic family (while keeping the multiplicative structure of the ring fixed). Here is another important example covered by my result: let $R$ be a perfectoid ring and let $R^\flat$ be its tilt and the multiplicative monoid $R^\flat$ of $R^\flat$. Then there exists a universal monoid formal group law for this monoid which interpolates the additive structures of untilts with tilt $R^\flat$. Thus in some sense one has a unified approach to various phenomenon which are well-known in anabelian geometry and in perfectoid geometry.


title: LPaintB: Learning to Paint from Self-Supervision

abstract: We present a novel reinforcement learning-based natural media painting algorithm. Our goal is to reproduce a reference image using brush strokes and we encode the objective through observations. Our formulation takes into account that the distribution of the reward in the action space is sparse and training a reinforcement learning algorithm from scratch can be difficult. We present an approach that combines self-supervised learning and reinforcement learning to effectively transfer negative samples into positive ones and change the reward distribution. We demonstrate the benefits of our painting agent to reproduce reference images with brush strokes. The training phase takes about one hour and the runtime algorithm takes about 30 seconds on a GTX1080 GPU reproducing a 1000x800 image with 20,000 strokes.


title: Spin-phonon interaction increased by compressive strain in antiferromagnetic MnO thin films

abstract: MnO thin films with various thicknesses and strains were grown on MgO substrates by pulsed laser deposition, then characterized using x-ray diffraction and infrared reflectance spectroscopy. Films grown on (001)-oriented MgO substrates exhibit homogenous biaxial compressive strain which increases as the film thickness is reduced. For that reason, the frequency of doubly-degenerate phonon increases with the strain, and splits below N\'eel temperature TN due to the magnetic-exchange interaction. Films grown on (110)-oriented MgO substrates exhibit a huge phonon splitting already at room temperature due to the anisotropic in-plane compressive strain. Below TN, additional phonon is activated in the IR spectra; this trend is evidence for a spin-order-induced structural phase transition from tetragonal to monoclinic phase. Total phonon splitting is 55 cm-1 in (110)-oriented MnO film, which is more than twice the value in bulk MnO. This result is evidence that the nearest neighbor exchange interaction, which is responsible for the magnetically driven phonon splitting, is greatly increased in compressively strained films.


title: Predicting Research Trends with Semantic and Neural Networks with an application in Quantum Physics

abstract: The vast and growing number of publications in all disciplines of science cannot be comprehended by a single human researcher. As a consequence, researchers have to specialize in narrow sub-disciplines, which makes it challenging to uncover scientific connections beyond the own field of research. Thus access to structured knowledge from a large corpus of publications could help pushing the frontiers of science. Here we demonstrate a method to build a semantic network from published scientific literature, which we call SemNet. We use SemNet to predict future trends in research and to inspire new, personalized and surprising seeds of ideas in science. We apply it in the discipline of quantum physics, which has seen an unprecedented growth of activity in recent years. In SemNet, scientific knowledge is represented as an evolving network using the content of 750,000 scientific papers published since 1919. The nodes of the network correspond to physical concepts, and links between two nodes are drawn when two physical concepts are concurrently studied in research articles. We identify influential and prize-winning research topics from the past inside SemNet thus confirm that it stores useful semantic knowledge. We train a deep neural network using states of SemNet of the past, to predict future developments in quantum physics research, and confirm high quality predictions using historic data. With the neural network and theoretical network tools we are able to suggest new, personalized, out-of-the-box ideas, by identifying pairs of concepts which have unique and extremal semantic network properties. Finally, we consider possible future developments and implications of our findings.

categories: cs.DL cs.LG physics.hist-ph physics.soc-ph quant-ph


title: Imaging Antiferromagnetic Domains in Nickel-oxide Thin Films by Magneto-optical Voigt Effect

abstract: Recent demonstrations of electrical detection and manipulation of antiferromagnets (AFMs) have opened new opportunities towards robust and ultrafast spintronics devices. However, it is difficult to establish the connection between the spin-transport behavior and the microscopic AFM domain states due to the lack of the real-time AFM domain imaging technique under the electric field. Here we report a significant Voigt rotation up to 60 mdeg in thin NiO(001) films at room temperature. Such large Voigt rotation allows us to directly observe AFM domains in thin-film NiO by utilizing a wide-field optical microscope. Further complementary XMLD-PEEM measurement confirms that the Voigt contrast originates from the NiO AFM order. We examine the domain pattern evolution at a wide range of temperature and with the application of external magnetic field. Comparing to large-scale-facility techniques such as the X-ray photoemission electron microscopy, the use with a wide-field, tabletop optical imaging method enables straightforward access to domain configurations of single-layer AFMs.


title: The cos$\phi_R$ asymmetry of dihadron production in double longitudinally polarized SIDIS

abstract: We present a study on the double longitudinal-spin asymmetry of dihadron production in semi-inclusive deep inelastic scattering, in which the total transverse momentum of the final-state hadron pairs is integrated out. In particular, we investigate the origin of the $\cos\phi_R$ azimuthal asymmetry for which we take into account the coupling of the helicity distribution $g_1$ and the twist-3 dihadron fragmentation function $\widetilde{D}^{\sphericalangle}$. We calculate the $s$-wave and $p$-wave interference term $\widetilde{D}^{\sphericalangle}_{ot}$ in a spectator model. We estimate the $\cos\phi_R$ asymmetry at the kinematics of COMPASS which is collecting data on dihadron production in polarized deep inelastic scattering. The prediction of the same asymmetry at JLab 12GeV and a future EIC are also presented. Our study indicates that measuring the $\cos\phi_R$ asymmetry in SIDIS may be a ideal way to probe the dihadron fragmentation function $\widetilde{D}^{\sphericalangle}$.


title: Exploring the properties of warm and cold atomic hydrogen in the Taurus and Gemini regions

abstract: We report Arecibo 21 cm absorption-emission observations to characterise the physical properties of neutral hydrogen (HI) in the proximity of five giant molecular clouds (GMCs): Taurus, California, Rosette, Mon OB1, and NGC 2264. Strong HI absorption was detected toward all 79 background continuum sources in the ~60x20 square degree region. Gaussian decompositions were performed to estimate temperatures, optical depths and column densities of the cold and warm neutral medium (CNM, WNM). The properties of individual CNM components are similar to those previously observed along random Galactic sightlines and in the vicinity of GMCs, suggesting a universality of cold HI properties. The CNM spin temperature (Ts) histogram peaks at ~50K. The turbulent Mach numbers of CNM vary widely, with a typical value of ~4, indicating that their motions are supersonic. About 60% of the total HI gas is WNM, and nearly 40% of the WNM lies in thermally unstable regime 500-5000K. The observed CNM fraction is higher around GMCs than in diffuse regions, and increases with increasing column density (NHI) to a maximum of ~75%. On average, the optically thin approximation (N*(HI)) underestimates the total N(HI) by ~21%, but we find large regional differences in the relationship between N(HI) and the required correction factor, f=N(HI)/N*(HI). We examine two different methods (linear fit of f vs log10(N*(HI)) and uniform Ts) to correct for opacity effects using emission data from the GALFA-HI survey. We prefer the uniform Ts method, since the linear relationship does not produce convincing fits for all subregions.


title: Structured Pruning of Recurrent Neural Networks through Neuron Selection

abstract: Recurrent neural networks (RNNs) have recently achieved remarkable successes in a number of applications. However, the huge sizes and computational burden of these models make it difficult for their deployment on edge devices. A practically effective approach is to reduce the overall storage and computation costs of RNNs by network pruning techniques. Despite their successful applications, those pruning methods based on Lasso either produce irregular sparse patterns in weight matrices, which is not helpful in practical speedup. To address these issues, we propose structured pruning method through neuron selection which can reduce the sizes of basic structures of RNNs. More specifically, we introduce two sets of binary random variables, which can be interpreted as gates or switches to the input neurons and the hidden neurons, respectively. We demonstrate that the corresponding optimization problem can be addressed by minimizing the L0 norm of the weight matrix. Finally, experimental results on language modeling and machine reading comprehension tasks have indicated the advantages of the proposed method in comparison with state-of-the-art pruning competitors. In particular, nearly 20 x practical speedup during inference was achieved without losing performance for language model on the Penn TreeBank dataset, indicating the promising performance of the proposed method


title: Well-posedness of compressible magneto-micropolar fluid equations

abstract: We are concerned with compressible magneto-micropolar fluid equations (1.1)-(1.2). The global existence and large time behaviour of solutions near a constant state to the magneto-micropolar-Navier-Stokes-Poisson (MMNSP) system is investigated in $\mathbb{R}^3$. By a refined energy method, the global existence is established under the assumption that the $H^3$ norm of the initial data is small, but the higher order derivatives can be large. If the initial data belongs to homogeneous Sobolev spaces or homogeneous Besov spaces, we prove the optimal time decay rates of the solution and its higher order spatial derivatives. Meanwhile, we also obtain the usual $L^p-L^2$ $(1\leq p\leq2)$ type of the decay rates without requiring that the $L^p$ norm of initial data is small.


title: A Multi-Task Architecture on Relevance-based Neural Query Translation

abstract: We describe a multi-task learning approach to train a Neural Machine Translation (NMT) model with a Relevance-based Auxiliary Task (RAT) for search query translation. The translation process for Cross-lingual Information Retrieval (CLIR) task is usually treated as a black box and it is performed as an independent step. However, an NMT model trained on sentence-level parallel data is not aware of the vocabulary distribution of the retrieval corpus. We address this problem with our multi-task learning architecture that achieves 16% improvement over a strong NMT baseline on Italian-English query-document dataset. We show using both quantitative and qualitative analysis that our model generates balanced and precise translations with the regularization effect it achieves from multi-task learning paradigm.


title: Non-vanishing of tidal Love numbers

abstract: In this note, we analyze black holes solutions under R3 corrections, which is the leading correction induced by quantum corrections in four-dimensions. We showed that perturbations around this black hole background will lead to non-zero tidal Love number (TLN). This further accentuates the \unnaturalness" of the vanishing TLN for Schwarzschild black hole under Einstein-Hilbert action.


title: Near Optimal Coflow Scheduling in Networks

abstract: The coflow scheduling problem has emerged as a popular abstraction in the last few years to study data communication problems within a data center. In this basic framework, each coflow has a set of communication demands and the goal is to schedule many coflows in a manner that minimizes the total weighted completion time. A coflow is said to complete when all its communication needs are met. This problem has been extremely well studied for the case of complete bipartite graphs that model a data center with full bisection bandwidth and several approximation algorithms and effective heuristics have been proposed recently. In this work, we study a slightly different model of coflow scheduling in general graphs (to capture traffic between data centers) and develop practical and efficient approximation algorithms for it. Our main result is a randomized 2 approximation algorithm for the single path and free path model, significantly improving prior work. In addition, we demonstrate via extensive experiments that the algorithm is practical, easy to implement and performs well in practice.