title: Quantitative theoretical analysis of lifetimes and decay rates relevant in laser cooling BaH

abstract: Tiny radiative losses below the 0.1% level can prove ruinous to the effective laser cooling of a molecule. In this paper the laser cooling of a hydride is studied with rovibronic detail using ab initio quantum chemistry in order to document the decays to all possible electronic states (not just the vibrational branching within a single electronic transition) and to identify the most populated final quantum states. The effect of spin-orbit and associated couplings on the properties of the lowest excited states of BaH are analysed in detail. The lifetimes of the A$^2{\Pi}_{1/2}$, H$^2{\Delta}_{3/2}$ and E$^2{\Pi}_{1/2}$ states are calculated (136 ns, 5.8 {\mu}s and 46 ns respectively) for the first time, while the theoretical value for B$^2{\Sigma}^+_{1/2}$ is in good agreement with experiments. Using a simple rate model the numbers of absorption-emission cycles possible for both one- and two-colour cooling on the competing electronic transitions are determined, and it is clearly demonstrated that the A$^2{\Pi}$ - X$^2{\Sigma}^+$ transition is superior to B$^2{\Sigma}^+$ - X$^2{\Sigma}^+$, where multiple tiny decay channels degrade its efficiency. Further possible improvements to the cooling method are proposed.


title: Three-dimensional Imaging for Large LArTPCs

abstract: High-performance event reconstruction is critical for current and future massive liquid argon time projection chambers (LArTPCs) to realize their full scientific potential. LArTPCs with readout using wire planes provide a limited number of 2D projections. In general, without a pixel-type readout it is challenging to achieve unambiguous 3D event reconstruction. As a remedy, we present a novel 3D imaging method, Wire-Cell, which incorporates the charge and sparsity information in addition to the time and geometry through simple and robust mathematics. The resulting 3D image of ionization density provides an excellent starting point for further reconstruction and enables the true power of 3D tracking calorimetry in LArTPCs.


title: Reaction rate theory for supramolecular kinetics: application to protein aggregation

abstract: Probing the reaction mechanisms of supramolecular processes in soft- and biological matter, such as protein aggregation, is inherently challenging. These processes emerge from the simultaneous action of multiple molecular mechanisms, each of which is associated with the rearrangement of a large number of weak bonds, resulting in a complex free energy landscape with many kinetic barriers. Reaction rate measurements of supramolecular processes at different temperatures can offer unprecedented insights into the underlying molecular mechanisms and their thermodynamic properties. However, to be able to interpret such measurements in terms of the underlying microscopic mechanisms, a key challenge is to establish which properties of the complex free energy landscapes are probed by the reaction rate. Here, we present a reaction rate theory for supramolecular kinetics based on Kramers rate theory for diffusive reactions over multiple kinetic barriers, and apply the results to protein aggregation. Using this framework and Monte Carlo simulations, we show that reaction rates for protein aggregation are of the Arrhenius-Eyring type and that the associated activation energies probe only one relevant barrier along the respective free energy landscapes. We apply this advancement to interpret, both in experiments and in coarse-grained computer simulations, reaction rate measurements of amyloid aggregation kinetics in terms of the underlying molecular mechanisms and associated thermodynamic signatures. Our results establish a general platform for probing the mechanisms and energetics of supramolecular phenomena in soft- and biological matter using the framework of chemical kinetics.


title: Local optimality of Zaks-Perles-Wills simplices

abstract: In 1982, Zaks, Perles and Wills discovered a d-dimensional lattice simplex S_{d,k} with k interior lattice points, whose volume is linear in k and doubly exponential in the dimension d. It is conjectured that, for all d \ge 3 and k \ge 1, the simplex S_{d,k} is a volume maximizer in the family P^d(k) of all d-dimensional lattice polytopes with k interior lattice points. To obtain a partial confirmation of this conjecture, one can try to verify it for a subfamily of P^d(k) that naturally contains S_{d,k} as one of the members. Currently, one does not even know whether S_{d,k} is optimal within the family S^d(k) of all d-dimensional lattice simplices with k interior lattice points. In view of this, it makes sense to look at even narrower families, for example, some subfamilies of S^d(k). The simplex S_{d,k} of Zaks, Perles and Wills has a facet with only one lattice point in the relative interior. We show that S_{d,k} is a volume maximizer in the family of simplices S \in S^d(k) that have a facet with one lattice point in its relative interior. We also show that, in the above family, the volume maximizer is unique up to unimodular transformations.


title: Regularized Bidimensional Estimation of the Hazard Rate

abstract: In epidemiological or demographic studies, with variable age at onset, a typical quantity of interest is the incidence of a disease (for example the cancer incidence). In these studies, the individuals are usually highly heterogeneous in terms of dates of birth (the cohort) and with respect to the calendar time (the period) and appropriate estimation methods are needed. In this article a new estimation method is presented which extends classical age-period-cohort analysis by allowing interactions between age, period and cohort effects. This paper introduces a bidimensional regularized estimate of the hazard rate where a penalty is introduced on the likelihood of the model. This penalty can be designed either to smooth the hazard rate or to enforce consecutive values of the hazard to be equal, leading to a parsimonious representation of the hazard rate. In the latter case, we make use of an iterative penalized likelihood scheme to approximate the L0 norm, which makes the computation tractable. The method is evaluated on simulated data and applied on breast cancer survival data from the SEER program.


title: M101: Spectral Observations of HII Regions and Their Physical Properties

abstract: By using the Hectospec 6.5 m Multiple Mirror Telescope (MMT) and the 2.16 m telescope of National Astronomical Observatories, Chinese Academy of Sciences (NAOC), we obtained 188 high signal-to-noise ratio (S/N) spectra of HII regions in the nearby galaxy M101, which are the largest spectroscopic sample of HII regions for this galaxy so far. These spectra cover a wide range of regions on M101, which enables us to analyze two dimensional distributions of its physical properties. The physical parameters are derived from emission lines or stellar continuum, including stellar population age, electron temperature, oxygen abundance and etc. The oxygen abundances are derived using two empirical methods based on O3N2 and R$_{23}$ indicators, as well as the direct Te method when OIII$\lambda4363$ is available. By applying the harmonic decomposition analysis to the velocity field, we obtained line-of-sight rotation velocity of 71 km s$^{-1}$ and a position angle of 36 degree. The stellar age profile shows an old stellar population in galaxy center and a relative young stellar population in outer regions, suggesting an old bulge and a young disk. Oxygen abundance profile exhibits a clear break at $\sim$18 kpc, with a gradient of $-$0.0364 dex kpc$^{-1}$ in the inner region and $-$0.00686 dex kpc$^{-1}$ in the outer region. Our results agree with the "inside-out" disk growth scenario of M101.


title: Atomization of correlated molecular-hydrogen chain: A fully microscopic Variational Monte-Carlo solution

abstract: We discuss electronic properties and their evolution for the linear chain of $H_2$ molecules in the presence of a uniform external force $f$ acting along the chain. The system is described by an extended Hubbard model within a fully microscopic approach. Explicitly, the microscopic parameters describing the intra- and inter-site Coulomb interactions are determined together with the hopping integrals by optimizing the system ground state energy and the single-particle wave functions in the correlated state. The many-body wave function is taken in the Jastrow form and the Variational Monte-Carlo (VMC) method is used in combination with an ab initio approach to determine the energy. Both the effective Bohr radii of the renormalized single-particle wave functions and the many-body wave function parameters are determined for each $f$. Hence, the evolution of the system can be analyzed in detail as a function of the equilibrium intermolecular distance, which in turn is determined for each $f$ value. The transition to the atomic state, including the Peierls distortion stability, can thus be studied in a systematic manner, particularly near the threshold of the dissociation of the molecular into atomic chain. The computational reliability of VMC approach is also estimated.


title: A System for the Generation of Synthetic Wide Area Aerial Surveillance Imagery

abstract: The development, benchmarking and validation of aerial Persistent Surveillance (PS) algorithms requires access to specialist Wide Area Aerial Surveillance (WAAS) datasets. Such datasets are difficult to obtain and are often extremely large both in spatial resolution and temporal duration. This paper outlines an approach to the simulation of complex urban environments and demonstrates the viability of using this approach for the generation of simulated sensor data, corresponding to the use of wide area imaging systems for surveillance and reconnaissance applications. This provides a cost-effective method to generate datasets for vehicle tracking algorithms and anomaly detection methods. The system fuses the Simulation of Urban Mobility (SUMO) traffic simulator with a MATLAB controller and an image generator to create scenes containing uninterrupted door-to-door journeys across large areas of the urban environment. This `pattern-of-life' approach provides three-dimensional visual information with natural movement and traffic flows. This can then be used to provide simulated sensor measurements (e.g. visual band and infrared video imagery) and automatic access to ground-truth data for the evaluation of multi-target tracking systems.


title: Efficient white noise sampling and coupling for multilevel Monte Carlo with non-nested meshes

abstract: When solving stochastic partial differential equations (SPDEs) driven by additive spatial white noise, the efficient sampling of white noise realizations can be challenging. Here, we present a new sampling technique that can be used to efficiently compute white noise samples in a finite element method and multilevel Monte Carlo (MLMC) setting. The key idea is to exploit the finite element matrix assembly procedure and factorize each local mass matrix independently, hence avoiding the factorization of a large matrix. Moreover, in a MLMC framework, the white noise samples must be coupled between subsequent levels. We show how our technique can be used to enforce this coupling even in the case of non-nested mesh hierarchies. We demonstrate the efficacy of our method with numerical experiments. We observe optimal convergence rates for the finite element solution of the elliptic SPDEs of interest in 2D and 3D and we show convergence of the sampled field covariances. In a MLMC setting, a good coupling is enforced and the telescoping sum is respected.


title: Expert identification of visual primitives used by CNNs during mammogram classification

abstract: This work interprets the internal representations of deep neural networks trained for classification of diseased tissue in 2D mammograms. We propose an expert-in-the-loop interpretation method to label the behavior of internal units in convolutional neural networks (CNNs). Expert radiologists identify that the visual patterns detected by the units are correlated with meaningful medical phenomena such as mass tissue and calcificated vessels. We demonstrate that several trained CNN models are able to produce explanatory descriptions to support the final classification decisions. We view this as an important first step toward interpreting the internal representations of medical classification CNNs and explaining their predictions.


title: On moments of integral exponential functionals of additive processes

abstract: For real-valued additive process $(X\_t)\_{t\geq 0}$ a recursive equation is derived for the entire positive moments of functionals $$I\_{s,t}= \int \_s^t\exp(-X\_u)du, \quad 0\leq s<t\leq\infty, $$ in case the Laplace exponent of $X\_t$ exists for positive values of the parameter. From the equation emergesan easy-to-apply sufficient condition for the finiteness of the moments. As an application we study first hitprocesses of diffusions.


title: NECTAR: Non-Interactive Smart Contract Protocol using Blockchain Technology

abstract: Blockchain-driven technologies are considered disruptive because of the availability of dis-intermediated, censorship-resistant and tamper-proof digital platforms of distributed trust. Among these technologies, smart contract platforms have the potential to take over functions usually done by intermediaries like banks, escrow or legal services. In this paper, we introduce a novel protocol aiming to execute smart contracts as part of a blockchain transaction validation. We enable extensions in the execution of smart contracts while guaranteeing their privacy, correctness and verifiability. Man-in-the-middle attacks are prevented, since no communication between participants is requested, and contract validations do not imply the re-execution of the code by all the nodes in the network. However, proofs of correct execution are stored on the blockchain and can be verified by multiple parties. Our solution is based on programming tools which optimize the time execution and the required memory while preserving the embedded functionality.


title: SHARVOT: secret SHARe-based VOTing on the blockchain

abstract: Recently, there has been a growing interest in using online technologies to design protocols for secure electronic voting. The main challenges include vote privacy and anonymity, ballot irrevocability and transparency throughout the vote counting process. The introduction of the blockchain as a basis for cryptocurrency protocols, provides for the exploitation of the immutability and transparency properties of these distributed ledgers. In this paper, we discuss possible uses of the blockchain technology to implement a secure and fair voting system. In particular, we introduce a secret share-based voting system on the blockchain, the so-called SHARVOT protocol. Our solution uses Shamir's Secret Sharing to enable on-chain, i.e. within the transactions script, votes submission and winning candidate determination. The protocol is also using a shuffling technique, Circle Shuffle, to de-link voters from their submissions.


title: Correlation Manipulating Circuits for Stochastic Computing

abstract: Stochastic computing (SC) is an emerging computing technique that promises high density, low power, and error tolerant solutions. In SC, values are encoded as unary bitstreams and SC arithmetic circuits operate on one or more bitstreams. In many cases, the input bitstreams must be correlated or uncorrelated for SC arithmetic to produce accurate results. As a result, a key challenge for designing SC accelerators is manipulating the impact of correlation across SC operations. This paper presents and evaluates a set of novel correlation manipulating circuits to manage correlation in SC computation: a synchronizer, desynchronizer, and decorrelator. We then use these circuits to propose improved SC maximum, minimum, and saturating adder designs. Compared to existing correlation manipulation techniques, our circuits are more accurate and up to 3x more energy efficient. In the context of an image processing pipeline, these circuits can reduce the total energy consumption by up to 24%.


title: Exact short-time height distribution for the flat Kardar-Parisi-Zhang interface

abstract: We determine the exact short-time distribution $-\ln \mathcal{P}_{\text{f}}\left(H,t\right)= S_{\text{f}} \left(H\right)/\sqrt{t}$ of the one-point height $H=h(x=0,t)$ of an evolving 1+1 Kardar-Parisi-Zhang (KPZ) interface for flat initial condition. This is achieved by combining (i) the optimal fluctuation method, (ii) a time-reversal symmetry of the KPZ equation in 1+1 dimension, and (iii) the recently determined exact short-time height distribution $-\ln \mathcal{P}_{\text{st}}\left(H,t\right)= S_{\text{st}} \left(H\right)/\sqrt{t}$ for \emph{stationary} initial condition. In studying the large-deviation function $S_{\text{st}} \left(H\right)$ of the latter, one encounters two branches: an analytic and a non-analytic. The analytic branch is non-physical beyond a critical value of $H$ where a second-order dynamical phase transition occurs. Here we show that, remarkably, it is the analytic branch of $S_{\text{st}} \left(H\right)$ which determines the large-deviation function $S_{\text{f}} \left(H\right)$ of the flat interface via a simple mapping $S_{\text{f}}\left(H\right)=2^{-3/2}S_{\text{st}}\left(2H\right)$.


title: Resource Allocation in Wireless Networks with Energy Constraints

abstract: This dissertation focuses on the development of novel scheduling and resource allocation schemes, which take into account and regulate the energy constraints imposed by the levels of harvested energy. To this direction, first, the optimal energy, time, and bandwidth allocation problem for the downlink of energy harvesting base stations (EHBSs) is investigated, with the main focus being on autonomous EHBSs. The presented analysis considers the impact of the energy constraint on users' preferences and the BS's revenue. In order to model the competitive nature of the problem, game theory is used. The next two chapters focus on wireless powered networks (WPNs) and simultaneous wireless information and power transfer (SWIPT) using radio frequency (RF) technology. One of the main contributions of these chapters is the introduction of both uplink and downlink non-orthogonal multiple access (NOMA) for WPNs. Moreover, the individual data rates and fairness are improved, while the formulated problems are optimally and efficiently solved. It is shown that, compared to orthogonal multiple access, NOMA offers a considerable improvement in throughput, fairness, and energy efficiency. Rather than this, proportional fairness is maximized and uplink/downlink of WPNs are jointly optimized, in which cases, except for NOMA, time division multiple access (TDMA) is also investigated. Furthermore, the role of interference is considered, which has been recognized as one of the main reasons of the asymmetric overall degradation of the users' performance, due to different path-loss values, called from now on as cascaded near-far problem. Moreover, SWIPT is investigated and efficiently optimized in the context of multicarrier cooperative communication networks. Finally, simultaneous lightwave information and power transfer (SLIPT) is introduced, while novel and fundamental techniques are proposed.


title: Mathematical models and search algorithms for the capacitated $p$-center problem

abstract: The capacitated p-center problem requires to select p facilities from a set of candidates to service a number of customers, subject to facility capacity constraints, with the aim of minimizing the maximum distance between a customer and its associated facility. The problem is well known in the field of facility location, because of the many applications that it can model. In this paper, we solve it by means of search algorithms that iteratively seek the optimal distance by solving tailored subproblems. We present different mathematical formulations for the subproblems and improve them by means of several valid inequalities, including an effective one based on a 0-1 disjunction and the solution of subset sum problems. We also develop an alternative search strategy that finds a balance between the traditional sequential search and binary search. This strategy limits the number of feasible subproblems to be solved and, at the same time, avoids large overestimates of the solution value, which are detrimental for the search. We evaluate the proposed techniques by means of extensive computational experiments on benchmark instances from the literature and new larger test sets. All instances from the literature with up to 402 vertices and integer distances are solved to proven optimality, including 13 open cases, and feasible solutions are found in 10 minutes for instances with up to 3038 vertices.


title: Dealing with Unknown Unknowns: Identification and Selection of Minimal Sensing for Fractional Dynamics with Unknown Inputs

abstract: This paper focuses on analysis and design of time-varying complex networks having fractional order dynamics. These systems are key in modeling the complex dynamical processes arising in several natural and man made systems. Notably, examples include neurophysiological signals such as electroencephalogram (EEG) that captures the variation in potential fields, and blood oxygenation level dependent (BOLD) signal, which serves as a proxy for neuronal activity. Notwithstanding, the complex networks originated by locally measuring EEG and BOLD are often treated as isolated networks and do not capture the dependency from external stimuli, e.g., originated in subcortical structures such as the thalamus and the brain stem. Therefore, we propose a paradigm-shift towards the analysis of such complex networks under unknown unknowns (i.e., excitations). Consequently, the main contributions of the present paper are threefold: (i) we present an alternating scheme that enables to determine the best estimate of the model parameters and unknown stimuli; (ii) we provide necessary and sufficient conditions to ensure that it is possible to retrieve the state and unknown stimuli; and (iii) upon these conditions we determine a small subset of variables that need to be measured to ensure that both state and input can be recovered, while establishing sub-optimality guarantees with respect to the smallest possible subset. Finally, we present several pedagogical examples of the main results using real data collected from an EEG wearable device.


title: Extreme field-sensitivity of the magnetic tunneling in Fe-doped Li$_3$N

abstract: The magnetic properties of dilute Li$_2$(Li$_{1-x}$Fe$_x$)N with $x \sim 0.001$ are dominated by the spin of single, isolated Fe atoms. Below $T = 10$ K the spin-relaxation times become temperature-independent indicating a crossover from thermal excitations to the quantum tunneling regime. We report on a strong increase of the spin-flip probability in $\textit{transverse}$ magnetic fields that proves the resonant character of this tunneling process. $\textit{Longitudinal}$ fields, on the other hand, lift the ground-state degeneracy and destroy the tunneling condition. An increase of the relaxation time by four orders of magnitude in applied fields of only a few milliTesla reveals exceptionally sharp tunneling resonances. Li$_2$(Li$_{1-x}$Fe$_x$)N represents a comparatively simple and clean model system that opens the possibility to study quantum tunneling of the magnetization at liquid helium temperatures.


title: Search-based optimal motion planning for automated driving

abstract: This paper presents a framework for fast and robust motion planning designed to facilitate automated driving. The framework allows for real-time computation even for horizons of several hundred meters and thus enabling automated driving in urban conditions. This is achieved through several features. Firstly, a convenient geometrical representation of both the search space and driving constraints enables the use of classical path planning approach. Thus, a wide variety of constraints can be tackled simultaneously (other vehicles, traffic lights, etc.). Secondly, an exact cost-to-go map, obtained by solving a relaxed problem, is then used by A*-based algorithm with model predictive flavour in order to compute the optimal motion trajectory. The algorithm takes into account both distance and time horizons. The approach is validated within a simulation study with realistic traffic scenarios. We demonstrate the capability of the algorithm to devise plans both in fast and slow driving conditions, even when full stop is required.


title: Rapidly evolving transients in the Dark Energy Survey

abstract: We present the results of a search for rapidly evolving transients in the Dark Energy Survey Supernova Programme. These events are characterized by fast light curve evolution (rise to peak in $\lesssim 10$ d and exponential decline in $\lesssim30$ d after peak). We discovered 72 events, including 37 transients with a spectroscopic redshift from host galaxy spectral features. The 37 events increase the total number of rapid optical transients by more than factor of two. They are found at a wide range of redshifts ($0.05<z<1.56$) and peak brightnesses ($-15.75>M_\mathrm{g}>-22.25$). The multiband photometry is well fit by a blackbody up to few weeks after peak. The events appear to be hot ($T\approx10000-30000$ K) and large ($R\approx 10^{14}-2\cdot10^{15}$ cm) at peak, and generally expand and cool in time, though some events show evidence for a receding photosphere with roughly constant temperature. Spectra taken around peak are dominated by a blue featureless continuum consistent with hot, optically thick ejecta. We compare our events with a previously suggested physical scenario involving shock breakout in an optically thick wind surrounding a core-collapse supernova (CCSNe), we conclude that current models for such a scenario might need an additional power source to describe the exponential decline. We find these transients tend to favor star-forming host galaxies, which could be consistent with a core-collapse origin. However, more detailed modeling of the light curves is necessary to determine their physical origin.


title: Narcissus: Deriving Correct-By-Construction Decoders and Encoders from Binary Formats

abstract: It is a neat result from functional programming that libraries of parser combinators can support rapid construction of decoders for quite a range of formats. With a little more work, the same combinator program can denote both a decoder and an encoder. Unfortunately, the real world is full of gnarly formats, as with the packet formats that make up the standard Internet protocol stack. Most past parser-combinator approaches cannot handle these formats, and the few exceptions require redundancy -- one part of the natural grammar needs to be hand-translated into hints in multiple parts of a parser program. We show how to recover very natural and nonredundant format specifications, covering all popular network packet formats and generating both decoders and encoders automatically. The catch is that we use the Coq proof assistant to derive both kinds of artifacts using tactics, automatically, in a way that guarantees that they form inverses of each other. We used our approach to reimplement packet processing for a full Internet protocol stack, inserting our replacement into the OCaml-based MirageOS unikernel, resulting in minimal performance degradation.


title: Anisotropic susceptibilities in the honeycomb Kitaev system $\alpha$-RuCl$_{3}$

abstract: The magnetic insulator $\alpha$-RuCl$_{3}$ is a promising candidate to realize Kitaev interactions on a quasi-2D honeycomb lattice. We perform extensive susceptibility measurements on single crystals of $\alpha$-RuCl$_{3}$, including angle-dependence of the in-plane longitudinal and transverse susceptibilities, which reveal a unidirectional anisotropy within the honeycomb plane. By comparing the experimental results to a high-temperature expansion of a Kitaev-Heisenberg-$\Gamma$ spin Hamiltonian with bond anisotropy, we find excellent agreement with the observed phase shift and periodicity of the angle-resolved susceptibilities. Within this model, we show that the pronounced difference between in-plane and out-of-plane susceptibilities as well as the finite transverse susceptibility are rooted in strong symmetric off-diagonal $\Gamma$ spin exchange. The $\Gamma$ couplings and relationships between other terms in the model Hamiltonian are quantified by extracting relevant Curie-Weiss intercepts from the experimental data.


title: Project Blue: Visible Light Imaging Search for Terrestrial-class Exoplanets in the Habitable Zones of Alpha Centauri A and B

abstract: Project Blue is designed to deliver a small coronagraphic space telescope mission to low-Earth orbit capable of detecting an Earth-like planet in the habitable zones of the Sun-like stars Alpha Centauri A and B within the next 4 years within a Mission of Opportunity budget envelope. The concept heavily leverages emerging commercial capabilities -- including the telescope optics, spacecraft bus and launch vehicle -- and emphasizes a public-private partnership approach.


title: Using Convolutional Neural Networks for Determining Reticulocyte Percentage in Cats

abstract: Recent advances in artificial intelligence (AI), specifically in computer vision (CV) and deep learning (DL), have created opportunities for novel systems in many fields. In the last few years, deep learning applications have demonstrated impressive results not only in fields such as autonomous driving and robotics, but also in the field of medicine, where they have, in some cases, even exceeded human-level performance. However, despite the huge potential, adoption of deep learning-based methods is still slow in many areas, especially in veterinary medicine, where we haven't been able to find any research papers using modern convolutional neural networks (CNNs) in medical image processing. We believe that using deep learning-based medical imaging can enable more accurate, faster and less expensive diagnoses in veterinary medicine. In order to do so, however, these methods have to be accessible to everyone in this field, not just to computer scientists. To show the potential of this technology, we present results on a real-world task in veterinary medicine that is usually done manually: feline reticulocyte percentage. Using an open source Keras implementation of the Single-Shot MultiBox Detector (SSD) model architecture and training it on only 800 labeled images, we achieve an accuracy of 98.7% at predicting the correct number of aggregate reticulocytes in microscope images of cat blood smears. The main motivation behind this paper is to show not only that deep learning can approach or even exceed human-level performance on a task like this, but also that anyone in the field can implement it, even without a background in computer science.


title: Robust Recursive Filtering and Smoothing

abstract: Using a perturbation technique, we derive a new approximate filtering and smoothing methodology generalizing along different directions several existing approaches to robust filtering based on the score and the Hessian matrix of the observation density. The main advantages of the methodology can be summarized as follows: (i) it relaxes the critical assumption of a Gaussian prior distribution for the latent states underlying such approaches; (ii) can be applied to a general class of state-space models including univariate and multivariate location, scale and count data models; (iii) has a very simple structure based on forward-backward recursions similar to the Kalman filter and smoother; (iv) allows a straightforward computation of confidence bands around the state estimates reflecting the combination of parameter and filtering uncertainty. We show through an extensive Monte Carlo study that the mean square loss with respect to exact simulation-based methods is small in a wide range of scenarios. We finally illustrate empirically the application of the methodology to the estimation of stochastic volatility and correlations in financial time-series.


title: Generating B\'ezout trees for Pythagorean pairs

abstract: Relatively prime pairs of integers can be represented as nodes in three way branching trees. We construct trees of B\'ezout coefficients which correspond to the relatively prime pairs in the aforementioned trees. As one application, we compare the B\'ezout coefficients in these trees to those returned by the gcd function in Matlab. As another application, we use these trees to decrease the computation time required to create computer generated hyperbolic wallpaper designs.


title: Hypergeometric continuation of divergent perturbation series. I. Critical exponents of the Bose-Hubbard model

abstract: We study the connection between the exponent of the order parameter of the Mott insulator-to-superfluid transition occurring in the two-dimensional Bose-Hubbard model, and the divergence exponents of its one- and two-particle correlation functions. We find that at the multicritical points all divergence exponents are related to each other, allowing us to express the critical exponent in terms of one single divergence exponent. This approach correctly reproduces the critical exponent of the three-dimensional $XY$ universality class. Because divergence exponents can be computed in an efficient manner by hypergeometric analytic continuation, our strategy is applicable to a wide class of systems.


title: A machine learning-based approach for estimating and testing associations with multivariate outcomes

abstract: We propose a method for summarizing the strength of association between a set of variables and a multivariate outcome. Classical summary measures are appropriate when linear relationships exist between covariates and outcomes, while our approach provides an alternative that is useful in situations where complex relationships may be present. We utilize ensemble machine learning to detect nonlinear relationships and covariate interactions and propose a measure of association that captures these relationships. A hypothesis test about the proposed associative measure can be used to test the strong null hypothesis of no association between a set of variables and a multivariate outcome. Simulations demonstrate that this hypothesis test has greater power than existing methods against alternatives where covariates have nonlinear relationships with outcomes. We additionally propose measures of variable importance for groups of variables, which summarize each groups' association with the outcome. We demonstrate our methodology using data from a birth cohort study on childhood health and nutrition in the Philippines.


title: Effect of network clustering on mutually cooperative coinfections

abstract: The spread of an infectious disease can be promoted by previous infections with other pathogens. This cooperative effect can give rise to violent outbreaks, reflecting the presence of an abrupt epidemic transition. As for other diffusive dynamics, the topology of the interaction pattern of the host population plays a crucial role. It was conjectured that a discontinuous transition arises when there are relatively few short loops and many long loops in the contact network. Here we focus on the role of local clustering in determining the nature of the transition. We consider two mutually cooperative pathogens diffusing in the same population: an individual already infected with one disease has an increased probability of getting infected by the other. We look at how a disease obeying the susceptible-infected-removed dynamics spreads on contact networks with tunable clustering. Using numerical simulations we show that for large cooperativity the epidemic transition is always abrupt, with the discontinuity decreasing as clustering is increased. For large clustering strong finite size effects are present and the discontinuous nature of the transition is manifest only in large networks. We also investigate the problem of influential spreaders for cooperative infections, revealing that both cooperativity and clustering strongly enhance the dependence of the spreading influence on the degree of the initial seed.


title: Grigorchuk-Gupta-Sidki groups as a source for Beauville surfaces

abstract: If $G$ is a Grigorchuk-Gupta-Sidki group defined over a $p$-adic tree, where $p$ is an odd prime, we study the existence of Beauville surfaces associated to the quotients of $G$ by its level stabilizers $\st_G(n)$. We prove that if $G$ is periodic then the quotients $G/\st_G(n)$ are Beauville groups for every $n\geq 2$ if $p\geq 5$ and $n\geq 3$ if $p=3$. On the other hand, if $G$ is non-periodic, then none of the quotients $G/\st_G(n)$ are Beauville groups.


title: An Efficient Data Protection Architecture Based on Fragmentation and Encryption

abstract: In this thesis, a completely revisited data protection scheme based on selective encryption is presented. First, this new scheme is agnostic in term of data format, second it has a parallel architecture using GPGPU allowing performance to be at least comparable to full encryption algorithms. Bitmap, as a special uncompressed multimedia format, is addressed as a first use case. Discrete Cosine Transform (DCT) is the first transformation for splitting fragments, getting data protection, and storing data separately on local device and cloud servers. This work has largely improved the previous published ones for bitmap protection by providing new designs and practical experimentations. General purpose graphic processing unit (GPGPU) is exploited as an accelerator to guarantee the efficiency of the calculation compared with traditional full encryption algorithms. Then, an agnostic selective encryption based on lossless Discrete Wavelet Transform (DWT) is presented. This design, with practical experimentations on different hardware configurations, provides strong level of protection and good performance at the same time plus flexible storage dispersion schemes. Therefore, our agnostic data protection and transmission solution combining fragmentation, encryption, and dispersion is made available for a wide range of end-user applications. Also a complete set of security analysis are deployed to test the level of provided protection.


title: Reviewing KLEE's Sonar-Search Strategy in Context of Greybox Fuzzing

abstract: Automatic test-case generation techniques of symbolic execution and fuzzing are the most widely used methods to discover vulnerabilities in, both, academia and industry. However, both these methods suffer from fundamental drawbacks that stop them from achieving high path coverage that may, consequently, lead to discovering vulnerabilities at the numerical scale of static analysis. In this presentation, we examine systems-under-test (SUTs) at the granularity level of functions and postulate that achieving higher function coverage (execution of functions in a program at least once) than, both, symbolic execution and fuzzing may be a necessary condition for discovering more vulnerabilities than both. We will start this presentation with the design of a targeted search strategy for KLEE, sonar-search, that prioritizes paths leading to a target function, rather than maximizing overall path coverage in the program. Then, we will show that examining SUTs at the level of functions (compositional analysis) leads to discovering more vulnerabilities than symbolic execution from a single entry point. Using this finding, we will, then, demonstrate a greybox fuzzing method that can achieve higher function coverage than symbolic execution. Finally, we will present a framework to effectively manage vulnerabilities and assess their severities.


title: Hypergeometric continuation of divergent perturbation series. II. Comparison with Shanks transformation and Pad\'e approximation

abstract: We explore in detail how analytic continuation of divergent perturbation series by generalized hypergeometric functions is achieved in practice. Using the example of strong-coupling perturbation series provided by the two-dimensional Bose-Hubbard model, we compare hypergeometric continuation to Shanks and Pad\'e techniques, and demonstrate that the former yields a powerful, efficient and reliable alternative for computing the phase diagram of the Mott insulator-to-superfluid transition. In contrast to Shanks transformations and Pad\'e approximations, hypergeometric continuation also allows us to determine the exponents which characterize the divergence of correlation functions at the transition points. Therefore, hypergeometric continuation constitutes a promising tool for the study of quantum phase transitions.


title: Melting probe technology for subsurface exploration of extraterrestrial ice - Critical refreezing length and the role of gravity

abstract: The 'Ocean Worlds' of our Solar System are covered with ice, hence the water is not directly accessible. Using melting probe technology is one of the promising technological approaches to reach those scientifically interesting water reservoirs. Melting probes basically consist of a heated melting head on top of an elongated body that contains the scientific payload. The traditional engineering approach to design such melting probes starts from a global energy balance around the melting head and quantifies the power necessary to sustain a specific melting velocity while preventing the probe from refreezing and stall in the channel. Though this approach is sufficient to design simple melting probes for terrestrial applications, it is too simplistic to study the probe's performance for environmental conditions found on some of the Ocean's Worlds, e.g. a lower value of the gravitational acceleration. This will be important, however, when designing exploration technologies for extraterrestrial purposes. We tackle the problem by explicitly modeling the physical processes in the thin melt film between the probe and the underlying ice. Our model allows to study melting regimes on bodies of different gravitational acceleration, and we explicitly compare melting regimes on Europa, Enceladus and Mars. In addition to that, our model allows to quantify the heat losses due to convective transport around the melting probe. We discuss to which extent these heat losses can be utilized to avoid the necessity of a side wall heating system to prevent stall, and introduce the notion of the 'Critical Refreezing Length'. Our results allow to draw important conclusions towards the design of melting probe technology for future missions to icy bodies in our Solar System.


title: IDEL: In-Database Entity Linking with Neural Embeddings

abstract: We present a novel architecture, In-Database Entity Linking (IDEL), in which we integrate the analytics-optimized RDBMS MonetDB with neural text mining abilities. Our system design abstracts core tasks of most neural entity linking systems for MonetDB. To the best of our knowledge, this is the first defacto implemented system integrating entity-linking in a database. We leverage the ability of MonetDB to support in-database-analytics with user defined functions (UDFs) implemented in Python. These functions call machine learning libraries for neural text mining, such as TensorFlow. The system achieves zero cost for data shipping and transformation by utilizing MonetDB's ability to embed Python processes in the database kernel and exchange data in NumPy arrays. IDEL represents text and relational data in a joint vector space with neural embeddings and can compensate errors with ambiguous entity representations. For detecting matching entities, we propose a novel similarity function based on joint neural embeddings which are learned via minimizing pairwise contrastive ranking loss. This function utilizes a high dimensional index structures for fast retrieval of matching entities. Our first implementation and experiments using the WebNLG corpus show the effectiveness and the potentials of IDEL.


title: First passage upwards for state dependent-killed spectrally negative L\'evy processes

abstract: For a spectrally negative L\'evy process (snLp) $X$, killed according to a rate that is a function $\omega$ of its position, we analyse the exit probability of the one-sided upwards-passage problem. When $\omega$ is strictly positive, this problem is related to the determination of the Laplace transform of the first passage time upwards for $X$ that has been time-changed by the inverse of the additive functional $\int_0^\cdot \omega(X_u)du$. In particular our findings thus shed extra light on related results concerning first passage times upwards (downwards) of spectrally negative positive self-similar Markov processes (continuous state branching processes).


title: Examples of hypergeometric twistor $\mathcal{D}$-modules

abstract: We show that certain one-dimensional hypergeometric differential systems underlie objects of the category of irregular mixed Hodge modules, which was recently introduced by Sabbah, and compute the irregular Hodge filtration for them. We also provide a comparison theorem between two different types of Fourier-Laplace transformation for algebraic integrable twistor $\mathcal{D}$-modules.


title: Sub-gap optical response across the structural phase transition in van der Waals layered \alpha-RuCl$_3$

abstract: We report magnetic, thermodynamic, thermal expansion, and on detailed optical experiments on the layered compound $\alpha$-RuCl$_3$ focusing on the THz and sub-gap optical response across the structural phase transition from the monoclinic high-temperature to the rhombohedral low-temperature structure, where the stacking sequence of the molecular layers is changed. This type of phase transition is characteristic for a variety of tri-halides crystallizing in a layered honeycomb-type structure and so far is unique, as the low-temperature phase exhibits the higher symmetry. One motivation is to unravel the microscopic nature of spin-orbital excitations via a study of temperature and symmetry-induced changes. We document a number of highly unusual findings: A characteristic two-step hysteresis of the structural phase transition, accompanied by a dramatic change of the reflectivity. An electronic excitation, which appears in a narrow temperature range just across the structural phase transition, and a complex dielectric loss spectrum in the THz regime, which could indicate remnants of Kitaev physics. Despite significant symmetry changes across the monoclinic to rhombohedral phase transition, phonon eigenfrequencies and the majority of spin-orbital excitations are not strongly influenced. Obviously, the symmetry of the single molecular layers determine the eigenfrequencies of most of these excitations. Finally, from this combined terahertz, far- and mid-infrared study we try to shed some light on the so far unsolved low energy (< 1eV) electronic structure of the ruthenium $4d^5$ electrons in $\alpha$-RuCl$_3$.


title: A look at possible microwave dust emission via AKARI infrared all-sky surveys

abstract: The anomalous microwave emission (AME) still lacks a conclusive explanation. This excess of emission, roughly between 10 and 50 GHz, correlates spatially with interstellar dust, prompting a "spinning dust" hypothesis: electric dipole emission by rapidly rotating, small dust grains. The typical peak frequency range of the AME profile implicates grains on the order of ~1 nm, suggesting polycyclic aromatic hydrocarbon molecules (PAHs). We compare AKARI/Infrared Camera (IRC), with its thorough PAH-band coverage, to AME intensity estimates from the Planck Collaboration, in the {\lambda} Orionis region. We look also at infrared dust emission from other mid IR and far-IR bands. The results and discussion contained here apply to an angular scale of approximately 1{\deg}. In general, our results support an AME-from-dust hypothesis. In {\lambda} Orionis, we find that certainly dust mass correlates with AME, and that PAH-related emission in the AKARI/IRC 9 {\mu}m band may correlate slightly more strongly.


title: Planar anti-Ramsey numbers of matchings

abstract: Given a positive integer $n$ and a planar graph $H$, let $\mathcal{T}_n(H)$ be the family of all plane triangulations $T$ on $n$ vertices such that $T$ contains a subgraph isomorphic to $H$. The planar anti-Ramsey number of $H$, denoted $ar_{_\mathcal{P}}(n, H)$, is the maximum number of colors in an edge-coloring of a plane triangulation $T\in \mathcal{T}_n(H)$ such that $T$ contains no rainbow copy of $H$. In this paper we study planar anti-Ramsey numbers of matchings. For all $t\ge1$, let $M_t$ denote a matching of size $t$. We prove that for all $t\ge6$ and $n\ge 3t-6$, $2n+3t-15\le ar_{_{\mathcal{P}}}(n, {M}_t)\le 2n+4t-14$, which significantly improves the existing lower and upper bounds for $ar_{_\mathcal{P}}(n, M_t)$. It seems that for each $t\ge6$, the lower bound we obtained is the exact value of $ar_{_{\mathcal{P}}}(n, {M}_t)$ for sufficiently large $n$. This is indeed the case for $M_6$. We prove that $ar_{_\mathcal{P}}(n, M_6)=2n+3$ for all $n\ge30$.


title: Complex symmetric differential operators on Fock space

abstract: The space of entire functions which are integrable with respect to the Gaussian weight, known also as the Fock space, is one of the preferred functional Hilbert spaces for modelling and experimenting harmonic analysis, quantum mechanics or spectral analysis phenomena. This space of entire functions carries a three parameter family of canonical isometric involutions. We characterize the linear differential operators acting on Fock space which are complex symmetric with respect to these conjugations. In parallel, as a basis of comparison, we discuss the structure of self-adjoint linear differential operators. The computation of the point spectrum of some of these operators is carried out in detail.


title: MHDSTS: a new explicit numerical scheme for simulations of partially ionised solar plasma

abstract: The interaction of plasma with magnetic field in the partially ionised solar atmosphere is frequently modelled via a single-fluid approximation, which is valid for the case of a strongly coupled collisional media, such as solar photosphere and low chromosphere. Under the single-fluid formalism the main non-ideal effects are described by a series of extra terms in the generalised induction equation and in the energy conservation equation. These effects are: Ohmic diffusion, ambipolar diffusion, the Hall effect, and the Biermann battery effect. From the point of view of the numerical solution of the single-fluid equations, when ambipolar diffusion or Hall effects dominate can introduce severe restrictions on the integration time step and can compromise the stability of the numerical scheme. In this paper we introduce two numerical schemes to overcome those limitations. The first of them is known as Super Time-Stepping (STS) and it is designed to overcome the limitations imposed when the ambipolar diffusion term is dominant. The second scheme is called the Hall Diffusion Scheme (HDS) and it is used when the Hall term becomes dominant. These two numerical techniques can be used together by applying Strang operator splitting. This paper describes the implementation of the STS and HDS schemes in the single-fluid code Mancha3D. The validation for each of these schemes is provided by comparing the analytical solution with the numerical one for a suite of numerical tests.


title: Theoretical and empirical analysis of trading activity

abstract: Understanding the structure of financial markets deals with suitably determining the functional relation between financial variables. In this respect, important variables are the trading activity, defined here as the number of trades $N$, the traded volume $V$, the asset price $P$, the squared volatility $\sigma^2$, the bid-ask spread $S$ and the cost of trading $C$. Different reasonings result in simple proportionality relations ("scaling laws") between these variables. A basic proportionality is established between the trading activity and the squared volatility, i.e., $N \sim \sigma^2$. More sophisticated relations are the so called 3/2-law $N^{3/2} \sim \sigma P V /C$ and the intriguing scaling $N \sim (\sigma P/S)^2$. We prove that these "scaling laws" are the only possible relations for considered sets of variables by means of a well-known argument from physics: dimensional analysis. Moreover, we provide empirical evidence based on data from the NASDAQ stock exchange showing that the sophisticated relations hold with a certain degree of universality. Finally, we discuss the time scaling of the volatility $\sigma$, which turns out to be more subtle than one might naively expect.


title: A single population of red globular clusters around the massive compact galaxy NGC 1277

abstract: Massive galaxies are thought to form in two phases: an initial, early collapse of gas and giant burst of central star formation, followed by the later accretion of material that builds up their stellar and dark matter haloes. The globular cluster systems of such galaxies are believed to form in a similar manner. The initial central burst forms metal-rich (red) clusters, while more metal-poor (blue) clusters are brought in by the later accretion of less massive satellites. This formation process is thought to lead the creation of the multimodal optical colour distributions seen in the globular cluster systems of massive galaxies. Here we report HST/ACS observations of the massive relic galaxy NGC 1277 and its globular clusters, a nearby unevolved example of a high redshift "red nugget". The g-z cluster colour distribution shows that the globular cluster system of the galaxy is unimodal and uniquely red. This is in strong contrast to normal galaxies of similar and larger stellar mass, whose cluster systems always exhibit (and are generally dominated by) blue clusters. We argue that the globular cluster system of NGC 1277 indicates that the galaxy has undergone little (if any) mass accretion after its initial collapse and use analytic merger trees to show that the total stellar mass accretion is likely less than ~ 10 %. These results confirm that NGC 1277 is a genuine relic galaxy and show that the blue, metal-poor globular clusters constitute an accreted population in present day massive galaxies.


title: A Score-Driven Conditional Correlation Model for Noisy and Asynchronous Data: an Application to High-Frequency Covariance Dynamics

abstract: The analysis of the intraday dynamics of correlations among high-frequency returns is challenging due to the presence of asynchronous trading and market microstructure noise. Both effects may lead to significant data reduction and may severely underestimate correlations if traditional methods for low-frequency data are employed. We propose to model intraday log-prices through a multivariate local-level model with score-driven covariance matrices and to treat asynchronicity as a missing value problem. The main advantages of this approach are: (i) all available data are used when filtering correlations, (ii) market microstructure noise is taken into account, (iii) estimation is performed through standard maximum likelihood methods. Our empirical analysis, performed on 1-second NYSE data, shows that opening hours are dominated by idiosyncratic risk and that a market factor progressively emerges in the second part of the day. The method can be used as a nowcasting tool for high-frequency data, allowing to study the real-time response of covariances to macro-news announcements and to build intraday portfolios with very short optimization horizons.


title: Inverse source problem in a forced network

abstract: We address the nonlinear inverse source problem of identifying a time-dependent source occurring in one node of a network governed by a wave equation. We prove that time records of the associated state taken at a strategic set of two nodes yield uniqueness of the two unknown elements: the source position and the emitted signal. We develop a non-iterative identification method that localizes the source node by solving a set of well posed linear systems. Once the source node is localized, we identify the emitted signal using a deconvolution problem or a Fourier expansion. Numerical experiments on a $5$ node graph confirm the effectiveness of the approach.


title: Sub-voxel perfusion modeling in terms of coupled 3d-1d problem

abstract: We study perfusion by a multiscale model coupling diffusion in the tissue and diffusion along the one-dimensional segments representing the vasculature. We propose a block-diagonal preconditioner for the model equations and demonstrate its robustness by numerical experiments. We compare our model to a macroscale model by [P. Tofts, Modelling in DCE MRI, 2012].


title: Explosion in weighted Hyperbolic Random Graphs and Geometric Inhomogeneous Random Graphs

abstract: In this paper we study weighted distances in scale-free spatial network models: hyperbolic random graphs (HRG), geometric inhomogeneous random graphs (GIRG) and scale-free percolation (SFP). In HRGs, $n=\Theta(\mathrm{e}^{R/2})$ vertices are sampled independently from the hyperbolic disk with radius $R$ and two vertices are connected either when they are within hyperbolic distance $R$, or independently with a probability depending on the hyperbolic distance. In GIRGs and SFP, each vertex is given an independent weight and location from an underlying measured metric space and $\mathbb{Z}^d$, respectively, and two vertices are connected independently with a probability that is a function of their distance and weights. We assign i.i.d. weights to the edges of the random graphs and study the weighted distance between two uniformly chosen vertices. In SFP, we study the weighted distance from the origin of vertex-sequences with norm tending to infinity. In particular, we study the case when the parameters are so that the degree distribution in the graph follows a power law with exponent $\tau\in(2,3)$ (infinite variance), and the edge-weight distribution is such that it produces an explosive age-dependent branching process with power-law offspring distribution. We show that in all three models, typical distances within the giant/infinite component converge in distribution, solving an open question in [Explosion and distances in scale-free percolation (2017)]. The main tools of our proof are to couple the models to infinite versions, to follow the shortest paths to infinity and to connect these paths using weight-dependent percolation on the graphs: delete edges attached to vertices with higher weight with higher probability. We realise this using the edge-weights: only short edges connected to high weight vertices will stay, yielding arbitrarily short upper bounds for the connections.


title: Magneto-optic dynamics in a ferromagnetic nematic liquid crystal

abstract: We investigate dynamic magneto-optic effects in a ferromagnetic nematic liquid crystal experimentally and theoretically. Experimentally we measure the magnetization and the phase difference of the transmitted light when an external magnetic field is applied. As a model we study the coupled dynamics of the magnetization, M, and the director field, n, associated with the liquid crystalline orientational order. We demonstrate that the experimentally studied macroscopic dynamic behavior reveals the importance of a dynamic cross-coupling between M and n. The experimental data are used to extract the value of the dissipative cross-coupling coefficient. We also make concrete predictions about how reversible cross-coupling terms between the magnetization and the director could be detected experimentally by measurements of the transmitted light intensity as well as by analyzing the azimuthal angle of the magnetization and the director out of the plane spanned by the anchoring axis and the external magnetic field. We derive the eigenmodes of the coupled system and study their relaxation rates. We show that in the usual experimental set-up used for measuring the relaxation rates of the splay-bend or twist-bend eigenmodes of a nematic liquid crystal one expects for a ferromagnetic nematic liquid crystal a mixture of at least two eigenmodes.


title: Optimal Transport for Multi-source Domain Adaptation under Target Shift

abstract: In this paper, we propose to tackle the problem of reducing discrepancies between multiple domains referred to as multi-source domain adaptation and consider it under the target shift assumption: in all domains we aim to solve a classification problem with the same output classes, but with labels' proportions differing across them. This problem, generally ignored in the vast majority papers on domain adaptation papers, is nevertheless critical in real-world applications, and we theoretically show its impact on the adaptation success. To address this issue, we design a method based on optimal transport, a theory that has been successfully used to tackle adaptation problems in machine learning. Our method performs multi-source adaptation and target shift correction simultaneously by learning the class probabilities of the unlabeled target sample and the coupling allowing to align two (or more) probability distributions. Experiments on both synthetic and real-world data related to satellite image segmentation task show the superiority of the proposed method over the state-of-the-art.


title: Microscopic approach to understand Brownian dynamics in viscoelastic fluid

abstract: We present an entirely microscopic formulation of viscoleasticity of a fluid starting from the microscopic Stokes-Oldroyd B Model assuming instantaneous hydrodynamic friction, and show that linearization leads to a form for the frequency dependent viscosity that can be directly applied to the Langevin equation. Interestingly, the calculated expression of viscosity can be directly mapped to the Jeffreys' model which is essentially macroscopic in nature with the bulk viscoelasticity of the fluid being characterized by a complex elastic modulus G(\omega). Further, we demonstrate that the concerned Green's function is same as that in an incompressible, low Reynold's number Newtonian fluid with the simple incorporation of frequency dependence in the viscosity term. We proceed to evaluate the trajectory of a free Brownian particle in a viscoelastic environment using our formalism and calculate parameters such as the power spectral density, the autocorrelation function and the mean-square displacement, which we then extend to the particle confined in a harmonic potential in the fluid.


title: Topology and the optical Dirac equation

abstract: Through understanding Maxwell's equations as an effective Dirac equation (the `optical Dirac equation'), we re--examine the relationship between electromagnetic interface states and topology. We illustrate a simple case where electromagnetic material parameters play the roles of `mass' and `energy' in an equivalent Dirac equation. The modes trapped between a gyrotropic medium and a mirror are then the counterpart of those at a `domain wall', where the mass of the Dirac particle changes sign. Considering the general case of arbitrary electromagnetic media, we provide an analytical proof relating the integral of the Berry curvature (the Chern number) to the number of interface states. We show that this reduces to the usual result for periodic media, and also that the Chern number can be computed without knowledge of how the material parameters depend on frequency.


title: Higher-order nonclassical properties of a shifted symmetric cat state and a one-dimensional continuous superposition of coherent states

abstract: Role of quantum interference in the origin of higher-order nonclassical characteristics of radiation field has been probed vis-a-vis a discrete and a continuous superposition of coherent states. Specifically, the possibilities of observing higher-order nonclassical properties (e.g., higher-order antibunching (HOA), higher-order sub-Poissonian photon statistics (HOSPS), higher-order squeezing (HOS) of Hong-Mandel type and Hillery type) have been investigated using a shifted symmetric cat state that reduces to Yurke-Stoler, even and odd coherent states at various limits. This shifted symmetric cat state which can be viewed as a discrete superposition of coherent states is found to show HOA and HOSPS. Similarly, higher-order nonclassical properties of a one-dimensional continuous superposition of coherent states is also studied here. The investigation has revealed the existence of HOS and HOSPS in the one-dimensional continuous superposition of coherent states studied here. Effect of non-Gaussianity inducing operations (e.g., photon addition and addition followed by subtraction) on these superposition states have also been investigated. Finally, some comparisons have been made between the higher-order nonclassical properties of discrete and continuous superposition of coherent states.


title: Global secondary bifurcation, symmetry breaking and period-doubling

abstract: In this paper we provide a criterion for global secondary bifurcation via symmetry breaking. As an application, the occurrence of period-doubling bifurcations for the Lugiato-Lefever equation is proved.


title: On cosmology in nonlinear multidimensional gravity with multiple factor spaces

abstract: Within the scope of multidimensional Kaluza--Klein gravity with nonlinear curvature terms and two spherical extra spaces of dimensions $m$ and $n$, we study the properties of an effective action for the scale factors of the extra dimensions. Dimensional reduction leads to an effective 4D multiscalar-tensor theory. Based on qualitative estimates of the Casimir energy contribution on a physically reasonable length scale, we demonstrate the existence of such sets of initial parameters of the theory in the case $m=n$ that provide a minimum of the effective potential that yield a fine-tuned value of the effective 4D cosmological constant. The corresponding size of extra dimensions depends of which conformal frame is interpreted as the observational one: it is about three orders of magnitude larger than the standard Planck length if we adhere to the Einstein frame, but it is $n$-dependent in the Jordan frame, and its invisibility requirement restricts the total dimension to values $D = 4+2n \leq 20$.


title: On the KKLT Goldstino

abstract: We construct general asymptotically Klebanov-Strassler solutions of a five dimensional $SU(2) \times SU(2) \times \mathbb{Z}_2\times \mathbb{Z}_{2R}$ truncation of IIB supergravity on $ T^{1,1} $, that break supersymmetry. This generalizes results in the literature for the $SU(2) \times SU(2) \times \mathbb{Z}_2\times U(1)_R$ case, to a truncation that is general enough to capture the deformation of the conifold in the IR. We observe that there are only two SUSY-breaking modes even in this generalized set up, and by holographically computing Ward identities, we confirm that only one of them corresponds to spontaneous breaking: this is the mode triggered by smeared anti-D3 branes at the tip of the warped throat. Along the way, we address some aspects of the holographic computation of one-point functions of marginal and relevant operators in the cascading gauge theory. Our results strengthen the evidence that {\it if} the KKLT construction is meta-stable, it is indeed a spontaneously SUSY-broken (and therefore bona fide) vacuum of string theory.


title: Diagnostic-Driven Nonstationary Emulators Using Kernel Mixtures

abstract: Weakly stationary Gaussian processes (GPs) are the principal tool in the statistical approaches to the design and analysis of computer experiments (or Uncertainty Quantification). Such processes are fitted to computer model output using a set of training runs to learn the parameters of the process covariance kernel. The stationarity assumption is often adequate, yet can lead to poor predictive performance when the model response exhibits nonstationarity, for example, if its smoothness varies across the input space. In this paper, we introduce a diagnostic-led approach to fitting nonstationary GP emulators by specifying finite mixtures of region-specific covariance kernels. Our method first fits a stationary GP and, if traditional diagnostics exhibit nonstationarity, those diagnostics are used to fit appropriate mixing functions for a covariance kernel mixture designed to capture the nonstationarity, ensuring an emulator that is continuous in parameter space and readily interpretable. We compare our approach to the principal nonstationary GP models in the literature and illustrate its performance on a number of idealised test cases and in an application to modelling the cloud parameterization of the French climate model.


title: Quantization of Fully Convolutional Networks for Accurate Biomedical Image Segmentation

abstract: With pervasive applications of medical imaging in health-care, biomedical image segmentation plays a central role in quantitative analysis, clinical diagno- sis, and medical intervention. Since manual anno- tation su ers limited reproducibility, arduous e orts, and excessive time, automatic segmentation is desired to process increasingly larger scale histopathological data. Recently, deep neural networks (DNNs), par- ticularly fully convolutional networks (FCNs), have been widely applied to biomedical image segmenta- tion, attaining much improved performance. At the same time, quantization of DNNs has become an ac- tive research topic, which aims to represent weights with less memory (precision) to considerably reduce memory and computation requirements of DNNs while maintaining acceptable accuracy. In this paper, we apply quantization techniques to FCNs for accurate biomedical image segmentation. Unlike existing litera- ture on quantization which primarily targets memory and computation complexity reduction, we apply quan- tization as a method to reduce over tting in FCNs for better accuracy. Speci cally, we focus on a state-of- the-art segmentation framework, suggestive annotation [22], which judiciously extracts representative annota- tion samples from the original training dataset, obtain- ing an e ective small-sized balanced training dataset. We develop two new quantization processes for this framework: (1) suggestive annotation with quantiza- tion for highly representative training samples, and (2) network training with quantization for high accuracy. Extensive experiments on the MICCAI Gland dataset show that both quantization processes can improve the segmentation performance, and our proposed method exceeds the current state-of-the-art performance by up to 1%. In addition, our method has a reduction of up to 6.4x on memory usage.


title: Low complexity algorithms in knot theory

abstract: We show that the genus problem for alternating knots with $n$ crossings has linear time complexity and is in Logspace$(n)$. Almost all alternating knots of given genus possess additional combinatorial structure, we call them standard. We show that the genus problem for these knots belongs to $TC^0$ circuit complexity class. We also show, that the equivalence problem for such knots with $n$ crossings has time complexity $n\log (n)$ and is in Logspace$(n)$ and $TC^{0}$ complexity classes.


title: Around logical perfection

abstract: In this article we present and describe a notion of "logical perfection". We extract the notion of "perfection" from the contemporary logical concept of categoricity. Categoricity (in power) has become in the past half century a main driver of ideas in model theory, both mathematically (stability theory may be regarded as a way of approximating categoricity) and philosophically. In the past two decades, categoricity notions have started to overlap with more classical notions of robustness and smoothness. These have been crucial in various parts of mathematics since the nineteenth century. We postulate and present the category of logical perfection. We draw on various notions of perfection from mathematics of the 19th and 20th centuries and then trace the relation to the concept of categoricity in power as a logical notion of what a "mathematically perfect" structure is.


title: Experimental evidence of robust acoustic valley Hall edge states in a non-resonant topological elastic waveguide

abstract: This paper presents experimental evidence of the existence of acoustic valley Hall (AVHE) edge states in topological elastic waveguides. The fundamental lattice is assembled based on a non-resonant unit where space inversion symmetry (SIS) is broken by simply perturbing the underlying lattice geometry. This aspect is in net contrast with existing elastic AVHE designs that exploit locally-resonant units and require the addition of masses in order to break SIS. The experimental results presented in this study validate findings so far presented only at theoretical and numerical level. In particular, it is found that edge modes can effectively propagate along domain walls between topologically dissimilar domains and that disorder-induced backscattering is substantially suppressed due to the weak coupling between oppositely valley-polarized modes. The coupling between valley modes is also further investigated and linked to an evident chiral flux of the mechanical energy. Finally, we show that the weak coupling between the valleys can be exploited to achieve selective mode injection at the domain wall, hence realizing a very effective excitation strategy of the chiral edge states.


title: Some overdetermined problems related to the anisotropic capacity

abstract: We characterize the Wulff shape of an anisotropic norm in terms of solutions to overdetermined problems for the Finsler $p$-capacity of a convex set $\Omega \subset \mathbb{R}^N$, with $1<p<N$. In particular we show that if the Finsler $p$-capacitary potential $u$ associated to $\Omega$ has two homothetic level sets then $\Omega$ is Wulff shape. Moreover, we show that the concavity exponent of $u$ is $q=-(p-1)/(N-p)$ if and only if $\Omega$ is Wulff shape.


title: Data-Driven Distributionally Robust Optimal Power Flow for Distribution Systems

abstract: Increasing penetration of distributed energy resources complicate operations of electric power distribution systems by amplifying volatility of nodal power injections. On the other hand, these resources can provide additional control means to the distribution system operator (DSO). This paper takes the DSO perspective and leverages a data-driven distributionally robust decision-making framework to overcome the uncertainty of these injections and its impact on the distribution system operations. We develop an AC OPF formulation for radial distribution systems based on the LinDistFlow AC power flow approximation and exploit distributionally robust optimization to immunize the optimized decisions against uncertainty in the probabilistic models of forecast errors obtained from the available observations. The model is reformulated to be computationally tractable and tested on multiple IEEE distribution test systems. We also release the code supplement that implements the proposed model in Julia and can be used to reproduce our numerical results.


title: On non-commutativity in quantum theory (I): from classical to quantum probability

abstract: A central feature of quantum mechanics is the non-commutativity of operators used to describe physical observables. In this article, we present a critical analysis on the role of non-commutativity in quantum theory, focusing on its consequences in the probabilistic description. Typically, a random phenomenon is described using the measure-theoretic formulation of probability theory. Such a description can also be done using algebraic methods, which are capable to deal with non-commutative random variables (like in quantum mechanics). Here we propose a method to construct a non-commutative probability theory starting from an ordinary measure-theoretic description of probability. This will be done using the entropic uncertainty relations between random variables, in order to evaluate the presence of non-commutativity in their algebraic description.


title: A probabilistic generalization of the Stirling numbers of the second kind

abstract: Associated to each random variable $Y$ having a finite moment generating function, we introduce a different generalization of the Stirling numbers of the second kind. Some characterizations and specific examples of such generalized numbers are provided. As far as their applications are concerned, attention is focused in extending in various ways the classical formula for sums of powers on arithmetic progressions. Illustrations involving rising factorials, Bell polynomials, polylogarithms, and a certain class of Appell polynomials, in connection with appropriate random variables $Y$ in each case, are discussed in detail.


title: Nutritionally recommended food for semi- to strict vegetarian diets based on large-scale nutrient composition data

abstract: Diet design for vegetarian health is challenging due to the limited food repertoire of vegetarians. This challenge can be partially overcome by quantitative, data-driven approaches that utilise massive nutritional information collected for many different foods. Based on large-scale data of foods' nutrient compositions, the recent concept of nutritional fitness helps quantify a nutrient balance within each food with regard to satisfying daily nutritional requirements. Nutritional fitness offers prioritisation of recommended foods using the foods' occurrence in nutritionally adequate food combinations. Here, we systematically identify nutritionally recommendable foods for semi- to strict vegetarian diets through the computation of nutritional fitness. Along with commonly recommendable foods across different diets, our analysis reveals favourable foods specific to each diet, such as immature lima beans for a vegan diet as an amino acid and choline source, and mushrooms for ovo-lacto vegetarian and vegan diets as a vitamin D source. Furthermore, we find that selenium and other essential micronutrients can be subject to deficiency in plant-based diets, and suggest nutritionally-desirable dietary patterns. We extend our analysis to two hypothetical scenarios of highly personalised, plant-based methionine-restricted diets. Our nutrient-profiling approach may provide a useful guide for designing different types of personalised vegetarian diets.


title: On non-commutativity in quantum theory (II): toy models for non-commutative kinematics

abstract: In this article, we continue our investigation on the role of non-commutativity in quantum theory. Using the method explained in "On non-commutativity in quantum theory (I): from classical to quantum probability", we analyze two toy models which exhibit non-commutativity between the corresponding position and velocity random variables. In particular, using ordinary probability theory, we study the kinematics of a point-like particle jumping at random over a discrete random space. We show that, after the removal of the random space from the model, the position and velocity of the particle do not commute, when represented as operators on the same Hilbert space.


title: Optical potential mapping with a levitated nanoparticle at sub-wavelength distances from a membrane

abstract: The controllable positioning of a vacuum-levitated object near a material surface is of importance for studying short-range forces, such as Casimir forces, interfacial friction forces, or gravity in yet unexplored parameter regimes. Here we optically levitate a nanoparticle in a laser beam strongly focused on a dielectric membrane. By investigating the motion of the trapped particle in vacuum, we map the position-dependent optical potential of the particle. We interferometrically measure the distance between the particle and the surface and demonstrate stable trapping in sub-wavelength proximity of the dielectric surface. Our work is important for the development of on-chip levitated optomechanics and for measuring short-range forces at sub-wavelength distances.


title: Uniform phases in fluids of hard isosceles triangles: one component and binary mixtures

abstract: We formulate the scaled particle theory for a general mixture of hard isosceles triangles and calculate different phase diagrams for the one-component fluid and for certain binary mixtures. The fluid of hard triangles exhibits a complex phase behavior: (i) the presence of a triatic phase with sixfold symmetry, (ii) the isotropic-uniaxial nematic transition is of first order for certain ranges of aspect ratios, and (iii) the one-component system exhibits nematic-nematic transitions ending in critical points. We found the triatic phase to be stable not only for equilateral triangles but also for triangles of similar aspect ratios. We focus the study of binary mixtures on the case of symmetric mixtures: equal particle areas with aspect ratios ($\kappa_i$) symmetric with respect to the equilateral one: $\kappa_1\kappa_2=3$. For these mixtures we found, aside from first-order isotropic-nematic and nematic-nematic transitions (the latter ending in a critical point): (i) A region of triatic phase stability even for mixtures made of particles that do not form this phase at the one-component limit, and (ii) the presence of a Landau point at which two isotropic-nematic first-order transitions and a nematic-nematic demixing transition coalesce. This phase behavior is analog to that of a symmetric three-dimensional mixture of rods and plates.


title: Implications of Neutron Star Mergers for Extraterrestrial Civilizations

abstract: The economy and fate of extraterrestrial civilizations should depend on the abundance of gold and uranium, made in neutron star mergers.


title: Unusual suppression of the superconducting energy gap and critical temperature in atomically thin NbSe2

abstract: It is well known that superconductivity in thin films is generally suppressed with decreasing thickness. This suppression is normally governed by either disorder-induced localization of Cooper pairs, weakening of Coulomb screening, or generation and unbinding of vortex-antivortex pairs as described by the Berezinskii-Kosterlitz-Thouless (BKT) theory. Defying general expectations, few-layer NbSe2 - an archetypal example of ultrathin superconductors - has been found to remain superconducting down to monolayer thickness. Here we report measurements of both the superconducting energy gap and critical temperature in high-quality monocrystals of few-layer NbSe2, using planar-junction tunneling spectroscopy and lateral transport. We observe a fully developed gap that rapidly reduces for devices with the number of layers N < 5, as does their ctitical temperature. We show that the observed reduction cannot be explained by disorder, and the BKT mechanism is also excluded by measuring its transition temperature that for all N remains very close to Tc. We attribute the observed behavior to changes in the electronic band structure predicted for mono- and bi- layer NbSe2 combined with inevitable suppression of the Cooper pair density at the superconductor-vacuum interface. Our experimental results for N > 2 are in good agreement with the dependences of the gap and Tc expected in the latter case while the effect of band-structure reconstruction is evidenced by a stronger suppression of the gap and the disappearance of its anisotropy for N = 2. The spatial scale involved in the surface suppression of the density of states is only a few angstroms but cannot be ignored for atomically thin superconductors.


title: On non-commutativity in quantum theory (III): determinantal point processes and non-relativistic quantum mechanics

abstract: This article concludes our critical analysis on the role of non-commutativity in quantum theory. After a brief introduction of the necessary notions on point processes, we re-analyse model B proposed in "On non-commutativity in quantum theory (II): toy models for non-commutative kinematics", using the point process theory. This viewpoint allows to generalize and modify the space process of model B in a simple manner, obtaining a new model (model C). This new model allows the recovery of non-relativistic quantum mechanics in a suitable limit.


title: A simple cure for numerical shock instability in HLLC Riemann solver

abstract: The Harten-Lax-van Leer with contact (HLLC) scheme is known to be plagued by various forms of numerical shock instabilities. In this paper, we propose a new framework for developing shock stable, contact and shear preserving approximate Riemann solvers based on the HLLC scheme for the Euler system of equations. The proposed framework termed as HLLC-SWM (\textbf{S}elective \textbf{W}ave \textbf{M}odified) scheme identifies and increases the magnitude of the inherent diffusive HLL component within the HLLC scheme in the vicinity of a shock wave while leaving its antidiffusive component unmodified to retain accuracy on linearly degenerate wavefields. We present two strategies to compute the requisite supplementary dissipation which results in HLLC-SWM-E and HLLC-SWM-P variants. Through a linear perturbation analysis of the HLLC-SWM framework, we clarify how the additional dissipation introduced in this way helps in damping of unphysical perturbations in primitive quantities under a derived CFL constraint. A matrix based stability analysis of a steady two-dimensional normal shock is used to show that both variants of the HLLC-SWM scheme are shock stable over a wide range of inlet Mach numbers. Results from standard test cases demonstrate that the HLLC-SWM schemes are capable of computing shock stable solutions on a variety of problems while retaining positivity and exact inviscid contact ability. On viscous flows, while the HLLC-SWM-P variant is quite accurate, the HLLC-SWM-E variant introduces slight inaccuracy which can be corrected through a simple Mach number based switching function.


title: Connectivity properties of the adjacency graph of SLE$_\kappa$ bubbles for $\kappa \in (4,8)$

abstract: We study the adjacency graph of bubbles---i.e., complementary connected components---of an SLE$_{\kappa}$ curve for $\kappa \in (4,8)$, with two such bubbles considered to be adjacent if their boundaries intersect. We show that this adjacency graph is a.s. connected for $\kappa \in (4,\kappa_0]$, where $\kappa_0 \approx 5.6158$ is defined explicitly. This gives a partial answer to a problem posed by Duplantier, Miller and Sheffield (2014). Our proof in fact yields a stronger connectivity result for $\kappa \in (4,\kappa_0]$, which says that there is a Markovian way of finding a path from any fixed bubble to $\infty$. We also show that there is a (non-explicit) $\kappa_1 \in (\kappa_0, 8)$ such that this stronger condition does not hold for $\kappa \in [\kappa_1,8)$. Our proofs are based on an encoding of SLE$_\kappa$ in terms of a pair of independent $\kappa/4$-stable processes, which allows us to reduce our problem to a problem about stable processes. In fact, due to this encoding, our results can be re-phrased as statements about the connectivity of the adjacency graph of loops when one glues together an independent pair of so-called $\kappa/4$-stable looptrees, as studied, e.g., by Curien and Kortchemski (2014). The above encoding comes from the theory of Liouville quantum gravity (LQG), but the paper can be read without any knowledge of LQG if one takes the encoding as a black box.


title: Dense Limit of the Dawid-Skene Model for Crowdsourcing and Regions of Sub-optimality of Message Passing Algorithms

abstract: Crowdsourcing is a strategy to categorize data through the contribution of many individuals. A wide range of theoretical and algorithmic contributions are based on the model of Dawid and Skene [1]. Recently it was shown in [2,3] that, in certain regimes, belief propagation is asymptotically optimal for data generated from the Dawid-Skene model. This paper is motivated by this recent progress. We analyze the dense limit of the Dawid-Skene model. It is shown that it belongs to a larger class of low-rank matrix estimation problems for which it is possible to express the asymptotic, Bayes-optimal, performance in a simple closed form. In the dense limit the mapping to a low-rank matrix estimation problem provides an approximate message passing algorithm that solves the problem algorithmically. We identify the regions where the algorithm efficiently computes the Bayes-optimal estimates. Our analysis refines the results of [2,3] about optimality of message passing algorithms by characterizing regions of parameters where these algorithms do not match the Bayes-optimal performance. We further study numerically the performance of approximate message passing, derived in the dense limit, on sparse instances and carry out experiments on a real world dataset.

categories: stat.ML cond-mat.dis-nn cond-mat.stat-mech physics.data-an


title: Ultra-efficient superconducting Dayem bridge field-effect transistor

abstract: Superconducting field-effect transitor (SuFET) and Josephson field-effect transistor (JoFET) technologies take advantage of electric field induced control of charge carrier concentration in order to modulate the channel superconducting properties. Despite field-effect is believed to be unaffective for superconducting metals, recent experiments showed electric field dependent modulation of the critical current (IC) in a fully metallic transistor. Yet, the grounding mechanism of this phenomenon is not completely understood. Here, we show the experimental realization of Ti-based Dayem bridge field-effect transistors (DB-FETs) able to control IC of the superconducting channel. Our easy fabrication process DB-FETs show symmetric full suppression of IC for an applied critical gate voltage as low as VCG~+-8V at temperatures reaching about the 85% of the record critical temperature TC~550mK for titanium. The gate-independent TC and normal state resistance (RN) coupled with the increase of resistance in the supercoducting state (RS) for gate voltages close to the critical value (VCG) suggest the creation of field-effect induced metallic puddles in the superconducting sea. Our devices show extremely high values of transconductance (gMAXm~15uA/V at VG~+-6.5V) and variations of Josephson kinetic inductance (LK) with VG of two orders of magnitude. Therefore, the DB-FET appears as an ideal candidate for the realization of superconducting electronics, superconducting qubits, tunable interferometers as well as photon detectors.


title: Active Reinforcement Learning with Monte-Carlo Tree Search

abstract: Active Reinforcement Learning (ARL) is a twist on RL where the agent observes reward information only if it pays a cost. This subtle change makes exploration substantially more challenging. Powerful principles in RL like optimism, Thompson sampling, and random exploration do not help with ARL. We relate ARL in tabular environments to Bayes-Adaptive MDPs. We provide an ARL algorithm using Monte-Carlo Tree Search that is asymptotically Bayes optimal. Experimentally, this algorithm is near-optimal on small Bandit problems and MDPs. On larger MDPs it outperforms a Q-learner augmented with specialised heuristics for ARL. By analysing exploration behaviour in detail, we uncover obstacles to scaling up simulation-based algorithms for ARL.


title: An Agent-Based Simulation of Residential Location Choice of Tenants in Tehran, Iran

abstract: Residential location choice modeling is one of the substantial components of land use and transportation models. While numerous aggregated mathematical and statistical approaches have been developed to model the residence choice behavior of households, disaggregated approaches such as the agent-based modeling have shown interesting capabilities. In this article, a novel agent-based approach is developed to simulate the residential location choice of tenants in Tehran, the capital of Iran. Tenants are considered as agents who select their desired residential alternatives according to their characteristics and preferences for various criteria such as the rent, accessibility to different services and facilities, environmental pollution, and distance from their workplace and former residence. The choice set of agents is limited to their desired residential alternatives by applying a constrained NSGA-II algorithm. Then, agents compete with each other to select their final residence among their alternatives. Results of the proposed approach are validated by comparing simulated and actual residences of a sample of tenants. Results show that the proposed approach is able to accurately simulate the residence of 59.3% of tenants at the traffic analysis zone level.


title: Post-recombination Dark Matter for the 21-cm Signal

abstract: We have no certain knowledge of the early history of dark matter (DM). In this paper we propose a scenario where DM is produced post-recombination but prior to the cosmic dawn. It helps to relax the bounds on DM interactions, in particular with baryons, from the CMB. It may be of interest in some circumstances, for example, to understand the recent cosmic dawn 21-cm signal anomaly. We argue that the cosmic gas cooling mechanism via the minicharged DM-baryon scattering may be viable even if it takes up the total DM budget. We also investigate the possibility of a gluon-philic mediator of a few 10 keV, to find that the most reliable exclusion is from the neutron scattering.


title: Structural Agnostic Modeling: Adversarial Learning of Causal Graphs

abstract: A new causal discovery method, Structural Agnostic Modeling (SAM), is presented in this paper. Leveraging both conditional independencies and distributional asymmetries, SAM aims to find the underlying causal structure from observational data. The approach is based on a game between different players estimating each variable distribution conditionally to the others as a neural net, and an adversary aimed at discriminating the generated data against the original data. A learning criterion combining distribution estimation, sparsity and acyclicity constraints is used to enforce the optimization of the graph structure and parameters through stochastic gradient descent. SAM is extensively experimentally validated on synthetic and real data.


title: Polynomials and degrees of maps in real normed algebras

abstract: Let $\cal{A}$ be the algebra of quaternions $\mathbb{H}$ or octonions $\mathbb{O}$. In this manuscript a new proof is given, based on ideas of Cauchy and D' Alembert, of the fact that an ordinary polynomial $f(t) \in {\cal{A}}\, [t]$ has a root in $\cal{A}$. As a consequence, the Jacobian determinant $|J(f)|$ is always non negative in $\cal{A}$. Moreover, using the idea of the topological degree we show that a regular polynomial $g(t)$ over $\cal{A}$ has also a root in $\cal{A}$. Finally, utilizing multiplication $(*)$ in $\cal{A}$, we prove various results on the topological degree of products of maps. In particular, if $S$ is the unit sphere in $\cal{A}$ and $h_1, h_2: S \to S$ are smooth maps, it is shown that $\hbox{deg} (h_1 * h_2)=\hbox{deg} (h_1) + \hbox{deg} (h_2)$.


title: A polynomial ideal associated to any $t$-$(v,k,\lambda)$ design

abstract: We consider ordered pairs $(X,\mathcal{B})$ where $X$ is a finite set of size $v$ and $\mathcal{B}$ is some collection of $k$-element subsets of $X$ such that every $t$-element subset of $X$ is contained in exactly $\lambda$ "blocks" $B\in \mathcal{B}$ for some fixed $\lambda$. We represent each block $B$ by a zero-one vector $\mathbf{c}_B$ of length $v$ and explore the ideal $\mathcal{I}(\mathcal{B})$ of polynomials in $v$ variables with complex coefficients which vanish on the set $\{ \mathbf{c}_B \mid B \in \mathcal{B}\}$. After setting up the basic theory, we investigate two parameters related to this ideal: $\gamma_1(\mathcal{B})$ is the smallest degree of a non-trivial polynomial in the ideal $\mathcal{I}(\mathcal{B})$ and $\gamma_2(\mathcal{B})$ is the smallest integer $s$ such that $\mathcal{I}(\mathcal{B})$ is generated by a set of polynomials of degree at most $s$. We first prove the general bounds $t/2 < \gamma_1(\mathcal{B}) \le \gamma_2(\mathcal{B}) \le k$. Examining important families of examples, we find that, for symmetric 2-designs and Steiner systems, we have $\gamma_2(\mathcal{B}) \le t$. But we expect $\gamma_2(\mathcal{B})$ to be closer to $k$ for less structured designs and we indicate this by constructing infinitely many triple systems satisfying $\gamma_2(\mathcal{B})=k$.


title: Impacts of transport development on residence choice of renter households: An agent-based evaluation

abstract: Because of improving accessibility, transport developments play an important role in residence choice of renter households. In this paper, an agent-based model is developed to investigate impacts of different transport developments on residence choice of renter households in Tehran, the capital of Iran. In the proposed model, renter households are considered as agents who make a multi-objective decision and compete with each other to rent a preferred residential zone. Then, three transport development scenarios including construction a new highway, subway and bus rapid transit (BRT) line are simulated and resulting changes in residence choice of agents are evaluated. Results show that transport development scenarios significantly affect residence choice behavior of different socio-economic categories of renter households and lead to considerable changes in the residential demand, composition of residents, mean income level and mean car ownership in their vicinities.


title: Approximate Quantum Fourier Transform with $O(n \log(n))$ T gates

abstract: The ability to implement the Quantum Fourier Transform (QFT) efficiently on a quantum computer facilitates the advantages offered by a variety of fundamental quantum algorithms, such as those for integer factoring, computing discrete logarithm over Abelian groups, solving systems of linear equations, and phase estimation, to name a few. The standard fault-tolerant implementation of an $n$-qubit unitary QFT approximates the desired transformation by removing small-angle controlled rotations and synthesizing the remaining ones into Clifford+T gates, incurring the T-count complexity of $O(n \log^2(n))$. In this paper, we show how to obtain approximate QFT with the T-count of $O(n \log(n))$. Our approach relies on quantum circuits with measurements and feedforward, and on reusing a special quantum state that induces the phase gradient transformation. We report asymptotic analysis as well as concrete circuits, demonstrating significant advantages in both theory and practice.