title: Curve counting in genus one: elliptic singularities & relative geometry

abstract: We construct and study the reduced, relative, genus one Gromov--Witten theory of very ample pairs. These invariants form the principal component contribution to relative Gromov--Witten theory in genus one and are relative versions of Zinger's reduced Gromov--Witten invariants. We relate the relative and absolute theories by degeneration of the tangency conditions, and the resulting formulas generalise a well-known recursive calculation scheme put forward by Gathmann in genus zero. The geometric input is a desingularisation of the principal component of the moduli space of genus one logarithmic stable maps to a very ample pair, using the geometry of elliptic singularities. Our study passes through general techniques for calculating integrals on logarithmic blowups of moduli spaces of stable maps, which may be of independent interest.


title: Angular separability of data clusters or network communities in geometrical space and its relevance to hyperbolic embedding

abstract: Analysis of 'big data' characterized by high-dimensionality such as word vectors and complex networks requires often their representation in a geometrical space by embedding. Recent developments in machine learning and network geometry have pointed out the hyperbolic space as a useful framework for the representation of this data derived by real complex physical systems. In the hyperbolic space, the radial coordinate of the nodes characterizes their hierarchy, whereas the angular distance between them represents their similarity. Several studies have highlighted the relationship between the angular coordinates of the nodes embedded in the hyperbolic space and the community metadata available. However, such analyses have been often limited to a visual or qualitative assessment. Here, we introduce the angular separation index (ASI), to quantitatively evaluate the separation of node network communities or data clusters over the angular coordinates of a geometrical space. ASI is particularly useful in the hyperbolic space - where it is extensively tested along this study - but can be used in general for any assessment of angular separation regardless of the adopted geometry. ASI is proposed together with an exact test statistic based on a uniformly random null model to assess the statistical significance of the separation. We show that ASI allows to discover two significant phenomena in network geometry. The first is that the increase of temperature in 2D hyperbolic network generative models, not only reduces the network clustering but also induces a 'dimensionality jump' of the network to dimensions higher than two. The second is that ASI can be successfully applied to detect the intrinsic dimensionality of network structures that grow in a hidden geometrical space.


title: Fine-tuning of orbital angular momentum in an optical parametric oscillator

abstract: We investigate the dynamical properties of a type-II optical parametric oscillator pumped by a structured light beam carrying orbital angular momentum. Different dynamical regimes are theoretically derived and experimentally demonstrated, depending on the anisotropy effects imposed by the nonlinear medium. In a weakly anisotropic regime, fine-tuning of the orbital angular momentum transfer is obtained by adjusting the crystal parameters. Under strong anisotropy, orbital angular momentum transfer is not possible and a sharp dynamical behavior is observed, with abrupt switching between different transverse modes.


title: Coping with Junk Radiation in Binary Black Hole Simulations

abstract: Spurious junk radiation in the initial data for binary black hole numerical simulations has been an issue of concern. The radiation affects the masses and spins of the black holes, modifying their orbital dynamics and thus potentially compromising the accuracy of templates used in gravitational wave analysis. Our study finds that junk radiation effects are localized to the vicinity of the black holes. Using insights from single black hole simulations, we obtain fitting formulas to estimate the changes from junk radiation on the mass and spin magnitude of the black holes in binary systems. We demonstrate how these fitting formulas could be used to adjust the initial masses and spin magnitudes of the black holes, so the resulting binary has the desired parameters after the junk radiation has left the computational domain. A comparison of waveforms from raw simulations with those from simulations that have been adjusted for junk radiation demonstrate that junk radiation could have an appreciable effect on the templates for LIGO sources with SNRs above 30.


title: Classification of glomerular hypercellularity using convolutional features and support vector machine

abstract: Glomeruli are histological structures of the kidney cortex formed by interwoven blood capillaries, and are responsible for blood filtration. Glomerular lesions impair kidney filtration capability, leading to protein loss and metabolic waste retention. An example of lesion is the glomerular hypercellularity, which is characterized by an increase in the number of cell nuclei in different areas of the glomeruli. Glomerular hypercellularity is a frequent lesion present in different kidney diseases. Automatic detection of glomerular hypercellularity would accelerate the screening of scanned histological slides for the lesion, enhancing clinical diagnosis. Having this in mind, we propose a new approach for classification of hypercellularity in human kidney images. Our proposed method introduces a novel architecture of a convolutional neural network (CNN) along with a support vector machine, achieving near perfect average results with the FIOCRUZ data set in a binary classification (lesion or normal). Our deep-based classifier outperformed the state-of-the-art results on the same data set. Additionally, classification of hypercellularity sub-lesions was also performed, considering mesangial, endocapilar and both lesions; in this multi-classification task, our proposed method just failed in 4\% of the cases. To the best of our knowledge, this is the first study on deep learning over a data set of glomerular hypercellularity images of human kidney.


title: Spontaneous critical layer formation and robustness beneath rotational waves

abstract: Non-stationary rotational surface waves are considered, where the underlying current has constant vorticity. A study is presented on the robustness of a critical layer in the presence of a bottom topography, as well as on its spontaneous formation for waves generated from rest. The restriction, from previous studies, to a traveling-wave formulation is removed leading to a non-stationary set of equations. In this setting streamlines are not necessarily pathlines. Particle-trajectories are found evolving the respective submarine dynamical system with a cloud of tracers. Pathlines are then visualized and the respective submarine structures identified. Robustness is illustrated through surface waves interacting with topographic undulations. The respective Kelvin cat eye structure dynamically adjusts itself over the bottom topography without loosing its integrity. On the spontaneous formation of a Kelvin cat eye structure, the surface is initially undisturbed and waves are generated from either the current-topography interaction or by a surface pressure distribution suddenly imposed. Under the pressure forcing, an isolated Kelvin cat eye spontaneously forms. The two extreme critical points of the cat eye structure are connected with a stagnation segment, thus exhibiting a new form of a critical layer.


title: Empirical Study of the Benefits of Overparameterization in Learning Latent Variable Models

abstract: One of the most surprising and exciting discoveries in supervised learning was the benefit of overparameterization (i.e. training a very large model) to improving the optimization landscape of a problem, with minimal effect on statistical performance (i.e. generalization). In contrast, unsupervised settings have been under-explored, despite the fact that it was observed that overparameterization can be helpful as early as Dasgupta & Schulman (2007). We perform an empirical study of different aspects of overparameterization in unsupervised learning of latent variable models via synthetic and semi-synthetic experiments. We discuss benefits to different metrics of success (recovering the parameters of the ground-truth model, held-out log-likelihood), sensitivity to variations of the training algorithm, and behavior as the amount of overparameterization increases. We find that across a variety of models (noisy-OR networks, sparse coding, probabilistic context-free grammars) and training algorithms (variational inference, alternating minimization, expectation-maximization), overparameterization can significantly increase the number of ground truth latent variables recovered.


title: The Thermodynamic Variational Objective

abstract: We introduce the thermodynamic variational objective (TVO) for learning in both continuous and discrete deep generative models. The TVO arises from a key connection between variational inference and thermodynamic integration that results in a tighter lower bound to the log marginal likelihood than the standard variational variational evidence lower bound (ELBO) while remaining as broadly applicable. We provide a computationally efficient gradient estimator for the TVO that applies to continuous, discrete, and non-reparameterizable distributions and show that the objective functions used in variational inference, variational autoencoders, wake sleep, and inference compilation are all special cases of the TVO. We use the TVO to learn both discrete and continuous deep generative models and empirically demonstrate state of the art model and inference network learning.


title: Cross-product Penalized Component Analysis (XCAN)

abstract: Matrix factorization methods are extensively employed to understand complex data. In this paper, we introduce the cross-product penalized component analysis (XCAN), a sparse matrix factorization based on the optimization of a loss function that allows a trade-off between variance maximization and structural preservation. The approach is based on previous developments, notably (i) the Sparse Principal Component Analysis (SPCA) framework based on the LASSO, (ii) extensions of SPCA to constrain both modes of the factorization, like co-clustering or the Penalized Matrix Decomposition (PMD), and (iii) the Group-wise Principal Component Analysis (GPCA) method. The result is a flexible modeling approach that can be used for data exploration in a large variety of problems. We demonstrate its use with applications from different disciplines.


title: Algorithms for weighted independent transversals and strong colouring

abstract: An independent transversal (IT) in a graph with a given vertex partition is an independent set consisting of one vertex in each partition class. Several sufficient conditions are known for the existence of an IT in a given graph with a given vertex partition, which have been used over the years to solve many combinatorial problems. Some of these IT existence theorems have algorithmic proofs, but there remains a gap between the best bounds given by nonconstructive results, and those obtainable by efficient algorithms. Recently, Graf and Haxell (2018) described a new (deterministic) algorithm that asymptotically closes this gap, but there are limitations on its applicability. In this paper we develop a randomized version of this algorithm that is much more widely applicable, and demonstrate its use by giving efficient algorithms for two problems concerning the strong chromatic number of graphs.


title: Observing cosmological binary mergers with next generation neutrino and gravitational wave detectors

abstract: We discuss the potential of detecting thermal neutrinos from matter-rich binary mergers, via a decades-long multi-messenger campaign involving a Mt-scale water Cherenkov neutrino detector and one or more next generation gravitational wave detectors, capable of observing mergers up to redshift $z\sim 2$. The search of neutrinos in time-coincidence with gravitational wave detections will allow to identify single neutrinos from individual mergers above the background, and to study their distributions in energy, redshift and type (double neutron-star or neutron-star-black hole merger) of the candidate sources. We find that, for merger rates consistent with current LIGO-Virgo constraints, and for a $100~{\rm Mt\cdot yr}$ exposure, between ${\mathcal O(10^{-1})}$ and ${\mathcal O(10)}$ neutrino events are expected. For extreme cases of mergers with more than $10^{52}$ ergs emitted in $\mathrel{{\bar \nu}_e}$, the number of events can be as large as $\sim 100$, with sensitivity to mergers up to redshift $z\sim 0.5$ or so. Such scenarios can already be tested with a $10~{\rm Mt\cdot yr}$ exposure, resulting in constraints on the post-merger evolution of the systems being considered.


title: Densest versus jammed packings of bent-core trimers

abstract: We identify putatively maximally dense packings of tangent-sphere trimers with fixed bond angles ($\theta = \theta_0$) using a novel method, and contrast them to the disordered jammed states they form under quasistatic and dynamic athermal compression. Incommensurability of $\theta_0$ with 3D close-packing does not by itself inhibit formation of dense 3D crystals; all $\theta_0$ allow formation of crystals with $\phi_{max}(\theta_0) > 0.97\phi_{cp}$. Trimers are always able to arrange into periodic structures composed of close-packed bilayers or trilayers of triangular-lattice planes, separated by ``gap layers'' that accomodate the incommensurability. All systems have $\phi_J$ significantly below the monomeric value, indicating that trimers' quenched bond-length and bond-angle constraints always act to promote jamming. $\phi_J$ varies strongly with $\theta_0$; straight ($\theta_0 = 0$) trimers minimize $\phi_J$ while closed ($\theta_0 = 120^\circ$) trimers maximize it. Marginally jammed states of trimers with lower $\phi_J(\theta_0)$ exhibit quantifiably greater disorder, and the lower $\phi_J$ for small $\theta_0$ is apparently caused by trimers' decreasing effective configurational freedom as they approach linearity.


title: Novel Suboptimal approaches for Hyperparameter Tuning of Deep Neural Network [under the shelf of Optical Communication]

abstract: Hyperparameter tuning is the main challenge of machine learning (ML) algorithms. Grid search is a popular method in hyperparameter tuning of simple ML algorithms; however, high computational complexity in complex ML algorithms such as Deep Neural Networks (DNN) is the main barrier towards its practical implementation. In this paper, two novel suboptimal grid search methods are presented, which search the grid marginally and alternating. In order to examine these methods, hyperparameter tuning is applied on two different DNN based Optical Communication (OC) systems (Fiber OC, and Free Space Optical (FSO) communication). The hyperparameter tuning of ML algorithms, despite its importance is ignored in ML for OC investigations. In addition, this is the first consideration of both FSO and Fiber OC systems in an ML for OC investigation. Results indicate that despite greatly reducing computation load, favorable performance could be achieved by the proposed methods. In addition, it is shown that the alternating search method has better performance than marginal grid search method. In sum, the proposed structures are cost-effective, and appropriate for real-time applications.


title: 3D Channel Modeling and Characterization for Hypersurface Empowered Indoor Environment at 60 GHz Millimeter-Wave Band

abstract: This paper proposes a three-dimensional (3D) communication channel model for an indoor environment considering the effect of the Hypersurface. The Hypersurface is a software controlled intelligent metasurface, which can be used to manipulate electromagnetic waves, as for example for non-specular reflection and full absorption. Thus it can control the impinging rays from a transmitter towards a receiver location in both LOS and NLOS paths, e.g. to combat distance and improve wireless connectivity. We focus on the 60 GHz mmWave frequency band due to its increasing significance in 5G/6G networks and evaluate the effect of Hypersurface in an indoor environment in terms of attenuation coefficients related to the Hypersurface reflection and absorption functionalities, using CST simulation, a 3D electromagnetic simulator of high frequency components. To highlight the benefits of Hypersurface coated walls versus plain walls, we use the derived Hypersurface 3D channel model and a custom 3D ray-tracing simulator for plain walls considering a typical indoor scenario for different Tx-Rx location and separation distances.


title: The Practical Challenges of Active Learning: Lessons Learned from Live Experimentation

abstract: We tested in a live setting the use of active learning for selecting text sentences for human annotations used in training a Thai segmentation machine learning model. In our study, two concurrent annotated samples were constructed, one through random sampling of sentences from a text corpus, and the other through model-based scoring and ranking of sentences from the same corpus. In the course of the experiment, we observed the effect of significant changes to the learning environment which are likely to occur in real-world learning tasks. We describe how our active learning strategy interacted with these events and discuss other practical challenges encountered in using active learning in the live setting.


title: The Branching-Course MPC Algorithm for Maritime Collision Avoidance

abstract: This article presents a new algorithm for short-term maritime collision avoidance (COLAV) named the branching-course MPC (BC-MPC) algorithm. The algorithm is designed to be robust with respect to noise on obstacle estimates, which is a significant source of disturbance when using exteroceptive sensors such as e.g. radars for obstacle detection and tracking. Exteroceptive sensors do not require vessel-to-vessel communication, which enables COLAV toward vessels not equipped with e.g. automatic identification system (AIS) transponders, in addition to increasing the robustness with respect to faulty information which may be provided by other vessels. The BC-MPC algorithm is compliant with rules 8 and 17 of the International Regulations for Preventing Collisions at Sea (COLREGs), and favors maneuvers following rules 13-15. This results in a COLREGs-aware algorithm which can ignore rules 13-15 when necessary. The algorithm is experimentally validated in several full-scale experiments in the Trondheimsfjord in 2017 using a radar-based system for obstacle detection and tracking. The COLAV experiments show good performance in compliance with the desired algorithm behavior.


title: Cavity dark mode of distant coupled atom-cavity systems

abstract: We report on a combined experimental and theoretical investigation into the normal modes of an all-fiber coupled cavity-quantum-electrodynamics system. The interaction between atomic ensembles and photons in the same cavities, and that between the photons in these cavities and the photons in the fiber connecting these cavities, generates five non-degenerate normal modes. We demonstrate our ability to excite each normal mode individually. We study particularly the `cavity dark mode', in which the two cavities coupled directly to the atoms do not exhibit photonic excitation. Through the observation of this mode, we demonstrate remote excitation and nonlocal saturation of atoms.


title: Exploring the structural , electronic and magnetic properties of cation ordered 3d-5d double perovskite Bi$_2$FeReO$_6$ and Bi$_2$FeIrO$_6$ thin-films from first-principles

abstract: We report a first-principles study of Bi-based 3$d$-5$d$ ordered double perovskite oxides (A$_2$BB$^\prime$O$_6$) with a 3$d$ atom (Fe) at the B-site and 5$d$ atoms (Re,Ir) at the B$^\prime$-site while keeping highly polarizable ions (Bi$^{3+}$) at the A-site. We find that, under coherent heteroepitaxy, Bi$_2$FeReO$_6$} exhibits a strain-driven anti-ferromagnetic insulator to ferrimagnetic semi-metal transition, while Bi$_2$FeIrO$_6$ shows correlation driven ferromagnetic insulator to ferrimagnetic half-metal transition with calculated magnetic moments of 5 $\mu_B$/f.u. and 3 $\mu_B$/f.u., respectively. These properties along with the low band gaps in the insulating phases make the compounds appealing for spintronics applications. Furthermore, in Bi$_2$FeIrO$_6$, the conduction and valence states are localized on different transition metal sublattices implying more efficient electron-hole separation upon photoexcitation, a desirable feature for photovoltaic applications.


title: Rhythm Dungeon: A Blockchain-based Music Roguelike Game

abstract: Rhythm Dungeon is a rhythm game which leverages the blockchain as a shared open database. During the gaming session, the player explores a roguelike dungeon by inputting specific sequences in time to music rhythm. By integrating smart contract to the game program, the enemies through the venture are generated from other games which share the identical blockchain. On the other hand, the player may upload their characters at the end of their journey, so that their own character may appear in other games and make an influence. Rhythm Dungeon is designed and implemented to show the potential of decentralized gaming experience, which utilizes the blockchain to provide asynchronous interactions among massive players.


title: The Reconstructed Power Spectrum in the Zeldovich Approximation

abstract: Density-reconstruction sharpens the baryon acoustic oscillations signal by undoing some of the smoothing incurred by nonlinear structure formation. In this paper we present an analytical model for reconstruction based on the Zeldovich approximation, which for the first time includes a complete set of counterterms and bias terms up to quadratic order and can fit real and redshift-space data pre- and post-reconstruction data in both Fourier and configuration space over a wide range of scales. We compare our model to n-body data at $z = 0$ from the {\tt DarkSky} simulation \cite{skillman14}, finding sub-percent agreement in both real space and in the redshift-space power spectrum monopole out to $k = 0.4\ h$ Mpc$^{-1}$, and out to $k = 0.2\ h$ Mpc$^{-1}$ in the quadrupole, with comparable agreement in configuration space. We compare our model with several popular existing alternatives, updating existing theoretical results for exponential damping in wiggle/no-wiggle splits of the BAO signal and discuss the usually-ignored effect of higher bias contributions on the reconstructed signal. In the appendices, we re-derive the former within our formalism, present exploratory results on higher-order corrections due to nonlinearities inherent to reconstruction, and present numerical techniques with which to calculate the redshift-space power spectrum of biased tracers within the Zeldovich approximation.


title: Free-Fermion entanglement and orthogonal polynomials

abstract: We present a simple construction for a tridiagonal matrix $T$ that commutes with the hopping matrix for the entanglement Hamiltonian ${\cal H}$ of open finite free-Fermion chains associated with families of discrete orthogonal polynomials. It is based on the notion of algebraic Heun operator attached to bispectral problems, and the parallel between entanglement studies and the theory of time and band limiting. As examples, we consider Fermionic chains related to the Chebychev, Krawtchouk and dual Hahn polynomials. For the former case, which corresponds to a homogeneous chain, the outcome of our construction coincides with a recent result of Eisler and Peschel; the latter cases yield commuting operators for particular inhomogeneous chains. Since $T$ is tridiagonal and non-degenerate, it can be readily diagonalized numerically, which in turn can be used to calculate the spectrum of ${\cal H}$, and therefore the entanglement entropy.


title: Pulse-Enhanced Two-Photon Interference with Solid State Quantum Emitters

abstract: The ability to entangle distant quantum nodes is essential for the construction of quantum networks and for quantum information processing. For solid-state quantum emitters used as qubits, it can be achieved by photon interference. When the emitter is subject to spectral diffusion, this process can become highly inefficient, impeding the achievement of scalable quantum technologies. We study two-photon interference in the context of a Hong-Ou-Mandel (HOM)-type experiment for two separate quantum emitters, with different detunings with respect to a specific target frequency. We evaluate the second order coherences that characterize photon indistinguishability between the two emitters. We find that the two-photon interference operation that is inefficient in the absence of a control protocol, when the two detunings are different and spectral overlap is lessened, can be highly improved by a periodic sequence of $\pi$ pulses at a set target frequency. Photon indistinguishability in solid state and other quantum emitters subject to spectral diffusion can thus be enhanced by the proposed pulse sequence and similar external control protocols despite the fluctuations in the environment.


title: Black holes and other exact spherical solutions in Quadratic Gravity

abstract: We study static, spherically symmetric vacuum solutions to Quadratic Gravity, extending considerably our previous Rapid Communication [Phys. Rev. D 98, 021502(R) (2018)] on this topic. Using a conformal-to-Kundt metric ansatz, we arrive at a much simpler form of the field equations in comparison with their expression in the standard spherically symmetric coordinates. We present details of the derivation of this compact form of two ordinary differential field equations for two metric functions. Next, we apply analytical methods and express their solutions as infinite power series expansions. We systematically derive all possible cases admitted by such an ansatz, arriving at six main classes of solutions, and provide recurrent formulas for all the series coefficients. These results allow us to identify the classes containing the Schwarzschild black hole as a special case. It turns out that one class contains only the Schwarzschild black hole, three classes admit the Schwarzschild solution as a special subcase, and two classes are not compatible with the Schwarzschild solution at all since they have strictly nonzero Bach tensor. In our analysis, we naturally focus on the classes containing the Schwarzschild spacetime, in particular on a new family of the Schwarzschild-Bach black holes which possesses one additional non-Schwarzschild parameter corresponding to the value of the Bach tensor invariant on the horizon. We study its geometrical and physical properties, such as basic thermodynamical quantities and tidal effects on free test particles induced by the presence of the Bach tensor. We also compare our results with previous findings in the literature obtained using the standard spherically symmetric coordinates.


title: Deformation and sorting of capsules in a T-junction

abstract: We study experimentally the motion and deformation of individual capsules transported by a constant volume-flux flow of low Reynolds number, through the T-junction of a channel with rectangular cross-section. We use millimetric ovalbumin-alginate capsules which we manufacture and characterise independently of the flow experiment. Centred capsules travel at constant velocity down the straight channel leading to the T-junction where they decelerate and expand in the spanwise direction before turning into one of the two identical daughter channels. There, non-inertial lift forces act to re-centre them and relax their shape until they reach a steady state of propagation. We find that the dynamics of fixed-size capsules within our channel geometry are governed by a capillary number Ca defined as the ratio of viscous shear forces to elastic restoring forces. We quantify the elastic forces by statically compressing the capsule to 50% of its initial diameter between parallel plates rather than by the Young's modulus of the encapsulating membrane, in order to account for different membrane thickness, pre-inflation and non-linear elastic deformation. We show that the maximum extension in the T-junction of capsules of different stiffness collapses onto a master curve in Ca. Thus, it provides a sensitive measure of the relative stiffness of capsules at constant flow rate, particularly for softer capsules. We also find that the T-junction can sort fixed-size capsules according to their stiffness because the position in the T-junction from which capsules are entrained into the daughter channel depends uniquely on Ca. We demonstrate that a T-junction can be used as a sorting device by enhancing this initial capsule separation through a diffuser.


title: Bridging the Architecture Gap: Abstracting Performance-Relevant Properties of Modern Server Processors

abstract: We describe a universal modeling approach for predicting single- and multicore runtime of steady-state loops on server processors. To this end we strictly differentiate between application and machine models: An application model comprises the loop code, problem sizes, and other runtime parameters, while a machine model is an abstraction of all performance-relevant properties of a CPU. We introduce a generic method for determining machine models and present results for relevant server-processor architectures by Intel, AMD, IBM, and Marvell/Cavium. Considering this wide range of architectures, the set of features required for adequate performance modeling is surprisingly small. To validate our approach, we compare performance predictions to empirical data for an OpenMP-parallel preconditioned CG algorithm, which includes compute- and memory-bound kernels. Both single- and multicore analysis shows that the model exhibits average and maximum relative errors of 5% and 10%. Deviations from the model and insights gained are discussed in detail.


title: Statistical properties of thermally expandable particles in soft Rayleigh-Benard convection

abstract: The dynamics of inertial particles in Rayleigh-B\'{e}nard convection, where both particles and fluid exhibit thermal expansion, is studied using direct numerical simulations (DNS). We consider the effect of particles with a thermal expansion coefficient larger than that of the fluid, causing particles to become lighter than the fluid near the hot bottom plate and heavier than the fluid near the cold top plate. Because of the opposite directions of the net Archimedes' force on particles and fluid, particles deposited at the plate now experience a relative force towards the bulk. The characteristic time for this motion towards the bulk to happen, quantified as the time particles spend inside the thermal boundary layers (BLs) at the plates, is shown to depend on the thermal response time, $\tau_T$, and the thermal expansion coefficient of particles relative to that of the fluid, $K = \alpha_p / \alpha_f$. In particular, the residence time is constant for small thermal response times, $\tau_T \lesssim 1$, and increasing with $\tau_T$ for larger thermal response times, $\tau_T \gtrsim 1$. Also, the thermal BL residence time is increasing with decreasing $K$. A one-dimensional (1D) model is developed, where particles experience thermal inertia and their motion is purely dependent on the buoyancy force. Although the values do not match one-to-one, this highly simplified 1D model does predict a regime of a constant thermal BL residence time for smaller thermal response times and a regime of increasing residence time with $\tau_T$ for larger response times, thus explaining the trends in the DNS data well.


title: State-of-the-Art on Query & Transaction Processing Acceleration

abstract: The vast amount of processing power and memory bandwidth provided by modern Graphics Processing Units (GPUs) make them a platform for data-intensive applications. The database community identified GPUs as effective co-processors for data processing. In the past years, there were many approaches to make use of GPUs at different levels of a database system. In this Internal Technical Report, based on the [1] and some other research papers, we identify possible research areas at LIP6 for GPU-accelerated database management systems. We describe some key properties, typical challenges of GPU-aware database architectures, and identify major open challenges.


title: Onset of rebound suppression in non Newtonian droplets post impact on superhydrophobic surfaces

abstract: Droplet deposition after impact on superhydrophobic surfaces has been an important area of study in recent years due to its potential application in reduction of pesticides usage. Minute amounts of long chain polymers added to water has been known to arrest the droplet rebound effect on superhydrophobic surfaces. Previous studies have attributed different reasons like extensional viscosity, dominance of elastic stresses or slowing down of contact line in retraction phase due to stretching of polymer chains. The present study attempts to unravel the existence of critical criteria of polymer concentration and impact velocity on the inhibition of droplet rebound. The impact velocity will indirectly influence the shear rate during the retraction phase, and the polymer concentration dictates the relaxation timescale of the elastic fluids. Finally we show that the Weissenberg number (at onset of retraction), which quantifies both the elastic effects of polymer chains and the hydrodynamics, is the critical parameter in determining the regime of onset of rebound suppression, and that there exists a critical value which determines the onset of bounce arrest. The previous three causes, which are manifestations of elastic effects in non-Newtonian fluids, can be related with the proposed Weissenberg number criterion.


title: Intrinsic dissipation mechanisms in metallic glass resonators

abstract: Micro- and nano-resonators have important applications including sensing, navigation, and biochemical detection. Their performance is quantified using the quality factor $Q$, which gives the ratio of the energy stored to the energy dissipated per cycle. Metallic glasses are a promising materials class for micro- and nano-scale resonators since they are amorphous and can be fabricated precisely into complex shapes on these lengthscales. To understand the intrinsic dissipation mechanisms that ultimately limit large $Q$-values in metallic glasses, we perform molecular dynamics simulations to model metallic glass resonators subjected to bending vibrations. We calculate the vibrational density of states, redistribution of energy from the fundamental mode of vibration, and $Q$ versus the kinetic energy per atom $K$ of the excitation. In the linear and nonlinear response regimes where there are no atomic rearrangements, we find that $Q \rightarrow \infty$ (since we do not consider coupling to the environment). We identify a characteristic $K_r$ above which atomic rearrangements occur, and there is significant energy leakage from the fundamental mode to higher frequencies, causing finite $Q$. Thus, $K_r$ is a critical parameter determining resonator performance. We show that $K_r$ decreases as a power-law, $K_r\sim N^{-k},$ with increasing system size $N$, where $k \approx 1.3$. We estimate the critical strain $\langle \gamma_r \rangle \sim 10^{-8}$ for micron-sized resonators below which atomic rearrangements do not occur, and thus large $Q$-values can be obtained when they are operated below $\gamma_r$. We find that $K_r$ for amorphous resonators is comparable to that for resonators with crystalline order.


title: Composable Rate-Independent Computation in Continuous Chemical Reaction Networks

abstract: Biological regulatory networks depend upon chemical interactions to process information. Engineering such molecular computing systems is a major challenge for synthetic biology and related fields. The chemical reaction network (CRN) model idealizes chemical interactions, allowing rigorous reasoning about the computational power of chemical kinetics. Here we focus on function computation with CRNs, where we think of the initial concentrations of some species as the input and the equilibrium concentration of another species as the output. Specifically, we are concerned with CRNs that are rate-independent (the computation must be correct independent of the reaction rate law) and composable ($f\circ g$ can be computed by concatenating the CRNs computing $f$ and $g$). Rate independence and composability are important engineering desiderata, permitting implementations that violate mass-action kinetics, or even "well-mixedness", and allowing the systematic construction of complex computation via modular design. We show that to construct composable rate-independent CRNs, it is necessary and sufficient to ensure that the output species of a module is not a reactant in any reaction within the module. We then exactly characterize the functions computable by such CRNs as superadditive, positive-continuous, and piecewise rational linear. Thus composability severely limits rate-independent computation unless more sophisticated input/output encodings are used.


title: Topological and noninertial effects in an Aharonov-Bohm ring

abstract: In this paper, we study the influence of topological and noninertial effects on a Dirac particle confined in an Aharonov-Bohm (AB) ring. Next, we explicitly determine the Dirac spinor and the energy spectrum for the relativistic bound states. We observe that this spectrum depends on the quantum number $n$, magnetic flux $\Phi$ of the ring, angular velocity $\omega$ associated to the noninertial effects of a rotating frame, and on the deficit angle $\eta$ associated to the topological effects of a cosmic string. We verified that this spectrum is a periodic function and grows in values as a function of $n$, $\Phi$, $\omega$, and $\eta$. In the nonrelativistic limit, we obtain the equation of motion for the particle, where now the topological effects are generated by a conic space. However, unlike relativistic case, the spectrum of this equation depends linearly on the velocity $\omega$ and decreases in values as a function of $\omega$. Comparing our results with other works, we note that our problem generalizes some particular cases of the literature. For instance, in the absence of the topological and noninertial effects ($\eta=1$ and $\omega=0$) we recover the usual spectrum of a particle confined in an AB ring ($\Phi\neq{0}$) or in an 1D quantum ring ($\Phi=0$).


title: A Survey for New Members of Taurus from Stellar to Planetary Masses

abstract: We present a large sample of new members of the Taurus star-forming region that extend from stellar to planetary masses. To identify candidate members at substellar masses, we have used color-magnitude diagrams and proper motions measured with several wide-field optical and infrared (IR) surveys. At stellar masses, we have considered the candidate members that were found in a recent analysis of high-precision astrometry from the Gaia mission. Using new and archival spectra, we have measured spectral types and assessed membership for these 161 candidates, 79 of which are classified as new members. Our updated census of Taurus now contains 519 known members. According to Gaia data, this census should be nearly complete for spectral types earlier than M6-M7 at A_J<1. For a large field encompassing ~72% of the known members, the census should be complete for K<15.7 at A_J<1.5, which corresponds to ~5-13 M_Jup for ages of 1-10 Myr based on theoretical evolutionary models. Our survey has doubled the number of known members at greater or equal to M9 and has uncovered the faintest known member in M_K, which should have a mass of ~3-10 M_Jup for ages of 1-10 Myr. We have used mid-IR photometry from the Spitzer Space Telescope and the Wide-field Infrared Survey Explorer to determine whether the new members exhibit excess emission that would indicate the presence of circumstellar disks. The updated disk fraction for Taurus is ~0.7 at less than or equal to M3.5 and ~0.4 at >M3.5.


title: Extending de Bruijn sequences to larger alphabets

abstract: A de Bruijn sequence of order n over a k-symbol alphabet is a circular sequence where each length-n sequence occurs exactly once. We present a way of extending de Bruijn sequences by adding a new symbol to the alphabet: the extension is performed by embedding a given de Bruijn sequence into another one of the same order, but over the alphabet with one more symbol, while ensuring that there are no long runs without the new symbol. Our solution is based on auxiliary graphs derived from the de Bruijn graph and solving a problem of maximum flow.


title: Fast and Exact Simulation of Multivariate Normal and Wishart Random Variables with Box Constraints

abstract: Models which include domain constraints occur in myriad contexts such as econometrics, genomics, and environmetrics, though simulating from constrained distributions can be computationally expensive. In particular, repeated sampling from constrained distributions is a common task in Bayesian inferential methods, where coping with these constraints can cause troublesome computational burden. Here, we introduce computationally efficient methods to make exact and independent draws from both the multivariate normal and Wishart distributions with box constraints. In both cases, these variables are sampled using a direct algorithm. By substantially reducing computing time, these new algorithms improve the feasibility of Monte Carlo-based inference for box-constrained, multivariate normal and Wishart distributions.


title: Explainable Anatomical Shape Analysis through Deep Hierarchical Generative Models

abstract: Quantification of anatomical shape changes currently relies on scalar global indexes which are largely insensitive to regional or asymmetric modifications. Accurate assessment of pathology-driven anatomical remodeling is a crucial step for the diagnosis and treatment of many conditions. Deep learning approaches have recently achieved wide success in the analysis of medical images, but they lack interpretability in the feature extraction and decision processes. In this work, we propose a new interpretable deep learning model for shape analysis. In particular, we exploit deep generative networks to model a population of anatomical segmentations through a hierarchy of conditional latent variables. At the highest level of this hierarchy, a two-dimensional latent space is simultaneously optimised to discriminate distinct clinical conditions, enabling the direct visualisation of the classification space. Moreover, the anatomical variability encoded by this discriminative latent space can be visualised in the segmentation space thanks to the generative properties of the model, making the classification task transparent. This approach yielded high accuracy in the categorisation of healthy and remodelled left ventricles when tested on unseen segmentations from our own multi-centre dataset as well as in an external validation set, and on hippocampi from healthy controls and patients with Alzheimer's disease when tested on ADNI data. More importantly, it enabled the visualisation in three-dimensions of both global and regional anatomical features which better discriminate between the conditions under exam. The proposed approach scales effectively to large populations, facilitating high-throughput analysis of normal anatomy and pathology in large-scale studies of volumetric imaging.


title: Fano-like spectral features in nanodiamond solutions for biometric applications

abstract: Fano resonance is a unique feature of interacting quantum systems, exhibiting resonance shapes distinctively different from conventional symmetric resonance curves. Recently, Fano resonances have been found in plasmonic nanoparticles, photonic crystals, and electromagnetic metamaterials. Here we report Fano-like photoluminiscence curves in nanodiamond solutions as a result of incoherent combination of two or more scattering and fluorescence processes. We argue that, analogously to Fano resonances, the steep asymmetric dispersion of the photoluminiscence profile in nanodiamond solutions, in combination with biologically-compatible spectral features characterizing nanodiamond fluorescence, can find promising biometric applications in several areas such as bio-sensors, bio-switches and bio-filters.


title: Singular Perturbation Theory for a Finite-Dimensional, Discrete-Time Chi Nonlinear System

abstract: This paper develops the singular perturbation theory for a particular discrete-time nonlinear system which models the saturating inductor buck converter using cycle-by-cycle digital control.


title: Complexity of acyclic colorings of graphs and digraphs with degree and girth constraints

abstract: We consider acyclic r-colorings in graphs and digraphs: they color the vertices in r colors, each of which induces an acyclic graph or digraph. (This includes the dichromatic number of a digraph, and the arboricity of a graph.) For any girth and sufficiently high degree, we prove the NP-completeness of acyclic r-colorings; our method also implies the known analogue for classical colorings. The proofs use high girth graphs with high arboricity and dichromatic numbers. High girth graphs and digraphs with high chromatic and dichromatic numbers have been well studied; we re-derive the results from a general result about relational systems, which also implies the similar fact about high girth and high arboricity used in the proofs. These facts concern graphs and digraphs of high girth and low degree; we contrast them by considering acyclic colorings of tournaments (which have low girth and high degree). We prove that even though acyclic two-colorability of tournaments is known to be NP-complete, random acyclically r-colorable tournaments allow recovering an acyclic r-coloring in deterministic linear time, with high probablity.


title: DIEL: Interactive Visualization Beyond the Here and Now

abstract: Interactive visualization design and research have primarily focused on local data and synchronous events. However, for more complex use cases---e.g., remote database access and streaming data sources---developers must grapple with distributed data and asynchronous events. Currently, constructing these use cases is difficult and time-consuming; developers are forced to operationally program low-level details like asynchronous database querying and reactive event handling. This approach is in stark contrast to modern methods for browser-based interactive visualization, which feature high-level declarative specifications. In response, we present DIEL, a declarative framework that supports asynchronous events over distributed data. Like many declarative visualization languages, DIEL developers need only specify what data they want, rather than procedural steps for how to assemble it; uniquely, DIEL models asynchronous events (e.g., user interactions or server responses) as streams of data that are captured in event logs. To specify the state of a user interface at any time, developers author declarative queries over the data and event logs; DIEL compiles and optimizes a corresponding dataflow graph, and synthesizes necessary low-level distributed systems details. We demonstrate DIEL's performance and expressivity through ex-ample interactive visualizations that make diverse use of remote data and coordination of asynchronous events. We further evaluate DIEL's usability using the Cognitive Dimensions of Notations framework, revealing wins such as ease of change, and compromises such as premature commitments.


title: Bayesian Nonparametric Boolean Factor Models

abstract: We build upon probabilistic models for Boolean Matrix and Boolean Tensor factorisation that have recently been shown to solve these problems with unprecedented accuracy and to enable posterior inference to scale to Billions of observation. Here, we lift the restriction of a pre-specified number of latent dimensions by introducing an Indian Buffet Process prior over factor matrices. Not only does the full factor-conditional take a computationally convenient form due to the logical dependencies in the model, but also the posterior over the number of non-zero latent dimensions is remarkably simple. It amounts to counting the number false and true negative predictions, whereas positive predictions can be ignored. This constitutes a very transparent example of sampling-based posterior inference with an IBP prior and, importantly, lets us maintain extremely efficient inference. We discuss applications to simulated data, as well as to a real world data matrix with 6 Million entries.


title: Existence and nonexistence of HOMO-LUMO excitations in Kohn-Sham density functional theory

abstract: In numerical computations of response properties of electronic systems, the standard model is Kohn-Sham density functional theory (KS-DFT). Here we investigate the mathematical status of the simplest class of excitations in KS-DFT, HOMO-LUMO excitations. We show using concentration-compactness arguments that such excitations, i.e. excited states of the Kohn-Sham Hamiltonian, exist for $Z>N$, where $Z$ is the total nuclear charge and $N$ is the number of electrons. The result applies under realistic assumptions on the exchange-correlation functional, which we verify explicitly for the widely used PZ81 and PW92 functionals. By contrast, and somewhat surprisingly, we find using a method of Glaser, Martin, Grosse, and Thirring \cite{glaser1976} that in case of the hydrogen and helium atoms, excited states do not exist in the neutral case $Z=N$ when the self-consistent KS ground state density is replaced by a realistic but easier to analyze approximation (in case of hydrogen, the true Schr\"{o}dinger ground state density). Implications for interpreting minus the HOMO eigenvalue as an approximation to the ionization potential are indicated.


title: Radiation Brightening from Virus-like Particles

abstract: Concentration quenching is a well-known challenge in many fluorescence imaging applications. Here we show that the optical emission from hundreds of chromophores confined onto the surface of a virus particle 28 nm diameter can be recovered under pulsed irradiation. We have found that, as one increases the number of chromophores tightly-bound to the virus surface, fluorescence quenching ensues at first, but when the number of chromophores per particle is nearing the maximum number of surface sites allowable, a sudden brightening of the emitted light and a shortening of the excited state lifetime are observed. This radiation brightening occurs only under short pulse excitation; steady-state excitation is characterized by conventional concentration quenching for any number of chromophores per particle. The observed suppression of fluorescence quenching is consistent with efficient, collective relaxation at room temperature. Interestingly, radiation brightening disappears when the emitters' spatial and/or dynamic heterogeneity is increased, suggesting that the template structural properties may play a role and opening a way towards novel, virus-enabled imaging vectors that have qualitatively different optical properties than state-of-the-art biophotonic agents.


title: Lectures on Factorization Homology, Infinity-Categories, and Topological Field Theories

abstract: These are notes from an informal mini-course on factorization homology, infinity-categories, and topological field theories. The target audience was imagined to be graduate students who are not homotopy theorists.


title: Passage through a resonance for a mechanical system, having time-varying parameters and possessing a single trapped mode. The principal term of the resonant solution

abstract: We consider a forced oscillation and passage through resonance for an infinite-length system, having time-varying parameters and possessing a single trapped mode. The system is a string, lying on the Winkler foundation and equipped with a discrete linear mass-spring oscillator of time-varying stiffness. We obtain the principal term of the asymptotic expansion for the resonant solution describing the motion of the inclusion (i.e., the mass-spring oscillator). The obtained result was verified by independent numerical calculations based on solution of the corresponding partial differential equation by means of the method of finite differences. The comparison demonstrates a good agreement in a neighbourhood of the instant of resonance.


title: On Reducing Negative Jacobian Determinant of the Deformation Predicted by Deep Registration Networks

abstract: Image registration is a fundamental step in medical image analysis. Ideally, the transformation that registers one image to another should be a diffeomorphism that is both invertible and smooth. Traditional methods like geodesic shooting approach the problem via differential geometry, with theoretical guarantees that the resulting transformation will be smooth and invertible. Most previous research using unsupervised deep neural networks for registration have used a local smoothness constraint (typically, a spatial variation loss) to address the smoothness issue. These networks usually produce non-invertible transformations with ``folding'' in multiple voxel locations, indicated by a negative determinant of the Jacobian matrix of the transformation. While using a loss function that specifically penalizes the folding is a straightforward solution, this usually requires carefully tuning the regularization strength, especially when there are also other losses. In this paper we address this problem from a different angle, by investigating possible training mechanisms that will help the network avoid negative Jacobians and produce smoother deformations. We contribute two independent ideas in this direction. Both ideas greatly reduce the number of folding locations in the predicted deformation, without making changes to the hyperparameters or the architecture used in the existing baseline registration network.


title: A 1d convolutional network for leaf and time series classification

abstract: In this paper, a 1d convolutional neural network is designed for classification tasks of plant leaves. This network based classifier is analyzed in two directions. In the forward direction, the proposed network can be used in two ways: a classifier and an automatic feature extractor. As a classifier, it takes the simple centroid contour distance curve as the single feature and achieves comparable performance with state-of-art methods that usually require multiple extracted features. As a feature extractor, it produces nearly linear separable features, hence can be used together with other classifiers such as support vector machines to provide better performance. The proposed network adopts simple 1d input and is generally applicable for other tasks such as classifying one dimensional time series in an end-to-end fashion without changes. Experiments on some benchmark datasets show this architecture can provide classification accuracies that are comparable or higher than some existing methods. In the backward direction, methods like gradient-weighted class activation mapping and maximum activation map of neurons in the classification layer with respect to inputs are performed to help investigate and further validate that hidden signatures helping trigger the trained classifier's specific decisions can be human interpretable. Code for the paper is available at https://github.com/dykuang/Leaf_Project.


title: Modeling Response Time Distributions with Generalized Beta Prime

abstract: We use Generalized Beta Prime distribution, also known as GB2, for fitting response time distributions. This distribution, characterized by one scale and three shape parameters, is incredibly flexible in that it can mimic behavior of many other distributions. GB2 exhibits power-law behavior at both front and tail ends and is a steady-state distribution of a simple stochastic differential equation. We apply GB2 in contrast studies between two distinct groups -- in this case children with dyslexia and a control group -- and show that it provides superior fitting. We compare aggregate response time distributions of the two groups for scale and shape differences (including several scale-independent measures of variability, such as Hoover index), which may in turn reflect on cognitive dynamics differences. In this approach, response time distribution of an individual can be considered as a random variate of that individual's group distribution.


title: A new estimator for gravitational lensing using galaxy and intensity mapping surveys

abstract: We introduce the Galaxy Intensity Mapping cross-COrrelation estimator (GIMCO), which is a new tomographic estimator for the gravitational lensing potential, based on a combination of intensity mapping (IM) and galaxy number counts. The estimator can be written schematically as IM$(z_f)\times$galaxy$(z_b)$ $-$ galaxy$(z_f)\times$IM$(z_b)$ for a pair of distinct redshifts $(z_f,z_b)$; this combination allows to greatly reduce the contamination by density-density correlations, thus isolating the lensing signal. As an estimator constructed only from cross-correlations, it is additionally less susceptible to systematic effects. We show that the new estimator strongly suppresses cosmic variance and consequently improves the signal-to-noise ratio (SNR) for the detection of lensing, especially on linear scales and intermediate redshifts. %This makes it particularly valuable for future studies of dark energy and modified gravity. For cosmic variance dominated surveys, the SNR of our estimator is a factor 30 larger than the SNR obtained from the correlation of galaxy number counts only. Shot noise and interferometer noise reduce the SNR. For the specific example of the Dark Energy Survey (DES) cross-correlated with the Hydrogen Intensity mapping and Real time Analysis eXperiment (HIRAX), the SNR is around 4, whereas for Euclid cross-correlated with HIRAX it reaches 52. This corresponds to an improvement of a factor 4-5 compared to the SNR from DES alone. For Euclid cross-correlated with HIRAX the improvement with respect to Euclid alone strongly depends on the redshift. We find that the improvement is particularly important for redshifts below 1.6, where it reaches a factor of 5. This makes our estimator especially valuable to test dark energy and modified gravity, that are expected to leave an impact at low and intermediate redshifts.


title: Polynomial Preconditioned GMRES to Reduce Communication in Parallel Computing

abstract: Polynomial preconditioning with the GMRES minimal residual polynomial has the potential to greatly reduce orthogonalization costs, making it useful for communication reduction. We implement polynomial preconditioning in the Belos package from Trilinos and show how it can be effective in both serial and parallel implementations. We further show it is a communication-avoiding technique and is a viable option to CA-GMRES for large-scale parallel computing.


title: Latent geometry and dynamics of proximity networks

abstract: Proximity networks are time-varying graphs representing the closeness among humans moving in a physical space. Their properties have been extensively studied in the past decade as they critically affect the behavior of spreading phenomena and the performance of routing algorithms. Yet, the mechanisms responsible for their observed characteristics remain elusive. Here, we show that many of the observed properties of proximity networks emerge naturally and simultaneously in a simple latent space network model, called dynamic-$\mathbb{S}^{1}$. The dynamic-$\mathbb{S}^{1}$ does not model node mobility directly, but captures the connectivity in each snapshot---each snapshot in the model is a realization of the $\mathbb{S}^{1}$ model of traditional complex networks, which is isomorphic to hyperbolic geometric graphs. By forgoing the motion component the model facilitates mathematical analysis, allowing us to prove the contact, inter-contact and weight distributions. We show that these distributions are power laws in the thermodynamic limit with exponents lying within the ranges observed in real systems. Interestingly, we find that network temperature plays a central role in network dynamics, dictating the exponents of these distributions, the time-aggregated agent degrees, and the formation of unique and recurrent components. Further, we show that paradigmatic epidemic and rumor spreading processes perform similarly in real and modeled networks. The dynamic-$\mathbb{S}^{1}$ or extensions of it may apply to other types of time-varying networks and constitute the basis of maximum likelihood estimation methods that infer the node coordinates and their evolution in the latent spaces of real systems.


title: Forensic Analysis of Third Party Location Applications in Android and iOS

abstract: Location sharing applications are becoming increasingly common. These applications allow users to share their own locations and view contacts' current locations on a map. Location applications are commonly used by friends and family members to view Global Positioning System (GPS) location of an individual, but valuable forensic evidence may exist in this data when stored locally on smartphones. This paper aims to discover forensic artifacts from two popular third-party location sharing applications on iOS and Android devices. Industry standard mobile forensic suites are utilized to discover if any locally stored data could be used to assist investigations reliant on knowing the past location of a suspect. Security issues raised regarding the artifacts found during our analysis is also discussed.


title: Programming with Timespans in Interactive Visualizations

abstract: Modern interactive visualizations are akin to distributed systems, where user interactions, background data processing, remote requests, and streaming data read and modify the interface at the same time. This concurrency is crucial to provide an interactive user experience---forbidding it can cripple responsiveness. However, it is notoriously challenging to program distributed systems, and concurrency can easily lead to ambiguous or confusing interface behaviors. In this paper, we present DIEL, a declarative programming model to help developers reason about and reconcile concurrency-related issues. Using DIEL, developers no longer need to procedurally describe how the interface should update based on different input events, but rather declaratively specify what the state of the interface should be as queries over event history. We show that resolving conflicts from concurrent processes in real-world interactive visualizations can be done in a few lines of DIEL code.


title: Equivariant Grothendieck-Riemann-Roch and localization in operational K-theory

abstract: We produce a Grothendieck transformation from bivariant operational $K$-theory to Chow, with a Riemann-Roch formula that generalizes classical Grothendieck-Verdier-Riemann-Roch. We also produce Grothendieck transformations and Riemann-Roch formulas that generalize the classical Adams-Riemann-Roch and equivariant localization theorems. As applications, we exhibit a projective toric variety $X$ whose equivariant $K$-theory of vector bundles does not surject onto its ordinary $K$-theory, and describe the operational $K$-theory of spherical varieties in terms of fixed-point data. In an appendix, Vezzosi studies operational $K$-theory of derived schemes and constructs a Grothendieck transformation from bivariant algebraic $K$-theory of relatively perfect complexes to bivariant operational $K$-theory.


title: Noncommutative unicellular LLT polynomials

abstract: It is known that unicellular LLT polynomials are related to the quasi-symmetric chromatic polynomials of certain graphs by the $(t-1)$-transform of symmetric functions. We investigate the extension of this transformation to various combinatorial Hopf algebras and prove a noncommutative version of this property.


title: High pO2 Floating Zone Crystal Growth of the Perovskite Nickelate PrNiO3

abstract: Single crystals of PrNiO3 were grown under an oxygen pressure of 295 bar using a unique high-pressure optical-image floating zone furnace. The crystals, with volume in excess of 1 mm3, were characterized structurally using single crystal and powder X-ray diffraction. Resistivity, specific heat, and magnetic susceptibility were measured, all of which evidenced an abrupt, first order metal-insulator transition (MIT) at ~130 K, in agreement with previous literature reports on polycrystalline specimens. Temperature-dependent single crystal diffraction was performed to investigate changes through the MIT. Our study demonstrates the opportunity space for high fugacity, reactive environments for single crystal growth specifically of perovskite nickelates but more generally to correlated electron oxides.


title: Injection of Deformable Capsules in a Reservoir, a Systematic Analysis

abstract: A computational study of capsule ejection from a narrow channel into a reservoir is undertaken for a combination of varying deformable capsule sizes and channel dimensions. A mass spring membrane model is coupled to an Immersed Boundary Lattice Boltzmann model solver. The aim of the present work is the description of the capsules motion, deformation and the response of the fluid due to the complex particle dynamics. The interactions between the capsules affect the local velocity field significantly and are responsible for the dynamics observed. Capsule membrane deformability is also seen to affect inter capsule interaction, and we observe that the train of three particles locally homogenizes the velocity field and the leading capsule travels faster than the other two trailing capsules. On the contrary, variations in size of the reservoir do not seem to be relevant, while the ratio of capsule diameter with respect to channel diameter plays a major role as well as the ratio of capsule diameter to inter capsule spacing. This flow set up has not been covered in the literature, and consequently we focus on describing capsule motion, membrane deformation and fluid dynamics, as a preliminary investigation in this field.


title: Vascular journey and adhesion mechanics of micro-sized carriers in narrow capillaries

abstract: In this work a Lattice Boltzmann Immersed Boundary method is used for predicting the dynamics of rigid and deformable adhesive micro carriers (1 um) navigating a capillary by the size of 10 um with 20% hematocrit. Red cells and particles are modeled as a collection of mass spring elements responding to a bending potential, an elastic potential and total enclosed area conservation constraint. Furthermore, particle surfaces are uniformly decorated with adhesive molecules (ligands) interacting with receptors disposed on the walls. Particle adhesion is modeled as a short range ligad receptor interaction and in term of formation and destruction probability functions that discriminate whether a chemical bond can be formed or destroyed. If a bond is established an attractive elastic force is activated. Particle transport and adhesion are characterized in terms of their ability to reach the capillary peripheries (margination rate) and firmly adhere the vasculature. This analysis is carried out systematically by varying particles' and cells' releasing positions and stiffness (Ca = 0 and 10-2). Moreover, three rigid and soft representative particles are transported on a finer mesh (Dx = 15 nm) and the chemical strength of their adhesive coating is varied (s = 0.5, 1.0, and 2.0) to precisely analyze the resulting adhesion mechanics. Stiffness is found to weakly influence the margination rate while significantly affect the ability of such constructs to efficiently interact with the endothelium by forming stable chemical bonds.


title: A Note on Shift Retrieval Problems

abstract: In this note, we discuss the shift retrieval problems, both classical and compressed, and provide connections between them using circulant matrices. We review the properties of circulant matrices necessary for our calculations and then show how shifts can be recovered from a single measurement.


title: Millemeter-Wave Fixed Wireless Access Using IEEE 802.11ay

abstract: IEEE 802.11ay defines new PHY and MAC specifications that enable 100 Gbps communications in the 60 GHz millimeter-wave (mmWave) band. Among the various use cases supported by IEEE 802.11ay, fixed wireless access, a cost-efficient high-performance alternative and/or complement to conventional fixed access, differentiates itself due to its unique requirements and characteristics. In this article, our goal is to identify and describe key elements incorporated into IEEE 802.11ay, including scheduling, beamforming, and link maintenance, that efficiently support fixed wireless access. IEEE 802.11ay is thus a viable and strong candidate to form the basis of future generations of standards-compliant (i.e., non-proprietary) mmWave fixed wireless access networks.


title: Extracting Novel Facts from Tables for Knowledge Graph Completion (Extended version)

abstract: We propose a new end-to-end method for extending a Knowledge Graph (KG) from tables. Existing techniques tend to interpret tables by focusing on information that is already in the KG, and therefore tend to extract many redundant facts. Our method aims to find more novel facts. We introduce a new technique for table interpretation based on a scalable graphical model using entity similarities. Our method further disambiguates cell values using KG embeddings as additional ranking method. Other distinctive features are the lack of assumptions about the underlying KG and the enabling of a fine-grained tuning of the precision/recall trade-off of extracted facts. Our experiments show that our approach has a higher recall during the interpretation process than the state-of-the-art, and is more resistant against the bias observed in extracting mostly redundant facts since it produces more novel extractions.


title: Constraint-preserving hybrid finite element methods for Maxwell's equations

abstract: Maxwell's equations describe the evolution of electromagnetic fields, together with constraints on the divergence of the magnetic and electric flux densities. These constraints correspond to fundamental physical laws: the nonexistence of magnetic monopoles and the conservation of charge, respectively. However, one or both of these constraints may be violated when one applies a finite element method to discretize in space. This is a well-known and longstanding problem in computational electromagnetics. We use domain decomposition to construct a family of primal hybrid finite element methods for Maxwell's equations, where the Lagrange multipliers are shown to correspond to a numerical trace of the magnetic field and a numerical flux of the electric flux density. Expressing the charge-conservation constraint in terms of this numerical flux, we show that both constraints are strongly preserved. As a special case, these methods include a hybridized version of N\'ed\'elec's method, implying that it preserves the constraints more strongly than previously recognized. These constraint-preserving properties are illustrated using numerical experiments in both the time domain and frequency domain. Additionally, we observe a superconvergence phenomenon, where hybrid post-processing yields an improved estimate of the magnetic field.


title: Large-scale inference with block structure

abstract: The detection of weak and rare effects in large amounts of data arises in a number of modern data analysis problems. Known results show that in this situation the potential of statistical inference is severely limited by the large-scale multiple testing that is inherent in these problems. Here we show that fundamentally more powerful statistical inference is possible when there is some structure in the signal that can be exploited, e.g. if the signal is clustered in many small blocks, as is the case in some relevant applications. We derive the detection boundary in such a situation where we allow both the number of blocks and the block length to grow polynomially with sample size. We derive these results both for the univariate and the multivariate settings as well as for the problem of detecting clusters in a network. These results recover as special cases the sparse mixture detection problem (Donoho and Jin, 2004) where there is no structure in the signal, as well as the scan problem (Chan and Walther, 2013) where the signal comprises a single interval. We develop methodology that allows optimal adaptive detection in the general setting, thus exploiting the structure if it is present without incurring a relevant penalty in the case where there is no structure. The advantage of this methodology can be considerable, as in the case of no structure the means need to increase at the rate $\sqrt{\log n}$ to ensure detection, while the presence of structure allows detection even if the means $decrease$ at a polynomial rate.


title: Mentoring Undergraduate Interdisciplinary Mathematics Research Students: Junior Faculty Experiences

abstract: To be successful, junior faculty must properly manage their time in the face of expanding responsibilities. One such responsibility is supervising undergraduate research projects. Student research projects (either single or multi-student) can be undertaken as a full-time summer experience, or as a part-time academic year commitment. With many potential undergraduate research formats, and with different types of students, junior faculty may find challenges in forming their research group, establishing a structure that promotes student productivity, picking an appropriate project, or in effectively mentoring their students. This article draws from the authors' experiences to help junior faculty navigate these complexities so that all parties reap the benefits of undergraduate research in interdisciplinary mathematical disciplines.


title: DRAT-based Bit-Vector Proofs in CVC4

abstract: Many state-of-the-art Satisfiability Modulo Theories (SMT) solvers for the theory of fixed-size bit-vectors employ an approach called bit-blasting, where a given formula is translated into a Boolean satisfiability (SAT) problem and delegated to a SAT solver. Consequently, producing bit-vector proofs in an SMT solver requires incorporating SAT proofs into its proof infrastructure. In this paper, we describe three approaches for integrating DRAT proofs generated by an off-the-shelf SAT solver into the proof infrastructure of the SMT solver CVC4 and explore their strengths and weaknesses. We implemented all three approaches using cryptominisat as the SAT back-end for its bit-blasting engine and evaluated performance in terms of proof-production and proof-checking.


title: Using Random Forest Machine Learning Algorithms in Binary Supernovae Classification

abstract: In the era of large all-sky surveys, there will be a need for rapid, automatic classifications of newly discovered transient objects. Our focus here is the classification of supernovae (SNe). We consider random forest machine learning algorithm applied to the classification of Type Ia and core-collapse supernovae (CCSNe) by full light curves. We also quantitatively show the potential of early-epoch classification using the same machine learning techniques. The algorithm uses the shape and magnitude of the light curve peak to determine the classification. This an initial study to essentially determine SN type with as few data points as possible. New all-sky surveys will discover new transients at a rapid rate and decisions will have to be made on very little data which transients to follow-up. Here we present an initial study where as few as five data points near peak, we an identify whether it is a Type Ia or CCSNe with better than 80\% accuracy, but with several biases in the data that will require addressing in future work.


title: Learning to Identify Patients at Risk of Uncontrolled Hypertension Using Electronic Health Records Data

abstract: Hypertension is a major risk factor for stroke, cardiovascular disease, and end-stage renal disease, and its prevalence is expected to rise dramatically. Effective hypertension management is thus critical. A particular priority is decreasing the incidence of uncontrolled hypertension. Early identification of patients at risk for uncontrolled hypertension would allow targeted use of personalized, proactive treatments. We develop machine learning models (logistic regression and recurrent neural networks) to stratify patients with respect to the risk of exhibiting uncontrolled hypertension within the coming three-month period. We trained and tested models using EHR data from 14,407 and 3,009 patients, respectively. The best model achieved an AUROC of 0.719, outperforming the simple, competitive baseline of relying prediction based on the last BP measure alone (0.634). Perhaps surprisingly, recurrent neural networks did not outperform a simple logistic regression for this task, suggesting that linear models should be included as strong baselines for predictive tasks using EHR


title: A computational proof of the linear Arithmetic Fundamental Lemma of GL$_4$

abstract: Let $K/F$ be an unramified quadratic extension of non-Archimedian local fields with residue character not equals to 2. We prove the linear Arithmetic Fundamental Lemma for GL$_4$ with the unit element in the spherical Hecke Algebra. In this article, all measures are normalized by its hyperspecial subgroup.


title: Quadrature by Two Expansions: Evaluating Laplace Layer Potentials using Complex Polynomial and Plane Wave Expansions

abstract: The recently developed quadrature by expansion (QBX) technique accurately evaluates the layer potentials with singular, weakly or nearly singular, or even hyper singular kernels in the integral equation reformulations of partial differential equations. The idea is to form a local complex polynomial or partial wave expansion centered at a point away from the boundary to avoid the singularity in the integrand, and then extrapolate the expansion at points near or even exactly on the boundary. In this paper, in addition to the local complex Taylor polynomial expansion, we derive new representations of the Laplace layer potentials using both the local complex polynomial and plane wave expansions. Unlike in the QBX, the local complex polynomial expansion in the new quadrature by two expansions (QB2X) method only collects the far-field contributions and its number of expansion terms can be analyzed using tools from the classical fast multipole method. The plane wave type expansion in the QB2X method better captures the layer potential features near the boundary. It is derived by applying the Fourier extension technique to the density and boundary geometry functions and then analytically utilizing the Residue Theorem for complex contour integrals. The internal connections of the layer potential with its density function and curvature on the boundary are explicitly revealed in the plane wave expansion and its error is bounded by the Fourier extension errors. We present preliminary numerical results to demonstrate the accuracy of the QB2X representations and to validate our analysis.


title: Neck-Pinching of $CP^1$-structures in the $PSL(2,C)$-character variety

abstract: We characterize a certain neck-pinching degeneration of (marked) $CP^1$- structures on a closed oriented surface S of genus at least two. Namely, we consider a path $C_t$ of $CP^1$-structures on S leaving every compact subset in the deformation space of (marked) $CP^1$-structures on S, such that its holonomy converges in the PSL(2, C)-character variety. In this setting, it is known that the complex structure $X_t$ of $C_t$ also leaves every compact subset in the Teichm\"uller space. In this paper, under the assumption that $X_t$ is pinched along a single loop m, we describe the limit of $C_t$ in terms of the developing maps, holomorphic quadratic differentials, and pleated surfaces. Moreover, we give an example of such a path $C_t$ whose limit holonomy is the trivial representation in the character variety.


title: Data integration for high-resolution, continental-scale estimation of air pollution concentrations

abstract: Air pollution constitutes the highest environmental risk factor in relation to heath. In order to provide the evidence required for health impact analyses, to inform policy and to develop potential mitigation strategies comprehensive information is required on the state of air pollution. Information on air pollution traditionally comes from ground monitoring (GM) networks but these may not be able to provide sufficient coverage and may need to be supplemented with information from other sources (e.g. chemical transport models; CTMs). However, these may only be available on grids and may not capture micro-scale features that may be important in assessing air quality in areas of high population. We develop a model that allows calibration between multiple data sources available at different levels of support by allowing the coefficients of calibration equations to vary over space and time, enabling downscaling where the data is sufficient to support it. The model is used to produce high-resolution (1km $\times$ 1km) estimates of NO$_2$ and PM$_{2.5}$ across Western Europe for 2010-2016. Concentrations of both pollutants are decreasing during this period, however there remain large populations exposed to levels exceeding the WHO Air Quality Guidelines and thus air pollution remains a serious threat to health.


title: Fusion products of twisted modules in permutation orbifolds

abstract: Let $V$ be a vertex operator algebra, $k$ a positive integer and $\sigma$ a permutation automorphism of the vertex operator algebra $V^{\otimes k}$. In this paper, we determine the fusion product of any $V^{\otimes k}$-module with any $\sigma$-twisted $V^{\otimes k}$-module.


title: APOGEE/Kepler Overlap Yields Orbital Solutions for a Variety of Eclipsing Binaries

abstract: Spectroscopic Eclipsing Binaries (SEBs) are fundamental benchmarks in stellar astrophysics and today are observed in breathtaking detail by missions like TESS, Kepler, and APOGEE. We develop a methodology for simultaneous analysis of high precision Kepler light curves and high resolution near-IR spectra from APOGEE and present orbital solutions and evolutionary histories for a subset of SEBs within this overlap. Radial velocities extracted from APOGEE spectra using the Broadening Function technique are combined with Kepler light curves and to yield binary orbital solutions. The Broadening Function approach yields more precise radial velocities than the standard Cross-Correlation Function, which in turn yields more precise orbital parameters and enables the identification of tertiary stars. The orbital periods of these seven SEBs range from 4 to 40 days. Four of the systems (KIC 5285607, KIC 6864859, KIC 6778289, and KIC 4285087) are well-detached binaries. The remaining three systems have apparent tertiary companions, but each exhibits two eclipses along with at least one spectroscopically varying component (KIC 6449358, KIC 6131659, and KIC 6781535). Gaia distances are available for four targets which we use to estimate temperatures of both members of these SEBs. We explore evolutionary histories in H-R diagram space and estimate ages for this subset of our sample. Finally, we consider the implications for the formation pathways of close binary systems via interactions with tertiary companions. Our methodology combined with the era of big data and observation overlap opens up the possibility of discovering and analyzing large numbers of diverse SEBs, including those with high flux ratios and those in triple systems.


title: Solving Polynomial Systems with phcpy

abstract: The solutions of a system of polynomials in several variables are often needed, e.g.: in the design of mechanical systems, and in phase-space analyses of nonlinear biological dynamics. Reliable, accurate, and comprehensive numerical solutions are available through PHCpack, a FOSS package for solving polynomial systems with homotopy continuation. This paper explores new developments in phcpy, a scripting interface for PHCpack, over the past five years. For instance, phcpy is now available online through a JupyterHub server featuring Python2, Python3, and SageMath kernels. As small systems are solved in real-time by phcpy, they are suitable for interactive exploration through the notebook interface. Meanwhile, phcpy supports GPU parallelization, improving the speed and quality of solutions to much larger polynomial systems. From various model design and analysis problems in STEM, certain classes of polynomial system frequently arise, to which phcpy is well-suited.


title: Parallel Performance of Molecular Dynamics Trajectory Analysis

abstract: The performance of biomolecular molecular dynamics simulations has steadily increased on modern high performance computing resources but acceleration of the analysis of the output trajectories has lagged behind so that analyzing simulations is becoming a bottleneck. To close this gap, we studied the performance of parallel trajectory analysis with MPI and the Python MDAnalysis library on three different XSEDE supercomputers where trajectories were read from a Lustre parallel file system. Strong scaling performance was impeded by stragglers, MPI processes that were slower than the typical process. Stragglers were less prevalent for compute-bound workloads, thus pointing to file reading as a bottleneck for scaling. However, a more complicated picture emerged in which both the computation and the data ingestion exhibited close to ideal strong scaling behavior whereas stragglers were primarily caused by either large MPI communication costs or long times to open the single shared trajectory file. We improved overall strong scaling performance by either subfiling (splitting the trajectory into separate files) or MPI-IO with Parallel HDF5 trajectory files. The parallel HDF5 approach resulted in near ideal strong scaling on up to 384 cores (16 nodes), thus reducing trajectory analysis times by two orders of magnitude compared to the serial approach.


title: Robustness Guarantees for Deep Neural Networks on Videos

abstract: The widespread adoption of deep learning models places demands on their robustness. In this paper, we consider the robustness of deep neural networks on videos, which comprise both the spatial features of individual frames extracted by a convolutional neural network and the temporal dynamics between adjacent frames captured by a recurrent neural network. To measure robustness, we study the maximum safe radius problem, which computes the minimum distance from the optical flow sequence obtained from a given input to that of an adversarial example in the neighbourhood of the input. We demonstrate that, under the assumption of Lipschitz continuity, the problem can be approximated using finite optimisation via discretising the optical flow space, and the approximation has provable guarantees. We then show that the finite optimisation problem can be solved by utilising a two-player turn-based game in a cooperative setting, where the first player selects the optical flows and the second player determines the dimensions to be manipulated in the chosen flow. We employ an anytime approach to solve the game, in the sense of approximating the value of the game by monotonically improving its upper and lower bounds. We exploit a gradient-based search algorithm to compute the upper bounds, and the admissible A* algorithm to update the lower bounds. Finally, we evaluate our framework on the UCF101 video dataset.


title: Weighted $\mathsf{P}-$partitions enumerator

abstract: To an extended permutohedron we associate the weighted integer points enumerator, whose principal specialization is the $f$-polynomial. In the case of poset cones it refines Gessel's $\mathsf{P}$-partitions enumerator. We show that this enumerator is a quasisymmetric function obtained by universal morphism from the Hopf algebra of posets.


title: Gravito-optics and intensity correlations for binary inspiral signal detections

abstract: We examine the correlation functions associated with intensity interferometry and gravito-optics of gravitational wave signals from compact binary coalescences. Previous theoretical studies of the gravito-optics of gravitational waves has concentrated on the characterization of both the classical and the non-classical properties of signals from cosmological sources in the early Universe. These previous works assume a periodic signal similar to the signals studied widely in optics and quantum optics and do not apply to transient signals. We develop the gravito-optics of intensity correlations for descriptions of the detection of transient signals from compact binary coalescences and apply these methods to calculate the two-point intensity correlations for the gravitational wave discovery. We also discuss the necessary theoretical work required for the description of the quantum gravito-optics of intensity correlations in the detection of signals from binary inspirals.


title: Dislocation drag from phonon wind in an isotropic crystal at large velocities

abstract: The anharmonic interaction and scattering of phonons by a moving dislocation, the photon wind, imparts a drag force $v B(v, T, \rho)$ on the dislocation. In early studies the drag coefficient $B$ was computed and experimentally determined only for dislocation velocities $v$ much less than transverse sound speed, $c_t$. In this paper we derive analytic expressions for the velocity dependence of $B$ up to $c_t$ in terms of the third-order continuum elastic constants of an isotropic crystal, in the continuum Debye approximation, valid for dislocation velocities approaching the sound speed. In so doing we point out that the most general form of the third order elastic potential for such a crystal and the dislocation-phonon interaction requires two additional elastic constants involving asymmetric local rotational strains, which have been neglected previously. We compute the velocity dependence of the transverse phonon wind contribution to $B$ in the range 1%-90% $c_t$ for Al, Cu, Fe, and Nb in the isotropic Debye approximation. The drag coefficient for transverse phonons scattering from screw dislocations is finite as $v \rightarrow c_t$, whereas $B$ is divergent for transverse phonons scattering from edge dislocations in the same limit. This divergence indicates the breakdown of the Debye approximation and sensitivity of the drag coefficient at very high velocities to the microscopic crystalline lattice cutoff. We compare our results to experimental results wherever possible and identify ways to validate and further improve the theory of dislocation drag at high velocities with realistic phonon dispersion relations, inclusion of lattice cutoff effects, MD simulation data, and more accurate experimental measurements.


title: The Complexity of Tiling Problems

abstract: In this document, we collected the most important complexity results of tilings. We also propose a definition of a so-called deterministic set of tile types, in order to capture deterministic classes without the notion of games. We also pinpoint tiling problems complete for respectively LOGSPACE and NLOGSPACE.


title: Learning Effective Loss Functions Efficiently

abstract: We consider the problem of learning a loss function which, when minimized over a training dataset, yields a model that approximately minimizes a validation error metric. Though learning an optimal loss function is NP-hard, we present an anytime algorithm that is asymptotically optimal in the worst case, and is provably efficient in an idealized "easy" case. Experimentally, we show that this algorithm can be used to tune loss function hyperparameters orders of magnitude faster than state-of-the-art alternatives. We also show that our algorithm can be used to learn novel and effective loss functions on-the-fly during training.


title: On orbit performance of the GRACE Follow-On Laser Ranging Interferometer

abstract: The Laser Ranging Interferometer (LRI) instrument on the Gravity Recovery and Climate Experiment (GRACE) Follow-On mission has provided the first laser interferometric range measurements between remote spacecraft, separated by approximately 220 km. Autonomous controls that lock the laser frequency to a cavity reference and establish the 5 degree of freedom two-way laser link between remote spacecraft succeeded on the first attempt. Active beam pointing based on differential wavefront sensing compensates spacecraft attitude fluctuations. The LRI has operated continuously without breaks in phase tracking for more than 50 days, and has shown biased range measurements similar to the primary ranging instrument based on microwaves, but with much less noise at a level of $1\,{\rm nm}/\sqrt{\rm Hz}$ at Fourier frequencies above 100 mHz.


title: Tableau Stabilization and Rectangular Tableaux Fixed by Promotion Powers

abstract: We introduce tableau stabilization, a new phenomenon and statistic on Young tableaux based on jeu de taquin. We investigate bounds for tableau stabilization, the shape of stabilized tableaux, and tableau stabilization as a permutation statistic. We apply tableau stabilization to construct the sufficiently large rectangular tableaux fixed by powers of promotion, which were counted by Brendon Rhoades via the cyclic sieving phenomenon


title: Smart Charging Benefits in Autonomous Mobility on Demand Systems

abstract: In this paper, we study the potential benefits from smart charging for a fleet of electric vehicles (EVs) providing autonomous mobility-on-demand (AMoD) services. We first consider a profit-maximizing platform operator who makes decisions for routing, charging, rebalancing, and pricing for rides based on a network flow model. Clearly, each of these decisions directly influence the fleet's smart charging potential; however, it is not possible to directly characterize the effects of various system parameters on smart charging under a classical network flow model. As such, we propose a modeling variation that allows us to decouple the charging and routing problems faced by the operator. This variation allows us to provide closed-form mathematical expressions relating the charging costs to the maximum battery capacity of the vehicles as well as the fleet operational costs. We show that investing in larger battery capacities and operating more vehicles for rebalancing reduces the charging costs, while increasing the fleet operational costs. Hence, we study the trade-off the operator faces, analyze the minimum cost fleet charging strategy, and provide numerical results illustrating the smart charging benefits to the operator.


title: Adaptive Sequential Experiments with Unknown Information Arrival Processes

abstract: Sequential experiments are often characterized by an exploration-exploitation tradeoff that is captured by the multi-armed bandit (MAB) framework. This framework has been studied and applied, typically when at each time period feedback is received only on the action that was selected at that period. However, in many practical settings additional data may become available between decision epochs. We introduce a generalized MAB formulation, which considers a broad class of distributions that are informative about mean rewards, and allows observations from these distributions to arrive according to an arbitrary and a priori unknown arrival process. When it is known how to map auxiliary data to reward estimates, by obtaining matching lower and upper bounds we characterize a spectrum of minimax complexities for this class of problems as a function of the information arrival process, which captures how salient characteristics of this process impact achievable performance. In terms of achieving optimal performance, we establish that upper confidence bound and posterior sampling policies possess natural robustness with respect to the information arrival process without any adjustments, which uncovers a novel property of these popular policies and further lends credence to their appeal. When the mappings connecting auxiliary data and rewards are a priori unknown, we characterize necessary and sufficient conditions under which auxiliary information allows performance improvement. We devise a new policy that is based on two different upper confidence bounds (one that accounts for auxiliary observation and one that does not) and establish the near-optimality of this policy. We use data from a large media site to analyze the value that may be captured in practice by leveraging auxiliary data for designing content recommendations.


title: On a Center-of-Mass System of Coordinates for Symmetric Classical and Quantum Many-Body Problems

abstract: In the context of classical or quantum many-body problems involving identical bodies, a linear change of coordinates can be constructed with the properties that it includes the center-of-mass as one of the new coordinates and preserves the inherent permutation symmetry of both the Hamiltonian and the admissible states. This has advantages over the usual system of Jacobi coordinates in the study of many-body problems for which permutation symmetry of the bodies plays an important role. This paper contains the details of the construction of this system and the proof that these properties uniquely determine it, up to trivial modifications. Examples of applications to both classical and quantum problems are explored, including a generalization to problems involving groups of different species of bodies.