title: Improved Bounds for the Greedy Strategy in Optimization Problems with Curvature

abstract: Consider the problem of choosing a set of actions to optimize an objective function that is a real-valued polymatroid function subject to matroid constraints. The greedy strategy provides an approximate solution to the optimization problem, and it is known to satisfy some performance bounds in terms of the total curvature. The total curvature depends on the value of objective function on sets outside the constraint matroid. If we are given a function defined only on the matroid, the problem still makes sense, but the existing bounds involving the total curvature do not apply. This is puzzling: If the optimization problem is perfectly well defined, why should the bounds no longer apply? This motivates an alternative formulation of such bounding techniques. The first question that comes to mind is whether it is possible to extend a polymatroid function defined on a matroid to one on the entire power set. This was recently shown to be negative in general. Here, we provide necessary and sufficient conditions for the existence of an \emph{incremental} extension of a polymatroid function defined on the uniform matroid of rank $k$ to one defined on the uniform matroid of rank $k+1$, together with an algorithm for constructing the extension. Whenever a polymatroid objective function defined on a matroid can be extended to the entire power set, the greedy approximation bounds involving the total curvature of the extension apply. However, these bounds still depend on sets outside the constraint matroid. Motivated by this, we define a new notion of curvature called \emph{partial curvature}, involving only sets in the matroid. We derive necessary and sufficient conditions for an extension of the function to have a total curvature that is equal to the partial curvature. Moreover, we prove that the bounds in terms of the partial curvature are in general improved over the previous ones......


title: Plasma vs Drude modelling of the Casimir force: beyond the proximity force approximation

abstract: We calculate the Casimir force and its gradient between a spherical and a planar gold surface. Significant numerical improvements allow us to extend the range of accessible parameters into the experimental regime. We compare our numerically exact results with those obtained within the proximity force approximation (PFA) employed in the analysis of all Casimir force experiments reported in the literature so far. Special attention is paid to the difference between the Drude model and the dissipationless plasma model at zero frequency. It is found that the correction to PFA is too small to explain the discrepancy between the experimental data and the PFA result based on the Drude model. However, it turns out that for the plasma model, the corrections to PFA lie well outside the experimental bound obtained by probing the variation of the force gradient with the sphere radius [D. E. Krause et al., Phys. Rev. Lett. 98, 050403 (2007)]. The corresponding corrections based on the Drude model are significantly smaller but still in violation of the experimental bound for small distances between plane and sphere.


title: Iso-vector dipole resonance and shear viscosity in low energy heavy-ion collision

abstract: The ratio of shear viscosity over entropy density in low energy heavy-ion collision has been calculated by using the Green-Kubo method in the framework of an extended quantum molecular dynamics model. After the system almost reaching a local equilibration for a head-on $^{40}$Ca+$^{100}$Mo collision, thermodynamic and transport properties are extracted. Meanwhile, iso-vector giant dipole resonance (IVGDR) of the collision system is also studied. By the Gaussian fits to the IVGDR photon spectra, the peak energies of IVGDR are extracted at different incident energies. The result shows that the IVGDR peak energy has a positive correlation with the ratio of shear viscosity over entropy density. This is a quantum effect and indicates a difference between nuclear matter and classical fluid.


title: A matrix characterization of boundary representations of positive matrices in the Hardy space

abstract: Spectral measures give rise to a natural harmonic analysis on the unit disc via a boundary representation of a positive matrix arising from a spectrum of the measure. We consider in this paper the reverse: for a positive matrix in the Hardy space of the unit disc we consider which measures, if any, yield a boundary representation of the positive matrix. We introduce a potential characterization of those measures via a matrix identity and show that the characterization holds in several important special cases.


title: Prediction of Stable Ground-State Lithium Polyhydrides under High Pressures

abstract: Hydrogen-rich compounds are important for understanding the dissociation of dense molecular hydrogen, as well as searching for room temperature Bardeen-Cooper-Schrieffer (BCS) superconductors. A recent high pressure experiment reported the successful synthesis of novel insulating lithium polyhydrides when above 130 GPa. However, the results are in sharp contrast to previous theoretical prediction by PBE functional that around this pressure range all lithium polyhydrides (LiHn (n = 2-8)) should be metallic. In order to address this discrepancy, we perform unbiased structure search with first principles calculation by including the van der Waals interaction that was ignored in previous prediction to predict the high pressure stable structures of LiHn (n = 2-11, 13) up to 200 GPa. We reproduce the previously predicted structures, and further find novel compositions that adopt more stable structures. The van der Waals functional (vdW-DF) significantly alters the relative stability of lithium polyhydrides, and predicts that the stable stoichiometries for the ground-state should be LiH2 and LiH9 at 130-170 GPa, and LiH2, LiH8 and LiH10 at 180-200 GPa. Accurate electronic structure calculation with GW approximation indicates that LiH, LiH2, LiH7, and LiH9 are insulative up to at least 208 GPa, and all other lithium polyhydrides are metallic. The calculated vibron frequencies of these insulating phases are also in accordance with the experimental infrared (IR) data. This reconciliation with the experimental observation suggests that LiH2, LiH7, and LiH9 are the possible candidates for lithium polyhydrides synthesized in that experiment. Our results reinstate the credibility of density functional theory in description H-rich compounds, and demonstrate the importance of considering van der Waals interaction in this class of materials.

categories: cond-mat.mtrl-sci cond-mat.other physics.app-ph physics.chem-ph physics.comp-ph


title: On Benjamini--Schramm limits of congruence subgroups

abstract: Every sequence of orbifolds corresponding to pairwise non-conjugate congruence lattices in a higher rank semisimple group over local fields of zero characteristic is Benjamini--Schramm convergent to the universal cover.


title: What Can We Learn from Noise? -- Mesoscopic Nonequilibrium Statistical Physics --

abstract: Mesoscopic systems -- small electric circuits working in quantum regime -- offer us a unique experimental stage to explorer quantum transport in a tunable and precise way. The purpose of this Review is to show how they can contribute to statistical physics. We introduce the significance of fluctuation, or equivalently noise, as noise measurement enables us to address the fundamental aspects of a physical system. The significance of the fluctuation theorem (FT) in statistical physics is noted. We explain what information can be deduced from the current noise measurement in mesoscopic systems. As an important application of the noise measurement to statistical physics, we describe our experimental work on the current and current noise in an electron interferometer, which is the first experimental test of FT in quantum regime. Our attempt will shed new light in the research field of mesoscopic quantum statistical physics.


title: Temporal dynamics of online petitions

abstract: Online petitions are an important avenue for direct political action, yet the dynamics that determine when a petition will be successful are not well understood. Here we analyze the temporal characteristics of online-petition signing behavior in order to identify systematic differences between popular petitions, which receive a high volume of signatures, and unpopular ones. We find that, in line with other temporal characterizations of human activity, the signing process is typically non-Poissonian and non-homogeneous in time. However, this process exhibits anomalously high memory for human activity, possibly indicating that synchronized external influence or contagion play and important role. More interestingly, we find clear differences in the characteristics of the inter-event time distributions depending on the total number of signatures that petitions receive, independently of the total duration of the petitions. Specifically, popular petitions that attract a large volume of signatures exhibit more variance in the distribution of inter-event times than unpopular petitions with only a few signatures, which could be considered an indication that the former are more bursty. However, petitions with large signature volume are less bursty according to measures that consider the time ordering of inter-event times. Our results, therefore, emphasize the importance of accounting for time ordering to characterize human activity.


title: Investigation of giant dipole resonance in heavy deformed nuclei with the EQMD model

abstract: The deformation evolution of giant dipole resonance (GDR), in the chains of Sm and Nd isotopes, are investigated in the framework of an extended quantum molecular dynamics (EQMD) model. The mass number dependence of resonance peak position $(E_{m})$ in the major and minor axis directions of deformed nuclei as well as the difference $\Delta E_{m}$ between them are described in detail. The correlation between the splitting ($\Delta E_{m} /\bar{E}_m $) of the GDR spectra and the deformation($ \beta_{2}$) is further studied. The results confirm that $\Delta E_{m} /\bar{E}_m $ is proportional to $ \beta_{2}$. By comparing the calculation with the experimental data on photon absorption cross section $\sigma_{\gamma}$, it shows that the EQMD model can quite well reproduce the shape of GDR spectra from spherical to prolate shape. The GDR shapes in $^{134}$Sm, $^{136}$Sm, $^{138}$Sm, $^{130}$Nd, $^{132}$Nd and $^{134}$Nd are also predicted. In addition, the symmetry energy coefficient $(E_{sym})$ dependence of GDR spectra of $^{150}$Nd is also discussed. It is found that the calculated GDR spectrum in the EQMD model is perfectly consistent with the experimental results when $E_{sym}$ equals to 32 MeV.


title: Constraining magnetic field amplification in SN shocks using radio observations of SNe 2011fe and 2014J

abstract: We modeled the radio non-detection of two Type Ia supernovae (SNe) 2011fe and 2014J considering synchrotron emission from the interaction between SN ejecta and the circumstellar medium. For an ejecta with the outer part having a power law density structure we compare synchrotron emission with radio observations. Assuming that 20$\%$ of the bulk shock energy is being shared equally between electrons and magnetic fields we found a very low density medium around both the SNe. A less tenuous medium with particle density $\sim$ 1 $\rm cm^{-3}$, which could be expected around both SNe, can be estimated when the magnetic field amplification is less than that presumed for energy equipartition. This conclusion also holds if the progenitor of SN 2014J was a rigidly rotating white dwarf (WD) with a main sequence (MS) or red giant companion. For a He star companion, or a MS for SN 2014J, with 10$\%$ and 1$\%$ of bulk kinetic energy in magnetic fields, we obtain a mass loss rate $< 10^{-9}$ and $< \sim 4\times 10^{-9}$ M$_{\odot}$yr$^{-1}$ for a wind velocity of 100 km/s. The former requires a mass accretion efficiency $>$ 99$\%$ onto the WD, but is less restricted for the latter case. However, if the tenuous medium is due to a recurrent nova it is difficult from our model to predict synchrotron luminosities. Although the formation channels of SNe 2011fe and 2014J are not clear, the null detection in radio wavelengths could point toward a low amplification efficiency for magnetic fields in SN shocks.


title: Flipping out with many flips: hardness of testing k-monotonicity

abstract: A function f : {0, 1}^n -> {0, 1} is said to be k-monotone if it flips between 0 and 1 at most k times on every ascending chain. Such functions represent a natural generalization of (1-)monotone functions, and have been recently studied in circuit complexity, PAC learning, and cryptography. Our work is part of a renewed focus in understanding testability of properties characterized by freeness of arbitrary order patterns as a generalization of monotonicity. Recently, Canonne et al. (ITCS 2017) initiate the study of k-monotone functions in the area of property testing, and Newman et al. (SODA 2017) study testability of families characterized by freeness from order patterns on real-valued functions over the line [n] domain. We study k-monotone functions in the more relaxed parametrized property testing model, introduced by Parnas et al. (JCSS, 72(6), 2006). In this process we resolve a problem left open in previous work. Specifically, our results include the following. 1. Testing 2-monotonicity on the hypercube non-adaptively with one-sided error requires an exponential in \sqrt n number of queries. This behavior shows a stark contrast with testing (1-)monotonicity, which only needs O(\sqrt n) queries (Khot et al. (FOCS 2015)). Furthermore, even the apparently easier task of distinguishing 2-monotone functions from functions that are far from being n^.01 -monotone also requires an exponential number of queries. 2. On the hypercube [n]^d domain, there exists a testing algorithm that makes a constant number of queries and distinguishes functions that are k-monotone from functions that are far from being O(kd^2)-monotone. Such a dependency is likely necessary, given the lower bound above for the hypercube.


title: Nonlinear stability of Gardner breathers

abstract: We show that breather solutions of the Gardner equation, a natural generalization of the KdV and mKdV equations, are $H^2(\mathbb{R})$ stable. Through a variational approach, we characterize Gardner breathers as minimizers of a new Lyapunov functional and we study the associated spectral problem, through $(i)$ the analysis of the spectrum of explicit linear systems (\emph{spectral stability}), and $(ii)$ controlling degenerated directions by using low regularity conservation laws.


title: Bisector surfaces and circumscribed spheres of tetrahedra derived by translation curves in $\SOL$ geometry

abstract: In the present paper we study the $\SOL$ geometry that is one of the eight homogeneous Thurston 3-geomet\-ri\-es. We determine the equation of the translation-like bisector surface of any two points. We prove, that the isosceles property of a translation triangle is not equivalent to two angles of the triangle being equal and that the triangle inequalities do not remain valid for translation triangles in general. Moreover, we develop a method to determine the centre and the radius of the circumscribed translation sphere of a given {\it translation tetrahedron}. In our work we will use for computations and visualizations the projective model of $\SOL$ described by E. Moln\'ar in \cite{M97}.


title: Geometric Graph Manifolds with non-negative scalar curvature

abstract: We classify $n$-dimensional geometric graph manifolds with nonnegative scalar curvature, and first show that if $n>3$, the universal cover splits off a codimension 3 Euclidean factor. We then proceed with the classification of the 3-dimensional case by showing that such a manifold is either a lens space or a prism manifold with a very rigid metric. This allows us to also classify the moduli space of such metrics: it has infinitely many connected components for lens spaces, while it is connected for prism manifolds.


title: An Examination of Rotational Effects on Large Wind Turbine Blades

abstract: CFD studies have been performed to demonstrate the effects of rotation on large wind turbine blades. The 10 MW blade of the Advanced Aerodynamic Tools for Large Rotors (AVATAR) project was selected in the present investigation. In order to demonstrate the rotational augmentation, lift coefficients of the blade sections were compared with the two-dimensional simulations at consistent inflow conditions. The increase of the lift coefficient was observed at the inboard section of the blade, showing the effect of rotation. A significant delay of separation point as well as a massive reduction of the wake size occur at the inboard blade section. The results of the present manuscript provide more information of the rotational augmentation occurring in the near hub region of wind turbine blades.


title: Tracing HI Beyond the Local Universe

abstract: The SKA and its pathfinders will enable studies of HI emission at higher redshifts than ever before. In moving beyond the local Universe, this will require the use of cosmologically appropriate formulae that have traditionally been simplified to their low-redshift approximations. In this paper, we summarise some of the most important relations for tracing HI emission in the SKA era, and present an online calculator to assist in the planning and analysis of observations (hifi.icrar.org).


title: Numerical Investigation of Unsteady Aerodynamic Effects on Thick Flatback Airfoils

abstract: The unsteady characteristics of the flow over thick flatback airfoils have been investigated by means of CFD calculations. Sandia airfoils which have 35% maximum thickness with three different trailing edge thicknesses were selected. The calculations provided good results compared with available experimental data with regard to the lift curve and the impact of trailing edge thickness. Unsteady CFD simulations revealed that the Strouhal number is found to be independent of the lift coefficient before stall and increases with the trailing edge. The present work shows the dependency of the Strouhal number and the wake development on the trailing edge thickness. A recommendation of the Strouhal number definition is given for flatback airfoils by considering the trailing edge separation at low angle of attack. The detailed unsteady characteristics of thick flatback airfoils are discussed more in the present paper.


title: Competitive Equilibrium For Almost All Incomes: Existence and Fairness

abstract: Competitive equilibrium (CE) is a fundamental concept in market economics. Its efficiency and fairness properties make it particularly appealing as a rule for fair allocation of resources among agents with possibly different entitlements. However, when the resources are indivisible, a CE might not exist even when there is one resource and two agents with equal incomes. Recently, Babaioff and Nisan and Talgam-Cohen (2017) have suggested to consider the entire space of possible incomes, and check whether there exists a competitive equilibrium for almost all income-vectors --- all income-space except a subset of measure zero. They proved various existence and non-existence results, but left open the cases of four goods and three or four agents with monotonically-increasing preferences. This paper proves non-existence in both these cases, thus completing the characterization of CE existence for almost all incomes in the domain of monotonically increasing preferences. Additionally, the paper provides a complete characterization of CE existence in the domain of monotonically decreasing preferences, corresponding to allocation of chores. On the positive side, the paper proves that CE exists for almost all incomes when there are four goods and three agents with additive preferences. The proof uses a new tool for describing a CE, as a subgame-perfect equilibrium of a specific sequential game. The same tool also enables substantially simpler proofs to the cases already proved by Babaioff et al. Additionally, this paper proves several strong fairness properties that are satisfied by any CE allocation, illustrating its usefulness for fair allocation among agents with different entitlements.


title: Predicted novel insulating electride compound between alkali metals lithium and sodium under high pressure

abstract: The application of high pressure can fundamentally modify the crystalline and electronic structures of elements as well as their chemical reactivity, which could lead to the formation of novel materials. Here, we explore the reactivity of lithium with sodium under high pressure, using a swarm structure searching techniques combined with first-principles calculations, which identify a thermodynamically stable LiNa compound adopting an orthorhombic oP8 phase at pressure above 355 GPa. The formation of LiNa may be a consequence of strong concentration of electrons transfer from the lithium and the sodium atoms into the interstitial sites, which also leads to opening a relatively wide band gap for LiNa-op8. This is substantially different from the picture that share or exchange electrons in common compounds and alloys. In addition, lattice-dynamic calculations indicate that LiNa-op8 remains dynamically stable when pressure decompresses down to 70 GPa.

categories: physics.chem-ph cond-mat.mtrl-sci cond-mat.other physics.comp-ph


title: A phase-stable dual-comb interferometer

abstract: Improvements to dual-comb interferometers will benefit precision spectroscopy and sensing, distance metrology, tomography, telecommunications etc. A specific requirement of such interferometers is to enforce mutual coherence between the two combs over the measurement time. With feed-forward relative stabilization of the carrier-enveloppe offset frequencies, we experimentally realize such mutual coherence over times that exceed 300 seconds, two orders of magnitude longer than state-of-the-art systems. Illustration is given with near-infrared Fourier transform molecular spectroscopy, where two combs of slightly different repetition frequencies replace a scanning two-beam interferometer. Our technique can be implemented with any frequency comb generators including microresonators or quantum cascade lasers.


title: Improved convection cooling in steady channel flows

abstract: We find steady channel flows that are locally optimal for transferring heat from fixed-temperature walls, under the constraint of a fixed rate of viscous dissipation (enstrophy = $Pe^2$), also the power needed to pump the fluid through the channel. We generate the optima with net flux as a continuation parameter, starting from parabolic (Poiseuille) flow, the unique optimum at maximum net flux. Decreasing the flux, we eventually reach optimal flows that concentrate the enstrophy in boundary layers of thickness $\sim Pe^{-2/5}$ at the channel walls, and have a uniform flow with speed $\sim Pe^{4/5}$ outside the boundary layers. We explain the scalings using physical arguments with a unidirectional flow approximation, and mathematical arguments using a decoupled approximation. We also show that with channels of aspect ratio (length/height) $L$, the boundary layer thickness scales as $L^{3/5}$ and the outer flow speed scales as $L^{-1/5}$ in the unidirectional approximation. At the Reynolds numbers near the turbulent transition for 2D Poiseuille flow in air, we find a 60\% increase in heat transferred over that of Poiseuille flow.


title: Instability of the standing waves for the nonlinear Klein-Gordon equations in one dimension

abstract: In this paper, we consider the following nonlinear Klein-Gordon equation \begin{align*} \partial_{tt}u-\Delta u+u=|u|^{p-1}u,\qquad t\in \mathbb{R},\ x\in \mathbb{R}^d, \end{align*} with $1<p< 1+\frac{4}{d}$. The equation has the standing wave solutions $u_\omega=e^{i\omega t}\phi_{\omega}$ with the frequency $\omega\in(-1,1)$, where $\phi_{\omega}$ obeys \begin{align*} -\Delta \phi+(1-\omega^2)\phi-\phi^p=0. \end{align*} It was proved by Shatah (1983), and Shatah, Strauss (1985) that there exists a critical frequency $\omega_c\in (0,1)$ such that the standing waves solution $u_\omega$ is orbitally stable when $\omega_c<|\omega|<1$, and orbitally unstable when $|\omega|<\omega_c$. Further, the critical case $|\omega|=\omega_c$ in the high dimension $d\ge 2$ was considered by Ohta, Todorova (2007), who proved that it is strongly unstable, by using the virial identities and the radial Sobolev inequality. The one dimension problem was left after then. In this paper, we consider the one-dimension problem and prove that it is orbitally unstable when $|\omega|=\omega_c$.


title: The small-x gluon from forward charm production: implications for a 100 TeV proton collider

abstract: We review the constraints on the small-x gluon PDF that can be derived by exploiting the forward D meson production data from the LHCb experiment at $\sqrt{s}=5,7$ and 13 TeV. We then discuss the phenomenological implications of the resulting improved small-x gluon for ultra-high energy astrophysics, in particular neutrino telescopes, as well as for the proposed Future Circular Collider (FCC) with $\sqrt{s}=100$ TeV. We illustrate how at the FCC even electroweak scale cross-sections can become sensitive to the small-x region of the quark and gluon PDFs, and then demonstrate how the addition of the LHCb heavy meson production measurements leads to a reduction of PDF uncertainties for various benchmark cross-sections.


title: Vulnerability Assessment of Large-scale Power Systems to False Data Injection Attacks

abstract: This paper studies the vulnerability of large-scale power systems to false data injection (FDI) attacks through their physical consequences. Prior work has shown that an attacker-defender bi-level linear program (ADBLP) can be used to determine the worst-case consequences of FDI attacks aiming to maximize the physical power flow on a target line. This ADBLP can be transformed into a single-level mixed-integer linear program, but it is hard to solve on large power systems due to numerical difficulties. In this paper, four computationally efficient algorithms are presented to solve the attack optimization problem on large power systems. These algorithms are applied on the IEEE 118-bus system and the Polish system with 2383 buses to conduct vulnerability assessments, and they provide feasible attacks that cause line overflows, as well as upper bounds on the maximal power flow resulting from any attack.


title: A critical analysis of resampling strategies for the regularized particle filter

abstract: We analyze the performance of different resampling strategies for the regularized particle filter regarding parameter estimation. We show in particular, building on analytical insight obtained in the linear Gaussian case, that resampling systematically can prevent the filtered density from converging towards the true posterior distribution. We discuss several means to overcome this limitation, including kernel bandwidth modulation, and provide evidence that the resulting particle filter clearly outperforms traditional bootstrap particle filters. Our results are supported by numerical simulations on a linear textbook example, the logistic map and a non-linear plant growth model.


title: Weak-strong uniqueness in fluid dynamics

abstract: We give a survey of recent results on weak-strong uniqueness for compressible and incompressible Euler and Navier-Stokes equations, and also make some new observations. The importance of the weak-strong uniqueness principle stems, on the one hand, from the instances of non-uniqueness for the Euler equations exhibited in the past years; and on the other hand from the question of convergence of singular limits, for which weak-strong uniqueness represents an elegant tool.


title: Stochastic differential games with state constraints and Isaacs equations with nonlinear Neumann problems

abstract: We investigate a two-player zero-sum stochastic differential game problem with the state process being constrained in a connected bounded closed domain, and the cost functional described by the solution of a generalized backward stochastic differential equation (GBSDE for short). We show that the value functions enjoy a (strong) dynamic programming principle, and are the unique viscosity solution of the associated Hamilton-Jacobi-Bellman-Isaacs equations with nonlinear Neumann boundary problems. To obtain the existence for viscosity solutions, we provide a new approach utilizing the representation theorem for generators of the GBSDE, which is proved by a random time change method and is a novel result in its own right.


title: Efficient Dynamic Programming Solution to a Platoon Coordination Merge Problem With Stochastic Travel Times

abstract: The problem of maximizing the probability of two trucks being coordinated to merge into a platoon on a highway is considered. Truck platooning is a promising technology that allows heavy vehicles to save fuel by driving with small automatically controlled inter-vehicle distances. In order to leverage the full potential of platooning, platoons can be formed dynamically en route by small adjustments to their speeds. However, in heavily used parts of the road network, travel times are subject to random disturbances originating from traffic, weather and other sources. We formulate this problem as a stochastic dynamic programming problem over a finite horizon, for which solutions can be computed using a backwards recursion. By exploiting the characteristics of the problem, we derive bounds on the set of states that have to be explored at every stage, which in turn reduces the complexity of computing the solution. Simulations suggest that the approach is applicable to realistic problem instances.


title: Lower bounds on distances between a given quantum channel and certain classes of channels

abstract: The tight, in a sense, lower estimates of diamond-norm distance from a given quantum channel to the sets of degradable, antidegradable and entanglement-breaking channels are obtained via the tight continuity bounds for quantum mutual information and for relative entropy of entanglement in finite-dimensional case. As an auxiliary result there are established lower bounds of trace-norm distance from a given bipartite state to the set of all separable states.


title: Impact of an AGN featureless continuum on estimation of stellar population properties

abstract: The effect of the featureless power-law (PL) continuum of an active galactic nucleus (AGN) on the estimation of physical properties of galaxies with optical population spectral synthesis (PSS) remains largely unknown. With this in mind, we fit synthetic galaxy spectra representing a wide range of galaxy star formation histories (SFHs) and including distinct PL contributions of the form $F_{\nu} \propto \nu^{-\alpha}$ with the PSS code STARLIGHT to study to which extent various inferred quantities (e.g. stellar mass, mean age, and mean metallicity) match the input. The synthetic spectral energy distributions (SEDs) computed with our evolutionary spectral synthesis code include an AGN PL component with $0.5 \leq \alpha \leq 2$ and a fractional contribution $0.2 \leq x_{\mathrm{AGN}} \leq 0.8$ to the monochromatic flux at 4020 \AA. At the empirical AGN detection threshold $x_{\mathrm{AGN}}\simeq 0.26$ that we previously inferred in a pilot study on this subject, our results show that the neglect of a PL component in spectral fitting can lead to an overestimation by $\sim$2 dex in stellar mass and by up to $\sim$1 and $\sim$4 dex in the light- and mass-weighted mean stellar age, respectively, whereas the light- and mass-weighted mean stellar metallicity are underestimated by up to $\sim$0.3 and $\sim$0.6 dex, respectively. Other fitting set-ups including either a single PL or multiple PLs in the base reveal, on average, much lower unsystematic uncertainties of the order of those typically found when fitting purely stellar SEDs with stellar templates, however, reaching locally up to $\sim$1, 3 and 0.4 dex in mass, age and metallicity, respectively. Our results underscore the importance of an accurate modelling of the AGN spectral contribution in PSS fits as a minimum requirement for the recovery of the physical and evolutionary properties of stellar populations in active galaxies.


title: Universality and tails of long range interactions in one dimension

abstract: Long-range interactions and, in particular, two-body potentials with power-law long-distance tails are ubiquitous in nature. For two bosons or fermions in one spatial dimension, the latter case being formally equivalent to three-dimensional $s$-wave scattering, we show how generic asymptotic interaction tails can be accounted for in the long-distance limit of scattering wave functions. This is made possible by introducing a generalisation of the collisional phase shifts to include space dependence. We show that this distance dependence is universal, in that it does not depend on short-distance details of the interaction. The energy dependence is also universal, and is fully determined by the asymptotic tails of the two-body potential. As an important application of our findings, we describe how to eliminate finite-size effects with long-range potentials in the calculation of scattering phase shifts from exact diagonalisation. We show that even with moderately small system sizes it is possible to accurately extract phase shifts that would otherwise be plagued with finite-size errors. We also consider multi-channel scattering, focusing on the estimation of open channel asymptotic interaction strengths via finite-size analysis.


title: Robot Planning with Mathematical Models of Human State and Action

abstract: Robots interacting with the physical world plan with models of physics. We advocate that robots interacting with people need to plan with models of cognition. This writeup summarizes the insights we have gained in integrating computational cognitive models of people into robotics planning and control. It starts from a general game-theoretic formulation of interaction, and analyzes how different approximations result in different useful coordination behaviors for the robot during its interaction with people.


title: Improved Poincar\'e inequalities in fractional Sobolev spaces

abstract: We obtain improved fractional Poincar\'e and Sobolev Poincar\'e inequalities including powers of the distance to the boundary in John, $s$-John domains and H\"older-$\alpha$ domains, and discuss their optimality.


title: Incremental Learning Through Deep Adaptation

abstract: Given an existing trained neural network, it is often desirable to learn new capabilities without hindering performance of those already learned. Existing approaches either learn sub-optimal solutions, require joint training, or incur a substantial increment in the number of parameters for each added domain, typically as many as the original network. We propose a method called \emph{Deep Adaptation Networks} (DAN) that constrains newly learned filters to be linear combinations of existing ones. DANs precisely preserve performance on the original domain, require a fraction (typically 13\%, dependent on network architecture) of the number of parameters compared to standard fine-tuning procedures and converge in less cycles of training to a comparable or better level of performance. When coupled with standard network quantization techniques, we further reduce the parameter cost to around 3\% of the original with negligible or no loss in accuracy. The learned architecture can be controlled to switch between various learned representations, enabling a single network to solve a task from multiple different domains. We conduct extensive experiments showing the effectiveness of our method on a range of image classification tasks and explore different aspects of its behavior.


title: Stimulated Raman Adiabatic control of a nuclear spin in diamond

abstract: Coherent manipulation of nuclear spins is a highly desirable tool for both quantum metrology and quantum computation. However, most of the current techniques to control nuclear spins lack of being fast impairing their robustness against decoherence. Here, based on Stimulated Raman Adiabatic Passage, and its modification including shortcuts to adiabaticity, we present a fast protocol for the coherent manipulation of nuclear spins. Moreover, we show how to initialise a nuclear spin starting from a thermal state, and how to implement Raman control for performing Ramsey spectroscopy to measure the dynamical and geometric phases acquired by nuclear spins.


title: Blind Spots for Direct Detection with Simplified DM Models and the LHC

abstract: Using the existing simplified model framework, we build several dark matter models which have suppressed spin-independent scattering cross section. We show that the scattering cross section can vanish due to interference effects with models obtained by simple combinations of simplified models. For weakly interacting massive particle (WIMP) masses $\gtrsim$10 GeV, collider limits are usually much weaker than the direct detection limits coming from LUX or XENON100. However, for our model combinations, LHC analyses are more competitive for some parts of the parameter space. The regions with direct detection blind spots can be strongly constrained from the complementary use of several Large Hadron Collider (LHC) searches like mono-jet, jets + missing transverse energy, heavy vector resonance searches, etc. We evaluate the strongest limits for combinations of scalar + vector, "squark" + vector, and scalar + "squark" mediator, and present the LHC 14 TeV projections.


title: Generalized expressions for hyper-Rayleigh scattering from isotropic liquids

abstract: Generalized expressions for second harmonic scattering from isotropically distributed liquid molecules are derived for arbitrary scattering angles and polarization states.


title: A Redshift Survey of the Nearby Galaxy Cluster Abell 2199: Comparison of the Spatial and Kinematic Distributions of Galaxies with the Intracluster Medium

abstract: We present the results from an extensive spectroscopic survey of the central region of the nearby galaxy cluster Abell 2199 at $z=0.03$. By combining 775 new redshifts from the MMT/Hectospec observations with the data in the literature, we construct a large sample of 1624 galaxies with measured redshifts at $R<30^\prime$, which results in high spectroscopic completeness at $r_{\rm petro,0}<20.5$ (77%). We use these data to study the kinematics and clustering of galaxies focusing on the comparison with those of the intracluster medium (ICM) from Suzaku X-ray observations. We identify 406 member galaxies of A2199 at $R<30^\prime$ using the caustic technique. The velocity dispersion profile of cluster members appears smoothly connected to the stellar velocity dispersion profile of the cD galaxy. The luminosity function is well fitted with a Schechter function at $M_r<-15$. The radial velocities of cluster galaxies generally agree well with those of the ICM, but there are some regions where the velocity difference between the two is about a few hundred kilometer per second. The cluster galaxies show a hint of global rotation at $R<5^\prime$ with $v_{\rm rot}=300{-}600\,\textrm{km s}^{-1}$, but the ICM in the same region do not show such rotation. We apply a friends-of-friends algorithm to the cluster galaxy sample at $R<60^\prime$ and identify 32 group candidates, and examine the spatial correlation between the galaxy groups and X-ray emission. This extensive survey in the central region of A2199 provides an important basis for future studies of interplay among the galaxies, the ICM and the dark matter in the cluster.


title: Measurement of the properites of cosmic rays with the LOFAR radio telescope

abstract: High-energy cosmic rays, impinging on the atmosphere of the Earth initiate cascades of secondary particles, the extensive air showers. The electrons and positrons in the air shower emit electromagnetic radiation. This emission is detected with the LOFAR radio telescope in the frequency range from 10 to 240 MHz. The data are used to determine the properties of the incoming cosmic rays. The radio technique is now routinely used to measure the arrival direction, the energy, and the particle type (atomic mass) of cosmic rays in the energy range from $10^{17}$ to $10^{18}$ eV. This energy region is of particular astrophysical interest, since in this regime a transition from a Galactic to an extra-galactic origin of cosmic rays is expected. For illustration, the LOFAR results are used to set constraints on models to describe the origin of high-energy cosmic rays.


title: Generalized model for the diffusion of solvents in glassy polymers: From Fickian to super Case II

abstract: The diffusion of small solvent molecules in glassy polymers may take on a variety of different forms. Fickian, anomalous, Case II, and super Case II diffusion have all been observed, and theoretical models exist that describe each specific type of behavior. Here we present a single generalized kinetic model capable of yielding all these different types of diffusion on the basis of just two parameters. The principal determinant of the type of diffusion is observed to be a dimensionless parameter, $\gamma$, that describes the influence of solvent-induced swelling in lowering the potential barriers separating available solvent sites. A second parameter, $\eta$, which characterizes the effect of the solvent in reducing the potential energy of a solvent molecule when at rest at an available site, only influences the type of diffusion to a lesser extent. The theoretical analysis does not include any effects that are explicitly non-local in time, an example of which is the inclusion of polymer viscosity in the Thomas-Windle model; it thus represents a variant of Fick's second law utilizing a concentration-dependent diffusivity when $\eta$ is small. To check the significance of time-delayed swelling, a simulation was performed of a modified model that contained a history-dependent term. The results were found to be very similar to those in the time-local model.


title: Optical and X-ray luminosity of expanding nebulae around ultraluminous X-ray sources

abstract: We have performed a set of simulations of expanding, spherically symmetric nebulae inflated by winds from accreting black holes in ultraluminous X-ray sources (ULXs). We implemented a realistic cooling function to account for free-free and bound-free cooling. For all model parameters we considered, the forward shock in the interstellar medium becomes radiative at a radius $\sim $ 100 pc. The emission is primarily in the optical and UV, and the radiative luminosity is about 50% of the total kinetic luminosity of the wind. In contrast, the reverse shock in the wind is adiabatic so long as the terminal outflow velocity of the wind $v_{\rm w} \sim 0.003c$. The shocked wind in these models radiates in X-rays, but with a luminosity of only $\sim 10^{35} \rm\,erg\,s^{-1}$. For wind velocities $v_{\rm w} \sim 0.001c$, the shocked wind becomes radiative, but it is no longer hot enough to produce X-rays. Instead it emits in optical and UV, and the radiative luminosity is comparable to 100% of the wind kinetic luminosity. We suggest that measuring the optical luminosities and putting limits on the X-ray and radio emission from shock-ionized ULX bubbles may help in estimating the mass outflow rate of the central accretion disk and the velocity of the outflow.


title: Controlling scattering of light through topological transitions in all-dielectric metasurfaces

abstract: Topological phase transitions in condensed matter systems have shown extremely rich physics, unveiling such exotic states of matter as topological insulators, superconductors and superfluids. Photonic topological systems open a whole new realm of research exhibiting a number of important distinctions from their condensed matter counterparts. Photonic modes can couple to the continuum of free space modes which makes it feasible to control and manipulate scattering properties of the photonic structure via topology. At the same time, the direct connection of scattering and topological properties of the photonic states allows their probing by spectroscopic means via Fano resonances. Here we demonstrate that the radiative coupling of modes supported by an all-dielectric metasurface can be controlled and tuned under topological phase transitions due to band inversion, correspondingly inducing a distinct switching of the quality factors of the resonances associated with the bands. In addition, we develop a technique to retrieve the topological properties of all-dielectric metasurfaces from the measured far-field scattering characteristics. The collected angle-resolved transmission and reflection spectra allow extracting the momentum-dependent frequencies and lifetimes of the photonic modes. This enables retrieval of the effective photonic Hamiltonian, including the effects of a synthetic gauge field, and topological invariants~-- pseudo-spin Chern numbers. Our results thus open a new avenue to design a new class of metasurfaces with unique scattering characteristics controlled via topological effects. This work also demonstrates how topological states of open systems can be explored via far-field measurements.


title: MOCVD growth of N-polar GaN on on-axis sapphire substrate: impact of AlN nucleation layer on GaN surface hillock density

abstract: We report on the impact of growth conditions on surface hillock density of N-polar GaN grown on nominally on-axis (0001) sapphire substrate by metal organic chemical vapor deposition (MOCVD). Large reduction in hillock density was achieved by implementation of an optimized high temperature AlN nucleation layer and use of indium surfactant in GaN overgrowth. A reduction by more than a factor of five in hillock density from 1000 to 170 hillocks/cm$^{-2}$ was achieved as a result. Crystal quality and surface morphology of the resultant GaN films were characterized by high resolution x-ray diffraction and atomic force microscopy and found to be relatively unaffected by the buffer conditions. It is also shown that the density of smaller surface features is unaffected by AlN buffer conditions.


title: Braided module categories via quantum symmetric pairs

abstract: Let ${\mathfrak g}$ be a finite dimensional complex semisimple Lie algebra. The finite dimensional representations of the quantized enveloping algebra $U_q({\mathfrak g})$ form a braided monoidal category $O_{int}$. We show that the category of finite dimensional representations of a quantum symmetric pair coideal subalgebra $B_{c,s}$ of $U_q({\mathfrak g})$ is a braided module category over an equivariantization of $O_{int}$. The braiding for $B_{c,s}$ is realized by a universal K-matrix which lies in a completion of $B_{c,s}\otimes U_q({\mathfrak g})$. We apply these results to describe a distinguished basis of the center of $B_{c,s}$.


title: Shortcuts to adiabaticity in the presence of a continuum: applications to itinerant quantum state transfer

abstract: We present a method for accelerating adiabatic protocols for systems involving a coupling to a continuum, one that cancels both non-adiabatic errors as well as errors due to dissipation. We focus on applications to a generic quantum state transfer problem, where the goal is to transfer a state between a single level or mode, and a propagating temporal mode in a waveguide or transmission line. Our approach enables perfect adiabatic transfer protocols in this setup, despite a finite protocol speed and a finite waveguide coupling. Our approach even works in highly constrained settings, where there is only a single time-dependent control field available.


title: Resonant structure, formation and stability of the planetary system HD155358

abstract: Two Jovian-sized planets are orbiting the star HD155358 near exact mean motion resonance (MMR) commensurability. In this work we re-analyze the radial velocity (RV) data previously collected by Robertson et al. (2012). Using a Bayesian framework we construct two models - one that includes and one that excludes gravitational planet-planet interactions (PPI). We find that the orbital parameters from our PPI and noPPI models differ by up to $2\sigma$, with our noPPI model being statistically consistent with previous results. In addition, our new PPI model strongly favours the planets being in MMR while our noPPI model strongly disfavours MMR. We conduct a stability analysis by drawing samples from our PPI model's posterior distribution and simulating them for $10^9$ years, finding that our best-fit values land firmly in a stable region of parameter space. We explore a series of formation models that migrate the planets into their observed MMR. We then use these models to directly fit to the observed RV data, where each model is uniquely parameterized by only three constants describing its migration history. Using a Bayesian framework we find that a number of migration models fit the RV data surprisingly well, with some migration parameters being ruled out. Our analysis shows that planet-planet interactions are important to take into account when modelling observations of multi-planetary systems. The additional information that one can gain from interacting models can help constrain planet migration parameters.


title: Distributionally Robust Groupwise Regularization Estimator

abstract: Regularized estimators in the context of group variables have been applied successfully in model and feature selection in order to preserve interpretability. We formulate a Distributionally Robust Optimization (DRO) problem which recovers popular estimators, such as Group Square Root Lasso (GSRL). Our DRO formulation allows us to interpret GSRL as a game, in which we learn a regression parameter while an adversary chooses a perturbation of the data. We wish to pick the parameter to minimize the expected loss under any plausible model chosen by the adversary - who, on the other hand, wishes to increase the expected loss. The regularization parameter turns out to be precisely determined by the amount of perturbation on the training data allowed by the adversary. In this paper, we introduce a data-driven (statistical) criterion for the optimal choice of regularization, which we evaluate asymptotically, in closed form, as the size of the training set increases. Our easy-to-evaluate regularization formula is compared against cross-validation, showing good (sometimes superior) performance.


title: Repeating FRB 121102 : Eight-year Fermi-LAT Upper Limits and Implications

abstract: The repeating fast radio burst (FRB) source that produced FRB 121102 was recently localized in a star forming galaxy at $z=0.193$, which is associated with an extended radio source at the burst location. One possibility is that the repeating FRBs are produced by a new-born magnetar, which also powers the radio nebula. If so, the magnetar may produce $\gamma$-ray emission due to magnetic dipolar spin-down. The luminosity depends on the magnetar spin parameters and age. We process the eight-year Fermi LAT data at the position of FRB 121102 and place an energy flux upper limit of $\sim 10^{-11} \ {\rm erg \ cm^{-2} \ s^{-1}}$ in time bins with six-month intervals, and an accumulated energy flux upper limit of $\sim 4\times 10^{-12} \ {\rm erg \ cm^{-2} \ s^{-1}}$ over the eight year span. The corresponding $\gamma$-ray luminosity upper limits are $\sim 10^{45} \ {\rm erg \ s^{-1}}$ and $\sim 4\times 10^{44} \ {\rm erg \ s^{-1}}$ for the time-resolved and time-integrated analyses, respectively. We discuss the implications of these limits on the young magnetar model.


title: Spectral gap estimates in mean field spin glasses

abstract: We show that mixing for local, reversible dynamics of mean field spin glasses is exponentially slow in the low temperature regime. We introduce a notion of free energy barriers for the overlap, and prove that their existence imply that the spectral gap is exponentially small, and thus that mixing is exponentially slow. We then exhibit sufficient conditions on the equilibrium Gibbs measure which guarantee the existence of these barriers, using the notion of replicon eigenvalue and 2D Guerra Talagrand bounds. We show how these sufficient conditions cover large classes of Ising spin models for reversible nearest-neighbor dynamics and spherical models for Langevin dynamics. Finally, in the case of Ising spins, Panchenko's recent rigorous calculation [79] of the free energy for a system of "two real replica" enables us to prove a quenched LDP for the overlap distribution, which gives us a wider criterion for slow mixing directly related to the Franz-Parisi-Virasoro approach [43,60]. This condition holds in a wider range of temperatures.


title: Mid-infrared observations of O-type stars: spectral morphology

abstract: We present mid-infrared observations for a sample of 16 O-type stars. The data were acquired with the NASA Spitzer Space Telescope, using the IRS instrument at moderate resolution (R $\sim$ 600), covering the range $\sim 10-37$ microns. Our sample includes early, mid and late O supergiants and dwarfs. We explore for the first time their mid-IR spectral morphology in a quantitative way. We use NLTE expanding atmosphere models to help with line identifications, analyze profile contributions and line-formation regions. The O supergiants present a rich emission line spectra. The most intense features are from hydrogen - $6\alpha$, 7$\alpha$, and $8\alpha$ - which have non-negligible contributions of HeI or HeII lines, depending on the spectral type. The spectrum of early O supergiants is a composite of HI and HeII lines, HeI lines being absent. On the other hand, late O supergiants present features composed mainly by HI and HeI lines. All emission lines are formed throughout the stellar wind. We found that O dwarfs exhibit a featureless mid-IR spectrum. Two stars of our sample exhibit very similar mid-IR features, despite having a very different optical spectral classification. The analysis of O-type stars based on mid-IR spectra alone to infer spectral classes or to estimate physical parameters may thus be prone to substantial errors. Our results may therefore inform spectroscopic observations of massive stars located in heavily obscured regions and help establish an initial framework for observations of massive stars using the Mid-Infrared Instrument on the James Webb Space Telescope.


title: Theory of Exciton-Phonon Coupling

abstract: The effects of the electron-phonon interaction on optical excitations can be understood in terms of exciton-phonon coupling, and require a careful treatment in low-dimensional materials with strongly bound excitons or strong electron-hole interaction in general. Through phonon absorption and emission processes, the optically accessible excitons are scattered into otherwise optically dark finite-momentum exciton states. We derive a practical expression for the phonon-induced term of the exciton self-energy (denoted as the exciton-phonon self-energy) that gives the temperature dependence of the optical transition energies and their lifetime broadening resulting from the exciton's interaction with the phonons. We illustrate this theory on a two-dimensional model, and show that our expression for the exciton-phonon self-energy differs qualitatively from previous expressions found in the literature that neglect the exciton binding or electron-hole correlations.


title: Itinerant G-type antiferromagnetic order in SrCr$_2$As$_2$

abstract: Neutron diffraction and magnetic susceptibility studies of a polycrystalline SrCr$_2$As$_2$ sample reveal that this compound is an itinerant G-type antiferromagnet below the N${\rm \acute{e}}$el temperature $T_{\textrm N}$ = 590(5) K with the Cr magnetic moments aligned along the tetragonal $c$ axis. The system remains tetragonal to the lowest measured temperature ($\sim$12 K). The lattice parameter ratio $c/a$ and the magnetic moment saturate at about the same temperature below $\sim$ 200 K, indicating a possible magnetoelastic coupling. The ordered moment, $\mu=1.9(1)~\mu_{\rm B}$/Cr, measured at $T = 12$ K, is significantly reduced compared to its localized value ($4~\mu_{\rm B}$/Cr) due to the itinerant character brought about by the hybridization between the Cr $3d$ and As $4p$ orbitals.


title: Geodesic triangulations exist for cusped Platonic manifolds

abstract: We show that if a cusped hyperbolic manifold is Platonic, i.e., can be decomposed into isometric Platonic solids, it can also be decomposed into geodesic ideal tetrahedra.


title: Newton polyhedra, tropical geometry and the ring of condition for $(C^*)^n$

abstract: The ring of conditions defined by C. De Concini and C. Procesi is an intersection theory for algebraic cycles in a spherical homogeneous space. In the paper we consider the ring of conditions for the group $(C^*)^n$. Up to a big extend this ring can be reduced to the cohomology rings of smooth projective toric varieties. This ring also can be described using tropical geometry. We recall these known results and provide a new description of this ring in terms of convex integral polyhedra.


title: K-sets+: a Linear-time Clustering Algorithm for Data Points with a Sparse Similarity Measure

abstract: In this paper, we first propose a new iterative algorithm, called the K-sets+ algorithm for clustering data points in a semi-metric space, where the distance measure does not necessarily satisfy the triangular inequality. We show that the K-sets+ algorithm converges in a finite number of iterations and it retains the same performance guarantee as the K-sets algorithm for clustering data points in a metric space. We then extend the applicability of the K-sets+ algorithm from data points in a semi-metric space to data points that only have a symmetric similarity measure. Such an extension leads to great reduction of computational complexity. In particular, for an n * n similarity matrix with m nonzero elements in the matrix, the computational complexity of the K-sets+ algorithm is O((Kn + m)I), where I is the number of iterations. The memory complexity to achieve that computational complexity is O(Kn + m). As such, both the computational complexity and the memory complexity are linear in n when the n * n similarity matrix is sparse, i.e., m = O(n). We also conduct various experiments to show the effectiveness of the K-sets+ algorithm by using a synthetic dataset from the stochastic block model and a real network from the WonderNetwork website.


title: Maximal Rank Divisors on $\overline{\mathcal{M}}_{g,n}$

abstract: We compute the class of the effective divisors on $\overline{\mathcal{M}}_{g,n}$, which are set theoretically equal to the locus of moduli points $[C,p_1,\dots ,p_n]$ where $C$ lies on a quadric under the map given by the linear series $|K_C-p_1-\dots -p_n|$. Using this divisor class we show that the moduli spaces $\overline{\mathcal{M}}_{16,8}$ and $\overline{\mathcal{M}}_{17,8}$ are of general type. We also note that the divisor classes computed in the papers \cite{FV1} and \cite{FV2} can be used to show that $\overline{\mathcal{M}}_{12,10}$ is of general type as well.


title: Resolvent estimates for spacetimes bounded by Killing horizons

abstract: We show that the resolvent grows at most exponentially with frequency for the wave equation on a class of stationary spacetimes which are bounded by non-degenerate Killing horizons, without any assumptions on the trapped set. Correspondingly, there exists an exponentially small resonance-free region, and solutions of the Cauchy problem exhibit logarithmic energy decay.


title: Attenuation effect and neutrino oscillation tomography

abstract: Attenuation effect is the effect of weakening of contributions to the oscillation signal from remote structures of matter density profile. The effect is a consequence of integration over the neutrino energy within the energy resolution interval. Structures of a density profile situated at distances larger than the attenuation length, $\lambda_{att}$, are not "seen". We show that the origins of attenuation are (i) averaging of oscillations in certain layer(s) of matter, (ii) smallness of matter effect: $\epsilon \equiv 2EV/\Delta m^2 \ll 1$, where $V$ is the matter potential, and (iii) specific initial and final states on neutrinos. We elaborate on the graphic description of the attenuation which allows us to compute explicitly the effects in the $\epsilon^2$ order for various density profiles and oscillation channels. The attenuation in the case of partial averaging is described. The effect is crucial for interpretation of oscillation data and for the oscillation tomography of the Earth with low energy (solar, supernova, atmospheric, {\it etc.}) neutrinos.


title: Sketching Word Vectors Through Hashing

abstract: We propose a new fast word embedding technique using hash functions. The method is a derandomization of a new type of random projections: By disregarding the classic constraint used in designing random projections (i.e., preserving pairwise distances in a particular normed space), our solution exploits extremely sparse non-negative random projections. Our experiments show that the proposed method can achieve competitive results, comparable to neural embedding learning techniques, however, with only a fraction of the computational complexity of these methods. While the proposed derandomization enhances the computational and space complexity of our method, the possibility of applying weighting methods such as positive pointwise mutual information (PPMI) to our models after their construction (and at a reduced dimensionality) imparts a high discriminatory power to the resulting embeddings. Obviously, this method comes with other known benefits of random projection-based techniques such as ease of update.


title: Community Detection in Signed Networks: an Error-Correcting Code Approach

abstract: In this paper, we consider the community detection problem in signed networks, where there are two types of edges: positive edges (friends) and negative edges (enemies). One renowned theorem of signed networks, known as Harary's theorem, states that structurally balanced signed networks are clusterable. By viewing each cycle in a signed network as a parity-check constraint, we show that the community detection problem in a signed network with two communities is equivalent to the decoding problem for a parity-check code. We also show how one can use two renowned decoding algorithms in error- correcting codes for community detection in signed networks: the bit-flipping algorithm, and the belief propagation algorithm. In addition to these two algorithms, we also propose a new community detection algorithm, called the Hamming distance algorithm, that performs community detection by finding a codeword that minimizes the Hamming distance. We compare the performance of these three algorithms by conducting various experiments with known ground truth. Our experimental results show that our Hamming distance algorithm outperforms the other two.


title: Four-variable expanders over the prime fields

abstract: Let $\mathbb{F}_p$ be a prime field of order $p>2$, and $A$ be a set in $\mathbb{F}_p$ with very small size in terms of $p$. In this note, we show that the number of distinct cubic distances determined by points in $A\times A$ satisfies \[|(A-A)^3+(A-A)^3|\gg |A|^{8/7},\] which improves a result due to Yazici, Murphy, Rudnev, and Shkredov. In addition, we investigate some new families of expanders in four and five variables. We also give an explicit exponent of a problem of Bukh and Tsimerman, namely, we prove that \[\max \left\lbrace |A+A|, |f(A, A)|\right\rbrace\gg |A|^{6/5},\] where $f(x, y)$ is a quadratic polynomial in $\mathbb{F}_p[x, y]$ that is not of the form $g(\alpha x+\beta y)$ for some univariate polynomial $g$.


title: Characterizations of numerical semigroup complements via Ap\'ery sets

abstract: In this paper, we generalize the work of Tuenter to give an identity which completely characterizes the complement of a numerical semigroup in terms of its Ap\'ery sets. Using this result, we compute the $m$th power Sylvester and alternating Sylvester sums for free numerical semigroups. Explicit formulas are given for small $m$.


title: Differential ghost imaging in time domain

abstract: Differential ghost imaging was attempted in time domain, i.e., temporal differential ghost imaging (TDGI), using pseudo-randomized light pulses and a temporal object consisting of no-return-to-zero bit patterns of varying duty. Evaluation of the signal-to-noise characteristics by taking account of errors due to false cross-correlation between the reference and the bucket detector readings indicates that the TDGI outperforms its non-differential counterpart, i.e., time-domain GI, in terms of consistently high and even duty-independent signal-to-noise ratios that are achieved. Dynamic local averaging helps save data recording without compromising the essential features of the TDGI.


title: Probabilistic Image Colorization

abstract: We develop a probabilistic technique for colorizing grayscale natural images. In light of the intrinsic uncertainty of this task, the proposed probabilistic framework has numerous desirable properties. In particular, our model is able to produce multiple plausible and vivid colorizations for a given grayscale image and is one of the first colorization models to provide a proper stochastic sampling scheme. Moreover, our training procedure is supported by a rigorous theoretical framework that does not require any ad hoc heuristics and allows for efficient modeling and learning of the joint pixel color distribution. We demonstrate strong quantitative and qualitative experimental results on the CIFAR-10 dataset and the challenging ILSVRC 2012 dataset.


title: The Fast and the Fiducial: Augmented kludge waveforms for detecting extreme-mass-ratio inspirals

abstract: The extreme-mass-ratio inspirals (EMRIs) of stellar-mass compact objects into massive black holes are an important class of source for the future space-based gravitational-wave detector LISA. Detecting signals from EMRIs will require waveform models that are both accurate and computationally efficient. In this paper, we present the latest implementation of an augmented analytic kludge (AAK) model, publicly available at github.com/alvincjk/EMRI_Kludge_Suite as part of an EMRI waveform software suite. This version of the AAK model has improved accuracy compared to its predecessors, with two-month waveform overlaps against a more accurate fiducial model exceeding 0.97 for a generic range of sources; it also generates waveforms 5-15 times faster than the fiducial model. The AAK model is well suited for scoping out data analysis issues in the upcoming round of mock LISA data challenges. A simple analytic argument shows that it might even be viable for detecting EMRIs with LISA through a semi-coherent template bank method, while the use of the original analytic kludge in the same approach will result in around 90% fewer detections.


title: Highly sensitive atomic based MW interferometry

abstract: We theoretically study a scheme to develop an atomic based MW interferometry using the Rydberg states in Rb. Unlike the traditional MW interferometry, this scheme is not based upon the electrical circuits, hence the sensitivity of the phase and the amplitude/strength of the MW field is not limited by the Nyquist thermal noise. Further this system has great advantage due to its very high bandwidth, ranging from radio frequency (RF), micro wave (MW) to terahertz regime. In addition, this is \textbf{orders of magnitude} more sensitive to field strength as compared to the prior demonstrations on the MW electrometry using the Rydberg atomic states. However previously studied atomic systems are only sensitive to the field strength but not to the phase and hence this scheme provides a great opportunity to characterize the MW completely including the propagation direction and the wavefront. This study opens up a new dimension in the Radar technology such as in synthetic aperture radar interferometry. The MW interferometry is based upon a six-level loopy ladder system involving the Rydberg states in which two sub-systems interfere constructively or destructively depending upon the phase between the MW electric fields closing the loop.


title: Spatio-temporal evolution of the nonresonant instability in shock precursors of young supernova remnants

abstract: A nonresonant cosmic-ray-current-driven instability may operate in the shock precursors of young supernova remnants and be responsible for magnetic-field amplification, plasma heating and turbulence. Earlier simulations demonstrated magnetic-field amplification, and in kinetic studies a reduction of the relative drift between cosmic rays and thermal plasma was observed as backreaction. However, all published simulations used periodic boundary conditions, which do not account for mass conservation in decelerating flows and only allow the temporal development to be studied. Here we report results of fully kinetic Particle-In-Cell simulations with open boundaries that permit inflow of plasma on one side of the simulation box and outflow at the other end, hence allowing an investigation of both the temporal and the spatial development of the instability. Magnetic-field amplification proceeds as in studies with periodic boundaries and, observed here for the first time, the reduction of relative drifts causes the formation of a shock-like compression structure at which a fraction of the plasma ions are reflected. Turbulent electric field generated by the nonresonant instability inelastically scatters cosmic rays, modifying and anisotropizing their energy distribution. Spatial CR scattering is compatible with Bohm diffusion. Electromagnetic turbulence leads to significant nonadiabatic heating of the background plasma maintaining bulk equipartition between ions and electrons. The highest temperatures are reached at sites of large-amplitude electrostatic fields. Ion spectra show supra-thermal tails resulting from stochastic scattering in the turbulent electric field. Together, these modifications in the plasma flow will affect the properties of the shock and particle acceleration there.


title: Hidden Warming Effects of Cloud Cycle Errors in Climate Models

abstract: Clouds' efficiency at reflecting solar radiation and trapping the terrestrial one is strongly modulated by their diurnal cycle. Much attention has been paid to mean cloud properties due to their critical role in climate projections; however, less research has been devoted to their diurnal cycle. Here, we quantify the mean, amplitude, and phase of the cloud cycles in current climate models and compare them with satellite observations and reanalysis data. While the mean values appear to be reliable, the amplitude and phase of the diurnal cycles of clouds show marked inconsistencies, especially over land. We show that, to compensate for the increased net radiation input implied by such errors, an overestimation of the cloud liquid water path may be introduced during calibration of climate models to keep temperatures in line with observations. These discrepancies are likely to be related to cloud parametrization problems in relation to atmospheric convection.


title: Improving Resilience of Autonomous Moving Platforms by Real Time Analysis of Their Cooperation

abstract: Environmental changes, failures, collisions or even terrorist attacks can cause serious malfunctions of the delivery systems. We have presented a novel approach improving resilience of Autonomous Moving Platforms AMPs. The approach is based on multi-level state diagrams describing environmental trigger specifications, movement actions and synchronization primitives. The upper level diagrams allowed us to model advanced interactions between autonomous AMPs and detect irregularities such as deadlocks live-locks etc. The techniques were presented to verify and analyze combined AMPs' behaviors using model checking technique. The described system, Dedan verifier, is still under development. In the near future, a graphical form of verified system representation is planned.


title: N-body simulations of planet formation via pebble accretion I: First Results

abstract: Context. Planet formation with pebbles has been proposed to solve a couple of long-standing issues in the classical formation model. Some sophisticated simulations have been done to confirm the efficiency of pebble accretion. However, there has not been any global N-body simulations that compare the outcomes of planet formation via pebble accretion with observed extrasolar planetary systems. Aims. In this paper, we study the effects of a range of initial parameters of planet formation via pebble accretion, and present the first results of our simulations. Methods. We incorporate the pebble accretion model by Ida et al. (2016) in the N-body code SyMBA (Duncan et al. 1998), along with the effects of gas accretion, eccentricity and inclination damping and planet migration in the disc. Results. We confirm that pebble accretion leads to a variety of planetary systems, but have difficulty in reproducing observed properties of exoplanetary systems, such as planetary mass, semimajor axis, and eccentricity distributions. The main reason behind this is a too-efficient type I migration, which sensitively depends on the disc model. However, our simulations also lead to a few interesting predictions. First, we find that formation efficiencies of planets depend on the stellar metallicities, not only for giant planets, but also for Earths (Es) and Super-Earths (SEs). The dependency for Es/SEs is subtle. Although higher metallicity environments lead to faster formation of a larger number of Es/SEs, they also tend to be lost later via dynamical instability. Second, our results indicate that a wide range of bulk densities observed for Es and SEs is a natural consequence of dynamical evolution of planetary systems. Third, the ejection trend of our simulations suggest that one free-floating E/SE may be expected for two smaller-mass planets.


title: Dark Photon Search in the Mass Range Between 1.5 and 3.4 GeV/$c^2$

abstract: Using a data set of 2.93 fb$^{-1}$ taken at a center-of-mass energy $\sqrt{s}$ = 3.773 GeV with the BESIII detector at the BEPCII collider, we perform a search for an extra U(1) gauge boson, also denoted as a dark photon. We examine the initial state radiation reactions $e^+e^-\rightarrow e^+e^-\gamma_{\rm ISR}$ and $e^+e^-\rightarrow \mu^+\mu^-\gamma_{\rm ISR}$ for this search, where the dark photon would appear as an enhancement in the invariant mass distribution of the leptonic pairs. We observe no obvious enhancement in the mass range between 1.5 and 3.4 GeV/$c^{2}$ and set a 90% confidence level upper limit on the mixing strength of the dark photon and the Standard Model photon. We obtain a competitive limit in the tested mass range.


title: The Character Field Theory and Homology of Character Varieties

abstract: We construct an extended oriented $(2+\epsilon)$-dimensional topological field theory, the character field theory $X_G$ attached to a affine algebraic group in characteristic zero, which calculates the homology of character varieties of surfaces. It is a model for a dimensional reduction of Kapustin-Witten theory ($N=4$ $d=4$ super-Yang-Mills in the GL twist), and a universal version of the unipotent character field theory introduced in arXiv:0904.1247. Boundary conditions in $X_G$ are given by quantum Hamiltonian $G$-spaces, as captured by de Rham (or strong) $G$-categories, i.e., module categories for the monoidal dg category $D(G)$ of $D$-modules on $G$. We show that the circle integral $X_G(S^1)$ (the center and trace of $D(G)$) is identified with the category $D(G/G)$ of "class $D$-modules", while for an oriented surface $S$ (with arbitrary decorations at punctures) we show that $X_G(S)\simeq{\rm H}_*^{BM}(Loc_G(S))$ is the Borel-Moore homology of the corresponding character stack. We also describe the "Hodge filtration" on the character theory, a one parameter degeneration to a TFT whose boundary conditions are given by classical Hamiltonian $G$-spaces, and which encodes a variant of the Hodge filtration on character varieties.


title: A Cascaded Convolutional Neural Network for X-ray Low-dose CT Image Denoising

abstract: Image denoising techniques are essential to reducing noise levels and enhancing diagnosis reliability in low-dose computed tomography (CT). Machine learning based denoising methods have shown great potential in removing the complex and spatial-variant noises in CT images. However, some residue artifacts would appear in the denoised image due to complexity of noises. A cascaded training network was proposed in this work, where the trained CNN was applied on the training dataset to initiate new trainings and remove artifacts induced by denoising. A cascades of convolutional neural networks (CNN) were built iteratively to achieve better performance with simple CNN structures. Experiments were carried out on 2016 Low-dose CT Grand Challenge datasets to evaluate the method's performance.


title: On the order sequence of an embedding of the Ree curve

abstract: In this paper we compute the Weierstrass order-sequence associated with a certain linear series on the Deligne-Lusztig curve of Ree type. As a result, we determine that the set of Weierstrass points of this linear series consists entirely of $\mathbb F_q$-rational points.


title: Positioning for the Internet of Things: A 3GPP Perspective

abstract: Many use cases in the Internet of Things (IoT) will require or benefit from location information, making positioning a vital dimension of the IoT. The 3rd Generation Partnership Project (3GPP) has dedicated a significant effort during its Release 14 to enhance positioning support for its IoT technologies to further improve the 3GPP-based IoT eco-system. In this article, we identify the design challenges of positioning support in Long-Term Evolution Machine Type Communication (LTE-M) and Narrowband IoT (NB-IoT), and overview the 3GPP's work in enhancing the positioning support for LTE-M and NB-IoT. We focus on Observed Time Difference of Arrival (OTDOA), which is a downlink based positioning method. We provide an overview of the OTDOA architecture and protocols, summarize the designs of OTDOA positioning reference signals, and present simulation results to illustrate the positioning performance.


title: The stochastic growth of metal whiskers

abstract: The phenomenon of spontaneously growing metal whiskers (MW) raises significant reliability concerns due to its related arcing and shorting in electric equipment. The growth kinetics of MW remains poorly predictable. Here we present a theory describing the earlier observed intermittent growth of MW as caused by local energy barriers related to variations in the random electric fields generated by surface imperfections. We find the probabilistic distribution of MW stopping times, during which MW growth halts, which is important for reliability projections.


title: Lifting in Besov Spaces

abstract: Let $\Omega$ be a smooth bounded domain in $\mathbb R^n$ and u be a measurable function on $\Omega$ such that $|u(x)|=1$ almost everywhere in $\Omega$. Assume that u belongs to the $B^s_{p,q}(\Omega)$ Besov space. We investigate whether there exists a real-valued function $\varphi \in B^s_{p,q}$ such that $u=e^{i\varphi}$. This extends the corresponding study in Sobolev spaces due to Bourgain, Brezis and the first author. The analysis of this lifting problem leads us to prove some interesting new properties of Besov spaces, in particular a non restriction property when $q>p$.


title: Improved underwater image enhancement algorithms based on partial differential equations (PDEs)

abstract: The experimental results of improved underwater image enhancement algorithms based on partial differential equations (PDEs) are presented in this report. This second work extends the study of previous work and incorporating several improvements into the revised algorithm. Experiments show the evidence of the improvements when compared to previously proposed approaches and other conventional algorithms found in the literature.


title: Dual attainment for the martingale transport problem

abstract: We investigate existence of dual optimizers in one-dimensional martingale optimal transport problems. While [BNT16] established such existence for weak (quasi-sure) duality, [BHP13] showed existence for the natural stronger pointwise duality may fail even in regular cases. We establish that (pointwise) dual maximizers exist when $y\mapsto c(x,y)$ is convex, or equivalent to a convex function. It follows that when marginals are compactly supported, the existence holds when the cost $c(x,y)$ is twice continuously differentiable in $y$. Further, this may not be improved as we give examples with $c(x,\cdot)\in C^{2-\epsilon}$, $\epsilon >0$, where dual attainment fails. Finally, when measures are compactly supported, we show that dual optimizers are Lipschitz if $c$ is Lipschitz.


title: Plants to Polyelectrolytes: Theophylline Polymers and their Microsphere Synthesis

abstract: To extend fossil oil supplies, sustainable feed stocks for the production of useful reagents and polymers should be harnessed. In this regard, chemicals derived from plants are excellent candidates. While the vast majority of plant sources used for polymer science only contain CxHyOz, alkaloids such as caffeine, nicotine, and theophylline possess nitrogen functionality that can provide new functions for bio-derived polymers and their synthesis. In this context, we exploited the chemistry of theophylline, a natural product found in chocolate and tea, to create a cationic poly(theophylline) in a straightforward fashion for the first time. We demonstrate how this new polymer can be synthesized and used for the creation of narrowly disperse cationic microspheres.


title: Main-chain Polyimidazolium Polymers by One-pot Synthesis and Application as Nitrogen-doped Carbon Precursors

abstract: This paper reports on the one-pot synthesis of main-chain imidazolium-containing polymers, some of which show unusually high thermal stability. The imidazolium polymers were obtained by modified Debus-Radziszewski reactions for the chain build-up from simple organic compounds, here pyruvaldehyde, formaldehyde, acetic acid, and a variety of diamines. The reactions were performed in aqueous media at ambient conditions, being synthetically elegant, convenient and highly efficient. Finally, a simple anion-metathesis reaction was conducted to replace the acetate anion with dicyanamide, and the thermal properties of the main-chain polyimidazoliums before and after anion exchange were studied in detail, which demonstrated chain cross-linking by the counterion and a coupled unusually high carbonization yield of up to 66 wt% at 900 oC.


title: Realizable sets of catenary degrees of numerical monoids

abstract: The catenary degree is an invariant that measures the distance between factorizations of elements within an atomic monoid. In this paper, we classify which finite subsets of $\mathbb Z_{\ge 0}$ occur as the set of catenary degrees of a numerical monoid (i.e., a co-finite, additive submonoid of $\mathbb Z_{\ge 0}$). In particular, we show that, with one exception, every finite subset of $\mathbb Z_{\ge 0}$ that can possibly occur as the set of catenary degrees of some atomic monoid is actually achieved by a numerical monoid.


title: Latest Results from T2K

abstract: The T2K (Tokai to Kamioka) experiment is a long-baseline accelerator neutrino oscillation experiment. Intense muon neutrino and antineutrino beams are produced at the J-PARC accelerator complex situated in Tokai. After 280 metres the beam goes through the set of two near detectors: INGRID and ND280, where the beam parameters before oscillations are measured and neutrino cross-sections are studied. After another 295 km, the distance corresponding to the first oscillation maximum, the neutrino beam passes through the Super-Kamiokande far detector. Comparison of the beam characteristics in the near and far detectors makes it possible to determine the oscillation parameters. Based on muon neutrino disappearance, T2K has delivered the world's leading measurement of the $\theta_{23}$ mixing angle. It has been also the first experiment to observe electron neutrino appearance, with a significance of 7.3$\sigma$. These results made it possible to determine the $\theta_{13}$ mixing angle and to provide the first hint of a non-zero value of the $\delta_{CP}$ phase. The T2K experiment has also delivered several neutrino cross-section measurements at neutrino energies around 1 GeV. Currently, T2K is collecting data with a muon antineutrino beam, for the $\delta_{CP}$ and antineutrino cross-section measurements. A summary of the most recent oscillation measurements, including the new electron antineutrino appearance results, is presented.


title: Charting the space of 3D CFTs with a continuous global symmetry

abstract: We study correlation functions of a conserved spin-1 current $J_\mu$ in three dimensional Conformal Field Theories (CFTs). We investigate the constraints imposed by permutation symmetry and current conservation on the form of three point functions $\langle J_\mu J_\nu \mathcal O_{\Delta,\ell}\rangle $ and the four point function $\langle J_\mu J_\nu J_\rho J_\sigma \rangle $ and identify the minimal set of independent crossing symmetry conditions. We obtain a recurrence relation for conformal blocks for generic spin-1 operators in three dimensions. In the process, we improve several technical points, facilitating the use of recurrence relations. By applying the machinery of the numerical conformal bootstrap we obtain universal bounds on the dimensions of certain light operators as well as the central charge. Highlights of our results include numerical evidence for the conformal collider bound and new constraints on the parameters of the critical $O(2)$ model. The results obtained in this work apply to any unitary, parity-preserving three dimensional CFT.


title: Image-based immersed boundary model of the aortic root

abstract: Each year, approximately 300,000 heart valve repair or replacement procedures are performed worldwide, including approximately 70,000 aortic valve replacement surgeries in the United States alone. This paper describes progress in constructing anatomically and physiologically realistic immersed boundary (IB) models of the dynamics of the aortic root and ascending aorta. This work builds on earlier IB models of fluid-structure interaction (FSI) in the aortic root, which previously achieved realistic hemodynamics over multiple cardiac cycles, but which also were limited to simplified aortic geometries and idealized descriptions of the biomechanics of the aortic valve cusps. By contrast, the model described herein uses an anatomical geometry reconstructed from patient-specific computed tomography angiography (CTA) data, and employs a description of the elasticity of the aortic valve leaflets based on a fiber-reinforced constitutive model fit to experimental tensile test data. Numerical tests show that the model is able to resolve the leaflet biomechanics in diastole and early systole at practical grid spacings. The model is also used to examine differences in the mechanics and fluid dynamics yielded by fresh valve leaflets and glutaraldehyde-fixed leaflets similar to those used in bioprosthetic heart valves. Although there are large differences in the leaflet deformations during diastole, the differences in the open configurations of the valve models are relatively small, and nearly identical hemodynamics are obtained in all cases considered.


title: Cofinite Submodule Closed Categories and the Weyl group

abstract: We consider hereditary Artin algebras over arbitrary fields and prove that there is a natural bijection between the Weyl groups and the sets of full additive cofinite submodule closed subcategories of the module categories. While Oppermann, Reiten and Thomas have shown this for algebraically closed fields and finite fields, we give a different method of proof that holds independently of the field. In particular, we show a relatively simple way to construct all modules that contain a given preinjective module as a submodule.


title: SEAGLE: Sparsity-Driven Image Reconstruction under Multiple Scattering

abstract: Multiple scattering of an electromagnetic wave as it passes through an object is a fundamental problem that limits the performance of current imaging systems. In this paper, we describe a new technique-called Series Expansion with Accelerated Gradient Descent on Lippmann-Schwinger Equation (SEAGLE)-for robust imaging under multiple scattering based on a combination of a new nonlinear forward model and a total variation (TV) regularizer. The proposed forward model can account for multiple scattering, which makes it advantageous in applications where linear models are inaccurate. Specifically, it corresponds to a series expansion of the scattered wave with an accelerated-gradient method. This expansion guarantees the convergence even for strongly scattering objects. One of our key insights is that it is possible to obtain an explicit formula for computing the gradient of our nonlinear forward model with respect to the unknown object, thus enabling fast image reconstruction with the state-of-the-art fast iterative shrinkage/thresholding algorithm (FISTA). The proposed method is validated on both simulated and experimentally measured data.


title: Learning to see people like people

abstract: Humans make complex inferences on faces, ranging from objective properties (gender, ethnicity, expression, age, identity, etc) to subjective judgments (facial attractiveness, trustworthiness, sociability, friendliness, etc). While the objective aspects of face perception have been extensively studied, relatively fewer computational models have been developed for the social impressions of faces. Bridging this gap, we develop a method to predict human impressions of faces in 40 subjective social dimensions, using deep representations from state-of-the-art neural networks. We find that model performance grows as the human consensus on a face trait increases, and that model predictions outperform human groups in correlation with human averages. This illustrates the learnability of subjective social perception of faces, especially when there is high human consensus. Our system can be used to decide which photographs from a personal collection will make the best impression. The results are significant for the field of social robotics, demonstrating that robots can learn the subjective judgments defining the underlying fabric of human interaction.


title: Completely $p$-primitive binary quadratic forms

abstract: Let $f(x,y)=ax^2+bxy+cy^2$ be a binary quadratic form with integer coefficients. For a prime $p$ not dividing the discriminant of $f$, we say $f$ is completely $p$-primitive if for any non-zero integer $N$, the diophantine equation $f(x,y)=N$ has always an integer solution $(x,y)=(m,n)$ with $(m,n,p)=1$ whenever it has an integer solution. In this article, we study various properties of completely $p$-primitive binary quadratic forms. In particular, we give a necessary and sufficient condition for a definite binary quadratic form $f$ to be completely $p$-primitive.


title: Dynamical Functional Theory for Compressed Sensing

abstract: We introduce a theoretical approach for designing generalizations of the approximate message passing (AMP) algorithm for compressed sensing which are valid for large observation matrices that are drawn from an invariant random matrix ensemble. By design, the fixed points of the algorithm obey the Thouless-Anderson-Palmer (TAP) equations corresponding to the ensemble. Using a dynamical functional approach we are able to derive an effective stochastic process for the marginal statistics of a single component of the dynamics. This allows us to design memory terms in the algorithm in such a way that the resulting fields become Gaussian random variables allowing for an explicit analysis. The asymptotic statistics of these fields are consistent with the replica ansatz of the compressed sensing problem.


title: Non-commutative Courant algebroids and Quiver algebras

abstract: In this paper, we develop a differential-graded symplectic (Batalin-Vilkovisky) version of the framework of Crawley-Boevey, Etingof and Ginzburg on noncommutative differential geometry based on double derivations to construct non-commutative analogues of the Courant algebroids introduced by Liu, Weinstein and Xu. Adapting geometric constructions of \v{S}evera and Roytenberg for (commutative) graded symplectic supermanifolds, we express the BRST charge, given in our framework by a `homological double derivation', in terms of Van den Bergh's double Poisson algebras for graded bi-symplectic non-commutative 2-forms of weight 1, and in terms of our non-commutative Courant algebroids for graded bi-symplectic non-commutative 2-forms of weight 2 (here, the grading, or ghost degree, is called weight). We then apply our formalism to obtain examples of exact non-commutative Courant algebroids, using appropriate graded quivers equipped with bi-symplectic forms of weight 2, with a possible twist by a closed Karoubi-de Rham non-commutative differential 3-form.


title: Phase recovery and holographic image reconstruction using deep learning in neural networks

abstract: Phase recovery from intensity-only measurements forms the heart of coherent imaging techniques and holography. Here we demonstrate that a neural network can learn to perform phase recovery and holographic image reconstruction after appropriate training. This deep learning-based approach provides an entirely new framework to conduct holographic imaging by rapidly eliminating twin-image and self-interference related spatial artifacts. Compared to existing approaches, this neural network based method is significantly faster to compute, and reconstructs improved phase and amplitude images of the objects using only one hologram, i.e., requires less number of measurements in addition to being computationally faster. We validated this method by reconstructing phase and amplitude images of various samples, including blood and Pap smears, and tissue sections. These results are broadly applicable to any phase recovery problem, and highlight that through machine learning challenging problems in imaging science can be overcome, providing new avenues to design powerful computational imaging systems.

categories: cs.CV cs.IR cs.LG physics.app-ph physics.optics


title: Homodyne monitoring of post-selected decay

abstract: We use homodyne detection to monitor the radiative decay of a superconducting qubit. According to the classical theory of conditional probabilities, the excited state population differs from an exponential decay law if it is conditioned upon a later projective qubit measurement. Quantum trajectory theory accounts for the expectation values of general observables, and we use experimental data to show how a homodyne detection signal is conditioned upon both the initial state and the finally projected state of a decaying qubit. We observe, in particular, how anomalous weak values occur in continuous weak measurement for certain pre- and post-selected states. Subject to homodyne detection, the density matrix evolves in a stochastic manner, but it is restricted to a specific surface in the Bloch sphere. We show that a similar restriction applies to the information associated with the post-selection, and thus bounds the predictions of the theory.


title: Hardware-Software Codesign of Accurate, Multiplier-free Deep Neural Networks

abstract: While Deep Neural Networks (DNNs) push the state-of-the-art in many machine learning applications, they often require millions of expensive floating-point operations for each input classification. This computation overhead limits the applicability of DNNs to low-power, embedded platforms and incurs high cost in data centers. This motivates recent interests in designing low-power, low-latency DNNs based on fixed-point, ternary, or even binary data precision. While recent works in this area offer promising results, they often lead to large accuracy drops when compared to the floating-point networks. We propose a novel approach to map floating-point based DNNs to 8-bit dynamic fixed-point networks with integer power-of-two weights with no change in network architecture. Our dynamic fixed-point DNNs allow different radix points between layers. During inference, power-of-two weights allow multiplications to be replaced with arithmetic shifts, while the 8-bit fixed-point representation simplifies both the buffer and adder design. In addition, we propose a hardware accelerator design to achieve low-power, low-latency inference with insignificant degradation in accuracy. Using our custom accelerator design with the CIFAR-10 and ImageNet datasets, we show that our method achieves significant power and energy savings while increasing the classification accuracy.


title: Joint Spectrum Allocation and Structure Optimization in Green Powered Heterogeneous Cognitive Radio Networks

abstract: We aim at maximizing the sum rate of secondary users (SUs) in OFDM-based Heterogeneous Cognitive Radio (CR) Networks using RF energy harvesting. Assuming SUs operate in a time switching fashion, each time slot is partitioned into three non-overlapping parts devoted for energy harvesting, spectrum sensing and data transmission. The general problem of joint resource allocation and structure optimization is formulated as a Mixed Integer Nonlinear Programming task which is NP-hard and intractable. Thus, we propose to tackle it by decomposing it into two subproblems. We first propose a sub-channel allocation scheme to approximately satisfy SUs' rate requirements and remove the integer constraints. For the second step, we prove that the general optimization problem is reduced to a convex optimization task. Considering the trade-off among fractions of each time slot, we focus on optimizing the time slot structures of SUs that maximize the total throughput while guaranteeing the rate requirements of both real-time and non-real-time SUs. Since the reduced optimization problem does not have a simple closed-form solution, we thus propose a near optimal closed-form solution by utilizing Lambert-W function. We also exploit iterative gradient method based on Lagrangian dual decomposition to achieve near optimal solutions. Simulation results are presented to validate the optimality of the proposed schemes.