title: The Dynamics of Norm Change in the Cultural Evolution of Language

abstract: What happens when a new social convention replaces an old one? While the possible forces favoring norm change - such as institutions or committed activists - have been identified since a long time, little is known about how a population adopts a new convention, due to the difficulties of finding representative data. Here we address this issue by looking at changes occurred to 2,541 orthographic and lexical norms in English and Spanish through the analysis of a large corpora of books published between the years 1800 and 2008. We detect three markedly distinct patterns in the data, depending on whether the behavioral change results from the action of a formal institution, an informal authority or a spontaneous process of unregulated evolution. We propose a simple evolutionary model able to capture all the observed behaviors and we show that it reproduces quantitatively the empirical data. This work identifies general mechanisms of norm change and we anticipate that it will be of interest to researchers investigating the cultural evolution of language and, more broadly, human collective behavior.


title: Privacy-Utility Management of Hypothesis Tests

abstract: The trade-off of hypothesis tests on the correlated privacy hypothesis and utility hypothesis is studied. The error exponent of the Bayesian composite hypothesis test on the privacy or utility hypothesis can be characterized by the corresponding minimal Chernoff information rate. An optimal management protects the privacy by minimizing the error exponent of the privacy hypothesis test and meanwhile guarantees the utility hypothesis testing performance by satisfying a lower bound on the corresponding minimal Chernoff information rate. The asymptotic minimum error exponent of the privacy hypothesis test is shown to be characterized by the infimum of corresponding minimal Chernoff information rates subject to the utility guarantees.


title: NbTi/Nb/Cu multilayer shield for the superconducting shield (SuShi) septum

abstract: A passive superconducting shield was proposed earlier to realize a high-field (3-4 T) septum magnet for the Future Circular Collider. This paper presents the experimental results of a potential shield material, a NbTi/Nb/Cu multilayer sheet. A cylindrical shield was constructed from two halves, each consisting of 4 layers with a total thickness of 3.2~mm, and inserted into the bore of a spare LHC dipole corrector magnet (MCBY). At 4.2~K, up to about 3.1~T at the shield's surface only a leakage field of 12.5~mT was measured inside the shield. This can be attributed to the mis-alignment of the two half cylinders, as confirmed by finite element simulations. With a better configuration we estimate the shield's attenuation to be better than $\mathbf{4\times 10^{-5}}$, acceptable for the intended application. Above 3.1~T the field penetrated smoothly. Below that limit no flux jumps were observed even at the highest achievable ramp rate of more than 50~mT/s at the shield's surface. A 'degaussing' cycle was used to eliminate the effects of the field trapped in the thick wall of the shield, which could otherwise distort the homogeneous field pattern at the extracted beam's position. At 1.9~K the shield's performance was superior to that at 4.2~K, but it suffered from flux jumps.


title: Controlled creation of nanometric skyrmions using external magnetic fields

abstract: To exploit nanometric magnetic skyrmions as information carriers in high-density storage devices, a method is needed that creates intended number of skyrmions at specified places in the device preferably at a low energy cost. We theoretically propose that using a system with a fabricated hole or notch, the controlled creation of individual skyrmions can be achieved even when using an external magnetic field applied to the entire specimen. The fabricated defect turns out to work like a catalyst to reduce the energy barrier for the skyrmion creation.


title: Deep Learning in Information Security

abstract: Machine learning has a long tradition of helping to solve complex information security problems that are difficult to solve manually. Machine learning techniques learn models from data representations to solve a task. These data representations are hand-crafted by domain experts. Deep Learning is a sub-field of machine learning, which uses models that are composed of multiple layers. Consequently, representations that are used to solve a task are learned from the data instead of being manually designed. In this survey, we study the use of DL techniques within the domain of information security. We systematically reviewed 77 papers and presented them from a data-centric perspective. This data-centric perspective reflects one of the most crucial advantages of DL techniques -- domain independence. If DL-methods succeed to solve problems on a data type in one domain, they most likely will also succeed on similar data from another domain. Other advantages of DL methods are unrivaled scalability and efficiency, both regarding the number of examples that can be analyzed as well as with respect of dimensionality of the input data. DL methods generally are capable of achieving high-performance and generalize well. However, information security is a domain with unique requirements and challenges. Based on an analysis of our reviewed papers, we point out shortcomings of DL-methods to those requirements and discuss further research opportunities.


title: A Malin 1 'cousin' with counter-rotation: internal dynamics and stellar content of the giant low surface brightness galaxy UGC 1922

abstract: The formation scenario for giant low surface brightness (gLSB) galaxies with discs as large as 100 kpc still remains unclear. These stellar systems are rare and very hard to observe, therefore a detailed insight on every additional object helps to understand their nature. Here we present a detailed observational study of the gLSB UGC 1922 performed using deep optical imaging and spectroscopic observations combined with archival ultraviolet data. We derived spatially resolved properties of stellar and ionized gas kinematics and characteristics of stellar populations and interstellar medium. We reveal the presence of a kinematically decoupled central component, which counter rotates with respect to the main disc of UGC 1922. The radial metallicity gradient of the ionised gas is in agreement with that found for moderate-size LSB galaxies. At the same time, a slowly rotating and dynamically hot central region of the galaxy hosts a large number of old metal-rich stars, which creates an appearance of a giant elliptical galaxy, that grew an enormous star forming disc. We reproduce most of the observed features of UGC 1922 in N-body/hydrodynamical simulations of an in-plane merger of giant Sa and Sd galaxies. We also discuss alternative formation scenarios of this unusual system.


title: Blind prediction of protein B-factor and flexibility

abstract: Debye-Waller factor, a measure of X-ray attenuation, can be experimentally observed in protein X-ray crystallography. Previous theoretical models have made strong inroads in the analysis of B-factors by linearly fitting protein B-factors from experimental data. However, the blind prediction of B-factors for unknown proteins is an unsolved problem. This work integrates machine learning and advanced graph theory, namely, multiscale weighted colored graphs (MWCGs), to blindly predict B-factors of unknown proteins. MWCGs are local features that measure the intrinsic flexibility due to a protein structure. Global features that connect the B-factors of different proteins, e.g., the resolution of X-ray crystallography, are introduced to enable the cross-protein B-factor predictions. Several machine learning approaches, including ensemble methods and deep learning, are considered in the present work. The proposed method is validated with hundreds of thousands of experimental B-factors. Extensive numerical results indicate that the blind B-factor predictions obtained from the present method are more accurate than the least squares fittings using traditional methods.


title: Graphyne as the anode material of magnesium-ion batteries: ab initio study

abstract: Graphyne, a single atomic layer structure of the carbon six-member rings connected by one acetilenic linkage, is a promising anode of rechargeable batteries. In this paper, a first-principle study has been carried out on graphyne as a new candidate for the anode material of magnesium-ion batteries, using density functional theory calculations. The main focus is on the magnesium adsorption on graphyne surface. The structural properties such as adsorption height and energy, the most stable adsorption sites, the Band structure and DOS of the pristine graphyne the diverse Mg-decorated graphyne structures, and energy barrier against Mg diffusion are also calculated. As a consequence of the band structure and DOS of graphyne structures, it is found that the pristine graphyne and the Mg-decorated graphyne structures show a semiconducting nature and metallic behavior, respectively. Moreover, the migration behavior of Mg on graphyne for the main diffusion paths is determined by the Nudged Elastic Band (NEB) method.


title: SDSS-IV MaNGA: The Formation Sequence of S0 Galaxies

abstract: Gas stripping of spiral galaxies or mergers are thought to be the formation mechanisms of lenticular galaxies. In order to determine the conditions in which each scenario dominates, we derive stellar populations of both the bulge and disk regions of 279 lenticular galaxies in the MaNGA survey. We find a clear bimodality in stellar age and metallicity within the population of S0s and this is strongly correlated with stellar mass. Old and metal-rich bulges and disks belong to massive galaxies, and young and metal-poor bulges and disks are hosted by low-mass galaxies. From this we conclude that the bulges and disks are co-evolving. When the bulge and disk stellar ages are compared, we find that the bulge is almost always older than the disk for massive galaxies ($\textrm{M}_{\star} > 10^{10}~\textrm{M}_{\odot}$). The opposite is true for lower mass galaxies. We conclude that we see two separate populations of lenticular galaxies. The old, massive, and metal-rich population possess bulges that are predominantly older than their disks, which we speculate may have been caused by morphological or inside-out quenching. In contrast, the less massive and more metal-poor population have bulges with more recent star formation than their disks. We postulate they may be undergoing bulge rejuvenation (or disk fading), or compaction. Environment doesn't play a distinct role in the properties of either population. Our findings give weight to the notion that while the faded spiral scenario likely formed low-mass S0s, other processes, such as mergers, may be responsible for high-mass S0s.


title: Adaptive Newton-Type Schemes Based on Projections

abstract: In this work we present and discuss a possible globalization concept for Newton-type methods. We consider nonlinear problems $f(x)=0$ in $\mathbb{R}^{n}$ using the concepts from ordinary differential equations as a basis for the proposed numerical solution procedure. Thus, the starting point of our approach is within the framework of solving ordinary differential equations numerically. Accordingly, we are able to reformulate general Newton-type iteration schemes using an adaptive step size control procedure. In doing so, we derive and discuss a discrete adaptive solution scheme thereby trying to mimic the underlying continuous problem numerically without losing the famous quadratic convergence regime of the classical Newton method in a vicinity of a regular solution. The derivation of the proposed adaptive iteration scheme relies on a simple orthogonal projection argument taking into account that, sufficiently close to regular solutions, the vector field corresponding to the Newton scheme is approximately linear. We test and exemplify our adaptive root-finding scheme using a few low-dimensional examples. Based on the presented examples, we finally show some performance data.


title: Game time: statistical contests in the classroom

abstract: We describe a contest in variable selection which was part of a statistics course for graduate students. In particular, the possibility to create a contest themselves offered an additional challenge for more advanced students. Since working with data is becoming more important in teaching statistics, we greatly encourage other instructors to try the same.


title: Concurrent Robin Hood Hashing

abstract: In this paper we examine the issues involved in adding concurrency to the Robin Hood hash table algorithm. We present a non-blocking obstruction-free K-CAS Robin Hood algorithm which requires only a single word compare-and-swap primitive, thus making it highly portable. The implementation maintains the attractive properties of the original Robin Hood structure, such as a low expected probe length, capability to operate effectively under a high load factor and good cache locality, all of which are essential for high performance on modern computer architectures. We compare our data-structures to various other lock-free and concurrent algorithms, as well as a simple hardware transactional variant, and show that our implementation performs better across a number of contexts.


title: On framed simple purely real Hurwitz numbers

abstract: We present a study of real Hurwitz numbers enumerating a special kind of real meromorphic functions, which we call simple framed purely real functions. We deduce partial differential equations of cut-and-join type for generating functions for these numbers. We also construct a topological field theory for them.


title: Stochastic Integral Representation for the Dynamics of Disordered Systems

abstract: The dynamics of interacting quantum systems in the presence of disorder is studied and an exact representation for disorder-averaged quantities via Ito stochastic calculus is obtained. The stochastic integral representation affords many advantages, including amenability to analytic approximation, applicability to interacting systems, and compatibility with existing tensor network methods. The integral may be expanded to produce a series of approximations, the first of which already includes all diffusive corrections and, further, is manifestly completely positive. The addition of fluctuations leads to a convergent series of systematic corrections. As examples, expressions for the density of states, spectral form factor, and out-of-time-order correlators for the Anderson model are obtained.


title: An asymptotic expansion for the error term in the Brent-McMillan algorithm for Euler's constant

abstract: The Brent-McMillan algorithm is the fastest known procedure for the high-precision computation of Euler's constant $\gamma$ and is based on the modified Bessel functions $I_0(2x)$ and $K_0(2x)$. An error estimate for this algorithm relies on the optimally truncated asymptotic expansion for the product $I_0(2x) K_0(2x)$ when $x$ assumes large positive integer values. An asymptotic expansion for this optimal error term is derived by exploiting the techniques developed in hyperasymptotics, thereby enabling more precise information on the error term than recently obtained bounds and estimates.


title: Compact Optimization Algorithms with Re-sampled Inheritance

abstract: Compact optimization algorithms are a class of Estimation of Distribution Algorithms (EDAs) characterized by extremely limited memory requirements (hence they are called "compact"). As all EDAs, compact algorithms build and update a probabilistic model of the distribution of solutions within the search space, as opposed to population-based algorithms that instead make use of an explicit population of solutions. In addition to that, to keep their memory consumption low, compact algorithms purposely employ simple probabilistic models that can be described with a small number of parameters. Despite their simplicity, compact algorithms have shown good performances on a broad range of benchmark functions and real-world problems. However, compact algorithms also come with some drawbacks, i.e. they tend to premature convergence and show poorer performance on non-separable problems. To overcome these limitations, here we investigate a possible algorithmic scheme obtained by combining compact algorithms with a non-disruptive restart mechanism taken from the literature, named Re-Sampled Inheritance (RI). The resulting compact algorithms with RI are tested on the CEC 2014 benchmark functions. The numerical results show on the one hand that the use of RI consistently enhances the performances of compact algorithms, still keeping a limited usage of memory. On the other hand, our experiments show that among the tested algorithms, the best performance is obtained by compact Differential Evolution with RI.


title: The Wisdom of MaSSeS: Majority, Subjectivity, and Semantic Similarity in the Evaluation of VQA

abstract: We introduce MASSES, a simple evaluation metric for the task of Visual Question Answering (VQA). In its standard form, the VQA task is operationalized as follows: Given an image and an open-ended question in natural language, systems are required to provide a suitable answer. Currently, model performance is evaluated by means of a somehow simplistic metric: If the predicted answer is chosen by at least 3 human annotators out of 10, then it is 100% correct. Though intuitively valuable, this metric has some important limitations. First, it ignores whether the predicted answer is the one selected by the Majority (MA) of annotators. Second, it does not account for the quantitative Subjectivity (S) of the answers in the sample (and dataset). Third, information about the Semantic Similarity (SES) of the responses is completely neglected. Based on such limitations, we propose a multi-component metric that accounts for all these issues. We show that our metric is effective in providing a more fine-grained evaluation both on the quantitative and qualitative level.


title: Disparity among low first ionization potential elements

abstract: The elemental composition of the solar wind differs from the solar photospheric composition. Elements with low first ionization potential (FIP) appear enhanced compared to O in the solar wind relative to the respective photospheric abundances. This so-called FIP effect is different in the slow solar wind and the coronal hole wind. However, under the same plasma conditions, for elements with similar FIPs such as Mg, Si, and Fe, comparable enhancements are expected. We scrutinize the assumption that the FIP effect is always similar for different low FIP elements, namely Mg, Si, and Fe. We investigate the dependency of the FIP effect of low FIP elements on the O7+/O6+ charge state ratio depending on time and solar wind type. We order the observed FIP ratios with respect to the O7+/O6+ charge state ratio into bins and analyze separately the respective distributions of the FIP ratio of Mg, Si, and Fe for each O7+/O6+ charge state ratio bin. We observe that the FIP effect shows the same qualitative yearly behavior for Mg and Si, while Fe shows significant differences during the solar activity maximum and its declining phase. In each year, the FIP effect for Mg and Si always increases with increasing O7+/O6+ charge state ratio, but for high O7+/O6+ charge state ratios the FIP effect for Fe shows a qualitatively different behavior. During the years 2001-2006, instead of increasing with the O7+/O6+ charge state ratio, the Fe FIP ratio exhibits a broad peak. Also, the FIP distribution per O7+/O6+ charge state bin is significantly broader for Fe than for Mg and Si. These observations support the conclusion that the elemental fractionation is only partly determined by FIP. In particular, the qualitative difference behavior with increasing O7+/O6+ charge state ratio between Fe on the one hand and Mg and Si on the other hand is not yet well explained by models of fractionation.


title: Characterisation of SiPMs

abstract: Silicon photomultipliers, thanks to their excellent performance, robustness and relatively simple use, are the photon-detectors of choice for many present and future applications. This paper gives an overview of methods to characterise SiPMs. The different SiPM parameters are introduced and generic setups for their determination presented. Finally, ways to extract the parameters from the measurements are discussed and the results shown. If a parameter can be obtained from different measurements, the results are compared and recommendations given, which is considered to be the most reliable. The characterisation of SiPMs, in particular for high light intensities and in high radiation fields, is presently a field of intensive research with many open questions and problems which will be discussed.


title: High-dimensional Bayesian Fourier Analysis For Detecting Circadian Gene Expressions

abstract: In genomic applications, there is often interest in identifying genes whose time-course expression trajectories exhibit periodic oscillations with a period of approximately 24 hours. Such genes are usually referred to as circadian, and their identification is a crucial step toward discovering physiological processes that are clock-controlled. It is natural to expect that the expression of gene i at time j might depend to some degree on the expression of the other genes measured at the same time. However, widely-used rhythmicity detection techniques do not accommodate for the potential dependence across genes. We develop a Bayesian approach for periodicity identification that explicitly takes into account the complex dependence structure across time-course trajectories in gene expressions. We employ a latent factor representation to accommodate dependence, while representing the true trajectories in the Fourier domain allows for inference on period, phase, and amplitude of the signal. Identification of circadian genes is allowed through a carefully chosen variable selection prior on the Fourier basis coefficients. The methodology is applied to a novel mouse liver circadian dataset. Although motivated by time-course gene expression array data, the proposed approach is applicable to the analysis of dependent functional data at broad.


title: A Bayesian seamless phase I-II trial design with two stages for cancer clinical trials with drug combinations

abstract: The use of drug combinations in clinical trials is increasingly common during the last years since a more favorable therapeutic response may be obtained by combining drugs. In phase I clinical trials, most of the existing methodology recommends a one unique dose combination as "optimal", which may result in a subsequent failed phase II clinical trial since other dose combinations may present higher treatment efficacy for the same level of toxicity. We are particularly interested in the setting where it is necessary to wait a few cycles of therapy to observe an efficacy outcome and the phase I and II population of patients are different with respect to treatment efficacy. Under these circumstances, it is common practice to implement two-stage designs where a set of maximum tolerated dose combinations is selected in a first stage, and then studied in a second stage for treatment efficacy. In this article we present a new two-stage design for early phase clinical trials with drug combinations. In the first stage, binary toxicity data is used to guide the dose escalation and set the maximum tolerated dose combinations. In the second stage, we take the set of maximum tolerated dose combinations recommended from the first stage, which remains fixed along the entire second stage, and through adaptive randomization, we allocate subsequent cohorts of patients in dose combinations that are likely to have high posterior median time to progression. The methodology is assessed with extensive simulations and exemplified with a real trial.


title: Systematic effects from black hole-neutron star waveform model uncertainties on the neutron star equation of state

abstract: We identify various contributors of systematic effects in the measurement of the neutron star (NS) tidal deformability and quantify their magnitude for several types of neutron star - black hole (NSBH) binaries. Gravitational waves from NSBH mergers contain information about the components' masses and spins as well as the NS equation of state. Extracting this information requires comparison of the signal in noisy detector data with theoretical templates derived from some combination of post-Newtonian (PN) approximants, effective one-body (EOB) models and %analytic fits to numerical relativity (NR) simulations. The accuracy of these templates is limited by errors in the NR simulations, by the approximate nature of the PN/EOB waveforms, and by the hybridization procedure used to combine them. In this paper, we estimate the impact of these errors by constructing and comparing a set of PN-NR hybrid waveforms, for the first time with NR waveforms from two different codes, namely, SpEC and SACRA, for such systems. We then attempt to recover the parameters of the binary using two non-precessing template approximants. We find that systematic errors are too large for tidal effects to be accurately characterized for any realistic NS equation of state model. We conclude that NSBH waveform models must be significantly improved if they are to be useful for the extraction of NS equation of state information or even for distinguishing NSBH systems from binary black holes.


title: Mergers, tidal interactions, and mass exchange in a population of disc globular clusters

abstract: We present the results of a self-consistent $N$-body simulation following the evolution of a primordial population of thick disc globular clusters (GCs). We study how the internal properties of such clusters evolve under the action of mutual interactions, while they orbit a Milky Way-like galaxy. For the first time, through analytical and numerical considerations, we find that physical encounters between disc GCs are a crucial factor that contributed to the shape of the current properties of the Galactic GC system. Close passages or motion on similar orbits may indeed have a significant impact on the internal structure of clusters, producing multiple gravitationally bound sub-populations through the exchange of mass and even mergers. Our model produces two major mergers and a few small mass exchanges between pairs of GCs. Two of our GCs accrete stars from two companions, ending up with three internal sub-populations. We propose these early interactions and mergers between thick disc GCs with slightly different initial chemical compositions as a possible explanation for the presence of the spreads in metallicity observed in some of the massive Milky Ways GCs.


title: Many-body effects on Landau-level spectra and cyclotron resonance in graphene

abstract: Recently Russell et al. [Phys. Rev. Lett. 120, 047401 (2018)] have reported a clear signal of many-particle contributions to cyclotron resonance in high-mobility hBN-encapsulated graphene, observing significant variations of resonance energies as a function of the filling factor $\nu$ for a series of interband channels. To elucidate their results, Coulombic contributions to the Landau-level spectra and cyclotron resonance in graphene are examined with a possible band gap taken into account and with emphasis on revealing electron-hole ($eh$) conjugation symmetry underlying such level and resonance spectra. Theory, based on the single-mode approximation, gives a practically good account of the experimental data; the data suggest a band gap of ~ 10 meV and show a profile that apparently reflects $eh$ conjugation symmetry.


title: Divide-and-Conquer Strategy for Large-Scale Eulerian Solvent Excluded Surface

abstract: Motivation: Surface generation and visualization are some of the most important tasks in biomolecular modeling and computation. Eulerian solvent excluded surface (ESES) software provides analytical solvent excluded surface (SES) in the Cartesian grid, which is necessary for simulating many biomolecular electrostatic and ion channel models. However, large biomolecules and/or fine grid resolutions give rise to excessively large memory requirements in ESES construction. We introduce an out-of-core and parallel algorithm to improve the ESES software. Results: The present approach drastically improves the spatial and temporal efficiency of ESES. The memory footprint and time complexity are analyzed and empirically verified through extensive tests with a large collection of biomolecule examples. Our results show that our algorithm can successfully reduce memory footprint through a straightforward divide-and-conquer strategy to perform the calculation of arbitrarily large proteins on a typical commodity personal computer. On multi-core computers or clusters, our algorithm can reduce the execution time by parallelizing most of the calculation as disjoint subproblems. Various comparisons with the state-of-the-art Cartesian grid based SES calculation were done to validate the present method and show the improved efficiency. This approach makes ESES a robust software for the construction of analytical solvent excluded surfaces. Availability and implementation: http://weilab.math.msu.edu/ESES.


title: Self-adjoint local boundary problems on compact surfaces. II. Family index

abstract: The paper presents a first step towards a family index theorem for classical self-adjoint boundary value problems. We address here the simplest non-trivial case of manifolds with boundary, namely the case of two-dimensional manifolds. The first result of the paper is an index theorem for families of first order self-adjoint elliptic differential operators with local boundary conditions, parametrized by points of a compact topological space $X$. We compute the $K^1(X)$-valued index in terms of the topological data over the boundary. The second result is the universality of the index: we show that the index is a universal additive homotopy invariant for such families, if the vanishing on families of invertible operators is required.


title: Robust Beamforming for AN Aided MISO SWIPT System with Unknown Eavesdroppers and Non-linear EH Model

abstract: This work studies a beamforming design for downlink transmission of a multi-user multiple-input single-output (MISO) system where each legitimate user employs a power splitting (PS) based simultaneous wireless information and power transfer (SWIPT) technique. The transmitter intends to send confidential information to its legitimate users in the presence of purely unknown eavesdroppers. Since the transmitter does not have any knowledge of the eavesdroppers' channel state information (CSI), we consider an artificial noise (AN) approach to establishing secure communication. This beamforming design is developed by maximizing the AN power to interfere with the eavesdropper as much as possible. Based on the assumption of imperfect CSI of legitimate users at the transmitter, two robust design approaches for the joint beamforming and PS ratio have been studied to maximize the AN power under both energy harvesting (EH) and signal-to-interference-plus-noise ratio (SINR) requirements at each legitimate user. In the first robust design, we consider the bounded channel uncertainties, and employ semidefinite relaxation (SDR) and a linear matrix inequality (LMI) representation to transform the original problem into a semidefinite program (SDP). In the second robust design, we consider the statistical channel uncertainties, and show that the proposed problem can be reformulated into another form of SDP through both SDR and Bernstein-type inequality. In addition, the non-linear energy harvesting (EH) model is incorporated in this work as it could reflect the characteristics of practical radio frequency(RF)-EH conversion circuit. Simulation results have been provided to demonstrate the performance of our proposed robust designs.


title: Packing Sporadic Real-Time Tasks on Identical Multiprocessor Systems

abstract: In real-time systems, in addition to the functional correctness recurrent tasks must fulfill timing constraints to ensure the correct behavior of the system. Partitioned scheduling is widely used in real-time systems, i.e., the tasks are statically assigned onto processors while ensuring that all timing constraints are met. The decision version of the problem, which is to check whether the deadline constraints of tasks can be satisfied on a given number of identical processors, has been known ${\cal NP}$-complete in the strong sense. Several studies on this problem are based on approximations involving resource augmentation, i.e., speeding up individual processors. This paper studies another type of resource augmentation by allocating additional processors, a topic that has not been explored until recently. We provide polynomial-time algorithms and analysis, in which the approximation factors are dependent upon the input instances. Specifically, the factors are related to the maximum ratio of the period to the relative deadline of a task in the given task set. We also show that these algorithms unfortunately cannot achieve a constant approximation factor for general cases. Furthermore, we prove that the problem does not admit any asymptotic polynomial-time approximation scheme (APTAS) unless ${\cal P}={\cal NP}$ when the task set has constrained deadlines, i.e., the relative deadline of a task is no more than the period of the task.


title: Deep learning for time series classification: a review

abstract: Time Series Classification (TSC) is an important and challenging problem in data mining. With the increase of time series data availability, hundreds of TSC algorithms have been proposed. Among these methods, only a few have considered Deep Neural Networks (DNNs) to perform this task. This is surprising as deep learning has seen very successful applications in the last years. DNNs have indeed revolutionized the field of computer vision especially with the advent of novel deeper architectures such as Residual and Convolutional Neural Networks. Apart from images, sequential data such as text and audio can also be processed with DNNs to reach state-of-the-art performance for document classification and speech recognition. In this article, we study the current state-of-the-art performance of deep learning algorithms for TSC by presenting an empirical study of the most recent DNN architectures for TSC. We give an overview of the most successful deep learning applications in various time series domains under a unified taxonomy of DNNs for TSC. We also provide an open source deep learning framework to the TSC community where we implemented each of the compared approaches and evaluated them on a univariate TSC benchmark (the UCR/UEA archive) and 12 multivariate time series datasets. By training 8,730 deep learning models on 97 time series datasets, we propose the most exhaustive study of DNNs for TSC to date.


title: Charged-particle decays of highly excited states in $^{19}$F

abstract: Neutron-capture reactions on $^{18}$F in the helium-burning shell play an important role in the production of $^{15}$N during core-collapse supernovae. The competition between the $^{18}$F($n,p/\alpha$)$^{18}$O/$^{15}$N reactions controls the amount of $^{15}$N produced. The strengths of these reactions depend on the decay branching ratios of states in $^{19}$F above the neutron threshold. We report on an experiment investigating the decay branching ratios of these states in order to better constrain the strengths of the reactions.


title: Effects of rotation and magnetic field on the revival of a stalled shock in supernova explosions

abstract: We investigate axisymmetric steady solutions of (magneto)hydrodynamics equations that describe approximately accretion flows through a standing shock wave and discuss the effects of rotation and magnetic field on the revival of the stalled shock wave in supernova explosions. We develop a new powerful numerical method to calculate the 2-dimensional (2D) steady accretion flows self-consistently. We first confirm the results of preceding papers that there is a critical luminosity of irradiating neutrinos, above which there exists no steady solution in spherical models. If a collapsing star has rotation and/or magnetic field, the accretion flows are no longer spherical owing to the centrifugal force and/or Lorentz force and the critical luminosity is modified.In fact we find that the critical luminosity is reduced by about 50% - 70% for rapid rotations and about 20% - 50% for strong toroidal magnetic fields, depending on the mass accretion rate. These results may be also interpreted as an existence of the critical specific angular momentum or critical magnetic field, above which there exists no steady solution and the standing shock wave will revive for a given combination of mass accretion rate and neutrino luminosity.


title: Training Deep Neural Networks with Different Datasets In-the-wild: The Emotion Recognition Paradigm

abstract: A novel procedure is presented in this paper, for training a deep convolutional and recurrent neural network, taking into account both the available training data set and some information extracted from similar networks trained with other relevant data sets. This information is included in an extended loss function used for the network training, so that the network can have an improved performance when applied to the other data sets, without forgetting the learned knowledge from the original data set. Facial expression and emotion recognition in-the-wild is the test bed application that is used to demonstrate the improved performance achieved using the proposed approach. In this framework, we provide an experimental study on categorical emotion recognition using datasets from a very recent related emotion recognition challenge.


title: From Solar to Stellar Brightness Variations: The Effect of Metallicity

abstract: Context. Comparison studies of Sun-like stars with the Sun suggest an anomalously low photometric variability of the Sun compared to Sun-like stars with similar magnetic activity. Comprehensive understanding of stellar variability is needed, to find a physical reasoning for this observation. Aims. We investigate the effect of metallicity and effective temperature on the photometric brightness change of Sun-like stars seen at different inclinations. The considered range of fundamental stellar parameters is sufficiently small so the stars, investigated here, still count as Sun-like or even as solar twins. Methods. To model the brightness change of stars with solar magnetic activity, we extend a well established model of solar brightness variations, SATIRE (which stands for Spectral And Total Irradiance Reconstruction), which is based on solar spectra, to stars with different fundamental parameters. For that we calculate stellar spectra for different metallicities and effective temperature using the radiative transfer code ATLAS9. Results. We show that even a small change (e.g. within the observational error range) of metallicity or effective temperature significantly affects the photometric brightness change compared to the Sun. We find that for Sun-like stars, the amplitude of the brightness variations obtained for Str\"omgren (b + y)/2 reaches a local minimum for fundamental stellar parameters close to the solar metallicity and effective temperature. Moreover, our results show that the effect of inclination decreases for metallicity values greater than the solar metallicity. Overall, we find that an exact determination of fundamental stellar parameters is crucially important for understanding stellar brightness changes.


title: Affleck-Dine baryogenesis in the SUSY Dine-Fischler-Srednicki-Zhitnitsky axion model without R-parity

abstract: We investigate the baryon asymmetry in the supersymmetry Dine-Fischler-Srednicki-Zhitnitsky axion model without R-parity. It turns out that the R-parity violating terms economically explain the atmospheric mass-squared difference of neutrinos and the appropriate amount of baryon asymmetry through the Affleck-Dine mechanism. In this model, the axion is a promising candidate for the dark matter and the axion isocurvature perturbation is suppressed due to the large field values of Peccei-Quinn fields. Remarkably, in some parameter regions explaining the baryon asymmetry and the axion dark matter abundance, the proton decay will be explored in future experiments.


title: Iterative Delegations in Liquid Democracy with Restricted Preferences

abstract: In this paper, we study liquid democracy, a collective decision making paradigm which lies between direct and representative democracy. One main feature of liquid democracy is that voters can delegate their votes in a transitive manner so that: A delegates to B and B delegates to C leads to A delegates to C. Unfortunately, this process may not converge as there may not even exist a stable state (also called equilibrium). In this paper, we investigate the stability of the delegation process in liquid democracy when voters have restricted types of preference on the agent representing them (e.g., single-peaked preferences). We show that various natural structures of preferences guarantee the existence of an equilibrium and we obtain both tractability and hardness results for the problem of computing several equilibria with some desirable properties.


title: Optimality conditions and complete description of polytopes in combinatorial optimization

abstract: A combinatorial optimization problem (COP) has a finite groundset $E(\left|E\right|=N$), a weight vector $c=\left(c^e:e\in E\right)$ and a family $T\in E$ of feasible subsets with objective to find $t\in T$ with maximal weight: ${max}\{\sum_{e\in t}c^e$: $t\in T\}$. Polyhedral combinatorics reformulates combinatorial optimization as linear program: $T$ is mapped into the set $X\in R^N$ of 0/1 incidence vectors and $c\in R^N$ is maximized over the convex hull of $X$: ${max}\{cx: x\in conv(X)\}$. In theory, complementary slackness conditions for the induced linear program provide optimality conditions for the COP. However, in general case, optimality conditions for combinatorial optimization have not been formulated analytically as for many problems complete description of the induced polytopes is available only as a convex hull of extreme points rather than a system of linear inequalities. Here, we formulate optimality conditions for a COP in general case: $x_k\in X\ \ $is optimal if and only if $c\in cone\left(H_k\right)\ $ where $H_k$ $=\{h\in V:\ {hx}_k\ \ge hx$ for all $x\in X\}$ and $V$ is the set of all -1/0/1 valued vectors in $R^N$. This provides basis to get, in theory, a complete description of a combinatorial polytope induced by any COP: all facet inducing inequalities for $conv(X)$ can be written as $hx\le l$ where $h\in V$ and $l$ is integer. A vector $h\in V$ induces a nonredundent facet if and only if $h\in H_k^o$ $\in H_k$ for at least one $x_k\in X\ $(where $H_k^o$$=\{h\in H_k:h\notin cone(H_k\setminus \{h\})\}$ ) and $l=\ $ $x_kh$.


title: Thermal Features for Presentation Attack Detection in Hand Biometrics

abstract: This paper proposes a method for utilizing thermal features of the hand for the purpose of presentation attack detection (PAD) that can be employed in a hand biometrics system's pipeline. By envisaging two different operational modes of our system, and by employing a DCNN-based classifiers fine-tuned with a dataset of real and fake hand representations captured in both visible and ther- mal spectrum, we were able to bring two important deliverables. First, a PAD method operating in an open-set mode, capable of correctly discerning 100% of fake thermal samples, achieving Attack Presentation Classification Error Rate (APCER) and Bona-Fide Presentation Classification Error Rate (BPCER) equal to 0%, which can be easily implemented into any existing system as a separate component. Second, a hand biometrics system operating in a closed-set mode, that has PAD built right into the recognition pipeline, and operating simultaneously with the user-wise classification, achieving rank-1 recognition accuracy of up to 99.75%. We also show that thermal images of the human hand, in addition to liveness features they carry, can also improve classification accuracy of a biometric system, when coupled with visible light images. To follow the reproducibility guidelines and to stimulate further research in this area, we share the trained model weights, source codes, and a newly created dataset of fake hand representations with interested researchers.


title: Predicting citation counts based on deep neural network learning techniques

abstract: With the growing number of published scientific papers world-wide, the need to evaluation and quality assessment methods for research papers is increasing. Scientific fields such as scientometrics, informetrics and bibliometrics establish quantified analysis methods and measurements for scientific papers. In this area, an important problem is to predict the future influence of a published paper. Particularly, early discrimination between influential papers and insignificant papers may find important applications. In this regard, one of the most important metrics is the number of citations to the paper, since this metric is widely utilized in the evaluation of scientific publications and moreover, it serves as the basis for many other metrics such as h-index. In this paper, we propose a novel method for predicting long-term citations of a paper based on the number of its citations in the first few years after publication. In order to train a citations prediction model, we employed artificial neural networks which is a powerful machine learning tool with recently growing applications in many domains including image and text processing. The empirical experiments show that our proposed method out-performs state-of-the-art methods with respect to the prediction accuracy in both yearly and total prediction of the number of citations.


title: Photodynamical analysis of the triply eclipsing hierarchical triple system EPIC 249432662

abstract: Using Campaign 15 data from the K2 mission, we have discovered a triply-eclipsing triple star system: EPIC 249432662. The inner eclipsing binary system has a period of 8.23 days, with shallow $\sim$3% eclipses. During the entire 80-day campaign, there is also a single eclipse event of a third-body in the system that reaches a depth of nearly 50% and has a total duration of 1.7 days, longer than for any previously known third-body eclipse involving unevolved stars. The binary eclipses exhibit clear eclipse timing variations. A combination of photodynamical modeling of the lightcurve, as well as seven follow-up radial velocity measurements, has led to a prediction of the subsequent eclipses of the third star with a period of 188 days. A campaign of follow-up ground-based photometry was able to capture the subsequent pair of third-body events as well as two further 8-day eclipses. A combined photo-spectro-dynamical analysis then leads to the determination of many of the system parameters. The 8-day binary consists of a pair of M stars, while most of the system light is from a K star around which the pair of M stars orbits.


title: Non-equilibrium fluctuations for the SSEP with a slow bond

abstract: We prove the non-equilibrium fluctuations for the one-dimensional symmetric simple exclusion process with a slow bond. This generalizes a result of T. Franco, A. Neumann and P. Gon\c{c}alves (2013), which dealt with the equilibrium fluctuations. The foundation stone of our proof is a precise estimate on the correlations of the system, and that is by itself one of the main novelties of this paper. To obtain these estimates, we first deduce a spatially discrete PDE for the covariance function and we relate it to the local times of a random walk in a non-homogeneous environment via Duhamel's principle. Projection techniques and coupling arguments reduce the analysis to the problem of studying the local times of the classical random walk. We think that the method developed here can be applied to a variety of models, and we provide a discussion on this matter.


title: Symmetries of special 2-flags

abstract: This work is a continuation of authors' research interrupted in the year 2010. Derived are recursive relations describing for the first time all infinitesimal symmetries of special 2-flags (sometimes also misleadingly called `Goursat 2-flags'). When algorithmized to the software level, they will give an answer filling in the gap in knowledge as of 2010: on one side the local finite classification of special 2-flags known in lengths not exceeding four, on the other side the existence of a continuous numerical modulus of that classification in length seven.


title: Thick points of random walk and the Gaussian free field

abstract: We consider the thick points of random walk, i.e. points where the local time is a fraction of the maximum. In two dimensions, we answer a question of Dembo, Peres, Rosen and Zeitouni and compute the number of thick points of planar random walk, assuming that the increments are symmetric and have a finite moment of order two. The proof provides a streamlined argument based on the connection to the Gaussian free field and works in a very general setting including isoradial graphs. In higher dimensions, we study the scaling limit of the thick points. In particular, we show that the rescaled number of thick points converges to a nondegenerate random variable and that the centred maximum of the local times converges to a randomly shifted Gumbel distribution.


title: First steps towards the reconstruction of the squark flavour structure

abstract: Assuming the observation of a squark at the Large Hadron Collider, we investigate methods to access its flavour content and thus gain information on the underlying flavour structure of the theory. Based on simple observables, we apply a likelihood inference method to determine the top-flavour content of the observed particle. In addition, we employ a multivariate analysis in order to classify different flavour hypotheses. Both methods are discussed within a simplified model and the more general Minimal Supersymmetric Standard Model including most general squark mixing. We conclude that the likelihood inference may provide an estimation of the top-flavour content if additional knowledge, especially on the gaugino sector is available, while the multivariate analysis identifies different flavour patterns and can accommodate a more minimalistic set of observables.


title: Exceptional points for photon pairs bound by nonlinear dissipation in cavity arrays

abstract: We study theoretically the dissipative Bose-Hubbard model describing array of tunneling-coupled cavities with non-conservative photon-photon interaction. Our calculation of the complex energy spectrum for the photon pairs reveals exceptional points where the two-photon states bound by nonlinear dissipation are formed. This improves fundamental understanding of the interplay of non-Hermiticity and interactions in the quantum structures and can be potentially used for on-demand nonlinear light generation in photonic lattices.


title: Anomalous increasing of the intensity of field dependence optical mode ferromagnetic resonance in the exchange coupled bilayer system

abstract: Acoustic and optical ferromagnetic resonance (FMR) in the interlayer exchange coupled Fe$_{20}$Ni$_{80}$/Co bilayer have been investigated. In the optical mode, unexpected increasing tendencies of peak value at the resonance frequency has been observed under an increasing magnetic field. We presented analytical calculations with which the exchange coupling between Co and Fe$_{20}$Ni$_{80}$ layers, the magnetization and the in-plane uniaxial anisotropy are taken into account, to interpret the increasing of the maximum values of the optical permeability. Both experimental measurements and theoretical calculation show that such tendencies are dependent on the layer thickness t, and that there is a critical field above which the optical peak value begins to decrease. These results might help us to understand the mechanism of interlayer exchange coupling induced optical FMR and might enlighten us to find new possibility of high frequency applications of magnetic materials.


title: Global smoothness for a 1D supercritical transport model with nonlocal velocity

abstract: We are concerned with a nonlocal transport 1D-model with supercritical dissipation $\gamma\in(0,1)$ in which the velocity is coupled via the Hilbert transform, namely the so-called CCF model. This model arises as a lower dimensional model for the famous 2D dissipative quasi-geostrophic equation and in connection with vortex-sheet problems. It is known that its solutions can blow up in finite time when $\gamma\in(0,1/2)$. On the other hand, as stated by Kiselev (2010), in the supercritical subrange $\gamma\in\lbrack1/2,1)$ it is an open problem to know whether its solutions are globally regular. We show global existence of non-negative $H^{3/2}$-strong solutions in a supercritical subrange (close to 1) that depends on the initial data norm. Then, for each arbitrary smooth non-negative initial data, the model has a unique global smooth solution provided that $\gamma\in\lbrack\gamma_{1},1)$ where $\gamma_{1}$ depends on the $H^{3/2}$-initial data norm. Our approach is inspired on that of Coti Zelati and Vicol (IUMJ, 2016).


title: A Collaborative Multi-agent Reinforcement Learning Anti-jamming Algorithm in Wireless Networks

abstract: In this letter, we investigate the anti-jamming defense problem in multi-user scenarios, where the coordination among users is taken into consideration. The Markov game framework is employed to model and analyze the anti-jamming defense problem, and a collaborative multi-agent anti-jamming algorithm (CMAA) is proposed to obtain the optimal anti-jamming strategy. In sweep jamming scenarios, on the one hand, the proposed CMAA can tackle the external malicious jamming. On the other hand, it can effectively cope with the mutual interference among users. Simulation results show that the proposed CMAA is superior to both sensing based method and independent Q-learning method, and has the highest normalized rate.


title: The PAU Survey: Early demonstration of photometric redshift performance in the COSMOS field

abstract: The PAU Survey (PAUS) is an innovative photometric survey with 40 narrow bands at the William Herschel Telescope (WHT). The narrow bands are spaced at 100\AA\ intervals covering the range 4500\AA\ to 8500\AA\ and, in combination with standard broad bands, enable excellent redshift precision. This paper describes the technique, galaxy templates and additional photometric calibration used to determine early photometric redshifts from PAUS. Using BCNz2, a new photometric redshift code developed for this purpose, we characterise the photometric redshift performance using PAUS data on the COSMOS field. Comparison to secure spectra from zCOSMOS DR3 shows that PAUS achieves $\sigma_{68} /(1+z) = 0.0037$ to $i_{\mathrm{AB}} < 22.5$ when selecting the best 50% of the sources based on a photometric redshift quality cut. Furthermore, a higher photo-z precision ($\sigma_{68}/(1+z) \sim 0.001$) is obtained for a bright and high quality selection, which is driven by the identification of emission lines. We conclude that PAUS meets its design goals, opening up a hitherto uncharted regime of deep, wide, and dense galaxy survey with precise redshifts that will provide unique insights into the formation, evolution and clustering of galaxies, as well as their intrinsic alignments.


title: A note on three-point functions of unprotected operators

abstract: Given the recent progress in computing three-point functions in N=4 SYM via integrability, I provide here a novel direct calculation of some structure constants at weak coupling. The main focus is on correlators involving more than one unprotected operator, at two-loop order in the perturbative expansion.


title: Detection of time reversibility in time series by ordinal patterns analysis

abstract: Time irreversibility is a common signature of nonlinear processes, and a fundamental property of non-equilibrium systems driven by non-conservative forces. A time series is said to be reversible if its statistical properties are invariant regardless of the direction of time. Here we propose the Time Reversibility from Ordinal Patterns method (TiROP) to assess time-reversibility from an observed finite time series. TiROP captures the information of scalar observations in time forward, as well as its time-reversed counterpart by means of ordinal patterns. The method compares both underlying information contents by quantifying its (dis)-similarity via Jensen-Shannon divergence. The statistic is contrasted with a population of divergences coming from a set of surrogates to unveil the temporal nature and its involved time scales. We tested TiROP in different synthetic and real, linear and non linear time series, juxtaposed with results from the classical Ramsey's time reversibility test. Our results depict a novel, fast-computation, and fully data-driven methodology to assess time-reversibility at different time scales with no further assumptions over data. This approach adds new insights about the current non-linear analysis techniques, and also could shed light on determining new physiological biomarkers of high reliability and computational efficiency.


title: On the Boris solver in particle-in-cell simulation

abstract: A simple form of the Boris solver in particle-in-cell (PIC) simulation is proposed. It employs an exact solution of the Lorentz-force part, and it is equivalent to the Boris solver with a gyrophase correction. As a favorable property for stable schemes, this form preserves a volume in the phase space. Numerical tests of the Boris solvers are conducted by test-particle simulations and by PIC simulations. The proposed form provides better accuracy than the popular form, while it only requires few additional computation time.

categories: physics.comp-ph astro-ph.HE physics.plasm-ph physics.space-ph


title: Bayesian Semi-supervised Learning with Graph Gaussian Processes

abstract: We propose a data-efficient Gaussian process-based Bayesian approach to the semi-supervised learning problem on graphs. The proposed model shows extremely competitive performance when compared to the state-of-the-art graph neural networks on semi-supervised learning benchmark experiments, and outperforms the neural networks in active learning experiments where labels are scarce. Furthermore, the model does not require a validation data set for early stopping to control over-fitting. Our model can be viewed as an instance of empirical distribution regression weighted locally by network connectivity. We further motivate the intuitive construction of the model with a Bayesian linear model interpretation where the node features are filtered by an operator related to the graph Laplacian. The method can be easily implemented by adapting off-the-shelf scalable variational inference algorithms for Gaussian processes.


title: Binary MDS Array Codes with Optimal Repair

abstract: Consider a binary maximum distance separable (MDS) array code composed of an $m\times (k+r)$ array of bits with $k$ information columns and $r$ parity columns, such that any $k$ out of $k+r$ columns suffice to reconstruct the $k$ information columns. Our goal is to provide {\em optimal repair access} for binary MDS array codes, meaning that the bandwidth triggered to repair any single failed information or parity column is minimized. In this paper, we propose a generic transformation framework for binary MDS array codes, using EVENODD codes as a motivating example, to support optimal repair access for $k+1\le d \le k+r-1$, where $d$ denotes the number of non-failed columns that are connected for repair; note that when $d<k+r-1$, some of the chosen $d$ columns in repairing a failed column are specific. In addition, we show how our transformation framework applies to an example of binary MDS array codes with asymptotically optimal repair access of any single information column and enables asymptotically or exactly optimal repair access for any column. Furthermore, we present a new transformation for EVENODD codes with two parity columns such that the existing efficient repair property of any information column is preserved and the repair access of parity column is optimal.


title: Fermionic degeneracy and non-local contributions in flag-dipole spinors and mass dimension one fermions

abstract: We construct a mass dimension one fermionic field associated with flag-dipole spinors. These spinors are related to Elko (flag-pole spinors) by a one-parameter matrix transformation $\mathcal{Z}(z)$ where $z$ is a complex number. The theory is non-local and non-covariant. While it is possible to obtain a Lorentz-invariant theory via $\tau$-deformation, we choose to study the effects of non-locality and non-covariance. Our motivation for doing so is explained. We show that a fermionic field with $|z|\neq1$ and $|z|=1$ are physically equivalent. But for fermionic fields with more than one value of $z$, their interactions are $z$-dependent thus introducing an additional fermionic degeneracy that is absent in the Lorentz-invariant theory. We study the fermionic self-interaction and the local $U(1)$ interaction. In the process, we obtained non-local contributions for fermionic self-interaction that have previously been neglected. For the local $U(1)$ theory, the interactions contain time derivatives that renders the interacting density non-commutative at space-like separation. We show that this problem can be resolved by working in the temporal gauge. This issue is also discussed in the context of gravity.


title: A Framework for Approval-based Budgeting Methods

abstract: We define and study a general framework for approval-based budgeting methods and compare certain methods within this framework by their axiomatic and computational properties. Furthermore, we visualize their behavior on certain Euclidean distributions and analyze them experimentally.


title: On convergence of Chorin's projection method to a Leray-Hopf weak solution

abstract: The projection method to solve the incompressible Navier-Stokes equations was first studied by Chorin [Math. Comp., 1969] in the framework of a finite difference method and Temam [Arch. Rational Mech. and Anal., 1969] in the framework of a finite element method. Chorin showed convergence of approximation and its error estimates in problems with the periodic boundary condition assuming existence of a $C^5$-solution, while Temam demonstrated an abstract argument to obtain a Leray-Hopf weak solution in problems on a bounded domain with the no-slip boundary condition. In the present paper, the authors extend Chorin's result with full details to obtain convergent finite difference approximation of a Leray-Hopf weak solution to the incompressible Navier-Stokes equations on an arbitrary bounded Lipschitz domain of $\mathbb{R}^3$ with the no-slip boundary condition and an external force. We prove unconditional solvability of our implicit scheme and strong $L^2$-convergence (up to subsequence) under the scaling condition $ h^{3-\alpha}\le\tau$ (no upper bound is necessary), where $h,\tau$ are space, time discretization parameters, respectively, and $\alpha\in(0,2]$ is any fixed constant. The results contain a compactness method based on a new interpolation inequality for step functions.


title: Hybrid matrix compression for high-frequency problems

abstract: Boundary element methods for the Helmholtz equation lead to large dense matrices that can only be handled if efficient compression techniques are used. Directional compression techniques can reach good compression rates even for high-frequency problems. Currently there are two approaches to directional compression: analytic methods approximate the kernel function, while algebraic methods approximate submatrices. Analytic methods are quite fast and proven to be robust, while algebraic methods yield significantly better compression rates. We present a hybrid method that combines the speed and reliability of analytic methods with the good compression rates of algebraic methods.


title: A new beamline for laser spin-polarization at ISOLDE

abstract: A beamline dedicated to the production of laser-polarized radioactive beams has been constructed at ISOLDE, CERN. We present here different simulations leading to the design and construction of the setup, as well as technical details of the full setup and examples of the achieved polarizations for several radioisotopes. Beamline simulations show a good transmission through the entire line, in agreement with observations. Simulations of the induced nuclear spin-polarization as a function of atom-laser interaction length are presented for $^{26,28}$Na, [1] and for $^{35}$Ar, which is studied in this work. Adiabatic spin rotation of the spin-polarized ensemble of atoms, and how this influences the observed nuclear ensemble polarization, are also performed for the same nuclei. For $^{35}$Ar, we show that multiple-frequency pumping enhances the ensemble polarization by a factor 1.85, in agreement with predictions from a rate equations model. [1] J. Phys. G: Nucl. Part. Phys./1744084005


title: On the regular-convexity of Ricci shrinker limit spaces

abstract: In this paper, we study the structure of the pointed-Gromov-Hausdorff limits of sequences of Ricci shrinkers. We define a regular-singular decomposition following the work of Cheeger-Colding for manifolds with a uniform Ricci curvature lower bound, and prove that the regular part of any Ricci shrinker limit space is convex, inspired by Colding-Naber's original idea of parabolic smoothing of the distance functions.


title: Quantum version of Wielandt's Inequality revisited

abstract: Consider a linear space L of complex D-dimensional linear operators, and assume that some power L^k of L is the whole space of DxD matrices. Perez-Garcia, Verstraete, Wolf and Cirac conjectured that the sequence L^1,L^2,... stablilizes after O(D^2) terms; we prove that this happens after O(D^2 log(D)) terms, improving the previously known bound of O(D^4).


title: Weak convergence of a measure-valued process for social networks

abstract: This article formalizes the problem of modeling social networks into an interacting particle system on random geometric graphs. Each vertex of the graph is associated with a geometric position in a latent space describing the unobserved affinity between the network members and characterizes the strength of the interaction between them. We endow the system with two recruitment mechanisms that depend on the positions of the particles in the latent space and one departure mechanism independent of the latent position. We characterize each spatial position by a Dirac measure and the system state by a measure-valued process which is the sum of all Dirac masses. Therefore, we investigate the large-scale behavior of the system. In particular, using a renormalization technique, we study the system's behavior when the initial number of particles goes to infinity. We thus show the weak convergence of the rescaled measure-valued process towards the solution of a deterministic integrodifferential equation. Finally, we present some Monte Carlo simulations with different parameter sets to illustrate the merit of our approach in modeling some phenomena encountered in real-world networks.


title: Spatio-Temporal Data Fusion for Massive Sea Surface Temperature Data from MODIS and AMSR-E Instruments

abstract: Remote sensing data have been widely used to study various geophysical processes. With the advances in remote-sensing technology, massive amount of remote sensing data are collected in space over time. Different satellite instruments typically have different footprints, measurement-error characteristics, and data coverages. To combine datasets from different satellite instruments, we propose a dynamic fused Gaussian process (DFGP) model that enables fast statistical inference such as filtering and smoothing for massive spatio-temporal datasets in a data-fusion context. Based upon a spatio-temporal-random-effects model, the DFGP methodology represents the underlying true process with two components: a linear combination of a small number of basis functions and random coefficients with a general covariance matrix, together with a linear combination of a large number of basis functions and Markov random coefficients. To model the underlying geophysical process at different spatial resolutions, we rely on the change-of-support property, which also allows efficient computations in the DFGP model. To estimate model parameters, we devise a computationally efficient stochastic expectation-maximization (SEM) algorithm to ensure its scalability for massive datasets. The DFGP model is applied to a total of 3.7 million sea surface temperature datasets in the tropical Pacific Ocean for a one-week time period in 2010 from MODIS and AMSR-E instruments.


title: A Hand Combining Two Simple Grippers to Pick up and Arrange Objects for Assembly

abstract: This paper proposes a novel robotic hand design for assembly tasks. The idea is to combine two simple grippers -- an inner gripper which is used for precise alignment, and an outer gripper which is used for stable holding. Conventional robotic hands require complicated compliant mechanisms or complicated control strategy and force sensing to conduct assemble tasks, which makes them costly and difficult to pick and arrange small objects like screws or washers. Compared to the conventional hands, the proposed design provides a low-cost solution for aligning, picking up, and arranging various objects by taking advantages of the geometric constraints of the positioning fingers and gravity. It is able to deal with small screws and washers, and eliminate the position errors of cylindrical objects or objects with cylindrical holes. In the experiments, both real-world tasks and quantitative analysis are performed to validate the aligning, picking, and arrangements abilities of the design.


title: Quantum effects of operator time ordering in the nonlinear Jaynes-Cummings model

abstract: Recently, in [Phys. Rev. A 97, 043806 (2018)], the detuned and nonlinear Jaynes-Cummings model describing the quantized motion of a trapped ion was introduced and its corresponding dynamics was solved via considering the driving laser in a quantized manner. In this work we reconsider this model and show that it can likewise be solved with a classical driving laser field. Using the exact solution we investigate the quantum time-ordering effects of the system with respect to nonclassicality of the motional states of the ion. Furthermore, we use the Magnus expansion to analyze the impact of certain orders of the time ordering and derive and exact radius of convergence beyond the established and only sufficient criterion. Finally, the differences of the solution derived here and the previously found one using a quantized pump, are discussed.


title: Machine Learning tools for global PDF fits

abstract: The use of machine learning algorithms in theoretical and experimental high-energy physics has experienced an impressive progress in recent years, with applications from trigger selection to jet substructure classification and detector simulation among many others. In this contribution, we review the machine learning tools used in the NNPDF family of global QCD analyses. These include multi-layer feed-forward neural networks for the model-independent parametrisation of parton distributions and fragmentation functions, genetic and covariance matrix adaptation algorithms for training and optimisation, and closure testing for the systematic validation of the fitting methodology.


title: Maximizing the Diversity of Exposure in a Social Network

abstract: Social-media platforms have created new ways for citizens to stay informed and participate in public debates. However, to enable a healthy environment for information sharing, social deliberation, and opinion formation, citizens need to be exposed to sufficiently diverse viewpoints that challenge their assumptions, instead of being trapped inside filter bubbles. In this paper, we take a step in this direction and propose a novel approach to maximize the diversity of exposure in a social network. We formulate the problem in the context of information propagation, as a task of recommending a small number of news articles to selected users. We propose a realistic setting where we take into account content and user leanings, and the probability of further sharing an article. This setting allows us to capture the balance between maximizing the spread of information and ensuring the exposure of users to diverse viewpoints. The resulting problem can be cast as maximizing a monotone and submodular function subject to a matroid constraint on the allocation of articles to users. It is a challenging generalization of the influence maximization problem. Yet, we are able to devise scalable approximation algorithms by introducing a novel extension to the notion of random reverse-reachable sets. We experimentally demonstrate the efficiency and scalability of our algorithm on several real-world datasets.


title: Induced path factors of regular graphs

abstract: An induced path factor of a graph $G$ is a set of induced paths in $G$ with the property that every vertex of $G$ is in exactly one of the paths. The induced path number $\rho(G)$ of $G$ is the minimum number of paths in an induced path factor of $G$. We show that if $G$ is a connected cubic graph on $n>6$ vertices, then $\rho(G)\le(n-1)/3$. Fix an integer $k\ge3$. For each $n$, define $\mathcal{M}_n$ to be the maximum value of $\rho(G)$ over all connected $k$-regular graphs $G$ on $n$ vertices. As $n\rightarrow\infty$ with $nk$ even, we show that $c_k=\lim(\mathcal{M}_n/n)$ exists. We prove that $5/18\le c_3\le1/3$ and $3/7\le c_4\le1/2$ and that $c_k=\frac12-O(k^{-1})$ for $k\rightarrow\infty$.


title: On the relation between pole and running heavy quark masses beyond the four-loop approximation

abstract: The effective charges motivated method is applied to the relation between pole and $\rm{\overline{MS}}$-scheme heavy quark masses to study high order perturbative QCD corrections in the observable quantities proportional to the running quark masses. The non-calculated five- and six-loop perturbative QCD coefficients are estimated. This approach predicts for these terms the sign-alternating expansion in powers of number of lighter flavors $n_l$, while the analyzed recently infrared renormalon asymptotic expressions do not reproduce the same behavior. We emphasize that coefficients of the quark mass relation contain proportional to $\pi^2$ effects, which result from analytical continuation from the Euclidean region, where the scales of the running masses and QCD coupling constant are initially fixed, to the Minkowskian region, where the pole masses and the running QCD parameters are determined. For the $t$-quark the asymptotic nature of the non-resummed PT mass relation does not manifest itself at six-loops, while for the $b$-quark the minimal PT term appears at the probed by direct calculations four-loop level. The recent infrared renormalon based studies support these conclusions.


title: Finding Cheeger Cuts in Hypergraphs via Heat Equation

abstract: Cheeger's inequality states that a tightly connected subset can be extracted from a graph $G$ using an eigenvector of the normalized Laplacian associated with $G$. More specifically, we can compute a subset with conductance $O(\sqrt{\phi_G})$, where $\phi_G$ is the minimum conductance of a set in $G$. It has recently been shown that Cheeger's inequality can be extended to hypergraphs. However, as the normalized Laplacian of a hypergraph is no longer a matrix, we can only approximate to its eigenvectors; this causes a loss in the conductance of the obtained subset. To address this problem, we here consider the heat equation on hypergraphs, which is a differential equation exploiting the normalized Laplacian. We show that the heat equation has a unique solution and that we can extract a subset with conductance $\sqrt{\phi_G}$ from the solution. An analogous result also holds for directed graphs.

categories: cs.DS cs.DM cs.NA math.AP math.NA math.SP


title: Isolated and Ensemble Audio Preprocessing Methods for Detecting Adversarial Examples against Automatic Speech Recognition

abstract: An adversarial attack is an exploitative process in which minute alterations are made to natural inputs, causing the inputs to be misclassified by neural models. In the field of speech recognition, this has become an issue of increasing significance. Although adversarial attacks were originally introduced in computer vision, they have since infiltrated the realm of speech recognition. In 2017, a genetic attack was shown to be quite potent against the Speech Commands Model. Limited-vocabulary speech classifiers, such as the Speech Commands Model, are used in a variety of applications, particularly in telephony; as such, adversarial examples produced by this attack pose as a major security threat. This paper explores various methods of detecting these adversarial examples with combinations of audio preprocessing. One particular combined defense incorporating compressions, speech coding, filtering, and audio panning was shown to be quite effective against the attack on the Speech Commands Model, detecting audio adversarial examples with 93.5% precision and 91.2% recall.


title: Optimal Strong Convergence Rate of a Backward Euler Type Scheme for the Cox--Ingersoll--Ross Model Driven by Fractional Brownian Motion

abstract: In this paper, we investigate the optimal strong convergence rate of numerical approximations for the Cox--Ingersoll--Ross model driven by fractional Brownian motion with Hurst parameter $H\in(1/2,1)$. To deal with the difficulties caused by the unbounded diffusion coefficient, we study an auxiliary equation based on Lamperti transformation. By means of Malliavin calculus, we prove that the backward Euler scheme applied to this auxiliary equation ensures the positivity of the numerical solution, and is of strong order one. Furthermore, a numerical approximation for the original model is obtained and converges with the same order.


title: Artificial Intelligence for the Public Sector: Opportunities and challenges of cross-sector collaboration

abstract: Public sector organisations are increasingly interested in using data science and artificial intelligence capabilities to deliver policy and generate efficiencies in high uncertainty environments. The long-term success of data science and AI in the public sector relies on effectively embedding it into delivery solutions for policy implementation. However, governments cannot do this integration of AI into public service delivery on their own. The UK Government Industrial Strategy is clear that delivering on the AI grand challenge requires collaboration between universities and public and private sectors. This cross-sectoral collaborative approach is the norm in applied AI centres of excellence around the world. Despite their popularity, cross-sector collaborations entail serious management challenges that hinder their success. In this article we discuss the opportunities and challenges from AI for public sector. Finally, we propose a series of strategies to successfully manage these cross-sectoral collaborations.


title: Learning Deep Mixtures of Gaussian Process Experts Using Sum-Product Networks

abstract: While Gaussian processes (GPs) are the method of choice for regression tasks, they also come with practical difficulties, as inference cost scales cubic in time and quadratic in memory. In this paper, we introduce a natural and expressive way to tackle these problems, by incorporating GPs in sum-product networks (SPNs), a recently proposed tractable probabilistic model allowing exact and efficient inference. In particular, by using GPs as leaves of an SPN we obtain a novel flexible prior over functions, which implicitly represents an exponentially large mixture of local GPs. Exact and efficient posterior inference in this model can be done in a natural interplay of the inference mechanisms in GPs and SPNs. Thereby, each GP is -- similarly as in a mixture of experts approach -- responsible only for a subset of data points, which effectively reduces inference cost in a divide and conquer fashion. We show that integrating GPs into the SPN framework leads to a promising probabilistic regression model which is: (1) computational and memory efficient, (2) allows efficient and exact posterior inference, (3) is flexible enough to mix different kernel functions, and (4) naturally accounts for non-stationarities in time series. In a variate of experiments, we show that the SPN-GP model can learn input dependent parameters and hyper-parameters and is on par with or outperforms the traditional GPs as well as state of the art approximations on real-world data.


title: Mean-Field Leader-Follower Games with Terminal State Constraint

abstract: We analyze linear McKean-Vlasov forward-backward SDEs arising in leader-follower games with mean-field type control and terminal state constraints on the state process. We establish an existence and uniqueness of solutions result for such systems in time-weighted spaces as well as a {convergence} result of the solutions with respect to certain perturbations of the drivers of both the forward and the backward component. The general results are used to solve a novel single-player model of portfolio liquidation under market impact with expectations feedback as well as a novel Stackelberg game of optimal portfolio liquidation with asymmetrically informed players.


title: Torus orbifolds with two fixed points

abstract: The main objects of this paper are torus orbifolds that have exactly two fixed points. We study the equivariant topological type of these orbifolds and consider when we can use the results of the paper [DKS] (arXiv:1809.03678) to compute its integral equivariant cohomology, in terms of generators and relations, coming from the corresponding orbifold torus graph.


title: Label Denoising with Large Ensembles of Heterogeneous Neural Networks

abstract: Despite recent advances in computer vision based on various convolutional architectures, video understanding remains an important challenge. In this work, we present and discuss a top solution for the large-scale video classification (labeling) problem introduced as a Kaggle competition based on the YouTube-8M dataset. We show and compare different approaches to preprocessing, data augmentation, model architectures, and model combination. Our final model is based on a large ensemble of video- and frame-level models but fits into rather limiting hardware constraints. We apply an approach based on knowledge distillation to deal with noisy labels in the original dataset and the recently developed mixup technique to improve the basic models.


title: Reconstruction of the Real Quantum Channel via Convex Optimization

abstract: Quantum process tomography is often used to completely characterize an unknown quantum process. However, it may lead to an unphysical process matrix, which will cause the loss of information respect to the tomography result. Convex optimization, widely used in machine learning, is able to generate a global optimal model that best fits the raw data while keeping the process tomography in a legitimate region. Only by correctly revealing the original action of the process can we seek deeper into its properties like its phase transition and its Hamiltonian. Thus, we reconstruct the real quantum channel using convex optimization from our experimental result obtained in free-space seawater. In addition, we also put forward a criteria, state deviation, to evaluate how well the reconstructed process fits the tomography result. We believe that the crossover between quantum process tomography and convex optimization may help us move forward to machine learning of quantum channels.


title: High dimensional measurement device independent quantum key distribution on two dimensional subspaces

abstract: Quantum key distribution (QKD) provides ultimate cryptographic security based on the laws of quantum mechanics. For point-to-point QKD protocols, the security of the generated key is compromised by detector side channel attacks. This problem can be solved with measurement device independent QKD (mdi-QKD). However, mdi-QKD has shown limited performances in terms of the secret key generation rate, due to post-selection in the Bell measurements. We show that high dimensional (Hi-D) encoding (qudits) improves the performance of current mdi-QKD implementations. The scheme is proven to be unconditionally secure even for weak coherent pulses with decoy states, while the secret key rate is derived in the single photon case. Our analysis includes phase errors, imperfect sources and dark counts to mimic real systems. Compared to the standard bidimensional case, we show an improvement in the key generation rate.


title: The impact of the surface energy coefficient on the deformation properties of atomic nuclei as predicted by Skyrme energy density functionals

abstract: In the framework of nuclear energy density functional (EDF) methods, many nuclear phenomena are related to the deformation of intrinsic states. Their accurate modeling relies on the correct description of the change of nuclear binding energy with deformation. The two most important contributions to the deformation energy have their origin in shell effects that are correlated to the spectrum of single-particle states, and the deformability of nuclear matter, that can be characterized by a model-dependent surface energy coefficient a_{surf}. With the goal of improving the global performance of nuclear EDFs through fine-tuning of their deformation properties, the purpose of this study is threefold. First, to analyze the impact of systematic variations of a_{surf} on properties of nuclei; second, to identify observables that can be safely used to narrow down the range of appropriate values of a_{surf} to be targeted in future parameter fits; third, to analyze the interdependence of a_{surf} with other properties of a nuclear EDF. Results for a large variety of relevant observables of deformed nuclei obtained from self-consistent mean-field calculations with a set of purpose-built SLy5sX parameterizations of the Skyrme EDF are correlated with the value of a_{surf}. The three main conclusions are that there is an evident preference for a comparatively low value of a_{surf}, as expected from the performance of existing parameterizations; that the isospin dependence of the surface energy also needs further finetuning in order to describe trends across the chart of nuclei; and that a satisfying simultaneous description of fission barriers and superdeformed states requires a better description of the single-particle spectra. [Note: The abstract has been abbreviated because of length restrictions imposed by the arXiv. See the paper for the full abstract.]


title: Meta-analysis of few studies involving rare events

abstract: Meta-analyses of clinical trials targeting rare events face particular challenges when the data lack adequate numbers of events for all treatment arms. Especially when the number of studies is low, standard meta-analysis methods can lead to serious distortions because of such data sparsity. To overcome this, we suggest the use of weakly informative priors (WIP) for the treatment effect parameter of a Bayesian meta-analysis model, which may also be seen as a form of penalization. As a data model, we use a binomial-normal hierarchical model (BNHM) which does not require continuity corrections in case of zero counts in one or both arms. We suggest a normal prior for the log odds ratio with mean 0 and standard deviation 2.82, which is motivated (1) as a symmetric prior centred around unity and constraining the odds ratio to within a range from 1/250 to 250 with 95 % probability, and (2) as consistent with empirically observed effect estimates from a set of $\mbox{$37\,773$}$ meta-analyses from the Cochrane Database of Systematic Reviews. In a simulation study with rare events and few studies, our BNHM with a WIP outperformed a Bayesian method without a WIP and a maximum likelihood estimator in terms of smaller bias and shorter interval estimates with similar coverage. Furthermore, the methods are illustrated by a systematic review in immunosuppression of rare safety events following paediatric transplantation. A publicly available $\textbf{R}$ package, $\texttt{MetaStan}$, is developed to automate the $\textbf{Stan}$ implementation of meta-analysis models using WIPs.


title: Gravitational instability caused by the weight of heat

abstract: Thermal energy points toward a disordered, completely uniform state acting counter to gravity's tendency to generate order and structure through gravitational collapse. It is, therefore, expected to contribute to the stabilization of a self-gravitating, classical ideal gas over collapse. However, I identified an instability that always occurs at sufficiently high energies: the high-energy or relativistic gravothermal instability. I argue here that this instability presents an analogous core--halo structure as its Newtonian counterpart, the Antonov instability. The main difference is that in the former case the core is dominated by the gravitation of thermal energy and not rest mass energy. A relativistic generalization of Antonov's instability---the low-energy gravothermal instability---also occurs. The two turning points, which make themselves evident as a double spiral of the caloric curve, approach each other as relativistic effects become more intense and eventually merge in a single point. Thus, the high and low-energy cases may be realized as two aspects of a single phenomenon---the gravothermal instability---which involves a core--halo separation and an intrinsic heat flow. Finally, I argue that the core formed during a core-collapse supernova is subject to the relativistic gravothermal instability if it becomes sufficiently hot and compactified at the time of the bounce. In this case, it will continue to collapse towards the formation of a black hole.

categories: gr-qc astro-ph.HE cond-mat.stat-mech hep-th


title: Thermal behavior, entanglement entropy and parton distributions

abstract: The apparent thermalization of the particles produced in hadronic collisions can be obtained by quantum entanglement of the partons of the initial state once a fast hard collision is produced. The scale of the hard collision is related to the thermal temperature. As the probability distribution of these events is of the form $np(n)$, as a consequence, the von Neumann entropy is larger than in the minimum bias case. The leading contribution to this entropy comes from the logarithm of the number of partons $n$, all with equal probability, making maximal the entropy. In addition there is another contribution related to the width of the parton multiplicity. Asymptotically, the entanglement entropy becomes the logarithm of $\sqrt{n}$, indicating that the number of microstates changes with energy from $n$ to $\sqrt{n}$.


title: Optimal residence time control for stochastically perturbed prescription opioid epidemic models

abstract: In this paper, we consider an optimal control problem for a prescription opioid epidemic model that describes the interaction between the regular prescription or addictive use of opioid drugs, and the process of rehabilitation and that of relapsing into opioid drug use. In particular, our interest is in the situation, where the control appearing linearly in the opioid epidemics is interpreted as the rate at which the susceptible individuals are effectively removed from the population due to an opioid-related intervention policy or when the dynamics of the addicted is strategically influenced due to an accessible addiction treatment facility, while a small perturbing noise enters through the dynamics of the susceptible group in the population compartmental model. To this end, we introduce a mathematical apparatus that minimizes the asymptotic exit-rate with which the solution for such stochastically perturbed prescription opioid epidemics exits from a given bounded open domain. Moreover, under certain assumptions, we also provide an admissible optimal Markov control for the corresponding optimal control problem that optimally effected removal of the susceptible or recovered individuals from the population dynamics.


title: From Molecular to Metallic Gold Nanoparticles: The Role of Nanocrystal Symmetry in the Crossover Region

abstract: Bulk gold is a good metal, i.e., conductor of electricity and heat, due to its delocalized electron density that can respond to extremely small external perturbations such as electric field or temperature gradient. In energy space, the quantum states of the conduction electrons cross over the metal's Fermi energy continuously. But when gold is dispersed in finite nanometer-size particles or 'clusters', the delocalized electronic states are re-grouped in energy space to 'shells' according to allowed energy levels and the symmetry and shape of the atomic arrangement. This re-grouping generates also an energy gap in the vicinity of the Fermi energy akin to the energy gap between occupied and unoccupied electron orbitals in molecules. How does the formation of these shells and the energy gap affect the physico-chemical properties of the nanoparticle? How does one get from a 'molecule-like' gold nanoparticle to 'metal-like' gold nanoparticle? Here we analyse the electronic structure and optical and chiroptical properties of recently reported gold nanoparticles of 144, 146, and 246 gold atoms, that are made by wet chemistry methods and whose structures have been resolved to atomic precision. We demonstrate computationally how re-grouping of the quantum states of valence electrons can affect drastically the optical properties of nanoparticles in the crossover-size region, by either generating a multi-band molecule-like or a monotonous metallic optical absorption. The lower the symmetry of the gold core, the more metallic is the nanoparticle. The underlying mechanism arises from symmetry-sensitive distribution of the electronic levels of the nanoparticle close to Fermi energy. Overall, this work sheds lights on fundamental mechanisms on how molecular nanoparticle properties can change metallic in the crossover-size region.


title: Chondrocyte Heterogeneity; It Is the Time to Update the Understanding of Cartilage Histology

abstract: Chondrocytes were described as one cell populations in most cartilage literature. Two different chondrocyte populations; dark and light, were described in the articular cartilage and a third population, adipochondrocytes, was described in the elastic cartilage. The cur-rent literature of cartilage histology should be updated and highlight that three different populations of chondrocytes are existed in cartilage.


title: Comments on "Amplifying magneto-optical photonic crystal, Appl. Phys. Lett. 97, 061116 (2010)"

abstract: This report concerns analytical solutions predicting the complex refractive indices of doped crystals, based on empirical values of the refractive indices of the dopant and the undoped crystal. An analytical approach in calculating the refractive index of Erbium doped Gadolinium Gallium Garnet was reported by another author previously. My first comment is related to the proposed solution by the other author, where I have shown an alternative approach that leads to a different solution. Although my comments are made in relation to Erbium as the dopant and Gadolinium Gallium Garnet as the undoped lattice, the analytical solution I propose here may be applied to any doped crystals. My second comment concerns the limitations and mistreatment of McCumber relation in this context.


title: The Moon Illusion explained by the Projective Consciousness Model

abstract: The Moon often appears larger near the perceptual horizon and smaller high in the sky though the visual angle subtended is invariant. We show how this illusion results from the optimization of a projective geometrical frame for conscious perception through free energy minimization, as articulated in the Projective Consciousness Model. The model accounts for all documented modulations of the illusion without anomalies (e.g., the size-distance paradox), surpasses other theories in explanatory power, makes sense of inter- and intra-subjective variability vis-a-vis the illusion, and yields new quantitative and qualitative predictions.


title: Rethinking Location Privacy for Unknown Mobility Behaviors

abstract: Location Privacy-Preserving Mechanisms (LPPMs) in the literature largely consider that users' data available for training wholly characterizes their mobility patterns. Thus, they hardwire this information in their designs and evaluate their privacy properties with these same data. In this paper, we aim to understand the impact of this decision on the level of privacy these LPPMs may offer in real life when the users' mobility data may be different from the data used in the design phase. Our results show that, in many cases, training data does not capture users' behavior accurately and, thus, the level of privacy provided by the LPPM is often overestimated. To address this gap between theory and practice, we propose to use blank-slate models for LPPM design. Contrary to the hardwired approach, that assumes known users' behavior, blank-slate models learn the users' behavior from the queries to the service provider. We leverage this blank-slate approach to develop a new family of LPPMs, that we call Profile Estimation-Based LPPMs. Using real data, we empirically show that our proposal outperforms optimal state-of-the-art mechanisms designed on sporadic hardwired models. On non-sporadic location privacy scenarios, our method is only better if the usage of the location privacy service is not continuous. It is our hope that eliminating the need to bootstrap the mechanisms with training data and ensuring that the mechanisms are lightweight and easy to compute help fostering the integration of location privacy protections in deployed systems.