title: Fossil stellar streams and their globular cluster populations in the E-MOSAICS simulations

abstract: Stellar haloes encode a fossil record of a galaxy's accretion history, generally in the form of structures of low surface brightness, such as stellar streams. While their low surface brightness makes it challenging to determine their age, metallicity, kinematics and spatial structure, the infalling galaxies also deposit globular clusters (GCs) in the halo, which are bright and therefore easier to observe and characterise. To understand how GCs associated with stellar streams can be used to estimate the stellar mass and the infall time of their parent galaxy, we examine a subset of 15 simulations of galaxies and their star clusters from the E-MOSAICS project. E-MOSAICS is a suite of hydrodynamical simulations incorporating a sub-grid model for GC formation and evolution. We find that more massive accreted galaxies typically contribute younger and more metal rich GCs. This lower age results from a more extended cluster formation history in more massive galaxies. In addition, at fixed stellar mass, galaxies that are accreted later host younger clusters, because they can continue to form GCs without being subjected to environmental influences for longer. This explains the large range of ages observed for clusters associated with the Sagittarius dwarf galaxy in the halo of the Milky Way compared to clusters which are thought to have formed in satellites accreted early in the Milky Way's formation history. Using the ages of the GCs associated with the Sagittarius dwarf, we estimate a virial radius crossing lookback time (infall time) of $9.3 \pm 1.8 Gyr$.


title: Magnonic Analogue of Black/White Hole Horizon in Superfluid $^3$He-B

abstract: We report on theoretical model and experimental results of the experiment made in a limit of absolute zero temperature ($\sim$ 600\,$\mu$K) studying the spin wave analogue of black/white hole horizon using spin (magnonic) superfluidity in superfluid $^3$He-B. As an experimental tool simulating the properties of the black/white horizon we used the spin-precession waves propagating on the background of the spin super-currents between two Bose-Einstein condensates of magnons in form of homogeneously precessing domains. We provide experimental evidence of the white hole formation for spin precession waves in this system, together with observation of an amplification effect. Moreover, the estimated temperature of the spontaneous Hawking radiation in this system is about four orders of magnitude lower than the system's background temperature what makes it a promising tool to study the effect of spontaneous Hawking radiation.


title: Nonequispaced Fast Fourier Transform (NFFT) Interface for Julia

abstract: This report describes the newly added Julia interface to the NFFT3 library. We explain the multidimensional NFFT algorithm and basics of the interface. Furthermore, we go into detail about the different parameters and how to adjust them properly.


title: A Thermodynamical derivation of the quantum potential and the temperature of the wave function

abstract: In this paper a thermodynamical derivation of the quantum potential is pro- posed. Within the framework of Bohmian mechanics we show how the quantum potential can be derived, by adding an additional informational degree of freedom to the ordinary degrees of freedom of a physical system. Such a derivation uses the First Law of thermodynamics for this additional degree of freedom and basic equilibrium thermodynamics methods. By doing that, one may associate a temper- ature to each wave function. Features and behavior of this temperature in different situations is studied.


title: On nonsingularity of circulant matrices

abstract: In Communication theory and Coding, it is expected that certain circulant matrices having $k$ ones and $k+1$ zeros in the first row are nonsingular. We prove that such matrices are always nonsingular when $2k+1$ is either a power of a prime, or a product of two distinct primes. For any other integer $2k+1$ we construct circulant matrices having determinant $0$. The smallest singular matrix appears when $2k+1=45$. The possibility for such matrices to be singular is rather low, smaller than $10^{-4}$ in this case.


title: Heterogeneous large datasets integration using Bayesian factor regression

abstract: Two key challenges in modern statistical applications are the large amount of information recorded per individual, and that such data are often not collected all at once but in batches. These batch effects can be complex, causing distortions in both mean and variance. We propose a novel sparse latent factor regression model to integrate such heterogeneous data. The model provides a tool for data exploration via dimensionality reduction while correcting for a range of batch effects. We study the use of several sparse priors (local and non-local) to learn the dimension of the latent factors. Our model is fitted in a deterministic fashion by means of an EM algorithm for which we derive closed-form updates, contributing a novel scalable algorithm for non-local priors of interest beyond the immediate scope of this paper. We present several examples, with a focus on bioinformatics applications. Our results show an increase in the accuracy of the dimensionality reduction, with non-local priors substantially improving the reconstruction of factor cardinality, as well as the need to account for batch effects to obtain reliable results. Our model provides a novel approach to latent factor regression that balances sparsity with sensitivity and is highly computationally efficient.


title: Spin-locking metasurface for surface plasmon routing

abstract: The spin-dependent routing of surface plasmons on metal by use of mirror symmetric periodic metasurface is presented and analyzed. We incorporate the intrinsic transverse spin angular momentum of the plasmonic wave in order to efficiently couple the incident light carrying circular polarization to a desired direction. The combination of the oblique incidence of polarized light with the accurately designed momentum matching properties of the grating provides a perfect way to achieve a spin-locking metasurface that routes the plasmonic beams.


title: Neutrino recoil force in electron-capture decay of polarized nuclei: measurement prospects and potential applications

abstract: Due to a directional asymmetry of neutrino emission caused by parity violation, a sample of radioactive atoms experiences a small recoil force from neutrino radiation accompanying electron capture by polarized nuclei. An expression for this force is obtained for the case of allowed nuclear transitions. Prospects to measure this force by the use of modern micromechanical devices are considered. Numerical estimates for the force are presented for a number of most suitable radioactive isotopes. Potential applications for the weak interaction studies are discussed including the possibility to search for hypothetical Lorentz invariance violation.


title: Young massive clusters in the interacting LIRG Arp 299: evidence for the dependence of star cluster formation and evolution on environment

abstract: Archival WFC3/UVIS imaging of Arp 299 (NGC 3690E + NGC 3690W) is retrieved to investigate the young massive cluster (YMC) population of this ongoing merger. We extract 2182 cluster candidates, including 1323 high confidence photometric sources. Multiband photometry is matched with Yggdrasil models to estimate the age, mass, and extinction of each cluster. A Schechter fit of the truncated cluster mass function results in a characteristic mass ${\rm M_{*} = 1.6 \times 10^6 M_{\odot}}$. Our results confirm that intensely star-forming galaxies such as Arp 299 host more massive clusters than quiescent dwarf and normal spirals. In the case of NGC 3690E, we find that the cluster masses decrease with an increasing galactocentric radius likely due to the gas density distribution. On the other hand, the fitted age distributions of a mass-limited sample suggest that YMCs of the western component undergo stronger disruption than those hosted by the eastern galaxy. This is in agreement with the properties of the underlying cluster luminosity functions: a clear truncation at high luminosities with slopes generally shallower by $\sim 0.3$ dex than the ones of the NGC 3690E. Finally, the derived cluster formation efficiency, $\Gamma \sim 19$ percent, indicates that Arp 299 has $\sim 3-5$ times more star formation happening in bound clusters compared to the cases of gas-poor spirals like NGC 2997 and NGC 4395. The merger generally follows the $\Gamma - $ star formation rate density relation from the literature. The YMC photometric study of Arp 299 thus reveals that both formation and disruption mechanisms of the star cluster population are most likely environment-dependent.


title: Exploiting the Space Filling Curve Ordering of Particles in the Neighbour Search of Gadget3

abstract: Gadget3 is nowadays one of the most frequently used high performing parallel codes for cosmological hydrodynamical simulations. Recent analyses have shown t\ hat the Neighbour Search process of Gadget3 is one of the most time-consuming parts. Thus, a considerable speedup can be expected from improvements of the u\ nderlying algorithms. In this work we propose a novel approach for speeding up the Neighbour Search which takes advantage of the space-filling-curve particle ordering. Instead of performing Neighbour Search for all particles individually, nearby active particles can be grouped and one single Neighbour Search can be performed to obta\ in a common superset of neighbours. Thus, with this approach we reduce the number of searches. On the other hand, tree walks are performed within a larger searching radius. There is an optimal size of grouping that maximize the speedup, which we found by numerical experiments. We tested the algorithm within the boxes of the Magneticum project. As a result we obtained a speedup of $1.65$ in the Density and of $1.30$ in the Hydrodynamics computation, respectively, and a total speedup of $1.34.$


title: Dynamic Likelihood-free Inference via Ratio Estimation (DIRE)

abstract: Parametric statistical models that are implicitly defined in terms of a stochastic data generating process are used in a wide range of scientific disciplines because they enable accurate modeling. However, learning the parameters from observed data is generally very difficult because their likelihood function is typically intractable. Likelihood-free Bayesian inference methods have been proposed which include the frameworks of approximate Bayesian computation (ABC), synthetic likelihood, and its recent generalization that performs likelihood-free inference by ratio estimation (LFIRE). A major difficulty in all these methods is choosing summary statistics that reduce the dimensionality of the data to facilitate inference. While several methods for choosing summary statistics have been proposed for ABC, the literature for synthetic likelihood and LFIRE is very thin to date. We here address this gap in the literature, focusing on the important special case of time-series models. We show that convolutional neural networks trained to predict the input parameters from the data provide suitable summary statistics for LFIRE. On a wide range of time-series models, a single neural network architecture produced equally or more accurate posteriors than alternative methods.


title: From 2D to 3D: graphene moulding for transparent and flexible probes

abstract: Chemical vapor deposition (CVD) has been widely adopted as the most scalable method to obtain single layer graphene. Incorporating CVD graphene in planar devices can be performed via well established wet transfer methods or thermal adhesive release. Nevertheless, for applications involving 3D shapes, methods adopted for planar surface provide only a crude solution if a continuous, tear-free, wrinkle-free graphene layer is required. In this work, we present the fabrication and characterization of PDMS-supported 3D graphene probes. To accommodate 3D geometries, we perform CVD on catalysts possessing a non-trivial 3D topology, serving to mold the grown graphene to a final non-trivial 3D shape. This advance overcomes challenges observed in standard transfer processes that can result in uneven coverage, wrinkles and tears. To demonstrate the potential of our new transfer approach, we apply it to fabricate graphene electrical probes. Graphene, due to its flexibility, transparency and conductivity, is an ideal material with which to replace conventional metal based probes. In particular, with a contact resistance in the order of tens of kOhm, these graphene probes may find applications, such as in electrophysiology studies.


title: Calibration of the DAMPE Plastic Scintillator Detector and its on-orbit performance

abstract: DArk Matter Particle Explorer (DAMPE) is a space-borne apparatus for detecting the high-energy cosmic-rays like electrons, $\gamma$-rays, protons and heavy-ions. Plastic Scintillator Detector (PSD) is the top-most sub-detector of the DAMPE. The PSD is designed to measure the charge of incident high-energy particles and it also serves as a veto detector for discriminating $\gamma$-rays from charged particles. In this paper, PSD on-orbit calibration procedure is described, which includes five steps of pedestal, dynode correlation, response to minimum-ionizing particles (MIPs), light attenuation function and energy reconstruction. A method for reconstructing the charge of incident high energy cosmic-ray particles is introduced. The detection efficiency of each PSD strip is verified to be above 99.5%, the total efficiency of the PSD for charged particles is above 99.99%.


title: Ultra-low Hysteresis in Giant Magnetocaloric Mn1-xVxFe(P,Si,B) Compounds

abstract: Large thermal hysteresis in the MnFe(P, Si, B) system hinders the heat exchange rate and thus limits the magnetocaloric applications at high frequencies. Substitution of Mn by V in Mn1-xVxFe0.95P0.593Si0.33B0.077 and Mn1-xVxFe0.95P0.563Si0.36B0.077 alloys was found to reduce the thermal hysteresis due to a decrease in the latent heat. Introducing V increases both the field-induced transition temperature shift and the magnetic moment per formula unit. Thus, a decease in the thermal hysteresis is obtained without losing the giant magnetocaloric effect. In consequence, an ultralow hysteresis (0.7 K) and a giant adiabatic temperature change of 2.3 K were achieved, which makes these alloys promising candidates for commercial magnetic refrigerator using permanent magnets.


title: The evolved fast rotator Sargas. Stellar parameters and evolutionary status from VLTI/PIONIER and VLT/UVES

abstract: Gravity darkening (GD) and flattening are important consequences of stellar rotation. The precise characterization of these effects across the HRD is crucial to a deeper understanding of stellar structure and evolution. We seek to characterize such important effects on Sargas, an evolved, fast-rotating, intermediate-mass star, located in a region of the HRD where they have never been directly measured as far as we know. We use our numerical model CHARRON to analyze interferometric (VLTI/PIONIER) and spectroscopic (VLT/UVES) observations through a MCMC model-fitting procedure. The visibilities and closure phases from the PIONIER data are particularly sensitive to rotational flattening and GD. Adopting the Roche approximation, we investigate two GD models: (1) the beta-model (classical von Zeipel's law), and (2) the omega-model. Using this approach we measure several physical parameters of Sargas, namely, equatorial radius, mass, equatorial rotation velocity, mean Teff, inclination and position angle of the rotation axis, and beta. In particular, we show that the measured beta leads to a surface flux distribution equivalent to the one given by the omega-model. Thanks to our results, we also show that Sargas is most probably located in a rare and interesting region of the H-R diagram: within the Hertzsprung gap and over the hot edge of the instability strip. These results show once more the power of optical/IR long-baseline interferometry, combined with high-resolution spectroscopy, to directly measure fast-rotation effects and stellar parameters, in particular GD. As was the case for a few fast rotators previously studied by interferometry, the omega-model provides a physically more profound description of Sargas' GD, without the need of a beta exponent.


title: Singular robustly chain transitive sets are singular volume partial hyperbolic

abstract: For diffeomorphisms or for non-singular flows, there are many results relating properties persistent under C1 perturbations and global structures for the dynamics ( such as hyperbolicity, partial hyperbolicity, dominated splitting). However, a dif\'iculty appears when a robust property of a flow holds on a set containing recurrent orbits accumulating a singular point. In [BdL] with Christan Bonatti we propose a a general procedure for adapting the usual hyperbolic structures to the singularities. In this paper, using this tool, we recover the results in [BDP] for flows, showing that robustly chain transitive sets have a weak form of hyperbolicity. allowing us to conclude as well the kind of hyperbolicity carried by the examples in \cite{BLY} (a robust chain transitive singular attractor with periodic orbits of different indexes). Along with the results in [BdL], this shows that the way we propose to interpret the effect of singularities has the potential to adapt to other settings in which there is a coexistence of singularities and regular orbits with the goal of re-obtaining the results that we already know for diffeomorphisms.


title: A Defect-Free Naphthalene Diimide Bithiazole Copolymer via Regioselective Direct Arylation Polycondensation

abstract: An optimized direct arylation polycondensation (DAP) protocol for the synthesis of a novel naphthalene diimide (NDI) 2,2'-bithiazole (2-BTz) copolymer (PNDI-2-BTz) is presented. The regioselective C-H activation of 2-BTz at the 5-positions allows for the synthesis of fully regioregular and homocoupling-free PNDI-2-BTz of high molecular weight in less than 1 h in quantitative yield. Complete end group assignment shows functionalities according to monomer structures or to nucleophilic substitution, and allows for the reliable determination of absolute molecular weight. Compared to the well-known bithiophene analog PNDIT2, an exceptionally high thermal stability, a hypsochromically shifted charge transfer absorption band and a lower-lying LUMO energy level is found, making PNDI-2-BTz an interesting candidate for applications in organic electronic devices. In contrast to the selective and high yielding C-H activation of 2-BTz at the 5-position, the regioisomer 5,5'-bithiazole is inactive under a variety of conditions.


title: Approximating Quasiparticle and Excitation Energies from Ground State Generalized Kohn-Sham Calculations

abstract: Quasiparticle energies and fundamental band gaps in particular are critical properties of molecules and materials. It was rigorously established that the generalized Kohn-Sham HOMO and LUMO orbital energies are the chemical potentials of electron removal and addition and thus good approximations to band edges and fundamental gaps from a density functional approximation (DFA) with minimal delocalization error. For other quasiparticle energies, their connection to the generalized Kohn-Sham orbital energies has not been established but remains highly interesting. We provide the comparison of experimental quasiparticle energies for many finite systems with calculations from the GW Green's function and localized orbitals scaling correction (LOSC), a recently developed correction to semilocal DFAs, which has minimal delocalization error. Extensive results with over forty systems clearly show that LOSC orbital energies achieve slightly better accuracy than the GW calculations with little dependence on the semilocal DFA, supporting the use of LOSC DFA orbital energies to predict quasiparticle energies. This also leads to the calculations of excitation energies of the $N$-electron systems from the ground state DFA calculations of the $\left(N-1\right)$-electron systems. Results show good performance with accuracy similar to TDDFT and the delta SCF approach for valence excitations with commonly used DFAs with or without LOSC. For Rydberg states, good accuracy was obtained only with the use of LOSC DFA. This work highlights the pathway to quasiparticle and excitation energies from ground density functional calculations.


title: Linear Receivers in Non-stationary Massive MIMO Channels with Visibility Regions

abstract: In a massive MIMO system with large arrays, the channel becomes spatially non-stationary. We study the impact of spatial non-stationarity characterized by visibility regions (VRs) where the channel energy is significant on a portion of the array. Relying on a channel model based on VRs, we provide expressions of the signal-to-interference-plus-noise ratio (SINR) of conjugate beamforming (CB) and zero-forcing (ZF) precoders. We also provide an approximate deterministic equivalent of the SINR of ZF precoders. We identify favorable and unfavorable multi-user configurations of the VRs and compare the performance of both stationary and non-stationary channels through analysis and numerical simulations.


title: Interaction of Skyrmions and Pearl Vortices in Superconductor-Chiral Ferromagnet Heterostructures

abstract: We investigate a hybrid heterostructure with magnetic skyrmions (Sk) inside a chiral ferromagnet interfaced by a thin superconducting film via an insulating barrier. The barrier prevents the electronic transport between the superconductor and the chiral magnet, such that the coupling can only occur through the magnetic fields generated by these materials. We find that Pearl vortices (PV) are generated spontaneously in the superconductor within the skyrmion radius, while anti-Pearl vortices (\(\overline{\mathrm{PV}}\)) compensating the magnetic moment of the Pearl vortices are generated outside of the Sk radius, forming an energetically stable topological hybrid structure. Finally, we analyze the interplay of skyrmion and vortex lattices and their mutual feedback on each other. In particular, we argue that the size of the skyrmions will be greatly affected by the presence of the vortices offering another prospect of manipulating the skyrmionic size by the proximity to a superconductor.


title: Bayes Factor Asymptotics for Variable Selection in the Gaussian Process Framework

abstract: Although variable selection is one of the most popular areas of modern statistical research, much of its development has taken place in the classical paradigm compared to the Bayesian counterpart. Somewhat surprisingly, both the paradigms have focussed almost completely on linear models, in spite of the vast scope offered by the model liberation movement brought about by modern advancements in studying real, complex phenomena. In this article, we investigate general Bayesian variable selection in models driven by Gaussian processes, which allows us to treat linear, non-linear and nonparametric models, in conjunction with even dependent setups, in the same vein. We consider the Bayes factor route to variable selection, and develop a general asymptotic theory for the Gaussian process framework in the "large p, large n" settings even with p>>n, establishing almost sure exponential convergence of the Bayes factor under appropriately mild conditions. The fixed p setup is included as a special case. To illustrate, we apply our general result to variable selection in linear regression, Gaussian process model with squared exponential covariance function accommodating the covariates, and a first order autoregressive process with time-varying covariates. We also follow up our theoretical investigations with ample simulation experiments in the above regression contexts and variable selection in a real, riboflavin data consisting of 71 observations but 4088 covariates. For implementation of variable selection using Bayes factors, we develop a novel and effective general-purpose transdimensional, transformation based Markov chain Monte Carlo algorithm, which has played a crucial role in our simulated and real data applications.


title: Mechanical factors affecting the mobility of membrane proteins

abstract: The mobility of membrane proteins controls many biological functions. The application of the model of Saffman and Delbr\"uck to the diffusion of membrane proteins does not account for all the experimental measurements. These discrepancies have triggered a lot of studies on the role of the mechanical factors in the mobility. After a short review of the Saffman and Delbr\"uck model and of some key experiments, we explore the various ways to incorporate the effects of the different mechanical factors. Our approach focuses on the coupling of the protein to the membrane, which is the central element in the modelling. We present a general, polaron-like model, its recent application to the mobility of a curvature sensitive protein, and its various extensions to other couplings that may be relevant in future experiments.


title: Interpreting Frame Transformations as Diagonalization of Harmonic Transfer Functions

abstract: Analysis of ac electrical systems can be performed via frame transformations in the time-domain or via harmonic transfer functions (HTFs) in the frequency-domain. The two approaches each have unique advantages but are hard to reconcile because the coupling effect in the frequency-domain leads to infinite dimensional HTF matrices that need to be truncated. This paper explores the relation between the two representations and shows that applying a similarity transformation to an HTF matrix creates a direct equivalence to a frame transformation on the input-output signals. Under certain conditions, such similarity transformations have a diagonalizing effect which, essentially, reduces the HTF matrix order from infinity to two or one, making the matrix tractable mathematically without truncation or approximation. This theory is applied to a droop-controlled voltage source inverter as an illustrative example. A stability criterion is derived in the frequency-domain which agrees with the conventional state-space model but offers greater insights into the mechanism of instability in terms of the negative damping (non-passivity) under droop control. The paper not only establishes a unified view in theory but also offers an effective practical tool for stability assessment.


title: Efficient Bayesian Experimental Design for Implicit Models

abstract: Bayesian experimental design involves the optimal allocation of resources in an experiment, with the aim of optimising cost and performance. For implicit models, where the likelihood is intractable but sampling from the model is possible, this task is particularly difficult and therefore largely unexplored. This is mainly due to technical difficulties associated with approximating posterior distributions and utility functions. We devise a novel experimental design framework for implicit models that improves upon previous work in two ways. First, we use the mutual information between parameters and data as the utility function, which has previously not been feasible. We achieve this by utilising Likelihood-Free Inference by Ratio Estimation (LFIRE) to approximate posterior distributions, instead of the traditional approximate Bayesian computation or synthetic likelihood methods. Secondly, we use Bayesian optimisation in order to solve the optimal design problem, as opposed to the typically used grid search or sampling-based methods. We find that this increases efficiency and allows us to consider higher design dimensions.


title: Most Probable Evolution Trajectories in a Genetic Regulatory System Excited by Stable L\'evy Noise

abstract: We study the most probable trajectories of the concentration evolution for the transcription factor activator in a genetic regulation system, with non-Gaussian stable L\'evy noise in the synthesis reaction rate taking into account. We calculate the most probable trajectory by spatially maximizing the probability density of the system path, i.e., the solution of the associated nonlocal Fokker-Planck equation. We especially examine those most probable trajectories from low concentration state to high concentration state (i.e., the likely transcription regime) for certain parameters, in order to gain insights into the transcription processes and the tipping time for the transcription likely to occur. This enables us: (i) to visualize the progress of concentration evolution (i.e., observe whether the system enters the transcription regime within a given time period); (ii) to predict or avoid certain transcriptions via selecting specific noise parameters in particular regions in the parameter space. Moreover, we have found some peculiar or counter-intuitive phenomena in this gene model system, including (a) a smaller noise intensity may trigger the transcription process, while a larger noise intensity can not, under the same asymmetric L\'evy noise. This phenomenon does not occur in the case of symmetric L\'evy noise; (b) the symmetric L\'evy motion always induces transition to high concentration, but certain asymmetric L\'evy motions do not trigger the switch to transcription. These findings provide insights for further experimental research, in order to achieve or to avoid specific gene transcriptions, with possible relevance for medical advances.


title: Suggestion of Coherent Radio Reflections from an Electron-Beam Induced Particle Cascade

abstract: Testbeam experiment 576 (T576) at the SLAC National Accelerator Laboratory sought to make the first measurement of coherent radio reflections from the ionization produced in the wake of a high-energy particle shower. The >10 GeV electron beam at SLAC End Station A was directed into a large high-density polyethylene target to produce a shower analogous to that produced by an EeV neutrino interaction in ice. Continuous wave (CW) radio was transmitted into the target, and receiving antennas monitored for reflection of the transmitted signal from the ionization left in the wake of the shower. We detail the experiment and report on preliminary hints of a signal consistent with a radio reflection at a significance of 2.36$\sigma$. We recommend another test-beam measurement in order to verify the signal.


title: Spiderweb central configurations

abstract: In this paper we study spiderweb central configurations for the $N$-body problem, i.e configurations given by $N=n \times \ell+1$ masses located at the intersection points of $\ell$ concurrent equidistributed half-lines with $n$ circles and a central mass $m_0$, under the hypothesis that the $\ell$ masses on the $i$-th circle are equal to a positive constant $m_i$; we allow the particular case $m_0=0$. We focus on constructive proofs of the existence of spiderweb central configurations, which allow numerical implementation. Additionally, we prove the uniqueness of such central configurations when $\ell \in \{2,\dots,9\}$ and arbitrary $n$ and $m_i$; under the constraint $m_1\geq m_2\geq \ldots \geq m_n$ we also prove uniqueness for $\ell \in \{10,\dots,18\}$ and $n$ not too large. We also give an algorithm providing a rigorous proof of the existence and local unicity of such central configurations when given as input a choice of $n$, $\ell$ and $m_0, . . . ,m_n$. Finally, our numerical simulations highlight some interesting properties of the mass distribution.


title: On Fractional Annealing Process

abstract: The paper deals with a full annealing process, for example for a metal alloy perturbed by some fractional noise caused by conditions in furnace such as initial temperature, initial structure of material itsel. The energy of the heating matter system depending on its states at each time, we consider some steady state at full annealing range by introducing a fractional model for it. And we determine states of the system related to this steady state.


title: Dwarf spheroidal $J$-factor likelihoods for generalized NFW profiles

abstract: Indirect detection strategies of particle Dark Matter (DM) in Dwarf spheroidal satellite galaxies (dSphs) typically entail searching for annihilation signals above the astrophysical background. To robustly compare model predictions with the observed fluxes of product particles, most analyses of astrophysical data -- which are generally frequentist -- rely on estimating the abundance of DM by calculating the so-called $\textit{J-factor}$. This quantity is usually inferred from the kinematic properties of the stellar population of a dSph using Jeans equation, commonly by means of Bayesian techniques which entail the presence (and additional systematic uncertainty) of prior choice. Here, extending earlier work, we develop a scheme to derive the profile likelihood for $J$-factors of dwarf spheroidals for models with five or more free parameters. We validate our method on a publicly available simulation suite, released by the Gaia Challenge, finding satisfactory statistical properties for coverage and bias. We present the profile likelihood function and maximum likelihood estimates for the $J$-factor of ten dSphs . As an illustration, we apply these profiles likelihood to recently published analyses of gamma-ray data with the Fermi Large Area Telescope to derive new, consistent upper limits on the DM annihilation cross-section. We do this for a subset of systems, generally referred to as $\textit{classical dwarfs}$. The implications of these findings for DM searches are discussed, together with future improvements and extensions of this technique.

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


title: Ionic transport through a protein nanopore: a Coarse-Grained Molecular Dynamics Study

abstract: The MARTINI coarse-grained (CG) force field is used to test the ability of CG models to simulate ionic transport through protein nanopores. The ionic conductivity of CG ions in solution was computed and compared with experimental results. Next, we studied the electrostatic behavior of a solvated CG lipid bilayer in salt solution under an external electric field. We showed this approach correctly describes the experimental conditions under a potential bias. Finally, we performed CG molecular dynamics simulations of the ionic transport through a protein nanopore ($\alpha$-hemolysin) inserted in a lipid bilayer, under different electric fields, for 2-3 microseconds. The resulting $I-V$ curve is qualitatively consistent with experiments, although the computed current is one order of magnitude smaller. Current saturation was observed for potential biases over $\pm~350$~mV. We also discuss the time to reach a stationary regime and the role of the protein flexibility in our CG simulations.


title: Measuring Bilayer Surface Energy and Curvature in Asymmetric Droplet Interface Bilayers

abstract: For the past decade, droplet interface bilayers (DIBs) have had an increased prevalence in biomolecular and biophysical literature. However, much of the underlying physics of these platforms are poorly characterized. To further the understanding of these structures, a study of lipid membrane tension on DIB membranes is measured by analysing the equilibrium shape of asymmetric DIBs. To this end, the morphology of DIBs is explored for the first time using confocal laser scanning fluorescence microscopy (CLSM). The experimental results confirm that, in accordance with theory, the bilayer interface of a volume asymmetric DIB is curved toward the smaller droplet and a lipid asymmetric DIB is curved toward the droplet with the higher monolayer surface tension. Moreover, the DIB shape can be exploited to measure complex bilayer surface energies. In this study, the bilayer surface energy of DIBs composed of lipid mixtures of phosphatidylgylcerol (PG) and phosphatidylcholine (PC) are shown to increase linearly with PG concentrations up to 25%. The assumption that DIB bilayer area can be geometrically approximated as a spherical cap base is also tested, and it is discovered that the bilayer curvature is negligible for most practical symmetric or asymmetric DIB systems with respect to bilayer area.


title: Clustering Time Series with Nonlinear Dynamics: A Bayesian Non-Parametric and Particle-Based Approach

abstract: We propose a general statistical framework for clustering multiple time series that exhibit nonlinear dynamics into an a-priori-unknown number of sub-groups. Our motivation comes from neuroscience, where an important problem is to identify, within a large assembly of neurons, subsets that respond similarly to a stimulus or contingency. Upon modeling the multiple time series as the output of a Dirichlet process mixture of nonlinear state-space models, we derive a Metropolis-within-Gibbs algorithm for full Bayesian inference that alternates between sampling cluster assignments and sampling parameter values that form the basis of the clustering. The Metropolis step employs recent innovations in particle-based methods. We apply the framework to clustering time series acquired from the prefrontal cortex of mice in an experiment designed to characterize the neural underpinnings of fear.


title: Connectivity of inhomogeneous random K-out graphs

abstract: We propose inhomogeneous random K-out graphs $\mathbb{H}(n; \pmb{\mu}, \pmb{K}_n)$, where each of the $n$ nodes is assigned to one of $r$ classes independently with a probability distribution $\pmb{\mu} = \{\mu_1, \ldots, \mu_r\}$. In particular, each node is classified as class-$i$ with probability $\mu_i>0$, independently. Each class-$i$ node selects $K_{i,n}$ distinct nodes uniformly at random from among all other nodes. A pair of nodes are adjacent in $\mathbb{H}(n; \pmb{\mu}, \pmb{K}_n)$ if at least one selects the other. Without loss of generality, we assume that $K_{1,n} \leq K_{2,n} \leq \ldots \leq K_{r,n}$. Earlier results on homogeneous random K-out graphs $\mathbb{H}(n; K_n)$, where all nodes select the same number $K$ of other nodes, reveal that $\mathbb{H}(n; K_n)$ is connected with high probability (whp) if $K_n \geq 2$ which implies that $\mathbb{H}(n; \pmb{\mu}, \pmb{K}_n)$ is connected whp if $K_{1,n} \geq 2$. In this paper, we investigate the connectivity of inhomogeneous random K-out graphs $\mathbb{H}(n; \pmb{\mu}, \pmb{K}_n)$ for the special case when $K_{1,n}=1$, i.e., when each class-$1$ node selects only one other node. We show that $\mathbb{H}\left(n;\pmb{\mu},\pmb{K}_n\right)$ is connected whp if $K_{r,n}$ is chosen such that $\lim_{n \to \infty} K_{r,n} = \infty$. However, any bounded choice of the sequence $K_{r,n}$ gives a positive probability of $\mathbb{H}\left(n;\pmb{\mu},\pmb{K}_n\right)$ being not connected. Simulation results are provided to validate our results in the finite node regime.


title: Non-i.i.d. random holomorphic dynamical systems and the probability of tending to infinity

abstract: We consider random holomorphic dynamical systems on the Riemann sphere whose choices of maps are related to Markov chains. Our motivation is to generalize the facts which hold in i.i.d. random holomorphic dynamical systems. In particular, we focus on the function $T$ which represents the probability of tending to infinity. We show some sufficient conditions which make $T$ continuous on the whole space and we characterize the Julia sets in terms of the function $T$ under certain assumptions.


title: Learning Classical Planning Strategies with Policy Gradient

abstract: A common paradigm in classical planning is heuristic forward search. Forward search planners often rely on simple best-first search which remains fixed throughout the search process. In this paper, we introduce a novel search framework capable of alternating between several forward search approaches while solving a particular planning problem. Selection of the approach is performed using a trainable stochastic policy, mapping the state of the search to a probability distribution over the approaches. This enables using policy gradient to learn search strategies tailored to a specific distributions of planning problems and a selected performance metric, e.g. the IPC score. We instantiate the framework by constructing a policy space consisting of five search approaches and a two-dimensional representation of the planner's state. Then, we train the system on randomly generated problems from five IPC domains using three different performance metrics. Our experimental results show that the learner is able to discover domain-specific search strategies, improving the planner's performance relative to the baselines of plain best-first search and a uniform policy.


title: Mechanical Properties of Protomene: A Molecular Dynamics Investigation

abstract: Recently, a new class of carbon allotrope called protomene was proposed. This new structure is composed of sp2 and sp3 carbon-bonds. Topologically, protomene can be considered as an sp3 carbon structure (~80% of this bond type) doped by sp2 carbons. First-principles simulations have shown that protomene presents an electronic bandgap of ~3.4 eV. However, up to now, its mechanical properties have not been investigated. In this work, we have investigated protomene mechanical behavior under tensile strain through fully atomistic reactive molecular dynamics simulations using the ReaxFF force field, as available in the LAMMPS code. At room temperature, our results show that the protomene is very stable and the obtained ultimate strength and ultimate stress indicates an anisotropic behavior. The highest ultimate strength was obtained for the x-direction, with a value of ~110 GPa. As for the ultimate strain, the highest one was for the z-direction (~25% of strain) before protomene mechanical fracture.


title: Expansion of coset graphs of PSL_2(F_p)

abstract: Let $G$ be a finite group and let $H_1,H_2<G$ be two subgroups. In this paper, we are concerned with the bipartite graph whose vertices are $G/H_1\cup G/H_2$ and a coset $g_1H_1$ is connected with another coset $g_2H_2$ if and only if $g_1H_1\cap g_2 H_2\neq\varnothing$. The main result of the paper establishes the existence of such graphs with large girth and large spectral gap. Lubotzky, Manning and Wilton use such graphs to construct certain infinite groups of interest in geometric group theory.


title: Multi-linear Monogamy Relations for Three Qubit States

abstract: The monogamy of entanglement means that entanglement cannot be freely shared. In 2014, Oliveira et al. [ Oliveira et al., Phys. Rev. A. 89, 034303 (2014)] proposed a monogamy relation in the linear version and considered it in terms of entanglement of formation. Here we generalize the above version and consider a multi-linear monogamy relation for a multi-qubit system in terms of entanglement of formation and concurrence. Based on the above results, we present an entanglement criterion for genuine entangled states, also we consider the absolutely maximally entangled states and present what an absolutely maximally entangled state is for a three-qubit system. At last, we apply our results to a three-qubit pure state in terms of quantum discord.


title: Loschmidt echo of local dynamical processes in integrable and non integrable spin chains

abstract: The Loschmidt echo is investigated to track the effect of the local QDP. It is also quite sensitive to whether the background dynamics is integrable or not. For the integrable case, viz. the Heisenberg model, the Loschmidt echo depends on the parameters operators corresponding to the QDP as well as the time of QDP. The probability of reviving the system to its initial state is higher for incoherent QDPs occurring at large time intervals. Whereas each time coherent QDP occurs certain probability of reviving the state is always lost. For For the non-integrable case, viz. a kicked Harper model, it exhibits a decaying behaviour when contrasted with integrable dynamics. The decay rate is slower when the corresponding classical Hamiltonian is non chaotic. The Loschmidt echo also distinguishes the integrable and the nonintegrable dynamics when a QDP occurs.


title: Geometry of the Kahan discretizations of planar quadratic Hamiltonian systems

abstract: Kahan discretization is applicable to any quadratic vector field and produces a birational map which approximates the shift along the phase flow. For a planar quadratic Hamiltonian vector field, this map is known to be integrable and to preserve a pencil of cubic curves. Generically, the nine base points of this pencil include three points at infinity (corresponding to the asymptotic directions of cubic curves) and six finite points lying on a conic. We show that the Kahan discretization map can be represented in six different ways as a composition of two Manin involutions, corresponding to an infinite base point and to a finite base point. As a consequence, the finite base points can be ordered so that the resulting hexagon has three pairs of parallel sides which pass through the three base points at infinity. Moreover, this geometric condition on the base points turns out to be characteristic: if it is satisfied, then the cubic curves of the corresponding pencil are invariant under the Kahan discretization of a planar quadratic Hamiltonian vector field.


title: Teleoperated Robotic Arm Movement Using EMG Signal With Wearable MYO Armband

abstract: The main purpose of this research is to move the robotic arm (5DoF) in real-time, based on the surface Electromyography (sEMG) signals, as obtained from the wireless Myo gesture armband to distinguish seven hand movements. The sEMG signals are biomedical signals that estimate and record the electrical signals produced in muscles through their contraction and relaxation, representing neuromuscular activities. Therefore, controlling the robotic arm via the muscles of the human arm using sEMG signals is considered to be one of the most significant methods. The wireless Myo gesture armband is used to record sEMG signals from the forearm. In order to analyze these signals, the pattern recognition system is employed, which consists of three main parts: segmentation, feature extraction, and classification. Overlap technique is chosen for segmenting part of the signal. Six time domain features (MAV, WL, RMS, AR, ZC, and SSC) are extracted from each segment. The classifiers (SVM, LDA, and KNN) are employed to enable comparison between them in order to obtain optimum accuracy of the system. The results show that the SVM achieves higher system accuracy at 96.57 %, compared to LDA reaching 96.01 %, and 92.67 % accuracy achieved by KNN.


title: Machine Learning Accelerated Likelihood-Free Event Reconstruction in Dark Matter Direct Detection

abstract: Reconstructing the position of an interaction for any dual-phase time projection chamber (TPC) with the best precision is key to directly detecting Dark Matter. Using the likelihood-free framework, a new algorithm to reconstruct the 2-D (x; y) position and the size of the charge signal (e) of an interaction is presented. The algorithm uses the charge signal (S2) light distribution obtained by simulating events using a waveform generator. To deal with the computational effort required by the likelihood-free approach, we employ the Bayesian Optimization for Likelihood-Free Inference (BOLFI) algorithm. Together with BOLFI, prior distributions for the parameters of interest (x; y; e) and highly informative discrepancy measures to perform the analyses are introduced. We evaluate the quality of the proposed algorithm by a comparison against the currently existing alternative methods using a large-scale simulation study. BOLFI provides a natural probabilistic uncertainty measure for the reconstruction and it improved the accuracy of the reconstruction over the next best algorithm by up to 15% when focusing on events over a large radii (R > 30 cm, the outer 37% of the detector). In addition, BOLFI provides the smallest uncertainties among all the tested methods.

categories: astro-ph.IM hep-ex hep-ph physics.data-an stat.ML


title: Random Bernstein-Markov factors

abstract: For a polynomial $P_n$ of degree $n$, Bernstein's inequality states that $\|P_n'\| \le n \|P_n\|$ for all $L^p$ norms on the unit circle, $0<p\le\infty,$ with equality for $P_n(z)= c z^n.$ We study this inequality for random polynomials, and show that the expected (average) and almost sure value of $\Vert P_n' \Vert/\Vert P_n\Vert$ is often different from the classical deterministic upper bound $n$. In particular, for circles of radii less than one, the ratio $\Vert P_n' \Vert/\Vert P_n\Vert$ is almost surely bounded as $n$ tends to infinity, and its expected value is uniformly bounded for all degrees under mild assumptions on the random coefficients. For norms on the unit circle, Borwein and Lockhart mentioned that the asymptotic value of $\Vert P_n' \Vert/\Vert P_n\Vert$ in probability is $n/\sqrt{3},$ and we strengthen this to almost sure limit for $p=2.$ If the radius $R$ of the circle is larger than one, then the asymptotic value of $\Vert P_n' \Vert/\Vert P_n\Vert$ in probability is $n/R$, matching the sharp upper bound for the deterministic case. We also obtain bounds for the case $p=\infty$ on the unit circle.


title: The elixir phase of chain molecule

abstract: A phase of matter is a familiar notion for inanimate physical matter. The nature of a phase of matter transcends the microscopic material properties. For example, materials in the liquid phase have certain common properties independent of the chemistry of the constituents: liquids take the shape of the container; they flow; and they can be poured -- alcohol, oil and water as well as a Lennard-Jones computer model exhibit similar behavior when poised in the liquid phase. Here we identify a hitherto unstudied `phase' of matter, the elixir phase, in a simple model of a polymeric chain whose backbone has the correct local cylindrical symmetry induced by the tangent to the chain. The elixir phase appears on breaking the cylindrical symmetry by adding side spheres along the negative normal direction, as in proteins. This phase, nestled between other phases, has multiple ground states made up of building blocks of helices and almost planar sheets akin to protein native folds. We discuss the similarities of this `phase' of a finite size system to the liquid crystal and spin glass phases. Our findings are relevant for understanding proteins; the creation of novel bio-inspired nano-machines; and also may have implications for life elsewhere in the cosmos.


title: RFID based Automated Car Theft Detection and Arresting system

abstract: In this paper, we are presenting a new approach to car theft detection and arresting system using RFID technology. The main purpose of this paper is to establish the concept and the architecture of the whole system. The system which we are able to develop still now is basically a demonstration model or a prototype of the system. The key thing in this system is the passive RFID tags which will be hidden inside the car and act as a unique identification number for the car. The information of all such tags will be maintained by a centralized server. In case of any car being theft, a report should be logged into the server by an authorized user and a network of RFID scanners installed in different check post, traffic signal or toll plazas in the city will search for the reported tag. Once it is found necessary security system would be activated for arresting the car.


title: Inflationary Cosmology: From Theory to Observations

abstract: The main aim of this paper is to provide a qualitative introduction to the cosmic inflation and its relationship with current cosmological observations. The inflationary model solves many of the fundamental problems that challenge the Standard Big Bang cosmology i.e. Flatness, Horizon and Monopole problem, and additionally provides an explanation for the initial conditions observed throughout the Large-Scale Structure of the Universe, such as galaxies. In this review we describe the general solutions carry out by a single scalar field. Then with the use of current surveys, we show the constraints imposed on the inflationary parameters $(n_{\rm s},r)$ which allow us to make the connection between theoretical and observational cosmology. In this way, with the latest results, it is possible to choose or at least to constrain the right inflationary model, parameterised by a single scalar field potential $V(\phi)$.


title: Network Structural Dependency in the Human Connectome Across the Life-Span

abstract: Principles of network topology have been widely studied in the human connectome. Of particular interest is the modularity of the human brain, where the connectome is divided into subnetworks and subsequently changes with development, aging or disease are investigated. We present a weighted network measure, the Network Dependency Index (NDI), to identify an individual region's importance to the global functioning of the network. Importantly, we utilize NDI to differentiate four subnetworks (Tiers) in the human connectome following Gaussian Mixture Model fitting. We analyze the topological aspects of each subnetwork with respect to age and compare it to rich-club based subnetworks (rich-club, feeder and seeder). Our results first demonstrate the efficacy of NDI to identify more consistent, central nodes of the connectome across age-groups, when compared to the rich-club framework. Stratifying the connectome by NDI led to consistent subnetworks across the life-span revealing distinct patterns associated with age where, e.g., the key relay nuclei and cortical regions are contained in a subnetwork with highest NDI. The divisions of the human connectome derived from our data-driven NDI framework have the potential to reveal topological alterations described by network measures through the life-span.


title: Enhancing Stock Movement Prediction with Adversarial Training

abstract: This paper contributes a new machine learning solution for stock movement prediction, which aims to predict whether the price of a stock will be up or down in the near future. The key novelty is that we propose to employ adversarial training to improve the generalization of a neural network prediction model. The rationality of adversarial training here is that the input features to stock prediction are typically based on stock price, which is essentially a stochastic variable and continuously changed with time by nature. As such, normal training with static price-based features (e.g. the close price) can easily overfit the data, being insufficient to obtain reliable models. To address this problem, we propose to add perturbations to simulate the stochasticity of price variable, and train the model to work well under small yet intentional perturbations. Extensive experiments on two real-world stock data show that our method outperforms the state-of-the-art solution with 3.11% relative improvements on average w.r.t. accuracy, validating the usefulness of adversarial training for stock prediction task.


title: Non-data-aided SNR Estimation for QPSK Modulation in AWGN Channel

abstract: Signal-to-noise ratio (SNR) estimation is an important parameter that is required in any receiver or communication systems. It can be computed either by a pilot signal data-aided approach in which the transmitted signal would be known to the receiver, or without any knowledge of the transmitted signal, which is a non-data-aided (NDA) estimation approach. In this paper, a NDA SNR estimation algorithm for QPSK signal is proposed. The proposed algorithm modifies the existing Signal- to-Variation Ratio (SVR) SNR estimation algorithm in the aim to reduce its bias and mean square error in case of negative SNR values at low number of samples of it. We first present the existing SVR algorithm and then show the mathematical derivation of the new NDA algorithm. In addition, we compare our algorithm to two baselines estimation methods, namely the M2M4 and SVR algorithms, using different test cases. Those test cases include low SNR values, extremely high SNR values and low number of samples. Results showed that our algorithm had a better performance compared to second and fourth moment estimation (M2M4) and original SVR algorithms in terms of normalized mean square error (NMSE) and bias estimation while keeping almost the same complexity as the original algorithms


title: The Power of Low Frequencies: Faraday Tomography in the sub-GHz regime

abstract: Faraday tomography, the study of the distribution of extended polarized emission by strength of Faraday rotation, is a powerful tool for studying magnetic fields in the interstellar medium of our Galaxy and nearby galaxies. The strong frequency dependence of Faraday rotation results in very different observational strengths and limitations for different frequency regimes. I discuss the role these effects take in Faraday tomography below 1 GHz, emphasizing the 100-200 MHz band observed by the Low Frequency Array and the Murchison Widefield Array. With that theoretical context, I review recent Faraday tomography results in this frequency regime, and discuss expectations for future observations.


title: Hypergeometric function and Modular Curvature I. --Hypergeometric functions in Heat Coefficients

abstract: We give a new proof of the rearrangement lemma that works for all dimensions and all heat coefficients in the study of modular geometry on noncommutative tori. The building blocks of the spectral functions are landed in a hypergeometric family knowns as Lauricella functions of type $D$. We investigate the differential and recursive relations among the family and obtain a full reduction to Gauss hypergeometric functions. As for applications, we give phase one demonstration on how the hypergeometric features lead to new simplifications in computations involved in the modular geometry.


title: Health Monitoring of Critical Power System Equipments using Identifying Codes

abstract: High voltage power transformers are one of the most critical equipments in the electric power grid. A sudden failure of a power transformer can significantly disrupt bulk power delivery. Before a transformer reaches its critical failure state, there are indicators which, if monitored periodically, can alert an operator that the transformer is heading towards a failure. One of the indicators is the signal to noise ratio (SNR) of the voltage and current signals in substations located in the vicinity of the transformer. During normal operations, the width of the SNR band is small. However, when the transformer heads towards a failure, the widths of the bands increase, reaching their maximum just before the failure actually occurs. This change in width of the SNR can be observed by sensors, such as phasor measurement units (PMUs) located nearby. Identifying Code is a mathematical tool that enables one to uniquely identify one or more {\em objects of interest}, by generating a unique signature corresponding to those objects, which can then be detected by a sensor. In this paper, we first describe how Identifying Code can be utilized for detecting failure of power transformers. Then, we apply this technique to determine the fewest number of sensors needed to uniquely identify failing transformers in different test systems.


title: Brand > Logo: Visual Analysis of Fashion Brands

abstract: While lots of people may think branding begins and ends with a logo, fashion brands communicate their uniqueness through a wide range of visual cues such as color, patterns and shapes. In this work, we analyze learned visual representations by deep networks that are trained to recognize fashion brands. In particular, the activation strength and extent of neurons are studied to provide interesting insights about visual brand expressions. The proposed method identifies where a brand stands in the spectrum of branding strategy, i.e., from trademark-emblazoned goods with bold logos to implicit no logo marketing. By quantifying attention maps, we are able to interpret the visual characteristics of a brand present in a single image and model the general design direction of a brand as a whole. We further investigate versatility of neurons and discover "specialists" that are highly brand-specific and "generalists" that detect diverse visual features. A human experiment based on three main visual scenarios of fashion brands is conducted to verify the alignment of our quantitative measures with the human perception of brands. This paper demonstrate how deep networks go beyond logos in order to recognize clothing brands in an image.


title: Preprocessor Selection for Machine Learning Pipelines

abstract: Much of the work in metalearning has focused on classifier selection, combined more recently with hyperparameter optimization, with little concern for data preprocessing. Yet, it is generally well accepted that machine learning applications require not only model building, but also data preprocessing. In other words, practical solutions consist of pipelines of machine learning operators rather than single algorithms. Interestingly, our experiments suggest that, on average, data preprocessing hinders accuracy, while the best performing pipelines do actually make use of preprocessors. Here, we conduct an extensive empirical study over a wide range of learning algorithms and preprocessors, and use metalearning to determine when one should make use of preprocessors in ML pipeline design.


title: Dynamics of Order Parameters of Non-stoquastic Hamiltonians in the Adaptive Quantum Monte Carlo Method

abstract: We derive macroscopically deterministic flow equations with regard to the order parameters of the ferromagnetic $p$-spin model with infinite-range interactions. The $p$-spin model has a first-order phase transition for $p>2$. In the case of $p\geq5$ ,the $p$-spin model with anti-ferromagnetic XX interaction has a second-order phase transition in a certain region. In this case, however, the model becomes a non-stoqustic Hamiltonian, resulting in a negative sign problem. To simulate the $p$-spin model with anti-ferromagnetic XX interaction, we utilize the adaptive quantum Monte Carlo method. By using this method, we can regard the effect of the anti-ferromagnetic XX interaction as fluctuations of the transverse magnetic field. A previous study derived deterministic flow equations of the order parameters in the quantum Monte Carlo method. In this study, we derive macroscopically deterministic flow equations for the magnetization and transverse magnetization from the master equation in the adaptive quantum Monte Carlo method. Under the Suzuki-Trotter decomposition, we consider the Glauber-type stochastic process. We solve these differential equations by using the Runge-Kutta method and verify that these results are consistent with the saddle-point solution of mean-field theory. Finally, we analyze the stability of the equilibrium solutions obtained by the differential equations.


title: Data models for service failure prediction in supply-chain networks

abstract: We aim to predict and explain service failures in supply-chain networks, more precisely among last-mile pickup and delivery services to customers. We analyze a dataset of 500,000 services using (1) supervised classification with Random Forests, and (2) Association Rules. Our classifier reaches an average sensitivity of 0.7 and an average specificity of 0.7 for the 5 studied types of failure. Association Rules reassert the importance of confirmation calls to prevent failures due to customers not at home, show the importance of the time window size, slack time, and geographical location of the customer for the other failure types, and highlight the effect of the retailer company on several failure types. To reduce the occurrence of service failures, our data models could be coupled to optimizers, or used to define counter-measures to be taken by human dispatchers.


title: Analyzing Neuroimaging Data Through Recurrent Deep Learning Models

abstract: The application of deep learning (DL) models to neuroimaging data poses several challenges, due to the high dimensionality, low sample size and complex temporo-spatial dependency structure of these datasets. Even further, DL models act as as black-box models, impeding insight into the association of cognitive state and brain activity. To approach these challenges, we introduce the DeepLight framework, which utilizes long short-term memory (LSTM) based DL models to analyze whole-brain functional Magnetic Resonance Imaging (fMRI) data. To decode a cognitive state (e.g., seeing the image of a house), DeepLight separates the fMRI volume into a sequence of axial brain slices, which is then sequentially processed by an LSTM. To maintain interpretability, DeepLight adapts the layer-wise relevance propagation (LRP) technique. Thereby, decomposing its decoding decision into the contributions of the single input voxels to this decision. Importantly, the decomposition is performed on the level of single fMRI volumes, enabling DeepLight to study the associations between cognitive state and brain activity on several levels of data granularity, from the level of the group down to the level of single time points. To demonstrate the versatility of DeepLight, we apply it to a large fMRI dataset of the Human Connectome Project. We show that DeepLight outperforms conventional approaches of uni- and multivariate fMRI analysis in decoding the cognitive states and in identifying the physiologically appropriate brain regions associated with these states. We further demonstrate DeepLight's ability to study the fine-grained temporo-spatial variability of brain activity over sequences of single fMRI samples.


title: Lattice quantum gravity with scalar fields

abstract: We consider the four-dimensional Euclidean dynamical triangulations lattice model of quantum gravity based on triangulations of $S^{4}$. We couple it minimally to a scalar field in the quenched approximation. Our results suggest a multiplicative renormalization for the mass of the scalar field which is consistent with the shift symmetry of the discretized lattice action. We discuss the possibility of measuring the mass anomalous dimension and the gravitational binding energy between two scalar test particles, where a negative bound state energy would imply that this model has an attractive gravitational force.


title: Object-oriented lexical encoding of multiword expressions: Short and sweet

abstract: Multiword expressions (MWEs) exhibit both regular and idiosyncratic properties. Their idiosyncrasy requires lexical encoding in parallel with their component words. Their (at times intricate) regularity, on the other hand, calls for means of flexible factorization to avoid redundant descriptions of shared properties. However, so far, non-redundant general-purpose lexical encoding of MWEs has not received a satisfactory solution. We offer a proof of concept that this challenge might be effectively addressed within eXtensible MetaGrammar (XMG), an object-oriented metagrammar framework. We first make an existing metagrammatical resource, the FrenchTAG grammar, MWE-aware. We then evaluate the factorization gain during incremental implementation with XMG on a dataset extracted from an MWE-annotated reference corpus.


title: Further Towards Unambiguous Edge Bundling: Investigating Power-Confluent Drawings for Network Visualization

abstract: Bach et al. [1] recently presented an algorithm for constructing confluent drawings, by leveraging power graph decomposition to generate an auxiliary routing graph. We identify two issues with their method which we call the node split and short-circuit problems, and solve both by modifying the routing graph to retain the hierarchical structure of power groups. We also classify the exact type of confluent drawings that the algorithm can produce as 'power-confluent', and prove that it is a subclass of the previously studied 'strict confluent' drawing. A description and source code of our implementation is also provided, which additionally includes an improved method for power graph construction.


title: Stepwise Acquisition of Dialogue Act Through Human-Robot Interaction

abstract: A dialogue act (DA) represents the meaning of an utterance at the illocutionary force level (Austin 1962) such as a question, a request, and a greeting. Since DAs take charge of the most fundamental part of communication, we believe that the elucidation of DA learning mechanism is important for cognitive science and artificial intelligence. The purpose of this study is to verify that scaffolding takes place when a human teaches a robot, and to let a robot learn to estimate DAs and to make a response based on them step by step utilizing scaffolding provided by a human. To realize that, it is necessary for the robot to detect changes in utterance and rewards given by the partner and continue learning accordingly. Experimental results demonstrated that participants who continued interaction for a sufficiently long time often gave scaffolding for the robot. Although the number of experiments is still insufficient to obtain a definite conclusion, we observed that 1) the robot quickly learned to respond to DAs in most cases if the participants only spoke utterances that match the situation, 2) in the case of participants who builds scaffolding differently from what we assumed, learning did not proceed quickly, and 3) the robot could learn to estimate DAs almost exactly if the participants kept interaction for a sufficiently long time even if the scaffolding was unexpected.


title: Upper bounds for Courant-sharp Neumann and Robin eigenvalues

abstract: We consider the eigenvalues of the Laplacian on an open, bounded, connected set in $\mathbb{R}^n$ with $C^2$ boundary, with a Neumann boundary condition or a Robin boundary condition. We obtain upper bounds for those eigenvalues that have a corresponding eigenfunction which achieves equality in Courant's Nodal Domain theorem. In the case where the set is also assumed to be convex, we obtain explicit upper bounds in terms of some of the geometric quantities of the set.


title: GhostVLAD for set-based face recognition

abstract: The objective of this paper is to learn a compact representation of image sets for template-based face recognition. We make the following contributions: first, we propose a network architecture which aggregates and embeds the face descriptors produced by deep convolutional neural networks into a compact fixed-length representation. This compact representation requires minimal memory storage and enables efficient similarity computation. Second, we propose a novel GhostVLAD layer that includes {\em ghost clusters}, that do not contribute to the aggregation. We show that a quality weighting on the input faces emerges automatically such that informative images contribute more than those with low quality, and that the ghost clusters enhance the network's ability to deal with poor quality images. Third, we explore how input feature dimension, number of clusters and different training techniques affect the recognition performance. Given this analysis, we train a network that far exceeds the state-of-the-art on the IJB-B face recognition dataset. This is currently one of the most challenging public benchmarks, and we surpass the state-of-the-art on both the identification and verification protocols.


title: Agent-Based Modeling and Simulation of Connected and Automated Vehicles Using Game Engine: A Cooperative On-Ramp Merging Study

abstract: Agent-based modeling and simulation (ABMS) has been a popular approach to modeling autonomous and interacting agents in a multi-agent system. Specifically, ABMS can be applied to connected and automated vehicles (CAVs), since CAVs can be driven autonomously with the help of on-board sensors, and cooperate with each other through vehicle-to-everything (V2X) communications. In this work, we apply ABMS to CAVs using the game engine Unity3D, taking advantage of its visualization capability and other capabilities. Agent-based models of CAVs are built in the Unity3D environment, where vehicles are enabled with connectivity and autonomy by C#-based scripting API. We also build a simulation network in Unity3D based on the city of Mountain View, California. A case study of cooperative on-ramp merging has been carried out with the proposed distributed consensus-based protocol, and then compared with the human-in-the-loop simulation where the on-ramp vehicle is driven by four different human drivers on a driving simulator. The benefits of introducing the proposed protocol are evaluated in terms of travel time, energy consumption, and pollutant emissions. It is shown from the results that the proposed cooperative on-ramp merging protocol can reduce average travel time by 7%, reduce energy consumption and pollutant emissions by 8% and 58%, respectively, and guarantee the driving safety when compared to the human-in-the-loop scenario.


title: Dwarf Galaxies: Their Low Metallicity Interstellar Medium

abstract: This review describes where we are today in light of the dust and gas properties and their relation to star formation, in low metallicity galaxies of the local universe following recent surveys from sensitive infrared space telescopes, mainly Spitzer and Herschel space observatories as well as ground-based observations of the molecular gas reservoir. Models to interpret the ISM properties are gaining sophistication in order to account for the wide range of valuable observational diagnostics that we have today to trace the different gas phases, the broad range of photometry we have, from mid-infrared to submillimetre dust emission and the various galactic size scales that we can sample today. This review summarizes the rich multi-phase observations we can exploit today, and the multi-phase modeling approach to interpret the observations.


title: Sharply $k$-arc-transitive-digraphs: finite and infinite examples

abstract: A general method for constructing sharply $k$-arc-transitive digraphs, i.e. digraphs that are $k$-arc-transitive but not $(k+1)$-arc-transitive, is presented. Using our method it is possible to construct both finite and infinite examples. The infinite examples can have one, two or infinitely many ends. Among the one-ended examples there are also digraphs that have polynomial growth.


title: First demonstration of 3D optical readout of a TPC using a single photon sensitive Timepix3 based camera

abstract: The ARIADNE project is developing innovative optical readout technologies for two-phase liquid Argon time projection chambers (LArTPCs). Optical readout presents an exciting alternative to the current paradigm of charge readout. Optical readout is simple, scalable and cost effective. This paper presents first demonstration of 3D optical readout of TPC, using CF4 gas as a proof of principle. Both cosmic rays and an Americium-241 alpha source have been imaged in 100 mbar CF4. A single-photon sensitive camera was developed by combining a Timepix3 (TPX3) based camera with an image intensifier. When a pixel of TPX3 is hit, a packet containing all information about the hit is produced. This packet contains the x,y coordinates of the pixel, time of arrival (ToA) and time over threshold (ToT) information. The z position of the hit in the TPC is determined by combining drift velocity with ToA information. 3D event reconstruction is performed by combining the pixels x,y location with this calculated z position. Calorimetry is performed using time over threshold, a measure of the intensity of the hit.


title: Social Status and Communication Behavior in an Evolving Social Network

abstract: The degree to which individuals can exert influence on propagation of information and opinion dynamics in online communities is highly dependent on their social status. Therefore, there is a high demand for identifying influential users in a community by predicting their social position in that community. Moreover, understanding how people with various social status behave, can shed light on the dynamics of interaction in social networks. In this paper, I study an evolving online social network originated from an online community for university students and I tackle the problem of forecasting users' social status, represented as their PageRank, based on frequency of recurring temporal sequences of observed behavior, i.e. behavioral motifs. I show that individuals with different values of PageRank exhibit different behavior even in early weeks since the online community's inception and it is possible to forecast future PageRank values given frequency of behavioral motifs with high accuracy.


title: NSML: Meet the MLaaS platform with a real-world case study

abstract: The boom of deep learning induced many industries and academies to introduce machine learning based approaches into their concern, competitively. However, existing machine learning frameworks are limited to sufficiently fulfill the collaboration and management for both data and models. We proposed NSML, a machine learning as a service (MLaaS) platform, to meet these demands. NSML helps machine learning work be easily launched on a NSML cluster and provides a collaborative environment which can afford development at enterprise scale. Finally, NSML users can deploy their own commercial services with NSML cluster. In addition, NSML furnishes convenient visualization tools which assist the users in analyzing their work. To verify the usefulness and accessibility of NSML, we performed some experiments with common examples. Furthermore, we examined the collaborative advantages of NSML through three competitions with real-world use cases.


title: ISA Mapper: A Compute and Hardware Agnostic Deep Learning Compiler

abstract: Domain specific accelerators present new challenges and opportunities for code generation onto novel instruction sets, communication fabrics, and memory architectures. In this paper we introduce an intermediate representation (IR) which enables both deep learning computational kernels and hardware capabilities to be described in the same IR. We then formulate and apply instruction mapping to determine the possible ways a computation can be performed on a hardware system. Next, our scheduler chooses a specific mapping and determines the data movement and computation order. In order to manage the large search space of mappings and schedules, we developed a flexible framework that allows heuristics, cost models, and potentially machine learning to facilitate this search problem. With this system, we demonstrate the automated extraction of matrix multiplication kernels out of recent deep learning kernels such as depthwise-separable convolution. In addition, we demonstrate two to five times better performance on DeepBench sized GEMMs and GRU RNN execution when compared to state-of-the-art (SOTA) implementations on new hardware and up to 85% of the performance for SOTA implementations on existing hardware.


title: Doping-induced persistent spin helix with large spin splitting in monolayer SnSe

abstract: Finding a new class of materials supporting a long spin lifetime is essential in development of energy-saving spintronics, which is achievable by using a persistent spin helix (PSH) materials. However, for spintronic devices, the PSH states with large spin splitting is required for operation at room temperature. By employing first-principles calculations, we show that the PSH states with large spin splitting are achieved in the SnSe monolayer (ML) functionalized by a substitutional halogen impurity. We find the PSH states in the Fermi level where k-space Fermi surface is characterized by the shifted two loops, dominated by out-of-plane spin orientations. We clarify the PSH states in term of an effective k.p Hamiltonian obtained from symmetry consideration. Finally, large spin-orbit strength in the PSH states with a substantially small wavelength are found, rendering that this system is promising for the development of an efficient and high-density scalable spintronic devices operating at room temperatures.


title: Deficient and multiple points of maps into a manifold

abstract: For a map $f:X \to M$ into a manifold $M$, we study the sets of deficient and multiple points of $f$. In case of the set of deficient points, we estimate its dimension. For multiple points, we study its density in $X$, and we also provide examples where its complement is dense.


title: Global Well-posedness of the Two Dimensional Beris-Edwards System with General Laudau-de Gennes Free Energy

abstract: In this paper, we consider the Beris-Edwards system for incompressible nematic liquid crystal flows. The system under investigation consists of the Navier-Stokes equations for the fluid velocity $\mathbf{u}$ coupled with an evolution equation for the order parameter $Q$-tensor. One important feature of the system is that its elastic free energy takes a general form and in particular, it contains a cubic term that possibly makes it unbounded from below. In the two dimensional periodic setting, we prove that if the initial $L^\infty$-norm of the $Q$-tensor is properly small, then the system admits a unique global weak solution. The proof is based on the construction of a specific approximating system that preserves the $L^\infty$-norm of the $Q$-tensor along the time evolution.


title: Majorana-Weyl crossings in topological multi-terminal junctions

abstract: We analyze the Andreev spectrum in a four-terminal Josephson junction between one-dimensional topological superconductors in class D. We find that a topologically protected crossing in the space of three superconducting phase differences can occur between the two lowest Andreev bound states. This crossing can be detected through the transconductance quantization, in units of $2e^2/h$, between two voltage-biased terminals. Our prediction provides yet another example of topology in multi-terminal Josephson junctions. We discuss possible realizations of such junctions with semiconducting crossed nanowires and with quantum-spin Hall insulators.


title: Parabolic Jets from the Spinning Black Hole in M87

abstract: The M87 jet is extensively examined by utilizing general relativistic magnetohydrodynamic (GRMHD) simulations as well as the steady axisymmetric force-free electrodynamic (FFE) solution. Quasi-steady funnel jets are obtained in GRMHD simulations up to the scale of $\sim 100$ gravitational radius ($r_{\rm g}$) for various black hole (BH) spins. As is known, the funnel edge is approximately determined by the following equipartitions; i) the magnetic and rest-mass energy densities and ii) the gas and magnetic pressures. Our numerical results give an additional factor that they follow the outermost parabolic streamline of the FFE solution, which is anchored to the event horizon on the equatorial plane. We also identify the matter dominated, non-relativistic corona/wind play a dynamical role in shaping the funnel jet into the parabolic geometry. We confirm a quantitative overlap between the outermost parabolic streamline of the FFE jet and the edge of jet sheath in VLBI observations at $\sim 10^{1}$-$10^{5} \, r_{\rm g}$, suggesting that the M87 jet is likely powered by the spinning BH. Our GRMHD simulations also indicate a lateral stratification of the bulk acceleration (i.e., the spine-sheath structure) as well as an emergence of knotty superluminal features. The spin characterizes the location of the jet stagnation surface inside the funnel. We suggest that the limb-brightened feature could be associated with the nature of the BH-driven jet, if the Doppler beaming is a dominant factor. Our findings can be examined with (sub-)mm VLBI observations, giving a clue for the origin of the M87 jet.

categories: astro-ph.HE astro-ph.GA physics.flu-dyn physics.plasm-ph


title: Transient probability currents provide upper and lower bounds on non-equilibrium steady-state currents in the Smoluchowski picture

abstract: Probability currents are fundamental in characterizing the kinetics of non-equilibrium processes. Notably, the steady-state current $J_{ss}$ for a source-sink system can provide the exact mean-first-passage time (MFPT) for the transition from source to sink. Because transient non-equilibrium behavior is quantified in some modern path sampling approaches, such as the "weighted ensemble" strategy, there is strong motivation to determine bounds on $J_{ss}$ -- and hence on the MFPT -- as the system evolves in time. Here we show that $J_{ss}$ is bounded from above and below by the maximum and minimum, respectively, of the current as a function of the spatial coordinate at any time $t$ for one-dimensional systems undergoing over-damped Langevin (i.e., Smoluchowski) dynamics and for higher-dimensional Smoluchowski systems satisfying certain assumptions when projected onto a single dimension. These bounds become tighter with time, making them of potential practical utility in a scheme for estimating $J_{ss}$ and the long-timescale kinetics of complex systems. Conceptually, the bounds result from the fact that extrema of the transient currents relax toward the steady-state current.

categories: cond-mat.stat-mech physics.bio-ph physics.chem-ph physics.comp-ph


title: Using Deep Learning for price prediction by exploiting stationary limit order book features

abstract: The recent surge in Deep Learning (DL) research of the past decade has successfully provided solutions to many difficult problems. The field of quantitative analysis has been slowly adapting the new methods to its problems, but due to problems such as the non-stationary nature of financial data, significant challenges must be overcome before DL is fully utilized. In this work a new method to construct stationary features, that allows DL models to be applied effectively, is proposed. These features are thoroughly tested on the task of predicting mid price movements of the Limit Order Book. Several DL models are evaluated, such as recurrent Long Short Term Memory (LSTM) networks and Convolutional Neural Networks (CNN). Finally a novel model that combines the ability of CNNs to extract useful features and the ability of LSTMs' to analyze time series, is proposed and evaluated. The combined model is able to outperform the individual LSTM and CNN models in the prediction horizons that are tested.


title: Pseudovarieties of ordered completely regular semigroups

abstract: This paper is a contribution to the theory of finite semigroups and their classification in pseudovarieties, which is motivated by its connections with computer science. The question addressed is what role can play the consideration of an order compatible with the semigroup operation. In the case of unions of groups, so-called completely regular semigroups, the problem of which new pseudovarieties appear in the ordered context is solved. As applications, it is shown that the lattice of pseudovarieties of ordered completely regular semigroups is modular and that taking the intersection with the pseudovariety of bands defines a complete endomorphism of the lattice of all pseudovarieties of ordered semigroups.


title: Reconciling $\lambda$-Returns with Experience Replay

abstract: Modern deep reinforcement learning methods have departed from the incremental learning required for eligibility traces, rendering the implementation of the $\lambda$-return difficult in this context. In particular, off-policy methods that utilize experience replay remain problematic because their random sampling of minibatches is not conducive to the efficient calculation of $\lambda$-returns. Yet replay-based methods are often the most sample efficient, and incorporating $\lambda$-returns into them is a viable way to achieve new state-of-the-art performance. Towards this, we propose the first method to enable practical use of $\lambda$-returns in arbitrary replay-based methods without relying on other forms of decorrelation such as asynchronous gradient updates. By promoting short sequences of past transitions into a small cache within the replay memory, adjacent $\lambda$-returns can be efficiently precomputed by sharing Q-values. Computation is not wasted on experiences that are never sampled, and stored $\lambda$-returns behave as stable temporal-difference (TD) targets that replace the target network. Additionally, our method grants the unique ability to observe TD errors prior to sampling; for the first time, transitions can be prioritized by their true significance rather than by a proxy to it. Furthermore, we propose the novel use of the TD error to dynamically select $\lambda$-values that facilitate faster learning. We show that these innovations can enhance the performance of DQN when playing Atari 2600 games, even under partial observability. While our work specifically focuses on $\lambda$-returns, these ideas are applicable to any multi-step return estimator.


title: Schr\"odinger Functional of a Quantum Scalar Field in Static Space-Times from Precanonical Quantization

abstract: The functional Schr\"odinger representation of a scalar field on an $n$-dimensional static space-time background is argued to be a singular limiting case of the hypercomplex quantum theory of the same system obtained by the precanonical quantization based on the space-time symmetric De Donder-Weyl Hamiltonian theory. The functional Schr\"odinger representation emerges from the precanonical quantization when the ultraviolet parameter $\varkappa$ introduced by precanonical quantization is replaced by $\underline{\gamma}{}_0\delta^\mathrm{inv}(\mathbf{0})$, where $\underline{\gamma}{}_0$ is the time-like tangent space Dirac matrix and $\delta^\mathrm{inv}(\mathbf{0})$ is an invariant spatial $(n-1)$-dimensional Dirac's delta function whose regularized value at $\mathbf{x}=\mathbf{0}$ is identified with the cutoff of the volume of the momentum space. In this limiting case, the Schr\"odinger wave functional is expressed as the trace of the product integral of Clifford-algebra-valued precanonical wave functions restricted to a certain field configuration and the canonical functional derivative Schr\"odinger equation is derived from the manifestly covariant Dirac-like precanonical Schr\"odinger equation which is independent of a choice of a codimension-one foliation.


title: Simulative Suzuki-Gerghaty type contraction with $\mathcal{C}$-class functions and applications

abstract: The aim of this paper is to introduce the notion of a Suzuki-Gerghaty type contractive mapping via simulation function along with $\mathcal{C}$-class functions and prove the existence of fixed point result. An example is given to show the validity of our results given herein. Moreover, we prove the existence of solution of nonlinear Hammerstein integral equation.


title: Chirality meets Topology

abstract: By considering the topology of chiral crystals a new type of massless fermions, connected with giant arc-like surface states, are predicted. Such Kramers-Weyl fermions should manifest in a wide variety of chiral materials


title: Quantum Interference Control of Localized Carrier Distributions in the Brillouin Zone

abstract: Using transition-metal dichalcogenides as an example, we show that the quantum interference arising in two- and three-photon absorption processes can lead to controllable, highly localized carrier distributions in the Brillouin zone. We contrast this with the previously studied one- and two-photon absorption, and find qualitatively different features, including changes in the relevance of interband and intraband processes according to the excitation energy. Thus, the distribution of excitations arising under certain circumstances in two- and three-photon absorption can facilitate the study of far-from-equilibrium states that are initially well localized in crystal momentum space.


title: Experimental studies of the transmission of light through low coverage regular or random arrays of silica micropillars supported by a glass substrate

abstract: The transmission of light through low coverage regular and random arrays of glass supported silica micropillars of diameters 10 to 40 \micro\meter and height 10 \micro\meter is studied experimentally. Angle-resolved measurements of the transmitted intensity are performed at visible wavelengths by either a goniospectrophotometer or a multimodal imaging (Mueller) polarimetric microscope. It is demonstrated that for the regular arrays, the angle-resolved measurements are capable of resolving many of the densely packed diffraction orders that are expected for periodic structures of lattice constants 20 to 80 \micro\meter, but they also display features that are due to the scattering and guiding of light in individual micropillars or in the supporting glass slides. These latter features are also found in angle-resolved measurements on random arrays of micropillars of the same surface coverage. Finally we perform a comparison of direct measurements of haze in transmission for our patterned glass samples with what can be calculated from the angle-resolved transmitted intensity measurements. Good agreement between the two types of results are found which testifies to the accuracy of the angle-resolved measurements that we report.


title: Production of strange particles in hadronic interactions

abstract: The NA61/SHINE collaboration has recently published high precision data on production of $\pi^\pm$ and $K^\pm$ mesons, protons, antiprotons and $\Lambda$ hyperons in ${\rm pp}$ interactions at 20, 31, 40, 80 and 158 GeV/c, and in ${\rm pC}$ interactions at 31 GeV/c. The collaboration also presented experimental data on production of particles - $\pi^\pm$, $K^\pm$, $p^\pm$, $\rho^0$, $\omega$ and $K^{*0}$ in $\pi^-{\rm C}$ collisions at 158 and 350 GeV/c. The collaboration has compared these data with various Monte Carlo model calculations: UrQMD, EPOS, GiBUU, and others. All of the models have various problems. The latest version of the FTF (Fritiof) model of Geant4 solves most of these problems. In the FTF model, we have improved the fragmentation of quark-gluon strings with small masses and introduced dependencies of probabilities of strange mesons and baryon-antibaryon pair's creation on string masses. Due to these changes, we describe the data of the NA61/SHINE collaboration on particle production in ${\rm pp, pC}$, and $\pi^-{\rm C}$ interactions. The improved Geant4 FTF model also well reproduces experimental data on inclusive cross sections of $\Lambda, \bar{\Lambda}$ and $K^{0}$ production in antiproton-proton interactions at various energies. The modified FTF model allows one to simulate realistic processes with two particle productions - $\bar{p}p \rightarrow \Lambda \bar{\Lambda}$, $\bar{p}p \rightarrow K^{+} K^{-}$, $\bar{p}p \rightarrow \Lambda \bar{\Sigma}$, and $\bar{p}p \rightarrow\Sigma \bar{\Sigma}$, which will be studied in the future by the PANDA experiment at FAIR (GSI, Germany).


title: Algorithmic Traversals of Infinite Graphs

abstract: A traversal of a connected graph is a linear ordering of its vertices all of whose initial segments induce connected subgraphs. Traversals, and their refinements such as breadth-first and depth-first traversals, are computed by various graph searching algorithms. We extend the theory of generic search and breadth-first search from finite graphs to wellordered infinite graphs, recovering the notion of "search trees" in this context. We also prove tight upper bounds on the extent to which graph search and breadth-first search can modify the order type of the original graph, as well as characterize the traversals computed by these algorithms as lexicographically minimal.


title: Jump sets in local fields

abstract: We show how to use the combinatorial notion of jump sets to parametrize the possible structures of the group of principal units of local fields, viewed as filtered modules. We establish a natural bijection between the set of jump sets and the orbit space of a $p$-adic group of filtered automorphisms acting on a free filtered module. This, together with a Markov process on Eisenstein polynomials, culminates into a mass-formula for unit filtrations. As a bonus the proof leads in many cases to explicit invariants of Eisenstein polynomials, yielding a link between the filtered structure of the unit group and ramification theory. Finally, with the basic theory of filtered modules developed here, we recover, with a more conceptual proof, a classification, due to Miki, of the possible sets of upper jumps of a wild character: these are all jump sets, with a set of exceptions explicitly prescribed by the jump set of the local field and the size of its residue field.


title: Evidence of unconventional superconductivity in the Ni-doped NbB2 system

abstract: In this paper, a comprehensive study of the effects of Ni-doping on structural, electrical, thermal and magnetic properties of the NbB2 is presented. Low amounts (\leq 10 %) of Ni substitution on Nb sites cause structural distortions and induce drastic changes in the physical properties, such as the emergence of a bulk superconducting state with anomalous behaviors in the critical fields (lower and upper) and in the specific heat. Ni-doping at the 9 at.% level, for instance, is able to increase the critical temperature (TC) in stoichiometric NbB2 (< 1.3 K) to approximately 6.0 K. Bulk superconductivity is confirmed by magnetization, electronic transport, and specific heat measurements. Both Hc1 and Hc2 critical fields exhibit a linear dependence with reduced temperature (T/TC), and the specific heat deviates remarkably from the conventional exponential temperature dependence of the single-band BCS theory. These findings suggest multiband superconductivity in the composition range from 0.01 \leq x \leq 0.10.