title: Spatial Mapping with Gaussian Processes and Nonstationary Fourier Features

abstract: The use of covariance kernels is ubiquitous in the field of spatial statistics. Kernels allow data to be mapped into high-dimensional feature spaces and can thus extend simple linear additive methods to nonlinear methods with higher order interactions. However, until recently, there has been a strong reliance on a limited class of stationary kernels such as the Matern or squared exponential, limiting the expressiveness of these modelling approaches. Recent machine learning research has focused on spectral representations to model arbitrary stationary kernels and introduced more general representations that include classes of nonstationary kernels. In this paper, we exploit the connections between Fourier feature representations, Gaussian processes and neural networks to generalise previous approaches and develop a simple and efficient framework to learn arbitrarily complex nonstationary kernel functions directly from the data, while taking care to avoid overfitting using state-of-the-art methods from deep learning. We highlight the very broad array of kernel classes that could be created within this framework. We apply this to a time series dataset and a remote sensing problem involving land surface temperature in Eastern Africa. We show that without increasing the computational or storage complexity, nonstationary kernels can be used to improve generalisation performance and provide more interpretable results.


title: Strain gradient plasticity modeling of hydrogen diffusion to the crack tip

abstract: In this work hydrogen diffusion towards the fracture process zone is examined accounting for local hardening due to geometrically necessary dislocations (GNDs) by means of strain gradient plasticity (SGP). Finite element computations are performed within the finite deformation theory to characterize the gradient-enhanced stress elevation and subsequent diffusion of hydrogen towards the crack tip. Results reveal that GNDs, absent in conventional plasticity predictions, play a fundamental role on hydrogen transport ahead of a crack. SGP estimations provide a good agreement with experimental measurements of crack tip deformation and high levels of lattice hydrogen concentration are predicted within microns to the crack tip. The important implications of the results in the understanding of hydrogen embrittlement mechanisms are thoroughly discussed.


title: Reduction of electron heating by magnetizing ultracold neutral plasma

abstract: Electron heating in an ultracold neutral plasma is modeled using classical molecular dynamics simulations in the presence of an externally applied magnetic field. A sufficiently strong magnetic field is found to reduce disorder induced heating and three body recombination heating of electrons by constraining electron motion, and therefore heating, to the single dimension aligned with the magnetic field. A strong and long-lasting temperature anisotropy develops, and the overall kinetic electron temperature is effectively reduced by a factor of three. These results suggest that experiments may increase the effective electron coupling strength using an applied magnetic field.


title: P-spline smoothing for spatial data collected worldwide

abstract: Spatial data collected worldwide at a huge number of locations are frequently used in environmental and climate studies. Spatial modelling for this type of data presents both methodological and computational challenges. In this work we illustrate a computationally efficient non parametric framework to model and estimate the spatial field while accounting for geodesic distances between locations. The spatial field is modelled via penalized splines (P-splines) using intrinsic Gaussian Markov Random Field (GMRF) priors for the spline coefficients. The key idea is to use the sphere as a surrogate for the Globe, then build the basis of B-spline functions on a geodesic grid system. The basis matrix is sparse and so is the precision matrix of the GMRF prior, thus computational efficiency is gained by construction. We illustrate the approach on a real climate study, where the goal is to identify the Intertropical Convergence Zone using high-resolution remote sensing data.


title: An improved lower bound for the maximal length of a multivector

abstract: A new lower bound for the maximal length of a multivector is obtained. It is much closer to the best known upper bound than previously reported lower bound estimates. The maximal length appears to be unexpectedly large for $n$-vectors, with n>2, since the few exactly known values seem to grow linearly with vector space dimension, whereas the new lower bound has a polynomial order equal to n-1 like the best known upper bound. This result has implications for quantum chemistry.


title: Sharp non-asymptotic Concentration Inequalities for the Approximation of the Invariant Measure of a Diffusion

abstract: For an ergodic Brownian diffusion with invariant measure $\nu$, we consider a sequence of empirical distributions ($\nu$n) n$\ge$1 associated with an approximation scheme with decreasing time step ($\gamma$n) n$\ge$1 along an adapted regular enough class of test functions f such that f --$\nu$(f) is a coboundary of the infinitesimal generator A. Denote by $\sigma$ the diffusion coefficient and $\Phi$ the solution of the Poisson equation A$\Phi$ = f -- $\nu$(f). When the square norm of |$\sigma$ * $\Phi$| 2 lies in the same coboundary class as f , we establish sharp non-asymptotic concentration bounds for suitable normalizations of $\nu$n(f) -- $\nu$(f). Our bounds are optimal in the sense that they match the asymptotic limit obtained by Lamberton and Pag{\`e}s in [LP02], for a certain large deviation regime. In particular, this allows us to derive sharp non-asymptotic confidence intervals. We provide as well a Slutsky like Theorem, for practical applications, where the deviation bounds are also asymptotically independent of the corresponding Poisson problem. Eventually, we are able to handle, up to an additional constraint on the time steps, Lipschitz sources f in an appropriate non-degenerate setting.


title: Half-lives of neutron-rich $^{128-130}$Cd

abstract: The $\beta$-decay half-lives of $^{128\text{--}130}$Cd have been measured with the newly commissioned GRIFFIN $\gamma$-ray spectrometer at the TRIUMF-ISAC facility. The time structures of the most intense $\gamma$-rays emitted following the $\beta$-decay were used to determine the half-lives of $^{128}$Cd and $^{130}$Cd to be $T_{1/2}= 246.2(21)$ ms and $T_{1/2}= 126(4)$ ms, respectively. The half-lives of the 3/2$^+$ and 11/2$^-$ states of $^{129}$Cd were measured to be $T_{1/2}(3/2^+)= 157(8)$ ms and $T_{1/2}(11/2^-)= 147(3)$ ms. The half-lives of the Cd isotopes around the $N=82$ shell closure are an important ingredient in astrophysical simulations to derive the magnitude of the second $r$-process abundance peak in the $A\sim130$ region. Our new results are compared with recent literature values and theoretical calculations.


title: Exploring triplet-quadruplet fermionic dark matter at the LHC and future colliders

abstract: We study the signatures of the triplet-quadruplet dark matter model at the LHC and future colliders, including the 100 TeV Super Proton-Proton Collider and the 240 GeV Circular Electron Positron Collider. The dark sector in this model contains one fermionic electroweak triplet and two fermionic quadruplets, which have two kinds of Yukawa couplings to the Higgs doublet. Electroweak production signals of the dark sector fermions in the $\text{monojet}+/\!\!\!\!E_\mathrm{T}$, disappearing track, and $\text{multilepton}+/\!\!\!\!E_\mathrm{T}$ channels at the LHC and the Super Proton-Proton Collider are investigated. Moreover, we study the loop effects of this model on the Circular Electron Positron Collider precision measurements of $e^+e^-\to Zh$ and $h\to\gamma\gamma$. We find that most of the parameter regions allowed by the observed dark matter relic density will be well explored by such direct and indirect searches at future colliders.


title: Measurement of the ratio of branching fractions $\mathcal{B}(B_c^+\,\to\,J/\psi\tau^+\nu_\tau)$/$\mathcal{B}(B_c^+\,\to\,J/\psi\mu^+\nu_\mu)$

abstract: A measurement is reported of the ratio of branching fractions $\mathcal{R}(J/\psi)=\mathcal{B}(B_c^+\,\to\,J/\psi\tau^+\nu_\tau)/\mathcal{B}(B_c^+\,\to\,J/\psi \mu^+\nu_\mu)$, where the $\tau^+$ lepton is identified in the decay mode $\tau^+\,\to\,\mu^+\nu_\mu\overline{\nu}_\tau$. This analysis uses a sample of proton-proton collision data corresponding to 3.0$\mathrm{\,fb}^{-1}$ of integrated luminosity recorded with the LHCb experiment at center-of-mass energies $7\,\mathrm{TeV}$ and $8\,\mathrm{TeV}$. A signal is found for the decay $B_c^+\,\to\,J/\psi \tau^+\nu_\tau$ at a significance of 3 standard deviations, corrected for systematic uncertainty, and the ratio of the branching fractions is measured to be $\mathcal{R}(J/\psi) = 0.71 \pm 0.17 \mathrm{\,(stat)} \pm 0.18\mathrm{\,(syst)}$. This result lies within 2 standard deviations above the range of existing predictions in the Standard Model.


title: Gaussian width bounds with applications to arithmetic progressions in random settings

abstract: Motivated by problems on random differences in Szemer\'{e}di's theorem and on large deviations for arithmetic progressions in random sets, we prove upper bounds on the Gaussian width of point sets that are formed by the image of the $n$-dimensional Boolean hypercube under a mapping $\psi:\mathbb{R}^n\to\mathbb{R}^k$, where each coordinate is a constant-degree multilinear polynomial with 0-1 coefficients. We show the following applications of our bounds. Let $[\mathbb{Z}/N\mathbb{Z}]_p$ be the random subset of $\mathbb{Z}/N\mathbb{Z}$ containing each element independently with probability $p$. $\bullet$ A set $D\subseteq \mathbb{Z}/N\mathbb{Z}$ is $\ell$-intersective if any dense subset of $\mathbb{Z}/N\mathbb{Z}$ contains a proper $(\ell+1)$-term arithmetic progression with common difference in $D$. Our main result implies that $[\mathbb{Z}/N\mathbb{Z}]_p$ is $\ell$-intersective with probability $1 - o(1)$ provided $p \geq \omega(N^{-\beta_\ell}\log N)$ for $\beta_\ell = (\lceil(\ell+1)/2\rceil)^{-1}$. This gives a polynomial improvement for all $\ell \ge 3$ of a previous bound due to Frantzikinakis, Lesigne and Wierdl, and reproves more directly the same improvement shown recently by the authors and Dvir. $\bullet$ Let $X_k$ be the number of $k$-term arithmetic progressions in $[\mathbb{Z}/N\mathbb{Z}]_p$ and consider the large deviation rate $\rho_k(\delta) = \log\Pr[X_k \geq (1+\delta)\mathbb{E}X_k]$. We give quadratic improvements of the best-known range of $p$ for which a highly precise estimate of $\rho_k(\delta)$ due to Bhattacharya, Ganguly, Shao and Zhao is valid for all odd $k \geq 5$. We also discuss connections with error correcting codes (locally decodable codes) and the Banach-space notion of type for injective tensor products of $\ell_p$-spaces.


title: Charmonia production from $b$-hadron decays at LHC with $k_T$-factorization: $J/\psi$, $\psi(2S)$ and $J/\psi + Z$

abstract: We consider the production of $J/\psi$ and $\psi(2S)$ mesons originating from the decays of $b$-flavored hadrons at the LHC using the $k_T$-factorization approach. Our analysis covers both inclusive charmonia production and production of $J/\psi$ mesons in association with $Z$ bosons. We apply the transverse momentum dependent (or unintegrated) gluon density in a proton derived from Catani-Ciafaloni-Fiorani-Marchesini (CCFM) evolution equation and adopt fragmentation functions based on the non-relativistic QCD factorization to describe the inclusive $b$-hadron decays into the different charmonium states. Our predictions agree well with latest experimental data taken by the CMS, ATLAS and LHCb Collaborations at $\sqrt s = 7$, $8$ and $13$~TeV. The contributions from double parton scattering to the associated non-prompt $J/\psi + Z$ production are estimated and found to be small.


title: Deep Temporal-Recurrent-Replicated-Softmax for Topical Trends over Time

abstract: Dynamic topic modeling facilitates the identification of topical trends over time in temporal collections of unstructured documents. We introduce a novel unsupervised neural dynamic topic model named as Recurrent Neural Network-Replicated Softmax Model (RNNRSM), where the discovered topics at each time influence the topic discovery in the subsequent time steps. We account for the temporal ordering of documents by explicitly modeling a joint distribution of latent topical dependencies over time, using distributional estimators with temporal recurrent connections. Applying RNN-RSM to 19 years of articles on NLP research, we demonstrate that compared to state-of-the art topic models, RNNRSM shows better generalization, topic interpretation, evolution and trends. We also introduce a metric (named as SPAN) to quantify the capability of dynamic topic model to capture word evolution in topics over time.


title: Exploiting Layerwise Convexity of Rectifier Networks with Sign Constrained Weights

abstract: By introducing sign constraints on the weights, this paper proposes sign constrained rectifier networks (SCRNs), whose training can be solved efficiently by the well known majorization-minimization (MM) algorithms. We prove that the proposed two-hidden-layer SCRNs, which exhibit negative weights in the second hidden layer and negative weights in the output layer, are capable of separating any two (or more) disjoint pattern sets. Furthermore, the proposed two-hidden-layer SCRNs can decompose the patterns of each class into several clusters so that each cluster is convexly separable from all the patterns from the other classes. This provides a means to learn the pattern structures and analyse the discriminant factors between different classes of patterns.


title: Infinite order $\Psi$DOs: Composition with entire functions, new Shubin-Sobolev spaces, and index theorem

abstract: We study global regularity and spectral properties of power series of the Weyl quantisation $a^w$, where $a(x,\xi) $ is a classical elliptic Shubin polynomial. For a suitable entire function $P$, we associate two natural infinite order operators to $a^{w}$, $P(a^w)$ and $(P\circ a)^{w},$ and prove that these operators and their lower order perturbations are globally Gelfand-Shilov regular. They have spectra consisting of real isolated eigenvalues diverging to $\infty$ for which we find the asymptotic behaviour of their eigenvalue counting function. In the second part of the article, we introduce Shubin-Sobolev type spaces by means of $f$-$\Gamma^{*,\infty}_{A_p,\rho}$-elliptic symbols, where $f $ is a function of ultrapolynomial growth and $\Gamma^{*,\infty}_{A_p,\rho}$ is a class of symbols of infinite order studied in this and our previous papers. We study the regularity properties of these spaces, and show that the pseudo-differential operators under consideration are Fredholm operators on them. Their indices are independent on the order of the Shubin-Sobolev spaces; finally, we show that the index can be expressed via a Fedosov-H\"{o}rmander integral formula.


title: An example of explicit implementation strategy and preconditioning for the high order edge finite elements applied to the time-harmonic Maxwell's equations

abstract: In this paper we focus on high order finite element approximations of the electric field combined with suitable preconditioners, to solve the time-harmonic Maxwell's equations in waveguide configurations.The implementation of high order curl-conforming finite elements is quite delicate, especially in the three-dimensional case. Here, we explicitly describe an implementation strategy, which has been embedded in the open source finite element software FreeFem++ (http://www.freefem.org/ff++/). In particular, we use the inverse of a generalized Vandermonde matrix to build basis functions in duality with the degrees of freedom, resulting in an easy-to-use but powerful interpolation operator. We carefully address the problem of applying the same Vandermonde matrix to possibly differently oriented tetrahedra of the mesh over the computational domain. We investigate the preconditioning for Maxwell's equations in the time-harmonic regime, which is an underdeveloped issue in the literature, particularly for high order discretizations. In the numerical experiments, we study the effect of varying several parameters on the spectrum of the matrix preconditioned with overlapping Schwarz methods, both for 2d and 3d waveguide configurations.


title: Mathematical modeling of Zika disease in pregnant women and newborns with microcephaly in Brazil

abstract: We propose a new mathematical model for the spread of Zika virus. Special attention is paid to the transmission of microcephaly. Numerical simulations show the accuracy of the model with respect to the Zika outbreak occurred in Brazil.


title: Light cluster production at NICA

abstract: In a recent publication the importance of measuaring light clusters at future NICA and FAIR experiments as "rare probes" for in-medium characteristics was presented. Within theoretical conciderations is shown how clusters would be affected by medium affects as given collider energies.


title: Gaussian Parsimonious Clustering Models with Covariates and a Noise Component

abstract: We consider model-based clustering methods for continuous, correlated data that account for external information available in the presence of mixed-type fixed covariates by proposing the MoEClust suite of models. These models allow different subsets of covariates to influence the component weights and/or component densities by modelling the parameters of the mixture as functions of the covariates. A familiar range of constrained eigen-decomposition parameterisations of the component covariance matrices are also accommodated. This paper thus addresses the equivalent aims of including covariates in Gaussian parsimonious clustering models and incorporating parsimonious covariance structures into all special cases of the Gaussian mixture of experts framework. The MoEClust models demonstrate significant improvement from both perspectives in applications to both univariate and multivariate data sets. Novel extensions to include a uniform noise component for capturing outliers and to address initialisation of the EM algorithm, model selection, and the visualisation of results are also proposed.


title: Sub-cycle dynamics of electron-hole pairs and high-harmonic generation in bulk diamond subjected to intense femtosecond laser pulse

abstract: We present calculation of photoexcitation and high-harmonic generation in bulk diamond induced by intense near-infrared laser pulse of photon energy 1.55 eV, time duration of 15 fs and peak field strength $F=$ 0.4 V/{\AA}. Depending on the laser polarization direction, the pulsed irradiation creates electron-hole pairs and transient density fluctuations in a localized region of the crystal momentum space. As a consequence energetic inter- and intra-band harmonics are generated with alternating phase during each half-cycle of the driving pulse. The corresponding inter- and intra-band currents are in definite phase relation with the laser pulse, reflecting the build-up of coherent superposition of population between valence and conduction bands, and the accelerated motion of charge carriers in their respective bands.


title: Supressed splashing on elastic membranes

abstract: The dynamics of drop impact on solid surfaces can be changed significantly by tuning the elasticity of the solid. Most prominently, the substrate deformation causes an increase in the splashing threshold as compared to impact onto perfectly rigid surfaces, and can thus lead to splash suppression. Here, we experimentally determine the splashing threshold for impact on thin membranes as a function of the tension in the membrane and its elastic properties. The drop dynamics is correlated to the membrane deformation, which is simultaneously measured using a laser profilometry technique. The experimental results enable us to adapt current models for splashing, showing quantitatively how substrate deformation alters the splashing threshold.


title: An example of how false conclusions could be made with personalized health tracking and suggestions for avoiding similar situations

abstract: Personalizing interventions and treatments is a necessity for optimal medical care. Recent advances in computing, such as personal electronic devices, have made it easier than ever to collect and utilize vast amounts of personal data on individuals. This data could support personalized medicine; however, there are pitfalls that must be avoided. We discuss an example, longitudinal medical tracking, in which traditional methods of evaluating machine learning algorithms fail and present the opportunity for false conclusions. We then pose three suggestions for avoiding such opportunities for misleading results in medical applications, where reliability is essential.


title: Galaxy and Mass Assembly (GAMA): Small-scale anisotropic galaxy clustering and the pairwise velocity dispersion of galaxies

abstract: The galaxy pairwise velocity dispersion (PVD) can provide important tests of non-standard gravity and galaxy formation models. We describe measurements of the PVD of galaxies in the Galaxy and Mass Assembly (GAMA) survey as a function of projected separation and galaxy luminosity. Due to the faint magnitude limit ($r < 19.8$) and highly-complete spectroscopic sampling of the GAMA survey, we are able to reliably measure the PVD to smaller scales ($r_\bot = 0.01$ Mpc/h) than previous work. The measured PVD at projected separations $r_\bot <~ 1$ Mpc/h increases near-monotonically with increasing luminosity from $\sigma \approx 200$ km/s at $M_r = -17$ mag to $\sigma \approx 600$ km/s at $M_r \approx -22$ mag. Analysis of the Gonzalez-Perez (2014) GALFORM semi-analytic model yields no such trend of PVD with luminosity: the model over-predicts the PVD for faint galaxies. This is most likely a result of the model placing too many low-luminosity galaxies in massive halos.


title: First observation of $B^{+} \to D_s^{+}K^{+}K^{-}$ decays and a search for $B^{+} \to D_s^{+}\phi$ decays

abstract: A search for $B^{+} \to D_s^{+}K^{+}K^{-}$ decays is performed using $pp$ collision data corresponding to an integrated luminosity of 4.8 fb$^{-1}$, collected at centre-of-mass energies of 7, 8 and 13$\,$TeV with the LHCb experiment. A significant signal is observed for the first time and the branching fraction is determined to be \begin{equation*} \mathcal{B}(B^{+} \to D_s^{+}K^{+}K^{-} ) = (7.1 \pm 0.5 \pm 0.6 \pm 0.7) \times 10^{-6}, \end{equation*} where the first uncertainty is statistical, the second systematic and the third due to the uncertainty on the branching fraction of the normalisation mode $B^{+} \to D_s^{+} \overline{\kern -0.2em D}^{0}$. A search is also performed for the pure annihilation decay $B^{+} \to D_s^{+}\phi$. No significant signal is observed and a limit of \begin{equation*} \mathcal{B}(B^{+} \to D_s^{+}\phi) < 4.9 \times 10^{-7}~(4.2 \times 10^{-7}) \end{equation*} is set on the branching fraction at 95$\%$ (90$\%$) confidence level.


title: Geometric aspects of interpolating gauge fixing in Chern-Simons Theory

abstract: In this article we investigate an interpolating gauge fixing procedure in $(4l+3)$-dimensional abelian Chern-Simons theory. We show that this interpolating gauge is related to the covariant gauge in a constant anisotropic metric. We compute the corresponding propagators involved in various expressions of the linking number in various gauges. We comment on the geometric interpretations of these expressions, clarifying how to pass from one interpretation to another.


title: Self-consistent assessment of Englert-Schwinger model on atomic properties

abstract: Our manuscript investigates a self-consistent solution of the statistical atom model proposed by Berthold-Georg Englert and Julian Schwinger (the ES model) and benchmarks it against atomic Kohn-Sham and two orbital-free models of the Thomas-Fermi-Dirac (TFD)-$\lambda$vW family. Results show that the ES model generally offers the same accuracy as the well-known TFD-$\frac{1}{5}$vW model; however, the ES model corrects the failure in Pauli potential near-nucleus region. We also point to the inability of describing low-$Z$ atoms as the foremost concern in improving the present model.


title: Contraction-based Observers using non-Euclidean Norms with an Application to Traffic Networks

abstract: In this note, we study Luenberger-type full-state observers for nonlinear systems using contraction theory. We show that if the matrix measure of a suitably defined Jacobian matrix constructed from the dynamics of the system-observer interconnection is uniformly negative, then the state estimate converges exponentially to the actual state. This sufficient condition for convergence establishes that the distance between the estimate and state is infinitesimally contracting with respect to some norm on the state-space. In contrast to existing results for contraction-based observer design, we allow for contraction with respect to non-Euclidean norms. Such norms have proven useful in applications. To demonstrate our results, we study the problem of observing vehicular traffic density along a freeway modeled as interconnected, spatially homogenous compartments, and our approach relies on establishing contraction of the system-observer interconnection with respect to the one-norm.


title: Monotonicity-based inversion of the fractional Schr\"odinger equation I. Positive potentials

abstract: We consider the inverse problems of for the fractional Schr\"odinger equation by using monotonicity formulas. We provide if-and-only-if monotonicity relations between positive bounded potentials and their associated nonlocal Dirichlet-to-Neumann maps. Based on the monotonicity relation, we can prove uniqueness for the nonlocal Calder\'on problem in a constructive manner. Secondly, we offer a reconstruction method for an unknown obstacles in a given domain. Our method is independent of the dimension $n\geq 2$ and only requires the background solution of the fractional Schr\"odinger equation.


title: Quantitative Benchmarks and New Directions for Noise Power Estimation Methods in ISM Radio Environment

abstract: Noise power estimation is a key issue in modern wireless communication systems. It allows resource allocation by detecting white spectral spaces effectively, and gives control over the communication process by adjusting transmission power. Thus far, the proposed estimation methods in the literature are based on spectral averaging, eigenvalues of sample covariance matrix, information theory, and statistical signal analysis. Each method is characterized by certain stability, accuracy and complexity. However, the existing literature does not provide an appropriate comparison. In this paper, we evaluate the performance of the existing estimation techniques intensively in terms of stability and accuracy, followed by detailed complexity analysis. The basis for comparison is signal-to-noise ratio (SNR) estimation in simulated industrial, scientific and medical (ISM) band transmission. The source of used background distortions is complex noise measurement, recorded by USRP-2932 in an industrial production area. Based on the examined solutions, we also analyze the influence of noise samples separation techniques on the estimation process. As a response to the defects in the used methods, we propose a novel noise samples separation algorithm based on the adaptation of rank order filtering (ROF). In addition to simple implementation, the proposed method has a very good 0.5 dB root-mean-squared error (RMSE) and smaller than 0.1 dB resolution, thus achieving a performance that is comparable with the methods exploiting information theory concepts.


title: Hole transport across MgO-based magnetic tunnel junctions with high resistance-area product due to oxygen vacancies

abstract: The quantum mechanical tunnelling process conserves the quantum properties of the particle considered. As applied to solid-state tunnelling (SST), this physical law was verified, within the field of spintronics, regarding the electron spin in early experiments across Ge tunnel barriers, and in the 90s across Al2O3 barriers. The conservation of the quantum parameter of orbital occupancy, as grouped into electronic symmetries, was observed in the '00s across MgO barriers, followed by SrTiO3 (STO). Barrier defects, such as oxygen vacancies, partly conserve this electronic symmetry. In the solid-state, an additional subtlety is the sign of the charge carrier: are holes or electrons involved in transport? We demonstrate that SST across MgO magnetic tunnel junctions (MTJs) with a large resistance-area (RA) product involves holes by examining how shifting the MTJ's Fermi level alters the ensuing barrier heights defined by the barrier's oxygen vacancies. In the process, we consolidate the description of tunnel barrier heights induced by specific oxygen-vacancy induced localized states. Our work opens prospects to understand the concurrent observation of high TMR and spin transfer torque across MgO-based nanopillars.


title: Simulating optical coherence tomography for observing nerve activity: a finite difference time domain bi-dimensional model

abstract: We present a finite difference time domain (FDTD) model for computation of A line scans in time domain optical coherence tomography (OCT). By simulating only the end of the two arms of the interferometer and computing the interference signal in post processing, it is possible to reduce the computational time required by the simulations and, thus, to simulate much bigger environments. Moreover, it is possible to simulate successive A lines and thus obtaining a cross section of the sample considered. In this paper we present the model applied to two different samples: a glass rod filled with water-sucrose solution at different concentrations and a peripheral nerve. This work demonstrates the feasibility of using OCT for non-invasive, direct optical monitoring of peripheral nerve activity, which is a long-sought goal of neuroscience.


title: Relating the wave-function collapse with Euler's formula, with applications to Classical Statistical Field Theory

abstract: One attractive interpretation of quantum mechanics is the ensemble interpretation, where Quantum Mechanics merely describes a statistical ensemble of objects and not individual objects. But this interpretation does not address why the wave-function plays a central role in the calculations of probabilities, unlike most other interpretations of quantum mechanics. On the other hand, Classical Statistical Field Theory suffers from severe mathematical inconsistencies (specially for Hamiltonians which are non-polynomial in the fields, e.g. General relativistic statistical field theory). We claim that both problems are related to each other and we propose a solution to both. We prove: 1) the wave-function is a parametrization of any probability distribution of a statistical ensemble: there is a surjective map from an hypersphere to the set of all probability distributions; 2) for a quantum system defined in a 2-dimensional real Hilbert space, the role of the (2-dimensional real) wave-function is identical to the role of the Euler's formula in engineering, while the collapse of the wave-function is identical to selecting the real part of a complex number; 3) the collapse of the wave-function of any quantum system is a recursion of collapses of 2-dimensional real wave-functions; 4) the wave-function parametrization is key in the mathematical definition of Classical Statistical Field Theory we propose here. The same formalism is applied to Quantum Yang-Mills theory and Quantum Gravity in another article.


title: Spatial Joint Species Distribution Modeling using Dirichlet Processes

abstract: Species distribution models usually attempt to explain presence-absence or abundance of a species at a site in terms of the environmental features (socalled abiotic features) present at the site. Historically, such models have considered species individually. However, it is well-established that species interact to influence presence-absence and abundance (envisioned as biotic factors). As a result, there has been substantial recent interest in joint species distribution models with various types of response, e.g., presence-absence, continuous and ordinal data. Such models incorporate dependence between species response as a surrogate for interaction. The challenge we focus on here is how to address such modeling in the context of a large number of species (e.g., order 102) across sites numbering in the order of 102 or 103 when, in practice, only a few species are found at any observed site. Again, there is some recent literature to address this; we adopt a dimension reduction approach. The novel wrinkle we add here is spatial dependence. That is, we have a collection of sites over a relatively small spatial region so it is anticipated that species distribution at a given site would be similar to that at a nearby site. Specifically, we handle dimension reduction through Dirichlet processes joined with spatial dependence through Gaussian processes. We use both simulated data and a plant communities dataset for the Cape Floristic Region (CFR) of South Africa to demonstrate our approach. The latter consists of presence-absence measurements for 639 tree species on 662 locations. Through both data examples we are able to demonstrate improved predictive performance using the foregoing specification.


title: Regularity of the spectrum for expanding maps

abstract: In this short note, we propose to extend differentiability (with respect to a multidimensional parameter) of a normalized eigenfunction associated to the simple, dominating eigenvalue of the weighted transfer operator for a uniformly expanding map, to the whole discrete spectrum. We do so by studying directly the regularity (with respect to the parameter) of the resolvent operator, applying a general regularity result for fixed points of maps having loss of regularity.


title: Recent progress in applying lattice QCD to kaon physics

abstract: Standard lattice calculations in kaon physics are based on the evaluation of matrix elements of local operators between two single-hadron states or a single-hadron state and the vacuum. Recent progress in lattice QCD has gone beyond these standard observables. I will review the status and prospects of lattice kaon physics with an emphasis on non-leptonic $K\to\pi\pi$ decay and long-distance processes including $K^0$-$\overline{K^0}$ mixing and rare kaon decays.


title: Assessing the Usability of a Novel System for Programming Education

abstract: The authors present the results of a simple usability test performed on line_explorer, an innovative tool aimed at letting students explore programming. The system offers an interactive environment where students can learn, review, and practice programming independently or through step-by-step instruction. Students in Information Technology, Computer Science, and Information Systems were surveyed. The findings show that students have interest in this tool, whereas some groups find this tool more interesting and useful. The findings will help refine the user interface for the next phase of testing which include changes for simplicity, usability and expanded topic content. Overall the survey on line_explorer in its current design phase seem more useful for IT and CS majors, however significant changes are still needed.


title: A Tractable Product Channel Model for Line-of-Sight Scenarios

abstract: We present a general and tractable fading model for line-of-sight (LOS) scenarios, which is based on the product of two independent and non-identically distributed $\kappa$-$\mu$ shadowed random variables. Simple closed-form expressions for the probability density function, cumulative distribution function and moment-generating function are derived, which are as tractable as the corresponding expressions derived from a product of Nakagami-$m$ random variables. This model simplifies the challenging characterization of LOS product channels, as well as combinations of LOS channels with non-LOS ones. We leverage these results to analyze performance measures of interest in the contexts of wireless powered and backscatter communications, where both forward and reverse links are inherently of LOS nature, as well as in device-to-device communications subject to composite fading. In these contexts, the model shows a higher flexibility when fitting field measurements with respect to conventional approaches based on product distributions with deterministic LOS, together with a more complete physical interpretation of the underlying propagation characteristics.


title: Section sigma models coupled to symplectic duality bundles on Lorentzian four-manifolds

abstract: We give the global mathematical formulation of a class of generalized four-dimensional theories of gravity coupled to scalar matter and to Abelian gauge fields. In such theories, the scalar fields are described by a section of a surjective pseudo-Riemannian submersion $\pi$ over space-time, whose total space carries a Lorentzian metric making the fibers into totally-geodesic connected Riemannian submanifolds. In particular, $\pi$ is a fiber bundle endowed with a complete Ehresmann connection whose transport acts through isometries between the fibers. In turn, the Abelian gauge fields are "twisted" by a flat symplectic vector bundle defined over the total space of $\pi$. This vector bundle is endowed with a vertical taming which locally encodes the gauge couplings and theta angles of the theory and gives rise to the notion of twisted self-duality, of crucial importance to construct the theory. When the Ehresmann connection of $\pi$ is integrable, we show that our theories are locally equivalent to ordinary Einstein-Scalar-Maxwell theories and hence provide a global non-trivial extension of the universal bosonic sector of four-dimensional supergravity. In this case, we show using a special trivializing atlas of $\pi$ that global solutions of such models can be interpreted as classical "locally-geometric" U-folds. In the non-integrable case, our theories differ locally from ordinary Einstein-Scalar-Maxwell theories and may provide a geometric description of classical U-folds which are "locally non-geometric".


title: On the unit sphere of positive operators

abstract: Given a C$^*$-algebra $A$, let $S(A^+)$ denote the set of those positive elements in the unit sphere of $A$. Let $H_1$, $H_2,$ $H_3$ and $H_4$ be complex Hilbert spaces, where $H_3$ and $H_4$ are infinite-dimensional and separable. In this note we prove a variant of Tingley's problem by showing that every surjective isometry $\Delta : S(B(H_1)^+)\to S(B(H_2)^+)$ or (respectively, $\Delta : S(K(H_3)^+)\to S(K(H_4)^+)$) admits a unique extension to a surjective complex linear isometry from $B(H_1)$ onto $B(H_2))$ (respectively, from $K(H_3)$ onto $B(H_4)$). This provides a positive answer to a conjecture posed by G. Nagy [\emph{Publ. Math. Debrecen}, 2018].


title: Theoretical model of a finite force at the moving contact line

abstract: In theoretical analyses of the moving contact line, an infinite force along the solid wall has been reported based off the non-integrable stress along a single interface. In this investigation we demonstrate that the stress singularity is integrable and results in a finite force at the moving contact line if the contact line is treated as a one-dimensional manifold and all three interfaces that make up the moving contact line are taken into consideration. This is due to the dipole nature of the vorticity and pressure distribution around the moving contact line. Mathematically, this finite force is determined by summing all the forces that act over an infinitesimally small cylindrical control volume that encloses the entire moving contact line. With this finite force, we propose a new dynamic Young's equation for microscopic dynamic contact angle that is a function of known parameters only, specifically the interface velocity, surface tension, and fluid viscosity. We combine our model with Cox's model for apparent dynamic contact angle and find good agreement with published dynamic contact angle measurements.


title: Spectral stability of bi-frequency solitary waves in Soler and Dirac--Klein--Gordon models

abstract: We construct bi-frequency solitary waves of the nonlinear Dirac equation with the scalar self-interaction (the Soler model) and the Dirac--Klein--Gordon with Yukawa self-interaction. These solitary waves provide a natural implementation of qubit and qudit states in the theory of quantum computing. We show the relation of $\pm 2\omega\mathrm{i}$ eigenvalues of the linearization at a solitary wave, Bogoliubov $\mathbf{SU}(1,1)$ symmetry, and the existence of bi-frequency solitary waves. We show that the spectral stability of these waves reduces to spectral stability of usual (one-frequency) solitary waves.


title: A Stochastic Hybrid Framework for Driver Behavior Modeling Based on Hierarchical Dirichlet Process

abstract: Scalability is one of the major issues for real-world Vehicle-to-Vehicle network realization. To tackle this challenge, a stochastic hybrid modeling framework based on a non-parametric Bayesian inference method, i.e., hierarchical Dirichlet process (HDP), is investigated in this paper. This framework is able to jointly model driver/vehicle behavior through forecasting the vehicle dynamical time-series. This modeling framework could be merged with the notion of model-based information networking, which is recently proposed in the vehicular literature, to overcome the scalability challenges in dense vehicular networks via broadcasting the behavioral models instead of raw information dissemination. This modeling approach has been applied on several scenarios from the realistic Safety Pilot Model Deployment (SPMD) driving data set and the results show a higher performance of this model in comparison with the zero-hold method as the baseline.


title: Learning to Predict with Highly Granular Temporal Data: Estimating individual behavioral profiles with smart meter data

abstract: Big spatio-temporal datasets, available through both open and administrative data sources, offer significant potential for social science research. The magnitude of the data allows for increased resolution and analysis at individual level. While there are recent advances in forecasting techniques for highly granular temporal data, little attention is given to segmenting the time series and finding homogeneous patterns. In this paper, it is proposed to estimate behavioral profiles of individuals' activities over time using Gaussian Process-based models. In particular, the aim is to investigate how individuals or groups may be clustered according to the model parameters. Such a Bayesian non-parametric method is then tested by looking at the predictability of the segments using a combination of models to fit different parts of the temporal profiles. Model validity is then tested on a set of holdout data. The dataset consists of half hourly energy consumption records from smart meters from more than 100,000 households in the UK and covers the period from 2015 to 2016. The methodological approach developed in the paper may be easily applied to datasets of similar structure and granularity, for example social media data, and may lead to improved accuracy in the prediction of social dynamics and behavior.


title: Approximation of $m-$subharmonic functions on bounded domains in $\mathbb C^n$

abstract: Let $D$ be a bounded domain in $\mathbb C^n$. We study approximation of (not necessarily bounded from above) $m-$subharmonic function $D$ by continuous $m-$subharmonic ones defined on neighborhoods of $\overline{D}$. We also consider the existence of a $m-$subharmonic function on $D$ whose boundary values coincides with a given real valued continuous function on $\partial D$ except for a sufficiently small subset of $\partial D.$


title: The QCD chiral phase transition from non-integer numbers of flavors

abstract: Attempts to extract the order of the chiral transition of QCD at zero chemical potential, with two dynamical flavors of massless quarks, from simulations with progressively decreasing pion mass have remained inconclusive because of their increasing numerical cost. In an alternative approach to this problem, we consider the path integral as a function of continuous number $N_\text{f}$ of degenerate quarks. If the transition in the chiral limit is first-order for $N_\text{f} \ge 3$, a second-order transition for $N_\text{f}=2$ then requires a tricritical point in between. This in turn implies tricritical scaling of the critical boundary line between the first-order and crossover regions as the chiral limit is approached. Non-integer numbers of fermion flavors are easily implemented within the staggered fermion discretization. Exploratory simulations at $\mu=0$ and $N_\text{f} = 2.8, 2.6, 2.4, 2.2, 2.1$, on coarse $N_\tau = 4$ lattices, indeed show a smooth variation of the critical mass mapping out a critical line in the $(m,N_\text{f})$-plane. For the smallest masses the line appears consistent with tricritical scaling, allowing for an extrapolation to the chiral limit.


title: An inverse problem for Sturm-Liouville operators on trees with partial information given on the potentials

abstract: We consider Sturm-Liouville operators on geometrical graphs without cycles (trees) with singular potentials from the class $W_2^{-1}$. We suppose that the potentials are known on a part of the graph, and study the so-called partial inverse problem, which consists in recovering the potentials on the remaining part of the graph from some parts of several spectra. The main results of the paper are the uniqueness theorem and a constructive procedure for the solution of the partial inverse problem. Our method is based on the completeness and the Riesz-basis property of special systems of vector functions, and the reduction of the partial inverse problem to the complete one on a part of the graph.


title: A partial inverse problem for the Sturm-Liouville operator on the graph with a loop

abstract: The Sturm-Liouville operator with singular potentials on the lasso graph is considered. We suppose that the potential is known a priori on the boundary edge, and recover the potential on the loop from a part of the spectrum and some additional data. We prove the uniqueness theorem and provide a constructive algorithm for the solution of this partial inverse problem.


title: On the $D^*_s$ and charmonia leptonic decays

abstract: Among the different scenarios of New Physics, those with an extended Higgs sector are examined with a lot of attention. Recent experimental observations of several anomalies in flavour physics with respect to expectations of the Standard Model further motivate the effort of phenomenologists. First, informations about the $R_{D_s}$ ratio, a test of lepton flavour universality equivalent to $R_D$, already measured, but with the $s$ quark as spectator, are awaited in coming years to constrain the corner of an extended Higgs sector with charged doublets. On another side, leptonic widths of pseudoscalar quarkonia are particularly interesting to test an extended Higgs sector with a light CP-odd Higgs boson singlet, through the study of its mixing with quarkonia states. Hadronic parameters entering those processes have to be determined from lattice QCD with enough confidence on the control of systematic errors. We report on the very first step of a long-term program tackled with ${\rm N_f}=2$ Wilson-Clover fermions to put relevant constraints on extensions of the Higgs sector: extraction of decay constants of $D^*_s$, $\eta_c$, $\eta_c(2S)$, $J/\psi$ and $\psi(2S)$ with lattice ensembles provided by the CLS effort, considering 2 lattice spacings and a large range of pion masses to estimate cut-off effects and extrapolate results to the chiral limit.


title: Approached vectorial model for Fano resonances in guided mode resonance gratings

abstract: We propose a self-consistent vectorial method, based on a Green's function technique, to describe the Fano resonances that appear in guided mode resonance gratings. The model provides intuitive expressions of the reflectivity and transmittivity matrices of the structure, involving coupling integrals between the modes of a planar reference structure and radiative modes. These expressions are used to derive a physical analysis in configurations where the effect of the incident polarization is not trivial. We provide numerical validations of our model. On a technical point of view, we show how the Green's tensor of our planar reference structure can be expressed as two scalar Green's functions, and how to deal with the singularity of the Green's tensor.


title: Semi-Supervised Approaches to Efficient Evaluation of Model Prediction Performance

abstract: In many modern machine learning applications, the outcome is expensive or time-consuming to collect while the predictor information is easy to obtain. Semi-supervised learning (SSL) aims at utilizing large amounts of `unlabeled' data along with small amounts of `labeled' data to improve the efficiency of a classical supervised approach. Though numerous SSL classification and prediction procedures have been proposed in recent years, no methods currently exist to evaluate the prediction performance of a working regression model. In the context of developing phenotyping algorithms derived from electronic medical records (EMR), we present an efficient two-step estimation procedure for evaluating a binary classifier based on various prediction performance measures in the semi-supervised (SS) setting. In step I, the labeled data is used to obtain a non-parametrically calibrated estimate of the conditional risk function. In step II, SS estimates of the prediction accuracy parameters are constructed based on the estimated conditional risk function and the unlabeled data. We demonstrate that under mild regularity conditions the proposed estimators are consistent and asymptotically normal. Importantly, the asymptotic variance of the SS estimators is always smaller than that of the supervised counterparts under correct model specification. We also correct for potential overfitting bias in the SS estimators in finite sample with cross-validation and develop a perturbation resampling procedure to approximate their distributions. Our proposals are evaluated through extensive simulation studies and illustrated with two real EMR studies aiming to develop phenotyping algorithms for rheumatoid arthritis and multiple sclerosis.


title: Stationary Inviscid Limit to Shear Flows

abstract: In this note we establish a density result for certain stationary shear flows, $\mu(y)$, that vanish at the boundaries of a horizontal channel. We construct stationary solutions to 2D Navier-Stokes that are $\epsilon$-close in $L^\infty$ to the given shear flow. Our construction is based on a coercivity estimate for the Rayleigh operator, $R[v]$, which is based on a decomposition made possible by the vanishing of $\mu$ at the boundaries.


title: A characterization of Fuchsian actions by topological rigidity

abstract: We prove that any rigid representation of $\pi_1\Sigma_g$ in $\mathrm{Homeo}_+(S^1)$ with Euler number at least $g$ is necessarily semi-conjugate to a discrete, faithful representation into $\mathrm{PSL}(2,\mathbb{R})$. Combined with earlier work of Matsumoto, this precisely characterizes Fuchsian actions by a topological rigidity property. Though independent, this work can be read as an introduction to the companion paper {\em rigidity and geometricity for surface group actions on the circle}, by the same authors.


title: Quantum symmetries of the deformation quantization of $SU(3)$

abstract: We prove a criterion of when a coaction of a compact Lie group on an algebra of continuous functions on a compact manifold extends to a coaction of deformation quantizations of the Lie group and the algebra. We compute an explicit example of a compact quantum group $SU(3)_\theta$, which arises as a deformation quantization of the Lie group $SU(3)$ by an action of its maximal torus. Using the criterion, we determine exactly when the action of $SU(3)$ on $S^5$ extends to a coaction of $SU(3)_\theta$ on the noncommutative 5-sphere $S^5_{\theta'}$. Furthermore, this coaction is shown to be cotransitive. A coaction of $SU(3)_\theta$ on the product of two noncommutative 5-sphere $(S^5\times S^5)_{\theta'}$ with nontrivial algebra of coinvariant elements is given and associated projective modules are constructed.


title: A Friendly Smoothed Analysis of the Simplex Method

abstract: Explaining the excellent practical performance of the simplex method for linear programming has been a major topic of research for over 50 years. One of the most successful frameworks for understanding the simplex method was given by Spielman and Teng (JACM `04), who developed the notion of smoothed analysis. Starting from an arbitrary linear program with $d$ variables and $n$ constraints, Spielman and Teng analyzed the expected runtime over random perturbations of the LP (smoothed LP), where variance $\sigma^2$ Gaussian noise is added to the LP data. In particular, they gave a two-stage shadow vertex simplex algorithm which uses an expected $\widetilde{O}(d^{55} n^{86} \sigma^{-30})$ number of simplex pivots to solve the smoothed LP. Their analysis and runtime was substantially improved by Deshpande and Spielman (FOCS `05) and later Vershynin (SICOMP `09). The fastest current algorithm, due to Vershynin, solves the smoothed LP using an expected $O(d^3 \sigma^{-4} \log^3 n + d^9\log^7 n)$ number of pivots, improving the dependence on $n$ from polynomial to poly-logarithmic. While the original proof of Spielman and Teng has now been substantially simplified, the resulting analyses are still quite long and complex and the parameter dependencies far from optimal. In this work, we make substantial progress on this front, providing an improved and simpler analysis of shadow simplex methods, where our algorithm requires an expected \[ O(d^2 \sqrt{\log n} \sigma^{-2} + d^3 \log^{3/2} n) \] number of simplex pivots. We obtain our results via an improved \emph{shadow bound}, key to earlier analyses as well, combined with improvements on algorithmic techniques of Vershynin. As an added bonus, our analysis is completely modular, allowing us to obtain non-trivial bounds for perturbations beyond Gaussians, such as Laplace perturbations.


title: Accelerated horizons and Planck-scale kinematics

abstract: We extend the concept of accelerated horizons to the framework of deformed relativistic kinematics at the Planck scale. We show that the non-trivial effects due to symmetry deformation manifest in a finite blueshift for field modes as measured by a Rindler observer approaching the horizon. We investigate whether, at a field theoretic level, this effect could manifest in the possibility of a finite horizon contribution to the entropy, a sort of covariant brick-wall. In the specific model of symmetry deformation considered it will turn out that a non-diverging density of modes close to the horizon can be achieved only by introducing a momentum space measure which violates Lorentz invariance.


title: Analytic continuation of dimensions in supersymmetric localization

abstract: We compute the perturbative partition functions for gauge theories with eight supersymmetries on spheres of dimension $d\le5$, proving a conjecture by the second author. We apply similar methods to gauge theories with four supersymmetries on spheres with $d\le3$. The results are valid for non-integer $d$ as well. We further propose an analytic continuation from $d=3$ to $d=4$ that gives the perturbative partition function for an $\mathcal{N}=1$ gauge theory. The results are consistent with the free multiplets and the one-loop $\beta$-functions for general $\mathcal{N}=1$ gauge theories. We also consider the analytic continuation of an $\mathcal{N}=1$-preserving mass deformation of the maximally supersymmetric gauge theory and compare to recent holographic results for $\mathcal{N}=1^*$ super Yang-Mills. We find that the general structure for the real part of the free energy coming from the analytic continuation is consistent with the holographic results.


title: Anti-reflection coating design for metallic terahertz meta-materials

abstract: We demonstrate a silicon-based, single-layer anti-reflection coating that suppresses the reflectivity of metals at near-infrared frequencies, enabling optical probing of nano-scale structures embedded in highly reflective surroundings. Our design does not affect the interaction of terahertz radiation with metallic structures that can be used to achieve terahertz near-field enhancement. We have verified the functionality of the design by calculating and measuring the reflectivity of both infrared and terahertz radiation from a silicon/gold double layer as a function of the silicon thickness. We have also fabricated the unit cell of a terahertz meta-material, a dipole antenna comprising two 20-nm thick extended gold plates separated by a 2 $\mu$m gap, where the terahertz field is locally enhanced. We used the time-domain finite element method to demonstrate that such near-field enhancement is preserved in the presence of the anti-reflection coating. Finally, we performed magneto-optical Kerr effect measurements on a single 3-nm thick, 1-$\mu$m wide magnetic wire placed in the gap of such a dipole antenna. The wire only occupies 2\% of the area probed by the laser beam, but its magneto-optical response can be clearly detected. Our design paves the way for ultrafast time-resolved studies, using table-top femtosecond near-infrared lasers, of dynamics in nano-structures driven by strong terahertz radiation.


title: A spectral Szego theorem on the real line

abstract: We characterize even measures $\mu=wdx+\mu_s$ on the real line with finite entropy integral $\int_{R} \frac{\log w(t)}{1+t^2}dt>-\infty$ in terms of $2\times 2$ Hamiltonian generated by $\mu$ in the sense of inverse spectral theory. As a corollary, we obtain criterion for spectral measure of Krein string to have converging logarithmic integral.


title: Monotonicity of maximal equicontinuous factors and an application to toral flows

abstract: We show that for group actions on locally connected spaces the maximal equicontinuous factor map is always monotone, that is, the preimages of single points are connected. As an application, we obtain that if the maximal continuous factor of a homeomorphism of the two-torus is minimal, then it is either (i) an irrational translation of the two-torus, (ii) an irrational rotation on the circle or (iii) the identity on a singleton.


title: Counting factorisations of monomials over rings of integers modulo $N$

abstract: A sharp bound is obtained for the number of ways to express the monomial $X^n$ as a product of linear factors over $\mathbb{Z}/p^{\alpha}\mathbb{Z}$. The proof relies on an induction-on-scale procedure which is used to estimate the number of solutions to a certain system of polynomial congruences. The method also applies to more general systems of polynomial congruences that satisfy a non-degeneracy hypothesis.


title: On laws of large numbers in $L^2$ for supercritical branching Markov processes beyond $\lambda$-positivity

abstract: We give necessary and sufficient conditions for laws of large numbers to hold in $L^2$ for the empirical measure of a large class of branching Markov processes, including $\lambda$-positive systems but also some $\lambda$-transient ones, such as the branching Brownian motion with drift and absorption at $0$. This is a significant improvement over previous results on this matter, which had only dealt so far with $\lambda$-positive systems. Our approach is purely probabilistic and is based on spinal decompositions and many-to-few lemmas. In addition, we characterize when the limit in question is always strictly positive on the event of survival, and use this characterization to derive a simple method for simulating (quasi-)stationary distributions.


title: On joint distribution of range and terminal value of a Brownian motion

abstract: In this note, we present the closed form solution for the joint distribution of the range and terminal value of a Brownian motion. Based on this distribution we build a range scaled terminal value distribution and show the derivation steps of its density, further s-density. Finally, we sample the s-density from different groups of currency pairs and compare them with theoretical result.


title: Classification of kinematical Lie algebras

abstract: We summarise the classification of kinematical Lie algebras in arbitrary dimension and indicate which of the kinematical Lie algebras admit an invariant inner product.


title: PQSER: A Matlab package for spectral seriation

abstract: The seriation problem is an important ordering issue which consists of finding the best ordering of a set of units whose interrelationship is defined by a bipartite graph. It has important applications in, e.g., archaeology, anthropology, psychology, and biology. This paper presents a Matlab implementation of an algorithm for spectral seriation by Atkins et al., based on the use of the Fiedler vector of the Laplacian matrix associated to the problem, which encodes the set of admissible solutions into a PQ-tree. We introduce some numerical technicalities in the original algorithm to improve its performance, and point out that the presence of a multiple Fiedler value may have a substantial influence on the computation of an approximated solution, in the presence of inconsistent data sets. Practical examples and numerical experiments show how to use the toolbox to process data sets deriving from real-world applications.


title: Unsupervised Morphological Expansion of Small Datasets for Improving Word Embeddings

abstract: We present a language independent, unsupervised method for building word embeddings using morphological expansion of text. Our model handles the problem of data sparsity and yields improved word embeddings by relying on training word embeddings on artificially generated sentences. We evaluate our method using small sized training sets on eleven test sets for the word similarity task across seven languages. Further, for English, we evaluated the impacts of our approach using a large training set on three standard test sets. Our method improved results across all languages.


title: New measurements on water ice photodesorption and product formation under ultraviolet irradiation

abstract: The photodesorption of icy grain mantles has been claimed to be responsible for the abundance of gas-phase molecules toward cold regions. Being water a ubiquitous molecule, it is crucial to understand its role in photochemistry and its behavior under an ultraviolet field. We report new measurements on the UV-photodesorption of water ice and its H$_2$, OH, and O$_2$ photoproducts using a calibrated quadrupole mass spectrometer. Solid water was deposited under ultra-high-vacuum conditions and then UV-irradiated at various temperatures starting from 8 K with a microwave discharged hydrogen lamp. Deuterated water was used for confirmation of the results. We found a photodesorption yield of 1.3 $\times$ 10$^{-3}$ molecules per incident photon for water, and 0.7 $\times$ 10$^{-3}$ molecules per incident photon for deuterated water at the lowest irradiation temperature, 8 K. The photodesorption yield per absorbed photon is given and comparison with astrophysical scenarios, where water ice photodesorption could account for the presence of gas-phase water toward cold regions in the absence of a thermal desorption process is addressed.


title: An Unsupervised Approach for Mapping between Vector Spaces

abstract: We present a language independent, unsupervised approach for transforming word embeddings from source language to target language using a transformation matrix. Our model handles the problem of data scarcity which is faced by many languages in the world and yields improved word embeddings for words in the target language by relying on transformed embeddings of words of the source language. We initially evaluate our approach via word similarity tasks on a similar language pair - Hindi as source and Urdu as the target language, while we also evaluate our method on French and German as target languages and English as source language. Our approach improves the current state of the art results - by 13% for French and 19% for German. For Urdu, we saw an increment of 16% over our initial baseline score. We further explore the prospects of our approach by applying it on multiple models of the same language and transferring words between the two models, thus solving the problem of missing words in a model. We evaluate this on word similarity and word analogy tasks.


title: Forecasting dynamic return distributions based on ordered binary choice

abstract: We present a simple approach to forecasting conditional probability distributions of asset returns. We work with a parsimonious specification of ordered binary choice regression that imposes a connection on sign predictability across different quantiles. The model forecasts the future conditional probability distributions of returns quite precisely when using a past indicator and past volatility proxy as predictors. Direct benefits of the model are revealed in an empirical application to the 29 most liquid U.S. stocks. The forecast probability distribution is translated to significant economic gains in a simple trading strategy. Our approach can also be useful in many other applications where conditional distribution forecasts are desired.


title: The Gauss map of a free boundary minimal surface

abstract: In this paper, we study the Gauss map of a free boundary minimal surface. The main theorem asserts that if components of the Gauss map are eigenfunctions of the Jacobi-Steklov operator, then the surface must be rotationally symmetric.


title: Hydra: a C++11 framework for data analysis in massively parallel platforms

abstract: Hydra is a header-only, templated and C++11-compliant framework designed to perform the typical bottleneck calculations found in common HEP data analyses on massively parallel platforms. The framework is implemented on top of the C++11 Standard Library and a variadic version of the Thrust library and is designed to run on Linux systems, using OpenMP, CUDA and TBB enabled devices. This contribution summarizes the main features of Hydra. A basic description of the overall design, functionality and user interface is provided, along with some code examples and measurements of performance.

categories: cs.MS hep-ex physics.comp-ph physics.data-an


title: Reversible viscosity and Navier--Stokes fluids

abstract: Exploring the possibility of describing a fluid flow via a time-reversible equation and its relevance for the fluctuations statistics in stationary turbulent (or laminar) incompressible Navier-Stokes flows.


title: Hamiltonian Renormalisation I: Derivation from Osterwalder-Schrader Reconstruction

abstract: A possible avenue towards a non-perturbative Quantum Field Theory (QFT) on Minkowski space is the constructive approach which employs the Euclidian path integral formulation, in the presence of both ultraviolet (UV) and infrared (IR) regulators, as starting point. The UV regulator is to be taken away by renormalisation group techniques which in case of success leads to a measure on the space of generalised Euclidian fields in finite volume. The IR regulator corresponds to the thermodynamic limit of the system in the statistical physics sense. If the resulting measure obeys the Osterwalder-Schrader axioms, the actual QFT on Minkowski space is then obtained by Osterwalder-Schrader reconstruction. In this work we study the question whether it is possible to reformulate the renormalisation group non-perturbatively directly at the operator (Hamiltonian) level. Hamiltonian renormalisation would be the natural route to follow if one had easier access to an interacting Hamiltonian operator rather than to a path integral, at least in the presence of UV and/or IR cut-off, which is generically the case in complicated gauge theories such as General Relativity. Our guiding principle for the definition of the direct Hamiltonian renormalisation group is that it results in the same continuum theory as the covariant (path integral) renormalisation group. In order to achieve this, we invert the Osterwalder-Schrader reconstruction, which may be called Osterwalder-Schrader construction of a Wiener measure from the underlying Hamiltonian theory. The resulting correspondence between reflection positive measures and Osterwalder-Schrader data consisting of a Hilbert space, a Hamiltonian and a ground state vector allows us to monitor the effect of the renormalisation flow of measures in terms of their Osterwalder-Schrader data which motivates a natural direct Hamiltonian renormalisation scheme.

categories: gr-qc hep-lat hep-th math-ph math.MP quant-ph


title: Novel decision-theoretic and risk-stratification metrics of predictive performance: Application to deciding who should undergo genetic testing

abstract: Currently, women are referred for BRCA1/2 mutation-testing only if their family-history of breast/ovarian cancer implies that their risk of carrying a mutation exceeds 10\%. However, as mutation-testing costs fall, prominent voices have called for testing all women, which would strain clinical resources by testing millions of women, almost all of whom will test negative. To better evaluate risk-thresholds for BRCA1/2 testing, we introduce two broadly applicable, linked metrics: Mean Risk Stratification (MRS) and a decision-theoretic metric, Net Benefit of Information (NBI). MRS and NBI provide a range of risk thresholds at which a marker/model is "optimally informative", in the sense of maximizing both MRS and NBI. NBI is a function of only MRS and the risk-threshold for action, connecting decision-theory to risk-stratification and providing a decision-theoretic rationale for MRS. AUC and Youden's index reflect on both the fraction of maximum MRS, and of maximum NBI, attained by the marker/model, providing AUC and Youden's index with long-sought decision-theoretic and risk-stratification rationale. To evaluate risk-thresholds for BRCA1/2 testing, we propose an eclectic approach considering AUC, Net Benefit, and MRS/NBI. MRS/NBI interpret AUC in the context of mutation-prevalence and provide a range of risk thresholds for which the risk model is optimally informative.


title: Interstellar Interloper 1I/2017 U1: Observations from the NOT and WIYN Telescopes

abstract: We present observations of the interstellar interloper 1I/2017 U1 ('Oumuamua) taken during its 2017 October flyby of Earth. The optical colors B-V = 0.70$\pm$0.06, V-R = 0.45$\pm$0.05, overlap those of the D-type Jovian Trojan asteroids and are incompatible with the ultrared objects which are abundant in the Kuiper belt. With a mean absolute magnitude $H_V$ = 22.95 and assuming a geometric albedo $p_V$ = 0.1, we find an average radius of 55 m. No coma is apparent; we deduce a limit to the dust mass production rate of only $\sim$ 2$\times$10$^{-4}$ kg s$^{-1}$, ruling out the existence of exposed ice covering more than a few m$^2$ of the surface. Volatiles in this body, if they exist, must lie beneath an involatile surface mantle $\gtrsim$0.5 m thick, perhaps a product of prolonged cosmic ray processing in the interstellar medium. The lightcurve range is unusually large at $\sim$2.0$\pm$0.2 magnitudes. Interpreted as a rotational lightcurve the body has semi-axes $\sim$230 m $\times$ 35 m. A $\sim$6:1 axis ratio is extreme relative to most small solar system asteroids and suggests that albedo variations may additionally contribute to the variability. The lightcurve is consistent with a two-peaked period $\sim$8.26 hr but the period is non-unique as a result of aliasing in the data. Except for its unusually elongated shape, 1I/2017 U1 is a physically unremarkable, sub-kilometer, slightly red, rotating object from another planetary system. The steady-state population of similar, $\sim$100 m scale interstellar objects inside the orbit of Neptune is $\sim$10$^4$, each with a residence time $\sim$10 yr.


title: Hamiltonian Renormalization III. Renormalisation Flow of 1+1 dimensional free scalar fields: Properties

abstract: This is the third paper in a series of four in which a renormalisation flow is introduced which acts directly on the Osterwalder-Schrader data (OS data) without recourse to a path integral. Here the OS data consist of a Hilbert space, a cyclic vacuum vector therein and a Hamiltonian annihilating the vacuum which can be obtained from an OS measure, that is a measure respecting (a subset of) the OS axioms. In the previous paper we successfully tested our proposal for the two-dimensional massive Klein-Gordon model, that is, we could confirm that our framework finds the correct fixed point starting from a natural initial naive discretisation of the finite resolution Hamiltonians, in particular the underlying Laplacian on a lattice, and a natural coarse graining map that drives the renormalisation flow. However, several questions remained unanswered. How generic can the initial discretisation and the coarse graining map be in order that the fixed point is not changed or is at least not lost, in other words, how universal is the fixed point structure? Is the fix point in fact stable, that is, is the fixed point actually a limit of the renormalisation sequence? We will address these questions in the present paper.


title: Active particles in periodic lattices

abstract: Both natural and artificial small-scale swimmers may often self-propel in environments subject to complex geometrical constraints. While most past theoretical work on low-Reynolds number locomotion addressed idealised geometrical situations, not much is known on the motion of swimmers in heterogeneous environments. As a first theoretical model, we investigate numerically the behaviour of a single spherical micro-swimmer located in an infinite, periodic body-centred cubic lattice consisting of rigid inert spheres of the same size as the swimmer. Running a large number of simulations we uncover the phase diagram of possible trajectories as a function of the strength of the swimming actuation and the packing density of the lattice. We then use hydrodynamic theory to rationalise our computational results and show in particular how the far-field nature of the swimmer (pusher vs. puller) governs even the behaviour at high volume fractions.


title: Orthogonal foliations on riemannian manifolds

abstract: In this work, we find an equation that relates the Ricci curvature of a riemannian manifold $M$ and the second fundamental forms of two orthogonal foliations of complementary dimensions, $\mathcal{F}$ and $\mathcal{F}^{\bot}$, defined on $M$. Using this equation, we show a sufficient condition for the manifold M to be locally a riemannian product of the leaves of $\mathcal{F}$ and $\mathcal{F}^{\bot}$, if one of the foliations is totally umbilical. We also prove an integral formula for such foliations.


title: The Transit Light Source Effect: False Spectral Features and Incorrect Densities for M-dwarf Transiting Planets

abstract: Transmission spectra are differential measurements that utilize stellar illumination to probe transiting exoplanet atmospheres. Any spectral difference between the illuminating light source and the disk-integrated stellar spectrum due to starspots and faculae will be imprinted in the observed transmission spectrum. However, few constraints exist for the extent of photospheric heterogeneities in M dwarfs. Here, we model spot and faculae covering fractions consistent with observed photometric variabilities for M dwarfs and the associated 0.3-5.5 $\mu$m stellar contamination spectra. We find that large ranges of spot and faculae covering fractions are consistent with observations and corrections assuming a linear relation between variability amplitude and covering fractions generally underestimate the stellar contamination. Using realistic estimates for spot and faculae covering fractions, we find stellar contamination can be more than $10 \times$ larger than transit depth changes expected for atmospheric features in rocky exoplanets. We also find that stellar spectral contamination can lead to systematic errors in radius and therefore the derived density of small planets. In the case of the TRAPPIST-1 system, we show that TRAPPIST-1's rotational variability is consistent with spot covering fractions $f_{spot} = 8^{+18}_{-7}\%$ and faculae covering fractions $f_{fac} = 54^{+16}_{-46}\%$. The associated stellar contamination signals alter transit depths of the TRAPPIST-1 planets at wavelengths of interest for planetary atmospheric species by roughly 1-15 $\times$ the strength of planetary features, significantly complicating $JWST$ follow-up observations of this system. Similarly, we find stellar contamination can lead to underestimates of bulk densities of the TRAPPIST-1 planets of $\Delta(\rho) = -3^{+3}_{-8} \%$, thus leading to overestimates of their volatile contents.


title: Bundling of elastic filaments induced by hydrodynamic interactions

abstract: Peritrichous bacteria swim in viscous fluids by rotating multiple helical flagellar filaments. As the bacterium swims forward, all its flagella rotate in synchrony behind the cell in a helical bundle. When the bacterium changes its direction, the flagellar filaments unbundle and reorient the cell for a short period of time before returning to their bundled state and resuming swimming. This rapid bundling and unbundling is, at its heart, a mechanical process whereby hydrodynamic interactions balance with elasticity to determine the time-varying deformation of the filaments. Inspired by this biophysical problem, we present in this paper what is perhaps the simplest model of bundling whereby two, or more, straight elastic filaments immersed in a viscous fluid rotate about their centreline, inducing rotational flows which tend to bend the filaments around each other. We derive an integro-differential equation governing the shape of the filaments resulting from mechanical balance in a viscous fluid at low Reynolds number. We show that such equation may be evaluated asymptotically analytically in the long-wavelength limit, leading to a local partial differential equation governed by a single dimensionless Bundling number. A numerical study of the dynamics predicted by the model reveals the presence of two configuration instabilities with increasing Bundling numbers: first to a crossing state where filaments touch at one point and then to a bundled state where filaments wrap along each other in a helical fashion. We next provide an intuitive physical model for the crossing instability and show that it may be used to predict analytically its threshold, and adapted to address the transition to a bundling state. We then use a macro-scale experimental implementation of the two-filament configuration in order to validate our theoretical predictions and obtain excellent agreement.


title: Connecting Loop Quantum Gravity and String Theory via Quantum Geometry

abstract: We argue that String Theory and Loop Quantum Gravity can be thought of as describing different regimes of a single unified theory of quantum gravity. LQG can be thought of as providing the pre-geometric exoskeleton out of which macroscopic geometry emerges and String Theory then becomes the \emph{effective} theory which describes the dynamics of that exoskeleton. The core of the argument rests on the claim that the Nambu-Goto action of String Theory can be viewed as the expectation value of the LQG area operator evaluated on the string worldsheet. A concrete result is that the string tension of String Theory and the Barbero-Immirzi parameter of LQG turn out to be proportional to each other.


title: Extended I-Love relations for slowly rotating neutron stars

abstract: Observations of gravitational waves from inspiralling neutron star binaries---such as GW170817---can be used to constrain the nuclear equation of state by placing bounds on stellar tidal deformability. For slowly rotating neutron stars, the response to a weak quadrupolar tidal field is characterized by four internal-structure-dependent constants called "Love numbers." The tidal Love numbers $k_2^\text{el}$ and $k_2^\text{mag}$ measure the tides raised by the gravitoelectric and gravitomagnetic components of the applied field, and the rotational-tidal Love numbers $\mathfrak{f}^\text{o}$ and $\mathfrak{k}^\text{o}$ measure those raised by couplings between the applied field and the neutron star spin. In this work we compute these four Love numbers for perfect fluid neutron stars with realistic equations of state. We discover (nearly) equation-of-state independent relations between the rotational-tidal Love numbers and the moment of inertia, thereby extending the scope of I-Love-Q universality. We find that similar relations hold among the tidal and rotational-tidal Love numbers. These relations extend the applications of I-Love universality in gravitational-wave astronomy. As our findings differ from those reported in the literature, we derive general formulas for the rotational-tidal Love numbers in post-Newtonian theory and confirm numerically that they agree with our general-relativistic computations in the weak-field limit.


title: Hamiltonian Renormalisation IV. Renormalisation Flow of D+1 dimensional free scalar fields and Rotation Invariance

abstract: In this article we extend the test of Hamiltonian Renormalisation proposed in this series of articles to the D-dimensional case using a massive free scalar field. The concepts we introduce are explicitly computed for the D=2 case but transfer immediately to higher dimensions. In this article we define and verify a criterion that monitors, at finite resolution defined by a cubic lattice, whether the flow approaches a rotationally invariant fixed point.


title: An approach to evaluation of common DNS misconfigurations

abstract: DNS is a basic Internet service which almost all other user services depend on. However, what has been perceived in practice are a lot of inconsistencies and errors in the configuration of servers that cause different problems. The majority of such cases are included in this research with the aim of identifying and classifying the major problems of DNS availability, performance and security. In order to analyze these problems in correlation with DNS administrators working practice, we have developed a methodology and tool for testing, quantifying and analysis of DNS misconfigurations. The methodology and tool were applied on three heterogeneous domain categories - the most popular Internet domains, academic domains and one national top level domain. Our results confirm relatively high percentage of misconfigured domains, especially in the academic and national categories. However, we have shown that fixing the configuration on relatively small number of name servers can have significant impact to great number of domains. Proper domain management, permanent testing and collaboration with other administrators are identified as measures to improve domains operation, stability and security.


title: Motif-based Convolutional Neural Network on Graphs

abstract: This paper introduces a generalization of Convolutional Neural Networks (CNNs) to graphs with irregular linkage structures, especially heterogeneous graphs with typed nodes and schemas. We propose a novel spatial convolution operation to model the key properties of local connectivity and translation invariance, using high-order connection patterns or motifs. We develop a novel deep architecture Motif-CNN that employs an attention model to combine the features extracted from multiple patterns, thus effectively capturing high-order structural and feature information. Our experiments on semi-supervised node classification on real-world social networks and multiple representative heterogeneous graph datasets indicate significant gains of 6-21% over existing graph CNNs and other state-of-the-art techniques.


title: Efficient Verification of Multi-Property Designs (The Benefit of Wrong Assumptions) (Extended Version)

abstract: We consider the problem of efficiently checking a set of safety properties P1,....,Pk of one design. We introduce a new approach called JA-verification where JA stands for "Just-Assume" (as opposed to "assume-guarantee"). In this approach, when proving property Pi, one assumes that every property Pj for j!=i holds. The process of proving properties either results in showing that P1,....,Pk hold without any assumptions or finding a "debugging set" of properties. The latter identifies a subset of failed properties that cause failure of other properties. The design behaviors that cause the properties in the debugging set to fail must be fixed first. Importantly, in our approach, there is no need to prove the assumptions used. We describe the theory behind our approach and report experimental results that demonstrate substantial gains in performance, especially in the cases where a small debugging set exists.


title: Heisenberg-Langevin Formalism For Open Circuit-QED Systems

abstract: We present a Heisenberg-Langevin formalism to study the effective dynamics of a superconducting qubit coupled to an open multimode resonator, without resorting to the rotating wave, two level, Born or Markov approximations. Our effective equations are derived by eliminating resonator degrees of freedom while encoding their effect in the Green's function of the electromagnetic background. We account for the openness of the resonator exactly by employing a spectral representation for the Green's function in terms of a set of non-Hermitian modes. A well-behaved time domain perturbation theory is derived to systematically account for the nonlinearity of weakly nonlinear qubits like transmon. We apply this method to the problem of spontaneous emission, capturing accurately the non-Markovian features of the qubit dynamics, valid for any qubit-resonator coupling strength. Any discrete-level quantum system coupled to the electromagnetic continuum is subject to radiative decay and renormalization of its energy levels. When inside a cavity, these quantities can be strongly modified with respect to vacuum. Generally, this modification can be captured by including only the closest resonant cavity mode. In circuit-QED architecture, with substantial coupling strengths, it is however found that such rates are strongly influenced by far off-resonant modes. A multimode calculation over the infinite set of cavity modes leads to divergences unless an artificial cutoff is imposed. Previous studies have not pointed out what the source of this divergence is. We show that unless the effect of $A^2$ is accounted for up to all orders exactly, any multimode calculations of circuit-QED quantities is bound to diverge. Subsequently, we present the calculation of finite radiative corrections to qubit properties that is free of an artificially introduced high frequency cut-off.


title: Subdiffusion via dynamical localization induced by thermal equilibrium fluctuations

abstract: We reveal the mechanism of subdiffusion which emerges in a straightforward, one dimensional classical nonequilibrium dynamics of a Brownian ratchet driven by both a time-periodic force and Gaussian white noise. In a tailored parameter set for which the deterministic counterpart is in a non-chaotic regime, subdiffusion is a long-living transient whose lifetime can be many, many orders of magnitude larger than characteristic time scales of the setup thus being amenable to experimental observations. As a reason for this subdiffusive behaviour in the coordinate space we identify thermal noise induced dynamical localization in the velocity (momentum) space. This novel idea is distinct from existing knowledge and has never been reported for any classical or quantum systems. It suggests reconsideration of generally accepted opinion that subdiffusion is due to road distributions or strong correlations which reflect disorder, trapping, viscoelasticity of the medium or geometrical constraints.


title: Social Complex Contagion in Music Listenership: A Natural Experiment with 1.3 Million Participants

abstract: Can live music events generate complex contagion in music streaming? This paper finds evidence in the affirmative, but only for the most popular artists. We generate a novel dataset from Last.fm, a music tracking website, to analyse the listenership history of 1.3 million users over a two-month time horizon. We use daily play counts along with event attendance data to run a regression discontinuity analysis in order to show the causal impact of concert attendance on music listenership among attendees and their friends network. First, we show that attending a music artist's live concert increases that artist's listenership among the attendees of the concert by approximately 1 song per day per attendee (p-value<0.001). Moreover, we show that this effect is contagious and can spread to users who did not attend the event. However, the extent of contagion depends on the type of artist. We only observe contagious increases in listenership for well-established, popular artists (.06 more daily plays per friend of an attendee [p<0.001]), while the effect is absent for emerging stars. We also show that the contagion effect size increases monotonically with the number of friends who have attended the live event.


title: Brain Extraction from Normal and Pathological Images: A Joint PCA/Image-Reconstruction Approach

abstract: Brain extraction from images is a common pre-processing step. Many approaches exist, but they are frequently only designed to perform brain extraction from images without strong pathologies. Extracting the brain from images with strong pathologies, for example, the presence of a tumor or of a traumatic brain injury, is challenging. In such cases, tissue appearance may deviate from normal tissue and violates algorithmic assumptions for these approaches; hence, the brain may not be correctly extracted. This paper proposes a brain extraction approach which can explicitly account for pathologies by jointly modeling normal tissue and pathologies. Specifically, our model uses a three-part image decomposition: (1) normal tissue appearance is captured by principal component analysis, (2) pathologies are captured via a total variation term, and (3) non-brain tissue is captured by a sparse term. Decomposition and image registration steps are alternated to allow statistical modeling in a fixed atlas space. As a beneficial side effect, the model allows for the identification of potential pathologies and the reconstruction of a quasi-normal image in atlas space. We demonstrate the effectiveness of our method on four datasets: the IBSR and LPBA40 datasets which show normal images, the BRATS dataset containing images with brain tumors and a dataset containing clinical TBI images. We compare the performance with other popular models: ROBEX, BEaST, MASS, BET, BSE and a recently proposed deep learning approach. Our model performs better than these competing methods on all four datasets. Specifically, our model achieves the best median (97.11) and mean (96.88) Dice scores over all datasets. The two best performing competitors, ROBEX and MASS, achieve scores of 96.23/95.62 and 96.67/94.25 respectively. Hence, our approach is an effective method for high quality brain extraction on a wide variety of images.


title: Geometric parametric instability in periodically modulated GRIN multimode fibers

abstract: We present a theoretical and numerical study of light propagation in graded-index (GRIN) multimode fibers where the core diameter has been periodically modulated along the propagation direction. The additional degree of freedom represented by the modulation permits to modify the intrinsic spatiotemporal dynamics which appears in multimode fibers. More precisely, we show that modulating the core diameter at a periodicity close to the self-imaging distance allows to induce a Moir\'{e}-like pattern, which modifies the geometric parametric instability gain observed in homogeneous GRIN fibers.


title: Bayesian optimal designs for dose-response curves with common parameters

abstract: The issue of determining not only an adequate dose but also a dosing frequency of a drug arises frequently in Phase II clinical trials. This results in the comparison of models which have some parameters in common. Planning such studies based on Bayesian optimal designs offers robustness to our conclusions since these designs, unlike locally optimal designs, are efficient even if the parameters are misspecified. In this paper we develop approximate design theory for Bayesian $D$-optimality for nonlinear regression models with common parameters and investigate the cases of common location or common location and scale parameters separately. Analytical characterisations of saturated Bayesian $D$-optimal designs are derived for frequently used dose-response models and the advantages of our results are illustrated via a numerical investigation.


title: Contextual Object Detection with a Few Relevant Neighbors

abstract: A natural way to improve the detection of objects is to consider the contextual constraints imposed by the detection of additional objects in a given scene. In this work, we exploit the spatial relations between objects in order to improve detection capacity, as well as analyze various properties of the contextual object detection problem. To precisely calculate context-based probabilities of objects, we developed a model that examines the interactions between objects in an exact probabilistic setting, in contrast to previous methods that typically utilize approximations based on pairwise interactions. Such a scheme is facilitated by the realistic assumption that the existence of an object in any given location is influenced by only few informative locations in space. Based on this assumption, we suggest a method for identifying these relevant locations and integrating them into a mostly exact calculation of probability based on their raw detector responses. This scheme is shown to improve detection results and provides unique insights about the process of contextual inference for object detection. We show that it is generally difficult to learn that a particular object reduces the probability of another, and that in cases when the context and detector strongly disagree this learning becomes virtually impossible for the purposes of improving the results of an object detector. Finally, we demonstrate improved detection results through use of our approach as applied to the PASCAL VOC and COCO datasets.


title: Thermal transitions, pseudogap behavior and BCS-BEC crossover in Fermi-Fermi mixtures

abstract: We study the mass imbalanced Fermi-Fermi mixture within the framework of a two-dimensional lattice fermion model. Based on the thermodynamic and species dependent quasiparticle behavior we map out the finite temperature phase diagram of this system and show that unlike the balanced Fermi superfluid there are now two different pseudogap regimes as PG-I and PG-II. While within the PG-I regime both the fermionic species are pseudogapped, PG-II corresponds to the regime where pseudogap feature survives only in the light species. We believe that the single particle spectral features that we discuss in this paper are observable through the species resolved radio frequency spectroscopy and momentum resolved photo emission spectroscopy measurements on systems such as, 6$_{Li}$-40$_{K}$ mixture. We further investigate the interplay between the population and mass imbalances and report that at a fixed population imbalance the BCS-BEC crossover in a Fermi-Fermi mixture would require a critical interaction (U$_{c}$), for the realization of the uniform superfluid state. The effect of imbalance in mass on the exotic Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) superfluid phase has been probed in detail in terms of the thermodynamic and quasiparticle behavior of this phase. It has been observed that in spite of the s-wave symmetry of the pairing field a nodal superfluid gap is realized in the LO regime. Our results on the various thermal scales and regimes are expected to serve as benchmarks for the experimental observations on 6$_{Li}$-40$_{K}$ mixture.


title: Emission of primordial bosonic radiation during inflation

abstract: We study the emission of neutral massless $(1, 2)\hbar$-spin bosons during power-law inflation using unified spinor field theory. We shows that during inflation gravitons and photons were emitted with wavelengths (on physical coordinates) that increase as the Hubble radius: $\lambda_{Ph} \sim a/H$. The quantised action related to these bosons is calculated and results to be a fraction of the Planck constant.


title: Steinberg Summands in the free $\mathbb{F}_p$-module on the Equivariant Sphere Spectrum

abstract: Let $G$ be a finite $p$-group. The Eilenberg-Maclane spectrum of the constant Mackey functor $\underline{\mathbb{F}}_p$, denoted $H\underline{\mathbb{F}}_p$, is modeled by the free $\mathbb{F}_p$-module on the $G$-equivariant sphere spectrum. With this construction, one has a `word length' filtration $\{(H\underline{\mathbb{F}}_p)_n\}_{n\ge 1}$. Our main theorem is that the $k$-th layer $(H\underline{\mathbb{F}}_p)_{p^k}/(H\underline{\mathbb{F}}_p)_{p^{k-1}}$ is $p$-locally equivalent to the $k$-fold suspension of the Steinberg summand of the $G$-equivariant classifying space of $(\mathbb{Z}/p)^k$. This is a generalization of the main result of [21]. We also show that when one smashes this filtration with $H\underline{\mathbb{F}}_p$, the filtration splits into its layers. The future goal of this work is to compute the $C_p$-equivariant dual Steenrod algebra $H\underline{\mathbb{F}}_p\wedge H\underline{\mathbb{F}}_p$ when $p>2$, via explicit cellular constructions of equivariant classifying spaces.