title: Topological mixture estimation

abstract: Density functions that represent sample data are often multimodal, i.e. they exhibit more than one maximum. Typically this behavior is taken to indicate that the underlying data deserves a more detailed representation as a mixture of densities with individually simpler structure. The usual specification of a component density is quite restrictive, with log-concave the most general case considered in the literature, and Gaussian the overwhelmingly typical case. It is also necessary to determine the number of mixture components \emph{a priori}, and much art is devoted to this. Here, we introduce \emph{topological mixture estimation}, a completely nonparametric and computationally efficient solution to the one-dimensional problem where mixture components need only be unimodal. We repeatedly perturb the unimodal decomposition of Baryshnikov and Ghrist to produce a topologically and information-theoretically optimal unimodal mixture. We also detail a smoothing process that optimally exploits topological persistence of the unimodal category in a natural way when working directly with sample data. Finally, we illustrate these techniques through examples.

categories: stat.ME cs.IT math.AT math.IT math.ST stat.TH


title: Accuracy of the adiabatic-impulse approximation for closed and open quantum systems

abstract: We study the adiabatic-impulse approximation (AIA) as a tool to approximate the time evolution of quantum states, when driven through a region of small gap. The AIA originates from the Kibble-Zurek theory applied to continuous quantum phase transitions. The Kibble-Zurek mechanism was developed to predict the power-law scaling of the defect density across a continuous quantum phase transition. Instead here, we quantify the accuracy of the AIA via the trace norm distance with respect to the exact evolved state. As expected, we find that for short times/fast protocols, the AIA outperforms the simple adiabatic approximation. However, for large times/slow protocols, the situation is actually reversed and the AIA provides a worse approximation. Nevertheless, we found a variation of the AIA that can perform better than the adiabatic one. This counter-intuitive modification consists in crossing twice the region of small gap. Our findings are illustrated by several examples of driven closed and open quantum systems.


title: Fingerprint Spoof Buster

abstract: The primary purpose of a fingerprint recognition system is to ensure a reliable and accurate user authentication, but the security of the recognition system itself can be jeopardized by spoof attacks. This study addresses the problem of developing accurate, generalizable, and efficient algorithms for detecting fingerprint spoof attacks. Specifically, we propose a deep convolutional neural network based approach utilizing local patches centered and aligned using fingerprint minutiae. Experimental results on three public-domain LivDet datasets (2011, 2013, and 2015) show that the proposed approach provides state-of-the-art accuracies in fingerprint spoof detection for intra-sensor, cross-material, cross-sensor, as well as cross-dataset testing scenarios. For example, in LivDet 2015, the proposed approach achieves 99.03% average accuracy over all sensors compared to 95.51% achieved by the LivDet 2015 competition winners. Additionally, two new fingerprint presentation attack datasets containing more than 20,000 images, using two different fingerprint readers, and over 12 different spoof fabrication materials are collected. We also present a graphical user interface, called Fingerprint Spoof Buster, that allows the operator to visually examine the local regions of the fingerprint highlighted as live or spoof, instead of relying on only a single score as output by the traditional approaches.


title: SPIRou Input Catalogue: Global properties of 440 M dwarfs observed with ESPaDOnS at CFHT

abstract: Present and future high-precision radial-velocity spectrometers dedicated to the discovery of low-mass planets orbiting low-mass dwarfs need to focus on the best selected stars to make an efficient use of telescope time. In the framework of the preparation of the SPIRou Input Catalog, the CoolSnap program aims at screening M dwarfs in the solar neighborhood against binarity, rapid rotation, activity, ... To optimize the selection, the present paper describes the methods used to compute effective temperature, metallicity, projected rotation velocity of a large sample of 440 M dwarfs observed in the visible with the high-resolution spectro-polarimeter ESPaDOnS at CFHT. It also summarizes known and newly-discovered spectroscopic binaries, and stars known to belong to visual multiple systems. A calibration of the projected rotation velocity versus measured line widths for M dwarfs observed by the ESPaDOnS spectro-polarimeter is derived, and the resulting values are compared to equatorial rotation velocities deduced from rotation periods and radii. A comparison of the derived effective temperatures and metallicities with literature values is also conducted. Finally, the radial velocity uncertainty of each star in the sample is estimated, to narrow down the selection of stars to be included into the SPIRou Input Catalogue (SPIC).


title: Shintani theta lifts of harmonic Maass forms

abstract: We define a regularized Shintani theta lift which maps weight $2k+2$ ($k \in \Z, k \geq 0$) harmonic Maass forms for congruence subgroups to (sesqui-)harmonic Maass forms of weight $3/2+k$ for the Weil representation of an even lattice of signature $(1,2)$. We show that its Fourier coefficients are given by traces of CM values and regularized cycle integrals of the input harmonic Maass form. Further, the Shintani theta lift is related via the $\xi$-operator to the Millson theta lift studied in our earlier work. We use this connection to construct $\xi$-preimages of Zagier's weight $1/2$ generating series of singular moduli and of some of Ramanujan's mock theta functions.


title: Super-luminous Type II supernovae powered by magnetars

abstract: Magnetar power is believed to be at the origin of numerous super-luminous supernovae (SNe) of Type Ic, arising from compact, hydrogen-deficient, Wolf-Rayet type stars. Here, we investigate the properties that magnetar power would have on standard-energy SNe associated with 15-20Msun blue or red supergiant (BSG/RSG) stars. We use a combination of Eulerian grey radiation-hydrodynamics and non-LTE steady-state radiative transfer to study their dynamical, photometric, and spectroscopic properties. Adopting magnetar fields of 1, 3.5, 7 x 10^14G and rotational energies of 0.4, 1, and 3 x 10^51erg, we produce bolometric light curves with a broad maximum covering 50-150d and a magnitude of 10^43-10^44erg/s. The spectra at maximum light are analogous to those of standard SNe II-P but bluer. Although the magnetar energy is channelled roughly in equal proportion between SN kinetic energy and SN luminosity, the latter may be boosted by a factor 10-100 compared to a standard SN II. This influence breaks the observed relation between brightness and expansion rate of standard Type II SNe. Magnetar energy injection also delays recombination and may even cause re-ionization, with a reversal in photospheric temperature and velocity. Depositing the magnetar energy in a narrow mass shell at the ejecta base leads to the formation of a dense shell at a few 1000km/s, which causes a light-curve bump at the end of the photospheric phase. Depositing this energy over a broad range of mass in the inner ejecta, to mimic the effect of multi-dimensional fluid instabilities, prevents the formation of a dense shell and produces an earlier-rising and smoother light curve. The magnetar influence on the SN radiation is generally not visible prior to 20-30d, during which one may discern a BSG from a RSG progenitor. We propose a magnetar model for the super-luminous Type II SN OGLE-SN14-073.


title: A Data-Driven Sparse-Learning Approach to Model Reduction in Chemical Reaction Networks

abstract: In this paper, we propose an optimization-based sparse learning approach to identify the set of most influential reactions in a chemical reaction network. This reduced set of reactions is then employed to construct a reduced chemical reaction mechanism, which is relevant to chemical interaction network modeling. The problem of identifying influential reactions is first formulated as a mixed-integer quadratic program, and then a relaxation method is leveraged to reduce the computational complexity of our approach. Qualitative and quantitative validation of the sparse encoding approach demonstrates that the model captures important network structural properties with moderate computational load.


title: Self-consistency in the phonon space of the particle-phonon coupling model

abstract: In the paper the non-linear generalization of the time blocking approximation (TBA) is presented. The TBA is one of the versions of the extended random-phase approximation (RPA) developed within the Green-function method and the particle-phonon coupling model. In the generalized version of the TBA the self-consistency principle is extended onto the phonon space of the model. The numerical examples show that this non-linear version of the TBA leads to the convergence of the results with respect to enlarging the phonon space of the model.


title: Intra-node Memory Safe GPU Co-Scheduling

abstract: GPUs in High-Performance Computing systems remain under-utilised due to the unavailability of schedulers that can safely schedule multiple applications to share the same GPU. The research reported in this paper is motivated to improve the utilisation of GPUs by proposing a framework, we refer to as schedGPU, to facilitate intra-node GPU co-scheduling such that a GPU can be safely shared among multiple applications by taking memory constraints into account. Two approaches, namely a client-server and a shared memory approach are explored. However, the shared memory approach is more suitable due to lower overheads when compared to the former approach. Four policies are proposed in schedGPU to handle applications that are waiting to access the GPU, two of which account for priorities. The feasibility of schedGPU is validated on three real-world applications. The key observation is that a performance gain is achieved. For single applications, a gain of over 10 times, as measured by GPU utilisation and GPU memory utilisation, is obtained. For workloads comprising multiple applications, a speed-up of up to 5x in the total execution time is noted. Moreover, the average GPU utilisation and average GPU memory utilisation is increased by 5 and 12 times, respectively.


title: Computational Properties of Slime Trail

abstract: We investigate the combinatorial game Slime Trail.This game is played on a graph with a starting piece in a node. Each player's objective is to reach one of their own goal nodes. Every turn the current player moves the piece and deletes the node they came from. We show that the game is PSPACE-complete when played on a planar graph.


title: Nonunitary representations of the groups of $U(p,q)$-currents for $q\geq p>1$

abstract: The purpose of this paper is to give a construction of representations of the group of currents for semisimple groups of rank greater than one. Such groups have no unitary representations in the Fock space, since the semisimple groups of this form have no nontrivial cohomology in faithful irreducible representations. Thus we first construct cohomology of the semisimple groups in nonunitary representations. The principal method is to reduce all constructions to Iwasawa subgroups (solvable subgroups of the semisimple groups), with subsequent extension to the original group. The resulting representation is realized in the so-called quasi-Poisson Hilbert space associated with natural measures on infinite-dimensional spaces. Key words: Iwasawa subgroup, cohomology, group of currents, nonunitary representations.


title: Pixel-Level Statistical Analyses of Prescribed Fire Spread

abstract: Wildland fire dynamics is a complex turbulent dimensional process. Cellular automata (CA) is an efficient tool to predict fire dynamics, but the main parameters of the method are challenging to estimate. To overcome this challenge, we compute statistical distributions of the key parameters of a CA model using infrared images from controlled burns. Moreover, we apply this analysis to different spatial scales and compare the experimental results to a simple statistical model. By performing this analysis and making this comparison, several capabilities and limitations of CA are revealed.


title: Colloquium: Criticality and dynamical scaling in living systems

abstract: A celebrated and controversial hypothesis conjectures that some biological systems --parts, aspects, or groups of them-- may extract important functional benefits from operating at the edge of instability, halfway between order and disorder, i.e. in the vicinity of the critical point of a phase transition. Criticality has been argued to provide biological systems with an optimal balance between robustness against perturbations and flexibility to adapt to changing conditions, as well as to confer on them optimal computational capabilities, huge dynamical repertoires, unparalleled sensitivity to stimuli, etc. Criticality, with its concomitant scale invariance, can be conjectured to emerge in living systems as the result of adaptive and evolutionary processes that, for reasons to be fully elucidated, select for it as a template upon which higher layers of complexity can rest. This hypothesis is very suggestive as it proposes that criticality could constitute a general and common organizing strategy in biology stemming from the physics of phase transitions. However, despite its thrilling implications, this is still in its embryonic state as a well-founded theory and, as such, it has elicited some healthy skepticism. From the experimental side, the advent of high-throughput technologies has created new prospects in the exploration of biological systems, and empirical evidence in favor of criticality has proliferated, with examples ranging from endogenous brain activity and gene-expression patterns, to flocks of birds and insect-colony foraging, to name but a few...

categories: cond-mat.stat-mech nlin.AO physics.bio-ph q-bio.QM


title: A view of the H-band light-element chemical patterns in Globular Clusters under the AGB self-enrichment scenario

abstract: We discuss the self-enrichment scenario by AGB stars for the formation of multiple populations in globular clusters (GCs) by analyzing data set of giant stars observed in 9 Galactic GCs, covering a wide range of metallicities and for which the simultaneous measurements of C, N, O, Mg, Al, Si are available. To this aim we calculated 6 sets of AGB models, with the same chemical composition as the stars belonging to the first generation of each GC. We find that the AGB yields can reproduce the set of observations available, not only in terms of the degree of contamination shown by stars in each GC but, more important, also the observed trend with metallicity, which agrees well with the predictions from AGB evolution modelling. While further observational evidences are required to definitively fix the main actors in the pollution of the interstellar medium from which new generation of stars formed in GCs, the present results confirm that the gas ejected by stars of mass in the range $4~\rm M_{\odot} \leq \rm M \leq 8~\rm M_{\odot}$ during the AGB phase share the same chemical patterns traced by stars in GCs.


title: Maximum-Likelihood Power-Distortion Monitoring for GNSS Signal Authentication

abstract: We propose an extension to the so-called PD detector. The PD detector jointly monitors received power and correlation profile distortion to detect the presence of GNSS carry-off-type spoofing, jamming, or multipath. We show that classification performance can be significantly improved by replacing the PD detector's symmetric-difference-based distortion measurement with one based on the post-fit residuals of the maximum-likelihood estimate of a single-signal correlation function model. We call the improved technique the PD-ML detector. In direct comparison with the PD detector, the PD-ML detector exhibits improved classification accuracy when tested against an extensive library of recorded field data. In particular, it is (1) significantly more accurate at distinguishing a spoofing attack from a jamming attack, (2) better at distinguishing multipath-afflicted data from interference-free data, and (3) less likely to issue a false alarm by classifying multipath as spoofing. The PD-ML detector achieves this improved performance at the expense of additional computational complexity.


title: Incoherent diffractive photoproduction of $J/\psi$ and $\Upsilon$ on heavy nuclei in the color dipole approach

abstract: We calculate cross sections and transverse momentum distributions for the incoherent diffractive production of vector mesons $J/\psi$ and $\Upsilon$ on heavy nuclei. In distinction to coherent diffraction, the nucleus is allowed to break up, but except for the vector meson no new particles are produced in the reaction. Within the color dipole approach, we derive the multiple scattering expansion of the incoherent diffractive cross section as an expansion over quasielastic scatterings of the color dipole. We also compare our results to the measurement of the ALICE collaboration for incoherent $J/\psi$ production at $\sqrt{s_{\rm NN}} = 2.76 \, \rm {TeV}$ and show predictions at $\sqrt{s_{\rm NN}} = 5.02 \, \rm{TeV}$. We also briefly discuss a possible contribution to $J/\psi$ production in peripheral collisions in the $70 \div 90 \%$ centrality class.


title: The dp-rank of abelian groups

abstract: An equation to compute the dp-rank of any abelian group is given. It is also shown that its dp-rank, or more generally that of any one-based group, agrees with its Vapnik-Chervonenkis density. Furthermore, strong abelian groups are characterised to be precisely those abelian groups $A$ such that there is only finitely many primes $p$ such that the group $A/pA$ is infinite and for every prime $p$, there is only finitely many natural numbers $n$ such that $(p^nA)[p]/(p^{n+1}A)[p]$ is infinite. Finally, it is shown that an infinite stable field of finite dp-rank is algebraically closed.


title: Integral representation of channel flow with interacting particles

abstract: We construct a boundary integral representation for the low-Reynolds-number flow in a channel in the presence of freely-suspended particles (or droplets) of arbitrary size and shape. We demonstrate that lubrication theory holds away from the particles at horizontal distances exceeding the channel height and derive a multipole expansion of the flow which is dipolar to the leading approximation. We show that the dipole moment of an arbitrary particle is a weighted integral of the stress and the flow at the particle surface, which can be determined numerically. We introduce the equation of motion that describes hydrodynamic interactions between arbitrary, possibly different, distant particles, with interactions determined by the product of the mobility matrix and the dipole moment. Further, the problem of three identical interacting spheres initially aligned in the streamwise direction is considered and the experimentally observed "pair exchange" phenomenon is derived analytically and confirmed numerically. For non-aligned particles, we demonstrate the formation of a configuration with one particle separating from a stable pair. Our results suggest that in a dilute initially homogeneous particulate suspension flowing in a channel the particles will eventually separate into singlets and pairs.


title: Molecular outflow and feedback in the obscured Quasar XID2028 revealed by ALMA

abstract: We imaged with ALMA and ARGOS/LUCI the molecular gas and the dust and stellar continuum in XID2028, an obscured QSO at z=1.593, where the presence of a massive outflow in the ionized gas component traced by the [O III]5007 emission has been resolved up to 10 kpc. This target represents a unique test case to study QSO 'feedback in action' at the peak epoch of AGN-galaxy coevolution. The QSO has been detected in the CO(5-4) transition and in the 1.3mm continuum, at ~30 and ~20 {\sigma} significance respectively, with both emissions confined in the central (<4 kpc) radius area. Our analysis suggests the presence of a fast rotating molecular disc (v~400 km/s) on very compact scales, and well inside the galaxy extent seen in the rest-frame optical light (~10 kpc, as inferred from the LUCI data). Adding available measurements in additional two CO transitions, CO(2-1) and CO(3-2), we could derive a total gas mass of ~10$^{10}$ M$_\odot$, thanks to a critical assessment of CO excitation and the comparison with Rayleigh-Jeans continuum estimate. This translates into a very low gas fraction (<5%) and depletion time scales of 40-75 Myr, reinforcing the result of atypical gas consumption conditions in XID2028, possibly due to feedback effects on the host galaxy. Finally, we also detect at ~5{\sigma} the presence of high velocity CO gas, which we interpret as a signature of galaxy-scale molecular outflow, spatially coincident with the ionised gas outflow. XID2028 represents therefore a unique case where the measurement of total outflowing mass (~500-800 M$_\odot$/yr) including the molecular and atomic components, in both the ionised and neutral phases, has been attempted for a high-z QSO.


title: On Combinatorial Types of Periodic Orbits of the Map $x \mapsto kx$ (mod $\mathbb Z$)

abstract: We study the combinatorial types of periodic orbits of the standard covering endomorphisms ${\mathbf m}_k(x)=k x \ (\text{mod} \ {\mathbb Z})$ of the circle for integers $k \geq 2$ and the frequency with which they occur. For any $q$-cycle $\sigma$ in the permutation group ${\mathcal S}_q$, we give a full description of the set of period $q$ orbits of ${\mathbf m}_k$ that realize $\sigma$ and in particular count how many such orbits there are. The description is based on an invariant called the "fixed point distribution" vector and is achieved by reducing the realization problem to finding the stationary state of an associated Markov chain. Our results generalize earlier work on the special case where $\sigma$ is a rotation cycle, and can be viewed as a missing combinatorial ingredient for a proper understanding of the dynamics of complex polynomial maps of degree $\geq 3$ and the structure of their parameter spaces.


title: de Sitter geodesics in stereographic charts

abstract: It is shown that the stereographic coordinates on the $(1+3)$-dimensional de Sitter spacetime induce, in the parametric equations of all the timelike geodesics, artificial singularities which do not have any physical motivation.


title: Statistical Mechanics Basis of Macleod Formula

abstract: In this work the theoretical basis for the famous formula of Macleod, relating the surface tension of a liquid in equilibrium with its own vapor to the one-particle densities in the two phases of the system, is derived. Using the statistical- mechanical definition of the surface tension, it is proved that this property is, at the first approximation, given by the Macleod formula.


title: Camera Calibration for Daylight Specular-Point Locus

abstract: In this paper we present a new camera calibration method aimed at finding a straight-line locus, in a special colour feature space, that is traversed by daylights and as well also approximately followed by specular points. The aim of the calibration is to enable recovering the colour of the illuminant in a scene, using the calibrated camera. First we prove theoretically that any candidate specular points, for an image that is generated by a specific camera and taken under a daylight, must lie on a straight line in log-chromaticity space, for a chromaticity that is generated using a geometric-mean denominator. Use is made of the assumptions that daylight illuminants can be approximated using Planckians and that camera sensors are narrowband or can be made so by spectral sharpening. Then we show how a particular camera can be calibrated so as to discover this locus. As applications we use this curve for illuminant detection, and also for re-lighting of images to show they would appear under lighting having a different colour temperature.


title: A new holographic dark energy model in Brans-Dicke theory with logarithmic scalar field

abstract: We study a holographic dark energy model in the framework of Brans-Dicke (BD) theory with taking into account the interaction between dark matter and holographic dark energy. We use the recent observational data sets, namely SN Ia compressed Joint Light-Analysis(cJLA) compilation, Baryon Acoustic Oscillations (BAO) from BOSS DR12 and the Cosmic Microwave Background (CMB) of Planck 2015. After calculating the evolution of the equation of state as well as the deceleration parameters, we find that with a logarithmic form for the BD scalar field the phantom crossing can be achieved in the late time of cosmic evolution. Unlike the conventional theory of holographic dark energy in standard cosmology ($\omega_D=0$), our model results a late time accelerated expansion. It is also shown that the cosmic coincidence problem may be resolved in the proposed model. We execute the statefinder and Om diagnostic tools and demonstrate that interaction term does not play a significant role. Based on the observational data sets used in this paper it seems that the best value with $1\sigma$ and $2\sigma$ confidence interval are $\Omega_m=0.268^{+0.008~+0.010}_{-0.007~-0.009}$, $ \alpha=3.361^{+0.332~+0.483}_{-0.401~-0.522}$, $\beta=5.560^{+0.541~+0.780}_{-0.510~-0.729}$, $c=0.777^{+0.023~+0.029}_{-0.017~-0.023}$ and $b^2 =0.045$, according to which we find that the proposed model in the presence of interaction is compatible with the recent observational data.


title: X-Cube Fracton Model on Generic Lattices: Phases and Geometric Order

abstract: Fracton order is a new kind of quantum order characterized by topological excitations that exhibit remarkable mobility restrictions and a robust ground state degeneracy (GSD) which can increase exponentially with system size. In this paper, we present a generic lattice construction (in three dimensions) for a generalized X-cube model of fracton order, where the mobility restrictions of the subdimensional particles inherit the geometry of the lattice. This helps explain a previous result that lattice curvature can produce a robust GSD, even on a manifold with trivial topology. We provide explicit examples to show that the (zero temperature) phase of matter is sensitive to the lattice geometry. In one example, the lattice geometry confines the dimension-1 particles to small loops, which allows the fractons to be fully mobile charges, and the resulting phase is equivalent to (3+1)-dimensional toric code. However, the phase is sensitive to more than just lattice curvature; different lattices without curvature (e.g. cubic or stacked kagome lattices) also result in different phases of matter, which are separated by phase transitions. Unintuitively however, according to a previous definition of phase [Chen, Gu, Wen 2010], even just a rotated or rescaled cubic lattice results in different phases of matter, which motivates us to propose a new and coarser definition of phase for gapped ground states and fracton order. The new equivalence relation between ground states is given by the composition of a local unitary transformation and a quasi-isometry (which can rotate and rescale the lattice); equivalently, ground states are in the same phase if they can be adiabatically connected by varying both the Hamiltonian and the positions of the degrees of freedom (via a quasi-isometry). In light of the importance of geometry, we further propose that fracton orders should be regarded as a geometric order.


title: COLOSSUS: A python toolkit for cosmology, large-scale structure, and dark matter halos

abstract: This paper introduces Colossus, a public, open-source python package for calculations related to cosmology, the large-scale structure (LSS) of matter in the universe, and the properties of dark matter halos. The code is designed to be fast and easy to use, with a coherent, well-documented user interface. The cosmology module implements Friedman-Lemaitre-Robertson-Walker cosmologies including curvature, relativistic species, and different dark energy equations of state, and provides fast computations of the linear matter power spectrum, variance, and correlation function. The LSS module is concerned with the properties of peaks in Gaussian random fields and halos in a statistical sense, including their peak height, peak curvature, halo bias, and mass function. The halo module deals with spherical overdensity radii and masses, density profiles, concentration, and the splashback radius. To facilitate the rapid exploration of these quantities, Colossus implements more than 40 different fitting functions from the literature. I discuss the core routines in detail, with particular emphasis on their accuracy. Colossus is available at bitbucket.org/bdiemer/colossus.


title: Berry phase of adiabatic electronic configurations in fullerene anions

abstract: The selection rule on vibronic angular momentum of $t_{1u}^n \otimes h_g$ Jahn-Teller problem ($n = $ 1-5) is reinvestigated. It is shown that among three adiabatic orbitals only two have nonzero Berry phase. Thus, the Berry phase of adiabatic electronic configurations depends on the spin multiplicity as well as the number of electrons. On this basis, the general relation between the Berry phase and the angular momentum is described. It allows, in particular, to clarify the nature of vibronic states arising from high spin configurations. In comparison with the previous solution for the low-lying vibronic states for bimodal systems, the present solutions correctly fulfill all the symmetry requirement.


title: Large Higgs-electron Yukawa coupling in 2HDM

abstract: The present upper bound on $\kappa_e$, the ratio between the electron Yukawa coupling and its Standard Model value, is of ${\cal O}(600)$. We ask what would be the implications in case that $\kappa_e$ is close to this upper bound. The simplest extension that allows for such enhancement is that of two Higgs doublet models (2HDM) without natural flavor conservation. In this framework, we find the following consequences: (i) Under certain conditions, measuring $\kappa_e$ and $\kappa_V$ would be enough to predict values of Yukawa couplings for other fermions and for the $H$ and $A$ scalars. (ii) In the case that the scalar potential has a softly broken $Z_2$ symmetry, the second Higgs doublet must be light, but if there is hard breaking of the symmetry, the second Higgs doublet can be much heavier than the electroweak scale and still allow the electron Yukawa coupling to be very different from its SM value. (iii) CP must not be violated at a level higher than ${\cal O}(0.01/\kappa_e)$ in both the scalar potential and the Yukawa sector. (iv) LHC searches for $e^+e^-$ resonances constrain this scenario in a significant way. Finally, we study the implications for models where one of the scalar doublets couples only to the first generation, or only to the third generation.


title: New Distances to Four Supernova Remnants

abstract: Distances are found for four supernova remnants without previous distance measurements. H I spectra and H I channel maps are used to determine the maximum velocity of H I absorption for the four supernova remnants (SNRs). We examined $^{13}$CO emission spectra and channel maps to look for possible molecular gas associated with each SNR, but did not find any. The resulting distances for the SNRs are $3.5 \pm 0.2$ kpc (G24.7+0.6), $4.7 \pm 0.3$ kpc (G29.6+0.1), $4.1 \pm 0.5$ kpc (G41.5+0.4) and $4.5 \pm 0 .4 - 9.0 \pm 0.4$ kpc (G57.2+0.8).


title: Tuning the Molecular Weight Distribution from Atom Transfer Radical Polymerization Using Deep Reinforcement Learning

abstract: We devise a novel technique to control the shape of polymer molecular weight distributions (MWDs) in atom transfer radical polymerization (ATRP). This technique makes use of recent advances in both simulation-based, model-free reinforcement learning (RL) and the numerical simulation of ATRP. A simulation of ATRP is built that allows an RL controller to add chemical reagents throughout the course of the reaction. The RL controller incorporates fully-connected and convolutional neural network architectures and bases its decision upon the current status of the ATRP reaction. The initial, untrained, controller leads to ending MWDs with large variability, allowing the RL algorithm to explore a large search space. When trained using an actor-critic algorithm, the RL controller is able to discover and optimize control policies that lead to a variety of target MWDs. The target MWDs include Gaussians of various width, and more diverse shapes such as bimodal distributions. The learned control policies are robust and transfer to similar but not identical ATRP reaction settings, even under the presence of simulated noise. We believe this work is a proof-of-concept for employing modern artificial intelligence techniques in the synthesis of new functional polymer materials.


title: Two-shot measurement of spatial coherence

abstract: We propose and demonstrate an interferometric scheme for measuring the two-dimensional two-point cross-spectral density function in a two-shot manner. Our scheme comprises a Michelson interferometer with a converging lens in one of the arms of the interferometer, and the crosss-pectral density function of an input optical field gets encoded in the intensity distribution of the output interferograms. This scheme works for any cross-spectral density function that is real and that depends on the spatial coordinates only through their difference. Using this scheme, we report measurements of several lab-synthesized cross-spectral density functions with very good agreement with theory. Our measurement technique can be very important for applications that are based on utilizing the partial spatial coherence properties of optical fields.


title: Search for pulsars and transients with the GMRT

abstract: This paper details on the discovery of 21 pulsars using the Giant Metrewave Radio Telescope (GMRT) from targeted (Fermi directed search) and blind surveys (GMRT High Resolution Southern Sky - GHRSS) and results from the follow up studies. We discovered seven millisecond pulsars (MSPs) in the Fermi directed searches, which are the first Galactic MSPs discovered with the GMRT. We have discovered 13 pulsars (including a MSP and two mildly recycled pulsars) with the GHRSS survey, which is an off-Galactic-plane survey at 322 MHz with complementary target sky (declination range -40 deg to -54 deg) to other ongoing low-frequency surveys by GBT and LOFAR. The simultaneous time-domain and imaging study for localising pulsars and transients and efficient candidate investigation with machine learning are some of the features of the GHRSS survey, which are also finding application in the SKA design methodology.


title: Deciphering the 'fuzzy' interaction of FG nucleoporins and transport factors using SANS

abstract: The largely intrinsically disordered phenylalanine-glycine-rich nucleoporins (FG Nups) underline a selectivity mechanism, which enables the rapid translocation of transport factors (TFs) through the nuclear pore complexes (NPCs). Conflicting models of NPC transport have assumed that FG Nups undergo different conformational transitions upon interacting with TFs. To selectively characterize conformational changes in FG Nups induced by TFs we performed small-angle neutron scattering (SANS) with contrast matching. Conformational ensembles derived SANS data indicate an increase in the overall size of FG Nups is associated with TF interaction. Moreover, the organization of the FG motif in the interacting state is consistent with prior experimental analyses defining that FG motifs undergo conformational restriction upon interacting with TFs. These results provide structural insights into a highly dynamic interaction and illustrate how functional disorder imparts rapid and selective FG Nup / TF interactions.


title: Self-acceleration in scalar-bimetric theories

abstract: We describe scalar-bimetric theories where the dynamics of the Universe are governed by two separate metrics, each with an Einstein-Hilbert term. In this setting, the baryonic and dark matter components of the Universe couple to metrics which are constructed as functions of these two gravitational metrics. The scalar field, contrary to dark energy models, does not have a potential whose role is to mimic a late-time cosmological constant. The late-time acceleration of the expansion of the Universe can be easily obtained at the background level in these models by appropriately choosing the coupling functions appearing in the decomposition of the vierbeins for the baryonic and dark matter metrics. We explicitly show how the concordance model can be retrieved with negligible scalar kinetic energy. This requires the scalar coupling functions to show variations of order unity during the accelerated expansion era. This leads in turn to deviations of order unity for the effective Newton constants and a fifth force that is of the same order as Newtonian gravity, with peculiar features. The baryonic and dark matter self-gravities are amplified although the gravitational force between baryons and dark matter is reduced and even becomes repulsive at low redshift. This slows down the growth of baryonic density perturbations on cosmological scales, while dark matter perturbations are enhanced. In our local environment, the upper bound on the time evolution of Newton's constant requires an efficient screening mechanism that both damps the fifth force on small scales and decouples the local value of Newton constant from its cosmological value. This cannot be achieved by a quasi-static chameleon mechanism, and requires going beyond the quasi-static regime and probably using derivative screenings, such as Kmouflage or Vainshtein screening, on small scales.


title: Shaping Long-lived Electron Wavepackets for Customizable Optical Spectra

abstract: Electrons in atoms and molecules are versatile physical systems covering a vast range of light-matter interactions, enabling the physics of Rydberg states, photon-photon bound states, simulation of condensed matter Hamiltonians, and quantum sources of light. A limitation on the versatility of such electronic systems for optical interactions would appear to arise from the discrete nature of the electronic transitions and from the limited ionization energy, constraining the energy scale through which the rich physics of bound electrons can be accessed. In this work, we propose the concept of shaping spatially confined electron wavepackets as superpositions of extended states in the ionization continuum. These wavepackets enable customizable optical emission spectra transitions in the eV-keV range. We find that the specific shaping lengthens the diffraction lifetime of the wavepackets in exchange for increasing their spatial spreads. Finally, we study the spontaneous radiative transitions of these excited states, examining the influence of the wavepacket shape and size on the radiation rate of the excited state. We observe that the rate of radiative capture is primarily limited by the diffraction lifetime of the wavepacket. The approach proposed in this work could have applications towards developing designer emitters at tunable energy and length scales, potentially bringing the physics of Rydberg states to new energy scales and extending the range and versatility of emitters that can be developed using shaped electrons.


title: Decay of the de Sitter Vacuum

abstract: The decay rate of the Bunch-Davies state of a massive scalar field in the expanding flat spatial sections of de Sitter space is determined by an analysis of the particle pair creation process in real time. The Feynman definition of particle and antiparticle Fourier mode solutions of the scalar wave equation, and their adiabatic phase analytically continued to the complexified time domain, show conclusively that the Bunch-Davies state is not the vacuum state at late times. The closely analogous creation of charged particle pairs in a uniform electric field is reviewed and Schwinger's result for the vacuum decay rate is recovered by the real time analysis. The vacuum decay rate in each case is also calculated by switching the background field on adiabatically, allowing it to act for a very long time, and then adiabatically switching it off again. In both the uniform electric field and de Sitter cases the particles created while the field is switched on are verified to be real, in the sense that they persist in the final asymptotic flat zero-field region. In the de Sitter case there is an interesting residual dependence of the rate on how the de Sitter phase is ended, indicating a greater sensitivity to spatial boundary conditions. The electric current of the created particles in the E-field case and their energy density and pressure in the de Sitter case are also computed, and the magnitude of their backreaction effects on the background field estimated. Possible consequences of the Hubble scale instability of the de Sitter vacuum for cosmology, vacuum dark energy, and the cosmological `constant' problem are discussed.


title: Continuous Optimization of Adaptive Quadtree Structures

abstract: We present a novel continuous optimization method to the discrete problem of quadtree optimization. The optimization aims at achieving a quadtree structure with the highest mechanical stiffness, where the edges in the quadtree are interpreted as structural elements carrying mechanical loads. We formulate quadtree optimization as a continuous material distribution problem. The discrete design variables (i.e., to refine or not to refine) are replaced by continuous variables on multiple levels in the quadtree hierarchy. In discrete quadtree optimization, a cell is only eligible for refinement if its parent cell has been refined. We propose a continuous analogue to this dependency for continuous multi-level design variables, and integrate it in the iterative optimization process. Our results show that the continuously optimized quadtree structures perform much stiffer than uniform patterns and the heuristically optimized counterparts. We demonstrate the use of adaptive structures as lightweight infill for 3D printed parts, where uniform geometric patterns have been typically used in practice.


title: Conflict-Free Coloring of String Graphs

abstract: Conflict-free coloring (in short, CF-coloring) of a graph $G = (V,E)$ is a coloring of $V$ such that the neighborhood of each vertex contains a vertex whose color differs from the color of any other vertex in that neighborhood. Bounds on CF-chromatic numbers have been studied both for general graphs and for intersection graphs of geometric shapes. In this paper we obtain such bounds for several classes of string graphs, i.e., intersection graphs of curves in the plane: (i) We provide a general upper bound of $O(\chi(G)^2 \log n)$ on the CF-chromatic number of any string graph $G$ with $n$ vertices in terms of the classical chromatic number $\chi(G)$. This result stands in contrast to general graphs where the CF-chromatic number can be $\Omega(\sqrt{n})$ already for bipartite graphs. (ii) For some central classes of string graphs, the CF-chromatic number is as large as $\Theta(\sqrt{n})$, which is the upper bound for any graph even in the non-geometric context. For several such classes (e.g., intersection graphs of frames) we prove a tight bound of $\Theta(\log n)$ with respect to the notion of $k$-CF-coloring (in which the punctured neighborhood of each vertex contains a color that appears at most $k$ times), for a small constant $k$. (iii) We obtain a general upper bound on the $k$-CF-chromatic number of arbitrary hypergraphs: Any hypergraph with $m$ hyperedges can be $k$-CF colored with $\tilde{O}(m^{\frac{1}{k+1}})$ colors. This bound, which extends a bound of Pach and Tardos (2009), is tight for some string graphs, up to a logarithmic factor. (iv) Our fourth result concerns circle graphs in which coloring problems are motivated by VLSI designs. We prove a tight bound of $\Theta(\log n)$ on the CF-chromatic number of circle graphs, and an upper bound of $O(\log^{3} n)$ for a wider class that contains circle graphs, namely, intersection graphs of grounded L-shapes.


title: Scalable Graphene Aptasensors for Drug Quantification

abstract: Simpler and more rapid approaches for therapeutic drug-level monitoring are highly desirable to enable use at the point-of-care. We have developed an all-electronic approach for detection of the HIV drug tenofovir based on scalable fabrication of arrays of graphene field-effect transistors (GFETs) functionalized with a commercially available DNA aptamer. The shift in the Dirac voltage of the GFETs varied systematically with the concentration of tenofovir in deionized water, with a detection limit less than 1 ng/mL. Tests against a set of negative controls confirmed the specificity of the sensor response. This approach offers the potential for further development into a rapid and convenient point-of-care tool with clinically relevant performance.


title: The role of strangeness and isospin in low density expansions of hadronic matter

abstract: We compare relativistic mean field models with their low density expansion counterparts used to mimic non-relativistic models by consistently expanding the baryonic scalar density in powers of the baryonic number density up to ${\cal O}(13/3)$, which goes two orders beyond the order considered in previous works. We show that, due to the non-trivial density dependence of the Dirac mass, the convergence of the expansion is very slow, and the validity of the non-relativistic approximation {\bf is} questionable even at subsaturation densities. In order to analyze the roles played by strangeness and isospin we consider $n-\Lambda$ and $n-p$ matter separately. Our results indicate that these degrees of freedom play quite different roles in the expansion mechanism and $n-\Lambda$ matter can be better described by low density expansions than $n-p$ matter in general.


title: The Matrix Element Method at next-to-leading order QCD for hadronic collisions: Single top-quark production at the LHC as an example application

abstract: In a recent work the authors have presented a general algorithm to extend the Matrix Element Method (MEM) to the hadronic production of coloured partons taking into account next-to-leading-order (NLO) corrections in quantum chromodynamics (QCD). In this article, the general algorithm is applied to the production of single top quarks at the LHC. In particular, the generation of unweighted events following the NLO predictions is presented. Treating these events as the result of a toy experiment we illustrate the first application of the Matrix Element Method at NLO QCD for hadronic jet production. As a concrete example, we study the determination of the top-quark mass. We show that the inclusion of the NLO corrections can lead to sizeable effects compared to the Matrix Element Method relying on leading-order predictions only and that the incorporation of the NLO corrections is mandatory to obtain reliable estimates of the theoretical uncertainties. In addition, we find that measuring the top-quark mass using the MEM in single top-quark production offers an interesting alternative to mass measurements in top-quark pair production.


title: The constraint equations of Lovelock gravity theories: a new $\sigma_k$-Yamabe problem

abstract: This paper is devoted to the study of the constraint equations of the Lovelock gravity theories. In the case of an empty, compact, conformally flat, time-symmetric, and space-like manifold, we show that the Hamiltonian constraint equation becomes a generalisation of the $\sigma_k$-Yamabe problem. That is to say, the prescription of a linear combination of the $\sigma_k$-curvatures of the manifold. We search solutions in a conformal class for a compact manifold. Using the existing results on the $\sigma_k$-Yamabe problem, we describe some cases in which they can be extended to this new problem. This requires to study the concavity of some polynomial. We do it in two ways: regarding the concavity of an entire root of this polynomial, which is connected to algebraic properties of the polynomial; and seeking analytically a concavifying function. This gives several cases in which a conformal solution exists. At last we show an implicit function theorem in the case of a manifold with negative scalar curvature, and find a conformal solution when the Lovelock theories are close to General Relativity.


title: Controlling light emission with electron wave interference

abstract: It is a long standing question whether or not one can change the nature of spontaneous emission by a free electron through shaping the electron wavefunction. On one hand, shaping the electron wavefunction changes the respective charge and current densities of the electron. On the other hand, spontaneous emission of an electron is an incoherent process and can often be insensitive to the shape of the electron wavefunction. In this work, we arrive at an affirmative answer examining Bremsstrahlung radiation by free electron superposition states. We find that the radiation can be markedly different from an incoherent sum of the radiations of the two states because of interference of the radiation amplitudes from the two components of the superposition. The ability to control free electron spontaneous emission via interference may eventually result in a new degree of control over radiation over the entire electromagnetic spectrum in addition to the ability to deterministically introduce quantum behavior into normally classical light emission processes.


title: Precision Measurement of the Hadronic Contribution to the Muon Anomalous Magnetic Moment

abstract: We report on a precision measurement of the cross section for the reaction $e^+e^-\to\pi^+\pi^-$ in the mass range $0.30<M_{\pi\pi}<1.00$ GeV with the initial state radiation (ISR) method, using 817 pb$^{-1}$ of data at $e^+e^-$ center-of-mass energies near 3.77 GeV and 586 pb$^{-1}$ of data at $e^+e^-$ center-of-mass energies near 4.17 GeV, collected with the CLEO-c detector at the CESR $e^+e^-$ collider at Cornell University. The integrated cross sections in the range $0.30<M_{\pi\pi}<1.00$ GeV for the process $e^+e^-\to\pi^+\pi^-$ are determined with a statistical uncertainty of $0.7\%$ and a systematic uncertainty of $1.5\%$. The leading-order hadronic contribution to the muon anomalous magnetic moment calculated using these measured $e^+e^-\to\pi^+\pi^-$ cross sections in the range $M_{\pi\pi}=0.30$ to 1.00 GeV is calculated to be $(500.4\pm3.6 (\mathrm{stat})\pm 7.5(\mathrm{syst}))\times10^{-10}$.


title: Photodissociation dynamics in the first absorption band of pyrrole: II. Photofragment distributions for the $^1\!A_2(\pi \sigma^*) \leftarrow \tilde{X}^1\!A_1(\pi\pi)$ transition

abstract: The analysis of the total kinetic energy release (TKER) of the photofragments pyrrolyl + H-atom formed in the photodissociation of pyrrole in the low-lying state $^1\!A_2(\pi \sigma^*)$ is presented. The TKER distributions contain complementary and often more precise information on the fragmentation process than the broad diffuse absorption spectra. The distributions are calculated quantum mechanically for the diabatic state $^1\!A_2(\pi \sigma^*)$ either isolated or coupled to the ground electronic state at an exit channel conical intersection. The calculations use the novel ab initio quasi-diabatic potential energy matrix constructed in paper I. The approximate overlap integral-based adiabatic mapping approach is introduced with which the quantum mechanical TKER distributions can be efficiently and accurately reproduced. Finally, the calculated TKERs are compared with the experimental results. The main features of the measured vibrationally resolved distributions are reproduced, and the spectral peaks are assigned and interpreted in detail.


title: Supernova 1572 and other newly observed stars in the literature of the time

abstract: SN 1572, also known as the Tycho's Nova, was a supernova largely reported and discussed in the literature of the time. Here we talk about this literature. In the latest texts, we find also mentioned the Kepler's Nova, today known as SN 1604, and other newly observed stars. We discuss them too.


title: Quantum ergodicity and mixing and their classical limits with quantum kicked rotor

abstract: We study the ergodicity and mixing of quantum kicked rotor (QKR) with two distinct approaches. In one approach, we use the definitions of quantum ergodicity and mixing recently proposed in [Phys. Rev. E 94, 022150 (2016)], which involve only eigen-energies (Floquet quasi-energies for QKR). In the other approach, we study ergodicity and mixing with quantum Poincar\`e section, which is plotted with a method that maps a wave function unitarily onto quantum phase space composed of Planck cells. Classical Poincar\`e section can be recovered with the effective Planck constant gradually diminishing. We demonstrate that the two approaches can capture the quantum and classical characteristics of ergodicity and mixing of QKR, and give consistent results with classical model at semiclassical limit. Therefore, we establish a correspondence between quantum ergodicity (mixing) and classical ergodicity (mixing).


title: Selective light-induced mass transport in amorphous $As_{x}Se_{100-x}$ films driven by the composition tuning: effect of temperature on maximum acceleration

abstract: Under irradiation of amorphous $As_{x}Se_{100-x}$ films by band-gap light, it was observed that lateral mass transport of the film material changed the direction of movement in the region (4 < x < 5 at%.) of topological structural transition (crossover "from light to dark" to "from dark to light"). We propose a model that qualitatively describes the observed phenomenon. In experimental testing of the model by Kelvin probe force microscopy, for the first time to our knowledge, it was detected that with the increase of temperature (from ambient to glass-transition temperature, Tg) photo-induced mass transport in glassy semiconductors occurs with the rates, which significantly exceed those at room temperature.


title: Studying the Ultraviolet Spectrum of the First Spectroscopically Confirmed Supernova at redshift two

abstract: We present observations of DES16C2nm, the first spectroscopically confirmed hydrogen-free superluminous supernova (SLSN-I) at redshift z~2. DES16C2nm was discovered by the Dark Energy Survey (DES) Supernova Program, with follow-up photometric data from the Hubble Space Telescope, Gemini, and the European Southern Observatory Very Large Telescope supplementing the DES data. Spectroscopic observations confirm DES16C2nm to be at z=1.998, and spectroscopically similar to Gaia16apd (a SLSN-I at z=0.102), with a peak absolute magnitude of U=-22.26$\pm$0.06. The high redshift of DES16C2nm provides a unique opportunity to study the ultraviolet (UV) properties of SLSNe-I. Combining DES16C2nm with ten similar events from the literature, we show that there exists a homogeneous class of SLSNe-I in the UV (~2500A), with peak luminosities in the (rest-frame) U band, and increasing absorption to shorter wavelengths. There is no evidence that the mean photometric and spectroscopic properties of SLSNe-I differ between low (z<1) and high redshift (z>1), but there is clear evidence of diversity in the spectrum at <2000A, possibly caused by the variations in temperature between events. No significant correlations are observed between spectral line velocities and photometric luminosity. Using these data, we estimate that SLSNe-I can be discovered to z=3.8 by DES. While SLSNe-I are typically identified from their blue observed colors at low redshift (z<1), we highlight that at z>2 these events appear optically red, peaking in the observer-frame z-band. Such characteristics are critical to identify these objects with future facilities such as the Large Synoptic Survey Telescope, Euclid, and the Wide-Field Infrared Survey Telescope, which should detect such SLSNe-I to z=3.5, 3.7, and 6.6, respectively.

categories: astro-ph.HE astro-ph.CO astro-ph.GA astro-ph.SR


title: Bursting endemic bubbles in an adaptive network

abstract: The spread of an infectious disease is known to change people's behavior, which in turn affects the spread of disease. Adaptive network models that account for both epidemic and behavior change have found oscillations, but in an extremely narrow region of the parameter space, which contrasts with intuition and available data. In this paper we propose a simple SIS epidemic model on an adaptive network with time-delayed rewiring, and show that oscillatory solutions are now present in a wide region of the parameter space. Altering the transmission or rewiring rates reveals the presence of an endemic bubble - an enclosed region of the parameter space where oscillations are observed.


title: Hyperunified field theory and gravitational gauge-geometry duality

abstract: A hyperunified field theory is built in detail based on the postulates of gauge invariance and coordinate independence along with the conformal scaling symmetry. All elementary particles are merged into a single hyper-spinor field and all basic forces are unified into a fundamental interaction governed by the hyper-spin gauge symmetry SP(1,$D_h$-1). The dimension $D_h$ of hyper-spacetime is conjectured to have a physical origin in correlation with the hyper-spin charge of elementary particles. The hyper-gravifield fiber bundle structure of biframe hyper-spacetime appears naturally with the globally flat Minkowski hyper-spacetime as a base spacetime and the locally flat hyper-gravifield spacetime as a fiber that is viewed as a dynamically emerged hyper-spacetime characterized by a non-commutative geometry. The gravitational origin of gauge symmetry is revealed with the hyper-gravifield that plays an essential role as a Goldstone-like field. The gauge-gravity and gravity-geometry correspondences bring about the gravitational gauge-geometry duality. The basic properties of hyperunified field theory and the issue on the fundamental scale are analyzed within the framework of quantum field theory, which allows us to describe the laws of nature in deriving the gauge gravitational equation with the conserved current and the geometric gravitational equations of Einstein-like type and beyond.


title: Forecasting Extreme Events in the Complex Dynamics of a Semiconductor Laser with Feedback

abstract: Complex systems performing spiking dynamics are widespread in Nature. They cover from earthquakes, to neurons, variable stars, social networks, or stock markets. Understanding and characterizing their dynamics is relevant in order to detect transitions, or to predict unwanted extreme events. Here we study the output intensity of a semiconductor laser with feedback, in a regime where it develops a complex spiking behavior, under an ordinal patterns analysis. We unveil that the complex dynamics presents two competing behaviors that can be distinguished with a thresholding method, and we use temporal correlations to forecast the extreme events, and transitions between dynamics.


title: XMM-Newton Observations of the Southeastern Radio Relic in Abell 3667

abstract: Radio relics, elongated, non-thermal, structures located at the edges of galaxy clusters, are the result of synchrotron radiation from cosmic-ray electrons accelerated by merger-driven shocks at the cluster outskirts. However, X-ray observations of such shocks in some clusters suggest that they are too weak to efficiently accelerate electrons via diffusive shock acceleration to energies required to produce the observed radio power. We examine this issue in the merging galaxy cluster Abell 3667 (A3667), which hosts a pair of radio relics. While the Northwest relic in A3667 has been well studied in the radio and X-ray by multiple instruments, the Southeast relic region has only been observed so far by Suzaku, which detected a temperature jump across the relic, suggesting the presence of a weak shock. We present observations of the Southeastern region of A3667 with XMM-Newton centered on the radio relic. We confirm the existence of an X-ray shock with Mach number of about 1.8 from a clear detection of temperature jump and a tentative detection of a density jump, consistent with previous measurements by Suzaku. We discuss the implications of this measurement for diffusive shock acceleration as the main mechanism for explaining the origin of radio relics. We then speculate on the plausibility of alternative scenarios, including re-acceleration and variations in the Mach number along shock fronts.


title: Herschel and ALMA Observations of Massive SZE-selected Clusters

abstract: We present new Herschel observations of four massive, Sunyaev-Zel'dovich Effect (SZE)-selected clusters at $0.3 \leq z \leq 1.1$, two of which have also been observed with ALMA. We detect 19 Herschel/PACS counterparts to spectroscopically confirmed cluster members, five of which have redshifts determined via CO($4-3$) and [CI](${}^3P_1 - {}^3P_0$) lines. The mean [CI]/CO line ratio is $0.19 \pm 0.07$ in brightness temperature units, consistent with previous results for field samples. We do not detect significant stacked ALMA dust continuum or spectral line emission, implying upper limits on mean interstellar medium (H$_2$ + HI) and molecular gas masses. An apparent anticorrelation of $L_{IR}$ with clustercentric radius is driven by the tight relation between star formation rate and stellar mass. We find average specific star formation rate log(sSFR/yr$^{-1}$) = -10.36, which is below the SFR$-M_*$ correlation measured for field galaxies at similar redshifts. The fraction of infrared-bright galaxies (IRBGs; $\log (L_{IR}/L_\odot) > 10.6$) per cluster and average sSFR rise significantly with redshift. For CO detections, we find $f_{gas} \sim 0.2$, comparable to those of field galaxies, and gas depletion timescales of about 2 Gyr. We use radio observations to distinguish active galactic nuclei (AGNs) from star-forming galaxies. At least four of our 19 Herschel cluster members have $q_{IR} < 1.8$, implying an AGN fraction $f_{AGN} \gtrsim 0.2$ for our PACS-selected sample.


title: Analysis and Optimization of Aperture Design in Computational Imaging

abstract: There is growing interest in the use of coded aperture imaging systems for a variety of applications. Using an analysis framework based on mutual information, we examine the fundamental limits of such systems---and the associated optimum aperture coding---under simple but meaningful propagation and sensor models. Among other results, we show that when thermal noise dominates, spectrally-flat masks, which have 50% transmissivity, are optimal, but that when shot noise dominates, randomly generated masks with lower transmissivity offer greater performance. We also provide comparisons to classical pinhole cameras.


title: Learning Sparse Graphs for Prediction and Filtering of Multivariate Data Processes

abstract: We address the problem of prediction of multivariate data process using an underlying graph model. We develop a method that learns a sparse partial correlation graph in a tuning-free and computationally efficient manner. Specifically, the graph structure is learned recursively without the need for cross-validation or parameter tuning by building upon a hyperparameter-free framework. Our approach does not require the graph to be undirected and also accommodates varying noise levels across different nodes.Experiments using real-world datasets show that the proposed method offers significant performance gains in prediction, in comparison with the graphs frequently associated with these datasets.


title: A Mathematical Programming Approach for Integrated Multiple Linear Regression Subset Selection and Validation

abstract: Subset selection for multiple linear regression aims to construct a regression model that minimizes errors by selecting a small number of explanatory variables. Once a model is built, various statistical tests and diagnostics are conducted to validate the model and to determine whether the regression assumptions are met. Most traditional approaches require human decisions at this step. For example, the user adding or removing a variable until a satisfactory model is obtained. However, this trial-and-error strategy cannot guarantee that a subset that minimizes the errors while satisfying all regression assumptions will be found. In this paper, we propose a fully automated model building procedure for multiple linear regression subset selection that integrates model building and validation based on mathematical programming. The proposed model minimizes mean squared errors while ensuring that the majority of the important regression assumptions are met. We also propose an efficient constraint to approximate the constraint for the coefficient t-test. When no subset satisfies all of the considered regression assumptions, our model provides an alternative subset that satisfies most of these assumptions. Computational results show that our model yields better solutions (i.e., satisfying more regression assumptions) compared to the state-of-the-art benchmark models while maintaining similar explanatory power.


title: Hierarchical Bloom Filter Trees for Approximate Matching

abstract: Bytewise approximate matching algorithms have in recent years shown significant promise in de- tecting files that are similar at the byte level. This is very useful for digital forensic investigators, who are regularly faced with the problem of searching through a seized device for pertinent data. A common scenario is where an investigator is in possession of a collection of "known-illegal" files (e.g. a collection of child abuse material) and wishes to find whether copies of these are stored on the seized device. Approximate matching addresses shortcomings in traditional hashing, which can only find identical files, by also being able to deal with cases of merged files, embedded files, partial files, or if a file has been changed in any way. Most approximate matching algorithms work by comparing pairs of files, which is not a scalable approach when faced with large corpora. This paper demonstrates the effectiveness of using a "Hierarchical Bloom Filter Tree" (HBFT) data structure to reduce the running time of collection-against-collection matching, with a specific focus on the MRSH-v2 algorithm. Three experiments are discussed, which explore the effects of different configurations of HBFTs. The proposed approach dramatically reduces the number of pairwise comparisons required, and demonstrates substantial speed gains, while maintaining effectiveness.


title: Macroscopic stability and simplicial norms of hypersurfaces

abstract: We introduce a $\mathbb{Z}$--coefficient version of Guth's macroscopic stability inequality for almost-minimizing hypersurfaces. In manifolds with a lower bound on macroscopic scalar curvature, we use the inequality to prove a lower bound on areas of hypersurfaces in terms of the Gromov simplicial norm of their homology classes. We give examples to show that a very positive lower bound on macroscopic scalar curvature does not necessarily imply an upper bound on the areas of minimizing hypersurfaces.


title: Encoding DNA sequences by integer chaos game representation

abstract: DNA sequences are fundamental for encoding genetic information. The genetic information may not only be understood by symbolic sequences but also from the hidden signals inside the sequences. The symbolic sequences need to be transformed into numerical sequences so the hidden signals can be revealed by signal processing techniques. All current transformation methods encode DNA sequences into numerical values of the same length. These representations have limitations in the applications of genomic signal compression, encryption, and steganography. We propose an integer chaos game representation (iCGR) of DNA sequences and a lossless encoding method DNA sequences by the iCGR. In the iCGR method, a DNA sequence is represented by the iterated function of the nucleotides and their positions in the sequence. Then the DNA sequence can be uniquely encoded and recovered using three integers from iCGR. One integer is the sequence length and the other two integers represent the accumulated distributions of nucleotides in the sequence. The integer encoding scheme can compress a DNA sequence by 2 bits per nucleotide. The integer representation of DNA sequences provides a prospective tool for sequence compression, encryption, and steganography. The Python programs in this study are freely available to the public at https://github.com/cyinbox/iCGR


title: Photoacoustic spectroscopy of NO$_2$ using a mid-infrared pulsed optical parametric oscillator as light source

abstract: A photoacoustic (PA) sensor for spectroscopic measurements of NO$_2$-N$_2$ at ambient pressure and temperature is demonstrated. The PA sensor is pumped resonantly by a nanosecond pulsed single-mode mid-infrared (MIR) optical parametric oscillator (OPO). Spectroscopic measurements of NO$_2$-N$_2$ in the 3.25 $\mu$m to 3.55 $\mu$m wavelength region with a resolution bandwidth of 5 cm$^{-1}$ and with a single shot detection limit of 1.6 ppmV ($\mu$mol/mol) is demonstrated. The measurements were conducted with a constant flow rate of 300 ml/min, thus demonstrating the suitability of the gas sensor for real time trace gas measurements. The acquired spectra is compared with data from the Hitran database and good agreement is found. An Allan deviation analysis shows that the detection limit at optimum integration time for the PAS sensor is 14 ppbV (nmol/mol) at 170 seconds of integration time, corresponding to a normalized noise equivalent absorption (NNEA) coefficient of 3.3$\times 10^{-7}$ W cm$^{-1}$ Hz$^{-1/2}$.


title: Accessiblility Percolation with Crossing Valleys on $n$-ary Trees

abstract: In this paper we study a variation of the accessibility percolation model, this is also motivated by evolutionary biology and evolutionary computation. Consider a tree whose vertices are labeled with random numbers. We study the probability of having a monotone subsequence of a path from the root to a leaf, where any $k$ consecutive vertices in the path contain at least one vertex of the subsequence. An $n$-ary tree, with height $h$, is a tree whose vertices at distance at most $h-1$ to the root have $n$ children. For the case of $n$-ary trees, we prove that, as $h$ tends to infinity the probability of having such subsequence: tends to 1, if $n$ grows significantly faster than $\sqrt[k]{h/(ek)}$; and tends to 0, if $n$ grows significantly slower than $\sqrt[k]{h/(ek)}$.


title: Minors of two-connected graphs of large path-width

abstract: Let $P$ be a graph with a vertex $v$ such that $P\backslash v$ is a forest, and let $Q$ be an outerplanar graph. We prove that there exists a number $p=p(P,Q)$ such that every 2-connected graph of path-width at least $p$ has a minor isomorphic to $P$ or $Q$. This result answers a question of Seymour and implies a conjecture of Marshall and Wood. The proof is based on a new property of tree-decompositions.


title: Material phase classification by means of Support Vector Machines

abstract: The pixel's classification of images obtained from random heterogeneous materials is a relevant step to compute their physical properties, like Effective Transport Coefficients (ETC), during a characterization process as stochastic reconstruction. A bad classification will impact on the computed properties; however, the literature on the topic discusses mainly the correlation functions or the properties formulae, giving little or no attention to the classification; authors mention either the use of a threshold or, in few cases, the use of Otsu's method. This paper presents a classification approach based on Support Vector Machines (SVM) and a comparison with the Otsu's-based approach, based on accuracy and precision. The data used for the SVM training are the key for a better classification; these data are the grayscale value, the magnitude and direction of pixels gradient. For instance, in the case study, the accuracy of the pixel's classification is 77.6% for the SVM method and 40.9% for Otsu's method. Finally, a discussion about the impact on the correlation functions is presented in order to show the benefits of the proposal.


title: BMS Supertranslation Symmetry Implies Faddeev-Kulish Amplitudes

abstract: We show explicitly that, among the scattering amplitudes constructed from eigenstates of the BMS supertranslation charge, the ones that conserve this charge, are equal to those constructed from Faddeev-Kulish states. Thus, Faddeev-Kulish states naturally arise as a consequence of the asymptotic symmetries of perturbative gravity and all charge conserving amplitudes are infrared finite. In the process we show an important feature of the Faddeev-Kulish clouds dressing the external hard particles: these clouds can be moved from the incoming states to the outgoing ones, and vice-versa, without changing the infrared finiteness properties of S matrix elements. We also apply our discussion to the problem of the decoherence of momentum configurations of hard particles due to soft boson effects.


title: From a $(p,2)$-Theorem to a Tight $(p,q)$-Theorem

abstract: A family $F$ of sets is said to satisfy the $(p,q)$-property if among any $p$ sets of $F$ some $q$ intersect. The celebrated $(p,q)$-theorem of Alon and Kleitman asserts that any family of compact convex sets in $\mathbb{R}^d$ that satisfies the $(p,q)$-property for some $q \geq d+1$, can be pierced by a fixed number $f_d(p,q)$ of points. The minimum such piercing number is denoted by $HD_d(p,q)$. Already in 1957, Hadwiger and Debrunner showed that whenever $q>\frac{d-1}{d}p+1$ the piercing number is $HD_d(p,q)=p-q+1$; no exact values of $HD_d(p,q)$ were found ever since. While for an arbitrary family of compact convex sets in $\mathbb{R}^d$, $d \geq 2$, a $(p,2)$-property does not imply a bounded piercing number, such bounds were proved for numerous specific families. The best-studied among them is axis-parallel rectangles in the plane. Wegner and (independently) Dol'nikov used a $(p,2)$-theorem for axis-parallel rectangles to show that $HD_{\mathrm{rect}}(p,q)=p-q+1$ holds for all $q>\sqrt{2p}$. These are the only values of $q$ for which $HD_{\mathrm{rect}}(p,q)$ is known exactly. In this paper we present a general method which allows using a $(p,2)$-theorem as a bootstrapping to obtain a tight $(p,q)$-theorem, for families with Helly number 2, even without assuming that the sets in the family are convex or compact. To demonstrate the strength of this method, we obtain a significant improvement of an over 50 year old result by Wegner and Dol'nikov. Namely, we show that $HD_{\mathrm{d-box}}(p,q)=p-q+1$ holds for all $q > c' \log^{d-1} p$, and in particular, $HD_{\mathrm{rect}}(p,q)=p-q+1$ holds for all $q \geq 7 \log_2 p$ (compared to $q \geq \sqrt{2p}$ of Wegner and Dol'nikov). In addition, for several classes of families, we present improved $(p,2)$-theorems, some of which can be used as a bootstrapping to obtain tight $(p,q)$-theorems.


title: MaxDim of some simple groups

abstract: MaxDim, along with it relatives m and i, are group invariants analogous to the concept of dimension in vector spaces. In this paper we give the following result on the MaxDim of sporadic groups: 1. An inductive argument for the lower bounds of the MaxDim of those groups in the Suzuki tower, utilizing the fact that those groups act on strongly regular graphs. 2. A detailed computation of the MaxDim of Mathieu 24 group, utilizing the fact that it acts as automorphisms of Steiner system. 3. Estimations of the MaxDim for McLaughlin and Higmann-Sims Group.


title: Incommensurate phonon anomaly and the nature of charge density waves in cuprates

abstract: While charge density wave (CDW) instabilities are ubiquitous to superconducting cuprates, the different ordering wavevectors in various cuprate families have hampered a unified description of the CDW formation mechanism. Here we investigate the temperature dependence of the low energy phonons in the canonical CDW ordered cuprate La$_{1.875}$Ba$_{0.125}$CuO$_{4}$. We discover that the phonon softening wavevector associated with CDW correlations becomes temperature dependent in the high-temperature precursor phase and changes from a wavevector of 0.238 reciprocal space units (r.l.u.) below the ordering transition temperature up to 0.3~r.l.u. at 300~K. This high-temperature behavior shows that "214"-type cuprates can host CDW correlations at a similar wavevector to previously reported CDW correlations in non-"214"-type cuprates such as YBa$_{2}$Cu$_{3}$O$_{6+\delta}$. This indicates that cuprate CDWs may arise from the same underlying instability despite their apparently different low temperature ordering wavevectors.


title: Classification vs. Regression in Supervised Learning for Single Channel Speaker Count Estimation

abstract: The task of estimating the maximum number of concurrent speakers from single channel mixtures is important for various audio-based applications, such as blind source separation, speaker diarisation, audio surveillance or auditory scene classification. Building upon powerful machine learning methodology, we develop a Deep Neural Network (DNN) that estimates a speaker count. While DNNs efficiently map input representations to output targets, it remains unclear how to best handle the network output to infer integer source count estimates, as a discrete count estimate can either be tackled as a regression or a classification problem. In this paper, we investigate this important design decision and also address complementary parameter choices such as the input representation. We evaluate a state-of-the-art DNN audio model based on a Bi-directional Long Short-Term Memory network architecture for speaker count estimations. Through experimental evaluations aimed at identifying the best overall strategy for the task and show results for five seconds speech segments in mixtures of up to ten speakers.


title: Tracking Filament Evolution in the Low Solar Corona using Remote-Sensing and In-situ Observations

abstract: In the present work, we analyze a filament eruption associated with an ICME that arrived at L1 on August 5th, 2011. In multi-wavelength SDO/AIA images, three plasma parcels within the filament were tracked at high-cadence along the solar corona. A novel absorption diagnostic technique was applied to the filament material travelling along the three chosen trajectories to compute the column density and temperature evolution in time. Kinematics of the filamentary material were estimated using STEREO/EUVI and STEREO/COR1 observations. The Michigan Ionization Code used inputs of these density, temperature, and speed profiles for the computation of ionization profiles of the filament plasma. Based on these measurements we conclude the core plasma was in near ionization equilibrium, and the ionization states were not frozen-in at the altitudes where they were visible in absorption in AIA images. Additionally, we report that the filament plasma was heterogeneous, and the filamentary material was continuously heated as it expanded in the low solar corona.


title: Long Range Particle Dynamics and the Linear Boltzmann Equation

abstract: This paper gives the first full proof of the justification of the linear Boltzmann equation from an underlying long range particle evolution. We suppose that a tagged particle is interacting with a background via a two body potential that is decaying faster than $ C\exp\left(-C|x|^{\frac{3}{2}}\right) $, and that the background is initially distributed according to a function in $L^1((\mathbb{R}^3,(1+|v|^2)\mathrm{d}v)$ in velocity and uniformly in space. Under finite mass and energy assumptions on the initial density, the tagged particle density converges weak-$\star$ in $L^\infty$ to a solution of the linear Boltzmann equation. The proof uses estimates on two body scattering and on the relationship between long range dynamics and dynamics with a truncated interaction potential to explicitly estimate the error between densities for long and short range dynamics. To compare the difference between the short range dynamics and the linear Boltzmann equation, we use a tree based structure to encode the collisional history of the tagged particle.


title: Thermal radiation and entanglement in proton-proton collisions at the LHC

abstract: The origin of the apparent thermalization in high-energy collisions is investigated using the data of the ATLAS and CMS Collaborations at the LHC. For this purpose, we analyze the transverse momentum distributions in the following proton-proton collision processes, all at $\sqrt{s} = 13$ TeV: i) inclusive inelastic $pp$ collisions; ii) single- and double-diffractive Drell-Yan production $pp \to \mu^+ \mu^- X$; and iii) Higgs boson production. We confirm the relation between the effective temperature and the hard scattering scale observed at lower energies, and find that it extends even to the Higgs boson production process. In addition we find that the thermal component disappears in diffractive events (even though many charged hadrons are still produced). We discuss the implications of our study for the mechanism of multi-particle production -- in particular, we test the hypothesis about the link between quantum entanglement and thermalization in high-energy collisions.


title: Two dimensional collective electron magnetotransport, oscillations and chaos in a semiconductor superlattice

abstract: When quantized, traces of classically chaotic single particle systems include eigenvalue statistics and scars in eigenfuntions. Since 2001, many theoretical and experimental works have argued that classically chaotic single electron dynamics influences and controls collective electron transport. For transport in semiconductor superlattices under tilted magnetic and electric fields, these theories rely on a reduction to a one-dimensional self-consistent drift model. A two-dimensional theory based on self-consistent Boltzmann transport does not support that single electron chaos influences collective transport. This theory agrees with existing experimental evidence of current self-oscillations, predicts spontaneous collective chaos via a period doubling scenario and it could be tested unambiguously by measuring the electric potential inside the superlattice under a tilted magnetic field.


title: IT Risk Assessment for Group6 Healthcare Clinic Report

abstract: Information security and privacy in the healthcare sector is an issue of growing importance. The adoption of digital patient records, increased regulation, provider consolidation and the increasing need for information exchange between patients, providers and payers, all point towards the need for better information security. We critically survey the literature on information security and privacy in healthcare, published in information systems journals as well as many other related disciplines including health informatics, public health, law, medicine, and the trade press and industry reports


title: Scientific Polarization

abstract: Contemporary societies are often "polarized", in the sense that sub-groups within these societies hold stably opposing beliefs, even when there is a fact of the matter. Extant models of polarization do not capture the idea that some beliefs are true and others false. Here we present a model, based on the network epistemology framework of Bala and Goyal ["Learning from neighbors", \textit{Rev. Econ. Stud.} \textbf{65}(3), 784-811 (1998)], in which polarization emerges even though agents gather evidence about their beliefs, and true belief yields a pay-off advantage. The key mechanism that generates polarization involves treating evidence generated by other agents as uncertain when their beliefs are relatively different from one's own.


title: A bandwidth theorem for approximate decompositions

abstract: We provide a degree condition on a regular $n$-vertex graph $G$ which ensures the existence of a near optimal packing of any family $\mathcal H$ of bounded degree $n$-vertex $k$-chromatic separable graphs into $G$. In general, this degree condition is best possible. Here a graph is separable if it has a sublinear separator whose removal results in a set of components of sublinear size. Equivalently, the separability condition can be replaced by that of having small bandwidth. Thus our result can be viewed as a version of the bandwidth theorem of B\"ottcher, Schacht and Taraz in the setting of approximate decompositions. More precisely, let $\delta_k$ be the infimum over all $\delta\ge 1/2$ ensuring an approximate $K_k$-decomposition of any sufficiently large regular $n$-vertex graph $G$ of degree at least $\delta n$. Now suppose that $G$ is an $n$-vertex graph which is close to $r$-regular for some $r \ge (\delta_k+o(1))n$ and suppose that $H_1,\dots,H_t$ is a sequence of bounded degree $n$-vertex $k$-chromatic separable graphs with $\sum_i e(H_i) \le (1-o(1))e(G)$. We show that there is an edge-disjoint packing of $H_1,\dots,H_t$ into $G$. If the $H_i$ are bipartite, then $r\geq (1/2+o(1))n$ is sufficient. In particular, this yields an approximate version of the tree packing conjecture in the setting of regular host graphs $G$ of high degree. Similarly, our result implies approximate versions of the Oberwolfach problem, the Alspach problem and the existence of resolvable designs in the setting of regular host graphs of high degree.


title: Structural Relations of Symmetry among Players in Strategic Games

abstract: The notions of symmetry and anonymity in strategic games have been formalized in different ways in the literature. We propose a combinatorial framework to analyze these notions, using group actions. Then, the same framework is used to define partial symmetries in payoff matrices. With this purpose, we introduce the notion of the role a player plays with respect to another one, and combinatorial relations between roles are studied. Building on them, we define relations directly between players, which provide yet another characterization of structural symmetries in the payoff matrices of strategic games.


title: Approximate Convex Hull of Data Streams

abstract: Given a finite set of points $P \subseteq \mathbb{R}^d$, we would like to find a small subset $S \subseteq P$ such that the convex hull of $S$ approximately contains $P$. More formally, every point in $P$ is within distance $\epsilon$ from the convex hull of $S$. Such a subset $S$ is called an $\epsilon$-hull. Computing an $\epsilon$-hull is an important problem in computational geometry, machine learning, and approximation algorithms. In many real world applications, the set $P$ is too large to fit in memory. We consider the streaming model where the algorithm receives the points of $P$ sequentially and strives to use a minimal amount of memory. Existing streaming algorithms for computing an $\epsilon$-hull require $O(\epsilon^{-(d-1)/2})$ space, which is optimal for a worst-case input. However, this ignores the structure of the data. The minimal size of an $\epsilon$-hull of $P$, which we denote by $\text{OPT}$, can be much smaller. A natural question is whether a streaming algorithm can compute an $\epsilon$-hull using only $O(\text{OPT})$ space. We begin with lower bounds that show that it is not possible to have a single-pass streaming algorithm that computes an $\epsilon$-hull with $O(\text{OPT})$ space. We instead propose three relaxations of the problem for which we can compute $\epsilon$-hulls using space near-linear to the optimal size. Our first algorithm for points in $\mathbb{R}^2$ that arrive in random-order uses $O(\log n\cdot \text{OPT})$ space. Our second algorithm for points in $\mathbb{R}^2$ makes $O(\log(\frac{1}{\epsilon}))$ passes before outputting the $\epsilon$-hull and requires $O(\text{OPT})$ space. Our third algorithm for points in $\mathbb{R}^d$ for any fixed dimension $d$ outputs an $\epsilon$-hull for all but $\delta$-fraction of directions and requires $O(\text{OPT} \cdot \log \text{OPT})$ space.


title: Intra-axonal Diffusivity in Brain White Matter

abstract: Biophysical modeling is the mediator of evaluating the cellular structure of biological tissues using diffusion-weighted MRI. It is however the bottleneck of microstructural MRI. Beyond the complexity of diffusion, the current development is hindered by the fact that biophysical models heavily rely on diffusion-specific properties of diverse cellular compartments that are still unknown and must be measured in vivo. Obtaining such parameters by straightforward fitting is hindered by the degenerated landscape of the likelihood functions, in particular, the signal obtained for multiple diffusion directions and moderate diffusion weighting strength is not enough to estimate these parameters: different parameter constellations explain the signal equally well. The aim of this study is to measure the central parameter of white matter models, namely the intra-axonal water diffusivity in the normal human brain. Proper estimation of this parameter is complicated due to (i) the presence of both intra- and extra-axonal water compartments and (ii) the orientation dispersion of axons. Our measurement involves an efficient suppression of extra-axonal space and all cellular processes oriented outside a narrow cone around the principal fiber direction. This is achieved using a planar water mobility filter -- a strong diffusion weighting that suppresses signal from all molecules that are mobile in the plane transverse to the fiber bundle. Following the planar filter, the diffusivity in the remaining compartment is measured using linear and isotropic weighting. We find the specifically averaged intra-axonal diffusivity $D_0 = 2.25\pm 0.03{\,\rm \mu m^2/ms}$ for the timing of the applied gradients. Extrapolation to the infinite diffusion time gives $D_\infty \approx 2.0{\,\rm \mu m^2/ms}$. This result imposes a strong limitation on the parameter selection for biophysical modeling of diffusion-weighted MRI.


title: Individual Vortex Manipulation and Stick-Slip Motion in Periodic Pinning Arrays

abstract: We numerically examine the manipulation of vortices interacting with a moving trap representing a magnetic force tip translating across a superconducting sample containing a periodic array of pinning sites. As a function of the tip velocity and coupling strength, we find five distinct dynamic phases, including a decoupled regime where the vortices are dragged a short distance within a pinning site, an intermediate coupling regime where vortices in neighboring pinning sites exchange places, an intermediate trapping regime where individual vortices are dragged longer distances and exchange modes of vortices occur in the surrounding pins, an intermittent multiple trapping regime where the trap switches between capturing one or two vortices, and a strong couping regime in which the trap permanently captures and drags two vortices. In some regimes we observe the counterintuitive behavior that slow moving traps couple less strongly to vortices than faster moving traps; however, the fastest moving traps are generally decoupled. The different phases can be characterized by the distances the vortices are displaced and the force fluctuations exerted on the trap. We find different types of stick-slip motion depending on whether vortices are moving into and out of pinning sites, undergoing exchange, or performing correlated motion induced by vortices outside of the trap. Our results are general to the manipulation of other types of particle-based systems interacting with periodic trap arrays, such as colloidal particles or certain types of frictional systems.


title: Practical Bayesian optimization in the presence of outliers

abstract: Inference in the presence of outliers is an important field of research as outliers are ubiquitous and may arise across a variety of problems and domains. Bayesian optimization is method that heavily relies on probabilistic inference. This allows outstanding sample efficiency because the probabilistic machinery provides a memory of the whole optimization process. However, that virtue becomes a disadvantage when the memory is populated with outliers, inducing bias in the estimation. In this paper, we present an empirical evaluation of Bayesian optimization methods in the presence of outliers. The empirical evidence shows that Bayesian optimization with robust regression often produces suboptimal results. We then propose a new algorithm which combines robust regression (a Gaussian process with Student-t likelihood) with outlier diagnostics to classify data points as outliers or inliers. By using an scheduler for the classification of outliers, our method is more efficient and has better convergence over the standard robust regression. Furthermore, we show that even in controlled situations with no expected outliers, our method is able to produce better results.


title: Local existence of strong solutions and weak-strong uniqueness for the compressible Navier-Stokes system on moving domains

abstract: We consider the compressible Navier-Stokes system on time-dependent domains with prescribed motion of the boundary. For both the no-slip boundary conditions as well as slip boundary conditions we prove local-in-time existence of strong solutions. These results are obtained using a transformation of the problem to a fixed domain and an existence theorem for Navier-Stokes like systems with lower order terms and perturbed boundary conditions. We also show the weak-strong uniqueness principle for slip boundary conditions which remained so far open question.


title: Im2Pano3D: Extrapolating 360 Structure and Semantics Beyond the Field of View

abstract: We present Im2Pano3D, a convolutional neural network that generates a dense prediction of 3D structure and a probability distribution of semantic labels for a full 360 panoramic view of an indoor scene when given only a partial observation (<= 50%) in the form of an RGB-D image. To make this possible, Im2Pano3D leverages strong contextual priors learned from large-scale synthetic and real-world indoor scenes. To ease the prediction of 3D structure, we propose to parameterize 3D surfaces with their plane equations and train the model to predict these parameters directly. To provide meaningful training supervision, we use multiple loss functions that consider both pixel level accuracy and global context consistency. Experiments demon- strate that Im2Pano3D is able to predict the semantics and 3D structure of the unobserved scene with more than 56% pixel accuracy and less than 0.52m average distance error, which is significantly better than alternative approaches.


title: Flare Activity of Wide Binary Stars with Kepler

abstract: We present an analysis of flare activity in wide binary stars using a combination of value-added data sets from the NASA Kepler mission. The target list contains a set of previously discovered wide binary star systems identified by proper motions in the Kepler field. We cross-matched these systems with estimates of flare activity for $\sim$200,000 stars in the Kepler field, allowing us to compare relative flare luminosity between stars in coeval binaries. From a sample of 184 previously known wide binaries in the Kepler field, we find 58 with detectable flare activity in at least 1 component, 33 of which are similar in mass (q $>$ 0.8). Of these 33 equal-mass binaries, the majority display similar ($\pm 1$ dex) flare luminosity between both stars, as expected for stars of equal mass and age. However, we find two equal-mass pairs where the secondary (lower mass) star is more active than its counterpart, and two equal-mass pairs where the primary star is more active. The stellar rotation periods are also anomalously fast for stars with elevated flare activity. Pairs with discrepant rotation and activity qualitatively seem to have lower mass ratios. These outliers may be due to tidal spin-up, indicating these wide binaries could be hierarchical triple systems. We additionally present high resolution adaptive optics images for two wide binary systems to test this hypothesis. The demographics of stellar rotation and magnetic activity between stars in wide binaries may be useful indicators for discerning formation scenarios of these systems.


title: Hadron Mass Effects: Kaons at HERMES vs. COMPASS

abstract: Experimental data for integrated kaon multiplicities taken at HERMES and COMPASS measurements look incompatible with each other. In this talk, we investigate the effects of hadron masses calculated at leading-order and leading twist at the kinematics of these two experiments. We present evidence that Hadron Mass Corrections can fully reconcile the data for the $K^+/K^-$ multiplicity ratio, and can also sizeably reduce the apparent large discrepancy in the case of $K^+ + K^-$ data. Residual differences in the shape of the latter one remains to be understood.


title: Quotients of $S^2\times{S^2}$

abstract: We consider closed topological 4-manifolds $M$ with universal cover ${S^2\times{S^2}}$ and Euler characteristic $\chi(M) = 1$. All such manifolds with $\pi=\pi_1(M)\cong {\mathbb Z}/4$ are homotopy equivalent. In this case, we show that there are four homeomorphism types, and propose a candidate for a smooth example which is not homeomorphic to the geometric quotient. If $\pi\cong {\mathbb Z}/2 \times {\mathbb Z}/2$, we show that there are three homotopy types (and between 6 and 24 homeomorphism types).