title: Unsupervised Learning for Real-World Super-Resolution

abstract: Most current super-resolution methods rely on low and high resolution image pairs to train a network in a fully supervised manner. However, such image pairs are not available in real-world applications. Instead of directly addressing this problem, most works employ the popular bicubic downsampling strategy to artificially generate a corresponding low resolution image. Unfortunately, this strategy introduces significant artifacts, removing natural sensor noise and other real-world characteristics. Super-resolution networks trained on such bicubic images therefore struggle to generalize to natural images. In this work, we propose an unsupervised approach for image super-resolution. Given only unpaired data, we learn to invert the effects of bicubic downsampling in order to restore the natural image characteristics present in the data. This allows us to generate realistic image pairs, faithfully reflecting the distribution of real-world images. Our super-resolution network can therefore be trained with direct pixel-wise supervision in the high resolution domain, while robustly generalizing to real input. We demonstrate the effectiveness of our approach in quantitative and qualitative experiments.


title: Manipulation Attacks in Local Differential Privacy

abstract: Local differential privacy is a widely studied restriction on distributed algorithms that collect aggregates about sensitive user data, and is now deployed in several large systems. We initiate a systematic study of a fundamental limitation of locally differentially private protocols: they are highly vulnerable to adversarial manipulation. While any algorithm can be manipulated by adversaries who lie about their inputs, we show that any non-interactive locally differentially private protocol can be manipulated to a much greater extent. Namely, when the privacy level is high or the input domain is large, an attacker who controls a small fraction of the users in the protocol can completely obscure the distribution of the users' inputs. We also show that existing protocols differ greatly in their resistance to manipulation, even when they offer the same accuracy guarantee with honest execution. Our results suggest caution when deploying local differential privacy and reinforce the importance of efficient cryptographic techniques for emulating mechanisms from central differential privacy in distributed settings.


title: POD-Galerkin Model Order Reduction for Parametrized Time Dependent Linear Quadratic Optimal Control Problems in Saddle Point Formulation

abstract: In this work we recast parametrized time dependent optimal control problems governed by partial differential equations in a saddle point formulation and we propose reduced order methods as an effective strategy to solve them. Indeed, on one hand parametrized time dependent optimal control is a very powerful mathematical model which is able to describe several physical phenomena; on the other hand, it requires a huge computational effort. Reduced order methods are a suitable approach to have rapid and accurate simulations. We rely on POD-Galerkin reduction over the physical and geometrical parameters of the optimality system in a space-time formulation. Our theoretical results and our methodology are tested on two examples: a boundary time dependent optimal control for a Graetz flow and a distributed optimal control governed by time dependent Stokes equations. With these two experiments the convenience of the reduced order modelling is further extended to the field of time dependent optimal control.


title: Photometric Redshift Calibration with Self Organising Maps

abstract: Accurate photometric redshift calibration is central to the robustness of all cosmology constraints from cosmic shear surveys. Analyses of the KiDS re-weighted training samples from all overlapping spectroscopic surveys to provide a direct redshift calibration. Using self-organising maps (SOMs) we demonstrate that this spectroscopic compilation is sufficiently complete for KiDS, representing $99\%$ of the effective 2D cosmic shear sample. We use the SOM to define a $100\%$ represented `gold' cosmic shear sample, per tomographic bin. Using mock simulations of KiDS and the spectroscopic training set, we estimate the uncertainty on the SOM redshift calibration, and find that photometric noise, sample variance, and spectroscopic selection effects (including redshift and magnitude incompleteness) induce a combined maximal scatter on the bias of the redshift distribution reconstruction ($\Delta \langle z \rangle=\langle z \rangle_{\rm est}-\langle z \rangle_{\rm true}$) of $\sigma_{\Delta \langle z \rangle} \leq 0.006$ in all tomographic bins. We show that the SOM calibration is unbiased in the cases of noiseless photometry and perfectly representative spectroscopic datasets, as expected from theory. The inclusion of both photometric noise and spectroscopic selection effects in our mock data introduces a maximal bias of $\Delta \langle z \rangle =0.013\pm0.006$, or $\Delta \langle z \rangle \leq 0.025$ at $97.5\%$ confidence, once quality flags have been applied to the SOM. The method presented here represents a significant improvement over the previously adopted direct redshift calibration implementation for KiDS, owing to its diagnostic and quality assurance capabilities. The implementation of this method in future cosmic shear studies will allow better diagnosis, examination, and mitigation of systematic biases in photometric redshift calibration.


title: Quantifying the Impact of Cognitive Biases in Question-Answering Systems

abstract: Crowdsourcing can identify high-quality solutions to problems; however, individual decisions are constrained by cognitive biases. We investigate some of these biases in an experimental model of a question-answering system. In both natural and controlled experiments, we observe a strong position bias in favor of answers appearing earlier in a list of choices. This effect is enhanced by three cognitive factors: the attention an answer receives, its perceived popularity, and cognitive load, measured by the number of choices a user has to process. While separately weak, these effects synergistically amplify position bias and decouple user choices of best answers from their intrinsic quality. We end our paper by discussing the novel ways we can apply these findings to substantially improve how high-quality answers are found in question-answering systems.


title: Refined symmetry indicators for topological superconductors in all space groups

abstract: Topological superconductors are exotic phases of matter featuring robust surface states that could be leveraged for topological quantum computation. A useful guiding principle for the search of topological superconductors is to relate the topological invariants with the behavior of the pairing order parameter on the normal-state Fermi surfaces. The existing formulas, however, become inadequate for the prediction of the recently proposed classes of topological crystalline superconductors. In this work, we advance the theory of symmetry indicators for topological (crystalline) superconductors to cover all space groups. Our main result is the exhaustive computation of the indicator groups for superconductors under a variety of symmetry settings. We further illustrate the power of this approach by analyzing four-fold symmetric superconductors with or without inversion symmetry, and show that the indicators can diagnose topological superconductors with surface states of different dimensionalities or dictate gaplessness in the bulk excitation spectrum.


title: On the classification of 5d SCFTs

abstract: We determine all 5d SCFTs upto rank three by studying RG flows of 5d KK theories. Our analysis reveals the existence of new rank one and rank two 5d SCFTs not captured by previous classifications. In addition to that, we provide for the first time a systematic and conjecturally complete classification of rank three 5d SCFTs. Our methods are based on a recently studied geometric description of 5d KK theories, thus demonstrating the utility of these geometric descriptions. It is straightforward, though computationally intensive, to extend this work and systematically classify 5d SCFTs of higher ranks (greater than or equal to four) by using the geometric description of 5d KK theories.


title: 1/4 BPS AdS$_3$/CFT$_2$

abstract: We discuss new solutions in massive Type IIA supergravity with AdS$_3\times$S$^2$ factors, preserving ${\cal N}=(0,4)$ SUSY. We propose a duality with a precise family of quivers that flow to ${\cal N}=(0,4)$ fixed points at low energies. These quivers consist on two families of linear quivers coupled by matter fields. Physical observables such as the central charges provide stringent checks of the proposed duality. A formal mapping is presented connecting our backgrounds with those dual to six dimensional ${\cal N}=(1,0)$ CFTs, suggesting the existence of a flow across dimensions between the CFTs.


title: The Conformal Equivalence of the half string 3-vertex and the Sciuto,Caneschi-schwimmer-veneziano 3-vertex

abstract: In this paper we construct the operator connecting the three vertex in half string field theory and the dual model vertex of Sciuto, Caneshi Schwimmer and Veneziano (SCSV). This construction results in an explicit conformal transformation linking the two interactions in the matter sector at all levels. Thus establishing the equivalence in the matter sector between the two theories at least for N = 3.


title: Accident Risk Prediction based on Heterogeneous Sparse Data: New Dataset and Insights

abstract: Reducing traffic accidents is an important public safety challenge, therefore, accident analysis and prediction has been a topic of much research over the past few decades. Using small-scale datasets with limited coverage, being dependent on extensive set of data, and being not applicable for real-time purposes are the important shortcomings of the existing studies. To address these challenges, we propose a new solution for real-time traffic accident prediction using easy-to-obtain, but sparse data. Our solution relies on a deep-neural-network model (which we have named DAP, for Deep Accident Prediction); which utilizes a variety of data attributes such as traffic events, weather data, points-of-interest, and time. DAP incorporates multiple components including a recurrent (for time-sensitive data), a fully connected (for time-insensitive data), and a trainable embedding component (to capture spatial heterogeneity). To fill the data gap, we have - through a comprehensive process of data collection, integration, and augmentation - created a large-scale publicly available database of accident information named US-Accidents. By employing the US-Accidents dataset and through an extensive set of experiments across several large cities, we have evaluated our proposal against several baselines. Our analysis and results show significant improvements to predict rare accident events. Further, we have shown the impact of traffic information, time, and points-of-interest data for real-time accident prediction.


title: A new class of galaxy models with a central BH -- I. The spherical case

abstract: The dynamical properties of spherically symmetric galaxy models, where a Jaffe (1983) stellar density profile is embedded in a total mass density decreasing as $r^{-3}$ at large radii, are presented. The orbital structure of the stellar component is described by the Osipkov--Merritt anisotropy; the dark matter halo is isotropic, and a black hole is added at the center of the galaxy. First, the conditions for a nowhere negative and monotonically decreasing dark matter halo density profile are derived; this profile can be made asymptotically coincident with a NFW profile at the center and at large radii. Then the minimum value of the anisotropy radius for phase-space consistency is derived as a function of the galaxy parameters. The Jeans equations for the stellar component are solved analytically; the projected velocity dispersion at the center and at large radii is also obtained, for generic values of the anisotropy radius. Finally, analytical expressions for the terms entering the Virial Theorem are derived, and the fiducial anisotropy limit required to prevent the onset of Radial Orbit Instability is determined as a function of the galaxy parameters. The presented models, built following an approach already adopted in our previous works, can be a useful starting point for a more advanced modeling of the dynamics of elliptical galaxies, and can be easily implemented in numerical simulations requiring a realistic dynamical model of a galaxy.


title: Beam-energy dependence of the azimuthal anisotropic flow from RHIC

abstract: Recent STAR measurements of the anisotropic flow coefficients (v$_{n}$) are presented for Au+Au collisions spanning the beam energy range $\sqrt{s_{NN}} = 7.7 - 200$ GeV. The measurements indicate dependences on the harmonic number ($n$), transverse momentum ($p_T$), pseudorapidity ($\eta$), collision centrality and beam energy ($\sqrt{s_{NN}}$) which could serve as important constraints to test different initial-state models and to aid precision extraction of the temperature dependence of the specific shear viscosity.


title: Productivity propagation with networks transformation

abstract: We model sectoral production by cascading binary compounding processes. The sequence of processes is discovered in a self-similar hierarchical structure stylized in the economy-wide networks of production. Nested substitution elasticities and Hicks-neutral productivity growth are measured such that the general equilibrium feedbacks between all sectoral unit cost functions replicate the transformation of networks observed as a set of two temporally distant input-output coefficient matrices. We examine this system of unit cost functions to determine how idiosyncratic sectoral productivity shocks propagate into aggregate macroeconomic fluctuations in light of potential network transformation. Additionally, we study how sectoral productivity increments propagate into the dynamic general equilibrium, thereby allowing network transformation and ultimately producing social benefits.


title: Zeros of primitive characters of finite groups

abstract: We classify finite non-solvable groups with a faithful primitive irreducible complex character that vanishes on a unique conjugacy class. Our results answer a question of Dixon and Rahnamai Barghi and suggest an extension of Burnside's classical theorem on zeros of characters.


title: Connected Baranyai's Theorem

abstract: Let $K_n^h=(V,\binom{V}{h})$ be the complete $h$-uniform hypergraph on vertex set $V$ with $|V|=n$. Baranyai showed that $K_n^h$ can be expressed as the union of edge-disjoint $r$-regular factors if and only if $h$ divides $rn$ and $r$ divides $\binom{n-1}{h-1}$. Using a new proof technique, in this paper we prove that $\lambda K_n^h$ can be expressed as the union $\mathcal G_1\cup \ldots \cup\mathcal G_k$ of $k$ edge-disjoint factors, where for $1\leq i\leq k$, $\mathcal G_i$ is $r_i$-regular, if and only if (i) $h$ divides $r_in$ for $1\leq i\leq k$, and (ii) $\sum_{i=1}^k r_i=\lambda \binom{n-1}{h-1}$. Moreover, for any $i$ ($1\leq i\leq k$) for which $r_i\geq 2$, this new technique allows us to guarantee that $\mathcal G_i$ is connected, generalizing Baranyai's theorem, and answering a question by Katona.


title: Decomposing complete equipartite multigraphs into cycles of variable lengths: the amalgamation-detachment approach

abstract: Using the technique of amalgamation-detachment, we show that the complete equipartite multigraph $\lambda K_{n\times m}$ can be decomposed into cycles of lengths $c_1m,\dots,c_km$ (plus a 1-factor if the degree is odd) whenever there exists a decomposition of $\lambda m K_n$ into cycles of lengths $c_1, \dots,c_k$ (plus a 1-factor if the degree is odd). In addition, we give sufficient conditions for the existence of some other, related cycle decompositions of the complete equipartite multigraph $\lambda K_{n\times m}$.


title: The Effect of Star--Disk Interactions on Highly Eccentric Stellar Orbits in Active Galactic Nuclei: A Disk Loss Cone and Implications for Stellar Tidal Disruption Events

abstract: Galactic center black holes appear to be nearly universally surrounded by dense stellar clusters. When these black holes go through an active accretion phase, the multiple components of the accretion disk, stellar cluster, and black hole system all coexist. We analyze the effect of drag forces on highly eccentric stellar orbits incurred as stars puncture through the disk plane. Disk crossings dissipate orbital energy, drawing eccentric stars into more circular orbits. For high surface density disks, such as those found around black holes accreting near the Eddington mass accretion limit, the magnitude of this energy dissipation can be larger than the mean scatterings that stars receive by two body relaxation. One implication of this is the presence of a disk "loss cone" for highly eccentric stellar orbits where the dissipation from disk interaction outweighs scatter via two body relaxation. The disk loss cone is larger than the tidal disruption loss cone for near-Eddington black hole accretion rates. Stellar orbits within the disk loss cone are lost from the overall cluster as stellar orbits are circularized and stars are potentially ablated by their high-velocity impacts with the disk. We find, however, that the presence of the disk loss cone has a minimal effect on the overall rate of stellar tidal disruptions. Stars are still efficiently fed to the black hole from more-distant stellar orbits that receive large-enough per orbit scatter to jump over the disk loss cone and end up tidally disrupted.


title: Topological semimetallic phase in PbO$_2$ promoted by temperature

abstract: Materials exhibiting topological order host exotic phenomena that could form the basis for novel developments in areas ranging from low-power electronics to quantum computers. The past decade has witnessed multiple experimental realizations and thousands of predictions of topological materials. However, it has been determined that increasing temperature destroys topological order, restricting many topological materials to very low temperatures and thus hampering practical applications. Here, we propose a material realization of temperature promoted topological order. We show that a semiconducting oxide that has been widely used in lead-acid batteries, $\beta$-PbO$_2$, hosts a topological semimetallic phase driven by both thermal expansion and electron-phonon coupling upon increasing temperature. We identify the interplay between the quasi-two-dimensional nature of the charge distribution of the valence band with the three-dimensional nature of the charge distribution of the conduction band as the microscopic mechanism driving this unconventional temperature dependence. Thus, we propose a general principle to search for and design novel topological materials whose topological order is stabilized by increasing temperature. This provides a clear roadmap for taking topological materials from the laboratory to technological devices.


title: Small-Scale HI Channel Map Structure is Cold: Evidence from Na I Absorption at High Galactic Latitudes

abstract: The spatial distribution of neutral hydrogen (H I) emission is a powerful probe of interstellar medium physics. The small-scale structure in H I channel maps is often assumed to probe the velocity field rather than real density structures. In this work we directly test this assumption, using high-resolution GALFA-H I observations and 50,985 quasar spectra from SDSS. We measure the equivalent widths of interstellar Na I D$_1$ and Na I D$_2$ absorption, and robustly conclude that they depend more strongly on the column density of small-scale structure in H I than on either the large-scale H I structure or the total H I column. This is inconsistent with the hypothesis that small-scale channel map structure is driven by velocity crowding. Instead, the data favor the interpretation that this emission structure predominantly originates in cold, dense interstellar material, consistent with a clumpy cold neutral medium.


title: Detecting the neutral IGM in filaments with the SKA

abstract: The intergalactic medium (IGM) plays an important role in the formation and evolution of galaxies. Recent developments in upcoming radio telescopes are starting to open up the possibility of making a first direct detection of the 21 cm signal of neutral hydrogen (HI) from the warm gas of the IGM in large-scale filaments. The cosmological hydrodynamical EAGLE simulation is used to estimate the typical IGM filament signal. Assuming the same average signal for all filaments, a prediction is made for the detectability of such a signal with the upcoming mid-frequency array of the Square Kilometer Array (SKA1-mid) or the future upgrade to SKA2. The signal-to-noise (S/N) then only depends on the size and orientation of each filament. With filament spines inferred from existing galaxy surveys as a proxy for typical real filaments, we find hundreds of filaments in the region of the sky accessible to the SKA that can be detected. Once the various phases of the SKA telescope become operational, their own surveys will be able to find the galaxies required to infer the position of even more filaments within the survey area. We find that in 120 h, SKA1-mid/SKA2 will detect HI emission from the strongest filaments in the field with a S/N of the order of 10 to $\sim$150 for the most pessimistic model considered here. Some of the brighter filaments can be detected with an integration time of a few minutes with SKA1-mid and a few seconds with SKA2. Therefore, SKA2 will be capable of not only detecting but also mapping a large part of the IGM in these filaments.


title: Hybrid kernel polynomial method

abstract: The kernel polynomial method allows to sample overall spectral properties of a quantum system, while sparse diagonalization provides accurate information about a few important states. We present a method combining these two approaches without loss of performance or accuracy. We apply this hybrid kernel polynomial method to improve the computation of thermodynamic quantities and the construction of perturbative effective models, in a regime where neither of the methods is sufficient on its own. To achieve this we develop a perturbative kernel polynomial method to compute arbitrary order series expansions of expectation values. We demonstrate the efficiency of our approach on three examples: the calculation of supercurrent and inductance in a Josephson junction, the interaction of spin qubits defined in a two dimensional electron gas, and the calculation of the effective band structure in a realistic model of a semiconductor nanowire.


title: Smoothed particle magnetohydrodynamics with the geometric density average force expression

abstract: We present a novel method of magnetohydrodynamics (MHD) within the smoothed particle hydrodynamics scheme using the Geometric Density average force expression (GDSPH). GDSPH has recently been shown to reduce the leading order errors and greatly improve the accuracy near density discontinuities, eliminating surface tension effects. Here, we extend the study to investigate how SPMHD benefits from this method. We implement ideal MHD in the Gasoline2 and Changa codes with both GDSPH and traditional SPH (TSPH) schemes. An constrained hyperbolic divergence cleaning scheme is employed to control the divergence error, and a switch for artificial resistivity with minimized dissipation is used. We test the codes with a large suite of MHD tests, and show that in all problems the results are comparable or improved over previous SPMHD implementations. While both GDSPH and TSPH perform well with relatively smooth or highly supersonic flows, GDSPH shows significant improvements in the presence of strong discontinuities and large dynamic scales. In particular, when applied to an astrophysical problem of the collapse of a magnetized cloud, GDSPH realistically captures the development of a magnetic tower and jet launching in the weak-field regime, and exhibit fast convergence with resolution, while TSPH failed to do so. Our new method shows qualitatively similar results to the ones from the meshless finite mass/volume (MFM/MFV) schemes within the Gizmo code, while remaining computationally less expensive.


title: WARPFIELD-EMP: The Self-Consistent Prediction of Emission Lines from Evolving HII Regions in Dense Molecular Clouds

abstract: We present the {\sc warpfield} emission predictor, {\sc warpfield-emp}, which couples the 1D stellar feedback code {\sc warpfield} with the {\sc cloudy} \hii region/PDR code and the {\sc polaris} radiative transfer code, in order to make detailed predictions for the time-dependent line and continuum emission arising from the H{\sc ii} region and PDR surrounding an evolving star cluster. {\sc warpfield-emp} accounts for a wide range of physical processes (stellar winds, supernovae, radiation pressure, gravity, thermal conduction, radiative cooling, dust extinction etc.) and yet runs quickly enough to allow us to explore broad ranges of different cloud parameters. We compare the results of an extensive set of models with SITELLE observations of a large sample of \hii regions in NGC~628 and find very good agreement, particularly for the highest signal-to-noise observations. We show that our approach of modeling individual clouds from first principles (instead of in terms of dimensionless quantities such as the ionization parameter) allows us to avoid long-standing degeneracies in the interpretation of \hii region diagnostics and enables us to relate these diagnostics to important physical parameters such as cloud mass or cluster age. Finally, we explore the implications of our models regarding the reliability of simple metallicity diagnostics, the properties of long-lived embedded clusters, and the role played by winds and supernovae in regulating \hii region and PDR line emission.


title: Do Rydberg chains yield Fibonacci anyons?

abstract: Recent experiments have focused attention on the properties of chains of atoms in which the atoms are either in their ground states or in highly excited Rydberg states which block similar excitations in their immediate neighbors. As the low energy Hilbert space of such chains is isomorphic to that of a chain of Fibonacci anyons, they have been proposed as a platform for topological quantum computation and for simulating anyon dynamics. We show that generic local operators in the Rydberg chain correspond to non-local anyonic operators that do not preserve a topological symmetry of the physical anyons. Consequently, we argue that Rydberg chains do not yield Fibonacci anyons and quantum computation with Rydberg atoms is not topologically protected.


title: Tuning entanglement by squeezing magnons in anisotropic magnets

abstract: We theoretically study the entanglement between two arbitrary spins in a magnetic material, where magnons naturally form a general squeezed coherent state, in the presence of an applied magnetic field and axial anisotropies. Employing concurrence as a measure of entanglement, we demonstrate that spins are generally entangled in thermodynamic equilibrium, with the amount of entanglement controlled by the external fields and anisotropies. As a result, the magnetic medium can serve as a resource to store and process quantum information. We, furthermore, show that the entanglement can jump discontinuously when decreasing the applied magnetic field. This tunable entanglement can be potentially used as an efficient switch in quantum-information processing tasks.


title: Entanglement and matrix elements of observables in interacting integrable systems

abstract: We study the bipartite von Neumann entanglement entropy and matrix elements of local operators in the eigenstates of an interacting integrable Hamiltonian (the paradigmatic spin-1/2 XXZ chain), and we contrast their behavior with that of quantum chaotic systems. We find that the leading term of the average (over all eigenstates in the zero magnetization sector) eigenstate entanglement entropy has a volume-law coefficient that is smaller than the universal (maximal entanglement) one in quantum chaotic systems. This establishes the entanglement entropy as a powerful measure to distinguish integrable models from generic ones. Remarkably, our numerical results suggest that the volume-law coefficient of the average entanglement entropy of eigenstates of the spin-1/2 XXZ Hamiltonian is very close to, or the same as, the one for translationally invariant quadratic fermionic models. We also study matrix elements of local operators in the eigenstates of the spin-1/2 XXZ Hamiltonian at the center of the spectrum. For the diagonal matrix elements, we show evidence that the support does not vanish with increasing system size, while the average eigenstate-to-eigenstate fluctuations vanish in a power-law fashion. For the off-diagonal matrix elements, we show that they follow a distribution that is close to (but not quite) log-normal, and that their variance is a well-defined function of $\omega=E_{\alpha}-E_{\beta}$ ($\{E_{\alpha}\}$ are the eigenenergies) proportional to $1/D$, where $D$ is the Hilbert space dimension.


title: Circumplanetary Disk Dynamics in the Isothermal and Adiabatic Limits

abstract: Circumplanetary disks (CPDs) may be essential to the formation of planets, regulating their spin and accretion evolution. We perform a series of 3D hydrodynamics simulations in both the isothermal and adiabatic limits to systematically measure the rotation rates, sizes, and masses of CPDs as functions of $q_{\rm thermal}$, the ratio of the planet mass to the disk thermal mass. Our $q_{\rm thermal}$ ranges from 0.1 to 4; for our various disk temperatures, this corresponds to planet masses between 1 Earth mass and 4 Jupiter masses. Within this parameter space, we find that isothermal CPDs are disky and bound within $\sim$10\% of the planet's Bondi radius $r_{\rm B}$, with the innermost $\sim0.05\,r_{\rm B}$ in full rotational support. Adiabatic CPDs are spherical (therefore not actually "disks"), bound within $\sim0.2\,r_{\rm B}$, and mainly pressure-supported with rotation rates scaling linearly with $q_{\rm thermal}$; extrapolation suggests full rotational support of adiabatic envelopes at $\sim10\,q_{\rm thermal}$. Fast rotation and 3D supersonic flow render isothermal CPDs significantly different in structure from --- and orders of magnitude less massive than --- their 1D isothermal hydrostatic counterparts. Inside a minimum-mass solar nebula, even a maximally cooled, isothermal CPD around a 10 Earth-mass core may have less than 1 Earth mass, suggesting that gas giant formation may hinge on angular momentum transport processes in CPDs. Our CPD sizes and masses appear consistent with the regular satellites orbiting solar system giants.


title: Understanding and Robustifying Differentiable Architecture Search

abstract: Differentiable Architecture Search (DARTS) has attracted a lot of attention due to its simplicity and small search costs achieved by a continuous relaxation and an approximation of the resulting bi-level optimization problem. However, DARTS does not work robustly for new problems: we identify a wide range of search spaces for which DARTS yields degenerate architectures with very poor test performance. We study this failure mode and show that, while DARTS successfully minimizes validation loss, the found solutions generalize poorly when they coincide with high validation loss curvature in the architecture space. We show that by adding one of various types of regularization we can robustify DARTS to find solutions with less curvature and better generalization properties. Based on these observations, we propose several simple variations of DARTS that perform substantially more robustly in practice. Our observations are robust across five search spaces on three image classification tasks and also hold for the very different domains of disparity estimation (a dense regression task) and language modelling.


title: The Information Geometry of Space-time

abstract: The method of maximum entropy is used to model curved physical space in terms of points defined with a finite resolution. Such a blurred space is automatically endowed with a metric given by information geometry. The corresponding space-time is such that the geometry of any embedded spacelike surface is given by its information geometry. The dynamics of blurred space, its geometrodynamics, is constructed by requiring that as space undergoes the deformations associated with evolution in local time, it sweeps a four-dimensional space-time. This reproduces Einstein's equations for vacuum gravity. We conclude with brief comments on some of the peculiar properties of blurred space: There is a minimum length and blurred points have a finite volume. There is a relativistic "blur dilation". The volume of space is a measure of its entropy.


title: Birational and noncommutative lifts of antichain toggling and rowmotion

abstract: The rowmotion action on order ideals or on antichains of a finite partially ordered set has been studied (under a variety of names) by many authors. Depending on the poset, one finds unexpectedly interesting orbit structures, instances of (small order) periodicity, cyclic sieving, and homomesy. Many of these nice features still hold when the action is extended to $[0,1]$-labelings of the poset or (via detropicalization) to labelings by rational functions (the birational setting). In this work, we parallel the birational lifting already done for order-ideal rowmotion to antichain rowmotion. We give explicit equivariant bijections between the birational toggle groups and between their respective liftings. We further extend all of these notions to labellings by noncommutative rational functions, setting an unpublished periodicity conjecture of Grinberg in a broader context.


title: Orientational ordering and layering of hard plates in narrow slit-like pores

abstract: We examine the ordering behavior of hard plate-like particle in a very narrow slit-like pore using the Parsons-Lee density functional theory and the restricted orientation approximation. We observe that the plates are orientationally ordered and align perpendicularly (face-on) to the walls at low densities, a first order layering transition occurs between uniaxial nematic structures having n and n+1 layers at intermediate densities and even a phase transition between a monolayer with parallel (edge-on) orientational order and n layers with perpendicular one can be detected at high densities. In addition to this, the edge-on monolayer is usually biaxial nematic and a uniaxial-biaxial nematic phase transition can be also seen at very high densities.


title: Constraining level densities through quantitative correlations with cross-section data

abstract: The adopted level densities (LD) for the nuclei produced through different reaction mechanisms significantly impact the calculation of cross sections for the many reaction channels. Common LD models make simplified assumptions regarding the overall behavior of the total LD and the intrinsic spin and parity distributions of the excited states. However, very few experimental constraints are taken into account: LD at neutron separation energy coming from average resonance spacings, whenever they have been previously measured, and the sometimes subjective extrapolation of discrete levels. These, however, constrain the LD only for very specific spins, parities and excitation energies. This work aims to establish additional experimental constraints on LD through quantitative correlations between cross sections and LD. This allows for the fitting and determination of detailed structures in LD. For this we use the microscopic Hartree-Fock-Bogoliubov (HFB) LD to associate variations predicted by the model with the structure observed in double-differential spectra at low outgoing neutron energy, which is dominated by the LD input. We also use \nuc{56}{Fe} ($n,p$) as an example cross sections are extremely sensitive to LD. For comparison purposes we also perform calculations with the GC model. With this approach we are able to perform fits of the LD based on actual experimental data, constraining the model and ensuring its consistency. This approach can be particularly useful in extrapolating the LD to nuclei for which high-excited discrete levels and/or resonance spacings are unknown. It also predicts inelastic gamma cross sections that can significantly differ from more standard phenomenological LD.


title: A Dwarf Disrupting -- Andromeda XXVII and the North West Stream

abstract: We present a kinematic and spectroscopic analysis of 38 red giant branch stars, in 7 fields, spanning the dwarf spheroidal galaxy Andromeda XXVII and the upper segment of the North West Stream. Both features are located in the outer halo of the Andromeda galaxy at a projected radius of 50-80 kpc, with the stream extending for $\sim$3$^{\circ}$ on the sky. Our data is obtained as part of the PAndAS survey and enables us to confirm that Andromeda XXVII's heliocentric distance is 827 $\pm$ 47 kpc and spectroscopic metallicity is -2.1$^{+0.4}_{-0.5}$. We also re-derive Andromeda XXVII's kinematic properties, measuring a systemic velocity = -526.1$^{+10.0}_{-11.0}${\kms} and a velocity dispersion that we find to be non-Gaussian but for which we derive a formal value of 27.0$^{+2.2}_{-3.9}${\kms}. In the upper segment of the North West Stream we measure mean values for the metallicity = -1.8$\pm$0.4, systemic velocity = -519.4 $\pm$4.0{\kms} and velocity dispersion = 10.0$\pm$4.0{\kms}. We also detect a velocity gradient of 1.7$\pm$0.3 {\kms} kpc$^{-1}$ on an infall trajectory towards M31. With a similar gradient, acting in the same direction, in the lower segment we suggest that the North West Stream is not a single structure. As the properties of the upper segment of the North West Stream and Andromeda XXVII are consistent within ~90\% confidence limits, it is likely that the two are related and plausible that Andromeda XXVII is the progenitor of this stream.


title: Mine Tunnel Exploration using Multiple Quadrupedal Robots

abstract: Robotic exploration of underground environments is a particularly challenging problem due to communication, endurance, and traversability constraints which necessitate high degrees of autonomy and agility. These challenges are further exacerbated by the need to minimize human intervention for practical applications. While legged robots have the ability to traverse extremely challenging terrain, they also engender new challenges for planning, estimation, and control. In this work, we describe a fully autonomous system for multi-robot mine exploration and mapping using legged quadrupeds, as well as a distributed database mesh networking system for reporting data. In addition, we show results from the DARPA Subterranean Challenge (SubT) Tunnel Circuit demonstrating localization of artifacts after traversals of hundreds of meters. These experiments describe fully autonomous exploration of an unknown Global Navigation Satellite System (GNSS)-denied environment undertaken by legged robots.


title: Exact Ground State of Lieb-Mattis Hamiltonian as a Superposition of N\'eel states

abstract: We show that the exact ground state of the Lieb-Mattis Hamiltonian is an equal-weight superposition of all possible classical N\'{e}el states, and provide an exact formulation of this superposition in the $z$-spin basis for both $S=1/2$ and general $S$ using Schwinger bosons. In general, a superposition of possible rotations on a general initial state is symmetric if and only if the initial state has a nonzero overlap with a singlet state and is otherwise made up of states that vanish due to the symmetrization. Most notably, $|s, m=0 \rangle$ states will vanish if symmetrized, which explains how a superposition of N\'{e}el states projects onto its singlet component.

categories: cond-mat.str-el cond-mat.other math-ph math.MP


title: Multidimensional quantum-enhanced target detection via spectro-temporal correlation measurements

abstract: In this work we investigate quantum-enhanced target detection in the presence of large background noise using multidimensional quantum correlations between photon pairs generated through spontaneous parametric down-conversion. Until now similar experiments have only utilized one of the photon pairs' many degrees of freedom such as temporal correlations and photon number correlations. Here, we utilized both temporal and spectral correlations of the photon pairs and achieved over an order of magnitude reduction to the background noise and in turn significant reduction to data acquisition time when compared to utilizing only temporal modes. We believe this work represents an important step in realizing a practical, real-time quantum-enhanced target detection system. The demonstrated technique will also be of importance in many other quantum sensing applications and quantum communications.


title: A note on additive twists, reciprocity laws and quantum modular forms

abstract: We prove that the central values of additive twists of a cuspidal $L$-function define a quantum modular form in the sense of Zagier, generalizing recent results of Bettin and Drappeau. From this we deduce a reciprocity law for the twisted first moment of multiplicative twists of cuspidal $L$-functions, similar to reciprocity laws discovered by Conrey for the twisted second moment of Dirichlet $L$-functions. Furthermore we give an interpretation of quantum modularity at infinity for additive twists of $L$-functions of weight 2 cusp forms in terms of the corresponding functional equations.


title: Bergman and Szego projections, Extremal Problems, and Square Functions

abstract: We study estimates for Hardy space norms of analytic projections. We first find a sufficient condition for the Bergman projection of a function in the unit disc to belong to the Hardy space $H^p$ for $1 < p < \infty$. We apply the result to prove a converse to an extension of Ryabykh's theorem about Hardy space regularity for Bergman space extremal functions. We also prove that the $H^q$ norm of the Szeg\"{o} projection of $F^{p/2} \overline{F}^{(p/2)-1}$ cannot be too small if $F$ is analytic, for certain values of $p$ and $q$. We apply this to show that the best analytic approximation in $L^p$ of a function in both $L^p$ and $L^q$ will also lie in $L^q$, for certain values of $p$ and $q$.


title: Online Hierarchical Clustering Approximations

abstract: Hierarchical clustering is a widely used approach for clustering datasets at multiple levels of granularity. Despite its popularity, existing algorithms such as hierarchical agglomerative clustering (HAC) are limited to the offline setting, and thus require the entire dataset to be available. This prohibits their use on large datasets commonly encountered in modern learning applications. In this paper, we consider hierarchical clustering in the online setting, where points arrive one at a time. We propose two algorithms that seek to optimize the Moseley and Wang (MW) revenue function, a variant of the Dasgupta cost. These algorithms offer different tradeoffs between efficiency and MW revenue performance. The first algorithm, OTD, is a highly efficient Online Top Down algorithm which provably achieves a 1/3-approximation to the MW revenue under a data separation assumption. The second algorithm, OHAC, is an online counterpart to offline HAC, which is known to yield a 1/3-approximation to the MW revenue, and produce good quality clusters in practice. We show that OHAC approximates offline HAC by leveraging a novel split-merge procedure. We empirically show that OTD and OHAC offer significant efficiency and cluster quality gains respectively over baselines.


title: Instability Zones in the Dynamics of a Quantum Mechanical Quasiperiodic Parametric Oscillator

abstract: Quasi-periodically driven quantum parametric oscillators have been the subject of several recent investigations. Here we show that for such oscillators, the instability zones of the mean position and variance (alternatively the mean energy) for a time developing wave packet are identical for the strongest resonance in the three-dimensional parameter space of the quasi-periodic modulation as it is for the two-dimensional parameter space of the periodic modulations.


title: Building a Library of Tactile Skills Based on FingerVision

abstract: Camera-based tactile sensors are emerging as a promising inexpensive solution for tactile-enhanced manipulation tasks. A recently introduced FingerVision sensor was shown capable of generating reliable signals for force estimation, object pose estimation, and slip detection. In this paper, we build upon the FingerVision design, improving already existing control algorithms, and, more importantly, expanding its range of applicability to more challenging tasks by utilizing raw skin deformation data for control. In contrast to previous approaches that rely on the average deformation of the whole sensor surface, we directly employ local deviations of each spherical marker immersed in the silicone body of the sensor for feedback control and as input to learning tasks. We show that with such input, substances of varying texture and viscosity can be distinguished on the basis of tactile sensations evoked while stirring them. As another application, we learn a mapping between skin deformation and force applied to an object. To demonstrate the full range of capabilities of the proposed controllers, we deploy them in a challenging architectural assembly task that involves inserting a load-bearing element underneath a bendable plate at the point of maximum load.


title: On the ATLAS Top Mass Measurements and the Potential for Stealth Stop Contamination

abstract: The discovery of the stop - the Supersymmetric partner of the top quark - is a key goal of the physics program enabled by the Large Hadron Collider. Although much of the accessible parameter space has already been probed, all current searches assume the top mass is known. This is relevant for the "stealth stop" regime, which is characterized by decay kinematics that force the final state top quark off its mass shell; such decays would contaminate the top mass measurements. We investigate the resulting bias imparted to the template method based ATLAS approach. A careful recasting of these results shows that effect can be as large as 2.0 GeV, comparable to the current quoted uncertainty on the top mass. Thus, a robust exploration of the stealth stop splinter requires the simultaneous consideration of the impact on the top mass. Additionally, we explore the robustness of the template technique, and point out a simple strategy for improving the methodology implemented for the semi-leptonic channel.


title: Angled crested like water waves with surface tension: Wellposedness of the problem

abstract: We consider the capillary-gravity water wave equation in two dimensions. We assume that the fluid is inviscid, incompressible, irrotational and the air density is zero. We construct an energy functional and prove a local wellposedness result without assuming the Taylor sign condition. When the surface tension $\sigma$ is zero, the energy reduces to a lower order version of the energy obtained by Kinsey and Wu [23] and allows angled crest interfaces. For positive surface tension, the energy does not allow angled crest interfaces but admits initial data with large curvature of the order of $\sigma^{-\frac{1}{3} + \epsilon} $ for any $\epsilon >0$.


title: Re-evaluating evidence for Hawking points in the CMB

abstract: We investigate recent claims for a detection of "Hawking points" (positions on the sky with unusually large temperature gradients between rings) in the cosmic microwave background (CMB) temperature maps at the 99.98% confidence level. We find that, after marginalization over the size of the rings, an excess is detected in Planck satellite maps at only an 87% confidence level (i.e., little more than 1$\sigma$). Therefore, we conclude that there is no statistically significant evidence for the presence of Hawking points in the CMB.


title: Perspectives on the distribution of orbits of distant Trans-Neptunian Objects

abstract: Looking at the orbits of small bodies with large semimajor axes, we are compelled to see patterns. Some of these patterns are noted as strong indicators of new or hidden processes in the outer Solar System, others are substantially generated by observational biases, and still others may be completely overlooked. We can gain insight into the current and past structure of the outer Solar System through a careful examination of these orbit patterns. In this chapter, we discuss the implications of the observed orbital distribution of distant trans-Neptunian objects (TNOs). We start with some cautions on how observational biases must affect the known set of TNO orbits. Some of these biases are intrinsic to the process of discovering TNOs, while others can be reduced or eliminated through careful observational survey design. We discuss some orbital element correlations that have received considerable attention in the recent literature. We examine the known TNOs in the context of the gravitational processes that the known Solar System induces in orbital distributions. We discuss proposed new elements of the outer Solar System, posited ancient processes, and the types of TNO orbital element distributions that they predict to exist. We conclude with speculation.


title: Controlling Assistive Robots with Learned Latent Actions

abstract: Assistive robotic arms enable users with physical disabilities to perform everyday tasks without relying on a caregiver. Unfortunately, the very dexterity that makes these arms useful also makes them challenging to teleoperate: the robot has more degrees-of-freedom than the human can directly coordinate with a handheld joystick. Our insight is that we can make assistive robots easier for humans to control by leveraging latent actions. Latent actions provide a low-dimensional embedding of high-dimensional robot behavior: for example, one latent dimension might guide the assistive arm along a pouring motion. In this paper, we design a teleoperation algorithm for assistive robots that learns latent actions from task demonstrations. We formulate the controllability, consistency, and scaling properties that user-friendly latent actions should have, and evaluate how different low-dimensional embeddings capture these properties. Finally, we conduct two user studies on a robotic arm to compare our latent action approach to both state-of-the-art shared autonomy baselines and a teleoperation strategy currently used by assistive arms. Participants completed assistive eating and cooking tasks more efficiently when leveraging our latent actions, and also subjectively reported that latent actions made the task easier to perform. The video accompanying this paper can be found at: https://youtu.be/wjnhrzugBj4.


title: Cross-Dataset Person Re-Identification via Unsupervised Pose Disentanglement and Adaptation

abstract: Person re-identification (re-ID) aims at recognizing the same person from images taken across different cameras. To address this challenging task, existing re-ID models typically rely on a large amount of labeled training data, which is not practical for real-world applications. To alleviate this limitation, researchers now targets at cross-dataset re-ID which focuses on generalizing the discriminative ability to the unlabeled target domain when given a labeled source domain dataset. To achieve this goal, our proposed Pose Disentanglement and Adaptation Network (PDA-Net) aims at learning deep image representation with pose and domain information properly disentangled. With the learned cross-domain pose invariant feature space, our proposed PDA-Net is able to perform pose disentanglement across domains without supervision in identities, and the resulting features can be applied to cross-dataset re-ID. Both of our qualitative and quantitative results on two benchmark datasets confirm the effectiveness of our approach and its superiority over the state-of-the-art cross-dataset Re-ID approaches.


title: Robust and high-resolution seismic complex trace analysis

abstract: Seismic attributes calculated by conventional methods are susceptible to noise. Conventional filtering reduces the noise in the cost of losing the spectral bandwidth. The challenge of having a high-resolution and robust signal processing tool motivated us to propose a sparse time-frequency decomposition while is stabilized for random noise. The procedure initiates by using Sparsity-based adaptive S-transform to regularize abrupt variations in frequency content of the nonstationary signals. Then, considering the fact that a higher amplitude of a frequency component results in a higher signal to noise ratio, an adaptive filter is applied to the time-frequency spectrum which is sparcified previously. The proposed zero adaptive filter enhances the high amplitude frequency components while suppresses the lower ones. The performance of the proposed method is compared to the sparse S-transform and the robust window Hilbert transform in estimation of instantaneous attributes by applying on synthetic and real data sets. Seismic attributes estimated by the proposed method is superior to the conventional ones in terms of its robustness and high resolution image. The proposed approach has a vast application in interpretation and identification of geological structures.


title: Gradual Network for Single Image De-raining

abstract: Most advances in single image de-raining meet a key challenge, which is removing rain streaks with different scales and shapes while preserving image details. Existing single image de-raining approaches treat rain-streak removal as a process of pixel-wise regression directly. However, they are lacking in mining the balance between over-de-raining (e.g. removing texture details in rain-free regions) and under-de-raining (e.g. leaving rain streaks). In this paper, we firstly propose a coarse-to-fine network called Gradual Network (GraNet) consisting of coarse stage and fine stage for delving into single image de-raining with different granularities. Specifically, to reveal coarse-grained rain-streak characteristics (e.g. long and thick rain streaks/raindrops), we propose a coarse stage by utilizing local-global spatial dependencies via a local-global subnetwork composed of region-aware blocks. Taking the residual result (the coarse de-rained result) between the rainy image sample (i.e. the input data) and the output of coarse stage (i.e. the learnt rain mask) as input, the fine stage continues to de-rain by removing the fine-grained rain streaks (e.g. light rain streaks and water mist) to get a rain-free and well-reconstructed output image via a unified contextual merging sub-network with dense blocks and a merging block. Solid and comprehensive experiments on synthetic and real data demonstrate that our GraNet can significantly outperform the state-of-the-art methods by removing rain streaks with various densities, scales and shapes while keeping the image details of rain-free regions well-preserved.


title: Sequential Dynamic Resource Allocation for Epidemic Control

abstract: Under the Dynamic Resource Allocation (DRA) model, an administrator has the mission to allocate dynamically a limited budget of resources to the nodes of a network in order to reduce a diffusion process (DP) (e.g. an epidemic). The standard DRA assumes that the administrator has constantly full information and instantaneous access to the entire network. Towards bringing such strategies closer to real-life constraints, we first present the Restricted DRA model extension where, at each intervention round, the access is restricted to only a fraction of the network nodes, called sample. Then, inspired by sequential selection problems such as the well-known Secretary Problem, we propose the Sequential DRA (SDRA) model. Our model introduces a sequential aspect in the decision process over the sample of each round, offering a completely new perspective to the dynamic DP control. Finally, we incorporate several sequential selection algorithms to SDRA control strategies and compare their performance in SIS epidemic simulations.


title: Good Lambda Inequalities and Variable Orlicz Hardy Spaces

abstract: We prove a general theorem showing that local good-$\lambda$ inequalities imply bounds in certain variable Orlicz spaces. We use this to prove results about variable Orlicz Hardy spaces in the unit disc.


title: Bogolyubov invariant via relative spectral invariants on manifolds

abstract: We introduce and study new spectral invariant of two elliptic partial differential operators of Laplace and Dirac type on compact smooth manifolds without boundary that depends on both the eigenvalues and the eigensections of the operators, which is a equal to the regularized number of created particles from the vacuum when the dynamical operator depends on time. We study the asymptotic expansion of this invariant for small adiabatic parameter and compute explicitly the first two coefficients of the asymptotic expansion.


title: The Locally Gaussian Partial Correlation

abstract: It is well known that the dependence structure for jointly Gaussian variables can be fully captured using correlations, and that the conditional dependence structure in the same way can be described using partial correlations. The partial orrelation does not, however, characterize conditional dependence in many non-Gaussian populations. This paper introduces the local Gaussian partial correlation (LGPC), a new measure of conditional dependence. It is a local version of the partial correlation coefficient that characterizes conditional dependence in a large class of populations. It has some useful and novel properties besides: The LGPC reduces to the ordinary partial correlation for jointly normal variables, and it distinguishes between positive and negative conditional dependence. Furthermore, the LGPC can be used to study departures from conditional independence in specific parts of the distribution. We provide several examples of this, both simulated and real, and derive estimation theory under a local likelihood framework. Finally, we indicate how the LGPC can be used to construct a powerful test for conditional independence, which, again, can be used to detect Granger causality in time series.


title: Do internal software quality tools measure validated metrics?

abstract: Internal software quality determines the maintainability of the software product and influences the quality in use. There is a plethora of metrics which purport to measure the internal quality of software, and these metrics are offered by static software analysis tools. To date, a number of reports have assessed the validity of these metrics. No data are available, however, on whether metrics offered by the tools are somehow validated in scientific studies. The current study covers this gap by providing data on which tools and how many validated metrics are provided. The results show that a range of metrics that the tools provided do not seem to be validated in the literature and that only a small percentage of metrics are validated in the provided tools.


title: Wavefront error tolerancing for direct imaging of exo-Earths with a large segmented telescope in space

abstract: Direct imaging of exo-Earths and search for life is one of the most exciting and challenging objectives for future space observatories. Segmented apertures in space will be required to reach the needed large diameters beyond the capabilities of current or planned launch vehicles. These apertures present additional challenges for high-contrast coronagraphy, not only in terms of static phasing but also in terms of their stability. The Pair-based Analytical model for Segmented Telescope Imaging from Space (PASTIS) was developed to model the effects of segment-level optical aberrations on the final image contrast. In this paper, we extend the original PASTIS propagation model from a purely analytical to a semi-analytical method, in which we substitute the use of analytical images with numerically simulated images. The inversion of this model yields a set of orthonormal modes that can be used to determine segment-level wavefront tolerances. We present results in the case of segment-level piston error applied to the baseline coronagraph design of LUVOIR A, with minimum and maximum wavefront error constraint between 56 pm and 290 pm per segment. The analysis is readily generalizable to other segment-level aberrations modes, and can also be expanded to establish stability tolerances for these missions.


title: From the Monster to Thompson to O'Nan

abstract: The commencement of monstrous moonshine is a connection between the largest sporadic simple group---the monster---and complex elliptic curves. Here we explain how a closer look at this connection leads, via the Thompson group, to recently observed relationships between the non-monstrous sporadic simple group of O'Nan and certain families of elliptic curves defined over the rationals. We also describe umbral moonshine from this perspective.


title: Arecibo 4.5/1.4/0.33-GHz Polarimetric Single-Pulse Emission Survey

abstract: We report on an Arecibo 4.5-GHz polarimetric single-pulse survey of the brightest pulsars at high frequency within its sky. The high frequency profiles are accompanied by a collection of both previously published and unpublished high quality 1.4- and 0.33-GHz observations. Here our analyses and discussion primarily involve the average and statistical properties of the 46 pulsar's polarimetric pulse sequences, profile classification and frequency evolution, and polarimetric profiles and peak-occurrence histograms. In most cases both the fractional linear polarization and profile widths decrease with frequency as expected, but there are some exceptions. Similarly, we were able to review and/or extend the profile classifications for this population of pulsars and work out their beaming characteristics quantitatively showing that almost all show properties compatible with the core/double-cone emission beam model. The entirety of these observation's average profiles are accessible for download.


title: The Symplectic Fueter-Sce Theorem

abstract: In this paper we present a symplectic analogue of the Fueter theorem. This allows the construction of special (polynomial) solutions for the symplectic Dirac operator $D_s$, which is defined as the first-order $\mathfrak{sp}(2n)$-invariant differential operator acting on functions on ${\mathbb R}^{2n}$ taking values in the metaplectic spinor representation.


title: The Hellan-Herrmann-Johnson method with curved elements

abstract: We study the finite element approximation of the Kirchhoff plate equation on domains with curved boundaries using the Hellan-Herrmann-Johnson (HHJ) method. We prove optimal convergence on domains with piecewise $C^{k+1}$ boundary for $k \geq 1$ when using a parametric (curved) HHJ space. Computational results are given that demonstrate optimal convergence and how convergence degrades when curved triangles of insufficient polynomial degree are used. Moreover, we show that the lowest order HHJ method on a polygonal approximation of the disk does not succumb to the classic Babu\v{s}ka paradox, highlighting the geometrically non-conforming aspect of the HHJ method.


title: Revisiting the Asymptotic Optimality of RRT$^*$

abstract: RRT* is one of the most widely used sampling-based algorithms for asymptotically-optimal motion planning. This algorithm laid the foundations for optimality in motion planning as a whole, and inspired the development of numerous new algorithms in the field, many of which build upon RRT* itself. In this paper, we first identify a logical gap in the optimality proof of RRT*, which was developed in Karaman and Frazzoli (2011). Then, we present an alternative and mathematically-rigorous proof for asymptotic optimality. Our proof suggests that the connection radius used by RRT* should be increased from $\gamma \left(\frac{\log n}{n}\right)^{1/d}$ to $\gamma' \left(\frac{\log n}{n}\right)^{1/(d+1)}$ in order to account for the additional dimension of time that dictates the samples' ordering. Here $\gamma$, $\gamma'$, are constants, and $n$, $d$, are the number of samples and the dimension of the problem, respectively.


title: Optical Enhancement of Superconductivity via Targeted Destruction of Charge Density Waves

abstract: It has been experimentally established that the occurrence of charge density waves is a common feature of various under-doped cuprate superconducting compounds. The observed states, which are often found in the form of bond density waves (BDW), often occur in a temperature regime immediately above the superconducting transition temperature. Motivated by recent optical experiments on superconducting materials, where it has been shown that optical irradiation can transiently improve the superconducting features, here, we propose a new approach for the enhancement of superconductivity by the targeted destruction of the BDW order. Since BDW states are usually found in competition with superconductivity, suppression of the BDW order enhances the tendency of electrons to form Cooper pairs after reaching a steady-state. By investigating the optical coupling of gapless, collective fluctuations of the BDW modes, we argue that the resonant excitation of these modes can melt the underlying BDW order parameter. We propose an experimental setup to implement such an optical coupling using 2D plasmon-polariton hybrid systems.


title: A Deep Learning-Based Approach for Measuring the Domain Similarity of Persian Texts

abstract: In this paper, we propose a novel approach for measuring the degree of similarity between categories of two pieces of Persian text, which were published as descriptions of two separate advertisements. We built an appropriate dataset for this work using a dataset which consists of advertisements posted on an e-commerce website. We generated a significant number of paired texts from this dataset and assigned each pair a score from 0 to 3, which demonstrates the degree of similarity between the domains of the pair. In this work, we represent words with word embedding vectors derived from word2vec. Then deep neural network models are used to represent texts. Eventually, we employ concatenation of absolute difference and bit-wise multiplication and a fully-connected neural network to produce a probability distribution vector for the score of the pairs. Through a supervised learning approach, we trained our model on a GPU, and our best model achieved an F1 score of 0.9865.


title: NSURL-2019 Shared Task 8: Semantic Question Similarity in Arabic

abstract: Question semantic similarity (Q2Q) is a challenging task that is very useful in many NLP applications, such as detecting duplicate questions and question answering systems. In this paper, we present the results and findings of the shared task (Semantic Question Similarity in Arabic). The task was organized as part of the first workshop on NLP Solutions for Under Resourced Languages (NSURL 2019) The goal of the task is to predict whether two questions are semantically similar or not, even if they are phrased differently. A total of 9 teams participated in the task. The datasets created for this task are made publicly available to support further research on Arabic Q2Q.


title: New Perspectives on the Schr{\"o}dinger-Pauli Theory of Electrons: Part I

abstract: Schr{\"o}dinger-Pauli (SP) theory is a description of electrons in the presence of a static electromagnetic field in which the interaction of the magnetic field with both the orbital and spin moments is explicitly considered. The theory is described from a new perspective, viz. that of the individual electron via its equation of motion or `Quantal Newtonian' first law. The law leads to new physical and mathematical insights into the theory. The law is in terms of `classical' fields whose sources are quantum mechanical expectation values of Hermitian operators taken with respect to the system wave function. The law states that each electron experiences an external and an internal field, the sum of which vanish. The external field is the sum of the electrostatic and a Lorentz field. The internal field is a sum of fields: the electron-interaction, differential density, kinetic, and internal magnetic fields. These fields are respectively representative of a property of the system: electron correlations due to the Pauli exclusion principle and Coulomb repulsion, the electron density, kinetic effects, and the physical current density. The energy components can be expressed in integral virial form in terms of these fields. The law leads to the further understanding that the Hamiltonian is an exactly known and universal functional of the wave function. This allows for the generalization of the SP equation, which then proves it to be intrinsically self-consistent. A Quantal density functional (local effective potential) theory of the SP system is developed. Further generalizations of the present work to the temporal case, and relativistic Dirac theory are proposed.


title: Formal Connections between Template and Anchor Models via Approximate Simulation

abstract: Reduced-order template models like the Linear Inverted Pendulum (LIP) and Spring-Loaded Inverted Pendulum (SLIP) are widely used tools for controlling high-dimensional humanoid robots. However, connections between templates and whole-body models have lacked formal underpinnings, preventing formal guarantees when it comes to integrated controller design. We take a small step towards addressing this gap by considering the notion of approximate simulation. Derived from simulation relations for discrete transition systems in formal methods, approximate similarity means that the outputs of two systems can remain $\epsilon$-close. In this paper, we consider the case of controlling a balancer via planning with the LIP model. We show that the balancer approximately simulates the LIP and derive linear constraints that are sufficient conditions for maintaining ground contact. This allows for rapid planning and replanning with the template model by solving a quadratic program that enforces contact constraints in the full model. We demonstrate the efficacy of this planning and control paradigm in a simulated push recovery scenario for a planar 4-link balancer.


title: An inversion formula with hypergeometric polynomials and application to singular integral operators

abstract: Given parameters $x \notin \mathbb{R}^- \cup \{1\}$ and $\nu$, $\mathrm{Re}(\nu) < 0$, and the space $\mathscr{H}_0$ of entire functions in $\mathbb{C}$ vanishing at $0$, we consider the family of operators $\mathfrak{L} = c_0 \cdot \delta \circ \mathfrak{M}$ with constant $c_0 = \nu(1-\nu)x/(1-x)$, $\delta = z \, \mathrm{d}/\mathrm{d}z$ and integral operator $\mathfrak{M}$ defined by $$ \mathfrak{M}f(z) = \int_0^1 e^{- \frac{z}{x}t^{-\nu}(1-(1-x)t)} \, f \left ( \frac{z}{x} \, t^{-\nu}(1-t) \right ) \, \frac{\mathrm{d}t}{t}, \qquad z \in \mathbb{C}, $$ for all $f \in \mathscr{H}_0$. Inverting $\mathfrak{L}$ or $\mathfrak{M}$ proves equivalent to solve a singular Volterra equation of the first kind. The inversion of operator $\mathfrak{L}$ on $\mathscr{H}_0$ leads us to derive a new class of linear inversion formulas $T = A(x,\nu) \cdot S \Leftrightarrow S = B(x,\nu) \cdot T$ between sequences $S = (S_n)_{n \in \mathbb{N}^*}$ and $T = (T_n)_{n \in \mathbb{N}^*}$, where the infinite lower-triangular matrix $A(x,\nu)$ and its inverse $B(x,\nu)$ involve Hypergeometric polynomials $F(\cdot)$, namely $$ \left\{ \begin{array}{ll} A_{n,k}(x,\nu) = \displaystyle (-1)^k\binom{n}{k}F(k-n,-n\nu;-n;x), B_{n,k}(x,\nu) = \displaystyle (-1)^k\binom{n}{k}F(k-n,k\nu;k;x) \end{array} \right. $$ for $1 \leqslant k \leqslant n$. Functional relations between the ordinary (resp. exponential) generating functions of the related sequences $S$ and $T$ are also given. These relations finally enable us to derive the integral representation $$ \mathfrak{L}^{-1}f(z) = \frac{1-x}{2i\pi x} \, e^{z} \int_{(0+)}^1 \frac{e^{-xtz}}{t(1-t)} \, f \left ( xz \, (-t)^{\nu}(1-t)^{1-\nu} \right ) \, \mathrm{d}t, \quad z \in \mathbb{C}, $$ for the inverse $\mathfrak{L}^{-1}$ of operator $\mathfrak{L}$ on $\mathscr{H}_0$, where the integration contour encircles the point 0.


title: Epidemic spreading on modular networks: The fear to declare a pandemic

abstract: In the past few decades, the frequency of pandemics has been increased due to the growth of urbanization and mobility among countries. Since a disease spreading in one country could become a pandemic with a potential worldwide humanitarian and economic impact, it is important to develop models to estimate the probability of a worldwide pandemic. In this paper, we propose a model of disease spreading in a structural modular complex network (having communities) and study how the number of bridge nodes $n$ that connect communities affects disease spread. We find that our model can be described at a global scale as an infectious transmission process between communities with global infectious and recovery time distributions that depend on the internal structure of each community and $n$. We find that near the critical point as $n$ increases, the disease reaches most of the communities, but each community has only a small fraction of recovered nodes. In addition, we obtain that in the limit $n \to \infty$, the probability of a pandemic increases abruptly at the critical point. This scenario could make the decision on whether to launch a pandemic alert or not more difficult. Finally, we show that link percolation theory can be used at a global scale to estimate the probability of a pandemic since the global transmissibility between communities has a weak dependence on the global recovery time.


title: A Gated Self-attention Memory Network for Answer Selection

abstract: Answer selection is an important research problem, with applications in many areas. Previous deep learning based approaches for the task mainly adopt the Compare-Aggregate architecture that performs word-level comparison followed by aggregation. In this work, we take a departure from the popular Compare-Aggregate architecture, and instead, propose a new gated self-attention memory network for the task. Combined with a simple transfer learning technique from a large-scale online corpus, our model outperforms previous methods by a large margin, achieving new state-of-the-art results on two standard answer selection datasets: TrecQA and WikiQA.


title: Continuous variable direct secure quantum communication using Gaussian states

abstract: Continuous variable one-way and controlled-two-way secure direct quantum communication schemes have been designed using Gaussian states. Specifically, a scheme for continuous variable quantum secure direct communication and another scheme for continuous variable controlled quantum dialogue are proposed using single-mode squeezed coherent states. The security of the proposed schemes against a set of attacks (e.g., Gaussian quantum cloning machine and intercept resend attacks) has been proved. Further, it is established that the proposed schemes do not require two-mode squeezed states which are essential for a set of existing proposals. The controlled two-way communication scheme is shown to be very general in nature as it can be reduced to schemes for various relatively simpler cryptographic tasks like controlled deterministic secure communication, quantum dialogue, quantum key distribution. In addition, it is briefly discussed that the proposed schemes can provide us tools to design quantum cryptographic solutions for several socioeconomic problems.


title: A LN$_2$ Based Cooling System for a Next Generation Liquid Xenon Dark Matter Detector

abstract: In recent years cooling technology for Liquid Xenon (LXe) detectors has advanced driven by the development of Dark Matter (DM) detectors with target mass in the 100 - 1,000 kg range. The next generation of DM detectors based on LXe will be in the 50,000 kg (50 t) range requiring more than 1 kW of cooling power. Most of the prior cooling methods become impractical at this level. For cooling a 50 t scale LXe detector, a method is proposed in which Liquid Nitrogen (LN$_2$) in a small local reservoir cools the xenon gas via a cold finger. The cold finger incorporates a heating unit to provide temperature regulation. The proposed cooling method is simple, reliable, and suitable for the required long-term operation for a rare event search. The device can be easily integrated into present cooling systems, e.g. the 'Cooling Bus' employed for the PandaX I and II experiments. It is still possible to cool indirectly with no part of the cooling or temperature control system getting in direct contact with the clean xenon in the detector. Also the cooling device can be mounted at a large distance, i.e. the detector is cooled remotely from a distance of 5 - 10 m. The method was tested in a laboratory setup at Columbia University to carry out different measurements with a small LXe detector and behaved exactly as predicted.


title: Induction and Reference of Entities in a Visual Story

abstract: We are enveloped by stories of visual interpretations in our everyday lives. The way we narrate a story often comprises of two stages, which are, forming a central mind map of entities and then weaving a story around them. A contributing factor to coherence is not just basing the story on these entities but also, referring to them using appropriate terms to avoid repetition. In this paper, we address these two stages of introducing the right entities at seemingly reasonable junctures and also referring them coherently in the context of visual storytelling. The building blocks of the central mind map, also known as entity skeleton are entity chains including nominal and coreference expressions. This entity skeleton is also represented in different levels of abstractions to compose a generalized frame to weave the story. We build upon an encoder-decoder framework to penalize the model when the decoded story does not adhere to this entity skeleton. We establish a strong baseline for skeleton informed generation and then extend this to have the capability of multitasking by predicting the skeleton in addition to generating the story. Finally, we build upon this model and propose a glocal hierarchical attention model that attends to the skeleton both at the sentence (local) and the story (global) levels. We observe that our proposed models outperform the baseline in terms of automatic evaluation metric, METEOR. We perform various analysis targeted to evaluate the performance of our task of enforcing the entity skeleton such as the number and diversity of the entities generated. We also conduct human evaluation from which it is concluded that the visual stories generated by our model are preferred 82% of the times. In addition, we show that our glocal hierarchical attention model improves coherence by introducing more pronouns as required by the presence of nouns.


title: Retrofitting Contextualized Word Embeddings with Paraphrases

abstract: Contextualized word embedding models, such as ELMo, generate meaningful representations of words and their context. These models have been shown to have a great impact on downstream applications. However, in many cases, the contextualized embedding of a word changes drastically when the context is paraphrased. As a result, the downstream model is not robust to paraphrasing and other linguistic variations. To enhance the stability of contextualized word embedding models, we propose an approach to retrofitting contextualized embedding models with paraphrase contexts. Our method learns an orthogonal transformation on the input space, which seeks to minimize the variance of word representations on paraphrased contexts. Experiments show that the retrofitted model significantly outperforms the original ELMo on various sentence classification and language inference tasks.


title: New Perspectives on the Schr{\"o}dinger-Pauli Theory of Electrons: Part II: Application to the Triplet State of a Quantum Dot in a Magnetic Field

abstract: The Schr\"odinger-Pauli theory of electrons in the presence of a static electromagnetic field can be described from the perspective of the individual electron via its equation of motion or `Quantal Newtonian' first law. The law is in terms of `classical' fields whose sources are quantum-mechanical expectation values of Hermitian operators taken with respect to the wave function. The law states that the sum of the external and internal fields experienced by each electron vanishes. The external field is the sum of the binding electrostatic and Lorentz fields. The internal field is the sum of fields representative of properties of the system: electron correlations due to the Pauli exclusion principle and Coulomb repulsion; the electron density; kinetic effects; the current density. Thus, the internal field is a sum of the electron-interaction, differential density, kinetic, and internal magnetic fields. The energy can be expressed in integral virial form in terms of these fields. Via this perspective, the Schr\"odinger-Pauli equation can be written in a generalized form which then shows it to be intrinsically self-consistent. This new perspective is explicated by application to the triplet $2^{3}S$ state of a 2-D 2-electron quantum dot in a magnetic field. The quantal sources of the density; the paramagnetic, diamagnetic, and magnetization current densities; pair-correlation density; the Fermi-Coulomb hole charge; and the single-particle density matrix are obtained, and from them the corresponding fields determined. The fields are shown to satisfy the `Quantal Newtonian' first law. The components of the energy too are determined from these fields. Finally, the example is employed to demonstrate the intrinsic self-consistent nature of the Schr\"odinger-Pauli equation.


title: Using Clinical Notes with Time Series Data for ICU Management

abstract: Monitoring patients in ICU is a challenging and high-cost task. Hence, predicting the condition of patients during their ICU stay can help provide better acute care and plan the hospital's resources. There has been continuous progress in machine learning research for ICU management, and most of this work has focused on using time series signals recorded by ICU instruments. In our work, we show that adding clinical notes as another modality improves the performance of the model for three benchmark tasks: in-hospital mortality prediction, modeling decompensation, and length of stay forecasting that play an important role in ICU management. While the time-series data is measured at regular intervals, doctor notes are charted at irregular times, making it challenging to model them together. We propose a method to model them jointly, achieving considerable improvement across benchmark tasks over baseline time-series model. Our implementation can be found at \url{https://github.com/kaggarwal/ClinicalNotesICU}.


title: Risk-sensitive safety specifications for stochastic systems using Conditional Value-at-Risk

abstract: This paper proposes a safety analysis method that facilitates a tunable balance between the worst-case and risk-neutral perspectives. First, we define a risk-sensitive safe set to specify the degree of safety attained by a stochastic system. This set is defined as a sublevel set of the solution to an optimal control problem that is expressed using the Conditional Value-at-Risk (CVaR) measure. This problem does not satisfy Bellman's Principle, thus our next contribution is to show how risk-sensitive safe sets can be under-approximated by the solution to a CVaR-Markov Decision Process. We adopt an existing value iteration algorithm to find an approximate solution to the reduced problem for a class of linear systems. Then, we develop a realistic numerical example of a stormwater system to show that this approach can be applied to non-linear systems. Finally, we compare the CVaR criterion to the exponential disutility criterion. The latter allocates control effort evenly across the cost distribution to reduce variance, while the CVaR criterion focuses control effort on a given worst-case quantile--where it matters most for safety.


title: Measuring Domain Portability and ErrorPropagation in Biomedical QA

abstract: In this work we present Google's submission to the BioASQ 7 biomedical question answering (QA) task (specifically Task 7b, Phase B). The core of our systems are based on BERT QA models, specifically the model of \cite{alberti2019bert}. In this report, and via our submissions, we aimed to investigate two research questions. We start by studying how domain portable are QA systems that have been pre-trained and fine-tuned on general texts, e.g., Wikipedia. We measure this via two submissions. The first is a non-adapted model that uses a public pre-trained BERT model and is fine-tuned on the Natural Questions data set \cite{kwiatkowski2019natural}. The second system takes this non-adapted model and fine-tunes it with the BioASQ training data. Next, we study the impact of error propagation in end-to-end retrieval and QA systems. Again we test this via two submissions. The first uses human annotated relevant documents and snippets as input to the model and the second predicted documents and snippets. Our main findings are that domain specific fine-tuning can benefit Biomedical QA. However, the biggest quality bottleneck is at the retrieval stage, where we see large drops in metrics -- over 10pts absolute -- when using non gold inputs to the QA model.


title: A Layered Architecture for Active Perception: Image Classification using Deep Reinforcement Learning

abstract: We propose a planning and perception mechanism for a robot (agent), that can only observe the underlying environment partially, in order to solve an image classification problem. A three-layer architecture is suggested that consists of a meta-layer that decides the intermediate goals, an action-layer that selects local actions as the agent navigates towards a goal, and a classification-layer that evaluates the reward and makes a prediction. We design and implement these layers using deep reinforcement learning. A generalized policy gradient algorithm is utilized to learn the parameters of these layers to maximize the expected reward. Our proposed methodology is tested on the MNIST dataset of handwritten digits, which provides us with a level of explainability while interpreting the agent's intermediate goals and course of action.


title: Do Compressed Representations Generalize Better?

abstract: One of the most studied problems in machine learning is finding reasonable constraints that guarantee the generalization of a learning algorithm. These constraints are usually expressed as some simplicity assumptions on the target. For instance, in the Vapnik-Chervonenkis (VC) theory the space of possible hypotheses is considered to have a limited VC dimension. In this paper, the constraint on the entropy $H(X)$ of the input variable $X$ is studied as a simplicity assumption. It is proven that the sample complexity to achieve an $\epsilon$-$\delta$ Probably Approximately Correct (PAC) hypothesis is bounded by $\frac{2^{ \left.6H(X)\middle/\epsilon\right.}+\log{\frac{1}{\delta}}}{\epsilon^2}$ which is sharp up to the $\frac{1}{\epsilon^2}$ factor. Morever, it is shown that if a feature learning process is employed to learn the compressed representation from the dataset, this bound no longer exists. These findings have important implications on the Information Bottleneck (IB) theory which had been utilized to explain the generalization power of Deep Neural Networks (DNNs), but its applicability for this purpose is currently under debate by researchers. In particular, this is a rigorous proof for the previous heuristic that compressed representations are exponentially easier to be learned. However, our analysis pinpoints two factors preventing the IB, in its current form, to be applicable in studying neural networks. Firstly, the exponential dependence of sample complexity on $\frac{1}{\epsilon}$, which can lead to a dramatic effect on the bounds in practical applications when $\epsilon$ is small. Secondly, our analysis reveals that arguments based on input compression are inherently insufficient to explain generalization of methods like DNNs in which the features are also learned using available data.


title: Integrable systems in planar robotics

abstract: The main purpose of this paper is to investigate commuting flows and integrable systems on the configuration spaces of planar linkages. Our study leads to the definition of a natural volume form on each configuration space of planar linkages, the notion of cross products of integrable systems, and also the notion of multi-Nambu integrable systems. The first integrals of our systems are functions of Bott-Morse type, which may be used to study the topology of configuration spaces. Dedicated to Anatoly T. Fomenko on the occasion of his 75th birthday


title: Measuring Conceptual Entanglement in Collections of Documents

abstract: Conceptual entanglement is a crucial phenomenon in quantum cognition because it implies that classical probabilities cannot model non--compositional conceptual phenomena. While several psychological experiments have been developed to test conceptual entanglement, this has not been explored in the context of Natural Language Processing. In this paper, we apply the hypothesis that words of a document are traces of the concepts that a person has in mind when writing the document. Therefore, if these concepts are entangled, we should be able to observe traces of their entanglement in the documents. In particular, we test conceptual entanglement by contrasting language simulations with results obtained from a text corpus. Our analysis indicates that conceptual entanglement is strongly linked to the way in which language is structured. We discuss the implications of this finding in the context of conceptual modeling and of Natural Language Processing.


title: SkyNet: a Hardware-Efficient Method for Object Detection and Tracking on Embedded Systems

abstract: Object detection and tracking are challenging tasks for resource-constrained embedded systems. While these tasks are among the most compute-intensive tasks from the artificial intelligence domain, they are only allowed to use limited computation and memory resources on embedded devices. In the meanwhile, such resource-constrained implementations are often required to satisfy additional demanding requirements such as real-time response, high-throughput performance, and reliable inference accuracy. To overcome these challenges, we propose SkyNet, a hardware-efficient neural network to deliver the state-of-the-art detection accuracy and speed for embedded systems. Instead of following the common top-down flow for compact DNN (Deep Neural Network) design, SkyNet provides a bottom-up DNN design approach with comprehensive understanding of the hardware constraints at the very beginning to deliver hardware-efficient DNNs. The effectiveness of SkyNet is demonstrated by winning the competitive System Design Contest for low power object detection in the 56th IEEE/ACM Design Automation Conference (DAC-SDC), where our SkyNet significantly outperforms all other 100+ competitors: it delivers 0.731 Intersection over Union (IoU) and 67.33 frames per second (FPS) on a TX2 embedded GPU; and 0.716 IoU and 25.05 FPS on an Ultra96 embedded FPGA. The evaluation of SkyNet is also extended to GOT-10K, a recent large-scale high-diversity benchmark for generic object tracking in the wild. For state-of-the-art object trackers SiamRPN++ and SiamMask, where ResNet-50 is employed as the backbone, implementations using our SkyNet as the backbone DNN are 1.60X and 1.73X faster with better or similar accuracy when running on a 1080Ti GPU, and 37.20X smaller in terms of parameter size for significantly better memory and storage footprint.


title: An Efficient Move Blocking Strategy for Multiple Shooting based Nonlinear Model Predictive Control

abstract: Move blocking (MB) is a widely used strategy to reduce the degrees of freedom of the Optimal Control Problem (OCP) arising in receding horizon control. The size of the OCP is reduced by forcing the input variables to be constant over multiple discretization steps. In this paper, we focus on developing computationally efficient MB schemes for multiple shooting based nonlinear model predictive control (NMPC). The degrees of freedom of the OCP is reduced by introducing MB in the shooting step, resulting in a smaller but sparse OCP. Therefore, the discretization accuracy and level of sparsity is maintained. A condensing algorithm that exploits the sparsity structure of the OCP is proposed, that allows to reduce the computation complexity of condensing from quadratic to linear in the number of discretization nodes. As a result, active-set methods with warm-start strategy can be efficiently employed, thus allowing the use of a longer prediction horizon. A detailed comparison between the proposed scheme and the nonuniform grid NMPC is given. Effectiveness of the algorithm in reducing computational burden while maintaining optimization accuracy and constraints fulfillment is shown by means of simulations with two different problems.


title: The complications of learning from Super Soft Source X-ray spectra

abstract: Super Soft X-ray Sources (SSS) are powered by nuclear burning on the surface of an accreting white dwarf, they are seen around 0.1-1 keV (thus in the soft X-ray regime), depending on effective temperature and the amount of intervening interstellar neutral hydrogen N_H. The most realistic model to derive physical parameters from observed SSS spectra would be an atmosphere model that simulates the radiation transport processes. However, observed SSS high-resolution grating spectra reveal highly complex details that cast doubts on the feasibility of achieving unique results from atmosphere modeling. In this article, I discuss two independent atmosphere model analyses of the same data set, leading to different results. I then show some of the details that complicate the analysis and conclude that we need to approach the interpretation of high-resolution SSS spectra differently. We need to focus more on the data than the models and to use more phenomenological approaches as is traditionally done with optical spectra.


title: Learning an Adaptive Learning Rate Schedule

abstract: The learning rate is one of the most important hyper-parameters for model training and generalization. However, current hand-designed parametric learning rate schedules offer limited flexibility and the predefined schedule may not match the training dynamics of high dimensional and non-convex optimization problems. In this paper, we propose a reinforcement learning based framework that can automatically learn an adaptive learning rate schedule by leveraging the information from past training histories. The learning rate dynamically changes based on the current training dynamics. To validate this framework, we conduct experiments with different neural network architectures on the Fashion MINIST and CIFAR10 datasets. Experimental results show that the auto-learned learning rate controller can achieve better test results. In addition, the trained controller network is generalizable -- able to be trained on one data set and transferred to new problems.


title: Beyond the one-dimensional configuration coordinate model of photoluminescence

abstract: The one-dimensional configuration coordinate model (1D-CCM) is widely used for the analysis of photoluminescence in molecules and doped solids, and relies on a linear combination of the equilibrium nuclear configurations of ground and excited states. It delivers an estimation of the energy barrier at which ground and excited state curves cross, semi-classically linked to non-radiative transition rate and thermal quenching. To assess its predictive power for the latter properties, we propose a new \textit{optimized configuration path (OCP) method in which} the ground-state and excited-state forces are mixed instead of their configurations. We also define another one-parameter model thanks a double energy parabola hypothesis (DEPH). We compare the OCP method and the DEPH reference with the 1D-CCM for three paradigmatic 4f-5d phosphors Y$_3$Al$_5$O$_{12}$:Ce, Lu$_2$SiO$_5$:Ce, and YAlO$_3$:Ce. We find that the OCP and DEPH methods yield similar results with geometries that have significantly lower ground-state energies than the 1D-CCM for the same 4f-5d energy difference. However, the OCP method suffers from the appearance of multiple local minima, rendering the clear determination of the optimal geometry very difficult in practice. Still the OCP method allows one to quantify the deviations from the 1D-CCM, therefore increasing confidence in the lower bound obtained from the DEPH for the 4f-5d crossing barrier, and its comparison with the energy of the auto-ionization thermal quenching mechanism. We expect the OCP approach to be applicable to other luminescent materials or molecules.


title: Recursive Calculation of the Optical Response of Multicomponent Metamaterials

abstract: We develop a recursive computational procedure to efficiently calculate the macroscopic dielectric function of multi-component metamaterials of arbitrary geometry and composition within the long wavelength approximation. Although the microscopic response of the system might correspond to non-Hermitian operators, we develop a representation of the microscopic fields and of the response, and we introduce an appropriate metric that makes all operators symmetric. This allows us to use a modified Haydock recursion, introducing complex Haydock coefficients that allow an efficient computation of the macroscopic response and the microscopic fields. We test our procedure comparing our results to analytical ones in simple systems, and verifying they obey a generalized multicomponent Keller's theorem and the Mortola and Stef\'e's theorem for four component metalic and dielectric systems.