title: Temperature Dependent Thermal Boundary Conductance of Monolayer MoS$_2$ by Raman Thermometry

abstract: The electrical and thermal behavior of nanoscale devices based on two-dimensional (2D) materials is often limited by their contacts and interfaces. Here we report the temperature-dependent thermal boundary conductance (TBC) of monolayer MoS$_2$ with AlN and SiO$_2$, using Raman thermometry with laser-induced heating. The temperature-dependent optical absorption of the 2D material is crucial in such experiments, which we characterize here for the first time above room temperature. We obtain TBC ~ 15 MWm$^-$$^2$K$^-$$^1$ near room temperature, increasing as ~ T$^0$$^.$$^6$$^5$ in the range 300 - 600 K. The similar TBC of MoS$_2$ with the two substrates indicates that MoS$_2$ is the "softer" material with weaker phonon irradiance, and the relatively low TBC signifies that such interfaces present a key bottleneck in energy dissipation from 2D devices. Our approach is needed to correctly perform Raman thermometry of 2D materials, and our findings are key for understanding energy coupling at the nanoscale.


title: Lifshitz transition from valence fluctuations in YbAl3

abstract: In Kondo lattice systems with mixed valence, such as YbAl3, interactions between localized electrons in a partially filled f shell and delocalized conduction electrons can lead to fluctuations between two different valence configurations with changing temperature or pressure. The impact of this change on the momentum-space electronic structure and Fermi surface topology is essential for understanding their emergent properties, but has remained enigmatic due to a lack of appropriate experimental probes. Here by employing a combination of molecular beam epitaxy (MBE) and in situ angle-resolved photoemission spectroscopy (ARPES) we show that valence fluctuations can lead to dramatic changes in the Fermi surface topology, even resulting in a Lifshitz transition. As the temperature is lowered, a small electron pocket in YbAl3 becomes completely unoccupied while the low-energy ytterbium (Yb) 4f states become increasingly itinerant, acquiring additional spectral weight, longer lifetimes, and well-defined dispersions. Our work presents the first unified picture of how local valence fluctuations connect to momentum space concepts including band filling and Fermi surface topology in the longstanding problem of mixed-valence systems.


title: A dense plasma globule in the solar neighborhood

abstract: The radio source J1819+3845 underwent a period of extreme interstellar scintillation between circa 1999 and 2007. The plasma structure responsible for this scintillation was determined to be just $1$-$3\,$pc from the solar system and to posses a density of $n_e\sim 10^2\,$cm$^{-3}$ that is three orders of magnitude higher than the ambient interstellar density (de Bruyn & Macquart 2015). Here we present radio-polarimetric images of the field towards J1819+3845 at wavelengths of 0.2, 0.92 and 2$\,$m. We detect an elliptical plasma globule of approximate size $1^\circ \times \gtrsim 2^\circ$ (major-axis position angle of $\approx -40^\circ$), via its Faraday-rotation imprint ($\approx 15\,$rad$\,$m$^{-2}$) on the diffuse Galactic synchrotron emission. The extreme scintillation of J1819+3845 was most likely caused at the turbulent boundary of the globule (J1819+3845 is currently occulted by the globule). The origin and precise nature of the globule remain unknown. Our observations are the first time plasma structures that likely cause extreme scintillation have been directly imaged.


title: Unusual thickness dependence of exciton characteristics in 2D perovskite quantum wells

abstract: Understanding the nature and energy distribution of optical resonances is of central importance in low-dimensional materials$^{1-4}$ and its knowledge is critical for designing efficient optoelectronic devices. Ruddlesden-Popper halide perovskites are 2D solution-processed quantum wells with a general formula A$_2$A'$_{n-1}$M$_n$X$_{3n+1}$, where optoelectronic properties can be tuned by varying the perovskite layer thickness (n value), and have recently emerged as efficient semiconductors with technologically relevant stability. However, fundamental questions concerning the nature of optical resonances (excitons or free-carriers) and the exciton reduced mass, and their scaling with quantum well thickness remains unresolved. Here, using optical spectroscopy and 60-Tesla magneto-absorption supported by modelling, we unambiguously demonstrate that the optical resonances arise from tightly bound excitons with unexpectedly high exciton reduced mass (0.20 m0) and binding energies varying from 470 meV to 125 meV with increasing thickness from n=1 to 5. Our work demonstrates the dominant role of Coulomb interactions in 2D solution-processed quantum wells and presents unique opportunities for next-generation optoelectronic and photonic devices.


title: Deep Voice 3: Scaling Text-to-Speech with Convolutional Sequence Learning

abstract: We present Deep Voice 3, a fully-convolutional attention-based neural text-to-speech (TTS) system. Deep Voice 3 matches state-of-the-art neural speech synthesis systems in naturalness while training ten times faster. We scale Deep Voice 3 to data set sizes unprecedented for TTS, training on more than eight hundred hours of audio from over two thousand speakers. In addition, we identify common error modes of attention-based speech synthesis networks, demonstrate how to mitigate them, and compare several different waveform synthesis methods. We also describe how to scale inference to ten million queries per day on one single-GPU server.


title: The infrared-radio correlation of spheroid- and disc-dominated star-forming galaxies to z $\sim$ 1.5 in the COSMOS field

abstract: Using infrared data from the Herschel Space Observatory and Karl G. Jansky Very Large Array (VLA) 3 GHz observations in the COSMOS field, we investigate the redshift evolution of the infrared-radio correlation (IRRC) for star-forming galaxies (SFGs) we classify as either spheroid- or disc-dominated based on their morphology. The sample predominantly consists of disc galaxies with stellar mass ${\gtrsim}10^{10}\,M_{\odot}$, and residing on the star-forming main sequence (MS). After the removal of AGN using standard approaches, we observe a significant difference between the redshift-evolution of the median IR/radio ratio $\overline{q}_{\mathrm{TIR}}$ of (i) a sample of ellipticals, plus discs with a substantial bulge component (`spheroid-dominated' SFGs) and, (ii) virtually pure discs and irregular systems (`disc-dominated' SFGs). The spheroid-dominated population follows a declining $\overline{q}_{\mathrm{TIR}}$ vs. $z$ trend similar to that measured in recent evolutionary studies of the IRRC. However, for disc-dominated galaxies, where radio and IR emission should be linked to star formation in the most straightforward way, we measure very little change in $\overline{q}_{\mathrm{TIR}}$. This suggests that low-redshift calibrations of radio emission as an SFR-tracer may remain valid out to at least $z\,{\simeq}\,1\,{-}\,1.5$ for pure star-forming systems. We find that the different redshift-evolution of $q_{\rm TIR}$ for the spheroid- and disc-dominated sample is mainly due to an increasing radio excess for spheroid-dominated galaxies at $z\,{\gtrsim}\,$0.8, hinting at some residual AGN activity in these systems. This finding demonstrates that in the absence of AGN the IRRC is independent of redshift, and that radio observations can therefore be used to estimate SFRs at all redshifts for genuinely star-forming galaxies.


title: Reconstruction of cosmological matter perturbations in Modified Gravity

abstract: The analysis of perturbative quantities is a powerful tool to distinguish between different Dark Energy models and gravity theories degenerated at the background level. In this work, we generalise the integral solution of the matter density contrast for General Relativity gravity to a wide class of Modified Gravity (MG) theories. To calculate this solution is necessary prior knowledge of the Hubble rate, the density parameter at the present epoch ($\Omega_{m0}$) and the functional form of the effective Newton's constant that characterises the gravity theory. We estimate in a model-independent way the Hubble expansion rate by applying a non-parametric reconstruction method to model-independent cosmic chronometer data and high-$z$ quasar data. In order to compare our generalised solution of the matter density contrast, using the non-parametric reconstruction of $H(z)$ from observational data, with purely theoretical one, we choose a parameterisation of the Screened MG and the $\Omega_{m0}$ from WMAP-9 collaborations. Finally, we calculate the growth index for the analysed cases, finding very good agreement between theoretical values and the obtained ones using the approach presented in this work.


title: Locally Optimal Control of Complex Networks

abstract: It has recently been shown that the minimum energy solution of the control problem for a linear system produces a control trajectory that is nonlocal. An issue then arises when the dynamics represents a linearization of the underlying nonlinear dynamics of the system where the linearization is only valid in a local region of the state space. Here we provide a solution to the problem of optimally controlling a linearized system by deriving a time-varying set that represents all possible control trajectories parameterized by time and energy. As long as the control action terminus is defined within this set, the control trajectory is guaranteed to be local. If the desired terminus of the control action is far from the initial state, a series of local control actions can be performed in series, re-linearizing the dynamics at each new position.


title: Branched covers of quasipositive links and L-spaces

abstract: Let $L$ be a oriented link such that $\Sigma_n(L)$, the $n$-fold cyclic cover of $S^3$ branched over $L$, is an L-space for some $n \geq 2$. We show that if either $L$ is a strongly quasipositive link other than one with Alexander polynomial a multiple of $(t-1)^{2g(L) + (|L|-1)}$, or $L$ is a quasipositive link other than one with Alexander polynomial divisible by $(t-1)^{2g_4(L) + (|L|-1)}$, then there is an integer $n(L)$, determined by the Alexander polynomial of $L$ in the first case and the Alexander polynomial of $L$ and the smooth $4$-genus of $L$, $g_4(L)$, in the second, such that $n \leq n(L)$. If $K$ is a strongly quasipositive knot with monic Alexander polynomial such as an L-space knot, we show that $\Sigma_n(K)$ is not an L-space for $n \geq 6$, and that the Alexander polynomial of $K$ is a non-trivial product of cyclotomic polynomials if $\Sigma_n(K)$ is an L-space for some $n = 2, 3, 4, 5$. Our results allow us to calculate the smooth and topological 4-ball genera of, for instance, quasi-alternating quasipositive links. They also allow us to classify strongly quasipositive alternating links and $3$-strand pretzel links.


title: Point Neurons with Conductance-Based Synapses in the Neural Engineering Framework

abstract: The mathematical model underlying the Neural Engineering Framework (NEF) expresses neuronal input as a linear combination of synaptic currents. However, in biology, synapses are not perfect current sources and are thus nonlinear. Detailed synapse models are based on channel conductances instead of currents, which require independent handling of excitatory and inhibitory synapses. This, in particular, significantly affects the influence of inhibitory signals on the neuronal dynamics. In this technical report we first summarize the relevant portions of the NEF and conductance-based synapse models. We then discuss a na\"ive translation between populations of LIF neurons with current- and conductance-based synapses based on an estimation of an average membrane potential. Experiments show that this simple approach works relatively well for feed-forward communication channels, yet performance degrades for NEF networks describing more complex dynamics, such as integration.


title: Verifying Equivalence of Database-Driven Applications

abstract: This paper addresses the problem of verifying equivalence between a pair of programs that operate over databases with different schemas. This problem is particularly important in the context of web applications, which typically undergo database refactoring either for performance or maintainability reasons. While web applications should have the same externally observable behavior before and after schema migration, there are no existing tools for proving equivalence of such programs. This paper takes a first step towards solving this problem by formalizing the equivalence and refinement checking problems for database-driven applications. We also propose a proof methodology based on the notion of bisimulation invariants over relational algebra with updates and describe a technique for synthesizing such bisimulation invariants. We have implemented the proposed technique in a tool called Mediator for verifying equivalence between database-driven applications written in our intermediate language and evaluate our tool on 21 benchmarks extracted from textbooks and real-world web applications. Our results show that the proposed methodology can successfully verify 20 of these benchmarks.


title: Finite element approximation of nonlocal fracture models

abstract: We consider nonlocal nonlinear potentials and estimate the rate of convergence of time stepping schemes to the peridynamic equation of motion. We begin by establishing the existence of $H^2$ solutions over any finite time interval. Here spatial approximation by finite element interpolations are considered. The energy stability of the associated semi-discrete time stepping scheme is established and the approximation of strong and weak formulations of the evolution using FE interpolations of $H^2$ solutions are investigated. The strong and weak form of approximations are shown to converge to the actual solution in the mean square norm at the rate $C_t\Delta t +C_s h^2/\epsilon^2$ where $h$ is the mesh size, $\epsilon$ is the size of nonlocal interaction and $\Delta t$ is the time step. The constants $C_t$ and $C_s$ are independent of $\Delta t$, and $h$. In the absence of nonlinearity a CFL like condition for the energy stability of the central difference time discretization scheme is developed.


title: Employing Fusion of Learned and Handcrafted Features for Unconstrained Ear Recognition

abstract: We present an unconstrained ear recognition framework that outperforms state-of-the-art systems in different publicly available image databases. To this end, we developed CNN-based solutions for ear normalization and description, we used well-known handcrafted descriptors, and we fused learned and handcrafted features to improve recognition. We designed a two-stage landmark detector that successfully worked under untrained scenarios. We used the results generated to perform a geometric image normalization that boosted the performance of all evaluated descriptors. Our CNN descriptor outperformed other CNN-based works in the literature, specially in more difficult scenarios. The fusion of learned and handcrafted matchers appears to be complementary as it achieved the best performance in all experiments. The obtained results outperformed all other reported results for the UERC challenge, which contains the most difficult database nowadays.


title: The Dawn of the Post-Naturalness Era

abstract: In an imaginary conversation with Guido Altarelli, I express my views on the status of particle physics beyond the Standard Model and its future prospects.


title: On the Tur\'an number of some ordered even cycles

abstract: A classical result of Bondy and Simonovits in extremal graph theory states that if a graph on $n$ vertices contains no cycle of length $2k$ then it has at most $O(n^{1+1/k})$ edges. However, matching lower bounds are only known for $k=2,3,5$. In this paper we study ordered variants of this problem and prove some tight estimates for a certain class of ordered cycles that we call bordered cycles. In particular, we show that the maximum number of edges in an ordered graph avoiding bordered cycles of length at most $2k$ is $\Theta(n^{1+1/k})$. Strengthening the result of Bondy and Simonovits in the case of 6-cycles, we also show that it is enough to forbid these bordered orderings of the 6-cycle to guarantee an upper bound of $O(n^{4/3})$ on the number of edges.


title: Entropy of real rational surface automorphisms

abstract: We compare real and complex dynamics for automorphisms of rational surfaces that are obtained by lifting \chg{some} quadratic birational maps of the plane. In particular, we show how to exploit the existence of an invariant cubic curve to understand how the real part of an automorphism acts on homology. We apply this understanding to give examples where the entropy of the full (complex) automorphism is the same as its real restriction. Conversely and by different methods, we exhibit different examples where the entropy is strictly decreased by restricting to the real part of the surface. Finally, we give an example of a rational surface automorphism with positive entropy whose periodic cycles are all real.


title: On the Relative Projective Space

abstract: Let $(\mathcal C,\otimes,\mathbb 1)$ be an abelian symmetric monoidal category satisfying certain exactness conditions. In this paper we define a presheaf $\mathbb P^{n}_{\mathcal C}$ on the category of commutative algebras in $\mathcal C$ and we prove that this functor is a $\mathcal C$-scheme in the sense of Toen and Vaqui\'e. This construction gives us a context of non-associative relative algebraic geometry. The most important example of the construction is the octonionic projective space.


title: The slowly variable star FY Lacertae

abstract: Photometric observations and analysis of FY Lac were performed in order to determine the variability characteristics, using measurements taken at Loiano site of Bologna Astronomical Observatory and AAVSO data. The star shows a Long Term variability with a preliminary period of about 274 +/- 28 days and photometric characteristics compatible with a red M5 III giant star, with a temperature of 3420 K degrees.


title: Endpoint $L^p \to L^q$ bounds for integration along certain polynomial curves

abstract: We establish strong-type endpoint $L^p(\mathbb R^d) \to L^q(\mathbb R^d)$ bounds for the operator given by convolution with affine arclength measure on polynomial curves for $d \geq 4$. The bounds established depend only on the dimension $d$ and the degree of the polynomial.


title: Black hole stability under odd-parity perturbations in Horndeski gravity

abstract: We study the stability under linear odd-parity perturbations of static spherically symmetric black holes in Horndeski gravity. We derive the master equation for these perturbations and obtain the conditions of no-ghost and Laplacian instability. In order for the black hole solutions to be stable, we study their generalized "Regge-Wheeler potential". It turns out that the problem is reduced to an algebraic problem where three functions characterizing the black hole should be positive outside the horizon to prove the stability. We found that these conditions are similar to the no-ghost and Laplacian instability conditions. We apply our results to various known solutions.


title: Search for nucleon decays with EXO-200

abstract: A search for instability of nucleons bound in $^{136}$Xe nuclei is reported with 223 kg$\cdot$yr exposure of $^{136}$Xe in the EXO-200 experiment. Lifetime limits of 3.3$\times 10^{23}$ and 1.9$\times 10^{23}$ yrs are established for nucleon decay to $^{133}$Sb and $^{133}$Te, respectively. These are the most stringent to date, exceeding the prior decay limits by a factor of 9 and 7, respectively.


title: Statistics of zero crossings in rough interfaces with fractional elasticity

abstract: We study numerically the distribution of zero crossings in one-dimensional elastic interfaces described by an overdamped Langevin dynamics with periodic boundary conditions. We model the %restoring elastic forces with a Riesz-Feller fractional Laplacian of order $z = 1 + 2\zeta$, such that the interfaces spontaneously relax, with a dynamical exponent $z$, to a self-affine geometry with roughness exponent $\zeta$. By continuously increasing from $\zeta=-1/2$ (macroscopically flat interface described by independent Ornstein--Uhlenbeck processes) to $\zeta=3/2$ (super-rough Mullins--Herring interface), three different regimes are identified: (I) $-1/2<\zeta<0$, (II) $0<\zeta<1$, and (III) $1<\zeta<3/2$. Starting from a flat initial condition, the mean number of zeros of the discretized interface (I) decays exponentially in time and reaches an extensive value in the system size, or decays as a power-law towards (II) a sub-extensive or (III) an intensive value. In the steady-state, the distribution of intervals between zeros changes from an exponential decay in (I) to a power-law decay $P(\ell) \sim \ell^{-\gamma}$ in (II) and (III). While in (II) $\gamma=1-\theta$ with $\theta=1-\zeta$ the steady-state persistence exponent, in (III) we obtain $\gamma = 3-2\zeta$, different from the exponent $\gamma=1$ expected from the prediction $\theta=0$ for infinite super-rough interfaces with $\zeta>1$. The effect on $P(\ell)$ of short-scale smoothening is also analyzed numerically and analytically. A tight relation between the mean interval, the mean width of the interface and the density of zeros is also reported. The results drawn from our analysis of rough interfaces subject to particular boundary conditions or constraints, along with discretization effects, are relevant for the practical analysis of zeros in interface imaging experiments or in numerical analysis.


title: Optimal cutting planes from the group relaxations

abstract: We study quantitative criteria for evaluating the strength of valid inequalities for Gomory and Johnson's finite and infinite group models and we describe the valid inequalities that are optimal for these criteria. We justify and focus on the criterion of maximizing the volume of the nonnegative orthant cut off by a valid inequality. For the finite group model of prime order, we show that the unique maximizer is an automorphism of the {\em Gomory Mixed-Integer (GMI) cut} for a possibly {\em different} finite group problem of the same order. We extend the notion of volume of a simplex to the infinite dimensional case. This is used to show that in the infinite group model, the GMI cut maximizes the volume of the nonnegative orthant cut off by an inequality.


title: $L^p$ improving multilinear Radon-like transforms

abstract: We characterize (up to endpoints) the $k$-tuples $(p_1,\ldots,p_k)$ for which certain $k$-linear generalized Radon transforms map $L^{p_1} \times \cdots \times L^{p_k}$ boundedly into $\mathbb R$. This generalizes a result of Tao and Wright.


title: Spin relaxation of a donor electron coupled to interface states

abstract: An electron spin qubit in a silicon donor atom is a promising candidate for quantum information processing because of its long coherence time. To be sensed with a single-electron transistor, the donor atom is usually located near an interface, where the donor states can be coupled with interface states. Here we study the phonon-assisted spin-relaxation mechanisms when a donor is coupled to confined (quantum-dot-like) interface states. We find that both Zeeman interaction and spin-orbit interaction can hybridize spin and orbital states, each contributing to phonon-assisted spin relaxation in addition to the spin relaxation for a bulk donor or a quantum dot. When the applied magnetic field $B$ is weak (compared to orbital spacing), the phonon assisted spin relaxation shows the $B^5$ dependence. We find that there are peaks (hot-spots) in the $B$-dependent and detuning dependent spin relaxation due to strong hybridization of orbital states with opposite spin. We also find spin relaxation dips (cool-spots) due to the interference of different relaxation channels. Qubit operations near spin relaxation hot-spots can be useful for the fast spin initialization and near cool-spots for the preservation of quantum information during the transfer of spin qubits.


title: Uniform bounds for convolution and restricted X-ray transforms along degenerate curves

abstract: We establish endpoint Lebesgue space bounds for convolution and restricted X-ray transforms along curves satisfying fairly minimal differentiability hypotheses, with affine and Euclidean arclengths. We also explore the behavior of certain natural interpolants and extrapolants of the affine and Euclidean versions of these operators.


title: Partial inversion of the 2D attenuated $X$-ray transform with data on an arc

abstract: In two dimensions, we consider the problem of inversion of the attenuated $X$-ray transform of a compactly supported function from data restricted to lines leaning on a given arc. We provide a method to reconstruct the function on the convex hull of this arc. The attenuation is assumed known. The method of proof uses the Hilbert transform associated with $A$-analytic functions in the sense of Bukhgeim.


title: Uniform $L^p$-improving for weighted averages on curves

abstract: We define variable parameter analogues of the affine arclength measure on curves and prove near-optimal $L^p$-improving estimates for associated multilinear generalized Radon transforms. Some of our results are new even in the convolution case.


title: $\mathcal{P}_m$ Interior Penalty Nonconforming Finite Element Methods for $2m$-th Order PDEs in $\mathbb{R}^n$

abstract: In this work, we propose a family of interior penalty nonconforming finite element methods for $2m$-th order partial differential equations in $\mathbb{R}^n$, for any $m \geq 0$, $n \geq 1$. For the nonconforming finite element, the shape function space consists of all polynomials with a degree not greater than $m$ and is hence minimal. This family of finite element spaces has some natural inclusion properties as in the corresponding Sobolev spaces in the continuous cases. By applying the interior penalty to the bilinear form, we establish quasi-optimal error estimates in the energy norm, provided that the corresponding conforming relatives exist. These theoretical results are further validated by the numerical tests.


title: Nonparametrically estimating dynamic bivariate correlation using visibility graph algorithm

abstract: Dynamic conditional correlation (DCC) is a method that estimates the correlation between two time series across time. Although used primarily in finance so far, DCC has been proposed recently as a model-based estimation method for quantifying functional connectivity during fMRI experiments. DCC could also be used to estimate the dynamic correlation between other types of time series such as local field potentials (LFP's) or spike trains recorded from distinct brain areas. DCC has very nice properties compared to other existing methods, but its applications for neuroscience are currently limited because of non-optimal performance in the presence of outliers. To address this issue, we developed a robust nonparametric version of DCC, based on an adaptation of the weighted visibility graph algorithm which converts a time series into a weighted graph. The modified DCC demonstrated better performance in the analysis of empirical data sets: one fMRI data set collected from a human subject performing a Stroop task; and one LFP data set recorded from an awake rat in resting state. Nonparametric DCC has the potential of enlarging the spectrum of analytical tools designed to assess the dynamic coupling and uncoupling of activity among brain areas.


title: Charged wormholes in f(R,T) extended theory of gravity

abstract: Wormholes are a solution for General Relativity field equations which characterize a passage or a tunnel that connects two different regions of space-time and is filled by some sort of exotic matter, that does not satisfy the energy conditions. On the other hand, it is known that in extended theories of gravity, the extra degrees of freedom once provided may allow the energy conditions to be obeyed and, consequently, the matter content of the wormhole to be non-exotic. In this work, we obtain, as a novelty in the literature, solutions for charged wormholes in the $f(R,T)$ extended theory of gravity. We show that the presence of charge in these objects may be a possibility to respect some stability conditions for their metric. Also, remarkably, the energy conditions are respected in the present approach.


title: Bulk-boundary correspondance for Sturmian Kohmoto like models

abstract: We consider one dimensional tight binding models on $\ell^2(\mathbb Z)$ whose spatial structure is encoded by a Sturmian sequence $(\xi_n)_n\in \{a,b\}^\mathbb Z$. An example is the Kohmoto Hamiltonian, which is given by the discrete Laplacian plus an onsite potential $v_n$ taking value $0$ or $1$ according to whether $\xi_n$ is $a$ or $b$. The only non-trivial topological invariants of such a model are its gap-labels. The bulk-boundary correspondence we establish here states that there is a correspondence between the gap label and a winding number associated to the edge states, which arises if the system is augmented and compressed onto half space $\ell^2(\mathbb N)$. This has been experimentally observed with polaritonic waveguides. A correct theoretical explanation requires, however, first a smoothing out of the atomic motion via phason flips. With such an interpretation at hand, the winding number corresponds to the mechanical work through a cycle which the atomic motion exhibits on the edge states.


title: Uniform estimates for Fourier restriction to polynomial curves in $\mathbb R^d$

abstract: We prove uniform $L^p \to L^q$ bounds for Fourier restriction to polynomial curves in $\mathbb R^d$ with affine arclength measure, in the conjectured range.


title: Linear and bilinear restriction to certain rotationally symmetric hypersurfaces

abstract: Conditional on Fourier restriction estimates for elliptic hypersurfaces, we prove optimal restriction estimates for polynomial hypersurfaces of revolution for which the defining polynomial has non-negative coefficients. In particular, we obtain uniform--depending only on the dimension and polynomial degree--estimates for restriction with affine surface measure, slightly beyond the bilinear range. The main step in the proof of our linear result is an (unconditional) bilinear adjoint restriction estimate for pieces at different scales.


title: The observing campaign on the deep-space debris WT1190F as a test case for short-warning NEO impacts

abstract: On 2015 November 13, the small artificial object designated WT1190F entered the Earth atmosphere above the Indian Ocean offshore Sri Lanka after being discovered as a possible new asteroid only a few weeks earlier. At ESA's SSA-NEO Coordination Centre we took advantage of this opportunity to organize a ground-based observational campaign, using WT1190F as a test case for a possible similar future event involving a natural asteroidal body.


title: Soft collinear effective theory for gravity

abstract: We present how to construct a Soft Collinear Effective Theory (SCET) for gravity at the leading and next-to-leading powers from the ground up. The soft graviton theorem and decoupling of collinear gravitons at the leading power are manifest from the outset in the effective symmetries of the theory. At the next-to-leading power, certain simple structures of amplitudes, which are completely obscure in Feynman diagrams of the full theory, are also revealed, which greatly simplifies calculations. The effective lagrangian is highly constrained by effectively multiple copies of diffeomorphism invariance that are inevitably present in gravity SCET due to mode separation, an essential ingredient of any SCET. Further explorations of effective theories of gravity with mode separation may shed light on lagrangian-level understandings of some of the surprising properties of gravitational scattering amplitudes. A gravity SCET with an appropriate inclusion of Glauber modes may serve as a powerful tool for studying gravitational scattering in the Regge limit.


title: Scale invariant Fourier restriction to a hyperbolic surface

abstract: This result sharpens the bilinear to linear deduction of Lee and Vargas for extension estimates on the hyperbolic paraboloid in $\mathbb R^3$ to the sharp line, leading to the first scale-invariant restriction estimates, beyond the Stein--Tomas range, for a hypersurface on which the principal curvatures have different signs.


title: Annama H chondrite - mineralogy, physical properties, cosmic ray exposure, and parent body history

abstract: The fall of the Annama meteorite occurred early morning (local time) on April 19, 2014 on the Kola Peninsula (Russia). Based on mineralogy and physical properties, Annama is a typical H chondrite. It has a high Ar-Ar age of 4.4 Ga. Its cosmic ray exposure history is atypical as it is not part of the large group of H chondrites with a prominent 7 - 8 Ma peak in the exposure age histograms. Instead, its exposure age is within uncertainty of a smaller peak at 30 \pm 4 Ma. The results from short-lived radionuclides are compatible with an atmosperic pre-entry radius of 30 - 40 cm. However, based on noble gas and cosmogenic radionuclide data, Annama must have been part of a larger body (radius >65 cm) for a large part of its cosmic ray exposure history. The 10Be concentration indicates a recent (3 - 5 Ma) breakup which may be responsible for the Annama parent body size reduction to 30 - 35 cm pre-entry radius.


title: Endpoint Lebesgue estimates for weighted averages on polynomial curves

abstract: We establish optimal Lebesgue estimates for a class of generalized Radon transforms defined by averaging functions along polynomial-like curves. The presence of an essentially optimal weight allows us to prove uniform estimates, wherein the Lebesgue exponents are completely independent of the curves and the operator norms depend only on the polynomial degree. Moreover, our weighted estimates possess rather strong diffeomorphism invariance properties, allowing us to obtain uniform bounds for averages on curves satisfying a natural nilpotency hypothesis.


title: Direct photon spectrum and elliptic flow produced from Pb+Pb collisions at $\sqrt{s_{NN}}=2.76$ TeV at the CERN Large Hadron Collider within an integrated hydrokinetic model

abstract: The photon transverse momentum spectrum and its anisotropy from Pb+Pb collisions at the CERN Large Hadron Collider energy $\sqrt {s_{NN}}=2.76$ TeV are investigated within the integrated hydrokinetic model (iHKM). Photon production is accumulated from the different processes at the various stages of relativistic heavy ion collisions: from the primary hard photons of very early stage of parton collisions to the thermal photons from equilibrated quark-gluon and hadron gas stages. Along the way a hadronic medium evolution is treated in two distinct, in a sense opposite, approaches: chemically equilibrated and chemically frozen system expansion. Studying the centrality dependence of the results obtained allows us to conclude that a relatively strong transverse momentum anisotropy of thermal radiation is suppressed by prompt photon emission which is an isotropic. We find out that this effect is getting stronger as centrality increases because of the simultaneous increase in the relative contribution of prompt photons in the soft part of the spectra. The substantial results obtained in iHKM with nonzero viscosity ($\eta/s=0.08$) for photon spectra and $v_2$ coefficients are mostly within the error bars of experimental data, but there is some systematic underestimation of both observables for the near central events. We claim that a situation could be significantly improved if an additional photon radiation that accompanies the presence of a deconfined environment is included. Since a matter of a space-time layer where hadronization takes place is actively involved in anisotropic transverse flow, both positive contributions to the spectra and $v_2$ are considerable, albeit such an argument needs further research and elaboration.


title: Tidal dissipation in rotating fluid bodies: the presence of a magnetic field

abstract: We investigate effects of the presence of a magnetic field on tidal dissipation in rotating fluid bodies. We consider a simplified model consisting of a rigid core and a fluid envelope, permeated by a background magnetic field (either a dipolar field or a uniform axial field). The wavelike tidal responses in the fluid layer are in the form of magnetic-Coriolis waves, which are restored by both the Coriolis force and the Lorentz force. Energy dissipation occurs through viscous damping and Ohmic damping of these waves. Our numerical results show that the tidal dissipation can be dominated by Ohmic damping even with a weak magnetic field. The presence of a magnetic field smooths out the complicated frequency-dependence of the dissipation rate, and broadens the frequency spectrum of the dissipation rate, depending on the strength of the background magnetic field. However, the frequency-averaged dissipation is independent of the strength and structure of the magnetic field, and of the dissipative parameters, in the approximation that the wave-like response is driven only by the Coriolis force acting on the non-wavelike tidal flow. Indeed, the frequency-averaged dissipation quantity is in good agreement with previous analytical results in the absence of magnetic fields. Our results suggest that the frequency-averaged tidal dissipation of the wavelike perturbations is insensitive to detailed damping mechanisms and dissipative properties.


title: The elliptic KZB connection and algebraic de Rham theory for unipotent fundamental groups of elliptic curves

abstract: In this paper, we develop an algebraic de Rham theory for unipotent fundamental groups of once punctured elliptic curves over a field of characteristic zero using the universal elliptic KZB connection of Calaque-Enriquez-Etingof and Levin-Racinet. We use it to give an explicit version of Tannaka duality for unipotent connections over an elliptic curve with a regular singular point at the identity.


title: Existence of extremizers for a model convolution operator

abstract: The operator $T$, defined by convolution with the affine arc length measure on the moment curve parametrized by $h(t)=(t,t^{2},...,t^{d})$ is a bounded operator from $L^{p}$ to $L^{q}$ if $(\frac{1}{p}, \frac{1}{q})$ lies on a line segment. In this article we prove that at non-end points there exist functions which extremize the associated inequality and any extremizing sequence is pre compact modulo the action of the symmetry of $T$. We also establish a relation between extremizers for $T$ at the end points and the extremizers of an X-ray transform restricted to directions along the moment curve. Our proof is based on the ideas of Michael Christ on convolution with the surface measure on the paraboloid.


title: Zero Variance and Hamiltonian Monte Carlo Methods in GARCH Models

abstract: In this paper, we develop Bayesian Hamiltonian Monte Carlo methods for inference in asymmetric GARCH models under different distributions for the error term. We implemented Zero-variance and Hamiltonian Monte Carlo schemes for parameter estimation to try and reduce the standard errors of the estimates thus obtaing more efficient results at the price of a small extra computational cost.


title: The MUSE Hubble Ultra Deep Field Survey: V. Spatially resolved stellar kinematics of galaxies at redshift $0.2\lesssim z \lesssim 0.8$

abstract: We present spatially resolved stellar kinematic maps, for the first time, for a sample of 17 intermediate redshift galaxies (0.2 < z < 0.8). We used deep MUSE/VLT integral field spectroscopic observations in the Hubble Deep Field South (HDFS) and Hubble Ultra Deep Field (HUDF), resulting from ~30h integration time per field, each covering 1'x1' field of view, with ~0.65" spatial resolution. We selected all galaxies brighter than 25mag in the I band and for which the stellar continuum is detected over an area that is at least two times larger than the spatial resolution. The resulting sample contains mostly late-type disk, main-sequence star-forming galaxies with 10^8.5 - 10^10.5 Msun. Using a full-spectrum fitting technique, we derive two-dimensional maps of the stellar and gas kinematics, including the radial velocity V and velocity dispersion sigma. We find that most galaxies in the sample are consistent with having rotating stellar disks with roughly constant velocity dispersions and that the Vrms=sqrt{V^2+sigma^2} of the gas and stars, a scaling proxy for the galaxy gravitational potential, compare well to each other. These spatially resolved observations of intermediate redshift galaxies suggest that the regular stellar kinematics of disk galaxies that is observed in the local Universe was already in place 4 - 7 Gyr ago and that their gas kinematics traces the gravitational potential of the galaxy, thus is not dominated by shocks and turbulent motions. Finally, we build dynamical axisymmetric Jeans models constrained by the derived stellar kinematics for two specific galaxies and derive their dynamical masses. These are in good agreement (within 25%) with those derived from simple exponential disk models based on the gas kinematics. The obtained mass-to-light ratios hint towards dark matter dominated systems within a few effective radii.


title: Verb Pattern: A Probabilistic Semantic Representation on Verbs

abstract: Verbs are important in semantic understanding of natural language. Traditional verb representations, such as FrameNet, PropBank, VerbNet, focus on verbs' roles. These roles are too coarse to represent verbs' semantics. In this paper, we introduce verb patterns to represent verbs' semantics, such that each pattern corresponds to a single semantic of the verb. First we analyze the principles for verb patterns: generality and specificity. Then we propose a nonparametric model based on description length. Experimental results prove the high effectiveness of verb patterns. We further apply verb patterns to context-aware conceptualization, to show that verb patterns are helpful in semantic-related tasks.


title: Quantum dynamics from a numerical linked cluster expansion

abstract: We demonstrate that a numerical linked cluster expansion method is a powerful tool to calculate quantum dynamics. We calculate the dynamics of the magnetization and spin correlations in the two-dimensional transverse field Ising and XXZ models evolved from a product state. Such dynamics are directly probed in ongoing experiments in ultracold atoms, molecules, and ions. We show that a numerical linked cluster expansion gives dramatically more accurate results at short-to-moderate times than exact diagonalization, and simultaneously requires fewer computational resources. More specifically, the cluster expansion frequently produces more accurate results than an exact diagonalization calculation that would require $10^{5}$--$10^{10}$ more computational operations and memory.


title: Radiative Heat Transfer Between Core-Shell Nanoparticles

abstract: Radiative heat transfer in systems with core-shell nanoparticles may exhibit not only a combination of disparate physical properties of its components but also further enhanced properties that arise from the synergistic properties of the core and shell components. We study the thermal conductance between two core-shell nanoparticles (CSNPs). The contribution of electric and magnetic dipole moments to the thermal conductance depend sensitively on the core and shell materials, and adjustable by core size and shell thickness. We predict that the radiative heat transfer in a dimer of Au@SiO2 CSNPs (i.e., silica-coated gold nanoparticles) could be enhanced several order of magnitude compared to bare Au nanoparticles. However, the reduction of several orders of magnitude in the heat transfer is possible between SiO2@Au CSNPs (i.e., silica as a core and gold as a shell) than that of uncoated SiO2 nanoparticles.


title: Shear-driven parametric instability in a precessing sphere

abstract: The present numerical study aims at shedding light on the mechanism underlying the precessional instability in a sphere. Precessional instabilities in the form of parametric resonance due to topographic coupling have been reported in a spheroidal geometry both analytically and numerically. We show that such parametric resonances can also develop in spherical geometry due to the conical shear layers driven by the Ekman pumping singularities at the critical latitudes. Scaling considerations lead to a stability criterion of the form, $|P_o|>O(E^{4/5})$, where $P_o$ represents the Poincar\'e number and $E$ the Ekman number. The predicted threshold is consistent with our numerical simulations as well as previous experimental results. When the precessional forcing is supercriticial, our simulations show evidence of an inverse cascade, i.e. small scale flows merging into large scale cyclones with a retrograde drift. Finally, it is shown that this instability mechanism may be relevant to precessing celestial bodies such as the Earth and Earth's moon.


title: An Ambarzumian type theorem on graphs with odd cycles

abstract: We consider an inverse problem for Schr\"odinger operators on a connected equilateral graph $G$ with standard matching conditions. The graph $G$ consists of at least two odd cycles glued together at a common vertex. We prove an Ambarzumian type result, i.e., if a specific part of the spectrum is the same as in the case of zero potential, then the potential has to be zero.


title: On the sharp upper bound related to the weak Muckenhoupt-Wheeden conjecture

abstract: We construct an example showing that the upper bound $[w]_{A_1}\log({\rm{e}}+[w]_{A_1})$ for the $L^1(w)\to L^{1,\infty}(w)$ norm of the Hilbert transform cannot be improved in general.


title: On production and asymmetric focusing of flat electron beams using rectangular capillary discharge plasmas

abstract: A method for the asymmetric focusing of electron bunches, based on the active plasma lensing technique is proposed. This method takes advantage of the strong inhomogeneous magnetic field generated inside the capillary discharge plasma to focus the ultrarelativistic electrons. The plasma and magnetic field parameters inside the capillary discharge are described theoretically and modeled with dissipative magnetohydrodynamic computer simulations enabling analysis of the capillaries of rectangle cross-sections. Large aspect ratio rectangular capillaries might be used to transport electron beams with high emittance asymmetries, as well as assist in forming spatially flat electron bunches for final focusing before the interaction point.


title: On the Consistency of Graph-based Bayesian Learning and the Scalability of Sampling Algorithms

abstract: A popular approach to semi-supervised learning proceeds by endowing the input data with a graph structure in order to extract geometric information and incorporate it into a Bayesian framework. We introduce new theory that gives appropriate scalings of graph parameters that provably lead to a well-defined limiting posterior as the size of the unlabeled data set grows. Furthermore, we show that these consistency results have profound algorithmic implications. When consistency holds, carefully designed graph-based Markov chain Monte Carlo algorithms are proved to have a uniform spectral gap, independent of the number of unlabeled inputs. Several numerical experiments corroborate both the statistical consistency and the algorithmic scalability established by the theory.


title: Collective behavior of light in vacuum

abstract: Under the action of light-by-light scattering, light beams show collective behaviors in vacuum. For instance, in the case of two counterpropagating laser beams with specific initial helicity, the polarization of each beam oscillates periodically between the left and right helicity. Furthermore, the amplitudes and the corresponding intensities of each polarization propagate like waves. Such polarization waves might be observationally accessible in future laser experiments, in a physical regime complementary to those explored by particle accelerators.


title: Predicativity and parametric polymorphism of Brouwerian implication

abstract: A common objection to the definition of intuitionistic implication in the Proof Interpretation is that it is impredicative. I discuss the history of that objection, argue that in Brouwer's writings predicativity of implication is ensured through parametric polymorphism of functions on species, and compare this construal with the alternative approaches to predicative implication of Goodman, Dummett, Prawitz, and Martin-L\"of.


title: Signatures of Parkinson's disease in complexity of low-frequency fluctuations of postural sway velocity

abstract: Posturography is routinely used to qualitatively assess one of the cardinal symptoms of Parkinson's disease -- postural instability. While most measures of balance control are derived from displacement of the center of pressure there is evidence that such control is more likely to be velocity-based. We performed static posturographic tests (eyes open and eyes closed) during quiet standing in narrow stance for n=30 patients with Parkinson's disease (PDP) in the ON state and compared the results with those of n=30 age-matched senior controls (HSC) and n=60 young controls (HYC). We used differentiator filters to generate time series of low-frequency fluctuations of sway velocity and calculated their Lempel-Ziv complexity (LZC). With eyes closed, the mediolateral LZC of HSC 0.21 (0.02) was significantly higher than those of HYC 0.19 (0.02) and PDP 0.18 (002). Thus aging and PD have opposite effects on mediolateral LZC which strongly differentiates between HSC and PDP (92% sensitivity and 87% specificity).


title: Low Precision RNNs: Quantizing RNNs Without Losing Accuracy

abstract: Similar to convolution neural networks, recurrent neural networks (RNNs) typically suffer from over-parameterization. Quantizing bit-widths of weights and activations results in runtime efficiency on hardware, yet it often comes at the cost of reduced accuracy. This paper proposes a quantization approach that increases model size with bit-width reduction. This approach will allow networks to perform at their baseline accuracy while still maintaining the benefits of reduced precision and overall model size reduction.


title: Anisotropic magnetic properties of the triangular plane lattice material TmMgGaO4

abstract: The crystal growth, structure, and basic magnetic properties of TmMgGaO4 are reported. The Tm ions are located in a planar triangular lattice consisting of distorted TmO6 octahedra, while the Mg and Ga atoms randomly occupy intermediary bilayers of M-O triangular bipyramids. The Tm ions are positionally disordered. The material displays an antiferromagnetic Curie Weiss theta of ~ -20 -25 K, with no clear ordering visible in the magnetic susceptibility down to 1.8 K; the structure and magnetic properties suggest that ordering of the magnetic moments is frustrated by both structural disorder and the triangular magnetic motif. Single crystal magnetization measurements indicate that the magnetic properties are highly anisotropic, with large moments measured perpendicular to the triangular planes. At 2 K, a broad step-like feature is seen in the field-dependent magnetization perpendicular to the plane on applied field near 2 Tesla.


title: The Effect of Crystal Symmetries on the Locality of Screw Dislocation Cores

abstract: In linearised continuum elasticity, the elastic strain due to a straight dislocation line decays as $O(r^{-1})$, where $r$ denotes the distance to the defect core. It is shown in Ehrlacher, Ortner, Shapeev (2016) that the core correction due to nonlinear and discrete (atomistic) effects decays like $O(r^{-2})$. In the present work, we focus on screw dislocations under pure anti-plane shear kinematics. In this setting we demonstrate that an improved decay $O(r^{-p})$, $p > 2$, of the core correction is obtained when crystalline symmetries are fully exploited and possibly a simple and explicit correction of the continuum far-field prediction is made. This result is interesting in its own right as it demonstrates that, in some cases, continuum elasticity gives a much better prediction of the elastic field surrounding a dislocation than expected, and moreover has practical implications for atomistic simulation of dislocations cores, which we discuss as well.


title: Solving the "false positives" problem in fraud prediction

abstract: In this paper, we present an automated feature engineering based approach to dramatically reduce false positives in fraud prediction. False positives plague the fraud prediction industry. It is estimated that only 1 in 5 declared as fraud are actually fraud and roughly 1 in every 6 customers have had a valid transaction declined in the past year. To address this problem, we use the Deep Feature Synthesis algorithm to automatically derive behavioral features based on the historical data of the card associated with a transaction. We generate 237 features (>100 behavioral patterns) for each transaction, and use a random forest to learn a classifier. We tested our machine learning model on data from a large multinational bank and compared it to their existing solution. On an unseen data of 1.852 million transactions, we were able to reduce the false positives by 54% and provide a savings of 190K euros. We also assess how to deploy this solution, and whether it necessitates streaming computation for real time scoring. We found that our solution can maintain similar benefits even when historical features are computed once every 7 days.


title: Direct determination of one-dimensional interphase structures using normalized crystal truncation rod analysis

abstract: A one-dimensional non-iterative direct method was employed for normalized crystal truncation rod. The non-iterative approach, utilizing Kramers-Kronig relation, avoids the ambiguities due to the improper initial model or the incomplete convergence in the conventional iterative methods. The validity and limitation of the present method are demonstrated through both numerical simulations and experiments with Pt (111) in 0.1 M CsF aqueous solution. The present method is compared to conventional iterative phase-retrieval methods.


title: Medical Implications of Space Radiation Exposure Due to Low Altitude Polar Orbits

abstract: Space radiation research has progressed rapidly in recent years, but there remain large uncertainties in predicting and extrapolating biological responses to humans. Exposure to cosmic radiation and Solar Particle Events may pose a critical health risk to future spaceflight crews and can have a serious impact to all biomedical aspects of space exploration. The relatively minimal shielding of the cancelled 1960's Manned Orbiting Laboratory program's space vehicle and the high inclination polar orbits would have left the crew susceptible to high exposures of cosmic radiation and high dose-rate SPEs that are mostly unpredictable in frequency and intensity. In this study, we have modeled the nominal and off-nominal radiation environment that a MOL-like spacecraft vehicle would be exposed to during a 30-day mission using high performance, multi-core computers. Projected doses from a historically large SPE (e.g. the August 1972 solar event) have been analyzed in the context of the MOL orbit profile, providing an opportunity to study its impact to crew health and subsequent contingencies.It is reasonable to presume that future commercial, government, and military spaceflight missions in low-Earth orbit will have vehicles with similar shielding and orbital profiles. Studying the impact of cosmic radiation to the mission's operational integrity and the health of MOL crewmembers provides an excellent surrogate and case-study for future commercial and military spaceflight missions.


title: Random Matrix Ensembles For Many-Body Quantum Systems

abstract: Classical random matrix ensembles were originally introduced in physics to approximate quantum many-particle nuclear interactions. However, there exists a plethora of quantum systems whose dynamics is explained in terms of few-particle (predominantly two-particle) interactions. The random matrix models incorporating the few-particle nature of interactions are known as embedded random matrix ensembles. In the present paper, we provide a brief overview of these two ensembles and illustrate how the embedded ensembles can be successfully used to study decoherence of a qubit interacting with an environment, both for fermionic and bosonic embedded ensembles. Numerical calculations show the dependence of decoherence on the nature of the environment.


title: Effect of PYTHIA8 tunes on event shapes and top-quark reconstruction in $e^+e^-$ annihilation at CLIC

abstract: This paper describes the effect of PYTHIA8 tunes on event simulation of $e^+e^-$ collisions with center-of-mass (CM) energies of 380 GeV and 3 TeV at the proposed CLIC collider. Event shapes, such as thrust, thrust major, thrust minor, oblateness, as well as particle multiplicities have been analyzed and relative differences with respect to the default PYTHIA8 tune were determined. The effect of tunes on top-mass reconstruction in the resolved and boosted regimes was analyzed. No statistically significant variation for reconstructed top masses using invariant masses of three jets was found for events with a CM energy of 380 GeV. For the fully boosted top reconstruction at a CM energy of 3 TeV, a significant shift in reconstructed top mass of about 700 MeV for the "Montull" tune was observed. This shift correlates with an increase in particle multiplicity compared to all other PYTHIA8 tunes.


title: Negative Ricci curvature on some non-solvable Lie groups II

abstract: We construct many examples of Lie groups with compact Levi factor admitting a left-invariant metric with negative Ricci curvature. We start with a Lie algebra with Levi factor su(n) or so(n) acting on an abelian nilradical via the representation on the space of homogeneous polynomials. In the case of su(2) we obtain a more general construction where the nilradical can be any nilpotent Lie algebra. We also prove a general result in the case when the Levi factor is a semisimple Lie algebra of non-compact type.


title: The Axion and the Goldstone Higgs

abstract: We consider the renormalizable $SO(5)/SO(4)$ $\sigma$-model, in which the Higgs particle has a pseudo-Nambu-Goldstone boson character, and explore what the minimal field extension required to implement the Peccei-Quinn symmetry (PQ) is, within the partial compositeness scenario. It turns out that the minimal model does not require the enlargement of the exotic fermionic sector, but only the addition of a singlet scalar: it is sufficient that the exotic fermions involved in partial compositeness and the singlet scalar become charged under Peccei-Quinn transformations. We explore the phenomenological predictions for photonic signals in axion searches for all models discussed. Because of the constraints imposed on the exotic fermion sector by the Standard Model fermion masses, the expected range of allowed axion-photon couplings turns out to be generically narrowed with respect to that of standard invisible axion models, impacting the experimental quest.


title: A Statistical Characterization of Localization Performance in Wireless Networks

abstract: Localization performance in wireless networks has been traditionally benchmarked using the Cramer-Rao lower bound (CRLB), given a fixed geometry of anchor nodes and a target. However, by endowing the target and anchor locations with distributions, this paper recasts this traditional, scalar benchmark as a random variable. The goal of this work is to derive an analytical expression for the distribution of this now random CRLB, in the context of Time-of-Arrival-based positioning. To derive this distribution, this work first analyzes how the CRLB is affected by the order statistics of the angles between consecutive participating anchors (i.e., internodal angles). This analysis reveals an intimate connection between the second largest internodal angle and the CRLB, which leads to an accurate approximation of the CRLB. Using this approximation, a closed-form expression for the distribution of the CRLB, conditioned on the number of participating anchors, is obtained. Next, this conditioning is eliminated to derive an analytical expression for the marginal CRLB distribution. Since this marginal distribution accounts for all target and anchor positions, across all numbers of participating anchors, it therefore statistically characterizes localization error throughout an entire wireless network. This paper concludes with a comprehensive analysis of this new network-wide-CRLB paradigm.


title: Magnetic field amplification in supernova remnants

abstract: Based on the new findings on the turbulent dynamo in \citet{XL16}, we examine the magnetic field amplification in the context of supernova remnants. Due to the strong ion-neutral collisional damping in the weakly-ionized interstellar medium, the dynamo in the preshock turbulence remains in the damping kinematic regime, which leads to {\it a linear-in-time growth of the magnetic field strength}. The resultant magnetic field structure enables effective diffusion upstream and shock acceleration of cosmic rays to energies above the "knee". Differently, the nonlinear dynamo in the postshock turbulence leads to {\it a linear-in-time growth of the magnetic energy} due to the turbulent magnetic diffusion. Given a weak initial field strength in the postshock region, the magnetic field saturates at a significant distance from the shock front as a result of the inefficiency of the nonlinear dynamo. This result is in a good agreement with existing numerical simulations and well explains the X-ray spots detected far behind the shock front.


title: Transverse instabilities of stripe domains in magnetic thin films with perpendicular magnetic anisotropy

abstract: Stripe domains are narrow, elongated, reversed regions that exist in magnetic materials with perpendicular magnetic anisotropy. Stripe domains appear as a pair of domain walls that can exhibit topology with a nonzero chirality. Recent experimental and numerical investigations identify an instability of stripe domains in the long direction as a means of nucleating isolated magnetic skyrmions. Here, the onset and nonlinear evolution of transverse instabilities for a dynamic stripe domain known as the bion stripe are investigated. Both non-topological and topological variants of the bion stripe are shown to exhibit a long-wavelength transverse instability with different characteristic features. In the former, small transverse variations in the stripe's width lead to a neck instability that eventually pinches the non-topological stripe into a chain of two-dimensional breathers composed of droplet soliton pairs. In the latter case, small variations in the stripe's center results in a snake instability whose topological structure leads to the nucleation of dynamic magnetic skyrmions and antiskyrmions as well as perimeter-modulated droplets. Quantitative, analytical predictions for both the early, linear evolution and the long-time, nonlinear evolution are achieved using an averaged Lagrangian approach that incorporates both exchange (dispersion) and anisotropy (nonlinearity). The method of analysis is general and can be applied to other filamentary structures.


title: Uniformly bounded regret in the multi-secretary problem

abstract: In the secretary problem of Cayley (1875) and Moser (1956), $n$ non-negative, independent, random variables with common distribution are sequentially presented to a decision maker who decides when to stop and collect the most recent realization. The goal is to maximize the expected value of the collected element. In the $k$-choice variant, the decision maker is allowed to make $k \leq n$ selections to maximize the expected total value of the selected elements. Assuming that the values are drawn from a known distribution with finite support, we prove that the best regret---the expected gap between the optimal online policy and its offline counterpart in which all $n$ values are made visible at time $0$---is uniformly bounded in the the number of candidates $n$ and the budget $k$. Our proof is constructive: we develop an adaptive Budget-Ratio policy that achieves this performance. The policy selects or skips values depending on where the ratio of the residual budget to the remaining time stands relative to multiple thresholds that correspond to middle points of the distribution. We also prove that being adaptive is crucial: in general, the minimal regret among non-adaptive policies grows like the square root of $n$. The difference is the value of adaptiveness.


title: The rotation problem

abstract: Any reasonable form of quantum gravity can explain (by phase interference) why on a large scale, inertial frames seem not to rotate relative to the average matter distribution in the universe without the need for absolute space, finely tuned initial conditions, or without giving up independent degrees of freedom for the gravitational field. A simple saddlepoint approximation to a path-integral calculation for a perfect fluid cosmology shows that only cosmologies with an average present relative rotation rate smaller than about $T^*H^2 \approx 10^{-71}$ radians per year could contribute significantly to a measurement of relative rotation rate in our universe, where $T^*\approx 10^{-51}$ years is the Planck time and $H \approx 10^{-10}$ yr$^{-1}$ is the present value of the Hubble parameter. A more detailed calculation (taking into account that with vorticity, flow lines are not normal to surfaces of constant global time, and approximating the action to second order in the mean square vorticity) shows that the saddlepoint at zero vorticity is isolated and that only cosmologies with an average present relative rotation rate smaller than about $T^*H^2 a_1^{1/2} \approx 10^{-73}$ radians per year could contribute significantly to a measurement of relative rotation rate in our universe, where $a_1 \approx 10^{-4}$ is the value of the cosmological scale factor at the time when matter became more significant than radiation in the cosmological expansion. This is consistent with measurements indicating a present relative rotation rate less than about $10^{-20}$ radians per year. The observed lack of relative rotation may be evidence for the existence of quantum gravity.


title: Data Mining for better material synthesis: the case of pulsed laser deposition of complex oxides

abstract: The pursuit of more advanced electronics, finding solutions to energy needs, and tackling a wealth of social issues often hinges upon the discovery and optimization of new functional materials that enable disruptive technologies or applications. However, the discovery rate of these materials is alarmingly low. Much of the information that could drive this rate higher is scattered across tens of thousands of papers in the extant literature published over several decades, and almost all of it is not collated and thus cannot be used in its entirety. Many of these limitations can be circumvented if the experimentalist has access to systematized collections of prior experimental procedures and results that can be analyzed and built upon. Here, we investigate the property-processing relationship during growth of oxide films by pulsed laser deposition. To do so, we develop an enabling software tool to (1) mine the literature of relevant papers for synthesis parameters and functional properties of previously studied materials, (2) enhance the accuracy of this mining through crowd sourcing approaches, (3) create a searchable repository that will be a community-wide resource enabling material scientists to leverage this information, and (4) provide through the Jupyter notebook platform, simple machine-learning-based analysis to learn the complex interactions between growth parameters and functional properties (all data and codes available on https://github.com/ORNL-DataMatls). The results allow visualization of growth windows, trends and outliers, and which can serve as a template for analyzing the distribution of growth conditions, provide starting points for related compounds and act as feedback for first-principles calculations. Such tools will comprise an integral part of the materials design schema in the coming decade.


title: Colorful coverings of polytopes and piercing numbers of colorful d-intervals

abstract: We prove a common strengthening of B\'ar\'any's colorful Carath\'eodory theorem and the KKMS theorem. In fact, our main result is a colorful polytopal KKMS theorem, which extends a colorful KKMS theorem due to Shih and Lee [Math. Ann. 296 (1993), no. 1, 35--61] as well as a polytopal KKMS theorem due to Komiya [Econ. Theory 4 (1994), no. 3, 463--466]. The (seemingly unrelated) colorful Carath\'eodory theorem is a special case as well. We apply our theorem to establish an upper bound on the piercing number of colorful d-interval hypergraphs, extending earlier results of Tardos [Combinatorica 15 (1995), no. 1, 123--134] and Kaiser [Discrete Comput. Geom. 18 (1997), no. 2, 195--203].


title: Generalized linear mixing model accounting for endmember variability

abstract: Endmember variability is an important factor for accurately unveiling vital information relating the pure materials and their distribution in hyperspectral images. Recently, the extended linear mixing model (ELMM) has been proposed as a modification of the linear mixing model (LMM) to consider endmember variability effects resulting mainly from illumination changes. In this paper, we further generalize the ELMM leading to a new model (GLMM) to account for more complex spectral distortions where different wavelength intervals can be affected unevenly. We also extend the existing methodology to jointly estimate the variability and the abundances for the GLMM. Simulations with real and synthetic data show that the unmixing process can benefit from the extra flexibility introduced by the GLMM.


title: The Weinstein conjecture for iterated planar contact structures

abstract: In this paper, we introduce the notions of an iterated planar Lefschetz fibration and an iterated planar open book decomposition and prove the Weinstein conjecture for contact manifolds supporting an open book that has iterated planar pages. For $n\geq 1$, we show that a $(2n+1)$-dimensional contact manifold $M$ supporting an iterated planar open book decomposition satisfies the Weinstein conjecture.


title: Optimal Placement and Patrolling of Autonomous Vehicles in Visibility-Based Robot Networks

abstract: In communication-denied or contested environments, Line-of-Sight (LoS) communication (e.g free space optical communication using infrared or visible light) becomes one of the most reliable and efficient ways to send information between geographically scattered mobile units. In this paper, we consider the problem of planning optimal locations and trajectories for a group of autonomous vehicles to see a set of units that are dispersed in an environment with obstacles. The contributions of the paper are the following: 1) We propose centralized and distributed algorithms to verify that the vehicles and units form a connected network through LoS; 2) We present an algorithm that can maintain visibility-based connectivity, if possible, by relocating a single vehicle; and 3) We study the computational


title: Continuous-wave room-temperature diamond maser

abstract: The maser, older sibling of the laser, has been confined to relative obscurity due to its reliance on cryogenic refrigeration and high-vacuum systems. Despite this it has found application in deep-space communications and radio astronomy due to its unparalleled performance as a low-noise amplifier and oscillator. The recent demonstration of a room-temperature solid- state maser exploiting photo-excited triplet states in organic pentacene molecules paves the way for a new class of maser that could find applications in medicine, security and sensing, taking advantage of its sensitivity and low noise. However, to date, only pulsed operation has been observed in this system. Furthermore, organic maser molecules have poor thermal and mechanical properties, and their triplet sub-level decay rates make continuous emission challenging: alternative materials are therefore required. Therefore, inorganic materials containing spin-defects such as diamond and silicon carbide have been proposed. Here we report a continuous-wave (CW) room-temperature maser oscillator using optically pumped charged nitrogen-vacancy (NV) defect centres in diamond. This demonstration unlocks the potential of room-temperature solid-state masers for use in a new generation of microwave devices.


title: Camera Based Two Factor Authentication Through Mobile and Wearable Devices

abstract: We introduce Pixie, a novel, camera based two factor authentication solution for mobile and wearable devices. A quick and familiar user action of snapping a photo is sufficient for Pixie to simultaneously perform a graphical password authentication and a physical token based authentication, yet it does not require any expensive, uncommon hardware. Pixie establishes trust based on both the knowledge and possession of an arbitrary physical object readily accessible to the user, called trinket. Users choose their trinkets similar to setting a password, and authenticate by presenting the same trinket to the camera. The fact that the object is the trinket, is secret to the user. Pixie extracts robust, novel features from trinket images, and leverages a supervised learning classifier to effectively address inconsistencies between images of the same trinket captured in different circumstances. Pixie achieved a false accept rate below 0.09% in a brute force attack with 14.3 million authentication attempts, generated with 40,000 trinket images that we captured and collected from public datasets. We identify master images, that match multiple trinkets, and study techniques to reduce their impact. In a user study with 42 participants over 8 days in 3 sessions we found that Pixie outperforms text based passwords on memorability, speed, and user preference. Furthermore, Pixie was easily discoverable by new users and accurate under field use. Users were able to remember their trinkets 2 and 7 days after registering them, without any practice between the 3 test dates.


title: A Computational Framework for Multi-Modal Social Action Identification

abstract: We create a computational framework for understanding social action and demonstrate how this framework can be used to build an open-source event detection tool with scalable statistical machine learning algorithms and a subsampled database of over 600 million geo-tagged Tweets from around the world. These Tweets were collected between April 1st, 2014 and April 30th, 2015, most notably when the Black Lives Matter movement began. We demonstrate how these methods can be used diagnostically-by researchers, government officials and the public-to understand peaceful and violent collective action at very fine-grained levels of time and geography.


title: Is space a word, too?

abstract: For words, rank-frequency distributions have long been heralded for adherence to a potentially-universal phenomenon known as Zipf's law. The hypothetical form of this empirical phenomenon was refined by Ben\^{i}ot Mandelbrot to that which is presently referred to as the Zipf-Mandelbrot law. Parallel to this, Herbet Simon proposed a selection model potentially explaining Zipf's law. However, a significant dispute between Simon and Mandelbrot, notable empirical exceptions, and the lack of a strong empirical connection between Simon's model and the Zipf-Mandelbrot law have left the questions of universality and mechanistic generation open. We offer a resolution to these issues by exhibiting how the dark matter of word segmentation, i.e., space, punctuation, etc., connect the Zipf-Mandelbrot law to Simon's mechanistic process. This explains Mandelbrot's refinement as no more than a fudge factor, accommodating the effects of the exclusion of the rank-frequency dark matter. Thus, integrating these non-word objects resolves a more-generalized rank-frequency law. Since this relies upon the integration of space, etc., we find support for the hypothesis that $all$ are generated by common processes, indicating from a physical perspective that space is a word, too.


title: Comment on "The Rotation-Vibration Spectrum of Diatomic Molecules with the Tietz-Hua Rotating Oscillator"

abstract: We present arguments demonstrating that the application of the Nikiforov-Uvarov polynomial method to solve the Schr\"odinger equation with the Tietz-Hua potential is valid only when $e^{-b_{h}r_{e}}\leq c_{h}<1$ and $r_{0}<r<+\infty $. In particular, it is pointed out that the numerical results with $c_{h}\neq 0$ for the diatomic molecules $\mathrm{HF}$, $\mathrm{N}_{2}$, $\mathrm{I}_{2}$, $\mathrm{H}_{2}$, $\mathrm{O}_{2}$ and $\mathrm{O}_{2}^{+}$ given in Tables 3-5 by Hamzavi and co-workers are wrong. When $-1\leq c_{h}<0$ or $0<c_{h}<e^{-b_{h}r_{e}}$, this approach is not suitable. In both cases, it is shown that the solutions of the Schr\"odinger equation are expressed in terms of the generalized hypergeometric functions $_{2}F_{1}\left( a,b,c;z\right) $. The determination of the energy levels requires the solution of transcendantal equations involving the hypergeometric function by means of the numerical procedure.


title: On the Annihilator Ideal of an Inverse Form

abstract: Let $K$ be a field. We simplify and extend work of Althaler \& D\"ur on finite sequences over $K$ by regarding $K[x^{-1},z^{-1}]$ as a $K[x,z]$ module, and studying forms in $K[x^{-1},z^{-1}]$ from first principles. Then we apply our results to finite sequences. First we define the annihilator ideal $I_F$ of a non-zero form $F\in K[x^{-1},z^{-1}]$, a homogeneous ideal. We inductively construct an ordered pair ($f_1$\,,\,$f_2$) of forms which generate $I_F$\,; our generators are special in that $z$ does not divide the leading grlex monomial of $f_1$ but $z$ divides $f_2$\,, and the sum of their total degrees is always $2-|F|$, where $|F|$ is the total degree of $F$. We show that $f_1,f_2$ is a maximal regular sequence for $I_F$, so that the height of $I_F$ is 2. The corresponding algorithm is $\sim |F|^2/2$. The row vector obtained by accumulating intermediate forms of the construction gives a minimal grlex Gr\"obner basis for $I_F$ for no extra computational cost other than storage and apply this to determining $\dim_K (K[x,z] /I_F)$\,. We show that either the form vector is reduced or a monomial of $f_1$ can be reduced by $f_2$\,. This enables us to efficiently construct the unique reduced Gr\"obner basis for $I_F$ from the vector extension of our algorithm. Then we specialise to the inverse form of a finite sequence, obtaining generator forms for its annihilator ideal and a corresponding algorithm which does not use the last 'length change' of Massey. We compute the intersection of two annihilator ideals using syzygies in $K[x,z]^5$. This improves a result of Althaler \& D\"ur. Finally, dehomogenisation induces a one-to-one correspondence ($f_1$\,,$f_2$) $\mapsto$ (minimal polynomial, auxiliary polynomial), the output of the author's variant of the Berlekamp-Massey algorithm. So we can also solve the LFSR synthesis problem via the corresponding algorithm for sequences.


title: A Tight Excess Risk Bound via a Unified PAC-Bayesian-Rademacher-Shtarkov-MDL Complexity

abstract: We present a novel notion of complexity that interpolates between and generalizes some classic existing complexity notions in learning theory: for estimators like empirical risk minimization (ERM) with arbitrary bounded losses, it is upper bounded in terms of data-independent Rademacher complexity; for generalized Bayesian estimators, it is upper bounded by the data-dependent information complexity (also known as stochastic or PAC-Bayesian, $\mathrm{KL}(\text{posterior} \operatorname{\|} \text{prior})$ complexity. For (penalized) ERM, the new complexity reduces to (generalized) normalized maximum likelihood (NML) complexity, i.e. a minimax log-loss individual-sequence regret. Our first main result bounds excess risk in terms of the new complexity. Our second main result links the new complexity via Rademacher complexity to $L_2(P)$ entropy, thereby generalizing earlier results of Opper, Haussler, Lugosi, and Cesa-Bianchi who did the log-loss case with $L_\infty$. Together, these results recover optimal bounds for VC- and large (polynomial entropy) classes, replacing localized Rademacher complexity by a simpler analysis which almost completely separates the two aspects that determine the achievable rates: 'easiness' (Bernstein) conditions and model complexity.


title: Hydrodynamic synchronization of externally driven colloids

abstract: The collective dynamics of externally driven $N_{\rm p}$-colloidal systems ($1\le N_{\rm p}\le 4$) in a confined viscous fluid have been investigated using three-dimensional direct numerical simulations with fully resolved hydrodynamics. The dynamical modes of collective particle motion are studied by changing the particle Reynolds number {as determined by the strength of the external driving force} and the confining wall distance. For a system with $N_{\rm p}= 3$, we found that at a critical Reynolds number, a dynamical mode transition occurs from the doublet-singlet mode to the triplet mode, which has not been reported experimentally. The dynamical mode transition was analyzed in detail from the following two viewpoints: (1) spectrum analysis of the time evolution of a tagged particle velocity and (2) the relative acceleration of the doublet cluster with respect to the singlet particle. For a system with $N_{\rm p}=4$, we found similar dynamical mode transitions from the doublet-singlet-singlet mode to the triplet-singlet mode and further to the quartet mode.


title: Optimally convergent hybridizable discontinuous Galerkin method for fifth-order Korteweg-de Vries type equations

abstract: We develop and analyze the first hybridizable discontinuous Galerkin (HDG) method for solving fifth-order Korteweg-de Vries (KdV) type equations. We show that the semi-discrete scheme is stable with proper choices of the stabilization functions in the numerical traces. For the linearized fifth-order equations, we prove that the approximations to the exact solution and its four spatial derivatives as well as its time derivative all have optimal convergence rates. The numerical experiments, demonstrating optimal convergence rates for both the linear and nonlinear equations, validate our theoretical findings.


title: ADA: A Game-Theoretic Perspective on Data Augmentation for Object Detection

abstract: The use of random perturbations of ground truth data, such as random translation or scaling of bounding boxes, is a common heuristic used for data augmentation that has been shown to prevent overfitting and improve generalization. Since the design of data augmentation is largely guided by reported best practices, it is difficult to understand if those design choices are optimal. To provide a more principled perspective, we develop a game-theoretic interpretation of data augmentation in the context of object detection. We aim to find an optimal adversarial perturbations of the ground truth data (i.e., the worst case perturbations) that forces the object bounding box predictor to learn from the hardest distribution of perturbed examples for better test-time performance. We establish that the game theoretic solution, the Nash equilibrium, provides both an optimal predictor and optimal data augmentation distribution. We show that our adversarial method of training a predictor can significantly improve test time performance for the task of object detection. On the ImageNet object detection task, our adversarial approach improves performance by over 16\% compared to the best performing data augmentation method


title: BigSparse: High-performance external graph analytics

abstract: We present BigSparse, a fully external graph analytics system that picks up where semi-external systems like FlashGraph and X-Stream, which only store vertex data in memory, left off. BigSparse stores both edge and vertex data in an array of SSDs and avoids random updates to the vertex data, by first logging the vertex updates and then sorting the log to sequentialize accesses to the SSDs. This newly introduced sorting overhead is reduced significantly by interleaving sorting with vertex reduction operations. In our experiments on a server with 32GB to 64GB of DRAM, BigSparse outperforms other in-memory and semi-external graph analytics systems for algorithms such as PageRank, BreadthFirst Search, and Betweenness-Centrality for terabyte-size graphs with billions of vertices. BigSparse is capable of highspeed analytics of much larger graphs, on the same machine configuration.


title: Dynamic mode decomposition for compressive system identification

abstract: Dynamic mode decomposition has emerged as a leading technique to identify spatiotemporal coherent structures from high-dimensional data, benefiting from a strong connection to nonlinear dynamical systems via the Koopman operator. In this work, we integrate and unify two recent innovations that extend DMD to systems with actuation [Proctor et al., 2016] and systems with heavily subsampled measurements [Brunton et al., 2015]. When combined, these methods yield a novel framework for compressive system identification [code is publicly available at: https://github.com/zhbai/cDMDc]. It is possible to identify a low-order model from limited input-output data and reconstruct the associated full-state dynamic modes with compressed sensing, adding interpretability to the state of the reduced-order model. Moreover, when full-state data is available, it is possible to dramatically accelerate downstream computations by first compressing the data. We demonstrate this unified framework on two model systems, investigating the effects of sensor noise, different types of measurements (e.g., point sensors, Gaussian random projections, etc.), compression ratios, and different choices of actuation (e.g., localized, broadband, etc.). In the first example, we explore this architecture on a test system with known low-rank dynamics and an artificially inflated state dimension. The second example consists of a real-world engineering application given by the fluid flow past a pitching airfoil at low Reynolds number. This example provides a challenging and realistic test-case for the proposed method, and results demonstrate that the dominant coherent structures are well characterized despite actuation and heavily subsampled data.

categories: cs.SY math.OC physics.data-an physics.flu-dyn


title: The Design and Implementation of Modern Online Programming Competitions

abstract: This paper presents a framework for the implementation of online programming competitions, including a set of principles for the design of the multiplayer game and a practical framework for the construction of the competition environment. The paper presents a successful example competition, the 2016-17 Halite challenge, and briefly mentions a second competition, the Halite II challenge, which launched in October 2017.


title: Learning Discrete Weights Using the Local Reparameterization Trick

abstract: Recent breakthroughs in computer vision make use of large deep neural networks, utilizing the substantial speedup offered by GPUs. For applications running on limited hardware, however, high precision real-time processing can still be a challenge. One approach to solving this problem is training networks with binary or ternary weights, thus removing the need to calculate multiplications and significantly reducing memory size. In this work, we introduce LR-nets (Local reparameterization networks), a new method for training neural networks with discrete weights using stochastic parameters. We show how a simple modification to the local reparameterization trick, previously used to train Gaussian distributed weights, enables the training of discrete weights. Using the proposed training we test both binary and ternary models on MNIST, CIFAR-10 and ImageNet benchmarks and reach state-of-the-art results on most experiments.


title: Program Synthesis using Abstraction Refinement

abstract: We present a new approach to example-guided program synthesis based on counterexample-guided abstraction refinement. Our method uses the abstract semantics of the underlying DSL to find a program $P$ whose abstract behavior satisfies the examples. However, since program $P$ may be spurious with respect to the concrete semantics, our approach iteratively refines the abstraction until we either find a program that satisfies the examples or prove that no such DSL program exists. Because many programs have the same input-output behavior in terms of their abstract semantics, this synthesis methodology significantly reduces the search space compared to existing techniques that use purely concrete semantics. While synthesis using abstraction refinement (SYNGAR) could be implemented in different settings, we propose a refinement-based synthesis algorithm that uses abstract finite tree automata (AFTA). Our technique uses a coarse initial program abstraction to construct an initial AFTA, which is iteratively refined by constructing a proof of incorrectness of any spurious program. In addition to ruling out the spurious program accepted by the previous AFTA, proofs of incorrectness are also useful for ruling out many other spurious programs. We implement these ideas in a framework called \tool. We have used the BLAZE framework to build synthesizers for string and matrix transformations, and we compare BLAZE with existing techniques. Our results for the string domain show that BLAZE compares favorably with FlashFill, a domain-specific synthesizer that is now deployed in Microsoft PowerShell. In the context of matrix manipulations, we compare BLAZE against Prose, a state-of-the-art general-purpose VSA-based synthesizer, and show that BLAZE results in a 90x speed-up over Prose.


title: On the Computational Complexity of the Bipartizing Matching Problem

abstract: We study the problem of determining whether a given graph~$G=(V,E)$ admits a matching~$M$ whose removal destroys all odd cycles of~$G$ (or equivalently whether~$G-M$ is bipartite). This problem is equivalent to determine whether~$G$ admits a~$(2,1)$-coloring, which is a~$2$-coloring of~$V(G)$ such that each color class induces a graph of maximum degree at most~$1$. We determine a dichotomy related to the~{\sf NP}-completeness of this problem, where we show that it is~{\sf NP}-complete even for $3$-colorable planar graphs of maximum degree~$4$, while it is known that the problem can be solved in polynomial time for graphs of maximum degree at most~$3$. In addition we present polynomial-time algorithms for some graph classes, including graphs in which every odd cycle is a triangle, graphs of small dominating sets, and~$P_5$-free graphs. Additionally, we show that the problem is fixed parameter tractable when parameterized by the clique-width, which implies polynomial-time solution for many interesting graph classes, such as distance-hereditary, outerplanar, and chordal graphs. Finally, an~$O\left(2^{O\left(vc(G)\right)} \cdot n\right)$-time algorithm and a kernel of at most~$2\cdot nd(G)$ vertices are presented, where~$vc(G)$ and~$nd(G)$ are the vertex cover number and the neighborhood diversity of~$G$, respectively.


title: Towards Black-box Iterative Machine Teaching

abstract: In this paper, we make an important step towards the black-box machine teaching by considering the cross-space machine teaching, where the teacher and the learner use different feature representations and the teacher can not fully observe the learner's model. In such scenario, we study how the teacher is still able to teach the learner to achieve faster convergence rate than the traditional passive learning. We propose an active teacher model that can actively query the learner (i.e., make the learner take exams) for estimating the learner's status and provably guide the learner to achieve faster convergence. The sample complexities for both teaching and query are provided. In the experiments, we compare the proposed active teacher with the omniscient teacher and verify the effectiveness of the active teacher model.


title: Ionosphere-thermosphere global time response to geomagnetic storms

abstract: In this study, we investigate thermospheric neutral mass density heating associated with 168 CME-driven geomagnetic storms in the period of May 2001 to September 2011. We use neutral density measured by two low-Earth orbit satellites: CHAMP and GRACE. For each storm, we superpose geomagnetic and density data for the time when the IMF B$_\mathrm{z}$ component turns sharply southward chosen as the zero epoch time. This indicates the storm main phase onset. We find that the average SYM-H index reaches the minimum of $-$42 nT near 12 hours after storm main phase onset. The Joule heating is enhanced by approximately 200\% in comparison to quiet values. In respect to thermosphere density, on average, high latitude regions (auroral zones and polar caps) of both hemispheres are highly heated in the first 1.5 hour of the storm. The equatorial response is presumably associated with direct equator-ward propagation of TADs (traveling atmospheric disturbances). A slight north-south asymmetry in thermosphere heating is found and is most likely due to a positive B$_\mathrm{y}$ component in the first hours of the storm main phase.


title: Symmetries of the hydrogen atom and algebraic families

abstract: We show how the Schr\"{o}dinger equation for the hydrogen atom in two dimensions gives rise to an algebraic family of Harish-Chandra pairs that codifies hidden symmetries. The hidden symmetries vary continuously between $SO(3)$, $SO(2,1)$ and the Euclidean group $O(2)\ltimes \mathbb{R}^2$. We show that solutions of the Schr\"{o}dinger equation may be organized into an algebraic family of Harish-Chandra modules. Furthermore, we use Jantzen filtration techniques to algebraically recover the spectrum of the Schr\"{o}dinger operator. This is a first application to physics of the algebraic families of Harish-Chandra pairs and modules developed in the work of Bernstein et al. [Int. Math. Res. Notices, rny147 (2018); rny146 (2018)].


title: High Performance Canny Edge Detector using Parallel Patterns for Scalability on Modern Multicore Processors

abstract: Canny Edge Detector (CED) is an edge detection operator commonly used by most Image Feature Extraction (IFE) Algorithms and Image Processing Applications. This operator involves the use of a multi-stage algorithm to detect edges in a wide range of images. Edge detection is at the forefront of image processing and hence, it is crucial to have at an up to scale level. Multicore Processors have emerged as the next solution for tackling compute intensive tasks that have a high demand for computational power. Having significant changes that restructured the microprocessor industry, it is evident that the best way to promote efficiency and improve performance is no longer by increasing the clock speeds on traditional monolithic processors but by adopting and utilizing Processors with Multicore architectures. In this paper we provide a high performance implementation of Canny Edge Detector using parallel patterns for improved performance and Scalability on Multicore Processors. The results show significant improvements in overall performance and this proves that our implementation using parallel patterns does not under utilize resources but scales well for multicore processors.


title: Stochastic Backward Euler: An Implicit Gradient Descent Algorithm for $k$-means Clustering

abstract: In this paper, we propose an implicit gradient descent algorithm for the classic $k$-means problem. The implicit gradient step or backward Euler is solved via stochastic fixed-point iteration, in which we randomly sample a mini-batch gradient in every iteration. It is the average of the fixed-point trajectory that is carried over to the next gradient step. We draw connections between the proposed stochastic backward Euler and the recent entropy stochastic gradient descent (Entropy-SGD) for improving the training of deep neural networks. Numerical experiments on various synthetic and real datasets show that the proposed algorithm provides better clustering results compared to $k$-means algorithms in the sense that it decreased the objective function (the cluster) and is much more robust to initialization.


title: Localization for MCMC: sampling high-dimensional posterior distributions with local structure

abstract: We investigate how ideas from covariance localization in numerical weather prediction can be used in Markov chain Monte Carlo (MCMC) sampling of high-dimensional posterior distributions arising in Bayesian inverse problems. To localize an inverse problem is to enforce an anticipated "local" structure by (i) neglecting small off-diagonal elements of the prior precision and covariance matrices; and (ii) restricting the influence of observations to their neighborhood. For linear problems we can specify the conditions under which posterior moments of the localized problem are close to those of the original problem. We explain physical interpretations of our assumptions about local structure and discuss the notion of high dimensionality in local problems, which is different from the usual notion of high dimensionality in function space MCMC. The Gibbs sampler is a natural choice of MCMC algorithm for localized inverse problems and we demonstrate that its convergence rate is independent of dimension for localized linear problems. Nonlinear problems can also be tackled efficiently by localization and, as a simple illustration of these ideas, we present a localized Metropolis-within-Gibbs sampler. Several linear and nonlinear numerical examples illustrate localization in the context of MCMC samplers for inverse problems.