title: The Pristine Dwarf-Galaxy survey - II. In-depth observational study of the faint Milky Way satellite Sagittarius II

abstract: We present an extensive study of the Sagittarius II (Sgr II) stellar system using MegaCam g and i photometry, narrow-band, metallicity-sensitive Calcium H&K doublet photometry, augmented with Keck II/DEIMOS multi-object spectroscopy. We are able to derive and refine the Sgr II structural and stellar properties: the colour-magnitude diagram implies Sgr II is old (12.0 +- 0.5 Gyr) and metal-poor. The CaHK photometry confirms the metal-poor nature of the satellite ([Fe/H]CaHK = -2.32 +- 0.04 dex) and suggests that Sgr II hosts more than one single stellar population ({\sigma} CaHK = 0.11+0.05-0.03 dex). From the deep spectroscopic data, the velocity dispersion of the system is found to be {\sigma}vr = 2.7+1.3-1.0 km s-1 after excluding two potential binary stars. Using the Ca infrared triplet measured from our highest signal-to-noise spectra, we are able to confirm the metallicity and dispersion inferred from the Pristine photometric metallicities: ([Fe/H] = -2.23 +- 0.05 dex, {\sigma}spectro = 0.10+0.06-0.04 dex). Sgr II's metallicity and absolute magnitude (MV = -5.7 +- 0.1 mag) place the system on the luminosity-metallicity relation of the Milky Way dwarf galaxies despite its small size. The low, but resolved metallicity and velocity dispersions paint the picture of a slightly dark matter-dominated satellite. Furthermore, using the Gaia Data Release 2, we constrain the orbit of the satellite and find an apocenter of 118.4+28.4-23.7 kpc and a pericenter of 54.8+3.3-6.1 kpc. The orbit of Sgr II is consistent with the trailing arm of the Sgr stream and indicates that it is possibly a satellite of the Sgr dSph that was tidally stripped from the dwarf's influence.


title: English Translation of Cournot's "Exposition de la th\'eorie des chances et des probabilit\'es" / "Exposition of the Theory of Chances and Probabilities"

abstract: In spite of the wide range of his book, Cournot did not know some essential discoveries in natural sciences (William Herschel, Daniel Bernoulli, Humboldt) and his deliberations about measurement were almost useless. But he introduced the median and a definition of probability that covered both the discrete and continuous cases; studied mixtures of densities; qualitatively explained randomness anticipating one of Poincar\`e's interpretations (small causes leading to considerable effects) and introducing an embryo of irregularity as a condition of randomness; introduced a pattern of stratified sampling; attempted to study dependence between the decisions of jurors; introduced (after J. F. Fries) philosophical, i.e. qualitative probabilities. Cournot stated that statistics ought to have its own theory, rules and principles and should penetrate into the essence of things (which was not yet a received opinion). Many of his pronouncements were thoughtful and echoed the statistical thought of the time. In 1849 his book was translated into German.


title: English Translation of Poisson's "Recherches sur la probabilit\'e des jugements en mati\`ere criminelle et en mati\`ere civile" / "Researches into the Probabilities of Judgements in Criminal and Civil Cases"

abstract: In spite of its title, the book mostly treats probability theory: the law of large numbers (regarded as a principle); formal definition of a random variable and law of distribution; the misnamed Cauchy distribution; functions now named after Dirac; difference between objective and subjective probabilities; approximation of the Bernoulli theorem for rare events (but not yet the Poisson distribution); investigation of the significance of empirical discrepancies. Following Laplace, he actually attributed probability to applied mathematics. Poisson calculated the optimal number of jurors and of their majority for conviction/exoneration of the accused and in general studied the conditions enabling to minimize miscarriage of justice. He replaced guilty and innocent by convictable or not.


title: Humor in Word Embeddings: Cockamamie Gobbledegook for Nincompoops

abstract: While humor is often thought to be beyond the reach of Natural Language Processing, we show that several aspects of single-word humor correlate with simple linear directions in Word Embeddings. In particular: (a) the word vectors capture multiple aspects discussed in humor theories from various disciplines; (b) each individual's sense of humor can be represented by a vector, which can predict differences in people's senses of humor on new, unrated, words; and (c) upon clustering humor ratings of multiple demographic groups, different humor preferences emerge across the different groups. Humor ratings are taken from the work of Engelthaler and Hills (2017) as well as from an original crowdsourcing study of 120,000 words. Our dataset further includes annotations for the theoretically-motivated humor features we identify.


title: A Survey of the Individual-Based Model applied in Biomedical and Epidemiology

abstract: Individual-based model (IBM) has been used to simulate and to design control strategies for dynamic systems that are subject to stochasticity and heterogeneity, such as infectious diseases. In the IBM, an individual is represented by a set of specific characteristics that may change dynamically over time. This feature allows a more realistic analysis of the spread of an epidemic. This paper presents a literature survey of IBM applied to biomedical and epidemiology research. The main goal is to present existing techniques, advantages and future perspectives in the development of the model. We evaluated 89 articles, which mostly analyze interventions aimed at endemic infections. In addition to the review, an overview of IBM is presented as an alternative to complement or replace compartmental models, such as the SIR (Susceptible-Infected-Recovered) model. Numerical simulations also illustrate the capabilities of IBM, as well as some limitations regarding the effects of discretization. We show that similar side-effects of discretization scheme for compartmental models may also occur in IBM, which requires careful attention.


title: Fundamental physics constraints from testing the stability of the fine-structure constant with the ELTs

abstract: The increased collecting area of the ELTs will bring fainter high-z targets within the reach of high-resolution ultra-stable spectrographs, thus enabling a new generation of precision consistency tests, including tests of the stability of nature's fundamental couplings. For example, the stability of the fine-structure constant can be tested by looking at metal absorption lines produced by the intervening clouds along the line of sight of distant quasars. In this contribution, we discuss the performance that can be expected from the ELTs in testing the stability of the fine-structure constant, based on the early ESPRESSO observations, and some comparative forecasts of the impact of these measurements for representative models of fundamental physics and cosmology.


title: The mean H$\alpha$ EW and Lyman-continuum photon production efficiency for faint $z\approx4-5$ galaxies

abstract: We present the first measurements of the Lyman-continuum photon production efficiency $\xi_{\textrm{ion,0}}$ at $z\sim4$-5 for galaxies fainter than 0.2 $L^*$ ($-$19 mag). $\xi_{\textrm{ion,0}}$ quantifies the production rate of ionizing photons with respect to the UV luminosity density assuming a fiducial escape fraction of zero. Extending previous measurements of $\xi_{\textrm{ion,0}}$ to the faint population is important, as ultra-faint galaxies are expected to contribute the bulk of the ionizing emissivity. We probe $\xi_{\textrm{ion,0}}$ to such faint magnitudes by taking advantage of 200-hour depth Spitzer/IRAC observations from the GREATS program and $\approx$300 3<$z$<6 galaxies with spectroscopic redshifts from the MUSE GTO Deep + Wide programs. Stacked IRAC [3.6]$-$[4.5] colors are derived and used to infer the H$\alpha$ rest-frame equivalent widths, which range from 403\r{A} to 2818\r{A}. The derived $\xi_{\textrm{ion,0}}$ is $\log_{10}(\xi_{\textrm{ion,0}} / \textrm{Hz erg}^{-1}) = 25.36 \pm 0.08$ over $-$20.5 < M$_{\textrm{UV}}$ < $-$17.5, similar to those derived for brighter galaxy samples at the same redshift and therefore suggesting that $\xi_{\textrm{ion}}$ shows no strong dependence on $M_{UV}$. The $\xi_{\textrm{ion,0}}$ values found in our sample imply that the Lyman-continuum escape fraction for $M_{\textrm{UV}} \approx -19$ star-forming galaxies cannot exceed $\approx$8-20\% in the reionization era.


title: The Hubble Space Telescope UV Legacy Survey of Galactic Globular Clusters. XVIII. Proper-motion kinematics of multiple stellar populations in the core regions of NGC 6352

abstract: We present the analysis of the radial distributions and kinematic properties of the multiple stellar populations (mPOPs) hosted in the globular cluster (GC) NGC 6352 as part of the Hubble Space Telescope (HST) "UV Legacy Survey of Galactic Globular Clusters" program. NGC 6352 is one of the few GCs for which the mPOP tagging in appropriate color-magnitude diagrams is clear in all evolutionary sequences. We computed high-precision stellar proper motions for the stars from the cluster's core out to 75 arcsec ($\sim$1.5 core radii, or $\sim$0.6 half-light radii). We find that, in the region explored, first- and second-generation stars share the same radial distribution and kinematic properties. Velocity dispersions, anisotropy radial profiles, differential rotation, and level of energy equipartition, all suggest that NGC 6352 is probably in an advanced evolutionary stage, and any possible difference in the structural and kinematic properties of its mPOPs have been erased by dynamical processes in the core of the cluster. We also provide an estimate of the mass of blue stragglers and of main-sequence binaries through kinematics alone. In general, in order to build a complete dynamical picture of this and other GCs, it will be essential to extend the analyses presented in this paper to the GCs' outer regions where some memories of the initial differences in the mPOP properties, and those imprinted by dynamical processes, might still be present.


title: Constraining Rapidly Oscillating Scalar Dark Matter Using Dynamic Decoupling

abstract: We propose and experimentally demonstrate a method for detection of light scalar Dark Matter (DM), through probing temporal oscillations of fundamental constants in an atomic optical transition. Utilizing the quantum information notion of Dynamic Decoupling (DD) in a table-top setting, we are able to obtain model-independent bounds on variations of $\alpha$ and $m_e$ at frequencies up to the MHz scale. We interpret our results to constrain the parameter space of light scalar DM field models. We consider the generic case, where the couplings of the DM field to the photon and to the electron are independent, as well as the case of a relaxion DM model, including the scenario of a DM boson star centered around Earth. Given the particular nature of DD, allowing to directly observe the oscillatory behaviour of coherent DM, and considering future experimental improvements, we conclude that our proposed method could be complimentary to, and possibly competitive with, gravitational probes of light scalar DM.


title: K2-288Bb: A Small Temperate Planet in a Low-mass Binary System Discovered by Citizen Scientists

abstract: Observations from the Kepler and K2 missions have provided the astronomical community with unprecedented amounts of data to search for transiting exoplanets and other astrophysical phenomena. Here, we present K2-288, a low-mass binary system (M2.0 +/- 1.0; M3.0 +/- 1.0) hosting a small (Rp = 1.9 REarth), temperate (Teq = 226 K) planet observed in K2 Campaign 4. The candidate was first identified by citizen scientists using Exoplanet Explorers hosted on the Zooniverse platform. Follow-up observations and detailed analyses validate the planet and indicate that it likely orbits the secondary star on a 31.39-day period. This orbit places K2-288Bb in or near the habitable zone of its low-mass host star. K2-288Bb resides in a system with a unique architecture, as it orbits at >0.1 au from one component in a moderate separation binary (aproj approximately 55 au), and further follow-up may provide insight into its formation and evolution. Additionally, its estimated size straddles the observed gap in the planet radius distribution. Planets of this size occur less frequently and may be in a transient phase of radius evolution. K2-288 is the third transiting planet system identified by the Exoplanet Explorers program and its discovery exemplifies the value of citizen science in the era of Kepler, K2, and the Transiting Exoplanet Survey Satellite.


title: Interacting multi-channel topological boundary modes in a quantum Hall valley system

abstract: Symmetry and topology play key roles in the identification of phases of matter and their properties. Both concepts are central to understanding quantum Hall ferromagnets (QHFMs), two-dimensional electronic phases with spontaneously broken spin or pseudospin symmetry whose wavefunctions also have topological properties. Domain walls between distinct broken symmetry QHFM phases are predicted to host gapless one-dimensional (1D) modes that emerge due to a topological change of the underlying electronic wavefunctions at such interfaces. Although a variety of QHFMs have been identified in different materials, probing interacting electronic modes at these domain walls has not yet been accomplished. Here we use a scanning tunneling microscope (STM) to directly visualize the spontaneous formation of boundary modes, within a sign-changing topological gap, at domain walls between different valley-polarized quantum Hall phases on the surface of bismuth. By changing the valley occupation and the corresponding number of modes at the domain wall, we can realize different regimes where the valley-polarized channels are either metallic or develop a spectroscopic gap. This behavior is a consequence of Coulomb interactions constrained by the symmetry-breaking valley flavor, which determines whether electrons in the topological modes can backscatter, making these channels a unique class of interacting Luttinger liquids.


title: Stability of the vacuum as constraint on $U$(1) extensions of the standard model

abstract: In the standard model the running quartic coupling becomes negative during its renormalization group flow, which destabilizes the vacuum. We consider U(1) extensions of the standard model, with an extra complex scalar field and a Majorana-type neutrino Yukawa coupling. These additional couplings affect the renormalization group flow of the quartic couplings. We compute the beta-functions of the extended model at one-loop order in perturbation theory and study how the parameter space of the new scalar couplings can be constrained by the requirement of stable vacuum and perturbativity up to the Planck scale.


title: Conditions for Reionizing the Universe with A Low Galaxy Ionizing Photon Escape Fraction

abstract: We explore scenarios for reionizing the intergalactic medium with low galaxy ionizing photon escape fractions. We combine simulation-based halo-mass dependent escape fractions with an extrapolation of the observed galaxy rest-ultraviolet luminosity functions to solve for the reionization history from z=20 to z=4. We explore the posterior distributions for key unknown quantities, including the limiting halo mass for star-formation, the ionizing photon production efficiency, and a potential contribution from active galactic nuclei (AGN). We marginalize over the allowable parameter space using a Markov Chain Monte Carlo method, finding a solution which satisfies the most model-independent constraints on reionization. Our fiducial model can match observational constraints with an average escape fraction of <5% throughout the bulk of the epoch of reionization if: i) galaxies form stars down to the atomic cooling limit before reionization and a photosuppression mass of log(M_h/Msol)~9 during/after reionization (-13<M_UV,lim<-11); ii) galaxies become more efficient producers of ionizing photons at higher redshifts and fainter magnitudes, and iii) there is a significant, but sub-dominant, contribution by AGN at z < 7. In this model the faintest galaxies (M_UV>-15) dominate the ionizing emissivity, leading to an earlier start to reionization and a smoother evolution of the ionized volume filling fraction than models which assume a single escape fraction at all redshifts and luminosities. The ionizing emissivity from this model is consistent with observations at z=4-5 (and below, when extrapolated), in contrast to some models which assume a single escape fraction. Our predicted ionized volume filling fraction at z=7 of Q_HII=78% (+\- 8%) is in ~1-2 sigma tension with observations of Lya emitters at z~7 and the damping wing analyses of the two known z>7 quasars, which prefer Q_HII,z=7~40-50%.


title: 'The Brick' is not a brick: A comprehensive study of the structure and dynamics of the Central Molecular Zone cloud G0.253+0.016

abstract: In this paper we provide a comprehensive description of the internal dynamics of G0.253+0.016 (a.k.a. 'the Brick'); one of the most massive and dense molecular clouds in the Galaxy to lack signatures of widespread star formation. As a potential host to a future generation of high-mass stars, understanding largely quiescent molecular clouds like G0.253+0.016 is of critical importance. In this paper, we reanalyse Atacama Large Millimeter Array cycle 0 HNCO $J=4(0,4)-3(0,3)$ data at 3 mm, using two new pieces of software which we make available to the community. First, scousepy, a Python implementation of the spectral line fitting algorithm scouse. Secondly, acorns (Agglomerative Clustering for ORganising Nested Structures), a hierarchical n-dimensional clustering algorithm designed for use with discrete spectroscopic data. Together, these tools provide an unbiased measurement of the line of sight velocity dispersion in this cloud, $\sigma_{v_{los}, {\rm 1D}}=4.4\pm2.1$ kms$^{-1}$, which is somewhat larger than predicted by velocity dispersion-size relations for the Central Molecular Zone (CMZ). The dispersion of centroid velocities in the plane of the sky are comparable, yielding $\sigma_{v_{los}, {\rm 1D}}/\sigma_{v_{pos}, {\rm 1D}}\sim1.2\pm0.3$. This isotropy may indicate that the line-of-sight extent of the cloud is approximately equivalent to that in the plane of the sky. Combining our kinematic decomposition with radiative transfer modelling we conclude that G0.253+0.016 is not a single, coherent, and centrally-condensed molecular cloud; 'the Brick' is not a \emph{brick}. Instead, G0.253+0.016 is a dynamically complex and hierarchically-structured molecular cloud whose morphology is consistent with the influence of the orbital dynamics and shear in the CMZ.


title: Line-Graph Lattices: Euclidean and Non-Euclidean Flat Bands, and Implementations in Circuit Quantum Electrodynamics

abstract: Materials science and the study of the electronic properties of solids are a major field of interest in both physics and engineering. The starting point for all such calculations is single-electron, or non-interacting, band structure calculations, and in the limit of strong on-site confinement this can be reduced to graph-like tight-binding models. In this context, both mathematicians and physicists have developed largely independent methods for solving these models. In this paper we will combine and present results from both fields. In particular, we will discuss a class of lattices which can be realized as line graphs of other lattices, both in Euclidean and hyperbolic space. These lattices display highly unusual features including flat bands and localized eigenstates of compact support. We will use the methods of both fields to show how these properties arise and systems for classifying the phenomenology of these lattices, as well as criteria for maximizing the gaps. Furthermore, we will present a particular hardware implementation using superconducting coplanar waveguide resonators that can realize a wide variety of these lattices in both non-interacting and interacting form.

categories: quant-ph cond-mat.other math-ph math.CO math.MP


title: Speeding past planets? Asteroids radiatively propelled by giant branch Yarkovsky effects

abstract: Understanding the fate of planetary systems through white dwarfs which accrete debris crucially relies on tracing the orbital and physical properties of exo-asteroids during the giant branch phase of stellar evolution. Giant branch luminosities exceed the Sun's by over three orders of magnitude, leading to significantly enhanced Yarkovsky and YORP effects on minor planets. Here, we place bounds on Yarkovsky-induced differential migration between asteroids and planets during giant branch mass loss by modelling one exo-Neptune with inner and outer exo-Kuiper belts. In our bounding models, the asteroids move too quickly past the planet to be diverted from their eventual fate, which can range from: (i) populating the outer regions of systems out to 10^4-10^5 au, (ii) being engulfed within the host star, or (iii) experiencing Yarkovsky-induced orbital inclination flipping without any Yarkovsky-induced semimajor axis drift. In these violent limiting cases, temporary resonant trapping of asteroids with radii of under about 10 km by the planet is insignificant, and capture within the planet's Hill sphere requires fine-tuned dissipation. The wide variety of outcomes presented here demonstrates the need to employ sophisticated structure and radiative exo-asteroid models in future studies. Determining where metal-polluting asteroids reside around a white dwarf depends on understanding extreme Yarkovsky physics.


title: An in-depth analysis of $b\to c(s)$ semileptonic observables with possible $\mu - \tau$ mixing

abstract: In a couple of recent publications ( arXiv:1706.08437 and arXiv:1712.01593 ), the authors attempted to achieve simultaneous explanation of the persistent flavor anomalies in $b\to s$ and $b\to c$ semileptonic decays with a minimal scheme by using only three unknown new parameters. The analysis was obtained with a handful of precise observables. Motivated by their proposal, in this paper we reanalyze the models proposed in the aforementioned papers with a total of 170 observables from those channels including newly available measurements, correlated theoretical results, and constraints. We validate our results by searching for the most influential points and outliers. By analyzing the parameter spaces and their relationship with the constraints, we gain new insight and statistical significance in those models. We also provide a new and precise calculation of $R(J/\Psi)$, obtained during the analysis.


title: Localization of Binary Black-Hole Mergers with Known Inclination

abstract: The localization of stellar-mass binary black hole mergers using gravitational waves is critical in understanding the properties of the binaries' host galaxies, observing possible electromagnetic emission from the mergers, or using them as a cosmological distance ladder. The precision of this localization can be substantially increased with prior astrophysical information about the binary system. In particular, constraining the inclination of the binary can reduce the distance uncertainty of the source. Here we present the first realistic set of localizations for binary black hole mergers, including different prior constraints on the binaries' inclinations. We find that prior information on the inclination can reduce the localization volume by a factor of 3. We discuss two astrophysical scenarios of interest: (i) follow-up searches for beamed electromagnetic/neutrino counterparts and (ii) mergers in the accretion disks of active galactic nuclei.


title: Triplet excitations in the frustrated spin ladder Li$_2$Cu$_2$O(SO$_4$)$_2$

abstract: Magnetic excitations of the recently discovered frustrated spin-1/2 two-leg ladder system Li$_2$Cu$_2$O(SO$_4$)$_2$ are investigated using inelastic neutron scattering, magnetic susceptibility and infrared absorption measurements. Despite the presence of a magnetic dimerization concomitant with the tetragonal-to-triclinic structural distortion occurring below 125 K, neutron scattering experiments reveal the presence of dispersive triplet excitations above a spin gap of $\Delta = 10.6$ meV at 1.5 K, a value consistent with the estimates extracted from magnetic susceptibility. The likely detection of these spin excitations in infrared spectroscopy is explained by invoking a dynamic Dzyaloshinskii-Moriya mechanism in which light is coupled to the dimer singlet-to-triplet transition through an optical phonon. These results are qualitatively explained by exact diagonalization and higher-order perturbation calculations carried out on the basis of the dimerized spin Hamiltonian derived from first-principles.


title: Classical Dimers on Penrose Tilings

abstract: We study the classical dimer model on rhombic Penrose tilings, whose edges and vertices may be identified with those of a bipartite graph. We find that Penrose tilings do not admit perfect matchings (defect-free dimer coverings). Instead, their maximum matchings have a monomer density of $81-50\varphi\approx 0.098$ in the thermodynamic limit, with $\varphi=\left(1+\sqrt{5}\right)/2$ the golden ratio. Maximum matchings divide the tiling into a fractal of nested closed regions bounded by loops that cannot be crossed by monomers. These loops connect second-nearest neighbour even-valence vertices, each of which lies on such a loop. Assigning a charge to each monomer with a sign fixed by its bipartite sublattice, we find that each bounded region has an excess of one charge, and a corresponding set of monomers, with adjacent regions having opposite net charge. The infinite tiling is charge neutral. We devise a simple algorithm for generating maximum matchings, and demonstrate that maximum matchings form a connected manifold under local monomer-dimer rearrangements. We show that dart-kite Penrose tilings feature an imbalance of charge between bipartite sub-lattices, leading to a minimum monomer density of $\left(7-4\varphi\right)/5\approx 0.106$ all of one charge.

categories: cond-mat.str-el cond-mat.stat-mech math-ph math.MP


title: Does AGN Fraction Depend on Redshift or Luminosity? An Extinction-Free Test by 18-band Mid-infrared SED Fitting in the AKARI NEP Wide Field

abstract: Revealing what fraction of galaxies harbor AGN is central in understanding black hole accretion history of the Universe. However, optical and soft X-ray surveys miss the most highly obscured AGNs. Infrared (IR), instead, is more robust against absorption. Previous IR photometric surveys, however, only had 4 or 5 filters in mid-IR. Our AKARI North Ecliptic Pole (NEP) wide field sample has 18 filters in mid-IR (9 from AKARI, 4 from WISE, and 5 from Spitzer), for the first time, allowing a sophisticated mid-IR SED fitting diagnosis for a statistical number of sources (89178 over 5.4 deg$^2$). By using a SED fitting technique, we investigate the evolution of AGN fraction as a function of redshift and IR (8-1000 $\mu$m) luminosity in an extinction-free way. We found that the AGN fraction (F$_{\rm AGN}$) shows no sign of strong redshift evolution. Instead, F$_{\rm AGN}$ increases with increasing IR luminosity in all redshifts bins ($0<z<2$).


title: Infrared luminosity functions based on 18 mid-infrared bands: revealing cosmic star formation history with AKARI and Hyper Suprime-Cam

abstract: Much of the star formation is obscured by dust. For the complete understanding of the cosmic star formation history (CSFH), infrared (IR) census is indispensable. AKARI carried out deep mid-infrared observations using its continuous 9-band filters in the North Ecliptic Pole (NEP) field (5.4 deg$^2$). This took significant amount of satellite's lifetime, $\sim$10\% of the entire pointed observations. By combining archival Spitzer (5 bands) and WISE (4 bands) mid-IR photometry, we have, in total, 18 band mid-IR photometry, which is the most comprehensive photometric coverage in mid-IR for thousands of galaxies. However previously, we only had shallow optical imaging ($\sim$25.9ABmag) in a small area of 1.0 deg$^2$. As a result, there remained thousands of AKARI's infrared sources undetected in optical. Using the new Hyper Suprime-Cam on Subaru telescope, we obtained deep enough optical images of the entire AKARI NEP field in 5 broad bands ($g\sim$27.5mag). These provided photometric redshift, and thereby IR luminosity for the previously undetected faint AKARI IR sources. Combined with the accurate mid-IR luminosity measurement, we constructed mid-IR LFs, and thereby performed a census of dust-obscured CSFH in the entire AKARI NEP field. We have measured restframe 8$\mu$m, 12$\mu$m luminosity functions (LFs), and estimated total infrared LFs at 0.35$<$z$<$2.2. Our results are consistent with our previous work, but with much reduced statistical errors thanks to the large area coverage of the new data. We have possibly witnessed the turnover of CSFH at $z\sim$2.


title: A strain tunable single-layer MoS2 photodetector

abstract: Strain engineering, which aims to tune the bandgap of a semiconductor by the application of strain, has emerged as an interesting way to control the electrical and optical properties of two-dimensional (2D) materials. Apart from the changes in the intrinsic properties of 2D materials, the application of strain can be also used to modify the characteristics of devices based on them. In this work, we study flexible and transparent photodetectors based on single-layer MoS2 under the application of biaxial strain. We find that by controlling the level of strain, we can tune the photoresponsivity (by 2-3 orders of magnitude), the response time (from <80 ms to 1.5 s) and the spectral bandwidth (with a gauge factor of 135 meV/% or 58 nm/%) of the device.


title: Revisiting the buckling metrology method to determine the Young's modulus of 2D materials

abstract: Measuring the mechanical properties of two-dimensional materials is a formidable task. While regular electrical and optical probing techniques are suitable even for atomically thin materials, conventional mechanical tests cannot be directly applied. Therefore, new mechanical testing techniques need to be developed. Up to now, the most widespread approaches require micro-fabrication to create freely suspended membranes, rendering their implementation complex and costly. Here, we revisit a simple yet powerful technique to measure the mechanical properties of thin films. The buckling metrology method, that does not require the fabrication of freely suspended structures, is used to determine the Young's modulus of several transition metal dichalcogenides (MoS2, MoSe2, WS2 and WSe2) with thicknesses ranging from 3 to 10 layers. We critically compare the obtained values for the Young's modulus and their uncertainty, finding that this simple technique provides results, which are in good agreement with those reported using other highly sophisticated testing methods. By comparing the cost, complexity and time required for the different methods reported in the literature, the buckling metrology method presents certain advantages that makes it an interesting mechanical test tool for 2D materials.


title: SiamVGG: Visual Tracking using Deeper Siamese Networks

abstract: Recently, we have seen a rapid development of Deep Neural Network (DNN) based visual tracking solutions. Some trackers combine the DNN-based solutions with Discriminative Correlation Filters (DCF) to extract semantic features and successfully deliver the state-of-the-art tracking accuracy. However, these solutions are highly compute-intensive, which require long processing time, resulting unsecured real-time performance. To deliver both high accuracy and reliable real-time performance, we propose a novel tracker called SiamVGG\footnote{https://github.com/leeyeehoo/SiamVGG}. It combines a Convolutional Neural Network (CNN) backbone and a cross-correlation operator, and takes advantage of the features from exemplary images for more accurate object tracking. The architecture of SiamVGG is customized from VGG-16 with the parameters shared by both exemplary images and desired input video frames. We demonstrate the proposed SiamVGG on OTB-2013/50/100 and VOT 2015/2016/2017 datasets with the state-of-the-art accuracy while maintaining a decent real-time performance of 50 FPS running on a GTX 1080Ti. Our design can achieve 2% higher Expected Average Overlap (EAO) compared to the ECO and C-COT in VOT2017 Challenge.


title: NOvA muon energy scale systematic

abstract: The systematic uncertainty on the correspondence between muon range and energy is developed for the NOvA neutrino experiment. NOvA consists of two detectors, the Near Detector at Fermilab and the Far Detector in northern Minnesota. Total errors are developed for the Near Detector, with its Muon Catcher treated separately, the Far Detector, and all combinations of correlated and uncorrelated errors between these three detectors. The absolute errors for the Near Detector (1.0%), the Far Detector (0.9%), and the fully correlated error shared by them (0.9%) are strongly dominated by Geant4's treatment of the Bethe density effect. At the Near Detector, the next biggest uncertainty is from stray hits caused by neutron capture pile-up. Other contributions are marginally significant, with the biggest, in descending order, being due to external measurements of the mean excitation energies of elements, detector mass accounting, and modification of energy loss by chemical binding. For the Muon Catcher, the absolute error is expressed as an offset instead of a percentage: 21 MeV. The density effect (at higher energies) and neutron capture pile-up (at lower energies) are the strongly dominant errors. The relative error between the Near and Far Detectors is 0.4% and is strongly dominated by neutron capture pile-up at the Near Detector, with a subdominant contribution from detector mass accounting.


title: Automated pebble mosaic stylization of images

abstract: Digital mosaics have usually used regular tiles, simulating the historical "tessellated" mosaics. In this paper, we present a method for synthesizing pebble mosaics, a historical mosaic style in which the tiles are rounded pebbles. We address both the tiling problem, where pebbles are distributed over the image plane so as to approximate the input image content, and the problem of geometry, creating a smooth rounded shape for each pebble. We adapt SLIC, simple linear iterative clustering, to obtain elongated tiles conforming to image content, and smooth the resulting irregular shapes into shapes resembling pebble cross-sections. Then, we create an interior and exterior contour for each pebble and solve a Laplace equation over the region between them to obtain height-field geometry. The resulting pebble set approximates the input image while presenting full geometry that can be rendered and textured for a highly detailed representation of a pebble mosaic.


title: The distance to NGC1042 in the context of its proposed association with the dark matter-deficient galaxies NGC1052-DF2 and NGC1052-DF4

abstract: It has been suggested that the dark matter-deficient galaxies NGC1052-DF2 and NGC1052-DF4 might not be members of the NGC1052 group but in the foreground at $\sim 13$ Mpc, and satellites of the bright spiral galaxy NGC1042. We previously showed that the CMDs of the galaxies are inconsistent with this hypothesis, and derived distances of 19-20 Mpc from their surface brightness fluctuation signals. Here we note that NGC1042 is almost certainly a member of the NGC1052 group as well, based on its radial velocity, the HI distribution in the NGC1052/NGC1042 system, and the Tully-Fisher relation.


title: ML Health: Fitness Tracking for Production Models

abstract: Deployment of machine learning (ML) algorithms in production for extended periods of time has uncovered new challenges such as monitoring and management of real-time prediction quality of a model in the absence of labels. However, such tracking is imperative to prevent catastrophic business outcomes resulting from incorrect predictions. The scale of these deployments makes manual monitoring prohibitive, making automated techniques to track and raise alerts imperative. We present a framework, ML Health, for tracking potential drops in the predictive performance of ML models in the absence of labels. The framework employs diagnostic methods to generate alerts for further investigation. We develop one such method to monitor potential problems when production data patterns do not match training data distributions. We demonstrate that our method performs better than standard "distance metrics", such as RMSE, KL-Divergence, and Wasserstein at detecting issues with mismatched data sets. Finally, we present a working system that incorporates the ML Health approach to monitor and manage ML deployments within a realistic full production ML lifecycle.


title: Kiloparsec-scale Variations in the Star Formation Efficiency of Dense Gas: the Antennae Galaxies (NGC 4038/39)

abstract: We study the relationship between dense gas and star formation in the Antennae galaxies by comparing ALMA observations of dense gas tracers (HCN, HCO$^+$, and HNC $\mathrm{J}=1-0$) to the total infrared luminosity ($\mathrm{L_{TIR}}$) calculated using data from the \textit{Herschel} Space Observatory and the \textit{Spitzer} Space Telescope. We compare the luminosities of our SFR and gas tracers using aperture photometry and employing two methods for defining apertures. We taper the ALMA dataset to match the resolution of our $\mathrm{L_{TIR}}$ maps and present new detections of dense gas emission from complexes in the overlap and western arm regions. Using OVRO CO $\mathrm{J}=1-0$ data, we compare with the total molecular gas content, $\mathrm{M(H_2)_{tot}}$, and calculate star formation efficiencies and dense gas mass fractions for these different regions. We derive HCN, HCO$^+$ and HNC upper limits for apertures where emission was not significantly detected, as we expect emission from dense gas should be present in most star-forming regions. The Antennae extends the linear $\mathrm{L_{TIR}-L_{HCN}}$ relationship found in previous studies. The $\mathrm{L_{TIR}-L_{HCN}}$ ratio varies by up to a factor of $\sim$10 across different regions of the Antennae implying variations in the star formation efficiency of dense gas, with the nuclei, NGC 4038 and NGC 4039, showing the lowest SFE$_\mathrm{dense}$ (0.44 and 0.70 $\times10^{-8}$ yr$^{-1}$). The nuclei also exhibit the highest dense gas fractions ($\sim 9.1\%$ and $\sim7.9\%$).


title: A new HLLD Riemann solver with Boris correction for reducing Alfv\'en speed

abstract: A new Riemann solver is presented for the ideal magnetohydrodynamics (MHD) equations with the so-called Boris correction. The Boris correction is applied to reduce wave speeds, avoiding an extremely small timestep in MHD simulations. The proposed Riemann solver, Boris-HLLD, is based on the HLLD solver. As done by the original HLLD solver, (1) the Boris-HLLD solver has four intermediate states in the Riemann fan when left and right states are given, (2) it resolves the contact discontinuity, Alfv\'en waves, and fast waves, and (3) it satisfies all the jump conditions across shock waves and discontinuities except for slow shock waves. The results of a shock tube problem indicate that the scheme with the Boris-HLLD solver captures contact discontinuities sharply and it exhibits shock waves without any overshoot when using the minmod limiter. The stability tests show that the scheme is stable when $|u| \lesssim 0.5c$ for a low Alfv\'en speed ($V_A \lesssim c$), where $u$, $c$, and $V_A$ denote the gas velocity, speed of light, and Alfv\'en speed, respectively. For a high Alfv\'en speed ($V_A \gtrsim c$), where the plasma beta is relatively low in many cases, the stable region is large, $|u| \lesssim (0.6-1) c$. We discuss the effect of the Boris correction on physical quantities using several test problems. The Boris-HLLD scheme can be useful for problems with supersonic flows in which regions with a very low plasma beta appear in the computational domain.

categories: physics.comp-ph astro-ph.IM astro-ph.SR physics.plasm-ph physics.space-ph


title: Classical and quantum integrable sigma models. Ricci flow, "nice duality" and perturbed rational conformal field theories

abstract: We consider classical and quantum integrable sigma models and their relations with the solutions of renormalization group equations. We say that an integrable sigma model possesses the "nice" duality property if the dual quantum field theory has the weak coupling region. As an example, we consider the deformed $CP(n-1)$ sigma model with additional quantum degrees of freedom. We formulate the dual integrable field theory and use perturbed conformal field theory, perturbation theory, $S$-matrix, Bethe Ansatz and renormalization group methods to show that this field theory has the "nice" duality property. We consider also an alternative approach to the analysis of sigma models on the deformed symmetric spaces, based on the perturbed rational conformal field theories.


title: Cooperative Training of Fast Thinking Initializer and Slow Thinking Solver for Conditional Learning

abstract: This paper studies the problem of learning the conditional distribution of a high-dimensional output given an input, where the output and input may belong to two different domains, e.g., the output is a photo image and the input is a sketch image. We solve this problem by cooperative training of a fast thinking initializer and slow thinking solver. The initializer generates the output directly by a non-linear transformation of the input as well as a noise vector that accounts for latent variability in the output. The slow thinking solver learns an objective function in the form of a conditional energy function, so that the output can be generated by optimizing the objective function, or more rigorously by sampling from the conditional energy-based model. We propose to learn the two models jointly, where the fast thinking initializer serves to initialize the sampling of the slow thinking solver, and the solver refines the initial output by an iterative algorithm. The solver learns from the difference between the refined output and the observed output, while the initializer learns from how the solver refines its initial output. We demonstrate the effectiveness of the proposed method on various conditional learning tasks, e.g., class-to-image generation, image-to-image translation, and image recovery. The advantage of our method over GAN-based methods is that our method is equipped with a slow thinking process that refines the solution guided by a learned objective function.


title: A Method to Search for Black Hole Candidates with Giant Companions by LAMOST

abstract: We propose a method to search for stellar-mass black hole (BH) candidates with giant companions from spectroscopic observations. Based on the stellar spectra of LAMOST Data Release 6, we obtain a sample of seven giants in binaries with large radial velocity variation $\Delta V_R > 80~{\rm km~s^{-1}}$. With the effective temperature, surface gravity, and metallicity provided by LAMOST, and the parallax given by {\it Gaia}, we can estimate the mass and radius of the giant, and therefore evaluate the possible mass of the optically invisible star in the binary. We show that the sources in our sample are potential BH candidates, and are worthy of dynamical measurement by further spectroscopic observations. Our method may be particularly valid for the selection of BH candidates in binaries with unknown orbital periods.


title: Hybrid magnetized stars within the Field Correlator Method

abstract: We present preliminary results of the study of intense magnetic fields effects on hybrid stars. For the description of the hadronic phase, we use the relativistic mean field approximation. For the quark matter phase, we employ the Field Correlator Method formalism. Once the the equation of state is built, we discuss the pressure anisotropy due the presence of the strong magnetic field. Finally, we calculate the structure of the compact stars using magnetized hybrid equations of state and their oscillation modes related with the emission of gravitational waves.


title: Discerning quantum memories based on electromagnetically-induced-transparency and Autler-Townes-splitting protocols

abstract: Electromagnetically induced transparency (EIT) and Autler-Townes splitting (ATS) are similar,but different quantum optical phenomena: EIT results from a Fano interference, whereas ATS is described by the AC-Stark effect. Likewise, despite their close resemblance, light-storage techniques based on the EIT memory protocol and the recently-proposed ATS memory protocol (E. Saglamyurek et al. Nature Photonics 12, 2018) are distinct: the EIT protocol relies on adiabatic elimination of absorption, whereas the ATS protocol is based on absorption. In this article, we elaborate on the distinction between EIT and ATS memory protocols through numerical analysis and experimental demonstrations in a cold rubidium ensemble. We find that their storage characteristics manifest opposite limits of the light-matter interaction due to their inherent adiabatic vs. non-adiabatic nature. Furthermore, we determine optimal memory conditions for each protocol and analyze ambiguous regimes in the case of broadband storage, where non-optimal memory implementations can possess characteristics of both EIT and ATS protocols. We anticipate that this investigation will lead to deeper understanding and improved technical development of quantum memories, while clarifying distinctions between the EIT and ATS protocols.


title: Revec: Program Rejuvenation through Revectorization

abstract: Modern microprocessors are equipped with Single Instruction Multiple Data (SIMD) or vector instructions which expose data level parallelism at a fine granularity. Programmers exploit this parallelism by using low-level vector intrinsics in their code. However, once programs are written using vector intrinsics of a specific instruction set, the code becomes non-portable. Modern compilers are unable to analyze and retarget the code to newer vector instruction sets. Hence, programmers have to manually rewrite the same code using vector intrinsics of a newer generation to exploit higher data widths and capabilities of new instruction sets. This process is tedious, error-prone and requires maintaining multiple code bases. We propose Revec, a compiler optimization pass which revectorizes already vectorized code, by retargeting it to use vector instructions of newer generations. The transformation is transparent, happening at the compiler intermediate representation level, and enables performance portability of hand-vectorized code. Revec can achieve performance improvements in real-world performance critical kernels. In particular, Revec achieves geometric mean speedups of 1.160$\times$ and 1.430$\times$ on fast integer unpacking kernels, and speedups of 1.145$\times$ and 1.195$\times$ on hand-vectorized x265 media codec kernels when retargeting their SSE-series implementations to use AVX2 and AVX-512 vector instructions respectively. We also extensively test Revec's impact on 216 intrinsic-rich implementations of image processing and stencil kernels relative to hand-retargeting.


title: A 30 GHz electro-optic frequency comb spanning 300 THz in the near infrared and visible

abstract: Beginning with a continuous wave laser at 1064 nm, we generate a 30 GHz electro-optic frequency comb which contains 100 lines spanning 3 THz. The initial comb is subsequently amplified, spectrally broadened in normal dispersion photonic crystal fiber, and then temporally compressed to provide 74 fs pulses with average power of up to 2.6 W. When launched into a second photonic crystal fiber with anomalous dispersion, a supercontinuum spanning 800-1350 nm is generated. Second harmonic generation allows for extension of the 30 GHz comb into the visible, yielding greater than 300 THz of total spectral bandwidth. Such a broad bandwidth, high repetition rate comb is a compelling source for astronomical spectrograph calibration.


title: An analytic semigroup generated by a fractional differential operator

abstract: We study a space-fractional diffusion problem, where the non-local diffusion flux involves the Caputo derivative of the diffusing quantity. We prove the unique existence of regular solutions to this problem by means of the semigroup theory. We show that the operator defined as divergence of product of a positive function with the Caputo derivative is a generator of analytic semigroup.


title: Appendix To Limits For Partial Maxima Of Gaussian Random Vectors

abstract: This appendix provides a short proof for sample path continuity of the Brownian motion induced by an arbitrary centered Gaussian measure on a separable Banach space, and also some perturbation results for the spectrum of compact self-adjoint operators on a Hilbert space.


title: Quenching low-mass satellite galaxies: evidence for a threshold ICM density

abstract: We compile a sample of SDSS galaxy clusters with high-quality Chandra X-ray data to directly study the influence of the dense intra-cluster medium (ICM) on the quenching of satellite galaxies. We study the quenched fractions of satellite galaxies as a function of ICM density for low- ($10^9 \lesssim M_\star \lesssim 10^{10}\,\mathrm{M_\odot}$), intermediate- ($10^{10} \lesssim M_\star \lesssim 10^{10.5}\,\mathrm{M_\odot}$), and high-mass ($M_\star \gtrsim 10^{10.5}\,\mathrm{M_\odot}$) satellite galaxies with $>\!3000$ satellite galaxies across 24 low-redshift ($z < 0.1$) clusters. For low-mass galaxies we find evidence for a broken powerlaw trend between satellite quenched fraction and local ICM density. The quenched fraction increases modestly at ICM densities below a threshold before increasing sharply beyond this threshold toward the cluster center. We show that this increase in quenched fraction at high ICM density is well matched by a simple, analytic model of ram pressure stripping. These results are consistent with a picture where low-mass cluster galaxies experience an initial, slow-quenching mode driven by steady gas depletion, followed by rapid quenching associated with ram pressure of cold-gas stripping near (one quarter of the virial radius, on average) the cluster center.


title: Sonine formulas and intertwining operators in Dunkl theory

abstract: Let $V_k$ denote Dunkl's intertwining operator associated with some root system $R$ and multiplicity function $k$. For two multiplicities $k, k^\prime$ on $R$, we study the operator $V_{k^\prime,k} = V_{k^\prime}\circ V_k^{-1}$, which intertwines the Dunkl operators for multiplicity $k$ with those for multiplicity $k^\prime.$ While it is well-known that the operator $V_k$ is positive for nonnegative $k$, it has been a long-standing conjecture that its generalizations $V_{k^\prime,k}$ are also positive if $k^\prime \geq k \geq 0,$ which is known to be true in rank one. In this paper, we disprove this conjecture by constructing examples for root system $B_n$ with multiplicites $k^\prime \geq k \geq 0$ for which $V_{k^\prime, k}$ is not positive. This matter is closely related to the existence of integral representations of Sonine type between the Dunkl kernels and Bessel functions associated with the relevant multiplicities. In our examples, such Sonine formulas do not exist. As a consequence, we obtain necessary conditions on Sonine-type integral formulas for Heckman-Opdam hypergeometric functions of type $BC_n$ as well as conditions on the existence of positive branching coefficients between systems of multivariable Jacobi polynomials.


title: CORSIKA 8 - Towards a modern framework for the simulation of extensive air showers

abstract: Current and future challenges in astroparticle physics require novel simulation tools to achieve higher precision and more flexibility. For three decades the FORTRAN version of CORSIKA served the community in an excellent way. However, the effort to maintain and further develop this complex package is getting increasingly difficult. To overcome existing limitations, and designed as a very open platform for all particle cascade simulations in astroparticle physics, we are developing CORSIKA 8 based on modern C++ and Python concepts. Here, we give a brief status report of the project.


title: Compatible Natural Gradient Policy Search

abstract: Trust-region methods have yielded state-of-the-art results in policy search. A common approach is to use KL-divergence to bound the region of trust resulting in a natural gradient policy update. We show that the natural gradient and trust region optimization are equivalent if we use the natural parameterization of a standard exponential policy distribution in combination with compatible value function approximation. Moreover, we show that standard natural gradient updates may reduce the entropy of the policy according to a wrong schedule leading to premature convergence. To control entropy reduction we introduce a new policy search method called compatible policy search (COPOS) which bounds entropy loss. The experimental results show that COPOS yields state-of-the-art results in challenging continuous control tasks and in discrete partially observable tasks.


title: Technical efficiency and inefficiency: SFA misspecification

abstract: The effect of external factors $z$ on technical inefficiency ($TI$) in stochastic frontier (SF) production models is often specified through the variance of inefficiency term $u$. In this setup the signs of marginal effects of $z$ on $TI$ and technical efficiency $TE$ identify how one should control $z$ to increase $TI$ or decrease $TE$. We prove that these signs for $TI$ and $TE$ are opposite for typical setups with normally distributed random error $v$ and exponentially or half-normally distributed $u$ for both conditional and unconditional cases. On the other hand, we give an example to show that signs of the marginal effects of $z$ on $TI$ and $TE$ may coincide, at least for some ranges of $z$. In our example, the distribution of $u$ is a mixture of two distributions, and the proportion of the mixture is a function of $z$. Thus if the real data comes from this mixture distribution, and we estimate model parameters with an exponential or half-normal distribution for $u$, the estimated efficiency and the marginal effect of $z$ on $TE$ would be wrong. Moreover, for a misspecified model, the rank correlations between the true and the estimated values of TE could be small and even negative for some subsamples of data. These results are demonstrated by simulations.


title: Best Practices for Administering Attitude and Beliefs Surveys

abstract: Physics faculty care about their students learning physics content. In addition, they usually hope that their students will learn some deeper lessons about thinking critically and scientifically. They hope that as a result of taking a physics class, students will come to appreciate physics as a coherent and logical method of understanding the world, and recognize that they can use reason and experimentation to figure things out about the world. Physics education researchers have created several surveys to assess one important aspect of thinking like a physicist: what students believe that learning physics is all about. In this article, we introduce attitudes and beliefs surveys; and give advice on how to choose, administer, and score them in your classes. This article is a companion to Best Practices for Administering Concept Inventories (The Physics Teacher, 2017), which introduces and answers common questions around concept inventories, which are research-based assessments of physics content topics.


title: Robustness Of Saak Transform Against Adversarial Attacks

abstract: Image classification is vulnerable to adversarial attacks. This work investigates the robustness of Saak transform against adversarial attacks towards high performance image classification. We develop a complete image classification system based on multi-stage Saak transform. In the Saak transform domain, clean and adversarial images demonstrate different distributions at different spectral dimensions. Selection of the spectral dimensions at every stage can be viewed as an automatic denoising process. Motivated by this observation, we carefully design strategies of feature extraction, representation and classification that increase adversarial robustness. The performances with well-known datasets and attacks are demonstrated by extensive experimental evaluations.


title: Modeling and Control of Soft Robots Using the Koopman Operator and Model Predictive Control

abstract: Controlling soft robots with precision is a challenge due in large part to the difficulty of constructing models that are amenable to model-based control design techniques. Koopman Operator Theory offers a way to construct explicit linear dynamical models of soft robots and to control them using established model-based linear control methods. This method is data-driven, yet unlike other data-driven models such as neural networks, it yields an explicit control-oriented linear model rather than just a "black-box" input-output mapping. This work describes this Koopman-based system identification method and its application to model predictive controller design. A model and MPC controller of a pneumatic soft robot arm was constructed via the method, and its performance was evaluated over several trajectory following tasks in the real-world. On all of the tasks, the Koopman-based MPC controller outperformed a benchmark MPC controller based on a linear state-space model of the same system.


title: All-electron fully relativistic Kohn-Sham theory for solids based on the Dirac-Coulomb Hamiltonian and Gaussian-type functions

abstract: We present the first full-potential method that solves the fully relativistic 4-component Dirac-Kohn-Sham equation for materials in the solid state within the framework of atom-centered Gaussian-type orbitals (GTOs). Our GTO-based method treats one-, two-, and three-dimensional periodic systems on an equal footing, and allows for a seamless transition to the methodology commonly used in studies of molecules with heavy elements. The scalar relativistic effects and spin-orbit interaction are handled variationally. The full description of the electron-nuclear potential in the core region of heavy nuclei is straightforward due to the local nature of the GTOs and does not pose any computational difficulties. We show how the time-reversal symmetry and a quaternion algebra-based formalism can be exploited to significantly reduce the increased methodological complexity and computational cost associated with multiple wave-function components coupled by the spin-orbit interaction. We provide a description of how to employ the matrix form of the multipole expansion and an iterative renormalization procedure to evaluate the conditionally convergent lattice sums arising for periodic systems. We investigate the problem of inverse variational collapse that arises if the Dirac operator containing a repulsive periodic potential is expressed in a basis containing diffuse functions, and suggest a possible solution. We demonstrate the validity of the method on silver halide crystals with large relativistic effects, and two-dimensional honeycomb structures (silicene and germanene) exhibiting the spin-orbit-driven quantum spin Hall effect. Our results are well-converged with respect to the basis limit using standard bases developed for molecular calculations, and indicate that the rule of removing basis functions with small exponents should not be applied when transferring the molecular basis to solids.


title: Low-cost Measurement of Industrial Shock Signals via Deep Learning Calibration

abstract: Special high-end sensors with expensive hardware are usually needed to measure shock signals with high accuracy. In this paper, we show that cheap low-end sensors calibrated by deep neural networks are also capable to measure high-g shocks accurately. Firstly we perform drop shock tests to collect a dataset of shock signals measured by sensors of different fidelity. Secondly, we propose a novel network to effectively learn both the signal peak and overall shape. The results show that the proposed network is capable to map low-end shock signals to its high-end counterparts with satisfactory accuracy. To the best of our knowledge, this is the first work to apply deep learning techniques to calibrate shock sensors.


title: Planarized Fabrication Process With Two Layers of SIS Josephson Junctions and Integration of SIS and SFS {\pi}-Junctions

abstract: We present our new fabrication Process for Superconductor Electronics (PSE2) that integrates two (2) layers of Josephson junctions in a fully planarized multilayer process on 200-mm wafers. The two junction layers can be, e.g., conventional Superconductor-Insulator-Superconductor (SIS) Nb/Al/AlO_x/Nb junctions with the same or different Josephson critical current densities, J_c. The process also allows integration of high-J_c Superconductor-Ferromagnet-Superconductor (SFS) or SFS'S JJs on the first junction layer with Nb/Al/AlO_x/Nb trilayer junctions on the second junction layer, or vice versa. In the present node, the SFS trilayer, Nb/Ni/Nb is placed below the standard SIS trilayer and separated by one niobium wiring layer. The main purpose of integrating the SFS and SIS junction layers is to provide compact {\pi}-phase shifters in logic cells of superconductor digital circuits and random access memories, and thereby increase the integration scale and functional density of superconductor electronics. The current node of the two-junction-layer process has six planarized niobium layers, two layers of resistors, and 350-nm minimum feature size. The target Josephson critical current densities for the SIS junctions are 100 {\mu}A/{\mu}m^2 and 200 {\mu}A/{\mu}m^2. We present the salient features of the new process, fabrication details, and characterization results on two layers of Josephson junctions integrated into one process, both for the conventional and {\pi}-junctions.


title: Evaluating Crowd Density Estimators via Their Uncertainty Bounds

abstract: In this work, we use the Belief Function Theory which extends the probabilistic framework in order to provide uncertainty bounds to different categories of crowd density estimators. Our method allows us to compare the multi-scale performance of the estimators, and also to characterize their reliability for crowd monitoring applications requiring varying degrees of prudence.


title: CO luminosity-Linewidth correlation of low and high redshift galaxies and its possible cosmological utilization

abstract: A linear correlation has been proposed between the CO luminosity ($\rm{L}^{\prime}_{\rm{CO}}$) and full-width at half maximum (FWHM) for high-redshift (z > 1) submillimeter galaxies. However, the controversy concerning the $\rm{L}^{\prime}_{\rm{CO}}$-FWHM correlation seems to have been caused by the use of heterogeneous samples (e.g., different transition lines) and/or data with large measurement uncertainties. In order to avoid the uncertainty caused by using different rotational transitions, in this work we make an extensive effort to select only CO($J = 1-0$) data from the literature. We separate these wide-ranging redshift data into two samples : the low-redshift (z < 1) and high-redshift (z > 1) samples. The samples are corrected for lensing magnification factors if gravitational-lensing effects appeared in the observations. The correlation analysis shows that there exists significant $\rm{L}^{\prime}_{\rm{CO}}$-FWHM correlations for both the low-redshift and high-redshift samples. A comparison of the low- and high-redshift $\rm{L}^{\prime}_{\rm{CO}}$-FWHM correlations does not show strong evolution with redshift. Assuming that there is no evolution, we can use this relation to determine the model independent distances of high-redshift galaxies. We then constrain cosmological models with the calibrated high-redshift CO data and the sample of Type Ia supernovae in the Union 2.1 compilation. In the constraint for wCDM with our samples, the derived values are w_{0} = -1.02 {\pm} 0.17, {\Omega}_{m0} = 0.30{\pm}0.02, and H_{0} = 70.00 {\pm}0.60 km\,s^{-1}\,Mpc^{-1}.


title: Exponential ergodicity for stochastic equations of nonnegative processes with jumps

abstract: In this work, we study ergodicity of continuous time Markov processes on state space $\mathbb{R}_{\geq 0} := [0,\infty)$ obtained as unique strong solutions to stochastic equations with jumps. Our first main result establishes exponential ergodicity in the Wasserstein distance, provided the stochastic equation satisfies a comparison principle and the drift is dissipative. In particular, it is applicable to continuous-state branching processes with immigration (shorted as CBI processes), possibly with nonlinear branching mechanisms or in L\'evy random environments. Our second main result establishes exponential ergodicity in total variation distance for subcritical CBI processes under a first moment condition on the jump measure for branching and a $\log$-moment condition on the jump measure for immigration.


title: The Long and the Short of It: Summarising Event Sequences with Serial Episodes

abstract: An ideal outcome of pattern mining is a small set of informative patterns, containing no redundancy or noise, that identifies the key structure of the data at hand. Standard frequent pattern miners do not achieve this goal, as due to the pattern explosion typically very large numbers of highly redundant patterns are returned. We pursue the ideal for sequential data, by employing a pattern set mining approach-an approach where, instead of ranking patterns individually, we consider results as a whole. Pattern set mining has been successfully applied to transactional data, but has been surprisingly under studied for sequential data. In this paper, we employ the MDL principle to identify the set of sequential patterns that summarises the data best. In particular, we formalise how to encode sequential data using sets of serial episodes, and use the encoded length as a quality score. As search strategy, we propose two approaches: the first algorithm selects a good pattern set from a large candidate set, while the second is a parameter-free any-time algorithm that mines pattern sets directly from the data. Experimentation on synthetic and real data demonstrates we efficiently discover small sets of informative patterns.


title: Threshold effects in heavy quarkonium spectroscopy

abstract: In this contribution, we discuss the possible importance of continuum-coupling (or threshold) effects in heavy quarkonium spectroscopy. Our calculations are carried out in a coupled-channel model, where meson-meson higher Fock (or molecular-type) components are introduced in $Q \bar Q$ bare meson wave functions by means of a pair-creation mechanism. After providing a quick resume of the main characteristics of the coupled-channel model, we briefly discuss its application to the calculation of the masses of heavy quarkonium-like $\chi_{\rm c}(2P)$ and $\chi_{\rm b}(3P)$ states with threshold corrections. We show that the introduction of pair-creation effects in the Quark Model (QM) formalism makes it possible to explain the deviations of $\chi_{\rm c}(2P)$ states' masses from the experimental data, without affecting the good QM description of the properties of $\chi_{\rm b}(3P)$ states.


title: A Novel Graph Analytic Approach to Monitor Terrorist Networks

abstract: Terrorist attacks all across the world have become a major source of concern for almost all national governments. The United States Department of State's Bureau of Counter-Terrorism, maintains a list of 66 terrorist organizations spanning the entire world. Actively monitoring a large number of organizations and their members, require considerable amounts of resources on the part of law enforcement agencies. Oftentimes, the law enforcement agencies do not have adequate resources to monitor these organizations and their members effectively. On multiple incidences of terrorist attacks in recent times across Europe, it has been observed that the perpetrators of the attack were in the suspect databases of the law enforcement authorities, but weren't under active surveillance at the time of the attack, due to resource limitations on the part of the authorities. As the suspect databases in various countries are very large, and it takes significant amount of technical and human resources to monitor a suspect in the database, monitoring all the suspects in the database may be an impossible task. In this paper, we propose a novel terror network monitoring approach that will significantly reduce the resource requirement of law enforcement authorities, but still provide the capability of uniquely identifying a suspect in case the suspect becomes active in planning a terrorist attack. The approach relies on the assumption that, when an individual becomes active in planning a terrorist attack, his/her friends/associates will have some inkling of the individuals plan. Accordingly, even if the individual is not under active surveillance by the authorities, but the individual's friends/associates are, then the individual planning the attack can be uniquely identified. We apply our techniques on various real-world terror network datasets and show the effectiveness of our approach.


title: Accuracy vs. Computational Cost Tradeoff in Distributed Computer System Simulation

abstract: Simulation is a fundamental research tool in the computer architecture field. These kinds of tools enable the exploration and evaluation of architectural proposals capturing the most relevant aspects of the highly complex systems under study. Many state-of-the-art simulation tools focus on single-system scenarios, but the scalability required by trending applications has shifted towards distributed computing systems integrated via complex software stacks. Web services with client-server architectures or distributed storage and processing of scale-out data analytics (Big Data) are among the main exponents. The complete simulation of a distributed computer system is the appropriate methodology to conduct accurate evaluations. Unfortunately, this methodology could have a significant impact on the already large computational effort derived from detailed simulation. In this work, we conduct a set of experiments to evaluate this accuracy/cost tradeoff. We measure the error made if client-server applications are evaluated in a single-node environment, as well as the overhead induced by the methodology and simulation tool employed for multi-node simulations. We quantify this error for different micro-architecture components, such as last-level cache and instruction/data TLB. Our findings show that accuracy loss can lead to completely wrong conclusions about the effects of proposed hardware optimizations. Fortunately, our results also demonstrate that the computational overhead of a multi-node simulation framework is affordable, suggesting multi-node simulation as the most appropriate methodology.


title: W algebra, Cosets and VOA for 4d N = 2 SCFT from M5 branes

abstract: We identify vertex operator algebras (VOAs) of a class of Argyres-Douglas (AD) matters with two types of non-abelian flavor symmetries. They are the $W$ algebra defined using nilpotent orbit with partition $[q^m,1^s]$. Gauging above AD matters, we can find VOAs for more general $\mathcal{N}=2$ SCFTs engineered from 6d $(2,0)$ theories. For example, the VOA for general $(A_{N-1}, A_{k-1})$ theory is found as the coset of a collection of above $W$ algebras. Various new interesting properties of 2d VOAs such as level-rank duality, conformal embedding, collapsing levels, coset constructions for known VOAs can be derived from 4d theory.


title: Higher Steenrod squares for Khovanov homology

abstract: We describe stable cup-i products on the cochain complex with $F^2$ coefficients of any augmented semi-simplicial object in the Burnside category. An example of such an object is the Khovanov functor of Lawson, Lipshitz and Sarkar. Thus we obtain explicit formulas for cohomology operations on the Khovanov homology of any link.


title: Higher order corks

abstract: It is shown that any finite list of smooth closed simply-connected 4-manifolds homeomorphic to a given one X can be obtained by removing a single compact contractible submanifold (or cork) from X, and then regluing it by powers of a boundary diffeomorphism. We then use this result to "separate" finite families of corks embedded in a fixed 4-manifold.


title: 3D Human Pose Estimation from Deep Multi-View 2D Pose

abstract: Human pose estimation - the process of recognizing a human's limb positions and orientations in a video - has many important applications including surveillance, diagnosis of movement disorders, and computer animation. While deep learning has lead to great advances in 2D and 3D pose estimation from single video sources, the problem of estimating 3D human pose from multiple video sensors with overlapping fields of view has received less attention. When the application allows use of multiple cameras, 3D human pose estimates may be greatly improved through fusion of multi-view pose estimates and observation of limbs that are fully or partially occluded in some views. Past approaches to multi-view 3D pose estimation have used probabilistic graphical models to reason over constraints, including per-image pose estimates, temporal smoothness, and limb length. In this paper, we present a pipeline for multi-view 3D pose estimation of multiple individuals which combines a state-of-art 2D pose detector with a factor graph of 3D limb constraints optimized with belief propagation. We evaluate our results on the TUM-Campus and Shelf datasets for multi-person 3D pose estimation and show that our system significantly out-performs the previous state-of-the-art with a simpler model of limb dependency.


title: Community Animation: Exploring a design space that leverages geosocial networking to increase community engagement

abstract: This paper explores a design study of a smartphone enabled meet-up app meant to inspire engagement in community innovation. Community hubs such as co-working spaces, incubators, and maker spaces attract community members with diverse interests. This paper presents these spaces as a design opportunity for an application that helps host community-centered meet-ups in smart and connected communities. Our design study explores three scenarios of use, inspired by previous literature, for organizing meet-ups and compares them by surveying potential users. Based on the results of our survey, we propose several design implications and implement them in the Community Animator geosocial networking application, which identifies nearby individuals that are willing to chat or perform community-centered activities. We present the results of both our survey and our prototype, discuss our design goals, and provide design implications for civic-minded, geosocial networking applications. Our contribution in this work is the development process, proposed design of a mobile application to support community-centered meet-ups, and insights for future work.


title: Stable maps, Q-operators and category O

abstract: Motivated by Maulik-Okounkov stable maps associated to quiver varieties, we define and construct algebraic stable maps on tensor products of representations in the category O of the Borel subalgebra of an arbitrary untwisted quantum affine algebra. Our representation-theoretical construction is based on the study of the action of Cartan-Drinfeld subalgebras. We prove the algebraic stable maps are invertible and depend rationally on the spectral parameter. As an application, we obtain new R-matrices in the category O and we establish that a large family of simple modules, including the prefundamental representations associated to Q-operators, generically commute as representations of the Cartan-Drinfeld subalgebra. We also establish categorified QQ*-systems in terms of the R-matrices we construct.


title: Entanglement Wedge Reconstruction using the Petz Map

abstract: At the heart of recent progress in AdS/CFT is the question of subregion duality, or entanglement wedge reconstruction: which part(s) of the boundary CFT are dual to a given subregion of the bulk? This question can be answered by appealing to the quantum error correcting properties of holography, and it was recently shown that robust bulk (entanglement wedge) reconstruction can be achieved using a universal recovery channel known as the twirled Petz map. In short, one can use the twirled Petz map to recover bulk data from a subset of the boundary. However, this map involves an averaging procedure over bulk and boundary modular time, and hence it can be somewhat intractable to evaluate in practice. We show that a much simpler channel, the Petz map, is sufficient for entanglement wedge reconstruction for any code space of fixed finite dimension - no twirling is required. Moreover, the error in the reconstruction will always be non-perturbatively small. From a quantum information perspective, we prove a general theorem extending the use of the Petz map as a general-purpose recovery channel to subsystem and operator algebra quantum error correction.


title: FaceSpoof Buster: a Presentation Attack Detector Based on Intrinsic Image Properties and Deep Learning

abstract: Nowadays, the adoption of face recognition for biometric authentication systems is usual, mainly because this is one of the most accessible biometric modalities. Techniques that rely on trespassing these kind of systems by using a forged biometric sample, such as a printed paper or a recorded video of a genuine access, are known as presentation attacks, but may be also referred in the literature as face spoofing. Presentation attack detection is a crucial step for preventing this kind of unauthorized accesses into restricted areas and/or devices. In this paper, we propose a novel approach which relies in a combination between intrinsic image properties and deep neural networks to detect presentation attack attempts. Our method explores depth, salience and illumination maps, associated with a pre-trained Convolutional Neural Network in order to produce robust and discriminant features. Each one of these properties are individually classified and, in the end of the process, they are combined by a meta learning classifier, which achieves outstanding results on the most popular datasets for PAD. Results show that proposed method is able to overpass state-of-the-art results in an inter-dataset protocol, which is defined as the most challenging in the literature.


title: Effects of non-pairwise repulsion on nanoparticle assembly

abstract: Electrostatic interactions provide a convenient way to modulate interactions between nanoparticles, colloids, and biomolecules because they can be adjusted by the solution pH or salt concentration. While the presence of salt provides an easy method to control the net interparticle interaction, the nonlinearities arising from electrostatic screening make it difficult to quantify the strength of the interaction. In particular, when charged particles assemble into clusters or aggregates, nonlinear effects render the interactions strongly non-pairwise. Here we report Brownian dynamics simulations to investigate the effect that the non-pairwise nature of electrostatic interactions has on nanoparticle assembly. We compare these simulations to a system in which the electrostatics are modeled by a strictly pairwise Yukawa potential. We find that both systems show a narrow range in parameter space where the particles form well-ordered crystals. Bordering this range are regions where the net interactions are too weak to stabilize aggregated structures, or strong enough that the system becomes kinetically trapped in a gel. The non-pairwise potential differs from the pairwise system in the appearance of an amorphous phase for strongly charged particles. This phase appears because the many-body electrostatic interactions limit the maximum density achievable in an assembly.


title: Regularity for Fully Nonlinear Parabolic Equations with Oblique Boundary Data

abstract: We obtain up to a flat boundary regularity results in parabolic H\"{o}lder spaces for viscosity solutions of fully nonlinear parabolic equations with oblique boundary conditions.


title: Conditionally Free Reduced Products of Hilbert Spaces

abstract: We present a product of pairs of pointed Hilbert spaces that, in the context of Boz\.ejko, Leinert and Speicher's theory of conditionally free probability, plays the role of the reduced free product of pointed Hilbert spaces, and thus gives a unified construction for the natural notions of independence defined by Muraki. We additionally provide important applications of this construction. We prove that, assuming minor restrictions, for any pair of conditionally free algebras there are copies of them that are conditionally free and also free, a property that is frequently assumed (as hypothesis) to prove several results in the literature. Finally, we give a short proof of the linearization property of the $^cR$-transform (the analog of Voiculescu's $R$-transform in the context of conditionally free probability).


title: f(R) and Brans-Dicke Theories and the Swampland

abstract: We discuss the viability of $f(R)$ and Brans-Dicke theories of gravity in light of the recent Swampland conjectures. We show that in the case of perfect fluid domination the swampland conjecture conditions may be easily satisfied and therefore we focus on the constraints in the vacuum theory. We derive the conditions for the swampland conjectures in $f(R)$ and Brans-Dicke framework, and find that a large portion of the parameter space is not allowed. Furthermore, we analyze particular $f(R)$ models of inflation and dark energy and in their case we set the allowed range of $R$.


title: Using Deep Q-learning To Prolong the Lifetime of Correlated Internet of Things Devices

abstract: Battery-powered sensors deployed in the Internet of Things (IoT) require energy-efficient solutions to prolong their lifetime. When these sensors observe a physical phenomenon distributed in space and evolving in time, the collected observations are expected to be correlated. We take advantage of the exhibited correlation and propose an updating mechanism that employs deep Q-learning. Our mechanism is capable of determining the frequency with which sensors should transmit their updates while taking into the consideration an ever-changing environment. We evaluate our solution using observations obtained in a real deployment, and show that our proposed mechanism is capable of significantly extending battery-powered sensors' lifetime without compromising the accuracy of the observations provided to the IoT service.


title: Towards Provably Not-at-Fault Control of Autonomous Robots in Arbitrary Dynamic Environments

abstract: As autonomous robots increasingly become part of daily life, they will often encounter dynamic environments while only having limited information about their surroundings. Unfortunately, due to the possible presence of malicious dynamic actors, it is infeasible to develop an algorithm that can guarantee collision-free operation. Instead, one can attempt to design a control technique that guarantees the robot is not-at-fault in any collision. In the literature, making such guarantees in real time has been restricted to static environments or specific dynamic models. To ensure not-at-fault behavior, a robot must first correctly sense and predict the world around it within some sufficiently large sensor horizon (the prediction problem), then correctly control relative to the predictions (the control problem). This paper addresses the control problem by proposing Reachability-based Trajectory Design for Dynamic environments (RTD-D), which guarantees that a robot with an arbitrary nonlinear dynamic model correctly responds to predictions in arbitrary dynamic environments. RTD-D first computes a Forward Reachable Set (FRS) offline of the robot tracking parameterized desired trajectories that include fail-safe maneuvers. Then, for online receding-horizon planning, the method provides a way to discretize predictions of an arbitrary dynamic environment to enable real-time collision checking. The FRS is used to map these discretized predictions to trajectories that the robot can track while provably not-at-fault. One such trajectory is chosen at each iteration, or the robot executes the fail-safe maneuver from its previous trajectory which is guaranteed to be not at fault. RTD-D is shown to produce not-at-fault behavior over thousands of simulations and several real-world hardware demonstrations on two robots: a Segway, and a small electric vehicle.


title: On Asymptotically Tight Tail Bounds for Sums of Geometric and Exponential Random Variables

abstract: In this note we prove bounds on the upper and lower probability tails of sums of independent geometric or exponentially distributed random variables. We also prove negative results showing that our established tail bounds are asymptotically tight.


title: CMC foliations of open spacetimes asymptotic to open Robertson-Walker spacetimes

abstract: We consider open globally hyperbolic spacetimes $N$ of dimension $n+1$, $n\ge 3$, which are spatially asymptotic to a Robertson-Walker spacetime or an open Friedmann universe with spatial curvature $\tilde\kappa = 0,-1$ and prove, under reasonable assumptions, that there exists a unique foliation by hypersurfaces of constant mean curvature and that the mean curvature function $\tau$ is a smooth time function if $N$ is smooth. Moreover, among the Friedmann universes which satisfy the necessary conditions are those that reflect the present assumptions of the development of the universe.


title: Static and semi-static hedging as contrarian or conformist bets

abstract: In this paper, we argue that, once the costs of maintaining the hedging portfolio are properly taken into account, semi-static portfolios should more properly be thought of as separate classes of derivatives, with non-trivial, model-dependent payoff structures. We derive new integral representations for payoffs of exotic European options in terms of payoffs of vanillas, different from Carr-Madan representation, and suggest approximations of the idealized static hedging/replicating portfolio using vanillas available in the market. We study the dependence of the hedging error on a model used for pricing and show that the variance of the hedging errors of static hedging portfolios can be sizably larger than the errors of variance-minimizing portfolios. We explain why the exact semi-static hedging of barrier options is impossible for processes with jumps, and derive general formulas for variance-minimizing semi-static portfolio. We show that hedging using vanillas only leads to larger errors than hedging using vanillas and first touch digitals. In all cases, efficient calculations of the weights of the hedging portfolios are in the dual space using new efficient numerical methods for calculation of the Wiener-Hopf factors and Laplace-Fourier inversion.


title: Integrated magnonic half-adder

abstract: Spin waves and their quanta magnons open up a promising branch of high-speed and low-power information processing. Several important milestones were achieved recently in the realization of separate magnonic data processing units including logic gates, a magnon transistor and units for non-Boolean computing. Nevertheless, the realization of an integrated magnonic circuit consisting of at least two logic gates and suitable for further integration is still an unresolved challenge. Here we demonstrate such an integrated circuit numerically on the example of a magnonic half-adder. Its key element is a nonlinear directional coupler serving as combined XOR and AND logic gate that utilizes the dependence of the spin wave dispersion on its amplitude. The circuit constitutes of only three planar nano-waveguides and processes all information within the magnon domain. Benchmarking of the proposed device is performed showing the potential for sub-aJ energy consumption per operation.


title: The ALMA Development Roadmap

abstract: The present document outlines a roadmap for future developments that will significantly expand ALMA's capabilities and enable it to produce even more exciting science in the coming decades. The proposed developments are motivated by the groundbreaking results achieved by ALMA during its first five years of operation. The roadmap described here is based on input on new scientific directions and technical feasibility of future developments from the ALMA Science Advisory Committee (ASAC), the community, and technical documents. The Working Group recommends that the top development priority, based on scientific merit and technical feasibility, is to broaden the receiver IF bandwidth by at least a factor two, and to upgrade the associated electronics and correlator. These developments will advance a wide range of scientific studies by significantly reducing the time required for blind redshift surveys, chemical spectral scans, and deep continuum surveys. In order of scientific priority, receiver upgrades are recommended for intermediate (200-425 GHz), low (< 200 GHz), and high (> 425 GHz) frequencies. The Working Group recommends that the receiver and throughput developments proceed as soon as fiscally and technically feasible. As a first step, a technical and scientific group should be formed to formalize the top-level requirements. A team of systems engineers should then be charged with flowing these requirements down to the subsystems to form a consistent new set of minimum requirements, which future development projects would have to meet. Given that upgrading the throughput will impact many ALMA subsystems, the Working Group recommends that a team within ALMA be charged with coordinating and monitoring these developments. (Abbreviated)


title: Pareto-optimal solution for the quantum battle of the sexes

abstract: Quantum games have gained much popularity in the last two decades. Many of these quantum games are a redefinition of iconic classical games to fit the quantum world, and they gain many different properties and solutions in this different view. In this letter, we attempt to find a solution to an asymmetric quantum game which still troubles quantum game researchers, the quantum battle of the sexes. To achieve that, we perform an analysis using the Eisert-Wilkens-Lewenstein's protocol for this asymmetric game. The protocol highlights two solutions for the game, which solve the dilemma and satisfy the Pareto-optimal definition, unlike previous reports that rely on Nash equilibrium. We perform an experimental implementation using the NMR technique in a two-qubit system. Our results eliminate dilemmas on the quantum battle of the sexes and provide us with arguments to elucidate that the Eisert-Wilkens-Lewenstein's protocol is not restricted to symmetric games when at the quantum regime.


title: Superionic UO2: A Model Anharmonic Crystalline Material

abstract: Crystalline materials at elevated temperatures and pressures can exhibit properties more reminiscent of simple liquids than ideal crystalline materials. Superionic crystalline materials having a liquid-like conductivity {\sigma} are particularly interesting for battery, fuel cell, and other energy applications, and we study UO2 as a prototypical superionic material since this material is widely studied given its commercial importance as reactor fuel. Using molecular dynamics, we first investigate basic thermodynamic and structural properties. We then quantify structural relaxation, dynamic heterogeneity, and average ions mobility. We find that the non-Arrhenius diffusion and structural relaxation time of this prototypical superionic material can be quantitatively described in terms of a generalized activated transport model ('string model') in which the activation energy varies in direct proportion to the average string length. Our transport data can also be described equally well by an Adam-Gibbs model in which the excess entropy density of the crystalline material is estimated from specific heat and thermal expansion data, consistent with the average scale of string-like collective motion scaling inversely with the excess entropy of the crystal. Strong differences in the temperature dependence of the interfacial mobility from non-ionic materials are observed, and we suggest that this difference is due to the relatively high cohesive interaction of ionic materials. In summary, the study of superionic UO2 provides insight into the role of cooperative motion in enhancing ion mobility in ionic materials and offers design principles for the development of new superionic materials for use in diverse energy applications.


title: Influence of Bonded Interactions on Structural Phases of Flexible Polymers

abstract: We introduce a novel coarse-grained bead-spring model for flexible polymers to systematically examine the effects of an adjusted bonded potential on the formation and stability of structural macrostates in a thermal environment. The density of states obtained in advanced replica-exchange Monte Carlo simulations is analyzed by employing the recently developed generalized microcanonical inflection-point analysis method, which enables the identification of diverse structural phases and the construction of a suitably parameterized hyperphase diagram. It reveals that icosahedral phases dominate for polymers with asymmetric and narrow bond potentials, whereas polymers with symmetric and more elastic bonds tend to form amorphous structures with non-icosahedral cores. We also observe a hierarchy in the freezing transition behavior associated with the formation of the surface layer after nucleation.


title: Crop Yield Prediction Using Deep Neural Networks

abstract: Crop yield is a highly complex trait determined by multiple factors such as genotype, environment, and their interactions. Accurate yield prediction requires fundamental understanding of the functional relationship between yield and these interactive factors, and to reveal such relationship requires both comprehensive datasets and powerful algorithms. In the 2018 Syngenta Crop Challenge, Syngenta released several large datasets that recorded the genotype and yield performances of 2,267 maize hybrids planted in 2,247 locations between 2008 and 2016 and asked participants to predict the yield performance in 2017. As one of the winning teams, we designed a deep neural network (DNN) approach that took advantage of state-of-the-art modeling and solution techniques. Our model was found to have a superior prediction accuracy, with a root-mean-square-error (RMSE) being 12% of the average yield and 50% of the standard deviation for the validation dataset using predicted weather data. With perfect weather data, the RMSE would be reduced to 11% of the average yield and 46% of the standard deviation. We also performed feature selection based on the trained DNN model, which successfully decreased the dimension of the input space without significant drop in the prediction accuracy. Our computational results suggested that this model significantly outperformed other popular methods such as Lasso, shallow neural networks (SNN), and regression tree (RT). The results also revealed that environmental factors had a greater effect on the crop yield than genotype.


title: Discovering Descriptive Tile Trees by Mining Optimal Geometric Subtiles

abstract: When analysing binary data, the ease at which one can interpret results is very important. Many existing methods, however, discover either models that are difficult to read, or return so many results interpretation becomes impossible. Here, we study a fully automated approach for mining easily interpretable models for binary data. We model data hierarchically with noisy tiles-rectangles with significantly different density than their parent tile. To identify good trees, we employ the Minimum Description Length principle. We propose STIJL, a greedy any-time algorithm for mining good tile trees from binary data. Iteratively, it finds the locally optimal addition to the current tree, allowing overlap with tiles of the same parent. A major result of this paper is that we find the optimal tile in only $\Theta(NM\min(N, M))$ time. STIJL can either be employed as a top-$k$ miner, or by MDL we can identify the tree that describes the data best. Experiments show we find succinct models that accurately summarise the data, and, by their hierarchical property are easily interpretable.


title: Lattices from tight frames and vertex transitive graphs

abstract: We show that real tight frames that generate lattices must be rational, and use this observation to describe a construction of lattices from vertex transitive graphs. In the case of irreducible group frames, we show that the corresponding lattice is always strongly eutactic. This is the case for the more restrictive class of distance transitive graphs. We show that such lattices exist in arbitrarily large dimensions and demonstrate examples arising from some notable families of graphs. In particular, some well-known root lattices and those related to them can be recovered this way. We discuss various properties of this construction and also mention some potential applications of lattices generated by incoherent systems of vectors.