title: Fractal iso-level sets in high-Reynolds-number scalar turbulence

abstract: We study the fractal scaling of iso-levels sets of a passive scalar mixed by three-dimensional homogeneous and isotropic turbulence at high Reynolds numbers. The Schmidt number is unity. A fractal box-counting dimension $D_F$ can be obtained for iso-levels below about $3$ standard deviations of the scalar fluctuation on either side of its mean value. The dimension varies systematically with the iso-level, with a maximum of about $8/3$ for the iso-level at the mean; this maximum dimension also follows as an upper bound from the geometric measure theory. We interpret this result to mean that mixing in turbulence is always incomplete. A unique box-counting dimension for all iso-levels results when we consider the spatial support of the steep cliffs of the scalar conditioned on local strain; that unique dimension is about $4/3$.


title: Set-based MPC for discrete-time LTI systems with maximal domain of attraction and minimal predictive control horizon

abstract: This paper presents a novel set-based model predictive control for tracking, which provides the largest domain of attraction, even with the minimal predictive/control horizon. The formulation - which consists of a single optimization problem - shows a dual behavior: one operating inside the maximal controllable set to the feasible equilibrium set, and the other operating at the N-controllable set to the same equilibrium set. Based on some finite-time convergence results, asymptotic stability of the resulting closed-loop is proved, while recursive feasibility is ensured for any change of the setpoint. The properties and advantages of the proposal have been tested on simulation models.


title: Tracking Dust Grains During Transport and Growth in Protoplanetary Disks

abstract: Protoplanetary disks are dynamic objects, within which dust grains and gas are expected to be redistributed over large distances. Evidence for this redistribution is seen both in other protoplanetary disks and in our own Solar System, with high-temperature materials thought to originate close to the central star found in the cold, outer regions of the disks. While models have shown this redistribution is possible through a variety of mechanisms, these models have generally ignored the possible growth of solids via grain-grain collisions that would occur during transit. Here we investigate the interplay of coagulation and radial and vertical transport of solids in protoplanetary disks, considering cases where growth is limited by bouncing or by fragmentation. We find that in all cases, growth effectively limits the ability for materials to be carried outward or preserved at large distances from the star. This is due to solids being incorporated into large aggregates which drift inwards rapidly under the effects of gas drag. We discuss the implications for mixing in protoplanetary disks, and how the preservation of high temperature materials in outer disks may require structures or outward flow patterns to avoid them being lost via radial drift.


title: DyKgChat: Benchmarking Dialogue Generation Grounding on Dynamic Knowledge Graphs

abstract: Data-driven, knowledge-grounded neural conversation models are capable of generating more informative responses. However, these models have not yet demonstrated that they can zero-shot adapt to updated, unseen knowledge graphs. This paper proposes a new task about how to apply dynamic knowledge graphs in neural conversation model and presents a novel TV series conversation corpus (DyKgChat) for the task. Our new task and corpus aids in understanding the influence of dynamic knowledge graphs on responses generation. Also, we propose a preliminary model that selects an output from two networks at each time step: a sequence-to-sequence model (Seq2Seq) and a multi-hop reasoning model, in order to support dynamic knowledge graphs. To benchmark this new task and evaluate the capability of adaptation, we introduce several evaluation metrics and the experiments show that our proposed approach outperforms previous knowledge-grounded conversation models. The proposed corpus and model can motivate the future research directions.


title: The hadronic vacuum polarization contribution to $(g_\mu-2)$: Lattice QCD+QED calculations

abstract: The anomalous magnetic moment of the muon $a_\mu$ is one of the most accurate quantities in Particle Physics. The long-standing discrepancy of about $3.7$ standard deviations between the experimental value and the prediction of the Standard Model could represent an intriguing indication of New Physics. The experiments at Fermilab (E989) and at J-PARC (E34) aim at reducing significantly the experimental uncertainty, thus making the theoretical one due to hadronic corrections the main limitation of this stringent test of the Standard Model. In this contribution we present the results of a first-principles lattice calculation of the hadronic vacuum polarization (HVP) contribution to $a_\mu$, including electromagnetic and $SU(2)$-breaking corrections. Our determination, $a_\mu^{\rm HVP} = 682.0 ~ (18.7) \cdot 10^{-10}$, turns out to be in agreement with recent theoretical determinations based on the dispersive analyses of the experimental cross section data for the annihilation process $e^+e^- \to hadrons$. Furthermore, we provide for the first time a lattice estimate for the missing part of $a_\mu^{\rm HVP}$ not covered in the MUonE experiment, $\left[a_\mu^{\rm HVP}\right]_> = 91.6 ~ (2.0) \cdot 10^{-10}$.


title: Multicomponent dark matter in extended $U(1)_{B-L}$: neutrino mass and high scale validity

abstract: Standard Model with right handed neutrinos charged under additional $U(1)_{B-L}$ gauge symmetry offer solutions to both dark matter (DM) problem and neutrino mass generation, although constrained severely from relic density, direct search and Higgs vacuum stability. We therefore investigate a multicomponent DM scenario augmented by an extra inert scalar doublet, that is neutral under $U(1)_{B-L}$, which aids to enlarge parameter space allowed by DM constraints and Higgs vacuum stability. The lightest right-handed neutrino and the $CP$-even inert scalar are taken as the dark matter candidates and constitute a two component dark matter framework as they are rendered stable by an unbroken $\mathbb{Z}_2 \times \mathbb{Z}_2^\prime$ symmetry. DM-DM conversion processes turn out crucial to render requisite relic abundance in mass regions of the RH neutrino that do not appear in the stand-alone $U(1)_{B-L}$ scenario. In addition, the one-loop renormalisation group (RG) equations in this model demonstrate that the electroweak (EW) vacuum can be stabilised till $\sim 10^{9}$ GeV in a parameter region compatible with the observed relic, the direct detection bound and other relevant constraints.


title: A family picture: tracing the dynamical path of the structural properties of multiple populations in globular clusters

abstract: We studied the spatial distributions of multiple stellar populations (MPs) in a sample of 20 globular clusters (GCs) spanning a broad range of dynamical ages. The differences between first-population (FP) and second-population (SP) stars were measured by means of the parameter $A^+$, defined as the area enclosed between their cumulative radial distributions. We provide the first purely observational evidence of the dynamical path followed by MPs from initial conditions toward a complete FP-SP spatial mixing. Less dynamically evolved clusters have SP stars more centrally concentrated than FPs, while in more dynamically evolved systems the spatial differences between FP and SP stars decrease and eventually disappear. By means of an appropriate comparison with a set of numerical simulations, we show that these observational results are consistent with the evolutionary sequence expected by the long-term dynamical evolution of clusters forming with an initially more centrally concentrated SP sub-system. This result is further supported by the evidence of a trend between $A^+$ and the stage of GC dynamical evolution inferred by the ratio between the present-day and the initial mass of the cluster.


title: The Choice Function Framework for Online Policy Improvement

abstract: There are notable examples of online search improving over hand-coded or learned policies (e.g. AlphaZero) for sequential decision making. It is not clear, however, whether or not policy improvement is guaranteed for many of these approaches, even when given a perfect evaluation function and transition model. Indeed, simple counter examples show that seemingly reasonable online search procedures can hurt performance compared to the original policy. To address this issue, we introduce the choice function framework for analyzing online search procedures for policy improvement. A choice function specifies the actions to be considered at every node of a search tree, with all other actions being pruned. Our main contribution is to give sufficient conditions for stationary and non-stationary choice functions to guarantee that the value achieved by online search is no worse than the original policy. In addition, we describe a general parametric class of choice functions that satisfy those conditions and present an illustrative use case of the framework's empirical utility.


title: A translation of Henri Joris' "Le chasseur perdu dans la for\^et" (1980)

abstract: This is an English translation of Henri Joris' article "Le chasseur perdu dans la for\^et (Un probl\`eme de g\'eom\'etrie plane)" that appeared in $Elemente\ der\ Mathematik$ v. 35 (1980) n. 1, 1--14. Given a point $P$ and a line $L$ in the plane, what is the shortest search path to find $L$, given its distance but not its direction from $P$? The shortest search path was described by Isbell (1957), but a complete and detailed proof was not published until Joris (1980). I am thankful to Natalya Pluzhnikov for her dedicated work and to the Swiss Mathematical Society for permission to post this translation on the arXiv.


title: Giant spontaneous exchange bias in an antiperovskite structure driven by a canted triangular magnetic structure

abstract: Exchange bias (EB) refers to a shift of the hysteresis loop along the field axis in materials consisting of ferromagnetic (FM) and antiferromagnetic (AFM) layers, generally after a cooling procedure in high magnetic field. This effect is highly desirable for technological applications ranging from spintronics to magnetic recording. Achieving giant EB effect near room temperature in a small cooling field is thus an on-going technologically relevant challenge for the materials science community. In this work, we present the experimental realization of such a fundamental goal by demonstrating the existence of giant EB after a zero field cooled (ZFC) procedure in antiperovskite Mn3.5Co0.5N below 256 K. We found that it exhibits an EB field of 0.28 T at 50 K after a ZFC procedure accompanied by a large vertical magnetization shift (VMS). Interestingly, this EB field can be elevated up to 1.2 T after a cooling procedure with a small applied field of just 500 Oe. Mn3:5Co0:5N bears the first intermetallic material showing a strong correlation between EB and VMS.We attribute the observed EB effect to a completely new canted triangular magnetic structure determined by neutron diffraction experiment. Finally, we discuss the striking effect of Co substitution on the physical properties of antiperovskite materials and put forward a new strategy for antiperovskite lattice to exploit new single phase materials showing large EB effect at room temperature.


title: Predicting materials properties without crystal structure: Deep representation learning from stoichiometry

abstract: Machine learning has the potential to accelerate materials discovery by accurately predicting materials properties at a low computational cost. However, the model inputs remain a key stumbling block. Current methods typically use descriptors constructed from knowledge of either the full crystal structure -- therefore only applicable to materials with already characterised structures -- or structure-agnostic fixed-length representations hand-engineered from the stoichiometry. We develop a machine learning approach that takes only the stoichiometry as input and automatically learns appropriate and systematically improvable descriptors from data. Our key insight is to treat the stoichiometric formula as a dense weighted graph between elements. Compared to the state of the art for structure-agnostic methods, our approach achieves lower errors with less data.


title: Omnipush: accurate, diverse, real-world dataset of pushing dynamics with RGB-D video

abstract: Pushing is a fundamental robotic skill. Existing work has shown how to exploit models of pushing to achieve a variety of tasks, including grasping under uncertainty, in-hand manipulation and clearing clutter. Such models, however, are approximate, which limits their applicability. Learning-based methods can reason directly from raw sensory data with accuracy, and have the potential to generalize to a wider diversity of scenarios. However, developing and testing such methods requires rich-enough datasets. In this paper we introduce Omnipush, a dataset with high variety of planar pushing behavior. In particular, we provide 250 pushes for each of 250 objects, all recorded with RGB-D and a high precision tracking system. The objects are constructed so as to systematically explore key factors that affect pushing -- the shape of the object and its mass distribution -- which have not been broadly explored in previous datasets, and allow to study generalization in model learning. Omnipush includes a benchmark for meta-learning dynamic models, which requires algorithms that make good predictions and estimate their own uncertainty. We also provide an RGB video prediction benchmark and propose other relevant tasks that can be suited with this dataset. Data and code are available at \url{https://web.mit.edu/mcube/omnipush-dataset/}.


title: The energy budget and figure of Earth during recovery from the Moon-forming giant impact

abstract: Quantifying the energy budget of Earth in the first few million years following the Moon-forming giant impact is vital to understanding Earth's initial thermal state and the dynamics of lunar tidal evolution. After the impact, the body was substantially vaporized and rotating rapidly, very different from the planet we know today. The subsequent evolution of Earth's energy budget, as the body cooled and angular momentum was transferred during lunar tidal recession, has not been accurately calculated with all relevant energy components included. Here, we use giant impact simulations and planetary structure models to calculate the energy budget at stages in Earth's evolution. We show that the figure and internal structure of Earth changed substantially during its post-impact evolution and that changes in kinetic, potential, and internal energy were all significant. These changes have important implications for the dynamics of tidal recession and the thermal structure of early Earth.


title: Holographic shear rheology of viscoelastic fluids

abstract: In this study, we report the use of digital holography microscopy (DHM) for 3D-resolved flow kinematics and shear rheometry of viscoelastic polymeric fluids. We computationally reconstruct the recorded holograms to visualize the tracer imbued flow volume in microchannels, followed by implementation of particle tracking velocimetry (PTV) to quantitate spatially-resolved velocity fields in 3D. In order to select optimal parameters for DHM-PTV characterization of complex fluids, we studied the effect of hologram recording distance, seeding density and particle size. Using the optimal parameters, we show quantitative characterization of the shear rheology from the velocity fields without any a-priori assumptions of wall boundary condition or constitutive equation. The viscosity versus shear rate data for Newtonian and polyethylene oxide solutions could be measured in the range of ~ 0.05 - 20,000 s-1 with just four input flow rates. This data from holographic shear rheometry was found to be in good agreement with computational fluid dynamics simulations and macrorheometry. The holographic shear rheology technique remained unaffected by wall-slip events and instead provided an avenue to quantitate slip severity. Finally, we discuss holographic visualization of particle migration in microfluidic flows which can limit flow field access while at the same time provide a fingerprint of the suspending fluid rheology.


title: Geometric percolation of hard nanorods: the interplay of spontaneous and externally induced uniaxial particle alignment

abstract: We present a numerical study on geometric percolation in liquid dispersions of hard slender colloidal particles subjected to an external orienting field. In the formulation and liquid-state processing of nanocomposite materials, the alignment of particles by external fields such as electric, magnetic or flow fields is practically inevitable, and often works against the emergence of large nanoparticle networks. Using continuum percolation theory in conjunction with Onsager theory, we investigate how the interplay between externally induced alignment and the spontaneous symmetry breaking of the uniaxial nematic phase affects cluster formation within nanoparticle dispersions. It is known that the enhancement of particle alignment by means of a density increase or an external field may result in the breakdown of an already percolating network. As a result, percolation can be limited to a small region of the phase diagram only. Here, we demonstrate that the existence and shape of such a "percolation island" in the phase diagram crucially depends on the connectivity length -- a critical distance defining direct connections between neighbouring particles. Deformations of this percolation island can lead to peculiar re-entrance effects, in which a system-spanning network forms and breaks down multiple times with increasing particle density.

categories: cond-mat.soft cond-mat.mtrl-sci cond-mat.stat-mech physics.app-ph


title: Generalized Non-Minimal Derivative Coupling: Application to Inflation and Primordial Black Hole Production

abstract: We introduce inflationary models where the inflaton features a field dependent non-minimal derivative coupling to the Einstein tensor, that we name GNMDC. This Horndeski term gives new and distinguishable inflationary predictions in a framework that ameliorates possible problems related with gradient instabilities during the reheating stage. We examine the inflationary phenomenology using power law monomial and exponential potentials. We further elaborate on the implications of the inflaton-modulated GNMDC and construct GNMDC terms that can amplify the power spectrum of primordial perturbations at small scales, triggering PBH production. An attractive feature of the GNMDC is that inflation, as well as PBH production, can be implemented utilizing the Higgs potential.


title: A Multiple Filter Based Neural Network Approach to the Extrapolation of Adsorption Energies on Metal Surfaces for Catalysis Applications

abstract: Computational catalyst discovery involves the development of microkinetic reactor models based on estimated parameters determined from density functional theory (DFT). For complex surface chemistries, the cost of calculating the adsorption energies by DFT for a large number of reaction intermediates can become prohibitive. Here, we have identified appropriate descriptors and machine learning models that can be used to predict part of these adsorption energies given data on the rest of them. Our investigations also included the case when the species data used to train the predictive model is of different size relative to the species the model tries to predict - an extrapolation in the data space which is typically difficult with regular machine learning models. We have developed a neural network based predictive model that combines an established model with the concepts of a convolutional neural network that, when extrapolating, achieves significant improvement over the previous models.


title: Discrete Laplacian in a half-space with a periodic surface potential I: Resolvent expansions, scattering matrix, and wave operators

abstract: We present a detailed study of the scattering system given by the Neumann Laplacian on the discrete half-space perturbed by a periodic potential at the boundary. We derive asymptotic resolvent expansions at thresholds and eigenvalues, we prove the continuity of the scattering matrix, and we establish new formulas for the wave operators. Along the way, our analysis puts into evidence a surprising relation between some properties of the potential, like the parity of its period, and the behaviour of the integral kernel of the wave operators.


title: Sub-nanosecond spin-torque switching of perpendicular magnetic tunnel junction nanopillars at cryogenic temperatures

abstract: Spin-transfer magnetic random access memory is of significant interest for cryogenic applications where a persistent, fast, low-energy consumption and high device density is needed. Here we report the low-temperature nanosecond duration spin-transfer switching characteristics of perpendicular magnetic tunnel junction (pMTJ) nanopillar devices (40 to 60 nm in diameter) and contrast them to their room temperature properties. Interestingly, at fixed pulse voltage overdrive the characteristic switching time decreases with temperature, in contrast to macrospin model predictions, with the largest reduction in switching time occurring between room temperature and 150 K. The switching energy increases with decreasing temperature, but still compares very favorably to other types of spin-transfer devices at 4 K, with < 300 fJ required per switch. Write error rate (WER) measurements show highly reliable (WER <= 5x10^-5 with 4 ns pulses at 4 K) demonstrating the promise of pMTJ devices for cryogenic applications and routes to further device optimization.


title: Pre- and post-processing in quantum-computational hydrologic inverse analysis

abstract: It was recently shown that certain subsurface hydrological inverse problems -- here framed as determining the composition of an aquifer from pressure readings -- can be solved on a quantum annealer. However, the quantum annealer performance suffered when solving problems where the aquifer was composed of materials with vastly different permeability, which is often encountered in practice. In this paper, we study why this regime is difficult and use several pre- and post-processing tools to address these issues. This study has three benefits: it improves quantum annealing performance for real-world problems in hydrology, it studies the scaling behavior for these problems (which were previously studied at a fixed size), and it elucidates a challenging class of problems that are amenable to quantum annealers.


title: Tropical Moduli Space of Rational Graphically Stable Curves

abstract: We study moduli spaces of rational graphically stable tropical curves and a refinement given by radial alignment. Given a complete multipartite graph $\Gamma$, the moduli space of radially aligned $\Gamma$-stable tropical curves can be given the structure of a balanced fan. This fan structure coincides with the Bergman fan of the cycle matroid of $\Gamma$.


title: Sensor Fusion: Gated Recurrent Fusion to Learn Driving Behavior from Temporal Multimodal Data

abstract: The Tactical Driver Behavior modeling problem requires understanding of driver actions in complicated urban scenarios from a rich multi modal signals including video, LiDAR and CAN bus data streams. However, the majority of deep learning research is focused either on learning the vehicle/environment state (sensor fusion) or the driver policy (from temporal data), but not both. Learning both tasks end-to-end offers the richest distillation of knowledge, but presents challenges in formulation and successful training. In this work, we propose promising first steps in this direction. Inspired by the gating mechanisms in LSTM, we propose gated recurrent fusion units (GRFU) that learn fusion weighting and temporal weighting simultaneously. We demonstrate it's superior performance over multimodal and temporal baselines in supervised regression and classification tasks, all in the realm of autonomous navigation. We note a 10% improvement in the mAP score over state-of-the-art for tactical driver behavior classification in HDD dataset and a 20% drop in overall Mean squared error for steering action regression on TORCS dataset.


title: Uniqueness of Optimal Point Sets Determining Two Distinct Triangles

abstract: In this paper, we show that the maximum number of points in $d\geq3$ dimensions determining exactly 2 distinct triangles is $2d$. We further show that this maximum is uniquely achieved by the vertices of the $d$-orthoplex. We build upon the work of Hirasaka and Shinohara who determined that the $d$-orthoplex is such an optimal configuration, but did not prove its uniqueness. Further, we present a more elementary argument for its optimality.


title: COHERENT Plans for D$_2$O at the Spallation Neutron Source

abstract: The Spallation Neutron Source (SNS) is a pulsed source of neutrons and, as a byproduct of this operation, an intense source of neutrinos via stopped-pion decay. The COHERENT collaboration uses this source to investigate coherent elastic neutrino-nucleus scattering (CEvNS) with a suite of detectors. To enable precise cross-section measurements, we must address an estimated 10\% uncertainty in our flux calculation associated with the lack of data for $\pi^\pm$ production from 1 GeV protons on an Hg target. We present here our Geant4 simulation of neutrino production at the SNS and our plans to experimentally normalize this flux with the development of a 670 kg D$_2$O detector. Using the precise cross section calculations for neutrino interactions on deuterium, we will dramatically reduce our flux uncertainty.


title: The correlation functions of certain random antiferromagnetic spin-1/2 critical chains

abstract: We study the spin-spin correlations in two distinct random critical XX spin-1/2 chain models via exact diagonalization. For the well-known case of uncorrelated random coupling constants, we study the non-universal numerical prefactors and relate them to the corresponding Lyapunov exponent of the underlying single-parameter scaling theory. We have also obtained the functional form of the correct scaling variables important for describing even the strongest finite-size effects. Finally, with respect to the distribution of the correlations, we have numerically determined that they converge to a universal (disorder-independent) non-trivial and narrow distribution when properly rescaled by the spin-spin separation distance in units of the Lyapunov exponent. With respect to the less known case of correlated coupling constants, we have determined the corresponding exponents and shown that both typical and mean correlations functions decay algebraically with the distance. While the exponents of the transverse typical and mean correlations are nearly equal, implying a narrow distribution of transverse correlations, the longitudinal typical and mean correlations critical exponents are quite distinct implying much broader distributions. Further comparisons between these models are given.


title: Kinematics of Hot Subdwarfs from the Gaia DR2 Catalogue

abstract: We have studied the kinematic properties of the candidates for hot subdwarfs (HSDs) selected by Geier et al. from the Gaia DR2 catalogue. We have used a total of 12 515 stars with relative trigonometric parallax errors less than 30\%. The HSDs are shown to have different kinematics, depending on their positions on the celestial sphere. For example, the sample of low-latitude $(|b|<20^\circ)$ HSDs rotates around the Galactic center with a linear velocity $V_0=221\pm5$ km s$^{-1}$. This suggests that they belong to the Galactic thin disk. At the same time, they lag behind the local standard of rest by $\Delta V_\odot\sim16$ km s$^{-1}$ due to the asymmetric drift. The high-latitude ($|b|\geq20^\circ$) HSDs rotate with a considerably lower velocity, $V_0=168\pm6$ km s$^{-1}$. Their lagging behind the local standard of rest is already $\Delta V_\odot\sim40$ km s$^{-1}$. Based on the entire sample of 12 515 HSDs, we have found a positive rotation around the $x$ axis significantly differing from zero with an angular velocity $\omega_1=1.36\pm0.24$ km s$^{-1}$ kpc$^{-1}$. We have studied the samples of HSDs that are complete within $r<1.5$ kpc. Based on them, we have traced the evolution of the parameters of the residual velocity ellipsoid as a function of both latitude $|b|$ and coordinate $|z|.$ The following vertical disk scale heights have been found: $h=180\pm6$ and $290\pm10$ pc from the low- and high-latitude HSDs, respectively. A new estimate of the local stellar density $\Sigma_{out}=53\pm4~M_{\odot}$/kpc$^2$ kpc$^{-2}$ has been obtained for $z_{out}=0.56$ kpc from the high-latitude HSDs.


title: Characterizing optimal point sets determining one distinct triangle

abstract: In this paper we determine the maximum number of points in $\mathbb{R}^d$ which form exactly $t$ distinct triangles, where we restrict ourselves to the case of $t = 1$. We denote this quantity by $F_d(t)$. It was known from the work of Epstein et al. that $F_2(1) = 4$. Here we show somewhat surprisingly that $F_3(1) = 4$ and $F_d(1) = d + 1$, whenever $d \geq 3$, and characterize the optimal point configurations. This is an extension of a variant of the distinct distance problem put forward by Erd\H{o}s and Fishburn.


title: Vortex-loop calculation of the specific heat of superfluid $^4$He under pressure

abstract: Vortex-loop renormalization is used to compute the specific heat of superfluid $^4$He near the lambda point at various pressures up to 26 bars. The input parameters are the the pressure dependence of T$_\lambda$ and the superfluid density, which determine the non-universal parameters of the vortex core energy and core size. The results for the specific heat are found to be in good agreement with experimental data, matching the expected universal pressure dependence to within about 5$\%$. The non-universal critical amplitude of the specific heat is found to be in reasonable agreement, a factor of four larger than the experiments. We point out problems with recent Gross-Pitaevskii simulations that claimed the vortex-loop percolation temperature did not match the critical temperature of the superfluid phase transition.


title: Extraction of Efficient Programs in $I\Sigma_1$-arithmetic

abstract: Clausal Language (CL) is a declarative programming and verifying system used in our teaching of computer science. CL is an implementation of, what we call, $\mathit{PR}{+}I\Sigma_1$ paradigm (primitive recursive functions with $I\Sigma_1$-arithmetic). This paper introduces an extension of $I\Sigma_1$-proofs called extraction proofs where one can extract from the proofs of $\Pi_2$-specifications primitive recursive programs as efficient as the hand-coded ones. This is achieved by having the programming constructs correspond exactly to the proof rules with the computational content.


title: Anomalous thermal expansion in one-dimensional transition-metal cyanides: Behavior of the trimetallic cyanide Cu$_{1/3}$Ag$_{1/3}$Au$_{1/3}$CN

abstract: The structural dynamics of a 1D mixed-metal cyanide, Cu$_{1/3}$Ag$_{1/3}$Au$_{1/3}$CN, with intriguing thermal properties is explored. All the current known related compounds with straight-chain structures, such as CuCN, AgCN, AuCN and M$_x$M'$_{1-x}$CN (M, M' = Cu, Ag, Au), exhibit 1D negative thermal expansion (NTE) along the chains and positive thermal expansion (PTE) perpendicular to them. Cu$_{1/3}$Ag$_{1/3}$Au$_{1/3}$CN exhibits similar PTE perpendicular to the chains, however PTE, rather than NTE, is also observed along the chains. In order to understand the origin of this unexpected behavior, inelastic neutron scattering (INS) measurements were carried out, underpinned by DFT calculations. Synchrotron-based PDF analysis and 13C SSNMR measurements were also performed to build an input structural model for the lattice dynamical study. The results indicate that transverse motions of the metal ions are responsible for the PTE perpendicular to the chains, as is the case for the related group 11 cyanides. However NTE along the chain due to the tension effect of these transverse motions is not observed. For this system, DFT-based phonon calculations predict small PTE along the chain due to low-energy chain-slipping modes induced by a bond-rotation effect on the weak metallophilic bonds. However the observed PTE is greater than that predicted with the straight-chain model. Small bends in the chain to accommodate truly co-planar metals provide an alternative explanation for thermal behavior. These would mitigate the tension effect induced by the transverse motions of the metals and, as temperature increases and the chains move further apart, a straightening could occur resulting in the observed PTE. This hypothesis is further supported by unusual evolution in the phonon spectra, which suggest small changes in local symmetry with temperature.


title: Essentia: Mining Domain-Specific Paraphrases with Word-Alignment Graphs

abstract: Paraphrases are important linguistic resources for a wide variety of NLP applications. Many techniques for automatic paraphrase mining from general corpora have been proposed. While these techniques are successful at discovering generic paraphrases, they often fail to identify domain-specific paraphrases (e.g., {staff, concierge} in the hospitality domain). This is because current techniques are often based on statistical methods, while domain-specific corpora are too small to fit statistical methods. In this paper, we present an unsupervised graph-based technique to mine paraphrases from a small set of sentences that roughly share the same topic or intent. Our system, Essentia, relies on word-alignment techniques to create a word-alignment graph that merges and organizes tokens from input sentences. The resulting graph is then used to generate candidate paraphrases. We demonstrate that our system obtains high-quality paraphrases, as evaluated by crowd workers. We further show that the majority of the identified paraphrases are domain-specific and thus complement existing paraphrase databases.


title: New symmetry-adapted ab initio approach to nuclear reactions for intermediate-mass nuclei

abstract: With a view toward describing reactions of intermediate-mass nuclei from first principles, we present first results for the norm and Hamiltonian overlaps (kernels) for the p-{\alpha}, p-16O and p-20Ne cluster systems using realistic nucleon-nucleon interactions. This is achieved in the framework of a new ab initio approach that combines the symmetry-adapted no-core shell model (SA-NCSM) with the resonating group method (RGM). In this model, a physically relevant basis based on the SU(3) symmetry is used. The structure of the clusters is provided by the ab initio SA-NCSM, which enables the description of spatially enhanced nuclear configurations and heavier nuclei, by exploiting symmetries known to dominate in nuclei. Here, we discuss the applicability and efficacy of this approach.


title: Ancient asymptotically cylindrical flows and applications

abstract: In this paper, we prove the mean-convex neighborhood conjecture for neck singularities of the mean curvature flow in $\mathbb{R}^{n+1}$ for all $n\geq 3$: we show that if a mean curvature flow $\{M_t\}$ in $\mathbb{R}^{n+1}$ has an $S^{n-1}\times \mathbb{R}$ singularity at $(x_0,t_0)$, then there exists an $\varepsilon=\varepsilon(x_0,t_0)>0$ such that $M_t\cap B(x_0,\varepsilon)$ is mean-convex for all $t\in(t_0-\varepsilon^2,t_0+\varepsilon^2)$. As in the case $n=2$, which was resolved by the first three authors in arXiv:1810.08467, the existence of such a mean-convex neighborhood follows from classifying a certain class of ancient Brakke flows that arise as potential blowup limits near a neck singularity. Specifically, we prove that any ancient unit-regular integral Brakke flow with a cylindrical blowdown must be either a round shrinking cylinder, a translating bowl soliton, or an ancient oval. In particular, combined with a prior result of the last two authors, we obtain uniqueness of mean curvature flow through neck singularities. The main difficulty in addressing the higher dimensional case is in promoting the spectral analysis on the cylinder to global geometric properties of the solution. Most crucially, due to the potential wide variety of self-shrinking flows with entropy lower than the cylinder when $n\geq 3$, smoothness does not follow from the spectral analysis by soft arguments. This precludes the use of the classical moving plane method to derive symmetry. To overcome this, we introduce a novel variant of the moving plane method, which we call "moving plane method without assuming smoothness" - where smoothness and symmetry are established in tandem.


title: On the Concavity of Expected Shortfall

abstract: It is well known that Expected Shortfall (also called Average Value-at-Risk) is a convex risk measure, i. e. Expected Shortfall of a convex linear combination of arbitrary risk positions is not greater than a convex linear combination with the same weights of Expected Shortfalls of the same risk positions. In this short paper we prove that Expected Shortfall is a concave risk measure with respect to probability distributions, i. e. Expected Shortfall of a finite mixture of arbitrary risk positions is not lower than the linear combination of Expected Shortfalls of the same risk positions (with the same weights as in the mixture).


title: An introduction to flexible methods for policy evaluation

abstract: This chapter covers different approaches to policy evaluation for assessing the causal effect of a treatment or intervention on an outcome of interest. As an introduction to causal inference, the discussion starts with the experimental evaluation of a randomized treatment. It then reviews evaluation methods based on selection on observables (assuming a quasi-random treatment given observed covariates), instrumental variables (inducing a quasi-random shift in the treatment), difference-in-differences and changes-in-changes (exploiting changes in outcomes over time), as well as regression discontinuities and kinks (using changes in the treatment assignment at some threshold of a running variable). The chapter discusses methods particularly suited for data with many observations for a flexible (i.e. semi- or nonparametric) modeling of treatment effects, and/or many (i.e. high dimensional) observed covariates by applying machine learning to select and control for covariates in a data-driven way. This is not only useful for tackling confounding by controlling for instance for factors jointly affecting the treatment and the outcome, but also for learning effect heterogeneities across subgroups defined upon observable covariates and optimally targeting those groups for which the treatment is most effective.


title: Even $A$-cycles have the edge-Erd\H{o}s-P\'osa property

abstract: I prove that even $A$-cycles have the edge-Erd\H{o}s-P\'osa property.


title: SlowMo: Improving Communication-Efficient Distributed SGD with Slow Momentum

abstract: Distributed optimization is essential for training large models on large datasets. Multiple approaches have been proposed to reduce the communication overhead in distributed training, such as synchronizing only after performing multiple local SGD steps, and decentralized methods (e.g., using gossip algorithms) to decouple communications among workers. Although these methods run faster than AllReduce-based methods, which use blocking communication before every update, the resulting models may be less accurate after the same number of updates. Inspired by the BMUF method of Chen & Huo (2016), we propose a slow momentum (SlowMo) framework, where workers periodically synchronize and perform a momentum update, after multiple iterations of a base optimization algorithm. Experiments on image classification and machine translation tasks demonstrate that SlowMo consistently yields improvements in optimization and generalization performance relative to the base optimizer, even when the additional overhead is amortized over many updates so that the SlowMo runtime is on par with that of the base optimizer. We provide theoretical convergence guarantees showing that SlowMo converges to a stationary point of smooth non-convex losses. Since BMUF can be expressed through the SlowMo framework, our results also correspond to the first theoretical convergence guarantees for BMUF.


title: On solvable factors of almost simple groups

abstract: Let $G$ be a finite almost simple group with socle $G_0$. A (nontrivial) factorization of $G$ is an expression of the form $G=HK$, where the factors $H$ and $K$ are core-free subgroups. There is an extensive literature on factorizations of almost simple groups, with important applications in permutation group theory and algebraic graph theory. In a recent paper, Li and Xia describe the factorizations of almost simple groups with a solvable factor $H$. Several infinite families arise in the context of classical groups and in each case a solvable subgroup of $G_0$ containing $H \cap G_0$ is identified. Building on this earlier work, in this paper we compute a sharp lower bound on the order of a solvable factor of every almost simple group and we determine the exact factorizations with a solvable factor. As an application, we describe the finite primitive permutation groups with a nilpotent regular subgroup, extending classical results of Burnside and Schur on cyclic regular subgroups, and more recent work of Li in the abelian case.


title: Constraints on terrestrial planet formation timescales and equilibration processes in the Grand Tack scenario from Hf-W isotopic evolution

abstract: We examine 141 N-body simulations of terrestrial planet late-stage accretion that use the Grand Tack scenario, coupling the collisional results with a hafnium-tungsten (Hf-W) isotopic evolution model. Accretion in the Grand Tack scenario results in faster planet formation than classical accretion models because of higher planetesimal surface density induced by a migrating Jupiter. Planetary embryos that grow rapidly experience radiogenic ingrowth of mantle $^{182}$W that is inconsistent with the measured terrestrial composition, unless much of the tungsten is removed by an impactor core that mixes thoroughly with the target mantle. For physically Earth-like surviving planets, we find that the fraction of equilibrating impactor core $k_\text{core} \geq 0.6$ is required to produce results agreeing with observed terrestrial tungsten anomalies (assuming equilibration with relatively large volumes of target mantle material; smaller equilibrating mantle volumes would require even larger $k_\text{core}$). This requirement of substantial core re-equilibration may be difficult to reconcile with fluid dynamical predictions and hydrocode simulations of mixing during large impacts, and hence this result does not favor the rapid planet building that results from Grand Tack accretion.


title: Assessment of Rydberg Atoms for Wideband Electric Field Sensing

abstract: Rydberg atoms have attracted significant interest recently as electric field sensors. In order to assess potential applications, detailed understanding of relevant figures of merit is necessary, particularly in relation to other, more mature, sensor technologies. Here we present a quantitative analysis of the Rydberg sensor's sensitivity to oscillating electric fields with frequencies between 1 kHz and 1 THz. Sensitivity is calculated using a combination of analytical and semi-classical Floquet models. Using these models, optimal sensitivity at arbitrary field frequency is determined. We validate the numeric Floquet model via experimental Rydberg sensor measurements over a range of 1-20 GHz. Using analytical models, we compare with two prominent electric field sensor technologies: electro-optic crystals and dipole antenna-coupled passive electronics.


title: Spontaneous Wave Function Collapse with Frame Dragging and Induced Gravity

abstract: I impose the Newtonian criteria of inertial frames on the c.o.m. trajectories of massive objects undergoing spontaneous collapse of their wave function. The corresponding modification of the so far used stochastic Schr\"odinger equation eliminates the Brownian motion of the c.o.m., restores the exact inertial motion for free masses. For the collapse of Schr\"odinger cat states the Born rule is satisfied invariably. The proposed machinery comes from the radical assumption that, in the vicinity of the spontaneously localized mass, the stochastic fluctuations of the c.o.m. ---inevitable in the collapse process--- would drag the physical inertial frame with themselves. The perspective of a general theory is presented where the spontaneous-collapse-caused breakdown of local energy-momentum conservation could be remedied by altering the metric, resulting in collapse-induced curvature of the space-time.


title: Unitarity of higher derivative spin-3 models

abstract: The fourth order in derivatives New Massive Gravity model NMG, describes a massive spin-2 particle in $D=2+1$. At the linearized level a proof of unitarity necessarily implies that the generalization to higher dimensions includes non-unitary massless spin-2 modes. The linearized version of NMG is dual to the Fierz-Pauli model FP. Here we examine the unitarity of higher derivative spin-3 models, analogues to the NMG, dual to the Singh-Hagen model SH. We find that the same kind of restriction on the dimension of the space also happens in this case, and the models are physical only in $D=2+1$.


title: Localization-based two-photon wave-function information encoding

abstract: In quantum communications, quantum states are employed for the transmission of information between remote parties. This usually requires sharing knowledge of the measurement bases through a classical public channel in the sifting phase of the protocol. Here, we demonstrate a quantum communication scheme where the information on the bases is shared "non-classically", by encoding this information in the same photons used for carrying the data. This enhanced capability is achieved by exploiting the localization of the photonic wave function, observed when the photons are prepared and measured in the same quantum basis. We experimentally implement our scheme by using a multi-mode optical fiber coupled to an adaptive optics setup. We observe a decrease in the error rate for higher dimensionality, indicating an improved resilience against noise.


title: pISTA-SENSE-ResNet for Parallel MRI Reconstruction

abstract: Magnetic resonance imaging has been widely applied in clinical diagnosis, however, is limited by its long data acquisition time. Although imaging can be accelerated by sparse sampling and parallel imaging, achieving promising reconstruction images with a fast reconstruction speed remains a challenge. Recently, deep learning approaches have attracted a lot of attention for its encouraging reconstruction results but without a proper interpretability. In this letter, to enable high-quality image reconstruction for the parallel magnetic resonance imaging, we design the network structure from the perspective of sparse iterative reconstruction and enhance it with the residual structure. The experimental results of a public knee dataset show that compared with the optimization-based method and the latest deep learning parallel imaging methods, the proposed network has less error in reconstruction and is more stable under different acceleration factors.


title: Memory Centric Characterization and Analysis of SPEC CPU2017 Suite

abstract: In this paper we provide a comprehensive, memory-centric characterization of the SPEC CPU2017 benchmark suite, using a number of mechanisms including dynamic binary instrumentation, measurements on native hardware using hardware performance counters and OS based tools. We present a number of results including working set sizes, memory capacity consumption and, memory bandwidth utilization of various workloads. Our experiments reveal that the SPEC CPU2017 workloads are surprisingly memory intensive, with approximately 50% of all dynamic instructions being memory intensive ones. We also show that there is a large variation in the memory footprint and bandwidth utilization profiles of the entire suite, with some benchmarks using as much as 16 GB of main memory and up to 2.3 GB/s of memory bandwidth. We also perform instruction execution and distribution analysis of the suite and find that the average instruction count for SPEC CPU2017 workloads is an order of magnitude higher than SPEC CPU2006 ones. In addition, we also find that FP benchmarks of the SPEC 2017 suite have higher compute requirements: on average, FP workloads execute three times the number of compute operations as compared to INT workloads.


title: Automated Weed Detection in Aerial Imagery with Context

abstract: In this paper, we demonstrate the ability to discriminate between cultivated maize plant and grass or grass-like weed image segments using the context surrounding the image segments. While convolutional neural networks have brought state of the art accuracies within object detection, errors arise when objects in different classes share similar features. This scenario often occurs when objects in images are viewed at too small of a scale to discern distinct differences in features, causing images to be incorrectly classified or localized. To solve this problem, we will explore using context when classifying image segments. This technique involves feeding a convolutional neural network a central square image along with a border of its direct surroundings at train and test times. This means that although images are labelled at a smaller scale to preserve accurate localization, the network classifies the images and learns features that include the wider context. We demonstrate the benefits of this context technique in the object detection task through a case study of grass (foxtail) and grass-like (yellow nutsedge) weed detection in maize fields. In this standard situation, adding context alone nearly halved the error of the neural network from 7.1% to 4.3%. After only one epoch with context, the network also achieved a higher accuracy than the network without context did after 50 epochs. The benefits of using the context technique are likely to particularly evident in agricultural contexts in which parts (such as leaves) of several plants may appear similar when not taking into account the context in which those parts appear.


title: Smart Palm: An IoT Framework for Red Palm Weevil Early Detection

abstract: Smart agriculture is an evolving trend in agriculture industry, where sensors are embedded into plants to collect vital data and help in decision making to ensure higher quality of crops and prevent pests, disease, and other possible threats. In Saudi Arabia, growing palms is the most important agricultural activity, and there is an increasing need to leverage smart agriculture technology to improve the production of dates and prevent diseases. One of the most critical diseases of palms if the red palm weevil, which is an insect that causes a lot of damage to palm trees and can devast large areas of palm trees. The most challenging problem is that the effect of the weevil is not visible by humans until the palm reaches an advanced infestation state. For this reason, there is a need to use advanced technology for early detection and prevention of infestation propagation. In this project, we have developed am IoT based smart palm monitoring prototype as a proof-of-concept that (1) allows to monitor palms remotely using smart agriculture sensors, (2) contribute to the early detection of red palm weevil. Users can use web/mobile application to interact with their palm farms and help them in getting early detection of possible infestations. We used Elm company IoT platform to interface between the sensor layer and the user layer. In addition, we have collected data using accelerometer sensors and we applied signal processing and statistical techniques to analyze collected data and determine a fingerprint of the infestation.


title: From Heisenberg-Euler Lagrangian to the discovery of chromomagnetic gluon condensation

abstract: I reexamine the phenomena of the chromomagnetic gluon condensation in Yang-Mills theory. The extension of the Heisenberg-Euler Lagrangian to the Yang-Mills theory allows to calculate the effective action, the energy-momentum tensor and demonstrate that the energy density curve crosses the zero energy level of the perturbative vacuum state at nonzero angle and continuously enters to the negative energy density region. At the crossing point and further down the effective coupling constant is small and demonstrate that the true vacuum state of the Yang-Mills theory is below the perturbative vacuum state and is described by the nonzero chromomagnetic gluon condensate. The renormalisation group analysis allows to express the energy momentum tensor, its trace and the first and second order derivatives in terms of Callan-Symanzik beta function and effective coupling constant. The derivatives define the convexity and the extremum of the energy density curve. In the vacuum the energy-momentum tensor is proportional to the space-time metric, and induces a negative contribution to the effective cosmological constant.


title: Descendant distributions for the impact of mutant contagion on networks

abstract: Contagion, broadly construed, refers to anything that can spread infectiously from peer to peer. Examples include communicable diseases, rumors, misinformation, ideas, innovations, bank failures, and electrical blackouts. Sometimes, as in the 1918 Spanish flu epidemic, a contagion mutates at some point as it spreads through a network. Here, using a simple susceptible-infected (SI) model of contagion, we explore the downstream impact of a single mutation event. Assuming that this mutation occurs at a random node in the contact network, we calculate the distribution of the number of "descendants," $d$, downstream from the initial "Patient Zero" mutant. We find that the tail of the distribution decays as $d^{-2}$ for complete graphs, random graphs, small-world networks, networks with block-like structure, and other infinite-dimensional networks. This prediction agrees with the observed statistics of memes propagating and mutating on Facebook, and is expected to hold for other effectively infinite-dimensional networks, such as the global human contact network. In a wider context, our approach suggests a possible starting point for a mesoscopic theory of contagion. Such a theory would focus on the paths traced by a spreading contagion, thereby furnishing an intermediate level of description between that of individual nodes and the total infected population. We anticipate that contagion pathways will hold valuable lessons, given their role as the conduits through which single mutations, innovations, or failures can sweep through a network as a whole.


title: Thresholds in the Lattice of Subspaces of $(\mathbb F_q)^n$

abstract: Let $Q$ be an ideal (downward-closed set) in the lattice of linear subspaces of $(\mathbb F_q)^n$, ordered by inclusion. For $0 \le k \le n$, let $\mu_k(Q)$ denote the fraction of $k$-dimensional subspaces that belong to $Q$. We show that these densities satisfy \[ \mu_k(Q) = \frac{1}{1+z} \quad\Longrightarrow\quad \mu_{k+1}(Q) \le \frac{1}{1+qz}. \] This implies a sharp threshold theorem: if $\mu_k(Q) \le 1-\varepsilon$, then $\mu_\ell(Q) \le \varepsilon$ for $\ell = k + O(\log_q(1/\varepsilon))$.


title: Anion and cation emission from water molecules after collisions with 6.6-keV $^{16}$O$^{+}$ ions

abstract: Anion and cation emission following water dissociation was studied for 6.6-keV $^{16}$O$^{+}$ + H$_{2}$O collisions. Absolute cross sections for the emission of all positively and negatively charged fragments, differential in both energy and observation angle, were measured. The fragments formed in hard, binary collisions appearing in peaks were distinguishable from those created in soft collisions with many-body dynamics that result in a broad energy spectrum. A striking feature is that anions and cations are emitted with similar energy and angular distributions, with a nearly constant ratio of about 1:100 for H$^{-}$ to H$^{+}$. Model calculations were performed at different levels of complexity. Four-body scattering simulations reproduce the measured fragment distributions if adequate kinetic-energy release of the target is taken into account. Providing even further insight into the underlying processes, predictions of a thermodynamic model indicate that transfer ionization at small impact parameters is the dominant mechanism for H$^{+}$ creation. The present findings confirm our earlier observation that in molecular fragmentation induced by slow, singly charged ions, the charge states of the emitted hydrogen fragments follow a simple statistical distribution independent of the way they are formed.


title: Binarized Gerchberg Saxton Algorithm for Hologram Generation Using a Digital Micromirror Device

abstract: In this paper, a modified Gerchberg Saxton algorithm for generating improved robust binary hologram is presented.


title: Forecasting Chaotic Systems with Very Low Connectivity Reservoir Computers

abstract: We explore the hyperparameter space of reservoir computers used for forecasting of the chaotic Lorenz '63 attractor with Bayesian optimization. We use a new measure of reservoir performance, designed to emphasize learning the global climate of the forecasted system rather than short-term prediction. We find that optimizing over this measure more quickly excludes reservoirs that fail to reproduce the climate. The results of optimization are surprising: the optimized parameters often specify a reservoir network with very low connectivity. Inspired by this observation, we explore reservoir designs with even simpler structure, and find well-performing reservoirs that have zero spectral radius and no recurrence. These simple reservoirs provide counterexamples to widely used heuristics in the field, and may be useful for hardware implementations of reservoir computers.


title: On fractional L\'evy processes: tempering, sample path properties and stochastic integration

abstract: We define two new classes of stochastic processes, called tempered fractional L\'{e}vy process of the first and second kinds (TFLP and TFLP $I\!I$, respectively). TFLP and TFLP $I\!I$ make up very broad finite-variance, generally non-Gaussian families of transient anomalous diffusion models that are constructed by exponentially tempering the power law kernel in the moving average representation of a fractional L\'{e}vy process. Accordingly, the increment processes of TFLP and TFLP $I\!I$ display semi-long range dependence. We establish the sample path properties of TFLP and TFLP $I\!I$. We further use a flexible framework of tempered fractional derivatives and integrals to develop the theory of stochastic integration with respect to TFLP and TFLP $I\!I$, which may not be semimartingales depending on the value of the memory parameter and choice of marginal distribution.


title: Beyond effective Hamiltonians: micromotion of Bose Einstein condensates in periodically driven optical lattices

abstract: We investigate a Bose Einstein condensate held in a 1D optical lattice whose phase undergoes a fast oscillation using a statistical analysis. The averaged potential experienced by the atoms boils down to a periodic potential having the same spatial period but with a renormalized depth. However, the atomic dynamics also contains a \emph{micromotion} whose main features are revealed by a Kolmorogov-Smirnov analysis of the experimental momentum distributions. We furthermore discuss the impact of the micromotion on a quench process corresponding to a proper sudden change of the driving amplitude which reverses the curvature of the averaged potential.


title: Enhancing high-content imaging for studying microtubule networks at large-scale

abstract: Given the crucial role of microtubules for cell survival, many researchers have found success using microtubule-targeting agents in the search for effective cancer therapeutics. Understanding microtubule responses to targeted interventions requires that the microtubule network within cells can be consistently observed across a large sample of images. However, fluorescence noise sources captured simultaneously with biological signals while using wide-field microscopes can obfuscate fine microtubule structures. Such requirements are particularly challenging for high-throughput imaging, where researchers must make decisions related to the trade-off between imaging quality and speed. Here, we propose a computational framework to enhance the quality of high-throughput imaging data to achieve fast speed and high quality simultaneously. Using CycleGAN, we learn an image model from low-throughput, high-resolution images to enhance features, such as microtubule networks in high-throughput low-resolution images. We show that CycleGAN is effective in identifying microtubules with 0.93+ AUC-ROC and that these results are robust to different kinds of image noise. We further apply CycleGAN to quantify the changes in microtubule density as a result of the application of drug compounds, and show that the quantified responses correspond well with known drug effects


title: A Computationally Efficient Pipeline Approach to Full Page Offline Handwritten Text Recognition

abstract: Offline handwriting recognition with deep neural networks is usually limited to words or lines due to large computational costs. In this paper, a less computationally expensive full page offline handwritten text recognition framework is introduced. This framework includes a pipeline that locates handwritten text with an object detection neural network and recognises the text within the detected regions using features extracted with a multi-scale convolutional neural network (CNN) fed into a bidirectional long short term memory (LSTM) network. This framework achieves comparable error rates to state of the art frameworks while using less memory and time. The results in this paper demonstrate the potential of this framework and future work can investigate production ready and deployable handwritten text recognisers.


title: Freeness and equivariant stable homotopy

abstract: We introduce a notion of freeness for $RO$-graded equivariant generalized homology theories, considering spaces or spectra $E$ such that the $R$-homology of $E$ splits as a wedge of the $R$-homology of induced virtual representation spheres. The full subcategory of these spectra is closed under all of the basic equivariant operations, and this greatly simplifies computation. Many examples of spectra and homology theories are included along the way. We refine this to a collection of spectra analogous to the pure and isotropic spectra considered by Hill--Hopkins--Ravenel. For these spectra, the $RO$-graded Bredon homology is extremely easy to compute, and if these spaces have additional structure, then this can also be easily determined. In particular, the homology of a space with this property naturally has the structure of a co-Tambara functor (and compatibly with any additional product structure). We work this out in the example of $BU_{\mathbb R}$ and coinduced versions of this. We finish by describing a readily computable bar and twisted bar spectra sequence, giving Bredon homology for various $E_{\infty}$ pushouts, and we apply this to describe the homology of $BBU_{\mathbb R}$.


title: On the Ejection of Dark Matter from Globular Clusters

abstract: We investigate analytically whether in a close encounter with a star, a Dark Matter particle can be accelerated above the escape speed of a Globular Cluster and be ejected. We find that this mechanism is not sufficient to eject a massive, extended Dark Matter halo by the present time. Combined with observations of isolated Globular Clusters that may not have had their halos tidally stripped, these results cast doubt on the scenario in which Globular Clusters formed in Dark Matter halos.


title: A Cepheid-Based Distance to the Seyfert Galaxy NGC 6814

abstract: We present a Cepheid-based distance to the nearby Seyfert galaxy NGC\,6814 from {\it Hubble Space Telescope} observations. We obtained F555W and F814W imaging over the course of 12 visits with logarithmic time spacing in 2013 August$-$October. We detected and made photometric measurements for 16,469 unique sources across all images in both filters, from which we identify 90 excellent Cepheid candidates spanning a range of periods of $13-84$ days. We find evidence for incompleteness in the detection of candidates at periods <21 days. Based on the analysis of Cepheid candidates above the incompleteness limit, we determine a distance modulus for NGC\,6814 relative to the LMC of $\mu_{\rm rel\,LMC}=13.200^{+0.031}_{-0.031}$ mag. Adopting the recent constraint of the distance modulus to the LMC determined by Pietrzynski et al., we find $m-M=31.677^{+0.041}_{-0.041}$ which gives a distance of $21.65 \pm 0.41$ Mpc to NGC 6814.


title: Covariance Matrix Estimation with Non Uniform and Data Dependent Missing Observations

abstract: In this paper we study covariance estimation with missing data. We consider missing data mechanisms that can be independent of the data, or have a time varying dependency. Additionally, observed variables may have arbitrary (non uniform) and dependent observation probabilities. For each mechanism, we construct an unbiased estimator and obtain bounds for the expected value of their estimation error in operator norm. Our bounds are equivalent, up to constant and logarithmic factors, to state of the art bounds for complete and uniform missing observations. Furthermore, for the more general non uniform and dependent cases, the proposed bounds are new or improve upon previous results. Our error estimates depend on quantities we call scaled effective rank, which generalize the effective rank to account for missing observations. All the estimators studied in this work have the same asymptotic convergence rate (up to logarithmic factors).


title: Wasserstein Neural Processes

abstract: Neural Processes (NPs) are a class of models that learn a mapping from a context set of input-output pairs to a distribution over functions. They are traditionally trained using maximum likelihood with a KL divergence regularization term. We show that there are desirable classes of problems where NPs, with this loss, fail to learn any reasonable distribution. We also show that this drawback is solved by using approximations of Wasserstein distance which calculates optimal transport distances even for distributions of disjoint support. We give experimental justification for our method and demonstrate performance. These Wasserstein Neural Processes (WNPs) maintain all of the benefits of traditional NPs while being able to approximate a new class of function mappings.


title: Redshift/blueshift inside the Schwarzschild black hole

abstract: We consider an observer who moves under the horizon of the Schwarzschild black hole and absorbs a photon. There are two different situations when (i) a photon comes from infinity, (ii) it is emitted by another observer under the horizon (say, by a surface of a collapsing star). We analyze the frequency change for absorption near the event horizon, in an intermediate region and near the singularity and compare the results for both scenarios. Near the singularity, in both scenarios infinite redshift (in the pure radial case) or infinite blueshift are possible, depending on angular momenta of a photon and an observer. The main difference between both scenarios manifests itself, if photons from the star surface are emitted near the horizon and received in the intermediate region or near the singularity by the observer with a nonzero angular momentum.


title: Algebraic structures on graph associahedra

abstract: M. Carr and S. Devadoss introduced in [7] the notion of tubing on a finite simple graph $\Gamma$, in the context of configuration spaces on the Hilbert plane. To any finite simple graph $\Gamma$ they associated a finite partially ordered set, whose elements are the tubings of $\Gamma$ and whose geometric realization is a convex polytope ${\mathcal K}\Gamma$, the graph-associahedron. For the complete graphs they recovered permutahedra, for linear graphs they got Stasheff's associahedra, while for simple graph they obtained the standard simplexes. The goal of the present work is to give an \emph{algebraic} description of graph associahedra. We introduce a substitution operation on tubings, which allows us to describe the set of faces of graph-associahedra as a free object, spanned by the set of all connected simple graphs, under operations given via connected subgraphs. The boundary maps of graph-associahedra defines natural derivations in this context. Along the way, we introduce a topological interpretation of the graph tubings and our new operations. In the last section, we show that substitution of tubings may be understood in the context of M. Batanin and M. Markl's operadic categories.


title: Planckian superconductor

abstract: The Planckian relaxation rate $\hbar/t_\mathrm{P} = 2\pi k_\mathrm{B} T$ sets a characteristic time scale for both equilibration of quantum critical systems and maximal quantum chaos. In this note, we show that at the critical coupling between a superconducting dot and the complex Sachdev-Ye-Kitaev model, known to be maximally chaotic, the pairing gap $\Delta$ behaves as $\eta \,\, \hbar/t_\mathrm{P}$ at low temperatures, where $\eta$ is an order one constant. The lower critical temperature emerges with a further increase of the coupling strength so that the finite $\Delta$ domain is settled between the two critical temperatures.


title: Satellite Constellation Pattern Optimization for Complex Regional Coverage

abstract: The use of regional coverage satellite constellations is on the rise, urging the need for an optimal constellation design method for complex regional coverage. Traditional constellations are often designed for continuous global coverage, and the few existing regional constellation design methods lead to suboptimal solutions for periodically time-varying or spatially-varying regional coverage requirements. This paper introduces a new general approach to design an optimal constellation pattern that satisfies such complex regional coverage requirements. To this end, the circular convolution nature of the repeating ground track orbit and common ground track constellation is formalized. This formulation enables a scalable constellation pattern analysis for multiple target areas and with multiple sub-constellations. The formalized circular convolution relationship is first used to derive a baseline constellation pattern design method with the conventional assumption of symmetry. Next, a novel method based on binary integer linear programming is developed, which aims to optimally design a constellation pattern with the minimum number of satellites. This binary integer linear programming method is shown to achieve optimal constellation patterns for general problem settings that the baseline method cannot achieve. Five illustrative examples are analyzed to demonstrate the value of the proposed new approach.


title: Learning to estimate label uncertainty for automatic radiology report parsing

abstract: Bootstrapping labels from radiology reports has become the scalable alternative to provide inexpensive ground truth for medical imaging. Because of the domain specific nature, state-of-the-art report labeling tools are predominantly rule-based. These tools, however, typically yield a binary 0 or 1 prediction that indicates the presence or absence of abnormalities. These hard targets are then used as ground truth to train image models in the downstream, forcing models to express high degree of certainty even on cases where specificity is low. This could negatively impact the statistical efficiency of image models. We address such an issue by training a Bidirectional Long-Short Term Memory Network to augment heuristic-based discrete labels of X-ray reports from all body regions and achieve performance comparable or better than domain-specific NLP, but with additional uncertainty estimates which enable finer downstream image model training.


title: An analytical solution of the Social Force Model for uni-directional flow

abstract: A function for the dependence of flow on pedestrian density is derived analytically from the Social Force Model (SFM) for the case of a homogeneous population walking in the same direction and being in steady state. Assuming that only nearest Voronoi neighbors effectuate forces the resulting function matches a variety of very different fundamental diagrams that were found empirically.


title: Series Representation of Jointly S$\alpha$S Distribution via Symmetric Covariations

abstract: We introduce the notion of symmetric covariation, which is a new measure of dependence between two components of a symmetric $\alpha$-stable random vector, where the stability parameter $\alpha$ measures the heavy-tailedness of its distribution. Unlike covariation that exists only when $\alpha\in(1,2]$, symmetric covariation is well defined for all $\alpha\in(0,2]$. We show that symmetric covariation can be defined using the proposed generalized fractional derivative, which has broader usages than those involved in this work. Several properties of symmetric covariation have been derived. These are either similar to or more general than those of the covariance functions in the Gaussian case. The main contribution of this framework is the representation of the characteristic function of bivariate symmetric $\alpha$-stable distribution via convergent series based on a sequence of symmetric covariations. This series representation extends the one of bivariate Gaussian.


title: Heliosheath Properties Measured from a Voyager 2 to Voyager 1 Transient

abstract: In mid-2012, a GMIR observed by Voyager 2 crossed through the heliosheath and collided with the heliopause, generating a pressure pulse that propagated into the very local interstellar medium. The effects of the transmitted wave were seen by Voyager 1 just 93 days after its own heliopause crossing. The passage of the transient was accompanied by long-lasting decreases in galactic cosmic ray intensities that occurred from ~2012.55 to ~2013.35 and ~2012.91 to ~2013.70 at Voyager 2 and Voyager 1, respectively. Omnidirectional (>20 MeV) proton-dominated measurements from each spacecraft's Cosmic Ray Subsystem reveal a remarkable similarity between these causally-related events, with a correlation coefficient of 91.2% and a time-lag of 130 days. Knowing the locations of the two spacecraft, we use the observed time-delay to calculate the GMIR's average speed through the heliosheath (inside the heliopause) as a function of temperature in the very local interstellar medium. This, combined with particle, field, and plasma observations enables us to infer previously unmeasured properties of the heliosheath, including a range of sound speeds and total effective pressures. For a nominal temperature of ~20,000 K just outside the heliopause, we find a sound speed of 314 (+/-) 32 km/s and total effective pressure of 267 (+/-) 55 fPa inside the heliopause. We compare these results with the Interstellar Boundary Explorer's data-driven models of heliosheath pressures derived from energetic neutral atom fluxes (the globally distributed flux) and present them as additional evidence that the heliosheath's dynamics are driven by suprathermal energetic processes.


title: NB-IoT Small Cell: A 3GPP Perspective

abstract: Narrowband Internet of Things (NB-IoT) technology was introduced in 3GPP Release 13 to accommodate device-generated traffic over cellular networks. There have been efforts in Release 14 to deliver further improvements to facilitate the deployment of Internet of Things (IoT) use cases, such as coverage enhancement, and support for low cost, low complexity and low power consumption device. Small cell NB-IoT is added as an attractive feature in 3GPP Release 15. In this article, we provide an overview of this feature and try to shed light on major aspects of small cell deployment in NB-IoT systems.


title: Implicit Trajectory Planning for Feedback Linearizable Systems: A Time-varying Optimization Approach

abstract: We develop an optimization-based framework for joint real-time trajectory planning and feedback control of feedback-linearizable systems. To achieve this goal, we define a target trajectory as the optimal solution of a time-varying optimization problem. In general, however, such trajectory may not be feasible due to , e.g., nonholonomic constraints. To solve this problem, we design a control law that generates feasible trajectories that asymptotically converge to the target trajectory. More precisely, for systems that are (dynamic) full-state linearizable, the proposed control law implicitly transforms the nonlinear system into an optimization algorithm of sufficiently high order. We prove global exponential convergence to the target trajectory for both the optimization algorithm and the original system. We illustrate the effectiveness of our proposed method on multi-target or multi-agent tracking problems with constraints.


title: TagSLAM: Robust SLAM with Fiducial Markers

abstract: TagSLAM provides a convenient, flexible, and robust way of performing Simultaneous Localization and Mapping (SLAM) with AprilTag fiducial markers. By leveraging a few simple abstractions (bodies, tags, cameras), TagSLAM provides a front end to the GTSAM factor graph optimizer that makes it possible to rapidly design a range of experiments that are based on tags: full SLAM, extrinsic camera calibration with non-overlapping views, visual localization for ground truth, loop closure for odometry, pose estimation etc. We discuss in detail how TagSLAM initializes the factor graph in a robust way, and present loop closure as an application example. TagSLAM is a ROS based open source package and can be found at https://berndpfrommer.github.io/tagslam_web.


title: Compactness by coarse-graining in long-range lattice systems

abstract: We consider energies on a periodic set ${\mathcal L}$ of ${\mathbb R}^d$ of the form $\sum_{i,j\in{\mathcal L}} a^\varepsilon_{ij}|u_i-u_j|$, defined on spin functions $u_i\in\{0,1\}$, and we suppose that the typical range of the interactions is $R_\varepsilon$ with $R_\varepsilon\to +\infty$, i.e., if $\|i-j\|\le R_\varepsilon$ then $a^\varepsilon_{ij}\ge c>0$. In a discrete-to-continuum analysis, we prove that the overall behaviour as $\varepsilon\to 0$ of such functionals is that of an interfacial energy. The proof is performed using a coarse-graining procedure which associates to scaled functions defined on $\varepsilon{\mathcal L}$ with equibounded energy a family of sets with equibounded perimeter. This agrees with the case of equibounded $R_\varepsilon$ and can be seen as an extension of coerciveness result for short-range interactions, but is different from that of other long-range interaction energies, whose limit exits the class of surface energies. A computation of the limit energy is performed in the case ${\mathcal L}={\mathbb Z}^d$.


title: An Iterative Quadratic Method for General-Sum Differential Games with Feedback Linearizable Dynamics

abstract: Iterative linear-quadratic (ILQ) methods are widely used in the nonlinear optimal control community. Recent work has applied similar methodology in the setting of multiplayer general-sum differential games. Here, ILQ methods are capable of finding local equilibria in interactive motion planning problems in real-time. As in most iterative procedures, however, this approach can be sensitive to initial conditions and hyperparameter choices, which can result in poor computational performance or even unsafe trajectories. In this paper, we focus our attention on a broad class of dynamical systems which are feedback linearizable, and exploit this structure to improve both algorithmic reliability and runtime. We showcase our new algorithm in three distinct traffic scenarios, and observe that in practice our method converges significantly more often and more quickly than was possible without exploiting the feedback linearizable structure.


title: Accelerated Robot Learning via Human Brain Signals

abstract: In reinforcement learning (RL), sparse rewards are a natural way to specify the task to be learned. However, most RL algorithms struggle to learn in this setting since the learning signal is mostly zeros. In contrast, humans are good at assessing and predicting the future consequences of actions and can serve as good reward/policy shapers to accelerate the robot learning process. Previous works have shown that the human brain generates an error-related signal, measurable using electroencephelography (EEG), when the human perceives the task being done erroneously. In this work, we propose a method that uses evaluative feedback obtained from human brain signals measured via scalp EEG to accelerate RL for robotic agents in sparse reward settings. As the robot learns the task, the EEG of a human observer watching the robot attempts is recorded and decoded into noisy error feedback signal. From this feedback, we use supervised learning to obtain a policy that subsequently augments the behavior policy and guides exploration in the early stages of RL. This bootstraps the RL learning process to enable learning from sparse reward. Using a robotic navigation task as a test bed, we show that our method achieves a stable obstacle-avoidance policy with high success rate, outperforming learning from sparse rewards only that struggles to achieve obstacle avoidance behavior or fails to advance to the goal.


title: Spin dependent ab initio nonlocal No-Core Shell-Model One-Body Densities

abstract: Constructing microscopic effective interactions (`optical potentials') for nucleon-nucleus (NA) elastic scattering requires in first order off-shell nucleon-nucleon (NN) scattering amplitudes between the projectile and the struck target nucleon and nonlocal one-body density matrices. While the NN amplitudes and the {\it ab intio} no-core shell-model (NCSM) calculations always contain the full spin structure of the NN problem, one-body density matrices used in traditional microscopic folding potential neglect spin contributions inherent in the one-body density matrix. Here we derive and show the expectation values of the spin-orbit contribution of the struck nucleon with respect to the rest of the nucleus for $^{4}$He, $^{6}$He, $^{12}$C, and $^{16}$O and compare them with the scalar one-body density matrix.


title: Orbit Characterization, Stabilization and Composition on 3D Underactuated Bipedal Walking via Hybrid Passive Linear Inverted Pendulum Model

abstract: A Hybrid passive Linear Inverted Pendulum (HLIP) model is proposed for characterizing, stabilizing and composing periodic orbits for 3D underactuated bipedal walking. Specifically, Period-1 (P1) and Period-2 (P2) orbits are geometrically characterized in the state space of the H-LIP. Stepping controllers are designed for global stabilization of the orbits. Valid ranges of the gains and their optimality are derived. The optimal stepping controller is used to create and stabilize the walking of bipedal robots. An actuated Spring-loaded Inverted Pendulum (aSLIP) model and the underactuated robot Cassie are used for illustration. Both the aSLIP walking with P1 or P2 orbits and the Cassie walking with all 3D compositions of the P1 and P2 orbits can be smoothly generated and stabilized from a stepping-in-place motion. This approach provides a perspective and a methodology towards continuous gait generation and stabilization for 3D underactuated walking robots.


title: Detection of new methanol maser transitions associated with G358.93-0.03

abstract: We report the detection of new 12.178, 12.229, 20.347, and 23.121 GHz methanol masers in the massive star-forming region G358.93-0.03, which are flaring on similarly short timescales (days) as the 6.668 GHz methanol masers also associated with this source. The brightest 12.178 GHz channel increased by a factor of over 700 in just 50 d. The masers found in the 12.229 and 20.347 GHz methanol transitions are the first ever reported and this is only the fourth object to exhibit associated 23.121 GHz methanol masers. The 12.178 GHz methanol maser emission appears to have a higher flux density than that of the 6.668 GHz emission, which is unusual. No associated near-infrared flare counterpart was found, suggesting that the energy source of the flare is deeply embedded.


title: Superconductivity in Nb5Ir3-xPtxO

abstract: We have investigated the superconducting critical temperature Tc of solid solution Nb5Ir3-xPtxO with the Ti5Ga4-type structure. The both end-members of Nb5Ir3O and Nb5Pt3O are reported to be superconductors with Tc of 10.5 K and 3.8 K, respectively. Particularly Nb5Ir3O is considered as a two-gap superconductor. The entire series of alloy hold the Ti5Ga4-type structure and show the linear x dependence of lattice parameters. On the other hand, a nonlinear Tc vs x plot was obtained, suggesting that Tc is not determined by only lattice parameters. The experimental result is discussed based on the Matthias rule and compared with those of the other two-gap superconductors such as MgB2, Mo8Ga41 and Nb3Sn.


title: Motion Decoupling and Composition via Reduced Order Model Optimization for Dynamic Humanoid Walking with CLF-QP based Active Force Control

abstract: In this paper, 3D humanoid walking is decoupled into periodic and transitional motion, each of which is decoupled into planar walking in the sagittal and lateral plane. Reduced order models (ROMs), i.e. actuated Spring-loaded Inverted Pendulum (aSLIP) models and Hybrid-Linear Inverted Pendulum (H-LIP) models, are utilized for motion generation on the desired center of mass (COM) dynamics for each type of planar motion. The periodic motion is planned via point foot (underactuated) ROMs for dynamic motion with minimum ankle actuation, while the transitional motion is planned via foot-actuated ROMs for fast and smooth transition. Composition of the planar COM dynamics yields the desired COM dynamics in 3D, which is embedded on the humanoid via control Lyapunov function based Quadratic programs (CLF-QPs). Additionally, the ground reaction force profiles of the aSLIP walking are used as desired references for ground contact forces in the CLF-QPs for smooth domain transitions. The proposed framework is realized on a lower-limb exoskeleton in simulation wherein different walking motions are achieved.


title: The Effects of Telluric Contamination in Iodine Calibrated Precise Radial Velocities

abstract: We characterized the effects of telluric absorption lines on the radial velocity (RV) precision of stellar spectra taken through an iodine cell. To isolate the effects induced by telluric contamination from other stellar, instrumental, or numerical systematic RV noise, we extracted RVs from simulated iodine calibrated spectra of three RV standard stars regularly observed by Keck/HIRES. We add in water absorption lines according to measured precipitable water vapor (PWV) contents over a one-year period. We concluded that telluric contamination introduces additional RV noise and spurious periodic signals on the level of 10-20 cm/s, consistent with similar previous studies. Our findings show that forward modeling the telluric lines effectively recovers the RV precision and accuracy, with no prior knowledge of the PWV needed. Such a recovery is less effective when the water absorption lines are relatively deep in the stellar template used in the forward modeling. Overall, telluric contamination plays an insignificant role for typical iodine-calibrated RV programs aiming at ~1-2 m/s, but we recommend adding modeling of telluric lines and taking stellar template observations on nights with low humidity for programs aiming to achieve sub-m/s precision.


title: Algebras from Congruences

abstract: We present a functorial construction which, starting from a congruence $\alpha$ of finite index in an algebra A, yields a new algebra C with the following properties: the congruence lattice of C is isomorphic to the interval of congruences between 0 and $\alpha$ on A, this isomorphism preserves higher commutators and TCT types, and C inherits all idempotent Maltsev conditions from A. As applications of this construction, we first show that supernilpotence is decidable for congruences of finite algebras in varieties that omit type 1. Secondly, we prove that the subpower membership problem for finite algebras with a cube term can be effectively reduced to membership questions in subdirect products of subdirectly irreducible algebras with central monoliths. As a consequence, we obtain a polynomial time algorithm for the subpower membership problem for finite algebras with a cube term in which the monolith of every subdirectly irreducible section has a supernilpotent centralizer.


title: Fabrication of a Cosmic Ray Veto System for the Mu2e Experiment

abstract: The Mu2e experiment at Fermilab will search for the charged-lepton flavor-violating process of a neutrino less muon-to-electron decay in the presence of a nucleus. The experiment expects a single-event sensitivity of $2.9\times10^{-17}$, which is four orders of magnitude below the current strongest limits on this process. This requires all backgrounds to sum to fewer than one event over the lifetime of the experiment. One major background is due to cosmic-ray muons producing electrons that fake a signal inside of the Mu2e apparatus. The Mu2e Cosmic Ray Veto (CRV) has been designed to veto these cosmic-ray backgrounds with an efficiency of 99.99%, while causing a low dead time and while operating in a high-intensity environment. The design and fabrication of the CRV is discussed.


title: Average number of solutions

abstract: Let $X$ be an $n$-dimensional manifold and $V_1,\ldots,V_n\subset C^\infty(X,\mathbb R)$ finite-dimensional vector spaces. For systems of equations $\{f_i = a_i\colon\: f_i\in V_i,\:a_i \in\mathbb R,\:i=1,\ldots,n\}$ we discover a relationship between the average number of their solutions and mixed volumes of convex bodies. To do this, we choose Banach metrics in the spaces $V_i$. Using these metrics, we construct 1) the measure in the space of systems, and 2) Banach convex bodies in $X$, i.e., collections of centrally symmetric convex bodies in the fibers of the cotangent bundle of $X$. It turns out that the average number of solutions is equal to the mixed symplectic volume of Banach convex bodies. Earlier this result was obtained for Euclidean metrics in spaces $V_i$. In Euclidean case, the Banach convex bodies are the collections of ellipsoids.


title: Retrieving Top Weighted Triangles in Graphs

abstract: Pattern counting in graphs is a fundamental primitive for many network analysis tasks, and a number of methods have been developed for scaling subgraph counting to large graphs. Many real-world networks carry a natural notion of strength of connection between nodes, which are often modeled by a weighted graph, but existing scalable graph algorithms for pattern mining are designed for unweighted graphs. Here, we develop a suite of deterministic and random sampling algorithms that enable the fast discovery of the 3-cliques (triangles) with the largest weight in a graph, where weight is measured by a generalized mean of a triangle's edges. For example, one of our proposed algorithms can find the top-1000 weighted triangles of a weighted graph with billions of edges in thirty seconds on a commodity server, which is orders of magnitude faster than existing "fast" enumeration schemes. Our methods thus open the door towards scalable pattern mining in weighted graphs.


title: Tracking Control by the Newton-Raphson Method with Output Prediction and Controller Speedup

abstract: This paper presents a control technique for output tracking of reference signals in continuous-time dynamical systems. The technique is comprised of the following three elements: (i) output prediction which has to track the reference signal, (ii) a controller based on an integrator with variable gain, and (iii) a speedup of the control action for enhancing the tracker's accuracy and, in some cases, guaranteeing stability of the closed-loop system. The technique is suitable for linear and nonlinear systems, implementable by simple algorithms, can track reference points as well as time-dependent reference signals, and may have large, even global domains of attraction. The derived theoretical results include convergence of the tracking controller and error analysis, and are supported by illustrative simulation and laboratory experiments.


title: RITnet: Real-time Semantic Segmentation of the Eye for Gaze Tracking

abstract: Accurate eye segmentation can improve eye-gaze estimation and support interactive computing based on visual attention; however, existing eye segmentation methods suffer from issues such as person-dependent accuracy, lack of robustness, and an inability to be run in real-time. Here, we present the RITnet model, which is a deep neural network that combines U-Net and DenseNet. RITnet is under 1 MB and achieves 95.3\% accuracy on the 2019 OpenEDS Semantic Segmentation challenge. Using a GeForce GTX 1080 Ti, RITnet tracks at $>$ 300Hz, enabling real-time gaze tracking applications. Pre-trained models and source code are available https://bitbucket.org/eye-ush/ritnet/.


title: Real-Time Thermospheric Density Estimation Via Two-Line-Element Data Assimilation

abstract: Inaccurate estimates of the thermospheric density are a major source of error in low Earth orbit prediction. To improve orbit prediction, real-time density estimation is required. In this work, we develop a reduced-order dynamic model for the thermospheric density by computing the main spatial modes of the atmosphere and deriving a linear model for the dynamics. The model is then used to estimate the density using two-line element (TLE) data by simultaneously estimating the reduced-order modes and the orbits and ballistic coefficients of several objects using an unscented Kalman filter. Accurate density estimation using the TLEs of 17 objects is demonstrated and validated against CHAMP and GRACE accelerometer-derived densities. Finally, the use of the model for density forecasting is shown.

categories: physics.space-ph astro-ph.EP math.DS physics.ao-ph


title: Improvement of Multiparametric MR Image Segmentation by Augmenting the Data with Generative Adversarial Networks for Glioma Patients

abstract: Every year thousands of patients are diagnosed with a glioma, a type of malignant brain tumor. Physicians use MR images as a key tool in the diagnosis and treatment of these patients. Neural networks show great potential to aid physicians in the medical image analysis. This study investigates the use of varying amounts of synthetic brain T1-weighted (T1), post-contrast T1-weighted (T1Gd), T2-weighted (T2), and T2 Fluid Attenuated Inversion Recovery (FLAIR) MR images created by a generative adversarial network to overcome the lack of annotated medical image data in training separate 2D U-Nets to segment enhancing tumor, peritumoral edema, and necrosis (non-enhancing tumor core) regions on gliomas. These synthetic MR images were assessed quantitively (SSIM=0.79) and qualitatively by a physician who found that the synthetic images seem stronger for delineation of structural boundaries but struggle more when gradient is significant, (e.g. edema signal in T2 modalities). Multiple 2D U-Nets were trained with original BraTS data and differing subsets of a quarter, half, three-quarters, and all synthetic MR images. There was not an obvious correlation between the improvement of values of the metrics in separate validation dataset for each structure and amount of synthetic data added, there is a strong correlation between the amount of synthetic data added and the number of best overall validation metrics. In summary, this study showed ability to generate high quality synthetic Flair, T2, T1, and T1CE MR images using the GAN. Using the synthetic MR images showed encouraging results to improve the U-Net segmentation performance which has the potential to address the scarcity of readily available medical images.


title: Graphs of bounded cliquewidth are polynomially $\chi$-bounded

abstract: We prove that if $\mathcal{C}$ is a hereditary class of graphs that is polynomially $\chi$-bounded, then the class of graphs that admit decompositions into pieces belonging to $\mathcal{C}$ along cuts of bounded rank is also polynomially $\chi$-bounded. In particular, this implies that for every positive integer $k$, the class of graphs of cliquewidth at most $k$ is polynomially $\chi$-bounded.


title: Re-balancing Variational Autoencoder Loss for Molecule Sequence Generation

abstract: Molecule generation is to design new molecules with specific chemical properties and further to optimize the desired chemical properties. Following previous work, we encode molecules into continuous vectors in the latent space and then decode the vectors into molecules under the variational autoencoder (VAE) framework. We investigate the posterior collapse problem of current RNN-based VAEs for molecule sequence generation. For the first time, we find that underestimated reconstruction loss leads to posterior collapse, and provide both theoretical and experimental evidence. We propose an effective and efficient solution to fix the problem and avoid posterior collapse. Without bells and whistles, our method achieves SOTA reconstruction accuracy and competitive validity on the ZINC 250K dataset. When generating 10,000 unique valid SMILES from random prior sampling, it costs JT-VAE1450s while our method only needs 9s. Our implementation is at https://github.com/chaoyan1037/Re-balanced-VAE.