title: The special linear group for nonassociative rings

abstract: We extend to arbitrary rings a definition of the octonion special linear group due to Baez. At the infinitesimal level we get a Lie ring, which we describe over some large classes of rings, including all associative rings and all algebras over a field. As a corollary we compute all the groups Baez defined.


title: Ab initio studies of the ground and first excited states of the Sr-H$_2$ and Yb-H$_2$ complexes

abstract: Accurate intermolecular potential-energy surfaces (IPESs) for the ground and first excited states of the Sr-H$_2$ and Yb-H$_2$ complexes were calculated. After an extensive methodological study, the CCSD(T) method with the Douglas-Kroll-Hess Hamiltonian and correlation-consistent basis sets of triple-$\zeta$ quality extended with 2 sets of diffuse functions and a set of midbond functions were chosen. The obtained ground-state IPESs are similar in both complexes, being relatively isotropic with two minima and two transition states (equivalent by symmetry). The global minima correspond to the collinear geometries with $R=$ 5.45 and 5.10~{\AA} and energies of $-$27.7 and $-$31.7~cm$^{-1}$ for the Sr-H$_2$ and Yb-H$_2$ systems, respectively. The calculated surfaces for the Sr($^3P$)-H$_2$ and Yb($^3P$)-H$_2$ states are deeper and more anisotropic and they exhibit similar patterns within both complexes. The deepest surfaces, where the singly occupied \textit{p}-orbital of the metal atom is perpendicular to the intermolecular axis, are characterised by the global minima of ca. $-$2053 and $-$2260~cm$^{-1}$ in the T-shape geometries at $R=$ 2.41 and 2.29~{\AA} for Sr-H$_2$ and Yb-H$_2$, respectively. Additional calculations for the complexes of Sr and Yb with the He atom revealed a similar, strong dependence of the interaction energy on the orientation of the \textit{p}-orbital in the the Sr($^3P$)-He and Yb($^3P$)-He states.


title: A Robust Tip of the Red Giant Branch Distance to the Fireworks Galaxy (NGC 6946)

abstract: Archival HST data taken in F606W+F814W of two different fields in the outer regions of NGC 6946 is used to measure a tip of the red giant branch (TRGB) distance to the galaxy. We employ a Bayesian maximum-likelihood modeling method that incorporates the completeness of the photometry, and allows us to model the luminosity function of the RGB population. Our two fields provide us with distances of 7.74 $\pm$ 0.42 Mpc and 7.69 $\pm$ 0.50 Mpc, respectively. Our final distance of 7.72 $\pm$ 0.32 Mpc is higher than most values published previously in the literature. This has important implications for supernova measurements, as NGC 6946 is host to the most observed supernovae (10) of any galaxy to date. We find that the supernovae in NGC 6946 are on average $\sim$ 2.3 times more luminous than previous estimates. Our distance gives NGC 6946 a peculiar velocity of $v_{pec}$ = $-229$ $\pm$ $29$ km/s in the Local Sheet frame. This velocity is the projected component of a substantial negative SGZ motion, indicating a coherent motion along with the other members of the M101 wall toward the plane of the Local Sheet. The M101 wall, a spur off the Local Sheet into the Local Void, is experiencing the expansion of the Local Void.


title: Studying the Milky Way Pulsar Population with Cosmic-Ray Leptons

abstract: Recent measurements of cosmic-ray electron and positron spectra at energies from a GeV to 5 TeV, as well as radio, X-ray and a wide range of gamma-ray observations of pulsar-wind nebulae, indicate that pulsars are significant sources of high-energy cosmic-ray electrons and positrons. Here we calculate the local cosmic-ray $e^\pm$ energy spectra from pulsars taking into account models for (a) the distribution of the pulsars spin-down properties; (b) the cosmic-ray source spectra; and (c) the physics of cosmic-ray propagation. We then use the measured cosmic-ray fluxes from AMS-02, CALET and DAMPE to constrain the space of pulsar and cosmic-ray-propagation models and in particular, local cosmic-ray diffusion and energy losses, the pulsars' energy-loss time-dependence, and the injected $e^{\pm}$ spectra. We find that the lower estimates for the local $e^{\pm}$ energy losses are inconsistent with the data. We also find that pulsar braking indexes of 2.5 or less for sources with ages more than 10 kyr are strongly disfavored. Moreover, the cosmic-ray data are consistent with a wide range of assumptions on the $e^{\pm}$ injection spectral properties and on the distribution of initial spin-down powers. Above a TeV in energy, we find that pulsars can easily explain the observed change in the $e^{+} + e^{-}$ spectral slope. These conclusions are valid as long as pulsars contribute $\gtrsim10\%$ of the observed cosmic-ray $e^\pm$ at energies $\gtrsim100$ GeV.


title: Carbon Chain Molecules Toward Embedded Low-Mass Protostars

abstract: Carbon chain molecules may be an important reservoir of reactive organics during star and planet formation. Carbon chains have been observed toward several low-mass young stellar objects (YSOs), but their typical abundances and chemical relationships in such sources are largely unconstrained. We present a carbon chain survey toward 16 deeply embedded (Class 0/I) low-mass protostars made with the IRAM 30 m telescope. Carbon chains are found to be common at this stage of protostellar evolution. We detect CCS, CCCS, HC$_3$N, HC$_5$N, l-C$_3$H, and C$_4$H toward 88%, 38%, 75%, 31%, 81%, and 88% of sources, respectively. Derived column densities for each molecule vary by one to two orders of magnitude across the sample. As derived from survival analysis, median column densities range between 1.2$\times 10^{11}$ cm$^{-2}$ (CCCS) and 1.5$\times 10^{13}$ cm$^{-2}$ (C$_4$H) and estimated fractional abundances with respect to hydrogen range between 2$\times 10^{-13}$ (CCCS) and 5$\times 10^{-11}$ (C$_4$H), which are low compared to cold cloud cores, warm carbon chain chemistry (WCCC) sources, and protostellar model predictions. We find significant correlations between molecules of the same carbon chain families, as well as between the cyanpolyynes (HC$_{\rm n}$N) and the pure hydrocarbon chains (C$_{\rm n}$H). This latter correlation is explained by a closely-related production chemistry of C$_{\rm{n}}$H and cyanpolyynes during low-mass star formation.


title: Magnon band topology in spin-orbital coupled magnets: classification and application to $\alpha$-RuCl$_3$

abstract: In spite of flourishing studies on the topology of spin waves, a generic framework to classify and compute magnon band topology in non-collinear magnets is still missing. In this work we provide such a theory framework, by mapping an arbitrary linear spin wave into a local free-fermion Hamiltonian with exactly the same spectrum, symmetry implementation and band topology, which allows for a full classification and calculation on any topological properties of magnon bands. We apply this fermionization approach to honeycomb Kitaev magnet $\alpha$-RuCl$_3$, and show the existence of topologically protected magnon band crossings, and field-induced magnon Chern bands under small magnetic fields.


title: Counting of States in Higgs Theories

abstract: We enumerate the micro-states in Higgs theories, addressing (i) the number of vacuum states and (ii) the appropriate measure in the quantum path integral. To address (i) we explicitly construct the set of ground state wave-functionals in the field basis focussing on scalar modes $\theta(x)$. Firstly, we show that in the limit in which the gauge coupling is zero, we obtain an infinite set of degenerate ground states at large volume distinguished by $\theta(x)\to\theta(x)+\theta_0$, spontaneously breaking the global symmetry, as is well known. We then show that at finite gauge coupling there is a unique ground state at large volume since the wave-functional only depends on $\nabla\theta$ in the IR, and we explain this at the level of the Lagrangian. Since gauge fields fall off exponentially from sources there are no conserved charges or symmetries in the Higgs phase; so the Higgs mechanism is the removal of symmetry from the theory. We show how physical features of defects, such as cosmic strings in the abelian Higgs model, are best understood in this context. Since there is a unique ground state, we address (ii) whether the path integral is a volume measure for the radial Higgs field $\mathcal{D}\rho\,\rho^{N-1}$ from the $N$ components of the Higgs multiplet, or a line measure $\mathcal{D}\rho$ as the $N-1$ would-be Goldstones can be removed in unitary gauge. We prove that the need to avoid quartic divergences demands a tower of counter terms that resum to exactly give the volume measure. So the size of the Hilbert space in the zero gauge coupling case and finite gauge coupling case are in one-to-one correspondence, despite the degeneracy of the ground state being lifted in the latter. As a cosmological application, we point out that the volume measure can make it exponentially more unlikely in $N(=4)$ for the Standard Model Higgs to relax to the electroweak vacuum in the early universe.


title: Optimal designs for frequentist model averaging

abstract: We consider the problem of designing experiments for the estimation of a target in regression analysis if there is uncertainty about the parametric form of the regression function. A new optimality criterion is proposed, which minimizes the asymptotic mean squared error of the frequentist model averaging estimate by the choice of an experimental design. Necessary conditions for the optimal solution of a locally and Bayesian optimal design problem are established. The results are illustrated in several examples and it is demonstrated that Bayesian optimal designs can yield a reduction of the mean squared error of the model averaging estimator up to $45\%$.


title: The Hawaii Infrared Parallax Program. III. 2MASS J0249-0557 c: A Wide Planetary-mass Companion to a Low-mass Binary in the beta Pic Moving Group

abstract: We have discovered a wide planetary-mass companion to the $\beta$ Pic moving group member 2MASSJ02495639-0557352 (M6 VL-G) using CFHT/WIRCam astrometry from the Hawaii Infrared Parallax Program. In addition, Keck laser guide star adaptive optics aperture-masking interferometry shows that the host is itself a tight binary. Altogether, 2MASSJ0249-0557ABc is a bound triple system with an $11.6^{+1.0}_{-1.3}$ $M_{\rm Jup}$ object separated by $1950\pm200$ AU (40") from a relatively close ($2.17\pm0.22$ AU, 0.04") pair of $48^{+12}_{-13}$ $M_{\rm Jup}$ and $44^{+11}_{-14}$ $M_{\rm Jup}$ objects. 2MASSJ0249-0557AB is one of the few ultracool binaries to be discovered in a young moving group and the first confirmed in the $\beta$ Pic moving group ($22\pm6$ Myr). The mass, absolute magnitudes, and spectral type of 2MASSJ0249-0557 c (L2 VL-G) are remarkably similar to those of the planet $\beta$ Pic b (L2, $13.0^{+0.4}_{-0.3}$ $M_{\rm Jup}$). We also find that the free-floating object 2MASSJ2208+2921 (L3 VL-G) is another possible $\beta$ Pic moving group member with colors and absolute magnitudes similar to $\beta$ Pic b and 2MASSJ0249-0557 c. $\beta$ Pic b is the first directly imaged planet to have a "twin," namely an object of comparable properties in the same stellar association. Such directly imaged objects provide a unique opportunity to measure atmospheric composition, variability, and rotation across different pathways of assembling planetary-mass objects from the same natal material.


title: Qualitative dynamics and inflationary attractors in loop cosmology

abstract: Qualitative dynamics of three different loop quantizations of spatially flat isotropic and homogeneous models is studied using effective spacetime description of the underlying quantum geometry. These include the standard loop quantum cosmology (LQC), its recently revived modification (referred to as mLQC-I), and another related modification of LQC (mLQC-II) whose dynamics is studied in detail for the first time. Various features of LQC, including quantum bounce and pre-inflationary dynamics, are found to be shared with the mLQC-I and mLQC-II models. We study universal properties of dynamics for chaotic inflation, fractional monodromy inflation, Starobinsky potential, non-minimal Higgs inflation, and an exponential potential. We find various critical points and study their stability, which reveal various qualitative similarities in the post-bounce phase for all these models. The pre-bounce qualitative dynamics of LQC and mLQC-II turns out to be very similar, but is strikingly different from that of mLQC-I. In the dynamical analysis, some of the fixed points turn out to be degenerate for which center manifold theory is used. For all these potentials, non-perturbative quantum gravitational effects always result in a non-singular inflationary scenario with a phase of super-inflation succeeded by the conventional inflation. We show the existence of inflationary attractors, and obtain scaling solutions in the case of the exponential potential. Since all of the models agree with general relativity at late times, our results are also of use in classical theory where qualitative dynamics of some of the potentials has not been studied earlier.


title: irbasis: Open-source database and software for intermediate-representation basis functions of imaginary-time Green's function

abstract: The open-source library, irbasis, provides easy-to-use tools for two sets of orthogonal functions named intermediate representation (IR). The IR basis enables a compact representation of the Matsubara Green's function and efficient calculations of quantum models. The IR basis functions are defined as the solution of an integral equation whose analytical solution is not available for this moment. The library consists of a database of pre-computed high-precision numerical solutions and computational code for evaluating the functions from the database. This paper describes technical details and demonstrates how to use the library.

categories: physics.comp-ph cond-mat.stat-mech cond-mat.str-el cond-mat.supr-con cs.LG


title: Efficient implementation of the continuous-time interaction-expansion quantum Monte Carlo method

abstract: We describe an open-source implementation of the continuous-time interaction-expansion quantum Monte Carlo method for cluster-type impurity models with onsite Coulomb interactions and complex Weiss functions. The code is based on the ALPS libraries.

categories: cond-mat.str-el cond-mat.stat-mech cond-mat.supr-con physics.comp-ph


title: Embedding Banach spaces into the space of bounded functions with countable support

abstract: We prove that a WLD subspace of the space $\ell_\infty^c(\Gamma)$ consisting of all bounded, countably supported functions on a set $\Gamma$ embeds isomorphically into $\ell_\infty$ if and only if it does not contain isometric copies of $c_0(\omega_1)$. Moreover, a subspace of $\ell_\infty^c(\omega_1)$ is constructed that has an unconditional basis, does not embed into $\ell_\infty$, and whose every weakly compact subset is separable (in particular, it cannot contain any isomorphic copies of $c_0(\omega_1)$).


title: Small-Field and Scale-Free: Inflation and Ekpyrosis at their Extremes

abstract: There is increasing evidence from string theory that effective field theories are only reliable over approximately sub-Planckian field excursions. The two most promising effective models for early universe cosmology, inflation and ekpyrosis, are mechanisms that, in order to address cosmological puzzles, must operate over vast expansion/energy ranges. This suggests that it might be appropriate to describe them using scaling laws. Here we combine these two ideas and demonstrate that they drive inflation and ekpyrosis to their extremes: inflation must start at ultra-slow-roll, and ekpyrosis at ultra-fast-roll. At face value, the implied spectra are overly tilted to the red, although in both cases minor departures from pure scale freedom bring the spectral indices within current observational bounds. These models predict a significant spectral running at a level detectable in the near future ($\alpha_s \approx -10^{-3}$). Ekpyrotic models with minimal coupling are nevertheless ruled out, as they lead to levels of non-Gaussianity that are at least an order of magnitude too large. No such restrictions apply to models with a kinetic coupling between the two ekpyrotic scalar fields, and these remain the most promising ekpyrotic models. An additional swampland criterion that was recently proposed for the slope of the scalar field potential would however rule out all ultra-slow-roll models of inflation. Finally, we speculate on the existence of corresponding restrictions on the slope at negative potential values, which might lead to similarly severe constraints on ekpyrotic models.


title: Approximation Algorithms for Clustering via Weighted Impurity Measures

abstract: An impurity measures $I:{R}^k \to {R}^+$ maps a $k$-dimensional vector ${\bf v}$ to a non-negative value $I({\bf v})$ so that the more homogeneous ${\bf v}$, the larger its impurity. We study clustering based on impurity measures: given a collection $V$ of $n$ many $k$-dimensional vectors and an impurity measure $I$, the goal is to find a partition ${\cal P}$ of $V$ into $L$ groups $V_1,\ldots,V_L$ that minimizes the total impurities of the groups in ${\cal P}$, i.e., $I({\cal P})= \sum_{m=1}^{L} I(\sum_{{\bf v} \in V_m}{\bf v}).$ Impurity minimization is widely used as quality assessment measure in probability distribution clustering and in categorical clustering where it is not possible to rely on geometric properties of the data set. However, in contrast to the case of metric based clustering, the current knowledge of impurity measure based clustering in terms of approximation and inapproximability results is very limited. Our research contributes to fill this gap. We first present a simple linear time algorithm that simultaneously achieves $3$-approximation for the Gini impurity measure and $O(\log(\sum_{{\bf v} \in V} \| {\bf v} \|_1))$-approximation for the Entropy impurity measure. Then, for the Entropy impurity measure---where we also show that finding the optimal clustering is strongly NP-hard---we are able to design a polynomial time $O(\log^2(\min\{k,L\}))$-approximation algorithm. Our algorithm relies on a nontrivial characterization of a class of clusterings that necessarily includes a partition achieving $O(\log^2(\min\{k,L\}))$--approximation of the impurity of the optimal partition. Remarkably, this is the first polynomial time algorithm with approximation guarantee independent of the number of points/vector and not relying on any restriction on the components of the vectors for producing clusterings with minimum entropy.


title: GreyFiber: A System for Providing Flexible Access to Wide-Area Connectivity

abstract: Access to fiber-optic connectivity in the Internet is traditionally offered either via lit circuits or dark fiber. Economic (capex vs. opex) and operational considerations (latency, capacity) dictate the choice between these two offerings, but neither may effectively address the specific needs of modern-day enterprises or service providers over a range of use scenarios. In this paper, we describe a new approach for fiber-optic connectivity in the Internet that we call GreyFiber. The core idea of GreyFiber is to offer flexible access to fiber-optic paths between end points (e.g., datacenters or colocation facilities) over a range of timescales. We identify and discuss operational issues and systems challenges that need to be addressed to make GreyFiber a viable and realistic option for offering flexible access to infrastructure (similar to cloud computing). We investigate the efficacy of GreyFiber with a prototype implementation deployed in the GENI and CloudLab testbeds. Our scaling experiments show that 50 circuits can be provisioned within a minute. We also show that backup paths can be provisioned 28 times faster than an OSPF-based solution during failure/maintenance events. Our experiments also examine GreyFiber overhead demands and show that the time spent in circuit creation is dependent on the network infrastructure, indicating avenues for future improvements.


title: Probing the Berry Curvature and Fermi Arcs of a Weyl Circuit

abstract: The Weyl particle is the massless fermionic cousin of the photon. While no fundamental Weyl particles have been identified, they arise in condensed matter and meta-material systems, where their spinor nature imposes topological constraints on low-energy dispersion and surface properties. Here we demonstrate a topological circuit with Weyl dispersion at low-momentum, realizing a 3D lattice that behaves as a half-flux Hofstadter model in all principal planes. The circuit platform provides access to the complete complex-valued spin-texture of all bulk- and surface- states, thereby revealing not only the presence of Weyl points and the Fermi arcs that connect their surface-projections, but also, for the first time, the Berry curvature distribution through the Brillouin zone and the associated quantized Chiral charge of the Weyl points. This work opens a path to exploration of interacting Weyl physics in superconducting circuits, as well as studies of how manifold topology impacts band topology in three dimensions.


title: Extending the D-Wave with support for Higher Precision Coefficients

abstract: D-Wave only guarantees to support coefficients with 4 to 5 bits of resolution or precision. This paper describes a method to extend the functionality of the D-Wave to solve problems that require the support of higher precision coefficients.


title: Performance of Humans in Iris Recognition: The Impact of Iris Condition and Annotation-driven Verification

abstract: This paper advances the state of the art in human examination of iris images by (1) assessing the impact of different iris conditions in identity verification, and (2) introducing an annotation step that improves the accuracy of people's decisions. In a first experimental session, 114 subjects were asked to decide if pairs of iris images depict the same eye (genuine pairs) or two distinct eyes (impostor pairs). The image pairs sampled six conditions: (1) easy for algorithms to classify, (2) difficult for algorithms to classify, (3) large difference in pupil dilation, (4) disease-affected eyes, (5) identical twins, and (6) post-mortem samples. In a second session, 85 of the 114 subjects were asked to annotate matching and non-matching regions that supported their decisions. Subjects were allowed to change their initial classification as a result of the annotation process. Results suggest that: (a) people improve their identity verification accuracy when asked to annotate matching and non-matching regions between the pair of images, (b) images depicting the same eye with large difference in pupil dilation were the most challenging to subjects, but benefited well from the annotation-driven classification, (c) humans performed better than iris recognition algorithms when verifying genuine pairs of post-mortem and disease-affected eyes (i.e., samples showing deformations that go beyond the distortions of a healthy iris due to pupil dilation), and (d) annotation does not improve accuracy of analyzing images from identical twins, which remain confusing for people.


title: Derangements, Ehrhart Theory, and Local h-polynomials

abstract: The Eulerian polynomials and derangement polynomials are two well-studied generating functions that frequently arise in combinatorics, algebra, and geometry. When one makes an appearance, the other often does so as well, and their corresponding generalizations are similarly linked. This is this case in the theory of subdivisions of simplicial complexes, where the Eulerian polynomial is an $h$-polynomial and the derangement polynomial is its local $h$-polynomial. Separately, in Ehrhart theory the Eulerian polynomials are generalized by the $h^\ast$-polynomials of $s$-lecture hall simplices. Here, we show that derangement polynomials are analogously generalized by the box polynomials, or local $h^\ast$-polynomials, of the $s$-lecture hall simplices, and that these polynomials are all real-rooted. We then connect the two theories by showing that the local $h$-polynomials of common subdivisions in algebra and topology are realized as local $h^\ast$-polynomials of $s$-lecture hall simplices. We use this connection to address some open questions on real-rootedness and unimodality of generating polynomials, some from each side of the story.


title: Channel Charting: Locating Users within the Radio Environment using Channel State Information

abstract: We propose channel charting (CC), a novel framework in which a multi-antenna network element learns a chart of the radio geometry in its surrounding area. The channel chart captures the local spatial geometry of the area so that points that are close in space will also be close in the channel chart and vice versa. CC works in a fully unsupervised manner, i.e., learning is only based on channel state information (CSI) that is passively collected at a single point in space, but from multiple transmit locations in the area over time. The method then extracts channel features that characterize large-scale fading properties of the wireless channel. Finally, the channel charts are generated with tools from dimensionality reduction, manifold learning, and deep neural networks. The network element performing CC may be, for example, a multi-antenna base-station in a cellular system and the charted area in the served cell. Logical relationships related to the position and movement of a transmitter, e.g., a user equipment (UE), in the cell can then be directly deduced from comparing measured radio channel characteristics to the channel chart. The unsupervised nature of CC enables a range of new applications in UE localization, network planning, user scheduling, multipoint connectivity, hand-over, cell search, user grouping, and other cognitive tasks that rely on CSI and UE movement relative to the base-station, without the need of information from global navigation satellite systems.


title: Domain-Specific Human-Inspired Binarized Statistical Image Features for Iris Recognition

abstract: Binarized statistical image features (BSIF) have been successfully used for texture analysis in many computer vision tasks, including iris recognition and biometric presentation attack detection. One important point is that all applications of BSIF in iris recognition have used the original BSIF filters, which were trained on image patches extracted from natural images. This paper tests the question of whether domain-specific BSIF can give better performance than the default BSIF. The second important point is in the selection of image patches to use in training for BSIF. Can image patches derived from eye-tracking experiments, in which humans perform an iris recognition task, give better performance than random patches? Our results say that (1) domain-specific BSIF features can out-perform the default BSIF features, and (2) selecting image patches in a task-specific manner guided by human performance can out-perform selecting random patches. These results are important because BSIF is often regarded as a generic texture tool that does not need any domain adaptation, and human-task-guided selection of patches for training has never (to our knowledge) been done. This paper follows the reproducible research requirements, and the new iris-domain-specific BSIF filters, the patches used in filter training, the database used in testing and the source codes of the designed iris recognition method are made available along with this paper to facilitate applications of this concept.


title: VTXO - Virtual Telescope for X-Ray Observations

abstract: The Virtual Telescope for X-Ray Observations (VTXO) is a conceptual mission under development to demonstrate a new instrument for astronomical observations in the X-ray band using a set of 6U CubeSats. VTXO will use a Phase Fresnel Lens, which has the potential to improve the imaging resolution several orders of magnitude over the current state-of-the-art X-ray telescopes. This technique requires long focal lengths (>100 m), which necessitates the lens and camera be on separate spacecraft, flying in precision formation. This work presents the results from a model developed to determine the {\Delta}V requirements to maintain formation, for both solar and galactic X-ray observations, from a Geostationary Transfer Orbit.


title: The Spin Structure of the Nucleon

abstract: We review the present understanding of the spin structure of protons and neutrons, the fundamental building blocks of nuclei collectively known as nucleons. The field of nucleon spin provides a critical window for testing Quantum Chromodynamics (QCD), the gauge theory of the strong interactions since it involves fundamental aspects of hadron structure, and it can be probed in detail in experiments, particularly deep inelastic lepton scattering on polarized targets. QCD was initially probed in high energy deep inelastic lepton scattering with unpolarized beams and targets. With time, interest shifted from testing perturbative QCD to illuminating the nucleon structure itself. In fact, the spin degrees of freedom of hadrons provide an essential and detailed verification of both perturbative and nonperturbative QCD dynamics. Nucleon spin was initially thought of coming mostly from the spin of its quark constituents, based on intuition from the parton model. However, the first experiments showed that this expectation was incorrect. It is now clear that nucleon physics is much more complex, involving quark orbital angular momenta as well as gluonic and sea quark contributions. Thus, the nucleon spin structure remains a most active aspect of QCD research, involving important advances such as the developments of generalized parton distributions (GPD) and transverse momentum distributions (TMD). Elastic and inelastic lepton-proton scattering, as well as photoabsorption experiments provide various ways to investigate non-perturbative QCD. Fundamental sum rules -- such as the Bjorken sum rule for polarized photoabsorption on polarized nucleons -- are also in the non-perturbative domain. This realization triggered a vigorous program to link the low energy effective hadronic description of the strong interactions to fundamental quarks and gluon degrees of freedom of...


title: Drifting Photons on Optical Paths: Mirrors, Sub-mm Resolution in 4 Dimensions, and Transverse/Longitudinal Phase Space: Exploiting Time Resolution

abstract: I discuss the status of MCP-based photo-detector amplification sections and Cherenkov light sources for precise timing measurements of charged particles and gamma rays. Sub-psec resolution is predicted for the large pulses such as those produced by a charged particle or electromagnetic shower traversing a photo-detector entrance window. Measuring events with sub-mm resolution in each of the 4 dimensions expands the optical phase space from 4 dimensions, allowing emittance transformations that can minimize expensive instrumented photo-sensitive area.

categories: physics.ins-det astro-ph.IM hep-ex hep-ph nucl-ex


title: The Dune Python Module

abstract: In this paper we present the new Dune-Python module which provides Python bindings for the Dune core, which is a C++ environment for solving partial differential equations. The aim of this new module is to firstly provide the general infrastructure for exporting realizations of statically polymorphic interfaces based on just-in-time compilation and secondly to provide bindings for the central interfaces of the dune core modules. In the first release we focus on the grid interface. Our aim is to only introduce a thin layer when passing objects into Python which can be removed when the object is passed back into a C++ algorithm. Thus no efficiency is lost and little additional code maintenance cost is incurred. To make the transition for Dune users to the Python environment straightforward the Python classes provide a very similar interface to their C++ counterparts. In addition, vectorized versions of many interfaces allow for more efficient code on the Python side. The infrastructure for exporting these interfaces and the resulting bindings for a Dune grid are explained in detail in this paper for both experienced Dune users and others interested in a flexible Python environment for implementing grid based schemes for solving partial differential equations.


title: Neutron inelastic scattering measurements on $^{136}$Xe at $E_{n}$ = 0.7 to 100 MeV

abstract: Experiments searching for neutrinoless double beta decay ($0\nu\beta\beta$) require precise energy calibration and extremely low backgrounds. One of the most popular isotopes for $0\nu\beta\beta$ experiments is $^{136}$Xe. In support of these experiments, the neutron inelastic scattering properties of this isotope have been measured at the GErmanium Array for Neutron Induced Excitations (GEANIE) at the Los Alamos Neutron Science Center. Time-of-flight techniques are utilized with high-purity germanium detectors to search for inelastic scattering $\gamma$ rays for neutron energies between 0.7 and 100 MeV. Limits are set on production of yet-unobserved $\gamma$ rays in the energy range critical for $0\nu\beta\beta$ studies, and measurements are made of multiple $\gamma$-ray production cross sections. In particular, we have measured the production of the 1313 keV $\gamma$ ray which comes from the transition of the first-excited to ground state of $^{136}$Xe. This neutron-induced $\gamma$ line may be useful for a novel energy calibration technique, described in this paper.


title: Cyclotron Damping along an Uniform Magnetic Field

abstract: We prove cyclotron damping for the collisionless Vlasov-Maxwell equations on $\mathbb{T}_{x}^{3}\times\mathbb{R}_{v}^{3}$ under the assumptions that the electric induction is zero and $(\mathcal{\mathbf{PSC}})$ holds. It is a crucial step to solve the stability problem of the Vlasov-Maxwell equations. Our proof is based on a new dynamical system of the plasma particles, originating from Faraday Law of Electromagnetic induction and Lenz's Law. On the basis of it, we use the improved Newton iteration scheme to show the damping mechanism.


title: Extremal primes of elliptic curves without complex multiplication

abstract: Fix an elliptic curve E over Q. An extremal prime for E is a prime p of good reduction such that the number of rational points on E modulo p is maximal or minimal in relation to the Hasse bound. Assuming that all the symmetric power L-functions associated to E are automorphic and satisfy the Generalized Riemann Hypothesis, we give the first non-trivial upper bounds for the number of such primes when E is a curve without complex multiplication. In order to obtain this bound, we use explicit equidistribution for the Sato-Tate measure as in the work of Rouse and Thorner (arXiv:1305.5283) and refine certain intermediate estimates taking advantage of the fact that extremal primes have a very small Sato-Tate measure.


title: A one-variable bracket polynomial for some Turk's head knots

abstract: We compute the Kauffman bracket polynomial of the three-lead Turk's head, the chain sinnet and the figure-eight chain shadow diagrams. Each of these knots can in fact be constructed by repeatedly concatenating the same 3-tangle, respectively, then taking the closure. The bracket is then evaluated by expressing the state diagrams of the concerned 3-tangle by means of the Kauffman monoid diagram's elements.


title: Complete monotonicity properties of a function involving the polygamma function

abstract: In this paper, we study completete monotonicity properties of the function $f_{a,k}(x)=\psi^{(k)}(x+a) - \psi^{(k)}(x) - \frac{ak!}{x^{k+1}}$, where $a\in(0,1)$ and $k\in \mathbb{N}_0$. Specifically, we consider the cases for $k\in \{ 2n: n\in \mathbb{N}_0 \}$ and $k\in \{ 2n+1: n\in \mathbb{N}_0 \}$. Subsequently, we deduce some inequalities involving the polygamma functions.


title: QR2: A Third-party Query Reranking Service Over Web Databases

abstract: The ranked retrieval model has rapidly become the de-facto way for search query processing in web databases. Despite the extensive efforts on designing better ranking mechanisms, in practice, many such databases fail to address the diverse and sometimes contradicting preferences of users. In this paper, we present QR2, a third-party service that uses nothing but the public search interface of a web database and enables the on-the-fly processing of queries with any user-specified ranking functions, no matter if the ranking function is supported by the database or not.


title: Dynamical Instability of Spherical Anisotropic Sources in $f(R,T,R_{\mu\nu}T^{\mu\nu})$ Gravity

abstract: In this paper, we study the effects of modification of gravity on the problem of dynamical instability of the spherical relativistic anisotropic interiors. We have considered non-zero influence of expansion scalar throughout during the evolutionary phases of spherical geometry that led to the use of fluid stiffness parameter. The modified hydrostatic equation for the stellar anisotropic matter distributions is constructed and then solved by using radial perturbation scheme. Such a differential equation can be further used to obtain instability constraints at both weak field and post-Newtonian approximations after considering a particular Harrison-Wheeler equation of state. This approach allows us to deal with the effects of usual and effective matter variables on the stability exotic stellar of self-gravitating structures.


title: Are NLS1s highly accreting low black hole mass AGNs?

abstract: In this work, we test the hypothesis that narrow-line Seyfert 1 galaxies (NSL1s) are active galactic nuclei in their early phase and are therefore younger and more active than the more common broad-line Seyfert 1 galaxies (BLS1s). If that is true, then NLS1s should, on average, have lower black hole masses and higher accretion rates than BLS1s. To test this, we use a sample of 35 NLS1s and 54 BLS1s with similar X-ray luminosity distributions and good XMM-Newton observations. To determine the black hole mass, we apply an X-ray scaling method that is independent of any assumptions on the broad-line region dynamics and the inclination of the objects. We find that, on average, NLS1s have lower black hole masses, but the difference between the average black hole masses of NLS1s and BLS1s in our sample is only marginally significant (at the 2.6 sigma level). According to a Kolmogorov-Smirnov test, the distribution of black hole mass values of NLS1s is different from that of BLS1s at the 99% confidence level. Even stronger differences between NLS1s and BLS1s are inferred when the accretion rate distributions of NLS1s are compared to BLS1s, suggesting that the two populations are indeed distinct. Our study also indicates that the black hole mass values (both for NLS1s and BLS1s) determined with the X-ray scaling method are fully consistent with those obtained using reverberation mapping.


title: Fulde-Ferrell-Larkin-Ovchinnikov superconductors near a surface

abstract: We show that the behaviour of the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) superconductors near a surface is considerably different from the usual case. The order parameter of the FF state is strongly deformed near the surface, which leads to a number of unusual features in the linear magnetic response, such as "anti-screening" or "over-screening" of the applied field. In a fully isotropic FF case, the Meissner effect is still present, despite the vanishing of the transverse superfluid density in the bulk. We also calculate the surface critical field Hc3, which exhibits a peculiar temperature dependence.


title: Analysing the role of entanglement in the three-qubit Vaidman's game

abstract: We analyse the role of degree of entanglement for Vaidman's game in a setting where the players share a set of partially entangled three-qubit states. Our results show that the entangled states combined with quantum strategies may not be always helpful in winning a game as opposed to the classical strategies. We further find the conditions under which quantum strategies are always helpful in achieving higher winning probability in the game in comparison to classical strategies. Moreover, we show that a special class of W states can always be used to win the game using quantum strategies irrespective of the degree of entanglement between the three qubits. Our analysis also helps us in comparing the Vaidman's game with the secret sharing protocol. Furthermore, we propose a new Vaidman-type game where the rule maker itself is entangled with the other two players and acts as a facilitator to share a secret key with the two players.


title: On local presentability of T/A

abstract: We prove that if $\mathcal{A}$ is a locally $\lambda$-presentable category and $T : \mathcal{A} \to \mathcal{A}$ is a $\lambda$-accessible functor then $T/\mathcal{A}$ is locally $\lambda$-presentable.


title: Stabilization and control for the biharmonic Schr\"odinger equation

abstract: The main purpose of this paper is to show the global stabilization and exact controllability properties for a fourth order nonlinear fourth order nonlinear Schr\"odinger system: $$i\partial_tu +\partial_x^2u-\partial_x^4u=\lambda |u|^2u,$$ on a periodic domain $\mathbb{T}$ with internal control supported on an arbitrary sub-domain of $\mathbb{T}$. More precisely, by certain properties of propagation of compactness and regularity in Bourgain spaces, for the solutions of the associated linear system, we show that the system is globally exponentially stabilizable. This property together with the local exact controllability ensures that fourth order nonlinear Schr\"odinger is globally exactly controllable.


title: Rebalancing Frequency Considerations for Kelly-Optimal Stock Portfolios in a Control-Theoretic Framework

abstract: In this paper, motivated by the celebrated work of Kelly, we consider the problem of portfolio weight selection to maximize expected logarithmic growth. Going beyond existing literature, our focal point here is the rebalancing frequency which we include as an additional parameter in our analysis. The problem is first set in a control-theoretic framework, and then, the main question we address is as follows: In the absence of transaction costs, does high-frequency trading always lead to the best performance? Related to this is our prior work on betting, also in the Kelly context, which examines the impact of making a wager and letting it ride. Our results on betting frequency can be interpreted in the context of weight selection for a two-asset portfolio consisting of one risky asset and one riskless asset. With regard to the question above, our prior results indicate that it is often the case that there are no performance benefits associated with high-frequency trading. In the present paper, we generalize the analysis to portfolios with multiple risky assets. We show that if there is an asset satisfying a new condition which we call dominance, then an optimal portfolio consists of this asset alone; i.e., the trader has "all eggs in one basket" and performance becomes a constant function of rebalancing frequency. Said another way, the problem of rebalancing is rendered moot. The paper also includes simulations which address practical considerations associated with real stock prices and the dominant asset condition.


title: Second Law and Non-Equilibrium Entropy of Schottky Systems -- Doubts and Verification

abstract: Meixner's historical remark in 1969 "... it can be shown that the concept of entropy in the absence of equilibrium is in fact not only questionable but that it cannot even be defined...." is investigated from today's insight. Several statements --such as the three laws of phenomenological thermodynamics, the embedding theorem and the adiabatical uniqueness-- are used to get rid of non-equilibrium entropy as a primitive concept. In this framework, Clausius inequality of open systems can be derived by use of the defining inequalities which establish the non-equilibrium quantities contact temperature and non-equilibrium molar entropy which allow to describe the interaction between the Schottky system and its controlling equilibrium environment.


title: Photosynthesis on a planet orbiting an M dwarf: enhanced effectiveness during flares

abstract: On planets near M dwarfs, photosynthesis (PS) will occur with an effectiveness which depends on the supply of visible photons with wavelengths between 400 and 700 nm. In this paper, we quantify the effectiveness of PS in two contexts which are relevant for M dwarfs. First, using photons from an M dwarf in its quiescent non-flaring state, we find that PS on an M dwarf planet in the HZ of its parent star is less effective than on Earth by a factor of 10 for a flare star with mid-M spectral type. For a flare star with late-M spectral type, PS effectiveness is smaller than on Earth by a factor of 100 or more. Second, using photons which are incident on the HZ planet during flares, we find that PS effectiveness can increase by factors of 5-20 above the quiescent values. In the case of a flare star with mid-M spectral type, we find that the PS effectiveness during a flare can increase up to as much as 50-60 percent of the values on Earth. However, for a late-M flare star, even during flares, the PS effectiveness remains almost one order of magnitude smaller than on Earth. We suggest that for biological processes on M dwarf planets, the stellar activity cycle may replace the orbital period as the year.


title: Tunable Bloch surface waves devices

abstract: This thesis is devoted to develop tunable devices on the base of one-dimensional photonic crystals (1DPhC) which can sustain Bloch surface waves (BSWs). First, we explore the possibilities to control the BSW propagation direction with polarization of incident light. In this case we manufacture additional passive structures such as gratings on the top of the 1DPhC, which are working both as a BSW launcher and polarization-controlled wave-splitters. We test this type of launcher in air and in water as an external medium. Then, we demonstrate the tunability of the BSW by adding an active layers into the multilayer stack. Here a crystalline X-cut thin film lithium niobate (TFLN) is used to introduce anisotropic properties to the whole 1DPhC. Different ways to manufacture 1D PhCs with LiNbO3 on the top would be described. Finally, we explore the concept of the electro-optically tuned BSW.


title: Non-commutative inspired black holes in Euler-Heisenberg non-linear electrodynamics

abstract: We find non-commutative inspired electrically and magnetically charged black hole solutions in Euler-Heisenberg non-linear electrodynamics. For these solutions, we determine the non-commutative corrections to the horizon radius for the general and extremal case. We also analyse the weak, dominant and strong energy conditions and the shadow associated with these metrics.


title: Curves, Jacobians, and Cryptography

abstract: The main purpose of this paper is to give an overview over the theory of abelian varieties, with main focus on Jacobian varieties of curves reaching from well-known results till to latest developments and their usage in cryptography. In the first part we provide the necessary mathematical background on abelian varieties, their torsion points, Honda-Tate theory, Galois representations, with emphasis on Jacobian varieties and hyperelliptic Jacobians. In the second part we focus on applications of abelian varieties on cryptography and treating separately, elliptic curve cryptography, genus 2 and 3 cryptography, including Diffie-Hellman Key Exchange, index calculus in Picard groups, isogenies of Jacobians via correspondences and applications to discrete logarithms. Several open problems and new directions are suggested.


title: Microscopic collectivity: The ridge and strangeness enhancement from string-string interactions

abstract: We present the microscopic model for collective effects, built on string-string interactions (string shoving and rope hadronization), recently implemented in the PYTHIA 8 and DIPSY event generators. Rope hadronization is shown to give a good description of strangeness enhancement across pp, pA and AA collisions systems, while string shoving qualitatively describes the ridge observed in pp collisions.


title: Galoisian and Qualitative Approaches to Linear Polyanin-Zaitsev Vector Fields

abstract: The analysis of dynamical systems has been a topic of great interest for researches mathematical sciences for a long times. The implementation of several devices and tools have been useful in the finding of solutions as well to describe common behaviors of parametric families of these systems. In this paper we study deeply a particular parametric family of differential equations, the so-called \emph{Linear Polyanin-Zaitsev Vector Field}, which has been introduced in a general case, in the previous paper by the same authors, as a correction of a family presented in the classical book of Polyanin and Zaitsev. Linear Polyanin-Zaitsev Vector Field is transformed into a Li\'enard equation and in particular we obtain the Van Der Pol equation. We present some algebraic and qualitative results to illustrate some interactions between algebra and the qualitative theory of differential equations in this parametric family.


title: High performing ionanofluid electrolyte with higher lithium salt concentration for safer lithium metal batteries

abstract: A new class of high-performance pyrrolidinium cation based ionanofluid electrolytes with higher lithium salt concentration are developed. The electrolytes are formed by dispersing imidazolium ionic liquid functionalized TiO2 nanoparticles in low conducting, 0.6 M lithium salt doped N-alkyl-N-butylpyrrolidinium bis(trifluoromethylsulfonyl)imide (Pyr14TFSI) ionic liquid (IL) hosted electrolyte. Viscosity, ionic conductivity and thermal properties of these electrolytes are compared with well studied 0.2 M salt doped Pyr14TFSI IL-based electrolyte. The highly crystalline 0.6 M lithium salt dissolved IL-based electrolytes gradually become amorphous with the increasing dispersion of surface functionalized nanoparticles within it. The ionic conductivity of the electrolytes shows unusual viscosity decoupled characteristics and at the 5.0 wt% nanoparticle dispersion it attains a maximum value, higher than that of pure IL host. As compared to pure IL-based electrolytes, the ionanofluid electrolyte also possesses a significantly higher value of lithium ion transference number. The Li/LiMn2O4 cell with the best conducting ionanofluid electrolyte delivers a discharge capacity of about 131 mAh g-1 at 25 degree C at a current density of 24 mA g-1, much higher than that obtained in 0.2 M Li salt dissociated Pyr14TFSI electrolyte (87 mAh g-1). Superior interfacial compatibility between ionanofluid electrolyte and electrodes as indicated by the excellent rate performance with outstanding capacity retention of the cell as compared to pure IL-based analogue, further establish great application potentiality of this optimized newly developed electrolyte for safer LMBs.


title: Sparse semiparametric canonical correlation analysis for data of mixed types

abstract: Canonical correlation analysis investigates linear relationships between two sets of variables, but often works poorly on modern data sets due to high-dimensionality and mixed data types such as continuous, binary and zero-inflated. To overcome these challenges, we propose a semiparametric approach for sparse canonical correlation analysis based on Gaussian copula. Our main contribution is a truncated latent Gaussian copula model for data with excess zeros, which allows us to derive a rank-based estimator of the latent correlation matrix for mixed variable types without the estimation of marginal transformation functions. The resulting canonical correlation analysis method works well in high-dimensional settings as demonstrated via numerical studies, as well as in application to the analysis of association between gene expression and micro RNA data of breast cancer patients.


title: LSTM-Based Zero-Velocity Detection for Robust Inertial Navigation

abstract: We present a method to improve the accuracy of a zero-velocity-aided inertial navigation system (INS) by replacing the standard zero-velocity detector with a long short-term memory (LSTM) neural network. While existing threshold-based zero-velocity detectors are not robust to varying motion types, our learned model accurately detects stationary periods of the inertial measurement unit (IMU) despite changes in the motion of the user. Upon detection, zero-velocity pseudo-measurements are fused with a dead reckoning motion model in an extended Kalman filter (EKF). We demonstrate that our LSTM-based zero-velocity detector, used within a zero-velocity-aided INS, improves zero-velocity detection during human localization tasks. Consequently, localization accuracy is also improved. Our system is evaluated on more than 7.5 km of indoor pedestrian locomotion data, acquired from five different subjects. We show that 3D positioning error is reduced by over 34% compared to existing fixed-threshold zero-velocity detectors for walking, running, and stair climbing motions. Additionally, we demonstrate how our learned zero-velocity detector operates effectively during crawling and ladder climbing. Our system is calibration-free (no careful threshold-tuning is required) and operates consistently with differing users, IMU placements, and shoe types, while being compatible with any generic zero-velocity-aided INS.


title: Robust Chauvenet Outlier Rejection

abstract: Sigma clipping is commonly used in astronomy for outlier rejection, but the number of standard deviations beyond which one should clip data from a sample ultimately depends on the size of the sample. Chauvenet rejection is one of the oldest, and simplest, ways to account for this, but, like sigma clipping, depends on the sample's mean and standard deviation, neither of which are robust quantities: Both are easily contaminated by the very outliers they are being used to reject. Many, more robust measures of central tendency, and of sample deviation, exist, but each has a tradeoff with precision. Here, we demonstrate that outlier rejection can be both very robust and very precise if decreasingly robust but increasingly precise techniques are applied in sequence. To this end, we present a variation on Chauvenet rejection that we call "robust" Chauvenet rejection (RCR), which uses three decreasingly robust/increasingly precise measures of central tendency, and four decreasingly robust/increasingly precise measures of sample deviation. We show this sequential approach to be very effective for a wide variety of contaminant types, even when a significant -- even dominant -- fraction of the sample is contaminated, and especially when the contaminants are strong. Furthermore, we have developed a bulk-rejection variant, to significantly decrease computing times, and RCR can be applied both to weighted data, and when fitting parameterized models to data. We present aperture photometry in a contaminated, crowded field as an example. RCR may be used by anyone at https://skynet.unc.edu/rcr, and source code is available there as well.


title: On an observation of Sibony

abstract: It is shown that if the boundary of a Reinhardt domain in $\mathbb{C}^n$ contains the origin, each holomorphic function on the domain which is infinitely many times differentiable up to the boundary extends holomorphically to a neighborhood of the origin.


title: The James Webb Space Telescope North Ecliptic Pole Time-Domain Field -- I: Field Selection of a JWST Community Field for Time-Domain Studies

abstract: We describe the selection of the James Webb Space Telescope (JWST) North Ecliptic Pole (NEP) Time-Domain Field (TDF), a ~14' diameter field located within JWST's northern Continuous Viewing Zone (CVZ) and centered at (RA, Dec)_J2000 = (17:22:47.896, +65:49:21.54). We demonstrate that this is the only region in the sky where JWST can observe a clean (i.e., free of bright foreground stars and with low Galactic foreground extinction) extragalactic deep survey field of this size at arbitrary cadence or at arbitrary orientation, and without a penalty in terms of a raised Zodiacal background. This will crucially enable a wide range of new and exciting time-domain science, including high redshift transient searches and monitoring (e.g., SNe), variability studies from Active Galactic Nuclei (AGN) to brown dwarf atmospheres, as well as proper motions of possibly extreme scattered Kuiper Belt and Inner Oort Cloud Objects, and of nearby Galactic brown dwarfs, low-mass stars, and ultracool white dwarfs. A JWST/NIRCam+NIRISS GTO program will provide an initial 0.8--5.0micron spectrophotometric characterization to m_AB ~ 28.8+/-0.3 mag of four orthogonal "spokes" within this field. The multi-wavelength (radio through X-ray) context of the field is in hand (ground-based near-UV--visible--near-IR), in progress (VLA 3GHz, VLBA 5GHz, HST UV--visible, Chandra X-ray, IRAM30m 1.3 and 2mm), or scheduled (JCMT 850micron). We welcome and encourage ground- and space-based follow-up of the initial GTO observations and ancillary data, to realize its potential as an ideal JWST time-domain community field.


title: The True Luminosities of Planetary Nebulae in M31's Bulge: Massive Central Stars from an Old Stellar Population

abstract: We measure the Balmer decrements of 23 of the brightest planetary nebulae (PNe) in the inner bulge ($r \lesssim 3$ arcmin) of M31 and de-redden the bright end of the region's [O III] $\lambda 5007$ planetary nebula luminosity function. We show that the most luminous PNe produce $\gtrsim 1{,}200 \, \rm{L}_{\odot}$ of power in their [O III] $\lambda 5007$ line, implying central star luminosities of at least $\sim 11{,}000 \, \rm{L}_{\odot}$. Even with the most recent accelerated-evolution post-AGB models, such luminosities require central star masses in excess of $0.66 \, \rm{M}_{\odot}$, and main sequence progenitors of at least $\sim 2.5 \, \rm{M}_{\odot}$. Since M31's bulge has very few intermediate-age stars, we conclude that conventional single-star evolution cannot be responsible for these extremely luminous objects. We also present the circumstellar extinctions for the region's bright PNe and demonstrate that the distribution is similar to that found for PNe in the Large Magellanic Cloud, with a median value of $A_{5007} = 0.71$. Finally, we compare our results to extinction measurements made for PNe in the E6 elliptical NGC 4697 and the interacting lenticular NGC 5128. We show that such extinctions are not unusual, and that the existence of very high-mass PN central stars is a general feature of old stellar populations. Our results suggest that single-star population synthesis models significantly underestimate the maximum luminosities and total integrated light of AGB stars.


title: Maximal Operators for cube skeletons

abstract: We study discretized maximal operators associated to averaging over (neighborhoods of) squares in the plane and, more generally, $k$-skeletons in $\mathbb{R}^n$. Although these operators are known not to be bounded on any $L^p$, we obtain nearly sharp $L^p$ bounds for every small discretization scale. These results are motivated by, and partially extend, recent results of T. Keleti, D. Nagy and P. Shmerkin, and of R. Thornton, on sets that contain a scaled $k$-sekeleton of the unit cube with center in every point of $\mathbb{R}^n$.


title: Numerical calculation of high-order QED contributions to the electron anomalous magnetic moment

abstract: This paper describes a method of numerical evaluating high-order QED contributions to the electron anomalous magnetic moment. The method is based on subtraction of infrared and ultraviolet divergences in Feynman-parametric space before integration and on nonadaptive Monte Carlo integration that is founded on Hepp sectors. A realization of the method on the graphics accelerator NVidia Tesla K80 is described. A method of removing round-off errors that emerge due to numerical subtraction of divergences without losing calculation speed is presented. The results of applying the method to all 2-loop, 3-loop, 4-loop QED Feynman graphs without lepton loops are presented. A detailed comparison of the 2-loop and 3-loop results with known analytical ones is given in the paper. A comparison of the contributions of 6 gauge invariant 4-loop graph classes with known analytical values is presented. Moreover, the contributions of 78 sets of 4-loop graphs for comparison with the direct subtraction on the mass shell are presented. Also, the contributions of the 5-loop and 6-loop ladder graphs are given as well as a comparison of these results with known analytical ones. The behavior of the generated Monte Carlo samples is described in detail, a method of the error estimation is presented. A detailed information about the graphics processor performance on these computations and about the Monte Carlo convergence is given in the paper.


title: The High Definition X-ray Imager (HDXI) Instrument on the Lynx X-Ray Surveyor

abstract: The Lynx X-ray Surveyor Mission is one of 4 large missions being studied by NASA Science and Technology Definition Teams as mission concepts to be evaluated by the upcoming 2020 Decadal Survey. By utilizing optics that couple fine angular resolution (<0.5 arcsec HPD) with large effective area (~2 m^2 at 1 keV), Lynx would enable exploration within a unique scientific parameter space. One of the primary soft X-ray imaging instruments being baselined for this mission concept is the High Definition X-ray Imager, HDXI. This instrument would achieve fine angular resolution imaging over a wide field of view (~ 22 x 22 arcmin, or larger) by using a finely-pixelated silicon sensor array. Silicon sensors enable large-format/small-pixel devices, radiation tolerant designs, and high quantum efficiency across the entire soft X-ray bandpass. To fully exploit the large collecting area of Lynx (~30x Chandra), without X-ray event pile-up, the HDXI will be capable of much faster frame rates than current X-ray imagers. The planned requirements, capabilities, and development status of the HDXI will be described.


title: When Are Two Gossips the Same? Types of Communication in Epistemic Gossip Protocols

abstract: We provide an in-depth study of the knowledge-theoretic aspects of communication in so-called gossip protocols. Pairs of agents communicate by means of calls in order to spread information---so-called secrets---within the group. Depending on the nature of such calls knowledge spreads in different ways within the group. Systematizing existing literature, we identify 18 different types of communication, and model their epistemic effects through corresponding indistinguishability relations. We then provide a classification of these relations and show its usefulness for an epistemic analysis in presence of different communication types. Finally, we explain how to formalise the assumption that the agents have common knowledge of a distributed epistemic gossip protocol.


title: Survey on Deep Learning Techniques for Person Re-Identification Task

abstract: Intelligent video-surveillance is currently an active research field in computer vision and machine learning techniques. It provides useful tools for surveillance operators and forensic video investigators. Person re-identification (PReID) is one among these tools. It consists of recognizing whether an individual has already been observed over a camera in a network or not. This tool can also be employed in various possible applications such as off-line retrieval of all the video-sequences showing an individual of interest whose image is given a query, and online pedestrian tracking over multiple camera views. To this aim, many techniques have been proposed to increase the performance of PReID. Among the systems, many researchers utilized deep neural networks (DNNs) because of their better performance and fast execution at test time. Our objective is to provide for future researchers the work being done on PReID to date. Therefore, we summarized state-of-the-art DNN models being used for this task. A brief description of each model along with their evaluation on a set of benchmark datasets is given. Finally, a detailed comparison is provided among these models followed by some limitations that can work as guidelines for future research.


title: Voltage Control of a van der Waals Spin-Filter Magnetic Tunnel Junction

abstract: Atomically thin chromium triiodide (CrI3) has recently been identified as a layered antiferromagnetic insulator, in which adjacent ferromagnetic monolayers are antiferromagnetically coupled. This unusual magnetic structure naturally comprises a series of anti-aligned spin filters which can be utilized to make spin-filter magnetic tunnel junctions with very large tunneling magnetoresistance (TMR). Here we report voltage control of TMR formed by four-layer CrI3 sandwiched by monolayer graphene contacts in a dual-gated structure. By varying the gate voltages at fixed magnetic field, the device can be switched reversibly between bistable magnetic states with the same net magnetization but drastically different resistance (by a factor of ten or more). In addition, without switching the state, the TMR can be continuously modulated between 17,000% and 57,000%, due to the combination of spin-dependent tunnel barrier with changing carrier distributions in the graphene contacts. Our work demonstrates new kinds of magnetically moderated transistor action and opens up possibilities for voltage-controlled van der Waals spintronic devices.


title: Schottky presentations of positive representations

abstract: We show that the notion of $3$-hyperconvexity on oriented flag manifolds defines a partial cyclic order. Using the notion of interval given by this partial cyclic order, we construct Schottky groups and show that they correspond to images of positive representations in the sense of Fock and Goncharov. We construct polyhedral fundamental domains for the domain of discontinuity that these groups admit in the projective space or the sphere, depending on the dimension.


title: The $\rho$ parameter at three loops and elliptic integrals

abstract: We describe the analytic calculation of the master integrals required to compute the two-mass three-loop corrections to the $\rho$ parameter. In particular, we present the calculation of the master integrals for which the corresponding differential equations do not factorize to first order. The homogeneous solutions to these differential equations are obtained in terms of hypergeometric functions at rational argument. These hypergeometric functions can further be mapped to complete elliptic integrals, and the inhomogeneous solutions are expressed in terms of a new class of integrals of combined iterative non-iterative nature.


title: Dark sector portal with vector-like leptons and flavor sequestering

abstract: We consider models with fermionic dark matter that transforms under a non-Abelian dark gauge group. Exotic, vector-like leptons that also transform under the dark gauge group can mix with standard model leptons after spontaneous symmetry breaking and serve as a portal between the dark and visible sectors. We show in an explicit, renormalizable model based on a dark SU(2) gauge group how this can lead to adequate dark matter annihilation to a standard model lepton flavor so that the correct relic density is obtained. We identify a discrete symmetry that allows mass mixing between the vector-like fermions and a single standard model lepton flavor, while preventing mixing between these fields and the remaining standard model leptons. This flavor sequestering avoids unwanted lepton-flavor-violating effects, substantially relaxing constraints on the mass scale of the vector-like states. We discuss aspects of the phenomenology of the model, including direct detection of the dark matter.


title: Transfer Learning for High-Precision Trajectory Tracking Through $\mathcal{L}_1$ Adaptive Feedback and Iterative Learning

abstract: Robust and adaptive control strategies are needed when robots or automated systems are introduced to unknown and dynamic environments where they are required to cope with disturbances, unmodeled dynamics, and parametric uncertainties. In this paper, we demonstrate the capabilities of a combined $\mathcal{L}_1$ adaptive control and iterative learning control (ILC) framework to achieve high-precision trajectory tracking in the presence of unknown and changing disturbances. The $\mathcal{L}_1$ adaptive controller makes the system behave close to a reference model; however, it does not guarantee that perfect trajectory tracking is achieved, while ILC improves trajectory tracking performance based on previous iterations. The combined framework in this paper uses $\mathcal{L}_1$ adaptive control as an underlying controller that achieves a robust and repeatable behavior, while the ILC acts as a high-level adaptation scheme that mainly compensates for systematic tracking errors. We illustrate that this framework enables transfer learning between dynamically different systems, where learned experience of one system can be shown to be beneficial for another different system. Experimental results with two different quadrotors show the superior performance of the combined $\mathcal{L}_1$-ILC framework compared with approaches using ILC with an underlying proportional-derivative controller or proportional-integral-derivative controller. Results highlight that our $\mathcal{L}_1$-ILC framework can achieve high-precision trajectory tracking when unknown and changing disturbances are present and can achieve transfer of learned experience between dynamically different systems. Moreover, our approach is able to achieve precise trajectory tracking in the first attempt when the initial input is generated based on the reference model of the adaptive controller.


title: Adaptive Model Predictive Control for High-Accuracy Trajectory Tracking in Changing Conditions

abstract: Robots and automated systems are increasingly being introduced to unknown and dynamic environments where they are required to handle disturbances, unmodeled dynamics, and parametric uncertainties. Robust and adaptive control strategies are required to achieve high performance in these dynamic environments. In this paper, we propose a novel adaptive model predictive controller that combines model predictive control (MPC) with an underlying $\mathcal{L}_1$ adaptive controller to improve trajectory tracking of a system subject to unknown and changing disturbances. The $\mathcal{L}_1$ adaptive controller forces the system to behave in a predefined way, as specified by a reference model. A higher-level model predictive controller then uses this reference model to calculate the optimal reference input based on a cost function, while taking into account input and state constraints. We focus on the experimental validation of the proposed approach and demonstrate its effectiveness in experiments on a quadrotor. We show that the proposed approach has a lower trajectory tracking error compared to non-predictive, adaptive approaches and a predictive, non-adaptive approach, even when external wind disturbances are applied.


title: LASSO: Large Adaptive optics Survey for Substellar Objects using the new SAPHIRA detector on Robo-AO

abstract: We report on initial results from the largest infrared AO direct imaging survey searching for wide orbit (>100 AU) massive exoplanets and brown dwarfs as companions around young nearby stars using Robo-AO at the 2.1-m telescope on Kitt Peak, Arizona. The occurrence rates of these rare substellar companions are critical to furthering our understanding of the origin of planetary-mass companions on wide orbits. The observing efficiency of Robo-AO allows us to conduct a survey an order of magnitude larger than previously possible. We commissioned a low-noise high-speed SAPHIRA near-infrared camera to conduct this survey and report on its sensitivity, performance, and data reduction process.


title: Neural Networks Regularization Through Representation Learning

abstract: Neural network models and deep models are one of the leading and state of the art models in machine learning. Most successful deep neural models are the ones with many layers which highly increases their number of parameters. Training such models requires a large number of training samples which is not always available. One of the fundamental issues in neural networks is overfitting which is the issue tackled in this thesis. Such problem often occurs when the training of large models is performed using few training samples. Many approaches have been proposed to prevent the network from overfitting and improve its generalization performance such as data augmentation, early stopping, parameters sharing, unsupervised learning, dropout, batch normalization, etc. In this thesis, we tackle the neural network overfitting issue from a representation learning perspective by considering the situation where few training samples are available which is the case of many real world applications. We propose three contributions. The first one presented in chapter 2 is dedicated to dealing with structured output problems to perform multivariate regression when the output variable y contains structural dependencies between its components. The second contribution described in chapter 3 deals with the classification task where we propose to exploit prior knowledge about the internal representation of the hidden layers in neural networks. Our last contribution presented in chapter 4 showed the interest of transfer learning in applications where only few samples are available. In this contribution, we provide an automatic system based on such learning scheme with an application to medical domain. In this application, the task consists in localizing the third lumbar vertebra in a 3D CT scan. This work has been done in collaboration with the clinic Rouen Henri Becquerel Center who provided us with data.


title: Markets Beyond Nash Welfare for Leontief Utilities

abstract: We study the allocation of divisible goods to competing agents via a market mechanism, focusing on agents with Leontief utilities. The majority of the economics and mechanism design literature has focused on \emph{linear} prices, meaning that the cost of a good is proportional to the quantity purchased. Equilibria for linear prices are known to be exactly the maximum Nash welfare allocations. \emph{Price curves} allow the cost of a good to be any (increasing) function of the quantity purchased. We show that price curve equilibria are not limited to maximum Nash welfare allocations with two main results. First, we show that an allocation can be supported by strictly increasing price curves if and only if it is \emph{group-domination-free}. A similarly characterization holds for weakly increasing price curves. We use this to show that given any allocation, we can compute strictly (or weakly) increasing price curves that support it (or show that none exist) in polynomial time. These results involve a connection to the \emph{agent-order matrix} of an allocation, which may have other applications. Second, we use duality to show that in the bandwidth allocation setting, any allocation maximizing a CES welfare function can be supported by price curves.


title: Weight distributions, zeta functions and Riemann hypothesis for linear and algebraic geometry codes

abstract: This is a survey on weight enumerators, zeta functions and Riemann hypothesis for linear and algebraic-geometry codes.


title: Non-abelian lattice gauge theory with a topological action

abstract: SU(2) gauge theory is investigated with a lattice action which is insensitive to small perturbations of the lattice gauge fields. Bare perturbation theory can not be defined for such actions at all. We compare non-perturbative continuum results with that obtained by the usual Wilson plaquette action. The compared observables span a wide range of interesting phenomena: zero temperature large volume behavior (topological susceptibility), finite temperature phase transition (critical exponents and critical temperature) and also the small volume regime (discrete beta-function or step-scaling function). In the continuum limit perfect agreement is found indicating that universality holds for these topological lattice actions as well.


title: Efficient Distribution Estimation and Uncertainty Quantification for Elliptic Problems on Domains with Stochastic Boundaries

abstract: We study the problem of uncertainty quantification for the numerical solution of elliptic partial differential equation boundary value problems posed on domains with stochastically varying boundaries. We also use the uncertainty quantification results to tackle the efficient solution of such problems. We introduce simple transformations that map a family of domains with stochastic boundaries to a fixed reference domain. We exploit the transformations to carry out a prior and a posteriori error analyses and to derive an efficient Monte Carlo sampling procedure.


title: On Lusztig-Dupont homology of flag complexes

abstract: Let $V$ be an $n$-dimensional vector space over the finite field of order $q$. The spherical building $X_V$ associated with $GL(V)$ is the order complex of the nontrivial linear subspaces of $V$. Let $\mathfrak{g}$ be the local coefficient system on $X_V$, whose value on the simplex $\sigma=[V_0 \subset \cdots \subset V_p] \in X_V$ is given by $\mathfrak{g}(\sigma)=V_0$. Following the work of Lusztig and Dupont, we study the homology module $D^k(V)=\tilde{H}_{n-k-1}(X_V;\mathfrak{g})$. Our results include a construction of an explicit basis of $D^1(V)$, and the following twisted analogue of a result of Smith and Yoshiara: For any $1 \leq k \leq n-1$, the minimal support size of a non-zero $(n-k-1)$-cycle in the twisted homology $\tilde{H}_{n-k-1}(X_V;\wedge^k \mathfrak{g})$ is $\frac{(n-k+2)!}{2}$.


title: Numerical Simulations of Polymer Flooding Process in Porous Media on Distributed-memory Parallel Computers

abstract: Polymer flooding is a mature enhanced oil recovery technique that has been successfully applied in many field projects. By injecting polymer into a reservoir, the viscosity of water is increased, and the efficiency of water flooding is improved. As a result, more oil can be recovered. This paper presents numerical simulations of a polymer flooding process using parallel computers, where the numerical modeling of polymer retention, inaccessible pore volumes, a permeability reduction and polymer absorption are considered. Darcy's law is employed to model the behavoir of a fluid in porous media, and the upstream finite difference (volume) method is applied to discretize the mass conservation equations. Numerical methods, including discretization schemes, linear solver methods, nonlinearization methods and parallel techniques are introduced. Numerical experiments show that, on one hand, computed results match those from the commercial simulator, Schlumberger-Eclipse, which is widely applied by the petroleum industry, and, on the other hand, our simulator has excellent scalability, which is demonstrated by field applications with up to 27 million grid blocks using up to 2048 CPU cores.


title: Orthogonalization of fermion k-body operators and representability

abstract: The reduced k-particle density matrix of a density matrix on finite-dimensional, fermion Fock space can be defined as the image under the orthogonal projection in the Hilbert-Schmidt geometry onto the space of k-body observables. A proper understanding of this projection is therefore intimately related to the representability problem, a long-standing open problem in computational quantum chemistry. Given an orthonormal basis in the finite-dimensional one-particle Hilbert space, we explicitly construct an orthonormal basis of the space of Fock space operators which restricts to an orthonormal basis of the space of k-body operators for all k.


title: A Bi-Directional Big Bang / Crunch Universe within a Two-State-Vector Quantum Mechanics?

abstract: A two boundary quantum mechanics incorporating a big bang / big crunch universe is carefully considered. After a short motivation of the concept we address the central question how a proposed a-causal quantum universe can be consistent with what is known about macroscopia and how it might find experimental support.


title: Robust parameter determination in epidemic models with analytical descriptions of uncertainties

abstract: Compartmental equations are primary tools in disease spreading studies. Their predictions are accurate for large populations but disagree with empirical and simulated data for finite populations, where uncertainties become a relevant factor. Starting from the agent-based approach, we investigate the role of uncertainties and autocorrelation functions in SIS epidemic model, including their relationship with epidemiological variables. We find new differential equations that take uncertainties into account. The findings provide improved predictions to the SIS model and it can offer new insights for emerging diseases.


title: High Energy Physics and Cosmology at the Unification Frontier: Opportunities and Challenges in the coming years

abstract: We give here an overview of recent developments in high energy physics and cosmology and their interconnections that relate to unification, and discuss prospects for the future. Thus there are currently three empirical data that point to supersymmetry as an underlying symmetry of particle physics: the unification of gauge couplings within supersymmetry, the fact that nature respects the supersymmetry prediction that the Higgs boson mass lie below 130 GeV, and vacuum stability up to the Planck scale with a Higgs boson mass at $\sim 125$ GeV while the standard model does not do that. Coupled with the fact that supersymmetry solves the big hierarchy problem related to the quadratic divergence to the Higgs boson mass square along with the fact that there is no alternative paradigm that allows us to extrapolate physics from the electroweak scale to the grand unification scale consistent with experiment, supersymmetry remains a compelling framework for new physics beyond the standard model. The large loop correction to the Higgs boson mass in supersymmetry to lift the tree mass to the experimentally observable value, indicates a larger value of the scale of weak scale supersymmetry, making the observation of sparticles more challenging but still within reach at the LHC for the lightest ones. Recent analyses show that a high energy LHC (HE-LHC) operating at 27 TeV running at its optimal luminosity of $2.5 \times 10^{35}$ cm$^{-2}$s$^{-1}$ can reduce the discovery period by several years relative to HL-LHC and significantly extend the reach in parameter space of models. In the coming years several experiments related to neutrino physics, searches for supersymmetry, on dark matter and dark energy will have direct impact on the unification frontier.


title: Suppression of superfluid stiffness near Lifshitz-point instability to finite momentum superconductivity

abstract: We derive the effective Ginzburg-Landau theory for finite momentum (FFLO/PDW) superconductivity without spin population imbalance from a model with local attraction and repulsive pair-hopping. We find that the GL free energy must include up to sixth order derivatives of the order parameter, providing a unified description of the interdependency of zero and finite momentum superconductivity. For weak pair-hopping the phase diagram contains a line of Lifshitz points where vanishing superfluid stiffness induces a continuous change to a long wavelength Fulde-Ferrell (FF) state. For larger pair-hopping there is a bicritical region where the pair-momentum changes discontinuously. Here the FF type state is near degenerate with the Larkin-Ovchinnikov (LO) or Pair-Density-wave (PDW) type state. At the intersection of these two regimes there is a "Super-Lifshitz" point with extra soft fluctuations. The instability to finite momentum superconductivity occurs for arbitrarily weak pair-hopping for sufficiently large attraction suggesting that even a small repulsive pair-hopping may be significant in a microscopic model of strongly correlated superconductivity. Several generic features of the model may have bearing on the cuprate superconductors, including the suppression of superfluid stiffness in proximity to a Lifshitz point as well as the existence of subleading FFLO order (or vice versa) in the bicritical regime.


title: Doped and Codoped Silicon Nanocrystals: the Role of Surfaces and Interfaces

abstract: Si nanocrystals have been extensively studied because of their novel properties and their potential applications in electronic, optoelectronic, photovoltaic, thermoelectric and biological devices. These new properties are achieved through the combination of the quantum confinement of carriers and the strong influence of surface chemistry. As in the case of bulk Si the tuning of the electronic, optical and transport properties is related to the possibility of doping, in a controlled way, the nanocrystals. This is a big challenge since several studies have revealed that doping in Si nanocrystals differs from the one of the bulk. Theory and experiments have underlined that doping and codoping are influenced by a large number of parameters such as size, shape, passivation and chemical environment of the silicon nanocrystals. However, the connection between these parameters and dopant localization as well as the occurrence of self-purification effects are still not clear. In this review we summarize the latest progress in this fascinating research field considering free-standing and matrix-embedded Si nanocrystals both from the theoretical and experimental point of view, with special attention given to the results obtained by ab-initio calculations and to size-, surface- and interface-induced effects.


title: The Origin of Heavy Element Content Trend in Giant Planets via Core Accretion

abstract: We explore the origin of the trend of heavy elements in observed massive exoplanets. Coupling of better measurements of mass ($M_p$) and radius of exoplanets with planet structure models enables estimating the total heavy element mass ($M_Z$) in these planets. The corresponding relation is characterized by a power-law profile, $M_Z \propto M_p^{3/5}$. We develop a simplified, but physically motivated analysis to investigate how the power-law profile can be produced under the current picture of planet formation. Making use of the existing semi-analytical formulae of accretion rates of pebbles and planetesimals, our analysis shows that the relation can be reproduced well if it traces the final stage of planet formation. In the stage, planets accrete solids from gapped planetesimal disks and gas accretion is limited by disk evolution. We also find that dust accretion accompanying with gas accretion does not contribute to $M_Z$ for planets with $M_p < 10^3 M_{\oplus}$. Our findings are broadly consistent with that of previous studies, yet we explicitly demonstrate how planetesimal dynamics is crucial for better understanding the relation. While our approach is simple, we can also reproduce the trend of a correlation between planet metallicity and $M_p$ that is obtained by detailed population synthesis calculations, when the same assumption is adopted. Our analysis suggests that pebble accretion would not play a direct role at the final stage of planet formation, whereas radial drift of pebbles might be important indirectly for metal enrichment of planets. Detailed numerical simulations and more observational data are required for confirming our analysis.


title: Generative Adversarial Privacy

abstract: We present a data-driven framework called generative adversarial privacy (GAP). Inspired by recent advancements in generative adversarial networks (GANs), GAP allows the data holder to learn the privatization mechanism directly from the data. Under GAP, finding the optimal privacy mechanism is formulated as a constrained minimax game between a privatizer and an adversary. We show that for appropriately chosen adversarial loss functions, GAP provides privacy guarantees against strong information-theoretic adversaries. We also evaluate GAP's performance on the GENKI face database.


title: How Do Classifiers Induce Agents To Invest Effort Strategically?

abstract: Algorithms are often used to produce decision-making rules that classify or evaluate individuals. When these individuals have incentives to be classified a certain way, they may behave strategically to influence their outcomes. We develop a model for how strategic agents can invest effort in order to change the outcomes they receive, and we give a tight characterization of when such agents can be incentivized to invest specified forms of effort into improving their outcomes as opposed to "gaming" the classifier. We show that whenever any "reasonable" mechanism can do so, a simple linear mechanism suffices.


title: TabulaROSA: Tabular Operating System Architecture for Massively Parallel Heterogeneous Compute Engines

abstract: The rise in computing hardware choices is driving a reevaluation of operating systems. The traditional role of an operating system controlling the execution of its own hardware is evolving toward a model whereby the controlling processor is distinct from the compute engines that are performing most of the computations. In this context, an operating system can be viewed as software that brokers and tracks the resources of the compute engines and is akin to a database management system. To explore the idea of using a database in an operating system role, this work defines key operating system functions in terms of rigorous mathematical semantics (associative array algebra) that are directly translatable into database operations. These operations possess a number of mathematical properties that are ideal for parallel operating systems by guaranteeing correctness over a wide range of parallel operations. The resulting operating system equations provide a mathematical specification for a Tabular Operating System Architecture (TabulaROSA) that can be implemented on any platform. Simulations of forking in TabularROSA are performed using an associative array implementation and compared to Linux on a 32,000+ core supercomputer. Using over 262,000 forkers managing over 68,000,000,000 processes, the simulations show that TabulaROSA has the potential to perform operating system functions on a massively parallel scale. The TabulaROSA simulations show 20x higher performance as compared to Linux while managing 2000x more processes in fully searchable tables.


title: No Strong Geometric Beaming in the Ultraluminous Neutron Star Binary NGC 300 ULX-1 (SN 2010da) from Swift and Gemini

abstract: We have obtained near-simultaneous Swift/XRT imaging and Gemini GMOS spectroscopy for the ultraluminous X-ray source (ULX) NGC~300 ULX-1 (formerly designated SN~2010da). The observed X-ray emission is consistent with an inhomogeneous wind that partially obscures a central, bright inner accretion disk. We simultaneously fit eleven 0.3-10 keV spectra obtained over a $\sim$1 year time period (2016 April to 2017 July) using the same partial covering model, and find that although the covering fraction varies significantly (from 78% to consistent with 0%), the unabsorbed luminosity remains essentially constant across all observations ($2-6\times10^{39}$ erg s$^{-1}$). A relatively high 0.3-10 keV fractional variability amplitude ($F_{\rm var}$) of $\sim$30% is observed in all eleven observations. Optical spectra from Gemini exhibit numerous emission lines (e.g., H$\alpha$, H$\beta$, He II $\lambda$4686) which suggest that the neutron star primary is photoionizing material in the immediate vicinity of the binary. We compare the He II $\lambda$4686 line luminosity ($\sim7-9\times10^{35}$ erg s$^{-1}$) to the contemporaneous soft X-ray emission and find the X-ray emission is broadly consistent with the observed He II line luminosity. The combination of our X-ray observations and optical spectroscopy suggest that geometric beaming effects in the ULX-1 system are minimal, making ULX-1 one of only a few bona fide ULXs to be powered by accretion onto a neutron star.


title: Data-constrained Model for Coronal Mass Ejections Using Graduated Cylindrical Shell Method

abstract: Coronal Mass Ejections (CMEs) are major drivers of extreme space weather conditions, this being a matter of serious concern for our modern technologically-dependent society. Development of numerical approaches that would simulate CME generation and propagation through the interplanetary space is an important step towards our capability to predict CME arrival times at Earth and their geo-effectiveness. In this paper, we utilize a data-constrained Gibson--Low (GL) flux rope model to generate CMEs. We derive the geometry of the initial GL flux rope using the Graduated Cylindrical Shell (GCS) method. This method uses multiple viewpoints from STEREO A & B Cor1/Cor2, and SOHO/LASCO C2/C3 coronagraphs to determine the size and orientation of a CME flux rope as it starts to erupt from the Sun. A flux rope generated in this way is inserted into a quasi-steady global magnetohydrodynamics (MHD) background solar wind flow driven by SDO/HMI line-of-sight magnetogram data, and erupts immediately. Numerical results obtained with the Multi-Scale Fluid-Kinetic Simulation Suite (MS-FLUKSS) code are compared with STEREO and SOHO/LASCO coronagraph observations in particular in terms of the CME speed, acceleration, and magnetic field structure.


title: A matching-based heuristic algorithm for school bus routing problems

abstract: School bus planning problem (SBPP) has drawn much research attention due to the huge costs of school transportation. In the literature, the SBPP is usually decomposed into the routing and scheduling subproblems due to its complexity. Because of the nature of the decomposition, even if all the subproblems are solved to optimality, the final solution may not be as good as the solution from the integrated model. In this paper, we present a new approach that incorporates the scheduling information (namely the trip compatibility) into the routing stage such that the interrelationship between the subproblems is still considered even in the decomposed problems. A novel two-step heuristic adopting the trip compatibility idea is presented to solve the school bus routing problem. The first step finds an initial solution using an iterative minimum cost matching-based insertion heuristic. Then, the initial trips are improved using a Simulated Annealing and Tabu Search hybrid method. Experiments were conducted on randomly generated problems and benchmark problems in the literature. The result shows that our two-step heuristic improves existing solutions up to 25% on the benchmark problems.


title: Comparing Redundant and Sky Model Based Interferometric Calibration: A First Look with Phase II of the MWA

abstract: Interferometric arrays seeking to measure the 21 cm signal from the Epoch of Reionization must contend with overwhelmingly bright emission from foreground sources. Accurate recovery of the 21 cm signal will require precise calibration of the array, and several new avenues for calibration have been pursued in recent years, including methods using redundancy in the antenna configuration. The newly upgraded Phase II of Murchison Widefield Array (MWA) is the first interferometer that has large numbers of redundant baselines while retaining good instantaneous UV-coverage. This array therefore provides a unique opportunity to compare redundant calibration with sky-model based algorithms. In this paper, we present the first results from comparing both calibration approaches with MWA Phase II observations. For redundant calibration, we use the package OMNICAL, and produce sky-based calibration solutions with the analysis package Fast Holographic Deconvolution (FHD). There are three principal results. (1) We report the success of OMNICAL on observations of ORBComm satellites, showing substantial agreement between redundant visibility measurements after calibration. (2) We directly compare OMNICAL calibration solutions with those from FHD, and demonstrate these two different calibration schemes give extremely similar results. (3) We explore improved calibration by combining OMNICAL and FHD. We evaluate these combined methods using power spectrum techniques developed for EoR analysis and find evidence for marginal improvements mitigating artifacts in the power spectrum. These results are likely limited by signal-to-noise in the six hours of data used, but suggest future directions for combining these two calibration schemes.


title: IV estimation of causal hazard ratio

abstract: Cox's proportional hazards model is one of the most popular statistical models to evaluate associations of exposure with a censored failure time outcome. When confounding factors are not fully observed, the exposure hazard ratio estimated using a Cox model is subject to unmeasured confounding bias. To address this, we propose a novel approach for the identification and estimation of the causal hazard ratio in the presence of unmeasured confounding factors. Our approach is based on a binary instrumental variable, and an additional no-interaction assumption in a first stage regression of the treatment on the IV and unmeasured confounders. We propose, to the best of our knowledge, the first consistent estimator of the (population) causal hazard ratio within an instrumental variable framework. A version of our estimator admits a closed-form representation. We derive the asymptotic distribution of our estimator, and provide a consistent estimator for its asymptotic variance. Our approach is illustrated via simulation studies and a data application.


title: Gamma Spaces and Information

abstract: We investigate the role of Segal's Gamma-spaces in the context of classical and quantum information, based on categories of finite probabilities with stochastic maps and density matrices with quantum channels. The information loss functional extends to the setting of probabilistic Gamma-spaces considered here. The Segal construction of connective spectra from Gamma-spaces can be used in this setting to obtain spectra associated to certain categories of gapped systems.