title: Electronic and phonon excitations in {\alpha}-RuCl$_3$

abstract: We report on THz, infrared reflectivity and transmission experiments for wave numbers from 10 to 8000 cm$^{-1}$ ($\sim$ 1 meV - 1 eV) and for temperatures from 5 to 295 K on the Kitaev candidate material {\alpha}-RuCl$_3$. As reported earlier, the compound under investigation passes through a first-order structural phase transition, from a monoclinic high-temperature to a rhombohedral low-temperature phase. The phase transition shows an extreme and unusual hysteretic behavior, which extends from 60 to 166 K. In passing this phase transition, in the complete frequency range investigated we found a significant reflectance change, which amounts almost a factor of two. We provide a broadband spectrum of dielectric constant, dielectric loss and optical conductivity from the THz to the mid infrared regime and study in detail the phonon response and the low-lying electronic density of states. We provide evidence for the onset of an optical energy gap, which is of order 200 meV, in good agreement with the gap derived from measurements of the DC electrical resistivity. Remarkably, the onset of the gap exhibits a strong blue shift on increasing temperatures.


title: Demystifying the Characteristics of 3D-Stacked Memories: A Case Study for Hybrid Memory Cube

abstract: Three-dimensional (3D)-stacking technology, which enables the integration of DRAM and logic dies, offers high bandwidth and low energy consumption. This technology also empowers new memory designs for executing tasks not traditionally associated with memories. A practical 3D-stacked memory is Hybrid Memory Cube (HMC), which provides significant access bandwidth and low power consumption in a small area. Although several studies have taken advantage of the novel architecture of HMC, its characteristics in terms of latency and bandwidth or their correlation with temperature and power consumption have not been fully explored. This paper is the first, to the best of our knowledge, to characterize the thermal behavior of HMC in a real environment using the AC-510 accelerator and to identify temperature as a new limitation for this state-of-the-art design space. Moreover, besides bandwidth studies, we deconstruct factors that contribute to latency and reveal their sources for high- and low-load accesses. The results of this paper demonstrates essential behaviors and performance bottlenecks for future explorations of packet-switched and 3D-stacked memories.


title: Short Spacing Considerations for the ngVLA

abstract: The next generation Very Large Array project (ngVLA) would represent a major step forward in sensitivity and resolution for radio astronomy, with ability to achieve 2 milli-arcsec resolution at 100 GHz (assuming a maximum baseline of 300 km). For science on spatial scales of >~ 1 arcsec, the ngVLA project should consider the use of a large single dish telescope to provide short-spacing data. Large single-dish telescopes are complementary to interferometers and are crucial to providing sensitivity to spatial scales lost by interferometry. Assuming the current vision of the ngVLA (300 18m dishes) and by studying possible array configurations, I argue that a single dish with a diameter of >= 45m with approximately 20 element receiver systems would be well matched to the ngVLA for mapping observations.


title: Targeted Nuclear Spallation from Moderator Block Design for a Ground-Based Space Radiation Analog

abstract: Current radiobiology studies on the effects of galactic cosmic ray radiation utilize mono-energetic beams, where the projected dose for an exploration mission is given using highly-acute exposures. This methodology does not replicate the multi-ion species and energies found in the space radiation environment, nor does it reflect the low dose-rate found in interplanetary space. In radiation biology studies as well as in the assessment of health risk to astronaut crews, the differences in the biological effectiveness of different ions is primarily attributed to differences in the linear energy transfer (LET) of the radiation spectrum. Here we show that the LET spectrum of the intravehicular environment of spaceflight vehicles can be simulated with a single particle, mono-energetic ion beam accelerated at target blocks constructed of one or more materials. The LET spectrum of the emerging field can then be moderated by the amount of mass or length of material the primary and secondary nuclei travels, thus preferentially producing specific nuclear spallation and fragmentation processes and allowing for a continuous generation of ionizing radiation that mimics the space radiation environment. This approach could allow more accurate simulation of not only intravehicular spaceflight conditions, but also could be used to simulate the external galactic cosmic ray field, planetary surface spectrum (e.g., Mars or Moon), and the local radiation environment of orbiting satellites, providing a much-needed ground-based space radiation analog for future experimentation.


title: The role of pair correlation function in the dynamical transition predicted by the mode coupling theory

abstract: In a recent study we have found that for a large number of systems the configuration entropy at pair level, $S_{c2}$, which is primarily determined by the structural information, vanishes at the mode coupling transition temperature $T_{c}$. Thus it appears that the information of the transition temperature is embedded in the structure of the liquid. In order to investigate this we describe the dynamics of the system at the mean field level and using the concepts of the dynamical density function theory show that the dynamics depends only on the structure of the liquid. Thus this theory is similar in spirit to the microscopic MCT. However unlike microscopic MCT, which predicts a very high transition temperature, the present theory predicts a transition temperature which is similar to $T_{c}$. Thus our study reveals that the information of the mode coupling transition temperature is embedded in the structure of the liquid.


title: Modeling of the outburst on July 29th, 2015 observed with OSIRIS cameras in the southern hemisphere of comet 67P/Churyumov-Gerasimenko

abstract: Images of the nucleus and the coma (gas and dust) of comet 67P/Churyumov- Gerasimenko have been acquired by the OSIRIS (Optical, Spectroscopic, and Infrared Remote Imaging System) cameras since March 2014 using both the Wide Angle Camera (WAC) and the Narrow Angle Camera (NAC). We use images from the NAC camera to study a bright outburst observed in the southern hemisphere on July 29, 2015. The high spatial resolution of the NAC is needed to localize the source point of the outburst on the surface of the nucleus. The heliocentric distance is 1.25 au and the spacecraft-comet distance is 186 km. Aiming to better understand the physics that led to the outgassing, we used the Direct Simulation Monte Carlo (DSMC) method to study the gas flow close to the nucleus and the dust trajectories. The goal is to understand the mechanisms producing the outburst. We reproduce the opening angle of the outburst in the model and constrain the outgassing ratio between the outburst source and the local region. The outburst is in fact a combination of both gas and dust, in which the active surface is approximately 10 times more active than the average rate found in the surrounding areas. We need a number of dust particles 7.83 $\times$ 10$^{11}$ - 6.90 $\times$ 10$^{15}$ (radius 1.97 - 185 {\mu}m), which corresponds to a mass of dust 220 - 21 $\times$ 10$^{3}$kg.


title: Random projections for trust region subproblems

abstract: The trust region method is an algorithm traditionally used in the field of derivative free optimization. The method works by iteratively constructing surrogate models (often linear or quadratic functions) to approximate the true objective function inside some neighborhood of a current iterate. The neighborhood is called "trust region in the sense that the model is trusted to be good enough inside the neighborhood. Updated points are found by solving the corresponding trust region subproblems. In this paper, we describe an application of random projections to solving trust region subproblems approximately.


title: Capacity Comparison between MIMO-NOMA and MIMO-OMA with Multiple Users in a Cluster

abstract: In this paper, the performance of multiple-input multiple-output non-orthogonal multiple access (MIMO-NOMA) is investigated when multiple users are grouped into a cluster. The superiority of MIMO-NOMA over MIMO orthogonal multiple access (MIMO-OMA) in terms of both sum channel capacity and ergodic sum capacity is proved analytically. Furthermore, it is demonstrated that the more users are admitted to a cluster, the lower is the achieved sum rate, which illustrates the tradeoff between the sum rate and maximum number of admitted users. On this basis, a user admission scheme is proposed, which is optimal in terms of both sum rate and number of admitted users when the signal-to-interference-plus-noise ratio thresholds of the users are equal. When these thresholds are different, the proposed scheme still achieves good performance in balancing both criteria. Moreover, under certain conditions,it maximizes the number of admitted users. In addition, the complexity of the proposed scheme is linear to the number of users per cluster. Simulation results verify the superiority of MIMO-NOMA over MIMO-OMA in terms of both sum rate and user fairness, as well as the effectiveness of the proposed user admission scheme.


title: The Grad-Shafranov Reconstruction of Toroidal Magnetic Flux Ropes: Method Development and Benchmark Studies

abstract: We develop an approach of Grad-Shafranov (GS) reconstruction for toroidal structures in space plasmas, based on in-situ spacecraft measurements. The underlying theory is the GS equation that describes two-dimensional magnetohydrostatic equilibrium as widely applied in fusion plasmas. The geometry is such that the arbitrary cross section of the torus has rotational symmetry about the rotation axis $Z$, with a major radius $r_0$. The magnetic field configuration is thus determined by a scalar flux function $\Psi$ and a functional $F$ that is a single-variable function of $\Psi$. The algorithm is implemented through a two-step approach: i) a trial-and-error process by minimizing the residue of the functional $F(\Psi)$ to determine an optimal $Z$ axis orientation, and ii) for the chosen $Z$, a $\chi^2$ minimization process resulting in the range of $r_0$. Benchmark studies of known analytic solutions to the toroidal GS equation with noise additions are presented to illustrate the two-step procedures and to demonstrate the performance of the numerical GS solver, separately. For the cases presented, the errors in $Z$ and $r_0$ are 9$^\circ$ and 22\%, respectively, and the relative percent error in the numerical GS solutions is less than 10\%. We also make public the computer codes for these implementations and benchmark studies.


title: Climbing a shaky ladder: Better adaptive risk estimation

abstract: We revisit the \emph{leaderboard problem} introduced by Blum and Hardt (2015) in an effort to reduce overfitting in machine learning benchmarks. We show that a randomized version of their Ladder algorithm achieves leaderboard error O(1/n^{0.4}) compared with the previous best rate of O(1/n^{1/3}). Short of proving that our algorithm is optimal, we point out a major obstacle toward further progress. Specifically, any improvement to our upper bound would lead to asymptotic improvements in the general adaptive estimation setting as have remained elusive in recent years. This connection also directly leads to lower bounds for specific classes of algorithms. In particular, we exhibit a new attack on the leaderboard algorithm that both theoretically and empirically distinguishes between our algorithm and previous leaderboard algorithms.


title: Rectifiability of line defects in liquid crystals with variable degree of orientation

abstract: In [AHL] Hardt, Lin and the author proved that the defect set of minimizers of the modified Ericksen energy for nematic liquid crystals consists locally of a finite union of isolated points and H\"older continuous curves with finitely many crossings. In this article, we show that each H\"older continuous curve in the defect set is of finite length. Hence, the defect set is locally rectifiable. For the most part, the proof follows the work of De Lellis, Marchese, Spadaro and Valtorta [DLMSV] on harmonic $\mathcal{Q}$-valued maps closely. The blow-up analysis in [AHL] allows us to simplify the covering arguments in [DLMSV] and locally estimate the length of line defects in a geometric fashion.


title: CortexNet: a Generic Network Family for Robust Visual Temporal Representations

abstract: In the past five years we have observed the rise of incredibly well performing feed-forward neural networks trained supervisedly for vision related tasks. These models have achieved super-human performance on object recognition, localisation, and detection in still images. However, there is a need to identify the best strategy to employ these networks with temporal visual inputs and obtain a robust and stable representation of video data. Inspired by the human visual system, we propose a deep neural network family, CortexNet, which features not only bottom-up feed-forward connections, but also it models the abundant top-down feedback and lateral connections, which are present in our visual cortex. We introduce two training schemes - the unsupervised MatchNet and weakly supervised TempoNet modes - where a network learns how to correctly anticipate a subsequent frame in a video clip or the identity of its predominant subject, by learning egomotion clues and how to automatically track several objects in the current scene. Find the project website at https://engineering.purdue.edu/elab/CortexNet/.


title: Higher-order symmetry energy and neutron star core-crust transition with Gogny forces

abstract: We study the symmetry energy and the core-crust transition in neutron stars using the finite-range Gogny nuclear interaction and examine the deduced crustal thickness and crustal moment of inertia. We start by analyzing the second-, fourth- and sixth-order coefficients of the Taylor expansion of the energy per particle in powers of the isospin asymmetry for Gogny forces. These coefficients provide information about the departure of the symmetry energy from the widely used parabolic law. The neutron star core-crust transition is evaluated by looking at the onset of thermodynamical instability of the liquid core. The calculation is performed with the exact (i.e., without Taylor expansion) Gogny EoS for the core, and also with its Taylor expansion in order to assess the influence of isospin expansions on locating the inner edge of neutron star crusts. It is found that the properties of the core-crust transition derived from the exact EoS differ from the predictions of the Taylor expansion even when the expansion is carried through sixth order in the isospin asymmetry. Gogny forces, using the exact EoS, predict the ranges $0.094 \text{ fm}^{-3} \lesssim \rho_t \lesssim 0.118\text{ fm}^{-3}$ for the transition density and $0.339 \text{ MeV fm}^{-3} \lesssim P_t \lesssim 0.665 \text{ MeV fm}^{-3}$ for the transition pressure. The transition densities show an anticorrelation with the slope parameter $L$ of the symmetry energy. The transition pressures are not found to correlate with $L$. Neutron stars obtained with Gogny forces have maximum masses below $1.74M_\odot$ and relatively small moments of inertia. The crustal mass and moment of inertia are evaluated and comparisons are made with the constraints from observed glitches in pulsars.


title: Advances in Joint CTC-Attention based End-to-End Speech Recognition with a Deep CNN Encoder and RNN-LM

abstract: We present a state-of-the-art end-to-end Automatic Speech Recognition (ASR) model. We learn to listen and write characters with a joint Connectionist Temporal Classification (CTC) and attention-based encoder-decoder network. The encoder is a deep Convolutional Neural Network (CNN) based on the VGG network. The CTC network sits on top of the encoder and is jointly trained with the attention-based decoder. During the beam search process, we combine the CTC predictions, the attention-based decoder predictions and a separately trained LSTM language model. We achieve a 5-10\% error reduction compared to prior systems on spontaneous Japanese and Chinese speech, and our end-to-end model beats out traditional hybrid ASR systems.


title: Origin of ultrastability in vapor-deposited glasses

abstract: Glass films created by vapor-depositing molecules onto a substrate can exhibit properties similar to those of ordinary glasses aged for thousands of years. It is believed that enhanced surface mobility is the mechanism that allows vapor deposition to create such exceptional glasses, but it is unclear how this effect is related to the final state of the film. Here we use molecular dynamics simulations to model vapor deposition and an efficient Monte Carlo algorithm to determine the deposition rate needed to create ultra-stable glassy films. We obtain a scaling relation that quantitatively captures the efficiency gain of vapor deposition over bulk annealing, and demonstrates that surface relaxation plays the same role in the formation of vapor-deposited glasses as bulk relaxation does in ordinary glass formation.


title: Cosmology at at Crossroads: Tension with the Hubble Constant

abstract: We are at an interesting juncture in cosmology. With new methods and technology, the accuracy in measurement of the Hubble constant has vastly improved, but a recent tension has arisen that is either signaling new physics or as-yet unrecognized uncertainties.


title: A New Class of Fully Discrete Sparse Fourier Transforms: Faster Stable Implementations with Guarantees

abstract: In this paper we consider Sparse Fourier Transform (SFT) algorithms for approximately computing the best $s$-term approximation of the Discrete Fourier Transform (DFT) $\mathbf{\hat{f}} \in \mathbb{C}^N$ of any given input vector $\mathbf{f} \in \mathbb{C}^N$ in just $\left( s \log N\right)^{\mathcal{O}(1)}$-time using only a similarly small number of entries of $\mathbf{f}$. In particular, we present a deterministic SFT algorithm which is guaranteed to always recover a near best $s$-term approximation of the DFT of any given input vector $\mathbf{f} \in \mathbb{C}^N$ in $\mathcal{O} \left( s^2 \log ^{\frac{11}{2}} (N) \right)$-time. Unlike previous deterministic results of this kind, our deterministic result holds for both arbitrary vectors $\mathbf{f} \in \mathbb{C}^N$ and vector lengths $N$. In addition to these deterministic SFT results, we also develop several new publicly available randomized SFT implementations for approximately computing $\mathbf{\hat{f}}$ from $\mathbf{f}$ using the same general techniques. The best of these new implementations is shown to outperform existing discrete sparse Fourier transform methods with respect to both runtime and noise robustness for large vector lengths $N$.


title: Probing the topology of density matrices

abstract: The mixedness of a quantum state is usually seen as an adversary to topological quantization of observables. For example, exact quantization of the charge transported in a so-called Thouless adiabatic pump is lifted at any finite temperature in symmetry-protected topological insulators. Here, we show that certain directly observable many-body correlators preserve the integrity of topological invariants for mixed Gaussian quantum states in one dimension. Our approach relies on the expectation value of the many-body momentum-translation operator, and leads to a physical observable --- the "ensemble geometric phase" (EGP) --- which represents a bona fide geometric phase for mixed quantum states, in the thermodynamic limit. In cyclic protocols, the EGP provides a topologically quantized observable which detects encircled spectral singularities ("purity-gap" closing points) of density matrices. While we identify the many-body nature of the EGP as a key ingredient, we propose a conceptually simple, interferometric setup to directly measure the latter in experiments with mesoscopic ensembles of ultracold atoms.


title: Twisted smooth Deligne cohomology

abstract: Deligne cohomology can be viewed as a differential refinement of integral cohomology, hence captures both topological and geometric information. On the other hand, it can be viewed as the simplest nontrivial version of a differential cohomology theory. While more involved differential cohomology theories have been explicitly twisted, the same has not been done to Deligne cohomology, although existence is known at a general abstract level. We work out what it means to twist Deligne cohomology, by taking degree one twists of both integral cohomology and de Rham cohomology. We present the main properties of the new theory and illustrate its use with examples and applications. Given how versatile Deligne cohomology has proven to be, we believe that this explicit and utilizable treatment of its twisted version will be useful.


title: On the intersection motive of certain Shimura varieties: the case of Siegel threefolds

abstract: In this article, we construct a Hecke-equivariant Chow motive whose realizations equal intersection cohomology of Siegel threefolds with regular algebraic coefficients. As a consequence, we are able to define Grothendieck motives for Siegel modular forms.


title: Avoiding Discrimination through Causal Reasoning

abstract: Recent work on fairness in machine learning has focused on various statistical discrimination criteria and how they trade off. Most of these criteria are observational: They depend only on the joint distribution of predictor, protected attribute, features, and outcome. While convenient to work with, observational criteria have severe inherent limitations that prevent them from resolving matters of fairness conclusively. Going beyond observational criteria, we frame the problem of discrimination based on protected attributes in the language of causal reasoning. This viewpoint shifts attention from "What is the right fairness criterion?" to "What do we want to assume about the causal data generating process?" Through the lens of causality, we make several contributions. First, we crisply articulate why and when observational criteria fail, thus formalizing what was before a matter of opinion. Second, our approach exposes previously ignored subtleties and why they are fundamental to the problem. Finally, we put forward natural causal non-discrimination criteria and develop algorithms that satisfy them.


title: A break in the high-redshift stellar mass Tully-Fisher relation

abstract: We investigate the stellar-mass Tully-Fisher relation (TFR) between the stellar mass and the integrated gas velocity dispersion, quantified by the kinematic estimator S_0.5 measured from strong emission lines in spectra of galaxies at 0<z<5. We combine luminosity-selected galaxies (`high-luminosity sample') with galaxies selected in other ways (`low-luminosity sample') to cover a range in stellar mass that spans almost five orders of magnitude: 7.0 < log M* < 11.5. We find that the logarithmic power-law slope and normalisation of the TFR are independent of redshift out to z~3. The scatter in the TFR is <0.5 dex such that the gas velocity dispersion can be used as a proxy for the stellar mass of a galaxy independently of its redshift. At z>3 the scatter increases and the existence of a correlation is not obvious. The high-luminosity sample exhibits a flatter slope of 1.5$\pm$0.2 at z<3 compared to the low-luminosity sample slope of 2.9$\pm$0.3, suggesting a turnover in the TFR. The combined sample is well fit with a break in the TFR at a characteristic stellar mass scale of M*~10$^{10}$ M$_{\odot}$, with no significant evolution out to z~3. We demonstrate that a break in the TFR with a steeper slope at the low-mass end is a natural consequence of galaxy models with a mass-dependent stellar to halo-mass ratio.


title: On the dimension of contact loci and the identifiability of tensors

abstract: Let $X\subset \mathbb{P}^r$ be an integral and non-degenerate variety. Set $n:= \dim (X)$. We prove that if the $(k+n-1)$-secant variety of $X$ has (the expected) dimension $(k+n-1)(n+1)-1<r$ and $X$ is not uniruled by lines, then $X$ is not $k$-weakly defective and hence the $k$-secant variety satisfies identifiability, i.e. a general element of it is in the linear span of a unique $S\subset X$ with $\sharp (S) =k$. We apply this result to many Segre-Veronese varieties and to the identifiability of Gaussian mixtures $G_{1,d}$. If $X$ is the Segre embedding of a multiprojective space we prove identifiability for the $k$-secant variety (assuming that the $(k+n-1)$-secant variety has dimension $(k+n-1)(n+1)-1<r$, this is a known result in many cases), beating several bounds on the identifiability of tensors.


title: Prospect for Measuring ${G_E^n}$ at High Momentum Transfers

abstract: Experiment E02-013, approved by PAC21, will measure the neutron electric form factor at \qsq~up to 3.4 \gvsq, which is twice that achieved to date. The main features of the new experiment will be the use of the electron spectrometer BigBite, a large array of neutron detectors, and a polarized \he3~target. We present the parameters and optimization of the experimental setup. A concept of an experiment for \GEN~where precision \GEP~data is used for calibration of the systematics of a Rosenbluth type measurement is also discussed. A concept of the circulating gas flow in the polarized He$^3$ target is presented.


title: The intricate Galaxy disk: velocity asymmetries in Gaia-TGAS

abstract: We use the Gaia-TGAS data to compare the transverse velocities in Galactic longitude (coming from proper motions and parallaxes) in the Milky Way disk for negative and positive longitudes as a function of distance. The transverse velocities are strongly asymmetric and deviate significantly from the expectations for an axisymmetric Galaxy. The value and sign of the asymmetry changes at spatial scales of several tens of degrees in Galactic longitude and about 0.5 kpc in distance. The asymmetry is statistically significant at 95% confidence level for 57% of the region probed, which extends up to ~1.2 kpc. A percentage of 24% of the region studied shows absolute differences at this confidence level larger than 5 km/s and 7% larger than 10 km/s. The asymmetry pattern shows mild variations in the vertical direction and with stellar type. A first qualitative comparison with spiral arm models indicates that the arms are unlikely to be the main source of the asymmetry. We briefly discuss alternative origins. This is the first time that global all-sky asymmetries are detected in the Milky Way kinematics, beyond the local neighbourhood, and with a purely astrometric sample.


title: Standard FITS template for simulated astrophysical scenes with the WFIRST coronagraph

abstract: The science investigation teams (SITs) for the WFIRST coronagraphic instrument have begun studying the capabilities of the instrument to directly image reflected light off from exoplanets at contrasts down to contrasts of ~10^-9 with respect to the stellar flux. Detection of point sources at these high contrasts requires yield estimates and detailed modeling of the image of the planetary system as it propagates through the telescope optics. While the SITs might generate custom astrophysical scenes, the integrated model, propagated through the internal speckle field, is typically done at JPL. In this white paper, we present a standard file format to ensure a single distribution system between those who produce the raw astrophysical scenes, and JPL modelers who incorporate those scenes into their optical modeling. At its core, our custom file format uses FITS files, and incorporates standards on packaging astrophysical scenes. This includes spectral and astrometric information for planetary and stellar point sources, zodiacal light and extragalactic sources that may appear as contaminants. Adhering to such a uniform data distribution format is necessary, as it ensures seamless work flow between the SITs and modelers at JPL for the goals of understanding limits of the WFIRST coronagraphic instrument.


title: ALMA and VLA observations of emission from the environment of Sgr A*

abstract: We present 44 and 226 GHz observations of the Galactic center within 20$"$ of Sgr A*. Millimeter continuum emission at 226 GHz is detected from eight stars that have previously been identified at near-IR and radio wavelengths. We also detect a 5.8 mJy source at 226 GHz coincident with the magnetar SGR~J1745-29 located 2.39$"$ SE of Sgr A* and identify a new 2.5$"\times1.5"$ halo of mm emission centered on Sgr A*. The X-ray emission from this halo has been detected previously and is interpreted in terms of a radiatively inefficient accretion flow. The mm halo surrounds an EW linear feature which appears to arise from Sgr A* and coincides with the diffuse X-ray emission and a minimum in the near-IR extinction. We argue that the millimeter emission is produced by synchrotron emission from relativistic electrons in equipartition with a $\sim 1.5$mG magnetic field. The origin of these is unclear but its coexistence with hot gas supports scenarios in which the gas is produced by the interaction of winds either from the fast moving S-stars, the photo-evaporation of low-mass YSO disks or by a jet-driven outflow from Sgr A*. The spatial anti-correlation of the X-ray, radio and mm emission from the halo and the low near-IR extinction provides compelling evidence for an outflow sweeping up the interstellar material, creating a dust cavity within 2$"$ of Sgr A*. Finally, the radio and mm counterparts to eight near-IR identified stars within $\sim$10\arcs\ of Sgr A* provide accurate astrometry to determine the positional shift between the peak emission at 44 and 226 GHz.


title: Radiative effects during the assembly of direct collapse black holes

abstract: We perform a post-processing radiative feedback analysis on a 3D ab initio cosmological simulation of an atomic cooling halo under the direct collapse black hole (DCBH) scenario. We maintain the spatial resolution of the simulation by incorporating native ray-tracing on unstructured mesh data, including Monte Carlo Lyman-alpha (Ly{\alpha}) radiative transfer. DCBHs are born in gas-rich, metal-poor environments with the possibility of Compton-thick conditions, $N_H \gtrsim 10^{24} {\rm cm}^{-2}$. Therefore, the surrounding gas is capable of experiencing the full impact of the bottled-up radiation pressure. In particular, we find that multiple scattering of Ly{\alpha} photons provides an important source of mechanical feedback after the gas in the sub-parsec region becomes partially ionized, avoiding the bottleneck of destruction via the two-photon emission mechanism. We provide detailed discussion of the simulation environment, expansion of the ionization front, emission and escape of Ly{\alpha} radiation, and Compton scattering. A sink particle prescription allows us to extract approximate limits on the post-formation evolution of the radiative feedback. Fully coupled Ly{\alpha} radiation hydrodynamics will be crucial to consider in future DCBH simulations.


title: Flow properties and hydrodynamic interactions of rigid spherical microswimmers

abstract: We analyze a minimal model for a rigid spherical microswimmer and explore the consequences of its extended surface on the interplay between its self-propulsion and flow properties. The model is the first order representation of microswimmers, such as bacteria and algae, with rigid bodies and flexible propelling appendages. The flow field of such a microswimmer at finite distances significantly differs from that of a point-force (Stokeslet) dipole. For a suspension of microswimmers, we derive the grand mobility matrix that connects the motion of an individual swimmer to the active and passive forces and torques acting on all the swimmers. Our investigation of the mobility tensors reveals that hydrodynamic interactions among rigid-bodied microswimmers differ considerably from those among the corresponding point-force dipoles. Our results are relevant for the study of collective behavior of hydrodynamically interacting microswimmers by means of Stokesian dynamics simulations at moderate concentrations.


title: Transport through correlated systems with density functional theory

abstract: We present recent advances in Density Functional Theory (DFT) for applications to the field of quantum transport, with particular emphasis on transport through strongly correlated systems. We review the foundations of the popular Landauer-B\"uttiker(LB)+DFT approach. This formalism, when using approximations to the exchange-correlation (xc) potential with steps at integer occupation, correctly captures the Kondo plateau in the zero bias conductance at zero temperature but completely fails to capture the transition to the Coulomb blockade (CB) regime as temperature increases. To overcome the limitations of LB+DFT the quantum transport problem is treated from a time-dependent (TD) perspective using TDDFT, an exact framework to deal with nonequilibrium situations. The steady-state limit of TDDFT shows that in addition to an xc potential in the junction, there also exists an xc correction to the applied bias. Open shell molecules in the CB regime provide the most striking examples of the importance of the xc bias correction. Using the Anderson model as guidance we estimate these corrections in the limit of zero bias. For the general case we put forward a steady-state DFT which is based on the one-to-one correspondence between the pair of basic variables steady density on and steady current across the junction and the pair local potential on and bias across the junction. Like TDDFT, this framework also leads to both an xc potential in the junction and an xc correction to the bias. Unlike in TDDFT, these potentials are independent of history. We highlight the universal features of both xc potential and xc bias corrections for junctions in the CB regime and provide an accurate parametrization for the Anderson model at arbitrary temperatures and interaction strengths thus providing a unified DFT description for both Kondo and CB regimes and the transition between them.


title: Statistically Characterizing the Electrical Parameters of the Grid Transformers and Transmission Lines

abstract: This paper presents a set of validation metrics for transmission network parameters that is applicable in both creation of synthetic power system test cases and validation of existing models. Using actual data from two real-world power grids, statistical analyses are performed to extract some useful statistics on transformers and transmission lines electrical parameters including per unit reactance, MVA rating, and their X/R ratio. It is found that conversion of per unit reactance calculated on system common base to transformer own power base will significantly stabilize its range and remove the correlation between per unit X and MVA rating. This is fairly consistent for transformers with different voltage levels and sizes and can be utilized as a strong validation metric for synthetic models. It is found that transmission lines exhibit different statistical properties than transformers with different distribution and range for the parameters. In addition, statistical analysis shows that the empirical PDF of transmission network electrical parameters can be approximated with mathematical distribution functions which would help appropriately characterize them in synthetic power networks. Kullback-Leibler divergence is used as a measure for goodness of fit for approximated distributions.


title: Visualization of Thomas-Wigner rotations

abstract: It is well known that a sequence of two non-collinear pure Lorentz transformations (boosts) is not a boost again, but involves a spatial rotation, the Wigner or Thomas-Wigner rotation. The formation of this rotation is visually illustrated by moving a Born-rigid object on a closed trajectory in several sections. Within each section the boost's proper time duration is assumed to be the same and the object's centre accelerates uniformly. Born-rigidity implies that the stern of this object accelerates faster than its bow. It is shown that at least five boosts are required to return the object's centre to its start position. With these assumptions, the Thomas-Wigner rotation angle depends on a single parameter only, the maximum speed reached within each boost section. The visualization highlights the close relationship between the Thomas-Wigner rotation and the relativity of simultaneity. Furthermore, it is illustrated that accelerated motion implies the formation of an event horizon. The event horizons associated with the five boosts constitute a boundary to the rotated Born-rigid object and ensure its finite size.


title: On Some Applications of Group Representation Theory to Algebraic Problems Related to the Congruence Principle for Equivariant Maps

abstract: Given a finite group $G$ and two unitary $G$-representations $V$ and $W$, possible restrictions on Brouwer degrees of equivariant maps between representation spheres $S(V)$ and $S(W)$ are usually expressed in a form of congruences modulo the greatest common divisor of lengths of orbits in $S(V)$ (denoted $\alpha(V)$). Effective applications of these congruences is limited by answers to the following questions: (i) under which conditions, is ${\alpha}(V)>1$? and (ii) does there exist an equivariant map with the degree easy to calculate? In the present paper, we address both questions. We show that ${\alpha}(V)>1$ for each irreducible non-trivial $C[G]$-module if and only if $G$ is solvable. For non-solvable groups, we use 2-transitive actions to construct complex representations with non-trivial ${\alpha}$-characteristic. Regarding the second question, we suggest a class of Norton algebras without 2-nilpotents giving rise to equivariant quadratic maps, which admit an explicit formula for the Brouwer degree.


title: Sympathy Begins with a Smile, Intelligence Begins with a Word: Use of Multimodal Features in Spoken Human-Robot Interaction

abstract: Recognition of social signals, from human facial expressions or prosody of speech, is a popular research topic in human-robot interaction studies. There is also a long line of research in the spoken dialogue community that investigates user satisfaction in relation to dialogue characteristics. However, very little research relates a combination of multimodal social signals and language features detected during spoken face-to-face human-robot interaction to the resulting user perception of a robot. In this paper we show how different emotional facial expressions of human users, in combination with prosodic characteristics of human speech and features of human-robot dialogue, correlate with users' impressions of the robot after a conversation. We find that happiness in the user's recognised facial expression strongly correlates with likeability of a robot, while dialogue-related features (such as number of human turns or number of sentences per robot utterance) correlate with perceiving a robot as intelligent. In addition, we show that facial expression, emotional features, and prosody are better predictors of human ratings related to perceived robot likeability and anthropomorphism, while linguistic and non-linguistic features more often predict perceived robot intelligence and interpretability. As such, these characteristics may in future be used as an online reward signal for in-situ Reinforcement Learning based adaptive human-robot dialogue systems.


title: Universal in vivo Textural Model for Human Skin based on Optical Coherence Tomograms

abstract: Currently, diagnosis of skin diseases is based primarily on visual pattern recognition skills and expertise of the physician observing the lesion. Even though dermatologists are trained to recognize patterns of morphology, it is still a subjective visual assessment. Tools for automated pattern recognition can provide objective information to support clinical decision-making. Noninvasive skin imaging techniques provide complementary information to the clinician. In recent years, optical coherence tomography has become a powerful skin imaging technique. According to specific functional needs, skin architecture varies across different parts of the body, as do the textural characteristics in OCT images. There is, therefore, a critical need to systematically analyze OCT images from different body sites, to identify their significant qualitative and quantitative differences. Sixty-three optical and textural features extracted from OCT images of healthy and diseased skin are analyzed and in conjunction with decision-theoretic approaches used to create computational models of the diseases. We demonstrate that these models provide objective information to the clinician to assist in the diagnosis of abnormalities of cutaneous microstructure, and hence, aid in the determination of treatment. Specifically, we demonstrate the performance of this methodology on differentiating basal cell carcinoma (BCC) and squamous cell carcinoma (SCC) from healthy tissue.


title: Local behavior of local times of super Brownian motion

abstract: For $x\in R^d- \{0\}$, in dimension $d=3$, we study the asymptotic behavior of the local time $L_t^x$ of super-Brownian motion $X$ starting from $\delta_0$ as $x \to 0$. Let $\psi(x)=((1/2\pi^2) \log (1/|x|))^{1/2}$ be a normalization, Theorem 1 implies that $(L_t^x-(1/2\pi|x|))/\psi(x)$ converges in distribution to a standard normal distributed random variable as $x\to 0$. For dimension $d=2$, Theorem 2 implies that $L^x_t-(1/\pi)\log(1/|x|)$ is $L^1$ bounded as $x\to 0$. To do this, we prove a Tanaka formula for the local time which refines a result in Barlow, Evans and Perkins.


title: Structure and interpolation of the turbulent velocity profile in parallel flow

abstract: The classical scaling theory of turbulent parallel flow provides a framework for the description of the mean velocity profile through two functions of one variable, traditionally named law of the wall and law of the wake, and a universal logarithmic law characterized by von K\'arm\'an's constant. Despite its widespread adoption in research and in teaching, discrepancies between this theory and both experiments and numerical simulations have been repeatedly observed in the literature. Recently we have shown that in the logarithmic layer such discrepancies can be physically interpreted and analytically accounted for through an equally universal correction caused by the pressure gradient. This finding opens the way to a likewise improvement in the description of the law of the wall and of the law of the wake, an analytical interpolation of either of which is often useful for practical applications.


title: Gated Orthogonal Recurrent Units: On Learning to Forget

abstract: We present a novel recurrent neural network (RNN) based model that combines the remembering ability of unitary RNNs with the ability of gated RNNs to effectively forget redundant/irrelevant information in its memory. We achieve this by extending unitary RNNs with a gating mechanism. Our model is able to outperform LSTMs, GRUs and Unitary RNNs on several long-term dependency benchmark tasks. We empirically both show the orthogonal/unitary RNNs lack the ability to forget and also the ability of GORU to simultaneously remember long term dependencies while forgetting irrelevant information. This plays an important role in recurrent neural networks. We provide competitive results along with an analysis of our model on many natural sequential tasks including the bAbI Question Answering, TIMIT speech spectrum prediction, Penn TreeBank, and synthetic tasks that involve long-term dependencies such as algorithmic, parenthesis, denoising and copying tasks.


title: Time-Series Analysis of Photovoltaic Distributed Generation Impacts on a Local Distributed Network

abstract: Increasing penetration level of photovoltaic (PV) distributed generation (DG) into distribution networks will have many impacts on nominal circuit operating conditions including voltage quality and reverse power flow issues. In U.S. most studies on PVDG impacts on distribution networks are performed for west coast and central states. The objective of this paper is to study the impacts of PVDG integration on local distribution network based on real-world settings for network parameters and time-series analysis. PVDG penetration level is considered to find the hosting capacity of the network without having major issues in terms of voltage quality and reverse power flow. Time-series analyses show that distributed installation of PVDGs on commercial buses has the maximum network energy loss reduction and larger penetration ratios for them. Additionally, the penetration ratio thresholds for which there will be no power quality and reverse power flow issues and optimal allocation of PVDG and penetration levels are identified for different installation scenarios.


title: On the Faithfulness of 1-dimensional Topological Quantum Field Theories

abstract: This paper explores 1-dimensional topological quantum field theories. We separately deal with strict and strong 1-dimensional topological quantum field theories. The strict one is regarded as a symmetric monoidal functor between the category of 1-cobordisms and the category of matrices, and the strong one is a symmetric monoidal functor between the category of 1-cobordisms and the category of finite dimensional vector spaces. It has been proved that both strict and strong 1-dimensional topological quantum field theories are faithful.


title: What modern vision science reveals about the awareness puzzle: Summary-statistic encoding plus decision limits underlie the richness of visual perception and its quirky failures

abstract: There is a fundamental puzzle in understanding our awareness of the visual world. On one hand, our subjective experience is one of a rich visual world, which we perceive effortlessly. However, when we actually test perception, observers know surprisingly little. A number of tasks, from search, through inattentional blindness, to change blindness, suggest that there is surprisingly little awareness or perception without attention. Meanwhile, another set of tasks, such as multiple object tracking, dual-task performance, and visual working memory tasks suggest that both attention and working memory have low capacity. These two components together - poor perception without attention, and greatly limited capacity for attention and memory - imply that perception is impoverished. How can we make sense of this awareness puzzle, of the riddle of our rich subjective experience coupled with poor performance on experimental tasks? I suggest that, looked at in the right way, there is in fact no awareness puzzle. In particular, I will argue that the tasks that show limits are inherently difficult tasks, and that there exists a unified explanation for both the rich subjective experience and the apparent limits.


title: Perfect Anomalous Reflection with an Aggressively Discretized Huygens' Metasurface

abstract: This paper investigates the discretization of a periodic metasurface and demonstrates how such a surface can achieve perfect anomalous reflection. Whilst most contemporary theoretical works on metasurfaces deal with continuous current or impedance distributions, we examine how discretization affects a metasurface, and show that in some cases one can discretize a metasurface aggressively --- to the extent of having only two cells per spatial period. Such aggressive discretization can lead to great simplifications in metasurface design, and perhaps more surprisingly, a possible performance improvement from continuous metasurfaces. Using this aggressive discretization technique, we report the design of a binary Huygens' metasurface which reflects an incident plane wave at 50$^\circ$ into a reflected direction of -22.5$^\circ$. Full-wave electromagnetic simulation shows the achievement of anomalous reflection with a power efficiency of 99.1%, which dramatically surpasses the performance of a corresponding passive continuous metasurface, for which the power efficiency is fundamentally limited to 69.6%.


title: Evolutionary Multitasking for Multiobjective Continuous Optimization: Benchmark Problems, Performance Metrics and Baseline Results

abstract: In this report, we suggest nine test problems for multi-task multi-objective optimization (MTMOO), each of which consists of two multiobjective optimization tasks that need to be solved simultaneously. The relationship between tasks varies between different test problems, which would be helpful to have a comprehensive evaluation of the MO-MFO algorithms. It is expected that the proposed test problems will germinate progress the field of the MTMOO research.


title: Extending the topological analysis and seeking the real-space subsystems in non-Coulombic systems with homogeneous potential energy functions

abstract: It is customary to conceive the interactions of all the constituents of a molecular system, i.e. electrons and nuclei, as Coulombic. However, in a more detailed analysis one may always find small but non-negligible non-Coulombic interactions in molecular systems originating from the finite size of nuclei, magnetic interactions, etc. While such small modifications of the Coulombic interactions do not seem to alter the nature of a molecular system in real world seriously, they are a serious obstacle for quantum chemical theories and methodologies which their formalism is strictly confined to the Coulombic interactions. Although the quantum theory of atoms in molecules (QTAIM) has been formulated originally for the Coulombic systems, some recent studies have demonstrated that apart from basin energy of an atom in a molecule, its theoretical ingredients are not sensitive to the explicit form of the potential energy operator. In this study, it is demonstrated that the basin energy may be defined not only for coulombic systems but for all real-space subsystems of those systems that are described by any member of the set of the homogeneous potential energy functions. On the other hand, this extension opens the door for seeking novel real-space subsystems, apart from atoms in molecules, in non-Coulombic systems. These novel real-space subsystems call for an extended formalism that goes beyond the orthodox QTAIM, which is not confined to the Coulombic systems nor to the atoms in molecules as the sole real-space subsystems. It is termed the quantum theory of real-space open subsystems (QTROS) and its potential applications are detailed. The harmonic trap model, containing non-interacting fermions or bosons, is considered as an example for the QTROS analysis. The QTROS analysis of bosonic systems is particularly quite unprecedented, not attempted before.


title: Random projections for linear programming

abstract: Random projections are random linear maps, sampled from appropriate distributions, that approx- imately preserve certain geometrical invariants so that the approximation improves as the dimension of the space grows. The well-known Johnson-Lindenstrauss lemma states that there are random ma- trices with surprisingly few rows that approximately preserve pairwise Euclidean distances among a set of points. This is commonly used to speed up algorithms based on Euclidean distances. We prove that these matrices also preserve other quantities, such as the distance to a cone. We exploit this result to devise a probabilistic algorithm to solve linear programs approximately. We show that this algorithm can approximately solve very large randomly generated LP instances. We also showcase its application to an error correction coding problem.


title: Source Forager: A Search Engine for Similar Source Code

abstract: Developers spend a significant amount of time searching for code: e.g., to understand how to complete, correct, or adapt their own code for a new context. Unfortunately, the state of the art in code search has not evolved much beyond text search over tokenized source. Code has much richer structure and semantics than normal text, and this property can be exploited to specialize the code-search process for better querying, searching, and ranking of code-search results. We present a new code-search engine named Source Forager. Given a query in the form of a C/C++ function, Source Forager searches a pre-populated code database for similar C/C++ functions. Source Forager preprocesses the database to extract a variety of simple code features that capture different aspects of code. A search returns the $k$ functions in the database that are most similar to the query, based on the various extracted code features. We tested the usefulness of Source Forager using a variety of code-search queries from two domains. Our experiments show that the ranked results returned by Source Forager are accurate, and that query-relevant functions can be reliably retrieved even when searching through a large code database that contains very few query-relevant functions. We believe that Source Forager is a first step towards much-needed tools that provide a better code-search experience.


title: Optimized Solar Photovoltaic Generation in a Real Local Distribution Network

abstract: Remarkable penetration of renewable energy in electric networks, despite its valuable opportunities, such as power loss reduction and loadability improvements, has raised concerns for system operators. Such huge penetration can lead to a violation of the grid requirements, such as voltage and current limits and reverse power flow. Optimal placement and sizing of Distributed Generation (DG) are one of the best ways to strengthen the efficiency of the power systems. This paper builds a simulation model for the local distribution network based on obtained load profiles, GIS information, solar insolation, feeder and voltage settings, and define the optimization problem of solar PVDG installation to determine the optimal siting and sizing for different penetration levels with different objective functions. The objective functions include voltage profile improvement and energy loss minimization and the considered constraints include the physical distribution network constraints (AC power flow), the PV capacity constraint, and the voltage and reverse power flow constraints.


title: On the global sup-norm of GL(3) cusp forms

abstract: Let $\phi$ be a spherical Hecke-Maass cusp form on the non-compact space $\mathrm{PGL}_3(\mathbb{Z})\backslash\mathrm{PGL}_3(\mathbb{R})$. We establish various pointwise upper bounds for $\phi$ in terms of its Laplace eigenvalue $\lambda_\phi$. These imply, for $\phi$ arithmetically normalized and tempered at the archimedean place, the bound $\|\phi\|_\infty\ll_\epsilon \lambda_{\phi}^{39/40+\epsilon}$ for the global sup-norm (without restriction to a compact subset). On the way, we derive a new uniform upper bound for the $\mathrm{GL}_3$ Jacquet-Whittaker function.


title: The Intensities of Cosmic Ray H and He Nuclei at ~250 MeV/nuc Measured by Voyagers 1 and 2 - Using these Intensities to Determine the Solar Modulation Parameter in the Inner Heliosphere and the Heliosheath Over a 40 Year Time Period

abstract: We have determined the solar modulation potential, phi, vs. time that is observed at Voyager 1 and 2 from measurements of the H and He nuclei intensities at a common energy of 250 MeVnuc. The H nuclei have a rigidity 0.7 GV, the He nuclei 1.4 GV. These measurements cover a 40 year time period, which includes almost 4 cycles of solar 11 year sunspot variations, throughout the inner heliosphere out to the HTS at distances of 95 AU and 85 AU, respectively at V1 and V2, and then beyond in the heliosheath. Inside the HTS the modulation potential vs. time curves at V1 and V2 show a very similar temporal structure to those observed at the Earth. During a later period of maximum solar modulation from 2000.0 to 2005.0 when V1 and V2 are in the outer heliosphere between 60-94 AU, the main temporal features of the modulation potential curves at all 3 locations match up with appropriate time delays at V1 and V2 if it is assumed that spatially coherent structures are moving outward past V1 and V2, with outward speeds of up to 700 Kms negative 1. After 2004.0 V1 and V2 are at latitudes of positive 35 and negative 30 respectively, placing lower limits on the latitude extent of these structures. Beyond the HTS in the heliosheath the modulation potential slowly decreases at both spacecraft with only a weak evidence of the unusual modulation minimum observed at the Earth in 2009, for example. A sudden decrease of the modulation potential 50 MV for both H and He nuclei occurs at V1 just before the heliopause crossing at about 122 AU. This decrease has not yet been observed at V2, which is now at 113 AU and still observing a modulation potential 60 MV.


title: The Raman fingerprint of rhombohedral graphite

abstract: Multi-layer graphene with rhombohedral stacking is a promising carbon phase possibly displaying correlated states like magnetism or superconductivity due to the occurrence of a flat surface band at the Fermi level. Recently, flakes of thickness up to 17 layers were tentatively attributed ABC sequences although the Raman fingerprint of rhombohedral multilayer graphene is currently unknown and the 2D resonant Raman spectrum of Bernal graphite not understood. We provide a first principles description of the 2D Raman peak in three and four layers graphene (all stackings) as well as in Bernal, rhombohedral and an alternation of Bernal and rhombohedral graphite. We give practical prescriptions to identify long range sequences of ABC multi-layer graphene. Our work is a prerequisite to experimental non-destructive identification and synthesis of rhombohedral graphite.


title: Experimental Demonstration of an Electride as a 2D Material

abstract: Because of their loosely bound electrons, electrides offer physical properties useful in chemical synthesis and electronics. For these applications and others, nano-sized electrides offer advantages, but to-date no electride has been synthesized as a nanomaterial. We demonstrate experimentally that Ca$_2$N, a layered electride in which layers of atoms are separated by layers of a 2D electron gas (2DEG), can be exfoliated into two-dimensional (2D) nanosheets using liquid exfoliation. The 2D flakes are stable in a nitrogen atmosphere or in select organic solvents for at least one month. Electron microscopy and elemental analysis reveal that the 2D flakes retain the crystal structure and stoichiometry of the parent 3D Ca$_2$N. In addition, the 2D flakes exhibit metallic character and an optical response that agrees with DFT calculations. Together these findings suggest that the 2DEG is preserved in the 2D material. With this work, we bring electrides into the nano-regime and experimentally demonstrate a 2D electride, Ca$_2$N.


title: A Mathematical Framework for Kinetochore-Driven Activation Feedback in the Mitotic Checkpoint

abstract: Proliferating cells properly divide into their daughter cells through a process that is mediated by kinetochores, protein-complexes that assemble at the centromere of each sister chromatid. Each kinetochore has to establish a tight bipolar attachment to the spindle apparatus before sister-chromatid separation is initiated. The Spindle Assembly Checkpoint (SAC) links the biophysical attachment status of the kinetochores to mitotic progression, and ensures that even a single misaligned kinetochore keeps the checkpoint active. The mechanism by which this is achieved is still elusive. Current computational models of the human SAC disregard important biochemical properties by omitting any kind of feedback loop, proper kinetochore signals, and other spatial properties such as the stability of the system and diffusion effects. To allow for more realistic in silico study of the dynamics of the SAC model, a minimal mathematical framework for SAC activation and silencing is introduced. A nonlinear ordinary differential equation model successfully reproduces bifurcation signaling switches with attachment of all 92 kinetochores and activation of APC/C by kinetochore-driven feedback. A partial differential equation model and mathematical linear stability analyses indicate the influence of diffusion and system stability. The conclusion is that quantitative models of the human SAC should account for the positive feedback on APC/C activation driven by the kinetochores which is essential for SAC silencing. Experimental diffusion coefficients for MCC sub-complexes are found to be insufficient for rapid APC/C inhibition. The presented analysis allows for systems-level understanding of mitotic control and the minimal new model can function as a basis for developing further quantitative-integrative models of the cell division cycle


title: Optimizing expected word error rate via sampling for speech recognition

abstract: State-level minimum Bayes risk (sMBR) training has become the de facto standard for sequence-level training of speech recognition acoustic models. It has an elegant formulation using the expectation semiring, and gives large improvements in word error rate (WER) over models trained solely using cross-entropy (CE) or connectionist temporal classification (CTC). sMBR training optimizes the expected number of frames at which the reference and hypothesized acoustic states differ. It may be preferable to optimize the expected WER, but WER does not interact well with the expectation semiring, and previous approaches based on computing expected WER exactly involve expanding the lattices used during training. In this paper we show how to perform optimization of the expected WER by sampling paths from the lattices used during conventional sMBR training. The gradient of the expected WER is itself an expectation, and so may be approximated using Monte Carlo sampling. We show experimentally that optimizing WER during acoustic model training gives 5% relative improvement in WER over a well-tuned sMBR baseline on a 2-channel query recognition task (Google Home).


title: Summary Analysis of the 2017 GitHub Open Source Survey

abstract: This report is a high-level summary analysis of the 2017 GitHub Open Source Survey dataset, presenting frequency counts, proportions, and frequency or proportion bar plots for every question asked in the survey.


title: Equality in Brascamp-Lieb-Luttinger Inequalities

abstract: An inequality of Brascamp-Lieb-Luttinger generalizes the Riesz-Sobolev inequality, stating that certain multilinear functionals, acting on nonnegative functions of one real variable with prescribed distribution functions, are maximized when these functions are symmetrized. It is shown that under certain hypotheses, when the functions are indicator functions of sets of prescribed measures, then up to the natural translation symmetries of the inequality, the maximum is attained only by intervals centered at the origin. Moreover, a quantitative form of this uniqueness is established, sharpening the inequality. The hypotheses include an auxiliary genericity assumption which may not be necessary.


title: A possibly universal red chromophore for modeling color variations on Jupiter

abstract: A new laboratory-generated chemical compound made from photodissociated ammonia (NH3) molecules reacting with acetylene (C2H2) was suggested as a possible coloring agent for Jupiter's Great Red Spot (GRS) by Carlson et al. (2016, Icarus 274, 106-115). Baines et al. (2016, Icarus, submitted) showed that the GRS spectrum measured by the visual channels of the Cassini VIMS instrument in 2000 could be accurately fit by a cloud model in which the chromophore appeared as a physically thin layer of small particles immediately above the main cloud layer of the GRS. Here we show that the same chromophore and same layer location can also provide close matches to the short wave spectra of many other cloud features on Jupiter, suggesting this material may be a nearly universal chromophore that could explain the various degrees of red coloration on Jupiter. This is a robust conclusion, even for 12% changes in VIMS calibration and large uncertainties in the refractive index of the main cloud layer due to uncertain fractions of NH4SH and NH3 in its cloud particles. The chromophore layer can account for color variations among north and south equatorial belts, equatorial zone, and the Great Red Spot, by varying particle size from 0.12 microns to 0.29 microns and 1-micron optical depth from 0.06 to 0.76. The total mass of the chromophore layer is much less variable, ranging from 18 to 30 micrograms/cm^2, except in the equatorial zone, where it is only 10-13 micrograms/cm^2. We also found a depression of the ammonia volume mixing ratio in the two belt regions, which averaged 0.4-0.5 X 10^{-4} immediately below the ammonia condensation level, while the other regions averaged twice that value.


title: Setting Players' Behaviors in World of Warcraft through Semi-Supervised Learning

abstract: Digital games are one of the major and most important fields on the entertainment domain, which also involves cinema and music. Numerous attempts have been done to improve the quality of the games including more realistic artistic production and computer science. Assessing the player's behavior, a task known as player modeling, is currently the need of the hour which leads to possible improvements in terms of: (i) better game interaction experience, (ii) better exploitation of the relationship between players, and (iii) increasing/maintaining the number of players interested in the game. In this paper we model players using the basic four behaviors proposed in \cite{BartleArtigo}, namely: achiever, explorer, socializer and killer. Our analysis is carried out using data obtained from the game "World of Warcraft" over 3 years (2006 $-$ 2009). We employ a semi-supervised learning technique in order to find out characteristics that possibly impact player's behavior.


title: Granger Causality Networks for Categorical Time Series

abstract: We present a new framework for learning Granger causality networks for multivariate categorical time series, based on the mixture transition distribution (MTD) model. Traditionally, MTD is plagued by a nonconvex objective, non-identifiability, and presence of many local optima. To circumvent these problems, we recast inference in the MTD as a convex problem. The new formulation facilitates the application of MTD to high-dimensional multivariate time series. As a baseline, we also formulate a multi-output logistic autoregressive model (mLTD), which while a straightforward extension of autoregressive Bernoulli generalized linear models, has not been previously applied to the analysis of multivariate categorial time series. We develop novel identifiability conditions of the MTD model and compare them to those for mLTD. We further devise novel and efficient optimization algorithm for the MTD based on the new convex formulation, and compare the MTD and mLTD in both simulated and real data experiments. Our approach simultaneously provides a comparison of methods for network inference in categorical time series and opens the door to modern, regularized inference with the MTD model.


title: $f(R)$ global monopole revisited

abstract: In this paper the $f(R)$ global monopole is reexamined. We provide an exact solution for the modified field equations in the presence of a global monopole for regions outside its core, generalizing previous results. Additionally, we discuss some particular cases obtained from this solution. We consider a setup consisting of a possible Schwarzschild black hole that absorbs the topological defect, giving rise to a static black hole endowed with a monopole's charge. Besides, we demonstrate how the asymptotic behavior of the Higgs field far from the monopole's core is shaped by a class of spacetime metrics which includes those ones analyzed here. In order to assess the gravitational properties of this system, we analyse the geodesic motion of both massive and massless test particles moving in the vicinity of such configuration. For the material particles we set the requirements they have to obey in order to experience stable orbits. On the other hand, for the photons we investigate how their trajectories are affected by the gravitational field of this black hole.


title: Linear Hashing is Awesome

abstract: We consider the hash function $h(x) = ((ax+b) \bmod p) \bmod n$ where $a,b$ are chosen uniformly at random from $\{0,1,\ldots,p-1\}$. We prove that when we use $h(x)$ in hashing with chaining to insert $n$ elements into a table of size $n$ the expected length of the longest chain is $\tilde{O}\!\left(n^{1/3}\right)$. The proof also generalises to give the same bound when we use the multiply-shift hash function by Dietzfelbinger et al. [Journal of Algorithms 1997].


title: The extra scalar degrees of freedom from the two Higgs doublet model for dark energy

abstract: In principle a minimal extension of the standard model of Particle Physics, the two Higgs doublet model, can be invoked to explain the scalar field responsible of dark energy. The two doublets are in general mixed. After diagonalization, the lightest CP-even Higgs and CP-odd Higgs are jointly taken to be the dark energy candidate. The dark energy obtained from Higgs fields in this case is indistinguishable from the cosmological constant.


title: Optimal parameters for bloom-filtered joins in Spark

abstract: In this paper, we present an algorithm that joins relational database tables efficiently in a distributed environment using Bloom filters of an optimal size. We propose not to use fixed-size bloom filters as in previous research, but to find an optimal size for the bloom filters, by creating a mathematical model of the join algorithm, and then finding the optimal parameters using traditional mathematical optimization. This algorithm with optimal parameters beats both previous approaches using bloom filters and the default SparkSQL engine not only on star-joins, but also on traditional database schema. The experiments were conducted on a standard TPC-H database stored as parquet files on a distributed file system.


title: The Evolution of Massive Stars: Bridging the Gap in the Local Group

abstract: The nearby galaxies of the Local Group can act as our laboratories in helping to bridge the gap between theory and observations. In this review we will describe the complications of identifying samples of OB stars, yellow and red supergiants, and Wolf-Rayet stars, and what we have so far learned from these studies.


title: Theoretical and experimental evidence of a site-selective Mott transition in Fe2O3 under pressure

abstract: We provide experimental and theoretical evidence for a novel type of pressure-induced insulator-metal transition characterized by site-selective delocalization of the electrons. M\"ossbauer spectroscopy, X-ray diffraction and electrical transport measurements on Fe$_2$O$_3$ to 100 GPa, along with dynamical mean-field theory (DFT+DMFT) calculations, reveal this site-selective Mott transition between 50 and 68 GPa, such that the metallization can be described by ($^\rm{VI}$Fe$^{3+\rm{HS}}$)$_2$O$_3$ [$R\bar{3}c$ structure] $\overrightarrow{\tiny\rm{50~GPa}}$ ($^\rm{VIII}$Fe$^{3+\rm{HS~VI}}$Fe$^\rm{M}$)O$_3$ [$P2_1/n$ structure] $\overrightarrow{\tiny\rm{68~GPa}}$ ($^\rm{VI}$Fe$^\rm{M}$)$_2$O$_3$ [$Aba2$ structure]. Within the $P2_1/n$ crystal structure, characterized by two distinct coordination sites (VI and VIII), we observe equal abundances of ferric ions (Fe$^{3+}$) and ions having delocalized electrons (Fe$^\rm{M}$), and only at higher pressures is a fully metallic $Aba2$ structure obtained, all at room temperature. The transition is characterized by delocalization/metallization of the $3d$ electrons on half the Fe sites, with a site-dependent collapse of local moments. Above $\sim$50 GPa, Fe$_2$O$_3$ is a strongly correlated metal with reduced electron mobility (large band renormalizations) of m*/m$\sim$4 and 6 near the Fermi level. Upon decompression, we observe a site-selective (metallic) to conventional Mott insulator phase transition ($^\rm{VIII}$Fe$^{3+\rm{HS~VI}}$Fe$^\rm{M}$)O$_3$ $\overrightarrow{\tiny\rm{50~GPa}}$ ($^\rm{VIII}$Fe$^{3+\rm{HS~VI}}$Fe$^{3+ \rm{HS}}$)O$_3$ within the same $P2_1/n$ structure, indicating a decoupling of the electronic and lattice degrees of freedom, characteristic of a true Mott transition. Our results show that the interplay of electronic correlations and lattice may result in rather complex behavior of the electronic structure and magnetic state.


title: Mass Formula for Light Nonstrange Mesons and Regge Trajectories in Quark Model

abstract: We study the Regge-like spectra of light mesons in a relativized quark model. An analytical mass formula is presented for the light unflavored mesons with the help of auxiliary field method, by which a quasi-linear Regge-Chew-Frautschi plot is predicted for the orbitally excited states. We show that the trajectory slope is proportional to the inverse confining parameter $1/a$ up to a factor depending on the strong coupling $\alpha_{s}$ when the orbital quantum number $L$ is large. The result is tested against the experimental data of the spectra of the meson families $\pi/b$,$\rho/a$,$\eta/h$ and $\omega/f$ in the ($L$,$M^{2}$) planes, with the fitted parameters consistent with that in the literatures.


title: On the Development of Intelligent Agents for MOBA Games

abstract: Multiplayer Online Battle Arena (MOBA) is one of the most played game genres nowadays. With the increasing growth of this genre, it becomes necessary to develop effective intelligent agents to play alongside or against human players. In this paper we address the problem of agent development for MOBA games. We implement a two-layered architecture agent that handles both navigation and game mechanics. This architecture relies on the use of Influence Maps, a widely used approach for tactical analysis. Several experiments were performed using {\em League of Legends} as a testbed, and show promising results in this highly dynamic real-time context.


title: The dangerous irrelevance of string theory

abstract: We comment on the relation between string theory and empirical science, grounding our discussion in cosmology, a subject with increasingly precise data in which this connection operates at several levels. It is important to take into account the phenomenon of dangerous irrelevance: over long times or large field ranges, physics can become sensitive to higher scales than the input energies. This pertains in inflationary cosmology (and possibly other aspects of horizon physics). String theory also contributes to our understanding of observational constraints and search strategies at the level of low energy field theory. We illustrate this with a current example concerning a new form of non-Gaussianity generated by very massive degrees of freedom coupling to the inflaton. New constraints on such fields and couplings can be obtained from existing data, increasing our empirical knowledge of the universe. This builds in part from the development of the string landscape, which is neither random nor an abdication of science as has sometimes been suggested. {\it Invited contribution to the proceedings of the conference `Why trust a theory'.}


title: Scattering of point particles by black holes: gravitational radiation

abstract: Gravitational waves can teach us not only about sources and the environment where they were generated, but also about the gravitational interaction itself. Here we study the features of gravitational radiation produced during the scattering of a point-like mass by a black hole. Our results are exact (to numerical error) at any order in a velocity expansion, and are compared against various approximations. At large impact parameter and relatively small velocities our results agree to within percent level with various post-Newtonian and weak-field results. Further, we find good agreement with scaling predictions in the weak-field/high-energy regime. Lastly, we achieve striking agreement with zero-frequency estimates.


title: The FastMap Algorithm for Shortest Path Computations

abstract: We present a new preprocessing algorithm for embedding the nodes of a given edge-weighted undirected graph into a Euclidean space. The Euclidean distance between any two nodes in this space approximates the length of the shortest path between them in the given graph. Later, at runtime, a shortest path between any two nodes can be computed with A* search using the Euclidean distances as heuristic. Our preprocessing algorithm, called FastMap, is inspired by the data mining algorithm of the same name and runs in near-linear time. Hence, FastMap is orders of magnitude faster than competing approaches that produce a Euclidean embedding using Semidefinite Programming. FastMap also produces admissible and consistent heuristics and therefore guarantees the generation of shortest paths. Moreover, FastMap applies to general undirected graphs for which many traditional heuristics, such as the Manhattan Distance heuristic, are not well defined. Empirically, we demonstrate that A* search using the FastMap heuristic is competitive with A* search using other state-of-the-art heuristics, such as the Differential heuristic.


title: From orbital measures to Littlewood-Richardson coefficients and hive polytopes

abstract: The volume of the hive polytope (or polytope of honeycombs) associated with a Littlewood- Richardson coefficient of SU(n), or with a given admissible triple of highest weights, is expressed, in the generic case, in terms of the Fourier transform of a convolution product of orbital measures. Several properties of this function -- a function of three non-necessarily integral weights or of three multiplets of real eigenvalues for the associated Horn problem-- are already known. In the integral case it can be thought of as a semi-classical approximation of Littlewood-Richardson coefficients. We prove that it may be expressed as a local average of a finite number of such coefficients. We also relate this function to the Littlewood-Richardson polynomials (stretching polynomials) i.e., to the Ehrhart polynomials of the relevant hive polytopes. Several SU(n) examples, for n=2,3,...,6, are explicitly worked out.


title: Rapid Randomized Restarts for Multi-Agent Path Finding Solvers

abstract: Multi-Agent Path Finding (MAPF) is an NP-hard problem well studied in artificial intelligence and robotics. It has many real-world applications for which existing MAPF solvers use various heuristics. However, these solvers are deterministic and perform poorly on "hard" instances typically characterized by many agents interfering with each other in a small region. In this paper, we enhance MAPF solvers with randomization and observe that they exhibit heavy-tailed distributions of runtimes on hard instances. This leads us to develop simple rapid randomized restart (RRR) strategies with the intuition that, given a hard instance, multiple short runs have a better chance of solving it compared to one long run. We validate this intuition through experiments and show that our RRR strategies indeed boost the performance of state-of-the-art MAPF solvers such as iECBS and M*.


title: A Deep Causal Inference Approach to Measuring the Effects of Forming Group Loans in Online Non-profit Microfinance Platform

abstract: Kiva is an online non-profit crowdsouring microfinance platform that raises funds for the poor in the third world. The borrowers on Kiva are small business owners and individuals in urgent need of money. To raise funds as fast as possible, they have the option to form groups and post loan requests in the name of their groups. While it is generally believed that group loans pose less risk for investors than individual loans do, we study whether this is the case in a philanthropic online marketplace. In particular, we measure the effect of group loans on funding time while controlling for the loan sizes and other factors. Because loan descriptions (in the form of texts) play an important role in lenders' decision process on Kiva, we make use of this information through deep learning in natural language processing. In this aspect, this is the first paper that uses one of the most advanced deep learning techniques to deal with unstructured data in a way that can take advantage of its superior prediction power to answer causal questions. We find that on average, forming group loans speeds up the funding time by about 3.3 days.


title: Dynamic Difficulty Adjustment on MOBA Games

abstract: This paper addresses the dynamic difficulty adjustment on MOBA games as a way to improve the player's entertainment. Although MOBA is currently one of the most played genres around the world, it is known as a game that offer less autonomy, more challenges and consequently more frustration. Due to these characteristics, the use of a mechanism that performs the difficulty balance dynamically seems to be an interesting alternative to minimize and/or avoid that players experience such frustrations. In this sense, this paper presents a dynamic difficulty adjustment mechanism for MOBA games. The main idea is to create a computer controlled opponent that adapts dynamically to the player performance, trying to offer to the player a better game experience. This is done by evaluating the performance of the player using a metric based on some game features and switching the difficulty of the opponent's artificial intelligence behavior accordingly. Quantitative and qualitative experiments were performed and the results showed that the system is capable of adapting dynamically to the opponent's skills. In spite of that, the qualitative experiments with users showed that the player's expertise has a greater influence on the perception of the difficulty level and dynamic adaptation.


title: Topological phase transitions in finite-size periodically driven translationally invariant systems

abstract: It is known that, in the thermodynamic limit, the Chern number of a translationally invariant system cannot change under unitary time evolutions that are smooth in momentum space. Yet a real-space counterpart of the Chern number, the Bott index, has been shown to change in periodically driven systems with open boundary conditions. Here we prove that the Bott index and the Chern number are identical in translationally invariant systems in the thermodynamic limit. Using the Bott index, we show that, in finite-size translationally invariant systems, a Fermi sea under a periodic drive that is turned on slowly can acquire a different topology from that of the initial state. This can happen provided that the gap-closing points in the thermodynamic limit are absent in the discrete Brillouin zone of the finite system. Hence, in such systems, a periodic drive can be used to dynamically prepare topologically nontrivial states starting from topologically trivial ones.


title: Renewal-Theoretic Packet Collision Modeling under Long-Tailed Heterogeneous Traffic

abstract: Internet-of-things (IoT), with the vision of billions of connected devices, is bringing a massively heterogeneous character to wireless connectivity in unlicensed bands. The heterogeneity in medium access parameters, transmit power and activity levels among the coexisting networks leads to detrimental cross-technology interference. The stochastic traffic distributions, shaped under CSMA/CA rules, of an interfering network and channel fading makes it challenging to model and analyze the performance of an interfered network. In this paper, to study the temporal interaction between the traffic distributions of two coexisting networks, we develop a renewal-theoretic packet collision model and derive a generic collision-time distribution (CTD) function of an interfered system. The CTD function holds for any busy- and idle-time distributions of the coexisting traffic. As the earlier studies suggest a long-tailed idle-time statistics in real environments, the developed model only requires the Laplace transform of long-tailed distributions to find the CTD. Furthermore, we present a packet error rate (PER) model under the proposed CTD and multipath fading of the interfering signals. Using this model, a computationally efficient PER approximation for interference-limited case is developed to analyze the performance of an interfered link.


title: Mathematical base of Difference Operator of Particle Method

abstract: Mathematical base of difference operators in Moving Particle Semi-implicit method (MPS) are not given sufficiently in contrast to Smooth Particle Hydrodynamics method (SPH). Iribe and Nakaza proposed a method to improve the accuracy of the gradient operator, and Khayyer and Gotoh gave an ingenuity also for gradient operator too. An extension to higher order difference operators of Iribe-Nakaza method is given in this paper. The proposed method is a special case of the author's method called Discrete Differential Operators on Irregular Nodes (DDIN).


title: Delocalization in infinite disordered 2D lattices of different geometry

abstract: The spectral approach to infinite disordered crystals is applied to an Anderson-type Hamiltonian to demonstrate the existence of extended states for nonzero disorder in 2D lattices of different geometries. The numerical simulations shown prove that extended states exist for disordered honeycomb, triangular, and square crystals. This observation stands in contrast to the predictions of scaling theory, and aligns with experiments in photonic lattices and electron systems. The method used is the only theoretical approach aimed at showing delocalization. A comparison of the results for the three geometries indicates that the triangular and honeycomb lattices experience transition in the transport behavior for same amount of disorder, which is to be expected from planar duality. This provides justification for the use of artificially-prepared triangular lattices as analogues for honeycomb materials, such as graphene. The analysis also shows that the transition in the honeycomb case happens more abruptly as compared to the other two geometries, which can be attributed to the number of nearest neighbors. We outline the advantages of the spectral approach as a viable alternative to scaling theory and discuss its applicability to transport problems in both quantum and classical 2D systems.


title: Semipullbacks of labelled Markov processes

abstract: A labelled Markov process (LMP) consists of a measurable space $S$ together with an indexed family of Markov kernels from $S$ to itself. This structure has been used to model probabilistic computations in Computer Science, and one of the main problems in the area is to define and decide whether two LMP $S$ and $S'$ "behave the same". There are two natural categorical definitions of sameness of behavior: $S$ and $S'$ are bisimilar if there exist an LMP $ T$ and measure preserving maps forming a diagram of the shape $ S\leftarrow T \rightarrow{S'}$; and they are behaviorally equivalent if there exist some $ U$ and maps forming a dual diagram $ S\rightarrow U \leftarrow{S'}$. These two notions differ for general measurable spaces but Doberkat (extending a result by Edalat) proved that they coincide for analytic Borel spaces, showing that from every diagram $S\rightarrow U \leftarrow{S'}$ one can obtain a bisimilarity diagram as above. Moreover, the resulting square of measure preserving maps is commutative (a semipullback). In this paper, we extend the previous result to measurable spaces $S$ isomorphic to a universally measurable subset of a Polish space with the trace of the Borel $\sigma$-algebra, using a version of Strassen's theorem on common extensions of finitely additive measures.


title: Rigidity questions for real half-classical manifolds

abstract: Let $X$ be a noncommutative real compact algebraic manifold, in the sense that $C(X)=C^*(x_1,\ldots,x_N|x_i=x_i^*,f_\alpha(x_1,\ldots,x_N)=0)$, with $f_\alpha\in\mathbb R<x_1,\ldots,x_N>$. Associated to $X$ are its classical version $X^\times$, obtained via the relations $ab=ba$, and its half-classical version $X^*$, obtained via the relations $abc=cba$. We discuss here some general questions regarding the inclusions $X^\times\subset X^*\subset X$, and notably the comparison of the corresponding quantum isometry groups. Our main results concern the half-classical case, $X=X^*$, and more specifically, the case of the submanifolds $X\subset S^{N-1}_{\mathbb R,*}$.


title: Scalable Kernel K-Means Clustering with Nystrom Approximation: Relative-Error Bounds

abstract: Kernel $k$-means clustering can correctly identify and extract a far more varied collection of cluster structures than the linear $k$-means clustering algorithm. However, kernel $k$-means clustering is computationally expensive when the non-linear feature map is high-dimensional and there are many input points. Kernel approximation, e.g., the Nystr\"om method, has been applied in previous works to approximately solve kernel learning problems when both of the above conditions are present. This work analyzes the application of this paradigm to kernel $k$-means clustering, and shows that applying the linear $k$-means clustering algorithm to $\frac{k}{\epsilon} (1 + o(1))$ features constructed using a so-called rank-restricted Nystr\"om approximation results in cluster assignments that satisfy a $1 + \epsilon$ approximation ratio in terms of the kernel $k$-means cost function, relative to the guarantee provided by the same algorithm without the use of the Nystr\"om method. As part of the analysis, this work establishes a novel $1 + \epsilon$ relative-error trace norm guarantee for low-rank approximation using the rank-restricted Nystr\"om approximation. Empirical evaluations on the $8.1$ million instance MNIST8M dataset demonstrate the scalability and usefulness of kernel $k$-means clustering with Nystr\"om approximation. This work argues that spectral clustering using Nystr\"om approximation---a popular and computationally efficient, but theoretically unsound approach to non-linear clustering---should be replaced with the efficient and theoretically sound combination of kernel $k$-means clustering with Nystr\"om approximation. The superior performance of the latter approach is empirically verified.


title: Network Dimensions in the Getty Provenance Index

abstract: In this article we make a case for a systematic application of complex network science to study art market history and more general collection dynamics. We reveal social, temporal, spatial, and conceptual network dimensions, i.e. network node and link types, previously implicit in the Getty Provenance Index (GPI). As a pioneering art history database active since the 1980s, the GPI provides online access to source material relevant for research in the history of collecting and art markets. Based on a subset of the GPI, we characterize an aggregate of more than 267,000 sales transactions connected to roughly 22,000 actors in four countries over 20 years at daily resolution from 1801 to 1820. Striving towards a deeper understanding on multiple levels we disambiguate social dynamics of buying, brokering, and selling, while observing a general broadening of the market, where large collections are split into smaller lots. Temporally, we find annual market cycles that are shifted by country and obviously favor international exchange. Spatially, we differentiate near-monopolies from regions driven by competing sub-centers, while uncovering asymmetries of international market flux. Conceptually, we track dynamics of artist attribution that clearly behave like product categories in a very slow supermarket. Taken together, we introduce a number of meaningful network perspectives dealing with historical art auction data, beyond the analysis of social networks within a single market region. The results presented here have inspired a Linked Open Data conversion of the GPI, which is currently in process and will allow further analysis by a broad set of researchers.


title: Compressed Stop Searches with Two Leptons and Two b-jets

abstract: In top squark (stop) searches with a compressed spectrum, it is very helpful to consider the stop production recoiling against a hard jet from the initial state radiation to obtain a significant amount of missing transverse energy. In particular, the kinematic variable $R_M$ which measures the ratio of the lightest neutralino mass and the stop mass proved to be crucial in separating the signals from the backgrounds in both the all-hadronic decay and the semileptonic decay of the stops. Here we generalize the search method to the dileptonic stop decays. In this case, due to the two missing neutrinos, there are not enough kinematic constraint equations to solve for the $R_M$ variable exactly, but only render an allowed interval consistent with the event. However, we show that the minimum and the maximum values of this allowed interval still provide useful variables in discriminating signals from the backgrounds. Although in the traditional stop decay to a top quark and the lightest neutralino, the dileptonic mode is not as competitive due to its small branching ratio, it becomes the main search mode if the stops decay through the charginos and sleptons with a compressed spectrum. We show that with the new variables, the dileptonic search of the stop can cover regions of the parameter space which have not been constrained before.


title: Controlling Spin-flips of Molecules in an Electromagnetic Trap

abstract: Doubly dipolar molecules exhibit complex internal spin-dynamics when electric and magnetic fields are both applied. Near magnetic trap minima, these spin-dynamics lead to enhancements in Majorana spin-flip transitions by many orders of magnitude relative to atoms, and are thus an important obstacle for progress in molecule trapping and cooling. We conclusively demonstrate and address this with OH molecules in a trap geometry where spin-flip losses can be tuned from over $200 \text{ s}^{-1} $ to below our $2\text{ s}^{-1}$ vacuum limited loss rate with only a simple external bias coil and with minimal impact on trap depth and gradient.


title: Learning to Embed Words in Context for Syntactic Tasks

abstract: We present models for embedding words in the context of surrounding words. Such models, which we refer to as token embeddings, represent the characteristics of a word that are specific to a given context, such as word sense, syntactic category, and semantic role. We explore simple, efficient token embedding models based on standard neural network architectures. We learn token embeddings on a large amount of unannotated text and evaluate them as features for part-of-speech taggers and dependency parsers trained on much smaller amounts of annotated data. We find that predictors endowed with token embeddings consistently outperform baseline predictors across a range of context window and training set sizes.


title: A randomized Halton algorithm in R

abstract: Randomized quasi-Monte Carlo (RQMC) sampling can bring orders of magnitude reduction in variance compared to plain Monte Carlo (MC) sampling. The extent of the efficiency gain varies from problem to problem and can be hard to predict. This article presents an R function rhalton that produces scrambled versions of Halton sequences. On some problems it brings efficiency gains of several thousand fold. On other problems, the efficiency gain is minor. The code is designed to make it easy to determine whether a given integrand will benefit from RQMC sampling. An RQMC sample of n points in $[0,1]^d$ can be extended later to a larger n and/or d.


title: Asymptotic behaviors in the homology of symmetric group and finite general linear group quandles

abstract: A quandle is an algebraic structure which attempts to generalize group conjugation. These structures have been studied extensively due to their connections with knot theory, algebraic combinatorics, and other fields. In this work, we approach the study of quandles from the perspective of the representation theory of categories. Namely, we look at collections of conjugacy classes of the symmetric groups and the finite general linear groups, and prove that they carry the structure of FI-quandles (resp. VIC$(q)$-quandles). As applications, we prove statements about the homology of these quandles, and construct FI-module and VIC$(q)$-module invariants of links.


title: A Discontinuous Differential Calculus in the Framework Colombeau's Full Algebra

abstract: Starting from the Colombeau's full generalized functions, the sharp topologies and the notion of generalized points, we introduce a new kind differential calculus (for functions between totally disconnected spaces). We study generalized pointvalues, Colombeau's differential algebra, holomorphic and analytic functions. We show that the Embedding Theorem and the Open Mapping Theorem hold in this framework. Moreover, we study some applications in differential equations.


title: Symmetric Contours and Convergent Interpolation

abstract: The essence of Stahl-Gonchar-Rakhmanov theory of symmetric contours as applied to the multipoint Pad\'e approximants is the fact that given a germ of an algebraic function and a sequence of rational interpolants with free poles of the germ, if there exists a contour that is "symmetric" with respect to the interpolation scheme, does not separate the plane, and in the complement of which the germ has a single-valued continuation with non-identically zero jump across the contour, then the interpolants converge to that continuation in logarithmic capacity in the complement of the contour. The existence of such a contour is not guaranteed. In this work we do construct a class of pairs interpolation scheme/symmetric contour with the help of hyperelliptic Riemann surfaces (following the ideas of Nuttall \& Singh and Baratchart \& the author. We consider rational interpolants with free poles of Cauchy transforms of non-vanishing complex densities on such contours under mild smoothness assumptions on the density. We utilize \( \bar\partial \)-extension of the Riemann-Hilbert technique to obtain formulae of strong asymptotics for the error of interpolation.