title: Effective drag of a rod in fluid-saturated granular beds

abstract: We measure the drag encountered by a vertically oriented rod moving across a sedimented granular bed immersed in a fluid under steady-state conditions. At low rod speeds, the presence of the fluid leads to a lower drag because of buoyancy, whereas a significantly higher drag is observed with increasing speeds. The drag as a function of depth is observed to decrease from being quadratic at low speeds to appearing more linear at higher speeds. By scaling the drag with the average weight of the grains acting on the rod, we obtain the effective friction $\mu_e$ encountered over six orders of magnitude of speeds. While a constant $\mu_e$ is found when the grain size, rod depth and fluid viscosity are varied at low speeds, a systematic increase is observed as the speed is increased. We analyze $\mu_e$ in terms of the inertial number $I$ and viscous number $J$ to understand the relative importance of inertia and viscous forces, respectively. For sufficiently large fluid viscosities, we find that the effect of varying the speed, depth, and viscosity can be described by the empirical function $\mu_e = \mu_o + k J^n$, where $\mu_o$ is the effective friction measured in the quasi-static limit, and $k$ and $n$ are material constants. The drag is then analyzed in terms of the effective viscosity $\eta_e$ and found to decrease systematically as a function of $J$. We further show that $\eta_e$ as a function of $J$ is directly proportional to the fluid viscosity and the $\mu_e$ encountered by the rod.


title: Tight Polynomial Worst-Case Bounds for Loop Programs

abstract: In 2008, Ben-Amram, Jones and Kristiansen showed that for a simple programming language - representing non-deterministic imperative programs with bounded loops, and arithmetics limited to addition and multiplication - it is possible to decide precisely whether a program has certain growth-rate properties, in particular whether a computed value, or the program's running time, has a polynomial growth rate. A natural and intriguing problem was to move from answering the decision problem to giving a quantitative result, namely, a tight polynomial upper bound. This paper shows how to obtain asymptotically-tight, multivariate, disjunctive polynomial bounds for this class of programs. This is a complete solution: whenever a polynomial bound exists it will be found. A pleasant surprise is that the algorithm is quite simple; but it relies on some subtle reasoning. An important ingredient in the proof is the forest factorization theorem, a strong structural result on homomorphisms into a finite monoid.


title: SurReal: Fr\'echet Mean and Distance Transform for Complex-Valued Deep Learning

abstract: We develop a novel deep learning architecture for naturally complex-valued data, which is often subject to complex scaling ambiguity. We treat each sample as a field in the space of complex numbers. With the polar form of a complex-valued number, the general group that acts in this space is the product of planar rotation and non-zero scaling. This perspective allows us to develop not only a novel convolution operator using weighted Fr\'echet mean (wFM) on a Riemannian manifold, but also a novel fully connected layer operator using the distance to the wFM, with natural equivariant properties to non-zero scaling and planar rotation for the former and invariance properties for the latter. Compared to the baseline approach of learning real-valued neural network models on the two-channel real-valued representation of complex-valued data, our method achieves surreal performance on two publicly available complex-valued datasets: MSTAR on SAR images and RadioML on radio frequency signals. On MSTAR, at 8% of the baseline model size and with fewer than 45,000 parameters, our model improves the target classification accuracy from 94% to 98% on this highly imbalanced dataset. On RadioML, our model achieves comparable RF modulation classification accuracy at 10% of the baseline model size.


title: Prompt ${J/\psi}$-pair production at the LHC: impact of loop-induced contributions and of the colour-octet mechanism

abstract: Prompt double-$J/\psi$ production at high-energy hadron colliders can be considered as a golden channel to probe double parton scatterings (DPS) --in particular to study gluon-gluon correlations inside the proton-- and, at the same time, to measure the distribution of linearly-polarised gluons inside the proton. Such studies however require a good control of both single and DPS in the respective regions where they are carried out. In this context, we have critically examined two mechanisms of single parton scatterings (SPS) that may be kinematically enhanced where DPS are thought to be dominant, even though they are either at higher orders in the strong-coupling or velocity expansion. First, we have considered a gauge-invariant and infrared-safe subset of the loop-induced contribution via Colour-Singlet (CS) transitions. We have found it to become the leading CS SPS contributions at large rapidity separation, yet too small to account for the data without invoking the presence of DPS yields. Second, we have surveyed the possible Colour-Octet (CO) contributions using both old and up-to-date non-perturbative long distance matrix elements (LDMEs). We have found that the pure CO yields crucially depend on the LDMEs. Among all the LDMEs we used, only two result into a visible modification of the NRQCD (CS+CO) yield, but only in two kinematical distributions measured by ATLAS, those of the rapidity separation and of the pair invariant mass. These modifications however do not impact the control region used for their DPS study.


title: On DICE-free Smart Cities, Particulate Matter, and Feedback-Enabled Access Control

abstract: The link between transport related emissions and human health is a major issue for city municipalities worldwide. PM emissions from exhaust and non-exhaust sources are one of the main worrying contributors to air-pollution. In this paper, we challenge the notion that a ban on internal combustion engine vehicles will result in clean and safe air in our cities, since emissions from tyres and other non-exhaust sources are expected to increase in the near future. To this end, we present data from the city of Dublin that document that the current amount of tyre-related PM emissions in the city might already be above or close to the levels deemed safe by the World Health Organization. As a solution to this problem, we present a feedback-enabled distributed access control mechanism and ride-sharing scheme to limit the number of vehicles in a city and therefore maintain the amount of transport-related PM to safe levels.


title: Conditional Expectation, Entropy, and Transport for Convex Gibbs Laws in Free Probability

abstract: Let $(X_1,\dots,X_m)$ be self-adjoint non-commutative random variables distributed according to the free Gibbs law given by a sufficiently regular convex and semi-concave potential $V$, and let $(S_1,\dots,S_m)$ be a free semicircular family. We show that conditional expectations and conditional non-microstates free entropy given $X_1$, \dots, $X_k$ arise as the large $N$ limit of the corresponding conditional expectations and entropy for the random matrix models associated to $V$. Then by studying conditional transport of measure for the matrix models, we construct an isomorphism $\mathrm{W}^*(X_1,\dots,X_m) \to \mathrm{W}^*(S_1,\dots,S_m)$ which maps $\mathrm{W}^*(X_1,\dots,X_k)$ to $\mathrm{W}^*(S_1,\dots,S_k)$ for each $k = 1, \dots, m$, and which also witnesses the Talagrand inequality for the law of $(X_1,\dots,X_m)$ relative to the law of $(S_1,\dots,S_m)$.


title: The Power Spectra of Polarized, Dusty Filaments

abstract: We develop an analytic model for the power spectra of polarized filamentary structures as a way to study the Galactic polarization foreground to the Cosmic Microwave Background. Our approach is akin to the cosmological halo-model framework, and reproduces the main features of the Planck 353 GHz power spectra. We model the foreground as randomly-oriented, three-dimensional, spheroidal filaments, accounting for their projection onto the sky. The main tunable parameters are the distribution of filament sizes, the filament physical aspect ratio, and the dispersion of the filament axis around the local magnetic field direction. The abundance and properties of filaments as a function of size determine the slopes of the foreground power spectra, as we show via scaling arguments. The filament aspect ratio determines the ratio of $B$-mode power to $E$-mode power, and specifically reproduces the Planck-observed dust ratio of one-half when the short axis is roughly one-fourth the length of the long axis. Filament misalignment to the local magnetic field determines the $TE$ cross-correlation, and to reproduce Planck measurements, we need a (three-dimensional) misalignment angle with a root mean squared dispersion of about 50 degrees. These parameters are not sensitive to the particular filament density profile. By artificially skewing the distribution of the misalignment angle, this model can reproduce the Planck-observed (and parity-violating) $TB$ correlation. The skewing of the misalignment angle necessary to explain $TB$ will cause a yet-unobserved, positive $EB$ dust correlation, a possible target for future experiments.


title: Block-coordinate and incremental aggregated proximal gradient methods for nonsmooth nonconvex problems

abstract: This paper analyzes block-coordinate proximal gradient methods for minimizing the sum of a separable smooth function and a (nonseparable) nonsmooth function, both of which are allowed to be nonconvex. The main tool in our analysis is the forward-backward envelope (FBE), which serves as a particularly suitable continuous and real-valued Lyapunov function. Global and linear convergence results are established when the cost function satisfies the Kurdyka-\L ojasiewicz property without imposing convexity requirements on the smooth function. Two prominent special cases of the investigated setting are regularized finite sum minimization and the sharing problem; in particular, an immediate byproduct of our analysis leads to novel convergence results and rates for the popular Finito/MISO algorithm in the nonsmooth and nonconvex setting with very general sampling strategies. This paper analyzes block-coordinate proximal gradient methods for minimizing the sum of a separable smooth function and a (nonseparable) nonsmooth function, both of which are allowed to be nonconvex. The main tool in our analysis is the forward-backward envelope (FBE), which serves as a particularly suitable continuous and real-valued Lyapunov function. Global and linear convergence results are established when the cost function satisfies the Kurdyka-\L ojasiewicz property without imposing convexity requirements on the smooth function. Two prominent special cases of the investigated setting are regularized finite sum minimization and the sharing problem; in particular, an immediate byproduct of our analysis leads to novel convergence results and rates for the popular Finito/MISO algorithm in the nonsmooth and nonconvex setting with very general sampling strategies.


title: Birth and decay of tensional wrinkles in hyperelastic sheets

abstract: We demonstrate with experiments that wrinkling in stretched latex sheets occur over finite strains, and that their amplitudes grow and then decay to zero over a greater range of applied strains compared with linear elastic materials. The wrinkles occur provided the sheet is sufficiently thin compared to its width, and only over a finite range of length-to-width ratios. We show with simulations that the Mooney-Rivlin hyperelastic model describes the observed growth and decay of the wrinkles in our experiments. The decrease of wavelength with applied tension is found to be consistent with a far-from-threshold scenario proposed by Cerda and Mahadevan in 2003. However, the amplitude is observed to decrease with increasing tensile load, in contrast with the prediction of their original model. We address the crucial assumption of {\it collapse of compressive stress}, as opposed to collapse of compressive strain, underlying the far-from-threshold analysis, and test it by measuring the actual arc-length of the stretched sheet in the transverse direction and its difference from the width of a planar projection of the wrinkled shape. Our experiments and numerical simulations indicate a complete {\it collapse of the compressive stress}, and reveal that a proper implementation of the far-from-threshold analysis is consistent with the non-monotonic dependence of the amplitude on applied tensile load observed in experiments and simulations. Thus, our work support and extend far-from-threshold analysis to the stretching problem of rectangular hyperelastic sheets.


title: On the spatial resolution of EBSD in magnesium

abstract: We measured the physical lateral resolution of the electron backscatter diffraction (EBSD) technique for the case of pure magnesium and tungsten. Spatial resolution, among other parameters, depends significantly on the accelerating voltage and the atomic number of the material. For the case of lighter metals, it is supposed to be lower than in the case of heavier metals for a given accelerating voltage. In the present work, lateral resolution was measured in dependence of accelerating voltage on a straight high angle grain boundary which was positioned parallel (horizontal boundary) and perpendicular (vertical boundary) to the tilt axis of the specimen. For magnesium the best lateral resolution of 240 nm was obtained at an accelerating voltage of 5 kV. The resolution dramatically worsened to values as high as 3500 nm as the voltage was increased from 15 kV to 30 kV. The aspect ratio of horizontal and vertical lateral resolution tended to 1.0 at the accelerating voltage of 5 kV and to 2.5 at the accelerating voltage of 30 kV. These values as function of accelerating voltages were compared with those obtained on the high atomic number metal tungsten. Here resolution at 5 kV was about a quarter of that of magnesium. With increasing voltage, the value almost didnt change. For all voltages the resolution aspect ratio stayed close to 1.0.


title: Parametric estimation for convolutionally observed diffusion processes

abstract: We propose a new statistical observation scheme of diffusion processes named convolutional observation, where it is possible to deal with smoother observation than ordinary diffusion processes by considering convolution of diffusion processes and some kernel functions with respect to time parameter. We discuss the estimation and test theories for the parameter determining the smoothness of the observation, as well as the least-square-type estimation for the parameters in the diffusion coefficient and the drift one of the latent diffusion process. In addition to the theoretical discussion, we also examine the performance of the estimation and the test with computational simulation, and show an example of real data analysis for one EEG data whose observation can be regarded as smoother one than ordinary diffusion processes with statistical significance.


title: Efficient and Effective Context-Based Convolutional Entropy Modeling for Image Compression

abstract: Precise estimation of the probabilistic structure of natural images plays an essential role in image compression. Despite the recent remarkable success of end-to-end optimized image compression, the latent codes are usually assumed to be fully statistically factorized in order to simplify entropy modeling. However, this assumption generally does not hold true and may hinder compression performance. Here we present context-based convolutional networks (CCNs) for efficient and effective entropy modeling. In particular, a 3D zigzag scanning order and a 3D code dividing technique are introduced to define proper coding contexts for parallel entropy decoding, both of which boil down to place translation-invariant binary masks on convolution filters of CCNs. We demonstrate the promise of CCNs for entropy modeling in both lossless and lossy image compression. For the former, we directly apply a CCN to the binarized representation of an image to compute the Bernoulli distribution of each code for entropy estimation. For the latter, the categorical distribution of each code is represented by a discretized mixture of Gaussian distributions, whose parameters are estimated by three CCNs. We then jointly optimize the CCN-based entropy model along with analysis and synthesis transforms for rate-distortion performance. Experiments on the Kodak and Tecnick datasets show that our methods powered by the proposed CCNs generally achieve comparable compression performance to the state-of-the-art while being much faster.


title: Comment on "Evidence for Supersymmetry in the Random-Field Ising Model at D=5''

abstract: In a recent letter, Fytas et al. [Phys. Rev. Lett. 122, 240603 (2019)] study the critical point of the equilibrium random-field Ising model (RFIM) in $D=5$ by means of state-of-art zero-temperature lattice simulations. We show that their numerical results concerning the virtual absence of violation of supersymmetry and the characterization of the main corrections to scaling are in full agreement with our predictions by means of a nonperturbative functional renormalization group and with the existence of a critical dimension $D_{DR}\approx 5.1$ above which supersymmetry and dimensional reduction are valid at the fixed point and below which they break down.


title: Efficient interpolation and evolution of parton distribution functions

abstract: We present an efficient numerical solution of the DGLAP equations for single and double parton distribution functions (PDFs and DPDs), based on the Chebyshev interpolation of these functions. For PDF evolution, our method allows for a higher numerical accuracy using a considerably smaller number of grid points compared to other methods. The DPD evolution is realized using an affordable number of grid points, and allows for two independent renormalization scales for the two partons. Both methods include NNLO DGLAP kernels and flavor matching.


title: Harmonic Analysis on the Positive Rationals. Computation of Character Sums

abstract: The authors review results implicit in their recent paper [2] on the product/quotient representation of rationals by rationals of the type $( an + b )/ ( An+ B )$ and give a detailed account of a particular related non-intuitive simultaneous ( i.e. two dimensional ) representation of pairs of rationals.


title: Isoperimetric relations between Dirichlet and Neumann eigenvalues

abstract: Inequalities between the Dirichlet and Neumann eigenvalues of the Laplacian have received much attention in the literature, but open problems abound. Here, we study the number of Neumann eigenvalues no greater than the first Dirichlet eigenvalue. Based on a combination of analytical and numerical results, we conjecture that this number is controlled by the isoperimetric ratio of the domain. This has applications to the nodal deficiency of eigenfunctions and is closely related to a long-standing conjecture of Yau on the Hausdorff measure of nodal sets.


title: Weak ferromagnetism in hexagonal Mn3Z (Z=Sn, Ge, Ga) alloys

abstract: We present combined spin model and first principles electronic structure calculations to study the weak ferromagnetism in bulk Mn$_3$Z (Z=Sn, Ge, Ga) compounds. The spin model parameters were determined from a spin-cluster expansion technique based on the relativistic disordered local moment formalism implemented in the screened Korringa--Kohn--Rostoker method. We describe the magnetic ground state of the system within a three-sublattice model and investigate the formation of the weak ferromagnetic states in terms of the relevant model parameters. First, we give a group-theoretical argument how the point-group symmetry of the lattice leads to the formation of weak ferromagnetic states. Then we study the ground states of the classical spin model and derive analytical expressions for the weak ferromagnetic distortions by recovering the main results of the group-theoretical analysis. As a third approach we obtain the weak ferromagnetic ground states from self-consistent density functional calculations and compare our results with previous first principles calculations and with available experimental data. In particular, we demonstrate that the orbital moments follow a decomposition predicted by group theory. For a deeper understanding of the formation of weak ferromagnetism we selectively trace the effect of the spin-orbit coupling at the Mn and Z sites. In addition, for the case of Mn$_3$Ga, we gain information on the role of the induced moment of Ga from constrained local density functional calculations.


title: All-Optical Bose-Einstein Condensates in Microgravity

abstract: We report on the all-optical production of Bose-Einstein condensates in microgravity using a combination of grey molasses cooling, light-shift engineering and optical trapping in a painted potential. Forced evaporative cooling in a 3-m high Einstein elevator results in $4 \times 10^4$ condensed atoms every 13.5 s, with a temperature as low as 35 nK. In this system, the atomic cloud can expand in weightlessness for up to 400 ms, paving the way for atom interferometry experiments with extended interrogation times and studies of ultra-cold matter physics at low energies on ground or in Space.


title: Variations on the Chebyshev-Lagrange Activation Function

abstract: We seek to improve the data efficiency of neural networks and present novel implementations of parameterized piece-wise polynomial activation functions. The parameters are the y-coordinates of n+1 Chebyshev nodes per hidden unit and Lagrangian interpolation between the nodes produces the polynomial on [-1, 1]. We show results for different methods of handling inputs outside [-1, 1] on synthetic datasets, finding significant improvements in capacity of expression and accuracy of interpolation in models that compute some form of linear extrapolation from either ends. We demonstrate competitive or state-of-the-art performance on the classification of images (MNIST and CIFAR-10) and minimally-correlated vectors (DementiaBank) when we replace ReLU or tanh with linearly extrapolated Chebyshev-Lagrange activations in deep residual architectures.


title: Modeling non-thermal emission from the jet-launching region of M 87 with adaptive mesh refinement

abstract: The galaxy M 87 harbors a kiloparsec-scale relativistic jet, whose origin coincides with a supermassive black hole. Observational mm-VLBI campaigns are capable of resolving the jet-launching region at the scale of the event horizon. In order to provide a context for interpreting these observations, realistic general-relativistic magnetohydrodynamical (GRMHD) models of the accretion flow are constructed. The characteristics of the observed spectral-energy distribution (SED) depend on the shape of the electrons' energy-distribution function (eDF). The dependency on the eDF is omitted in the modeling of the first Event Horizon Telescope results. In this work, we aim to model the M 87 SED from radio up to NIR/optical frequencies using a thermal-relativistic Maxwell- J\"uttner distribution, as well as a relativistic $\kappa$-distribution function. The electrons are injected based on sub-grid, particle-in-cell parametrizations for sub-relativistic reconnection. A GRMHD simulation in Cartesian-Kerr-Schild coordinates, using eight levels of adaptive mesh refinement (AMR), forms the basis of our model. To obtain spectra and images, the GRMHD data is post-processed with the ray-tracing code RAPTOR, which is capable of ray tracing through AMR GRMHD simulation data. We obtain radio spectra in both the thermal-jet and $\kappa$-jet models consistent with radio observations. Additionally, the $\kappa$-jet models also recover the NIR/optical emission. The models recover the observed source sizes and core shifts and obtain a jet power of $\approx 10^{43}$ ergs/s. In the $\kappa$-jet models, both the accretion rates and jet powers are approximately two times lower than the thermal-jet model. The frequency cut-off observed at $\nu \approx 10^{15}$ Hz is recovered when the accelerator size is $10^6$ - $10^8$ cm, this could potentially point to an upper limit for plasmoid sizes in the jet of M 87.


title: Inclusion modulo nonstationary

abstract: A classical theorem of Hechler asserts that the structure $\left(\omega^\omega,\le^*\right)$ is universal in the sense that for any $\sigma$-directed poset P with no maximal element, there is a ccc forcing extension in which $\left(\omega^\omega,\le^*\right)$ contains a cofinal order-isomorphic copy of P. In this paper, we prove a consistency result concerning the universality of the higher analogue $\left(\kappa^\kappa,\le^S\right)$: Theorem. Assume GCH. For every regular uncountable cardinal $\kappa$, there is a cofinality-preserving GCH-preserving forcing extension in which for every analytic quasi-order Q over $\kappa^\kappa$ and every stationary subset S of $\kappa$, there is a Lipschitz map reducing Q to $(\kappa^\kappa,\le^S)$.


title: Dynamics of pattern formation and emergence of swarming in C. elegans

abstract: Many animals in their natural habitat exhibit collective motion and form complex patterns to tackle environmental difficulties. Several physical and biological factors, such as animal motility, population densities, and chemical cues, play significant roles in this process. However, very little is known about how sensory information interplays with all these factors and controls the dynamics of collective response and pattern formation. Here, we use a model organism, Caenorhabditis elegans, to study the direct relation between oxygen sensing, pattern formation, and the emergence of swarming in active worm aggregates. We find that when thousands of animals gather on food, bacteria-mediated decrease in oxygen levels slowed down the animals and triggers motility-induced phase separation. Three coupled factors bacterial accumulation, aerotaxis, and population density act together and control the dynamics of pattern formation. Through several intermediate stages, aggregates converge to a large scale swarming phase and collectively move across the bacterial lawn. Additionally, our theoretical model captures behavioral differences resulting from the genetic variations and oxygen sensitivity. Altogether, our study provides many physical insights and a new platform for investigating the complex relationship between neural sensitivity, collective dynamics, and pattern formation.


title: Is It Worth the Attention? A Comparative Evaluation of Attention Layers for Argument Unit Segmentation

abstract: Attention mechanisms have seen some success for natural language processing downstream tasks in recent years and generated new State-of-the-Art results. A thorough evaluation of the attention mechanism for the task of Argumentation Mining is missing, though. With this paper, we report a comparative evaluation of attention layers in combination with a bidirectional long short-term memory network, which is the current state-of-the-art approach to the unit segmentation task. We also compare sentence-level contextualized word embeddings to pre-generated ones. Our findings suggest that for this task the additional attention layer does not improve upon a less complex approach. In most cases, the contextualized embeddings do also not show an improvement on the baseline score.


title: Nonequilibrium and Fluctuation Relation

abstract: A review on the fluctuation relation, fluctuation theorem and related topics.


title: Full next-to-leading-order calculations of Higgs boson decay rates in models with non-minimal scalar sectors

abstract: We present a complete set of decay rates of the Higgs boson with the mass of 125 GeV at the full next-to-leading order in a variety of extended Higgs models; i.e., a model with an additional real singlet scalar field, four types of two Higgs doublet models and the inert doublet model. All the one-loop contributions due to QCD and electroweak interactions as well as scalar interactions are taken into account, and the calculations are systematically performed. Branching ratios for all the decay modes are evaluated in these models, and patterns of deviations in each decay mode from the standard model predictions are comprehensively analyzed. We show how these models with extended Higgs sectors can be distinguished by using our calculation of the branching ratios and future precision measurements of the Higgs boson decays.


title: Black-hole scattering with general spin directions from minimal-coupling amplitudes

abstract: We study the link between classical scattering of spinning black holes and quantum amplitudes for massive spin-$s$ particles. Generic spin orientations of the black holes are considered, allowing their spins to be deflected on par with their momenta. We rederive the spin-exponentiated structure of the relevant tree-level amplitude from minimal coupling to Einstein's gravity, which in the $s\to\infty$ limit generates the black holes' complete series of spin-induced multipoles. The resulting scattering function is seen to encode in a simple way the known net changes in the black-hole momenta and spins at first post-Minkowskian order. We connect our findings to a rigorous framework developed elsewhere for computing such observables from amplitudes.


title: Principal Tensor Strikes Again: Separability of Vector Equations with Torsion

abstract: Many black hole spacetimes with a 3-form field exhibit a hidden symmetry encoded in a torsion generalization of the principal Killing--Yano tensor. This tensor determines basic properties of such black holes while also underlying the separability of the Hamilton--Jacobi, Klein--Gordon, and (torsion-modified) Dirac field equations in their background. As a specific example, we consider the Chong--Cveti\v{c}--L\"u--Pope black hole of $D=5$ minimal gauged supergravity and show that the torsion-modified vector field equations can also be separated, with the principal tensor playing a key role in the separability ansatz. For comparison, separability of the Proca field in higher-dimensional Kerr--NUT--AdS spacetimes (including new explicit formulae in odd dimensions) is also presented.


title: A Logic-Based Learning Approach to Explore Diabetes Patient Behaviors

abstract: Type I Diabetes (T1D) is a chronic disease in which the body's ability to synthesize insulin is destroyed. It can be difficult for patients to manage their T1D, as they must control a variety of behavioral factors that affect glycemic control outcomes. In this paper, we explore T1D patient behaviors using a Signal Temporal Logic (STL) based learning approach. STL formulas learned from real patient data characterize behavior patterns that may result in varying glycemic control. Such logical characterizations can provide feedback to clinicians and their patients about behavioral changes that patients may implement to improve T1D control. We present both individual- and population-level behavior patterns learned from a clinical dataset of 21 T1D patients.


title: Implementation of a Hybrid Classical-Quantum Annealing Algorithm for Logistic Network Design

abstract: The logistic network design is an abstract optimization problem that, under the assumption of minimal cost, seeks the optimal configuration of the supply chain's infrastructures and facilities based on customer demand. Key economic decisions are taken about the location, number, and size of manufacturing facilities and warehouses based on the optimal solution. Therefore, improvements in the methods to address this question, which is known to be in the NP-hard complexity class, would have relevant financial consequences. Here, we implement in the D-Wave quantum annealer a hybrid classical-quantum annealing algorithm. The cost function with constraints is translated to a spin Hamiltonian, whose ground state encodes the searched result. As a benchmark, we measure the accuracy of results for a set of paradigmatic problems against the optimal published solutions (the error is on average below $1\%$), and the performance is compared against the classical algorithm, showing a remarkable reduction in the number of iterations. This work shows that state-of-the-art quantum annealers may codify and solve relevant supply-chain problems even still far from useful quantum supremacy.


title: Distribution-Independent PAC Learning of Halfspaces with Massart Noise

abstract: We study the problem of {\em distribution-independent} PAC learning of halfspaces in the presence of Massart noise. Specifically, we are given a set of labeled examples $(\mathbf{x}, y)$ drawn from a distribution $\mathcal{D}$ on $\mathbb{R}^{d+1}$ such that the marginal distribution on the unlabeled points $\mathbf{x}$ is arbitrary and the labels $y$ are generated by an unknown halfspace corrupted with Massart noise at noise rate $\eta<1/2$. The goal is to find a hypothesis $h$ that minimizes the misclassification error $\mathbf{Pr}_{(\mathbf{x}, y) \sim \mathcal{D}} \left[ h(\mathbf{x}) \neq y \right]$. We give a $\mathrm{poly}\left(d, 1/\epsilon \right)$ time algorithm for this problem with misclassification error $\eta+\epsilon$. We also provide evidence that improving on the error guarantee of our algorithm might be computationally hard. Prior to our work, no efficient weak (distribution-independent) learner was known in this model, even for the class of disjunctions. The existence of such an algorithm for halfspaces (or even disjunctions) has been posed as an open question in various works, starting with Sloan (1988), Cohen (1997), and was most recently highlighted in Avrim Blum's FOCS 2003 tutorial.


title: Lower bounds on the radius of spatial analyticity for the Kawahara equation

abstract: In this paper we obtain lower bounds on the radius of spatial analyticity of solutions to the Kawahara equation $u_t + uu_x + \alpha u_{xxx} + \beta u_{xxxxx} = 0$, $\beta\neq0$, given initial data which is analytic with a fixed radius. It is shown that the uniform radius of spatial analyticity of solutions at later time $t$ can decay no faster than $1/|t|$ as $|t|\rightarrow\infty$.


title: Dynamic exciton funneling by local strain control in a monolayer semiconductor

abstract: The ability to control excitons in semiconductors underlies numerous proposed applications, from excitonic circuits for computing and communications to polariton condensates to energy transport in photovoltaics. 2D semiconductors are particularly promising for room-temperature applications due to their large exciton binding energy. Their enormous stretchability gives rise to a strain-engineerable bandgap that has been used to induce static exciton flux in predetermined structures. However, dynamic control of exciton flux represents an outstanding challenge. Here, we introduce a method to tune the bandgap of suspended 2D semiconductors by applying a local strain gradient with a nanoscale tip. This strain allows us to locally and reversibly shift the exciton energy and to steer the exciton flux over micron-scale distances, as observed by wide-field imaging and time-resolved photoluminescence spectroscopy. We anticipate that the ability to strongly and dynamically modulate the bandgap of a semiconductor at the nanoscale not only marks an important experimental tool but will also open a broad range of new applications from information processing to energy conversion.


title: Finding Optimal Solutions With Neighborly Help

abstract: Can we efficiently compute optimal solutions to instances of a hard problem from optimal solutions to neighboring (i.e., locally modified) instances? For example, can we efficiently compute an optimal coloring for a graph from optimal colorings for all one-edge-deleted subgraphs? Studying such questions not only gives detailed insight into the structure of the problem itself, but also into the complexity of related problems; most notably graph theory's core notion of critical graphs (e.g., graphs whose chromatic number decreases under deletion of an arbitrary edge) and the complexity-theoretic notion of minimality problems (also called criticality problems, e.g., recognizing graphs that become 3-colorable when an arbitrary edge is deleted). We focus on two prototypical graph problems, Colorability and Vertex Cover. For example, we show that it is NP-hard to compute an optimal coloring for a graph from optimal colorings for all its one-vertex-deleted subgraphs, and that this remains true even when optimal solutions for all one-edge-deleted subgraphs are given. In contrast, computing an optimal coloring from all (or even just two) one-edge-added supergraphs is in P. We observe that Vertex Cover exhibits a remarkably different behavior, demonstrating the power of our model to delineate problems from each other more precisely on a structural level. Moreover, we provide a number of new complexity results for minimality and criticality problems. For example, we prove that Minimal-3-UnColorability is complete for DP (differences of NP sets), which was previously known only for the more amenable case of deleting vertices rather than edges. For Vertex Cover, we show that recognizing beta-vertex-critical graphs is complete for Theta_2^p (parallel access to NP), obtaining the first completeness result for a criticality problem for this class.


title: First-passage percolation in random planar maps and Tutte's bijection

abstract: We consider large random planar maps and study the first-passage percolation distance obtained by assigning independent identically distributed lengths to the edges. We consider the cases of quadrangulations and of general planar maps. In both cases, the first-passage percolation distance is shown to behave in large scales like a constant times the usual graph distance. We apply our method to the metric properties of the classical Tutte bijection between quadrangulations with $n$ faces and general planar maps with $n$ edges. We prove that the respective graph distances on the quadrangulation and on the associated general planar map are in large scales equivalent when $n \to \infty$.


title: K\"ahler-Einstein metrics on symmetric general arrangement varieties

abstract: We calculate Chow quotients of some families of symmetric \(T\)-varieties. In complexity two we obtain new examples of K\"ahler-Einstein metrics by bounding the symmetric alpha invariant of their orbifold quotients. As an additional application we determine the homeomorpism class of the orbit space of the compact torus action.


title: EasyCrash: Exploring Non-Volatility of Non-Volatile Memory for High Performance Computing Under Failures

abstract: Emerging non-volatile memory (NVM) is promising for building future HPC. Leveraging the non-volatility of NVM as main memory, we can restart the application using data objects remaining on NVM when the application crashes. This paper explores this solution to handle HPC under failures, based on the observation that many HPC applications have good enough intrinsic fault tolerance. To improve the possibility of successful recomputation with correct outcomes and ignorable performance loss, we introduce EasyCrash, a framework to decide how to selectively persist application data objects during application execution. Our evaluation shows that EasyCrash transforms 54% of crashes that cannot correctly recompute into the correct computation while incurring a negligible performance overhead (1.5% on average). Using EasyCrash and application intrinsic fault tolerance, 82% of crashes can successfully recompute. When EasyCrash is used with a traditional checkpoint scheme, it enables up to 24% improvement (15% on average) in system efficiency.


title: Scalable Online Survey Framework: from Sampling to Analysis

abstract: With the advancement in technology, raw event data generated by the digital world have grown tremendously. However, such data tend to be insufficient and noisy when it comes to measuring user intention or satisfaction. One effective way to measure user experience directly is through surveys. In particular, with the popularity of online surveys, extensive work has been put in to study this field. Surveys at LinkedIn play a major role in influencing product and marketing decisions and supporting our sales efforts. We run an increasing number of surveys that help us understand shifts in awareness and perceptions with regards to our own products and also to peer companies. As the need to survey grows, both sampling and analysis of surveys have become more challenging. Instead of simply multiplying the number of surveys each user takes, we need a scalable approach to collect enough and representative samples for each survey analysis while maintaining good user experience. In this paper, we start with discussions on how we handle multiple email surveys under such constraints. We then shift our discussions to challenges of in-product surveys and how we address them at LinkedIn through a survey study conducted across two mobile apps. Finally, we share how in-product surveys can be utilized as monitoring tools and connect surveys with A/B testing.


abstract: Relativistic invariance in Euclidean formulations of quantum mechanics is discussed. Relativistic treatments of quantum theory are needed to study hadronic systems at sub-hadronic distance scales. Euclidean formulations of relativistic quantum mechanics have some computational advantages. In the Euclidean representation the physical Hilbert space inner product is expressed in terms of Euclidean space-time variables with no need for any analytic continuation. The identification of the complex Euclidean group with the complex Poincar\'e group relates the infinitesimal generators of both groups. In this work explicit representations of the Poincar\'e generators in Euclidean space-time variables for all positive-mass positive-energy irreducible representations of the Poincar\'e group are derived. The commutation relations are checked, both hermiticity and self-adjointness are established, and reflection positivity of the kernels is verified.


title: Long Run Feedback in the Broker Call Money Market

abstract: I unravel the basic long run dynamics of the broker call money market, which is the pile of cash that funds margin loans to retail clients (read: continuous time Kelly gamblers). Call money is assumed to supply itself perfectly inelastically, and to continuously reinvest all principal and interest. I show that the relative size of the money market (that is, relative to the Kelly bankroll) is a martingale that nonetheless converges in probability to zero. The margin loan interest rate is a submartingale that converges in mean square to the choke price $r_\infty:=\nu-\sigma^2/2$, where $\nu$ is the asymptotic compound growth rate of the stock market and $\sigma$ is its annual volatility. In this environment, the gambler no longer beats the market asymptotically a.s. by an exponential factor (as he would under perfectly elastic supply). Rather, he beats the market asymptotically with very high probability (think 98%) by a factor (say 1.87, or 87% more final wealth) whose mean cannot exceed what the leverage ratio was at the start of the model (say, $2:1$). Although the ratio of the gambler's wealth to that of an equivalent buy-and-hold investor is a submartingale (always expected to increase), his realized compound growth rate converges in mean square to $\nu$. This happens because the equilibrium leverage ratio converges to $1:1$ in lockstep with the gradual rise of margin loan interest rates.

categories: econ.GN econ.TH q-fin.CP q-fin.EC q-fin.GN q-fin.PM


title: Inflowing Gas in the Central Parsec of M81

abstract: Spectroscopic observations of the Seyfert 1/Liner nucleus of M81, obtained recently with the Space Telescope Imaging Spectrograph (STIS) aboard the Hubble Space Telescope (HST), have revealed a UV--visible spectrum rich with emission lines of a variety of widths, ionization potentials, and critical densities, including several in the UV that have not previously been reported. Even at the highest angular resolution currently achievable with HST, the broad-line region of M81 cannot be uniquely defined on the basis of commonly used observables such as the full-width at half maximum of the emission lines, or ratios of various emission lines. Numerous broad forbidden lines complicate interpretation of the spectra. At least three separate line-emitting components are inferred. A large, highly ionized, low density, low metallicity H${^+}$ region producing the broad Balmer lines. Located within the H${^+}$ region are smaller condensations spanning a wide-range in density, and the source of forbidden line emission through collisional excitation of the respective ions. Intermingled with the H${^+}$ region and the condensations is a curious extended source of time-variable CIV ${\lambda}$ 1548 emission. Collectively, these observations can be qualitatively understood in the context of a shock excited jet cavity within a large H${^+}$ region that is photoionized by the central UV--X-ray source. The H${^+}$ region contains ${\sim}$ 500 M${\odot}$ of low metallicity gas that is dynamically unstable to inflow. At the current rate, the available H${^+}$ gas can sustain the advection dominated accretion flow that powers the central UV--X-ray source for 10$^{5}$ years.


title: Analyzing CART

abstract: Decision trees with binary splits are popularly constructed using Classification and Regression Trees (CART) methodology. For binary classification and regression models, this approach recursively divides the data into two near-homogenous daughter nodes according to a split point that maximizes the reduction in sum of squares error (the impurity) along a particular variable. This paper aims to study the bias and adaptive properties of regression trees constructed with CART. In doing so, we derive an interesting connection between the bias and the mean decrease in impurity (MDI) measure of variable importance---a tool widely used for model interpretability---defined as the sum of impurity reductions over all non-terminal nodes in the tree. In particular, we show that the probability content of a terminal subnode for a variable is small when the MDI for that variable is large and that this relationship is exponential---confirming theoretically that decision trees with CART have small bias and are adaptive to signal strength and direction. Finally, we apply these individual tree bounds to tree ensembles and show consistency of Breiman's random forests. The context is surprisingly general and applies to a wide variety of multivariable data generating distributions and regression functions. The main technical tool is an exact characterization of the conditional probability content of the daughter nodes arising from an optimal split, in terms of the partial dependence function and reduction in impurity.


title: A Note on Clockability for Ordinal Turing Machines

abstract: We study clockability for Ordinal Turing Machines (OTMs). In particular, we show that, in contrast to the situation for ITTMs, admissible ordinals can be OTM-clockable, that $\Sigma_{2}$-admissible ordinals are never OTM-clockable and that gaps in the OTM-clockable ordinals are always started by admissible limits of admissible ordinals.


title: Binary icosahedral quasicrystals of hard spheres in spherical confinement

abstract: The influence of geometry on the local and global packing of particles is important to many fundamental and applied research themes such as the structure and stability of liquids, crystals and glasses. Here, we show by experiments and simulations that a binary mixture of hard-sphere-like particles crystallizing into the MgZn2 Laves phase in bulk, spontaneously forms 3D icosahedral quasicrystals in slowly drying droplets. Moreover, the local symmetry of 70-80% of the particles changes to that of the MgCu2 Laves phase. Both of these findings are significant for photonic applications. If the stoichiometry deviates from that of the Laves phase, our experiments show that the crystallization of MgZn2 is hardly affected by the spherical confinement. Our simulations show that the quasicrystals nucleate away from the spherical boundary and grow along five-fold symmetric structures. Our findings not only open the way for particle-level studies of nucleation and growth of 3D quasicrystals, but also of binary crystallization.

categories: cond-mat.soft cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.stat-mech


title: Image to Images Translation for Multi-Task Organ Segmentation and Bone Suppression in Chest X-Ray Radiography

abstract: Chest X-ray radiography is one of the earliest medical imaging technologies and remains one of the most widely-used for diagnosis, screening, and treatment follow up of diseases related to lungs and heart. The literature in this field of research reports many interesting studies dealing with the challenging tasks of bone suppression and organ segmentation but performed separately, limiting any learning that comes with the consolidation of parameters that could optimize both processes. This study, and for the first time, introduces a multitask deep learning model that generates simultaneously the bone-suppressed image and the organ-segmented image, enhancing the accuracy of tasks, minimizing the number of parameters needed by the model and optimizing the processing time, all by exploiting the interplay between the network parameters to benefit the performance of both tasks. The architectural design of this model, which relies on a conditional generative adversarial network, reveals the process on how the well-established pix2pix network (image-to-image network) is modified to fit the need for multitasking and extending it to the new image-to-images architecture. The developed source code of this multitask model is shared publicly on Github as the first attempt for providing the two-task pix2pix extension, a supervised/paired/aligned/registered image-to-images translation which would be useful in many multitask applications. Dilated convolutions are also used to improve the results through a more effective receptive field assessment. The comparison with state-of-the-art algorithms along with ablation study and a demonstration video are provided to evaluate efficacy and gauge the merits of the proposed approach.


title: Holographic Principle, Cosmological Constant and Cyclic Cosmology

abstract: The holographic principle provides a deep insight into quantum gravity and resolves the fine-tuning crisis concerning the cosmological constant. Holographic dark energy introduces new ultra-violet (UV) and infra-red (IR) cutoffs into quantum gravity which are necessarily strongly related. The equation of state for dark energy $\omega = p/\rho$ is discussed from the holographic point of view. The phantom option of $\omega < -1$ is resurrected, as in an earlier cyclic cosmology. Such a cyclic model can, however, equally use the cosmological constant with $\omega = -1$.


title: Fixed-time Control under Spatiotemporal and Input Constraints: A Quadratic Program Based Approach

abstract: In this paper, we present a control synthesis framework for a general class of nonlinear, control-affine systems under spatiotemporal and input constraints. First, we study the problem of fixed-time convergence in the presence of input constraints. The relation between the domain of attraction for fixed-time stability with respect to input constraints and the required time of convergence is established. It is shown that increasing the control authority or the required time of convergence can expand the domain of attraction for fixed-time stability. Then, we consider the problem of finding a control input that confines the closed-loop system trajectories in a safe set and steers them to a goal set within a fixed time. To this end, we present a Quadratic Program (QP) formulation to compute the corresponding control input. We use slack variables to guarantee feasibility of the proposed QP under input constraints. Furthermore, when strict complementary slackness holds, we show that the solution of the QP is a continuous function of the system states, and establish uniqueness of closed-loop solutions to guarantee forward invariance using Nagumo's theorem. We present two case studies, an example of adaptive cruise control problem and an instance of a two-robot motion planning problem, to corroborate our proposed methods.


title: On linear structure of non-commutative operator graphs

abstract: We continue the study of non-commutative operator graphs generated by resolutions of identity covariant with respect to unitary actions of the circle group and the Heisenber-Weyl group as well. It is shown that the graphs generated by the circle group has the system of unitary generators fulfilling permutations of basis vectors. For the graph generated by the Heisenberg-Weyl group the explicit formula for a dimension is given. Thus, we found a new description of the linear structure for the operator graphs introduced in our previous works.


title: Efficient Analysis of Unambiguous Automata Using Matrix Semigroup Techniques

abstract: We introduce a novel technique to analyse unambiguous B\"uchi automata quantitatively, and apply this to the model checking problem. It is based on linear-algebra arguments that originate from the analysis of matrix semigroups with constant spectral radius. This method can replace a combinatorial procedure that dominates the computational complexity of the existing procedure by Baier et al. We analyse the complexity in detail, showing that, in terms of the set $Q$ of states of the automaton, the new algorithm runs in time $O(|Q|^4)$, improving on an efficient implementation of the combinatorial algorithm by a factor of $|Q|$.


title: Density-based Clustering with Best-scored Random Forest

abstract: Single-level density-based approach has long been widely acknowledged to be a conceptually and mathematically convincing clustering method. In this paper, we propose an algorithm called "best-scored clustering forest" that can obtain the optimal level and determine corresponding clusters. The terminology "best-scored" means to select one random tree with the best empirical performance out of a certain number of purely random tree candidates. From the theoretical perspective, we first show that consistency of our proposed algorithm can be guaranteed. Moreover, under certain mild restrictions on the underlying density functions and target clusters, even fast convergence rates can be achieved. Last but not least, comparisons with other state-of-the-art clustering methods in the numerical experiments demonstrate accuracy of our algorithm on both synthetic data and several benchmark real data sets.


title: An Empirical Comparison of FAISS and FENSHSES for Nearest Neighbor Search in Hamming Space

abstract: In this paper, we compare the performances of FAISS and FENSHSES on nearest neighbor search in Hamming space--a fundamental task with ubiquitous applications in nowadays eCommerce. Comprehensive evaluations are made in terms of indexing speed, search latency and RAM consumption. This comparison is conducted towards a better understanding on trade-offs between nearest neighbor search systems implemented in main memory and the ones implemented in secondary memory, which is largely unaddressed in literature.


title: Audio-Visual Kinship Verification

abstract: Visual kinship verification entails confirming whether or not two individuals in a given pair of images or videos share a hypothesized kin relation. As a generalized face verification task, visual kinship verification is particularly difficult with low-quality found Internet data. Due to uncontrolled variations in background, pose, facial expression, blur, illumination and occlusion, state-of-the-art methods fail to provide high level of recognition accuracy. As with many other visual recognition tasks, kinship verification may benefit from combining visual and audio signals. However, voice-based kinship verification has received very little prior attention. We hypothesize that the human voice contains kin-related cues that are complementary to visual cues. In this paper, we address, for the first time, the use of audio-visual information from face and voice modalities to perform kinship verification. We first propose a new multi-modal kinship dataset, called TALking KINship (TALKIN), that contains several pairs of Internet-quality video sequences. Using TALKIN, we then study the utility of various kinship verification methods including traditional local feature based methods, statistical methods and more recent deep learning approaches. Then, early and late fusion methods are evaluated on the TALKIN dataset for the study of kinship verification with both face and voice modalities. Finally, we propose a deep Siamese fusion network with contrastive loss for multi-modal fusion of kinship relations. Extensive experiments on the TALKIN dataset indicate that by combining face and voice modalities, the proposed Siamese network can provide a significantly higher level of accuracy compared to baseline uni-modal and multi-modal fusion techniques. Experimental results also indicate that audio (vocal) information is complementary (to facial information) and useful for kinship verification.


title: Dirichlet energy-minimizers with analytic boundary

abstract: In this paper, we consider multi-valued graphs with a prescribed real analytic interface that minimize the Dirichlet energy. Such objects arise as a linearized model of area minimizing currents with real analytic boundaries and our main result is that their singular set is discrete in 2 dimensions. This confirms (and provides a first step to) a conjecture by B. White \cite{White97} that area minimizing $2$-dimensional currents with real analytic boundaries have a finite number of singularities. We also show that, in any dimension, Dirichlet energy-minimizers with a $C^1$ boundary interface are H\"older continuous at the interface.


title: Bayesian waveform-based calibration of high-pressure acoustic emission systems with ball drop measurements

abstract: Acoustic emission (AE) is a widely used technology to study source mechanisms and material properties during high-pressure rock failure experiments. It is important to understand the physical quantities that acoustic emission sensors measure, as well as the response of these sensors as a function of frequency. This study calibrates the newly built AE system in the MIT Rock Physics Laboratory using a ball-bouncing system. Full waveforms of multi-bounce events due to ball drops are used to infer the transfer function of lead zirconate titanate (PZT) sensors in high pressure environments. Uncertainty in the sensor transfer functions is quantified using a waveform-based Bayesian approach. The quantification of \textit{in situ} sensor transfer functions makes it possible to apply full waveform analysis for acoustic emissions at high pressures.


title: DynoPlan: Combining Motion Planning and Deep Neural Network based Controllers for Safe HRL

abstract: Many realistic robotics tasks are best solved compositionally, through control architectures that sequentially invoke primitives and achieve error correction through the use of loops and conditionals taking the system back to alternative earlier states. Recent end-to-end approaches to task learning attempt to directly learn a single controller that solves an entire task, but this has been difficult for complex control tasks that would have otherwise required a diversity of local primitive moves, and the resulting solutions are also not easy to inspect for plan monitoring purposes. In this work, we aim to bridge the gap between hand designed and learned controllers, by representing each as an option in a hybrid hierarchical Reinforcement Learning framework - DynoPlan. We extend the options framework by adding a dynamics model and the use of a nearness-to-goal heuristic, derived from demonstrations. This translates the optimization of a hierarchical policy controller to a problem of planning with a model predictive controller. By unrolling the dynamics of each option and assessing the expected value of each future state, we can create a simple switching controller for choosing the optimal policy within a constrained time horizon similarly to hill climbing heuristic search. The individual dynamics model allows each option to iterate and be activated independently of the specific underlying instantiation, thus allowing for a mix of motion planning and deep neural network based primitives. We can assess the safety regions of the resulting hybrid controller by investigating the initiation sets of the different options, and also by reasoning about the completeness and performance guarantees of the underpinning motion planners.


title: Kerr Black Holes as Elementary Particles

abstract: Long ago, Newman and Janis showed that a complex deformation $z\rightarrow z+i a$ of the Schwarzschild solution produces the Kerr solution. The underlying explanation for this relationship has remained obscure. The complex deformation has an electromagnetic counterpart: by shifting the Coloumb potential, we obtain the EM field of a certain rotating charge distribution which we term $\sqrt{\rm Kerr}$. In this note, we identify the origin of this shift as arising from the exponentiation of spin operators for the recently defined "minimally coupled" three-particle amplitudes of spinning particles coupled to gravity, in the large-spin limit. We demonstrate this by studying the impulse imparted to a test particle in the background of the heavy spinning particle. We first consider the electromagnetic case, where the impulse due to $\sqrt{\rm Kerr}$ is reproduced by a charged spinning particle; the shift of the Coloumb potential is matched to the exponentiated spin-factor appearing in the amplitude. The known impulse due to the Kerr black hole is then trivially derived from the gravitationally coupled spinning particle via the double copy.


title: LMVP: Video Predictor with Leaked Motion Information

abstract: We propose a Leaked Motion Video Predictor (LMVP) to predict future frames by capturing the spatial and temporal dependencies from given inputs. The motion is modeled by a newly proposed component, motion guider, which plays the role of both learner and teacher. Specifically, it {\em learns} the temporal features from real data and {\em guides} the generator to predict future frames. The spatial consistency in video is modeled by an adaptive filtering network. To further ensure the spatio-temporal consistency of the prediction, a discriminator is also adopted to distinguish the real and generated frames. Further, the discriminator leaks information to the motion guider and the generator to help the learning of motion. The proposed LMVP can effectively learn the static and temporal features in videos without the need for human labeling. Experiments on synthetic and real data demonstrate that LMVP can yield state-of-the-art results.


title: The effect of geometry parameters and flow characteristics on erosion and sedimentation in channels junction using finite volume method

abstract: One of the most critical problems in the river engineering field is scouring, sedimentation and morphology of a river bed. In this paper, a finite volume method FORTRAN code is provided and used. The code is able to model the sedimentation. The flow and sediment were modeled at the interception of the two channels. It is applied an experimental model to evaluate the results. Regarding the numerical model, the effects of geometry parameters such as proportion of secondary channel to main channel width and intersection angle and also hydraulic conditionals like secondary to main channel discharge ratio and inlet flow Froude number were studied on bed topographical and flow pattern. The numerical results show that the maximum height of bed increased to 32 percent as the discharge ratio reaches to 51 percent, on average. It is observed that the maximum height of sedimentation decreases by declining in main channel to secondary channel Froude number ratio. On the assessment of the channel width, velocity and final bed height variations have changed by given trend, in all the ratios. Also, increasing in the intersection angle accompanied by decreasing in flow velocity variations along the channel. The pattern of velocity and topographical bed variations are also constant in any studied angles.


title: Reciprocity and representation theorems for flux- and field-normalised decomposed wave fields

abstract: We consider wave propagation problems in which there is a preferred direction of propagation. To account for propagation in preferred directions, the wave equation is decomposed into a set of coupled equations for waves that propagate in opposite directions along the preferred axis. This decomposition is not unique. We discuss flux-normalised and field-normalised decomposition in a systematic way, analyse the symmetry properties of the decomposition operators and use these symmetry properties to derive reciprocity theorems for the decomposed wave fields, for both types of normalisation. Based on the field-normalised reciprocity theorems, we derive representation theorems for decomposed wave fields. In particular we derive double- and single-sided Kirchhoff-Helmholtz integrals for forward and backward propagation of decomposed wave fields. The single-sided Kirchhoff-Helmholtz integrals for backward propagation of field-normalised decomposed wave fields find applications in reflection imaging, accounting for multiple scattering.


title: Remote Estimation of Free-Flow Speeds

abstract: We propose an automated method to estimate a road segment's free-flow speed from overhead imagery and road metadata. The free-flow speed of a road segment is the average observed vehicle speed in ideal conditions, without congestion or adverse weather. Standard practice for estimating free-flow speeds depends on several road attributes, including grade, curve, and width of the right of way. Unfortunately, many of these fine-grained labels are not always readily available and are costly to manually annotate. To compensate, our model uses a small, easy to obtain subset of road features along with aerial imagery to directly estimate free-flow speed with a deep convolutional neural network (CNN). We evaluate our approach on a large dataset, and demonstrate that using imagery alone performs nearly as well as the road features and that the combination of imagery with road features leads to the highest accuracy.


title: Coded Distributed Computing: Performance Limits and Code Designs

abstract: We consider the problem of coded distributed computing where a large linear computational job, such as a matrix multiplication, is divided into $k$ smaller tasks, encoded using an $(n,k)$ linear code, and performed over $n$ distributed nodes. The goal is to reduce the average execution time of the computational job. We provide a connection between the problem of characterizing the average execution time of a coded distributed computing system and the problem of analyzing the error probability of codes of length $n$ used over erasure channels. Accordingly, we present closed-form expressions for the execution time using binary random linear codes and the best execution time any linear-coded distributed computing system can achieve. It is also shown that there exist good binary linear codes that attain, asymptotically, the best performance any linear code, not necessarily binary, can achieve. We also investigate the performance of coded distributed computing systems using polar and Reed-Muller (RM) codes that can benefit from low-complexity decoding, and superior performance, respectively, as well as explicit constructions. The proposed framework in this paper can enable efficient designs of distributed computing systems given the rich literature in the channel coding theory.


title: Implicitly Learning to Reason in First-Order Logic

abstract: We consider the problem of answering queries about formulas of first-order logic based on background knowledge partially represented explicitly as other formulas, and partially represented as examples independently drawn from a fixed probability distribution. PAC semantics, introduced by Valiant, is one rigorous, general proposal for learning to reason in formal languages: although weaker than classical entailment, it allows for a powerful model theoretic framework for answering queries while requiring minimal assumptions about the form of the distribution in question. To date, however, the most significant limitation of that approach, and more generally most machine learning approaches with robustness guarantees, is that the logical language is ultimately essentially propositional, with finitely many atoms. Indeed, the theoretical findings on the learning of relational theories in such generality have been resoundingly negative. This is despite the fact that first-order logic is widely argued to be most appropriate for representing human knowledge. In this work, we present a new theoretical approach to robustly learning to reason in first-order logic, and consider universally quantified clauses over a countably infinite domain. Our results exploit symmetries exhibited by constants in the language, and generalize the notion of implicit learnability to show how queries can be computed against (implicitly) learned first-order background knowledge.


title: Primal-dual fast gradient method with a model

abstract: In this work we consider a possibility to use the conception of $(\delta, L)$-model of a function for optimization tasks, whereby solving a primal problem there is a necessity to recover a solution of a dual problem. The conception of $(\delta, L)$-model is based on the conception of $(\delta, L)$-oracle which was proposed by Devolder-Glineur-Nesterov, herewith the authors proposed approximate a function with an upper bound using a convex quadratic function with some additive noise $\delta$. They managed to get convex quadratic upper bounds with noise even for nonsmooth functions. The conception of $(\delta, L)$-model continues this idea by using instead of a convex quadratic function a more complex convex function in an upper bound. Possibility to recover the solution of a dual problem gives great benefits in different problems, for instance, in some cases, it is faster to find a solution in a primal problem than in a dual problem. Note that primal-dual methods are well studied, but usually each class of optimization problems has its own primal-dual method. Our goal is to develop a method which can find solutions in different classes of optimization problems. This is realized through the use of the conception of $(\delta, L)$-model and adaptive structure of our methods. Thereby, we developed primal-dual adaptive gradient method and fast gradient method with $(\delta, L)$-model and proved convergence rates of the methods, moreover, for some classes of optimization problems the rates are optimal. The main idea is the following: we find a dual solution to an approximation of a primal problem using the conception of $(\delta, L)$-model. It is much easier to find a solution to an approximated problem, however, we have to do it in each step of our method, thereby the principle of "divide and conquer" is realized.


title: A lift of the Seiberg-Witten equations to Kaluza-Klein 5-manifolds

abstract: We consider Riemannian 4-manifolds $(X,g_X)$ with a Spin^c-structure and a suitable circle bundle $Y$ over $X$ such that the Spin^c-structure on $X$ lifts to a spin structure on $Y$. With respect to these structures a spinor $\phi$ on $X$ lifts to an untwisted spinor $\psi$ on $Y$ and a U(1)-gauge field $A$ for the Spin^c-structure can be absorbed into a Kaluza-Klein metric $g_Y^A$ on $Y$. We show that irreducible solutions $(A,\phi)$ to the Seiberg-Witten equations on $(X,g_X)$ for the given Spin^c-structure are equivalent to irreducible solutions $\psi$ of a Dirac equation with cubic non-linearity on the Kaluza-Klein circle bundle $(Y,g_Y^A)$. As an application we consider solutions to the equations in the case of Sasaki 5-manifolds which are circle bundles over Kaehler-Einstein surfaces.


title: CMRNet: Camera to LiDAR-Map Registration

abstract: In this paper we present CMRNet, a realtime approach based on a Convolutional Neural Network to localize an RGB image of a scene in a map built from LiDAR data. Our network is not trained in the working area, i.e. CMRNet does not learn the map. Instead it learns to match an image to the map. We validate our approach on the KITTI dataset, processing each frame independently without any tracking procedure. CMRNet achieves 0.27m and 1.07deg median localization accuracy on the sequence 00 of the odometry dataset, starting from a rough pose estimate displaced up to 3.5m and 17deg. To the best of our knowledge this is the first CNN-based approach that learns to match images from a monocular camera to a given, preexisting 3D LiDAR-map.


title: Modeling of CoRoT and Spitzer lightcurves in NGC 2264 caused by an optically thick warp

abstract: Aims: We present an analysis of simultaneously observed CoRoT and Spitzer lightcurves for $4$ systems in the stellar forming region NGC 2264: Mon-660, Mon-811, Mon-1140 and Mon-1308. These objects share in common a high resemblance between the optical and infrared lightcurves, such that the mechanism responsible to produce them is the same. The aim of this paper is to explain both lightcurves simultaneously with only one mechanism. Methods: We have modeled the infrared emission as coming from a warp composed of an optically thick wall and an optically thick asymmetric disk beyond this location. We have modeled the optical emission mainly by partial stellar occultation by the warp. Results: The magnitude amplitude of the CoRoT and Spitzer observations for all the objects can be described with the emission coming from the system components. The difference between them is the value of the disk flux compared with the wall flux and the azimuthal variations of the former. This result points out the importance of the hydrodynamical interaction between the stellar magnetic field and the disk. Conclusions: CoRoT and Spitzer lightcurves for the stellar systems Mon-660, Mon-811, Mon-1140 and Mon-1308 can be simultaneously explained using the emission coming from an asymmetric disk and emission with stellar occultation by an optically thick wall.


title: Cooling Rate Effects on the Structure of 45S5 Bioglass: Computational and Experimental Evidence of Si--P Avoidance

abstract: Due to its ability to bond with living tissues upon dissolution, 45S5 bioglass and related compositions are promising materials for the replacement, regeneration, and repair of hard tissues in the human body. However, the details of their atomic structure remain unclear. This is partially due to the non-equilibrium nature of glasses, as their non-crystalline structure is highly dependent on their thermal history, namely, the cooling rate used during quenching. Herein, using molecular dynamics (MD) simulations and magic angle spinning nuclear magnetic resonance (MAS-NMR) spectroscopy experiments, we investigate the structure of the nominal 45S5 bioglass composition prepared using cooling rates ranging over several orders of magnitude. We show that the simulations results are in very good agreement with experimental data, provided that they are extrapolated toward lower cooling rates achieved in experiments. These results highlight that previously reported inconsistencies between simulations and experiments stem from the difference in cooling rate, thereby addressing one of the longstanding questions on the structure of bioglass. Based on these results, we demonstrate the existence of a Si--P avoidance behavior, which may be key in controlling the bioactivity of 45S5 bioglass.


title: GANalyze: Toward Visual Definitions of Cognitive Image Properties

abstract: We introduce a framework that uses Generative Adversarial Networks (GANs) to study cognitive properties like memorability, aesthetics, and emotional valence. These attributes are of interest because we do not have a concrete visual definition of what they entail. What does it look like for a dog to be more or less memorable? GANs allow us to generate a manifold of natural-looking images with fine-grained differences in their visual attributes. By navigating this manifold in directions that increase memorability, we can visualize what it looks like for a particular generated image to become more or less memorable. The resulting ``visual definitions" surface image properties (like ``object size") that may underlie memorability. Through behavioral experiments, we verify that our method indeed discovers image manipulations that causally affect human memory performance. We further demonstrate that the same framework can be used to analyze image aesthetics and emotional valence. Visit the GANalyze website at http://ganalyze.csail.mit.edu/.


title: X-ray and gamma-ray limits on the primordial black hole abundance from Hawking radiation

abstract: The non-observation of Hawking radiation from primordial black holes of $10^{16}$g sets a conservative strong bound on their cosmological abundance. We revisit this bound and show how it can be improved (both in mass reach and strength) by an adequate modeling of the combined AGN and blazar emission in the MeV range. We also estimate the sensitivity to the primordial black hole abundance of a future X-ray experiment capable of identifying a significantly larger number of astrophysical sources contributing to the diffuse background in this energy range.


title: Trajectory of Alternating Direction Method of Multipliers and Adaptive Acceleration

abstract: The alternating direction method of multipliers (ADMM) is one of the most widely used first-order optimisation methods in the literature owing to its simplicity, flexibility and efficiency. Over the years, numerous efforts are made to improve the performance of the method, such as the inertial technique. By studying the geometric properties of ADMM, we discuss the limitations of current inertial accelerated ADMM and then present and analyze an adaptive acceleration scheme for the method. Numerical experiments on problems arising from image processing, statistics and machine learning demonstrate the advantages of the proposed acceleration approach.


title: Divide and Couple: Using Monte Carlo Variational Objectives for Posterior Approximation

abstract: Recent work in variational inference (VI) uses ideas from Monte Carlo estimation to tighten the lower bounds on the log-likelihood that are used as objectives. However, there is no systematic understanding of how optimizing different objectives relates to approximating the posterior distribution. Developing such a connection is important if the ideas are to be applied to inference-i.e., applications that require an approximate posterior and not just an approximation of the log-likelihood. Given a VI objective defined by a Monte Carlo estimator of the likelihood, we use a "divide and couple" procedure to identify augmented proposal and target distributions. The divergence between these is equal to the gap between the VI objective and the log-likelihood. Thus, after maximizing the VI objective, the augmented variational distribution may be used to approximate the posterior distribution.


title: Spectral and transport properties of a $\mathcal{PT}$-symmetric tight-binding chain with gain and loss

abstract: We derive a continuity equation to study transport properties in a $\mathcal{PT}$-symmetric tight-binding chain with gain and loss in symmetric configurations. This allows us to identify the density fluxes in the system, and to define a transport coefficient to characterize the efficiency of transport of each state. These quantities are studied explicitly using analytical expressions for the eigenvalues and eigenvectors of the system. We find that in states with broken $\mathcal{PT}$-symmetry, transport is inefficient, in the sense that either inflow exceeds outflow and density accumulates within the system, or outflow exceeds inflow, and the system becomes depleted. We also report the appearance of two subsets of interesting eigenstates whose eigenvalues are independent on the strength of the coupling to gain and loss. We call these opaque and transparent states. Opaque states are decoupled from the contacts and there is no transport; transparent states exhibit always efficient transport. Interestingly, the appearance of such eigenstates is connected with the divisors of the length of the system plus one and the position of the contacts. Thus the number of opaque and transparent states varies very irregularly.


title: The contact in the unitary Fermi gas across the superfluid phase transition

abstract: A quantity known as the contact plays a fundamental role in quantum many-body systems with short-range interactions. The determination of the temperature dependence of the contact for the unitary Fermi gas of infinite scattering length has been a major challenge, with different calculations yielding qualitatively different results. Here we use finite-temperature auxiliary-field quantum Monte Carlo (AFMC) methods on the lattice within the canonical ensemble to calculate the temperature dependence of the contact for the homogeneous spin-balanced unitary Fermi gas. We extrapolate to the continuum limit for 40, 66, and 114 particles. We observe a dramatic decrease in the contact as the superfluid critical temperature is approached from below, followed by a gradual weak decrease as the temperature increases in the normal phase. Our results are in excellent agreement with the most recent precision ultracold atomic gas experiments. We also present results for the energy of the unitary gas as a function of temperature in the continuum limit.

categories: cond-mat.quant-gas cond-mat.supr-con hep-lat nucl-th


title: Query-driven PAC-Learning for Reasoning

abstract: We consider the problem of learning rules from a data set that support a proof of a given query, under Valiant's PAC-Semantics. We show how any backward proof search algorithm that is sufficiently oblivious to the contents of its knowledge base can be modified to learn such rules while it searches for a proof using those rules. We note that this gives such algorithms for standard logics such as chaining and resolution.


title: Yukawa Hierarchies in Global F-theory Models

abstract: We argue that global F-theory compactifications to four dimensions generally exhibit higher rank Yukawa matrices from multiple geometric contributions known as Yukawa points. The holomorphic couplings furthermore have large hierarchies for generic complex structure moduli. Unlike local considerations, the compact setup realizes these features all through geometry, and requires no instanton corrections. As an example, we consider a concrete toy model with $SU(5) \times U(1)$ gauge symmetry. From the geometry, we find two Yukawa points for the ${\bf 10}_{-2} \, \bar{\bf 5}_6 \, \bar{\bf 5}_{-4}$ coupling, producing a rank two Yukawa matrix. Our methods allow us to track all complex structure dependencies of the holomorphic couplings and study the ratio numerically. This reveals hierarchies of ${\cal O}(10^5)$ and larger on a full-dimensional subspace of the moduli space.


title: Assembly line balancing with task division

abstract: In a commonly-used version of the Simple Assembly Line Balancing Problem (SALBP-1) tasks are assigned to stations along an assembly line with a fixed cycle time in order to minimize the required number of stations. It has traditionally been assumed that the total work needed for each product unit has been partitioned into economically indivisible tasks. However, in practice, it is sometimes possible to divide particular tasks in limited ways at additional time penalty cost. Despite the penalties, task division where possible, now and then leads to a reduction in the minimum number of stations. Deciding which allowable tasks to divide creates a new assembly line balancing problem, TDALBP (Task Division Assembly Line Balancing Problem). We propose a mathematical model of the TDALBP, an exact solution procedure for it and present promising computational results for the adaptation of some classical SALBP instances from the research literature. The results demonstrate that the TDALBP sometimes has the potential to significantly improve assembly line performance.


title: Metaheuristics optimized feedforward neural networks for efficient stock price prediction

abstract: The prediction of stock prices is an important task in economics, investment and making financial decisions. This has, for decades, spurred the interest of many researchers to make focused contributions to the design of accurate stock price predictive models; of which some have been utilized to predict the next day opening and closing prices of the stock indices. This paper proposes the design and implementation of a hybrid symbiotic organisms search trained feedforward neural network model for effective and accurate stock price prediction. The symbiotic organisms search algorithm is used as an efficient optimization technique to train the feedforward neural networks, while the resulting training process is used to build a better stock price prediction model. Furthermore, the study also presents a comparative performance evaluation of three different stock price forecasting models; namely, the particle swarm optimization trained feedforward neural network model, the genetic algorithm trained feedforward neural network model and the well-known ARIMA model. The system developed in support of this study utilizes sixteen stock indices as time series datasets for training and testing purpose. Three statistical evaluation measures are used to compare the results of the implemented models, namely the root mean squared error, the mean absolute percentage error and the mean absolution deviation. The computational results obtained reveal that the symbiotic organisms search trained feedforward neural network model exhibits outstanding predictive performance compared to the other models. However, the performance study shows that the three metaheuristics trained feedforward neural network models have promising predictive competence for solving problems of high dimensional nonlinear time series data, which are difficult to capture by traditional models.


title: X-ray photon scattering at a focused high-intensity laser pulse

abstract: We study x-ray photon scattering in the head-on collision of an XFEL pulse and a focused high-intensity laser pulse, described as paraxial Laguerre-Gaussian beam of arbitrary mode composition. For adequately chosen relative orientations of the polarization vectors of the colliding laser fields, this gives rise to a vacuum birefringence effect manifesting itself in polarization flipped signal photons. Throughout this article the XFEL is assumed to be mildly focused to a waist larger than that of the high-intensity laser beam. As previously demonstrated for the special case of a fundamental paraxial Gaussian beam, this scenario is generically accompanied by a scattering phenomenon of x-ray energy signal photons outside the forward cone of the XFEL beam, potentially assisting the detection of the effect in experiment. Here, we study the fate of the x-ray scattering signal under exemplary deformations of the transverse focus profile of the high-intensity pump.


title: A Schr\"odinger potential involving $x^\frac{2}{3}$ and centrifugal-barrier terms conditionally integrable in terms of the confluent hypergeometric functions

abstract: The solution of the one-dimensional Schr\"odinger equation for a potential involving an attractive $x^\frac{2}{3}$ and a repulsive centrifugal-barrier $\sim x^{-2}$ terms is presented in terms of the non-integer-order Hermite functions. The potential belongs to one of the five bi-confluent Heun families. This is a conditionally integrable potential in that the strength of the centrifugal-barrier term is fixed. The general solution of the problem is composed using fundamental solutions each of which presents an irreducible linear combination of two Hermite functions of a scaled and shifted argument. The potential presents an infinitely extended confining well defined on the positive semi-axis and sustains infinitely many bound states.


title: On Multi-Agent Learning in Team Sports Games

abstract: In recent years, reinforcement learning has been successful in solving video games from Atari to Star Craft II. However, the end-to-end model-free reinforcement learning (RL) is not sample efficient and requires a significant amount of computational resources to achieve superhuman level performance. Model-free RL is also unlikely to produce human-like agents for playtesting and gameplaying AI in the development cycle of complex video games. In this paper, we present a hierarchical approach to training agents with the goal of achieving human-like style and high skill level in team sports games. While this is still work in progress, our preliminary results show that the presented approach holds promise for solving the posed multi-agent learning problem.


title: A note on locally optimal designs for generalized linear models with restricted support

abstract: Optimal designs for generalized linear models require a prior knowledge of the regression parameters. At certain values of the parameters we propose particular assumptions which allow to derive a locally optimal design for a model without intercept from a locally optimal design for the corresponding model with intercept and vice versa. Applications to Poisson and logistic models and Extensions to nonlinear models are provided.


title: Evidence for Sub-Chandrasekhar Type Ia Supernovae from Stellar Abundances in Dwarf Galaxies

abstract: There is no consensus on the progenitors of Type Ia supernovae (SNe Ia) despite their importance for cosmology and chemical evolution. We address this question by using our previously published catalogs of Mg, Si, Ca, Cr, Fe, Co, and Ni abundances in dwarf galaxy satellites of the Milky Way to constrain the mass at which the white dwarf explodes during a typical SN Ia. We fit a simple bi-linear model to the evolution of [X/Fe] with [Fe/H], where X represents each of the elements mentioned above. We use the evolution of [Mg/Fe] coupled with theoretical supernova yields to isolate what fraction of the elements originated in SNe Ia. Then, we infer the [X/Fe] yield of SNe Ia for all of the elements except Mg. We compare these observationally inferred yields to recent theoretical predictions for two classes of Chandrasekhar-mass (M_Ch) SN Ia as well as sub-M_Ch SNe Ia. Most of the inferred SN Ia yields are consistent with all of the theoretical models, but [Ni/Fe] is consistent only with sub-M_Ch models. We conclude that the dominant type of SN Ia in ancient dwarf galaxies is the explosion of a sub-M_Ch white dwarf. The Milky Way and dwarf galaxies with extended star formation histories have higher [Ni/Fe] abundances, which could indicate that the dominant class of SN Ia is different for galaxies where star formation lasted for at least several Gyr.


title: Probing Proton Structure at the Large Hadron electron Collider

abstract: For the foreseeable future, the exploration of the high-energy frontier will be the domain of the Large Hadron Collider (LHC). Of particular significance will be its high-luminosity upgrade (HL-LHC), which will operate until the mid-2030s. In this endeavour, for the full exploitation of the HL-LHC physics potential an improved understanding of the parton distribution functions (PDFs) of the proton is critical. The HL-LHC program would be uniquely complemented by the proposed Large Hadron electron Collider (LHeC), a high-energy lepton-proton and lepton-nucleus collider based at CERN. In this work, we build on our recent PDF projections for the HL-LHC to assess the constraining power of the LHeC measurements of inclusive and heavy quark structure functions. We find that the impact of the LHeC would be significant, reducing PDF uncertainties by up to an order of magnitude in comparison to state-of-the-art global fits. In comparison to the HL-LHC projections, the PDF constraints from the LHeC are in general more significant for small and intermediate values of the momentum fraction x. At higher values of x, the impact of the LHeC and HL-LHC data is expected to be of a comparable size, with the HL-LHC constraints being more competitive in some cases, and the LHeC ones in others. Our results illustrate the encouraging complementarity of the HL-LHC and the LHeC in terms of charting the quark and gluon structure of the proton.


title: Discovery of an equal-mass "twin" binary population reaching 1000+ AU separations

abstract: We use a homogeneous catalog of 42,000 main-sequence wide binaries identified by Gaia to measure the mass ratio distribution, p(q), of binaries with primary masses $0.1<M_1/M_{\odot}<2.5$, mass ratios $0.1 \lesssim q<1$, and separations $50<s/{\rm AU}<50,000$. A well-understood selection function allows us to constrain p(q) in 35 independent bins of primary mass and separation, with hundreds to thousands of binaries in each bin. Our investigation reveals a sharp excess of equal-mass "twin" binaries that is statistically significant out to separations of 1,000 to 10,000 AU, depending on primary mass. The excess is narrow: a steep increase in p(q) at $0.95 \lesssim q<1$, with no significant excess at $q\lesssim 0.95$. A range of tests confirm the signal is real, not a data artifact or selection effect. Combining the Gaia constraints with those from close binaries, we show that the twin excess decreases with increasing separation, but its width ($q\gtrsim 0.95$) is constant over $0.01<a/{\rm AU}<10,000$. The wide twin population would be difficult to explain if the components of all wide binaries formed via core fragmentation, which is not expected to produce strongly correlated component masses. We conjecture that wide twins formed at closer separations ($a \lesssim 100$ AU), likely via accretion from circumbinary disks, and were subsequently widened by dynamical interactions in their birth environments. The separation-dependence of the twin excess then constrains the efficiency of dynamical widening and disruption of binaries in young clusters. We also constrain p(q) across $0.1 \lesssim q<1$. Besides changes in the twin fraction, p(q) is independent of separation at fixed primary mass over $100 \lesssim s/{\rm AU} < 50,000$. It is flatter than expected for random pairings from the IMF but more bottom-heavy for wide binaries than for binaries with $a\lesssim$100 AU.


title: The geometry of the gas surrounding the Central Molecular Zone: on the origin of localised molecular clouds with extreme velocity dispersions

abstract: Observations of molecular gas near the Galactic centre ($| l | < 10^\circ$, $| b | < 1^\circ$) reveal the presence of a distinct population of enigmatic compact clouds which are characterised by extreme velocity dispersions ($\Delta v > 100\, \rm km/s$). These Extended Velocity Features (EVFs) are very prominent in the datacubes and dominate the kinematics of molecular gas just outside the Central Molecular Zone (CMZ). The prototypical example of such a cloud is Bania Clump 2. We show that similar features are naturally produced in simulations of gas flow in a realistic barred potential. We analyse the structure of the features obtained in the simulations and use this to interpret the observations. We find that the features arise from collisions between material that has been infalling rapidly along the dust lanes of the Milky Way bar and material that belongs to one of the following two categories: (i) material that has `overshot' after falling down the dust lanes on the opposite side; (ii) material which is part of the CMZ. Both types of collisions involve gas with large differences in the line-of-sight velocities, which is what produces the observed extreme velocity dispersions. Examples of both categories can be identified in the observations. If our interpretation is correct, we are directly witnessing (a) collisions of clouds with relative speeds of $\sim 200\, \rm km/s$ and (b) the process of accretion of fresh gas onto the CMZ.


title: Evidence for Low Radiative Efficiency or Highly Obscured Growth of $z>7$ Quasars

abstract: The supermassive black holes (SMBHs) observed at the centers of all massive galaxies are believed to have grown via luminous accretion during quasar phases in the distant past. The fraction of inflowing rest mass energy emitted as light, the radiative efficiency, has been inferred to be 10\%, in agreement with expectations from thin disk accretion models. But the existence of billion solar-mass SMBHs powering quasars at $z > 7$ challenges this picture: provided they respect the Eddington limit, there is not enough time to grow $z>7$ SMBHs from stellar remnant seeds unless the radiative efficiency is below 10\%. Here we show that one can constrain the radiative efficiencies of the most distant quasars known using foreground neutral intergalactic gas as a cosmological-scale ionizing photon counter. From the Ly$\alpha$ absorption profiles of ULAS J1120+0641 ($z=7.09$) and ULAS J1342+0928 ($z=7.54$), we determine posterior median radiative efficiencies of 0.08\% and 0.1\%, respectively, and the combination of the two measurements rule out the canonical 10\% efficiency at 99.8\% credibility after marginalizing over the unknown obscured fraction. This low radiative efficiency implies rapid mass accretion for the earliest SMBHs, greatly easing the tension between the age of the Universe and the SMBH masses. However, our measured efficiency may instead reflect nearly complete obscuration by dusty gas in the quasar host galaxies over the vast majority of their SMBH growth. Assuming 10\% efficiency during unobscured phases, we find that the obscured fraction would be $>82\%$ at 95\% credibility, and imply a $25.7^{+49.6}_{-16.5}$ times larger obscured than unobscured luminous quasar population at $z>7$.