title: On stabilizing the variance of dynamic functional brain connectivity time series

abstract: Assessment of dynamic functional brain connectivity (dFC) based on fMRI data is an increasingly popular strategy to investigate temporal dynamics of the brain's large-scale network architecture. Current practice when deriving connectivity estimates over time is to use the Fisher transform which aims to stabilize the variance of correlation values that fluctuate around varying true correlation values. It is however unclear how well the stabilization of signal variance performed by the Fisher transform works for each connectivity time series, when the true correlation is assumed to be fluctuating. This is of importance because many subsequent analyses either assume or perform better when the time series have stable variance or adheres to an approximate Gaussian distribution. In this paper, using simulations and analysis of resting-state fMRI data, we analyze the effect of applying different variance stabilization strategies on connectivity time-series. We here focus our investigation on the Fisher transform, the Box Cox transform and an approach that combines both transforms. Our results show that, if the intention of stabilizing the variance is to use metrics on the time series where stable variance or a Gaussian distribution is desired (e.g. clustering), the Fisher transform is not optimal and may even skew connectivity time series away from being Gaussian. Further, we show that the suboptimal performance of the Fisher transform can be substantially improved by including an additional Box-Cox transformation after the dFC time series has been Fisher transformed.


title: Peculiar Lapse of Periodic Eclipsing Event at Low Mass X-ray Binary GRS 1747$-$312 during Suzaku Observation in 2009

abstract: GRS 1747$-$312 is a neutron star Low-Mass X-ray Binary in the globular cluster Terzan 6, located at a distance of 9.5 kpc from the Earth. During its outbursts, periodic eclipses were known to occur. Observations for the outbursts were performed with Chandra in 2004 and Swift in 2013. XMM-Newton observed its quiescent state in 2004. In addition, when Suzaku observed it in 2009 as a part of Galactic center mapping observations, GRS 1747$-$312 was found to be in a low luminosity state with $L_{\rm x} \sim 1.2 \times 10^{35}$ erg s$^{-1}$. All of the observations except for XMM-Newton included the time of the eclipses predicted. We analyzed archival data of these observations. During the Chandra and Swift observations, we found clear flux decreases at the expected time of the eclipses. During the Suzaku observation, however, there were no clear signs for the predicted eclipses. The lapse of the predicted eclipses during the Suzaku observation can be explained by a contaminant source quite close to GRS 1747$-$312. When GRS 1747$-$312 is in the quiescent state, we observe X-rays from the contaminant source rather than from GRS 1747$-$312. However, we have no clear evidence for the contaminant source in our data. The lapse might also be explained by thick material ($N_{\rm H} > 10^{24}$ cm$^{-2}$ ) between the neutron star and the companion star, though the origin of the thick material is not clear.


title: Heterogeneous Multi-Tier Networks: Improper Signaling For Joint Rate-Energy Optimization

abstract: Wireless nodes in future communication systems need to overcome three barriers when compared to their transitional counterparts, namely to support significantly higher data rates, have long-lasting energy supplies and remain fully operational in interference-limited heterogeneous networks. This could be partially achieved by providing three promising features, which are radio frequency (RF) energy harvesting, improper Gaussian signaling and operating in full-duplex communication mode, i.e., transmit and receive at the same time within the same frequency band. In this paper, we consider these aspects jointly in a multi-antenna heterogeneous two-tier network. In this network, the users in the femto- cell are sharing the scarce resources with the cellular users in the macro-cell and have to cope with the interference from the macro-cell base station as well as the transmitter noise and residual self- interference (RSI) due to imperfect full-duplex operation. Interestingly enough, while these impairments are detrimental from the achievable rate perspective, they are beneficial from the energy harvesting aspect as they carry RF energy. In this paper, we consider this natural trade-off jointly and propose appropriate optimization problems for beamforming and optimal resource allocation. Moreover, various receiver structures are employed for both information detection (ID) and energy harvesting (EH) capabilities. The paper aims at characterizing the trade-off between the achievable rates and harvested energies. Rate and energy maximization problems are thoroughly investigated. Finally, the numerical illustrations demonstrate the impact of energy harvesting on the achievable rate performance.


title: Existence of Self-Cheeger sets on Riemannian manifolds

abstract: Let $(\mathcal{M},g)$ be a compact Riemannian manifold of dimension $N\geq 2$. We prove the existence of a family $(\Omega_\varepsilon)_{\varepsilon\in (0,\varepsilon_0)}$ of self-Cheeger sets in $(\mathcal{M},g)$ . The domains $\Omega_\varepsilon\subset\mathcal{M}$ are perturbations of geodesic balls of radius $\varepsilon$ centered at $p \in \mathcal{M}$, and in particular, if $p_0$ is a non-degenerate critical point of the scalar curvature of $g$, then the family $(\partial \Omega_\varepsilon)_{\varepsilon \in (0,\varepsilon_0)}$ constitutes a smooth foliation of a neighborhood of $p_0$.


title: A star forming ring around Kappa Ori 250 pc from the Sun

abstract: X-rays are a powerful probe of activity in early stages of star formation. They allow us to identify young stars even after they have lost the IR signatures of circumstellar disks and provide constraints on their distance. Here we report on XMM-Newton observations which detect 121 young stellar objects (YSOs) in two fields between L1641S and $\kappa$ Ori. These observations extend the Survey of Orion A with XMM and Spitzer (SOXS). The YSOs are contained in a ring of gas and dust apparent at millimeter wavelengths, and in far-IR and near-IR surveys. The X-ray luminosity function of the young stellar objects detected in the two fields indicates a distance of 250-280 pc, much closer than the Orion A cloud and similar to distance estimates of $\kappa$ Ori. We propose that the ring is a 5-8 pc diameter shell that has been swept up by $\kappa$ Ori. This ring contains several groups of stars detected by Spitzer and WISE including one surrounding the Herbig Ae/Be stars V1818 Ori. In this interpretation, the $\kappa$ Ori ring is one of several shells swept up by massive stars within the Orion Eridanus Superbubble, and is unrelated to the southern portion of Orion A / L1641 S.


title: A new approach to the Tarry-Escott problem

abstract: In this paper we describe a new method of obtaining ideal solutions of the well-known Tarry-Escott problem, that is, the problem of finding two distinct sets of integers $x_i,\;i=1,\,2,\,\dots,\,k+1$ and $y_i,\;i=1,\,2,\,\dots,\,k+1$ such that $ \sum_{i=1}^{k+1} x_i^r = \sum_{i=1}^{k+1} y_i^r,\;\;\;r = 1,\,2,\,\dots,\,k$, where $k$ is a given positive integer. When $k > 3$, only a limited number of parametric/ numerical ideal solutions of the Tarry-Escott problem are known. In this paper, by applying the new method mentioned above, we find several new parametric ideal solutions of the problem when $k \leq 7$. The ideal solutions obtained by this new approach are more general and very frequently, simpler than the ideal solutions obtained by the earlier methods. We also obtain new parametric solutions of certain diophantine systems that are closely related to the Tarry-Escott problem. These solutions are also more general and simpler than the solutions of these diophantine systems published earlier.


title: Sets of bounded remainder for the continuous irrational rotation on $[0,1)^2$

abstract: We study sets of bounded remainder for the two-dimensional continuous irrational rotation $(\{x_1+t\}, \{x_2+t\alpha \})_{t \geq 0}$ in the unit square. In particular, we show that for almost all $\alpha$ and every starting point $(x_1, x_2)$, every polygon $S$ with no edge of slope $\alpha$ is a set of bounded remainder. Moreover, every convex set $S$ whose boundary is twice continuously differentiable with positive curvature at every point is a bounded remainder set for almost all $\alpha$ and every starting point $(x_1, x_2)$. Finally we show that these assertions are, in some sense, best possible.


title: Approximation of a free Poisson process by systems of freely independent particles

abstract: Let $\sigma$ be a non-atomic, infinite Radon measure on $\mathbb R^d$, for example, $d\sigma(x)=z\,dx$ where $z>0$. We consider a system of freely independent particles $x_1,\dots,x_N$ in a bounded set $\Lambda\subset\mathbb R^d$, where each particle $x_i$ has distribution $\frac1{\sigma(\Lambda)}\,\sigma$ on $\Lambda$ and the number of particles, $N$, is random and has Poisson distribution with parameter $\sigma(\Lambda)$. If the particles were classically independent rather than freely independent, this particle system would be the restriction to $\Lambda$ of the Poisson point process on $\mathbb R^d$ with intensity measure $\sigma$. In the case of free independence, this particle system is not the restriction of the free Poisson process on $\mathbb R^d$ with intensity measure $\sigma$. Nevertheless, we prove that this is true in an approximative sense: if bounded sets $\Lambda^{(n)}$ ($n\in\mathbb N$) are such that $\Lambda^{(1)}\subset\Lambda^{(2)}\subset\Lambda^{(3)}\subset\dotsm$ and $\bigcup_{n=1}^\infty \Lambda^{(n)}=\mathbb R^d$, then the corresponding particle system in $\Lambda^{(n)}$ converges (as $n\to\infty$) to the free Poisson process on $\mathbb R^d$ with intensity measure $\sigma$. We also prove the following $N/V$-limit: Let $N^{(n)}$ be a determinstic sequence of natural numbers such that $\lim_{n\to\infty}N^{(n)}/\sigma(\Lambda^{(n)})=1$. Then the system of $N^{(n)}$ freely independent particles in $\Lambda^{(n)}$ converges (as $n\to\infty$) to the free Poisson process. We finally extend these results to the case of a free L\'evy white noise (in particular, a free L\'evy process) without free Gaussian part.


title: Group C*-algebras without the completely bounded approximation property

abstract: It is proved that: (1) The Fourier algebra A(G) of a simple Lie group G of real rank at least 2 with finite center does not have a multiplier bounded approximate unit. (2) The reduced C*-algebra of any lattice in a non-compact simple Lie group of real rank at least 2 with finite center does not have the completely bounded approximation property. Hence, the results obtained by J. de Canniere and the author for SO(n,1), n at least 2, and by M. Cowling for SU(n,1) do not generalize to simple Lie groups of real rank at least 2.


title: Phase-only signal reconstruction by MagnitudeCut

abstract: In this paper, we present a new algorithm, called MagnitudeCut, for recovering a signal from the phase of its Fourier transform. We casted our recovering problem into a new convex optimization problem, and then solved it by the block coordinate descent algorithm and the interior point algorithm, in which the iteration process consists of matrix vector product and inner product. We used the new method for reconstruction of a set of signal/image. The simulation results reveal that the proposed MagnitudeCut method can reconstruct the original signal with fewer sampling number of the phase information than that of the Greedy algorithm and iterative method under the same reconstruction error. Moreover, our algorithm can also reconstruct the symmetric image from its Fourier phase.


title: On the Estimation Performance and Convergence Rate of the Generalized Power Method for Phase Synchronization

abstract: An estimation problem of fundamental interest is that of phase synchronization, in which the goal is to recover a collection of phases using noisy measurements of relative phases. It is known that in the Gaussian noise setting, the maximum likelihood estimator (MLE) has an expected squared $\ell_2$-estimation error that is on the same order as the Cram\'er-Rao lower bound. Moreover, even though the MLE is an optimal solution to a non-convex quadratic optimization problem, it can be found with high probability using semidefinite programming (SDP), provided that the noise power is not too large. In this paper, we study the estimation and convergence performance of a recently-proposed low-complexity alternative to the SDP-based approach, namely, the generalized power method (GPM). Our contribution is twofold. First, we bound the rate at which the estimation error decreases in each iteration of the GPM and use this bound to show that all iterates---not just the MLE---achieve an estimation error that is on the same order as the Cram\'er-Rao bound. Our result holds under the least restrictive assumption on the noise power and gives the best provable bound on the estimation error known to date. It also implies that one can terminate the GPM at any iteration and still obtain an estimator that has a theoretical guarantee on its estimation error. Second, we show that under the same assumption on the noise power as that for the SDP-based method, the GPM will converge to the MLE at a linear rate with high probability. This answers a question raised in [3] and shows that the GPM is competitive in terms of both theoretical guarantees and numerical efficiency with the SDP-based method. At the heart of our convergence rate analysis is a new error bound for the non-convex quadratic optimization formulation of the phase synchronization problem, which could be of independent interest.


title: Graphene with vacancies: supernumerary zero modes

abstract: The density of states, $\varrho(E)$, of graphene is investigated within the tight binding (H\"uckel) approximation in the presence of vacancies. They induce a non-vanishing density of zero modes, $n_\text{zm}$, that act as midgap states: $\varrho(E)=n_\text{zm}\delta(E) + \text{smooth}$. As is well known, the actual number of zero modes per sample can in principle exceed the sublattice imbalance: $N_\text{zm}\geq |N_\text{A}-N_\text{B}|$, where $N_\text{A}$, $N_\text{B}$ denote the number of carbon atoms in each sublattice. In this work, we establish a stronger relation that is valid in the thermodynamic limit and that involves the concentration of zero-modes: $n_\text{zm}>|c_\text{A}-c_\text{B}|$, where $c_\text{A}$ and $c_\text{B}$ denote the concentration of vacancies per sublattice; in particular, $n_\text{zm}$ is non-vanishing even in the case of balanced disorder, $N_\text{A}/N_\text{B}=1$. Adopting terminology from benzoid graph theory, the excess modes associated with the current carrying backbone (percolation cluster) are called supernumerary. In the simplest cases such modes can be associated with structural elements of internal boundaries like dangling bonds. Our result suggest that the continuum limit of bipartite hopping models supports nontrivial "supernumerary" terms that escape the present continuum descriptions.


title: An Optimal Algorithm for l1-Heavy Hitters in Insertion Streams and Related Problems

abstract: We give the first optimal bounds for returning the $\ell_1$-heavy hitters in a data stream of insertions, together with their approximate frequencies, closing a long line of work on this problem. For a stream of $m$ items in $\{1, 2, \dots, n\}$ and parameters $0 < \epsilon < \phi \leq 1$, let $f_i$ denote the frequency of item $i$, i.e., the number of times item $i$ occurs in the stream. With arbitrarily large constant probability, our algorithm returns all items $i$ for which $f_i \geq \phi m$, returns no items $j$ for which $f_j \leq (\phi -\epsilon)m$, and returns approximations $\tilde{f}_i$ with $|\tilde{f}_i - f_i| \leq \epsilon m$ for each item $i$ that it returns. Our algorithm uses $O(\epsilon^{-1} \log\phi^{-1} + \phi^{-1} \log n + \log \log m)$ bits of space, processes each stream update in $O(1)$ worst-case time, and can report its output in time linear in the output size. We also prove a lower bound, which implies that our algorithm is optimal up to a constant factor in its space complexity. A modification of our algorithm can be used to estimate the maximum frequency up to an additive $\epsilon m$ error in the above amount of space, resolving Question 3 in the IITK 2006 Workshop on Algorithms for Data Streams for the case of $\ell_1$-heavy hitters. We also introduce several variants of the heavy hitters and maximum frequency problems, inspired by rank aggregation and voting schemes, and show how our techniques can be applied in such settings. Unlike the traditional heavy hitters problem, some of these variants look at comparisons between items rather than numerical values to determine the frequency of an item.


title: Permuting Incomplete Paired Data: A Novel Exact and Asymptotic Correct Randomization Test

abstract: Various statistical tests have been developed for testing the equality of means in matched pairs with missing values. However, most existing methods are commonly based on certain distributional assumptions such as normality, 0-symmetry or homoscedasticity of the data. The aim of this paper is to develop a statistical test that is robust against deviations from such assumptions and also leads to valid inference in case of heteroscedasticity or skewed distributions. This is achieved by applying a novel randomization approach. The resulting test procedure is not only shown to be asymptotically correct but is also finitely exact if the distribution of the data is invariant with respect to the considered randomization group. Its small sample performance is further studied in an extensive simulation study and compared to existing methods. Finally, an illustrative data example is analyzed.


title: Vacuum Rabi Splitting in Nanomechanical QED System with Nonlinear Resonator

abstract: Considering the intrinsic nonlinearity in a nanomechanical resonator coupled to a charge qubit, vacuum Rabi splitting effect is studied in a nanomechanical QED (qubit-resonator) system. A driven nonlinear Jaynes-Cummings model describes the dynamics of this qubit-resonator system. Using quantum regression theorem and master equation approach, we have calculated the two-time correlation spectrum analytically. In the weak driving limit, these analytical results clarify the influence of the driving strength and nonlinearity parameter on the correlation spectrum. Also, numerical calculations confirm these analytical results.


title: Structural, electronic and optical properties of well-known primary explosive: Mercury fulminate

abstract: Mercury Fulminate (MF) is one of the well-known primary explosives since 17th century and it has rendered invaluable service over many years. However, the correct molecular and crystal structures are determined recently after 300 years of its discovery. In the present study, we report pressure dependent structural, elastic, electronic and optical properties of MF. Non-local correction methods have been employed to capture the weak van der Waals interactions in layered and molecular energetic MF. Among the non-local correction methods tested, optB88-vdW method works well for the investigated compound. The obtained equilibrium bulk modulus reveals that MF is softer than the well known primary explosives Silver Fulminate (SF), silver azide and lead azide. MF exhibits anisotropic compressibility (b>a>c) under pressure, consequently the corresponding elastic moduli decrease in the following order: C22>C11>C33. The structural and mechanical properties suggest that MF is more sensitive to detonate along c-axis (similar to RDX) due to high compressibility of Hg...O non-bonded interactions along that axis. Electronic structure and optical properties were calculated including spin-orbit (SO) interactions using full potential linearized augmented plane wave method within recently developed Tran-Blaha modified Becke-Johnson (TB-mBJ) potential. The calculated TB-mBJ electronic structures of SF and MF show that these compounds are indirect bandgap insulators. Also, SO coupling is found to be more pronounced for 4d and 5d-states of Ag and Hg atoms of SF and MF, respectively. Partial density of states and electron charge density maps were used to describe the nature of chemical bonding. Ag-C bond is more directional than Hg-C bond which makes SF to be more unstable than MF. The effect of SO coupling on optical properties has also been studied and found to be significant for both of the compounds.


title: A class of globally solvable Markovian quadratic BSDE systems and applications

abstract: We establish existence and uniqueness for a wide class of Markovian systems of backward stochastic differential equations (BSDE) with quadratic nonlinearities. This class is characterized by an abstract structural assumption on the generator, an a-priori local-boundedness property, and a locally-H\"older-continuous terminal condition. We present easily verifiable sufficient conditions for these assumptions and treat several applications, including stochastic equilibria in incomplete financial markets, stochastic differential games, and martingales on Riemannian manifolds.


title: Alternative multiplications and non-associativity in physics

abstract: Some physical systems like the quantum mechanics with magnetic charges or field theoretical models appearing in the context of string theory are formulated in terms of non-associative algebras. Hence, demand non-associative star products for its description. However, unlike its associative counterpart the situation with non-associative star products is quite unclear. At the moment it is not evident which condition should be used instead of the associativity in the definition of the non-associative star multiplication. In mathematics, the natural generalization of the associativity is the requirement that the multiplication should be alternative, the associator of three elements should vanish if each two of them are equal. We show that for the alternative closed star product the integrated associator vanishes. In particular, it means that on-shell there will be no violation of associativity in string scattering amplitudes. We discuss the examples of the alternative and closed star products. In particular, it is constructed the star product corresponding to the algebra of octonions.


title: The data processing pipeline for the Herschel SPIRE Fourier Transform Spectrometer

abstract: We present the data processing pipeline to generate calibrated data products from the Spectral and Photometric Imaging Receiver (SPIRE) imaging Fourier Transform Spectrometer on the Herschel Space Observatory. The pipeline processes telemetry from SPIRE observations and produces calibrated spectra for all resolution modes. The spectrometer pipeline shares some elements with the SPIRE photometer pipeline, including the conversion of telemetry packets into data timelines and calculation of bolometer voltages. We present the following fundamental processing steps unique to the spectrometer: temporal and spatial interpolation of the scan mechanism and detector data to create interferograms; Fourier transformation; apodization; and creation of a data cube. We also describe the corrections for various instrumental effects including first- and second-level glitch identification and removal, correction of the effects due to emission from the Herschel telescope and from within the spectrometer instrument, interferogram baseline correction, temporal and spatial phase correction, non-linear response of the bolometers, and variation of instrument performance across the focal plane arrays. Astronomical calibration is based on combinations of observations of standard astronomical sources and regions of space known to contain minimal emission.


title: Evolution of electron Fermi surface with doping in cobaltates

abstract: The notion of the electron Fermi surface is one of the characteristic concepts in the field of condensed matter physics, and it plays a crucial role in the understanding of the physical properties of doped Mott insulators. Based on the t-J model, we study the nature of the electron Fermi surface in the cobaltates, and qualitatively reproduce the essential feature of the evolution of the electron Fermi surface with doping. It is shown that the underlying hexagonal electron Fermi surface obeys Luttinger's theorem. The theory also predicts a Fermi-arc phenomenon at the low-doped regime, where the region of the hexagonal electron Fermi surface along the \Gamma-K direction is suppressed by the electron self-energy, and then six disconnected Fermi arcs located at the region of the hexagonal electron Fermi surface along the \Gamma-M direction emerge. However, this Fermi-arc phenomenon at the low-doped regime weakens with the increase of doping.


title: Bumblebees minimize control challenges by combining active and passive modes in unsteady winds

abstract: The natural wind environment that volant insects encounter is unsteady and highly complex, posing significant flight control and stability challenges. Unsteady airflows can range from structured chains of discrete vortices shed in the wake of an object to fully developed chaotic turbulence. It is critical to understand the flight control strategies insects employ to safely navigate in natural environments. We combined experiments on free flying bumblebees with high fidelity numerical simulations and lower order modeling to identify the salient mechanics that mediate insect flight in unsteady winds. We trained bumblebees to fly upwind towards an artificial flower in a wind tunnel under steady wind and in a von Karman street (23Hz) formed in the wake of a cylinder. The bees displayed significantly higher movement in the unsteady vortex street compared to steady winds. Correlation analysis revealed that at lower frequencies, less than 10 Hz, in both steady and unsteady winds the bees mediated lateral movement with body roll, typical casting motion. At higher frequencies in unsteady winds there was a negative correlation between body roll and lateral accelerations. Numerical simulations of a bumblebee in similar conditions permitted the separation of the passive and active components of the flight trajectories. Comparison between the free-flying and numerical bees revealed a novel mechanism that enables bees to passively ride out high frequency perturbations while performing active maneuvers and corrections at lower frequencies. The capacity of maintaining stability by combining passive and active modes at different timescales provides a viable means for volant animals and machines to tackle the control challenges posed by complex airflows.


title: Isotropic geometry of graph surfaces associated with product production functions in economics

abstract: A production function is a mathematical formalization in economics which denotes the relations between the output generated by a firm, an industry or an economy and the inputs that have been used in obtaining it. In this paper, we study the product production functions of 2 variables in terms of the geometry of their associated graph surfaces in the isotropic 3-space I^{3}. In particular, we derive several classification results for the graph surfaces of product production functions in I^{3} with constant curvature.


title: Segmental Recurrent Neural Networks for End-to-end Speech Recognition

abstract: We study the segmental recurrent neural network for end-to-end acoustic modelling. This model connects the segmental conditional random field (CRF) with a recurrent neural network (RNN) used for feature extraction. Compared to most previous CRF-based acoustic models, it does not rely on an external system to provide features or segmentation boundaries. Instead, this model marginalises out all the possible segmentations, and features are extracted from the RNN trained together with the segmental CRF. In essence, this model is self-contained and can be trained end-to-end. In this paper, we discuss practical training and decoding issues as well as the method to speed up the training in the context of speech recognition. We performed experiments on the TIMIT dataset. We achieved 17.3 phone error rate (PER) from the first-pass decoding --- the best reported result using CRFs, despite the fact that we only used a zeroth-order CRF and without using any language model.


title: Numerically Fitting The Electron Fermi Energy and The Electron Fraction in A Neutron Star

abstract: Based on the basic definition of Fermi energy of degenerate and relativistic electrons, we obtain a special solution to electron Fermi energy, $E_{\rm F}(e)$, and express $E_{\rm F}(e)$ as a function of electron fraction, $Y_{e}$, and matter density, $\rho$. Several useful analytical formulae for $Y_{e}$ and $\rho$ within classical models and the work of Dutra et al. 2014 (Type-2) in relativistic mean field theory are obtained using numerically fitting. When describing the mean-field Lagrangian, density, we adopt the TMA parameter set, which is remarkably consistent with with the updated astrophysical observations of neutron stars. Due to the importance of the density dependence of the symmetry energy, $S$, in nuclear astrophysics, a brief discussion on the symmetry parameters $S_v$ and $L$ (the slope of $S$) is presented. Combining these fit formulae with boundary conditions for different density regions, we can evaluate the value of $E_{\rm F}(e)$ in any given matter density, and obtain a schematic diagram of $E_{\rm F}(e)$ as a continuous function of $\rho$. Compared with previous study on the electron Fermi energy in other models, our methods of calculating $E_{\rm F}(e)$ are more simple and convenient, and can be universally suitable for the relativistic electron regions in the circumstances of common neutron stars. We have deduced a general expression of $E_{\rm F}(e)$ and $n_{e}$, which could be used to indirectly test whether one EoS of a NS is correct in our future studies on neutron star matter properties. Since URCA reactions are expected in the center of a massive star due to high-value electron Fermi energy and electron fraction, this study could be useful in the future studies on the NS thermal evolution.


title: Non equilibrium optical properties in semiconductors from first--principles: a combined theoretical and experimental study of bulk silicon

abstract: The calculation of the equilibrium optical properties of bulk silicon by using the Bethe--Salpeter equation solved in the Kohn--Sham basis represents a cornerstone in the development of an ab--initio approach to the optical and electronic properties of materials. Nevertheless calculations of the {\em transient} optical spectrum using the same efficient and successful scheme are scarce. We report, here, a joint theoretical and experimental study of the transient reflectivity spectrum of bulk silicon. Femtosecond transient reflectivity is compared to a parameter--free calculation based on the non--equilibrium Bethe--Salpeter equation. By providing an accurate description of the experimental results we disclose the different phenomena that determine the transient optical response of a semiconductor. We give a parameter--free interpretation of concepts like bleaching, photo--induced absorption and stimulated emission, beyond the Fermi golden rule. We also introduce the concept of optical gap renormalization, as a generalization of the known mechanism of band gap renormalization. The present scheme successfully describes the case of bulk silicon, showing its universality and accuracy.


title: Finding the Nucleoli of Large Cooperative Games: A Disproof with Counter-Example

abstract: Nguyen and Thomas (2016) claimed that they have found a method to compute the nucleoli of games with more than $50$ players using nested linear programs (LP). Unfortunately, this claim is false. They incorrectly applied the indirect proof by "$A \land \neg B$ implies $A \land \neg A$" to conclude that "if $A$ then $B$" is valid. In fact, they prove that a truth implies a falsehood. As established by Meinhardt (2015a), this is a wrong statement. Therefore, instead of giving a proof of their main Theorem 4b, they give a disproof. It comes as no surprise to us that the flow game example presented by these authors to support their arguments is obviously a counter-example of their algorithm. We show that the computed solution by this algorithm is neither the nucleolus nor a core element of the flow game. Moreover, the stopping criterion of all proposed methods is wrong, since it does not satisfy one of Kohlberg's properties (cf. Kohlberg (1971)). As a consequence, none of these algorithms is robust.


title: On the vortex filament conjecture for Euler flows

abstract: In this paper, we study the evolution of a vortex filament in an incompressible ideal fluid. Under the assumption that the vorticity is concentrated along a smooth curve in $\mathbb{R}^3$, we prove that the curve evolves to leading order by binormal curvature flow. Our approach combines new estimates on the distance of the corresponding Hamiltonian-Possion structures with stability estimates recently developed in Ref. 15.


title: Possible $B^{(\ast)} \bar{K}$ hadronic molecule state

abstract: In the present work, we estimate the decays of the $X(5568)$ and $X(5616)$ in a $B \bar{K}$ and a $B^\ast \bar{K}$ $S$-wave hadronic molecule scenarios, respectively, which may corresponding to the structure observed by D0 Collaboration. Our estimation indicates both $B\bar{K}$ and $B^\ast \bar{K}$ hadronic molecule decay widths could explain the experimental data in a proper model parameter range.


title: Fetishizing Food in Digital Age: #foodporn Around the World

abstract: What food is so good as to be considered pornographic? Worldwide, the popular #foodporn hashtag has been used to share appetizing pictures of peoples' favorite culinary experiences. But social scientists ask whether #foodporn promotes an unhealthy relationship with food, as pornography would contribute to an unrealistic view of sexuality. In this study, we examine nearly 10 million Instagram posts by 1.7 million users worldwide. An overwhelming (and uniform across the nations) obsession with chocolate and cake shows the domination of sugary dessert over local cuisines. Yet, we find encouraging traits in the association of emotion and health-related topics with #foodporn, suggesting food can serve as motivation for a healthy lifestyle. Social approval also favors the healthy posts, with users posting with healthy hashtags having an average of 1,000 more followers than those with unhealthy ones. Finally, we perform a demographic analysis which shows nation-wide trends of behavior, such as a strong relationship (r=0.51) between the GDP per capita and the attention to healthiness of their favorite food. Our results expose a new facet of food "pornography", revealing potential avenues for utilizing this precarious notion for promoting healthy lifestyles.


title: Neutrino masses and ordering via multimessenger astronomy

abstract: We define the theoretical framework and deduce the conditions under which multi-messenger astronomy can provide useful information about neutrino masses and their ordering. The framework uses time differences between the arrival of neutrinos and the other light messenger, i.e. the graviton, emitted in astrophysical catastrophes. We also provide a preliminary feasibility study elucidating the experimental reach and challenges for planned neutrino detectors such as Hyper-Kamiokande as well as future several megaton detectors. This study shows that future experiments can be useful in testing independently the cosmological bounds on absolute neutrino masses. Concretely the success of such measurements depends crucially on the available rate of astrophysical events and further requires development of high resolution timing besides the advocated need of megaton size detectors.


title: Classifying Unidentified Gamma-ray Sources

abstract: During its first 2 years of mission the Fermi-LAT instrument discovered more than 1,800 gamma-ray sources in the 100 MeV to 100 GeV range. Despite the application of advanced techniques to identify and associate the Fermi-LAT sources with counterparts at other wavelengths, about 40% of the LAT sources have no a clear identification remaining "unassociated". The purpose of my Ph.D. work has been to pursue a statistical approach to identify the nature of each Fermi-LAT unassociated source. To this aim, we implemented advanced machine learning techniques, such as logistic regression and artificial neural networks, to classify these sources on the basis of all the available gamma-ray information about location, energy spectrum and time variability. These analyses have been used for selecting targets for AGN and pulsar searches and planning multi-wavelength follow-up observations. In particular, we have focused our attention on the search of possible radio-quiet millisecond pulsar (MSP) candidates in the sample of the Fermi-LAT unidentified sources. These objects have not yet been detected but their discovery would have a formidable impact for our understanding of the MSP gamma-ray emission mechanism.


title: Cost Optimisation of an Instrument Suite at an Accelerator-Driven Spallation Source

abstract: This artcile presents an optimisation of performance and cost of neutron scattering instrumentation at the European Spallation Source. This is done by trading detailed cost functions against beam transmission functions in a multi-dimensional, yet simple, parameter space. On the one hand, the neutron guide cost increases as a power of the desired beam divergence, and inversely with the minimum wavelength, due to the supermirror coating needed. On the other hand, the more neutrons are transported to the instrument the greater are the shielding costs to deal with the gamma rays that result from the eventual absorption of the neutrons. There are additional factors in that many of the parameters defining the neutron guide geometry are continuous variables, and the straightness of the guide increases the transmission of high energy spallation products, which affect the specifications of particularly heavy hardware, such as heavy shutters and additional shielding, beam stops etc. Over the suite of 16 instruments, a maximum potential cost saving is identified in the vicinity of 20 M\euro \,, with negligible penalties on instrument performance and bandwidth, by trading minimum neutron wavelength $\lambda$ against total optical costs and shielding costs. Finally, these cost calculations are shown to be consistent with those at real world, existing spallation source facilities.


title: Damped vacuum states of light

abstract: We consider one-dimensional propagation of quantum light in the presence of a block of material, with a full account of dispersion and absorption. The electromagnetic zero-point energy for some frequencies is damped (suppressed) by the block below the free-space value, while for other frequencies it is increased. We also calculate the regularized (Casimir) zero-point energy at each frequency and find that it too is damped below the free-space value (zero) for some frequencies. The total Casimir energy is positive.


title: M-curves and symmetric products

abstract: Let $(X , \sigma)$ be a geometrically irreducible smooth projective M-curve of genus $g$ defined over the field of real numbers. We prove that the $n$-th symmetric product of $(X , \sigma)$ is an M-variety for $n=2 ,3$ and $n\geq 2g -1$.


title: Oracle Estimation of a Change Point in High Dimensional Quantile Regression

abstract: In this paper, we consider a high-dimensional quantile regression model where the sparsity structure may differ between two sub-populations. We develop $\ell_1$-penalized estimators of both regression coefficients and the threshold parameter. Our penalized estimators not only select covariates but also discriminate between a model with homogeneous sparsity and a model with a change point. As a result, it is not necessary to know or pretest whether the change point is present, or where it occurs. Our estimator of the change point achieves an oracle property in the sense that its asymptotic distribution is the same as if the unknown active sets of regression coefficients were known. Importantly, we establish this oracle property without a perfect covariate selection, thereby avoiding the need for the minimum level condition on the signals of active covariates. Dealing with high-dimensional quantile regression with an unknown change point calls for a new proof technique since the quantile loss function is non-smooth and furthermore the corresponding objective function is non-convex with respect to the change point. The technique developed in this paper is applicable to a general M-estimation framework with a change point, which may be of independent interest. The proposed methods are then illustrated via Monte Carlo experiments and an application to tipping in the dynamics of racial segregation.


title: Interference effect between $\phi$ and $\Lambda(1520)$ production channels in the $\gamma p \rightarrow K^+K^-p$ reaction near threshold

abstract: The $\phi$-$\Lambda(1520)$ interference effect in the $\gamma p\to K^+K^-p$ reaction has been measured for the first time in the energy range from 1.673 to 2.173 GeV. The relative phases between $\phi$ and $\Lambda(1520)$ production amplitudes were obtained in the kinematic region where the two resonances overlap. The measurement results support strong constructive interference when $K^+K^-$ pairs are observed at forward angles, but destructive interference for proton emission at forward angles. Furthermore, the observed interference effect does not account for the $\sqrt{s}=2.1$ GeV bump structure in forward differential cross sections for $\phi$ photoproduction. This fact suggests possible exotic structures such as a hidden-strangeness pentaquark state, a new Pomeron exchange and rescattering processes via other hyperon states.


title: Center of the quantum affine vertex algebra in type A

abstract: We consider the quantum vertex algebra associated with the double Yangian in type A as defined by Etingof and Kazhdan. We show that its center is a commutative associative algebra and construct algebraically independent families of topological generators of the center at the critical level.


title: Linearly degenerate hierarchies of quasiclassical SDYM type

abstract: We demonstrate that SDYM equations for the Lie algebra of one-dimensional vector fields represent a natural reduction in the framework of general linearly degenerate dispersionless hierarchy. We define the reduction in terms of wave functions, introduce generating relation, Lax-Sato equations and the dressing scheme for the reduced hierarchy. Multidimensional case is also discussed.


title: A Solution Theory for a General Class of SPDEs

abstract: In this article we present a way of treating stochastic partial differential equations with multiplicative noise by rewriting them as stochastically perturbed evolutionary equations in the sense of \cite{picardbook}, where a general solution theory for deterministic evolutionary equations has been developed. This allows us to present a unified solution theory for a general class of SPDEs which we believe has great potential for further generalizations. We will show that many standard stochastic PDEs fit into this class as well as many other SPDEs such as the stochastic Maxwell equation and time-fractional stochastic PDEs with multiplicative noise on sub-domains of $\Rd$. The approach is in spirit similar to the approach in \cite{dapratozabczyk}, but complementing it in the sense that it does not involve semi-group theory and allows for an effective treatment of coupled systems of SPDEs. In particular, the existence of a (regular) fundamental solution or Green's function is not required.


title: High pressure structural, electronic, and optical properties of polymorphic InVO4 phases

abstract: In the present work, we report a detailed density functional theory calculation on polymorphic InVO$_4$ phases by means of projector augmented wave method. The computed first-order structural phase transformation from orthorhombic \emph{(Cmcm)} to monoclinic \emph{(P2/c)} structure is found to occur around 5.6 GPa along with a large volume collapse of 16.6$\%$, which is consistent with previously reported experimental data. This transformation also leads to an increase in the coordination number of vanadium atom from 4 to 6. The computed equilibrium and high pressure structural properties of both InVO$_4$ phases, including unit cell parameters, equation of state, and bulk moduli, are in good agreement with the available experimental data. In addition, compressibility is found to be highly anisotropic and the \emph{b}-axis being more compressible than the other for both the structures. Electronic band structures for both the phases were calculated, and the band gap for orthorhombic and monoclinic InVO$_4$ are found to be 4.02 and 1.67 eV, respectively, within the Tran-Blaha Modified Becke-Johnson potential as implemented in linearized augmented planewave method. We further examined the optical properties such as dielectric function, refractive index, and absorption spectra for both the structures. From the implications of these results, it can be proposed that the high pressure InVO$_4$ phase can be more useful than orthorhombic phase for photo catalytic applications.


title: An Extension Theorem for convex functions of class $C^{1,1}$ on Hilbert spaces

abstract: Let $\mathbb{H}$ be a Hilbert space, $E \subset \mathbb{H}$ be an arbitrary subset and $f: E \rightarrow \mathbb{R}, \: G: E \rightarrow \mathbb{H}$ be two functions. We give a necessary and sufficient condition on the pair $(f,G)$ for the existence of a \textit{convex} function $F\in C^{1,1}(\mathbb{H})$ such that $F=f$ and $\nabla F =G$ on $E$. We also show that, if this condition is met, $F$ can be taken so that $\textrm{Lip}(\nabla F) = \textrm{Lip}(G)$. We give a geometrical application of this result, concerning interpolation of sets by boundaries of $C^{1,1}$ convex bodies in $\mathbb{H}$. Finally, we give a counterexample to a related question concerning smooth convex extensions of smooth convex functions with derivatives which are not uniformly continuous.


title: 3-nets realizing a diassociative loop in a projective plane

abstract: A \textit{$3$-net} of order $n$ is a finite incidence structure consisting of points and three pairwise disjoint classes of lines, each of size $n$, such that every point incident with two lines from distinct classes is incident with exactly one line from each of the three classes. The current interest around $3$-nets (embedded) in a projective plane $PG(2,K)$, defined over a field $K$ of characteristic $p$, arose from algebraic geometry. It is not difficult to find $3$-nets in $PG(2,K)$ as far as $0<p\le n$. However, only a few infinite families of $3$-nets in $PG(2,K)$ are known to exist whenever $p=0$, or $p>n$. Under this condition, the known families are characterized as the only $3$-nets in $PG(2,K)$ which can be coordinatized by a group. In this paper we deal with $3$-nets in $PG(2,K)$ which can be coordinatized by a diassociative loop $G$ but not by a group. We prove two structural theorems on $G$. As a corollary, if $G$ is commutative then every non-trivial element of $G$ has the same order, and $G$ has exponent $2$ or $3$. We also discuss the existence problem for such $3$-nets.


title: Precision Measurement of the Radiative $\Beta$ Decay of the Free Neutron

abstract: The standard model predicts that, in addition to a proton, an electron, and an antineutrino, a continuous spectrum of photons is emitted in the $\beta$ decay of the free neutron. We report on the RDK II experiment which measured the photon spectrum using two different detector arrays. An annular array of bismuth germanium oxide scintillators detected photons from 14 to 782~keV. The spectral shape was consistent with theory, and we determined a branching ratio of 0.00335 $\pm$ 0.00005 [stat] $\pm$ 0.00015 [syst]. A second detector array of large area avalanche photodiodes directly detected photons from 0.4 to 14~keV. For this array, the spectral shape was consistent with theory, and the branching ratio was determined to be 0.00582 $\pm$ 0.00023 [stat] $\pm$ 0.00062 [syst]. We report the first precision test of the shape of the photon energy spectrum from neutron radiative decay and a substantially improved determination of the branching ratio over a broad range of photon energies.


title: Caching Incentive Design in Wireless D2D Networks: A Stackelberg Game Approach

abstract: Caching in wireless device-to-device (D2D) networks can be utilized to offload data traffic during peak times. However, the design of incentive mechanisms is challenging due to the heterogeneous preference and selfish nature of user terminals (UTs). In this paper, we propose an incentive mechanism in which the base station (BS) rewards those UTs that share contents with others using D2D communication. We study the cost minimization problem for the BS and the utility maximization problem for each UT. In particular, the BS determines the rewarding policy to minimize his total cost, while each UT aims to maximize his utility by choosing his caching policy. We formulate the conflict among UTs and the tension between the BS and the UTs as a Stackelberg game. We show the existence of the equilibrium and propose an iterative gradient algorithm (IGA) to obtain the Stackelberg Equilibrium. Extensive simulations are carried out to evaluate the performance of the proposed caching scheme and comparisons are drawn with several baseline caching schemes with no incentives. Numerical results show that the caching scheme under our incentive mechanism outperforms other schemes in terms of the BS serving cost and the utilities of the UTs.


title: Landau-Zener-St\"uckelberg-Majorana lasing in circuit QED

abstract: We demonstrate amplification (and attenuation) of a probe signal by a driven two-level quantum system in the Landau-Zener-St\"{u}ckelberg-Majorana regime by means of an experiment, in which a superconducting qubit was strongly coupled to a microwave cavity, in a conventional arrangement of circuit quantum electrodynamics. Two different types of flux qubit, specifically a conventional Josephson junctions qubit and a phase-slip qubit, show similar results, namely, lasing at the working points where amplification takes place. The experimental data are explained by the interaction of the probe signal with Rabi-like oscillations. The latter are created by constructive interference of Landau-Zener-St\"{u}ckelberg-Majorana (LZSM) transitions during the driving period of the qubit. A detailed description of the occurrence of these oscillations and a comparison of obtained data with both analytic and numerical calculations are given.


title: Polynomial ideals from a nonlinear viewpoint

abstract: Classes of homogeneous polynomials between Banach spaces have been studied in the last three decades from the perspective of the so-called ideal property: if a polynomial P belongs to a class Q, then the composition u o P o v of P with linear operators u and v belongs to Q as well. In an attempt to explore the nonlinearity of the subject in a more consistent way, and taking into account recent results in the field, in this paper we propose the study of classes of homogeneous polynomials that are stable under the composition with homogeneous polynomials. Some importante classes justify the study of the intermediate concept of classes of polynomials Q such that if P belongs to Q, u is a linear operator and Q is a homogeneous polynomial, then u o P o Q belongs to Q.


title: Electrons diffraction scattering on a traveling wave - "Inelastic Kapitza-Dirac effect"

abstract: In this paper conceptual points regarding electrons elastic (Kapitza-Dirac effect) and inelastic diffraction effect on the different type slowed electromagnetic wave structures/light gratings are considered. From the unified point of view it is analyzed the main works on this subject for last four decades in chronological order, pointing out the essential peculiarity inherent in induced Cherenkov, Compton, and undulator/wiggler processes too. This review article has also a goal to resolve confusion in scientific literature connected with the recently appeared paper [6] relating the electrons diffraction on a traveling wave in a dielectric medium.


title: Class of exact memory-kernel master equations

abstract: A well-known situation in which a non-Markovian dynamics of an open quantum system $S$ arises is when this is coherently coupled to an auxiliary system $M$ in contact with a Markovian bath. In such cases, while the joint dynamics of $S$-$M$ is Markovian and obeys a standard (bipartite) Lindblad-type master equation (ME), this is in general not true for the reduced dynamics of $S$. Furthermore, there are several instances (\eg the dissipative Jaynes-Cummings model) in which a {\it closed} ME for the $S$'s state {\it cannot} even be worked out. Here, we find a class of bipartite Lindblad-type MEs such that the reduced ME of $S$ can be derived exactly and in a closed form for any initial product state of $S$-$M$. We provide a detailed microscopic derivation of our result in terms of a mapping between two collision models


title: Polarization angle swings in blazars: The case of 3C 279

abstract: Over the past few years, several occasions of large, continuous rotations of the electric vector position angle (EVPA) of linearly polarized optical emission from blazars have been reported. These events are often coincident with high energy gamma-ray flares and they have attracted considerable attention, as they could allow one to probe the magnetic field structure in the gamma-ray emitting region of the jet. The flat-spectrum radio quasar 3C279 is one of the most prominent examples showing this behaviour. Our goal is to study the observed EVPA rotations and to distinguish between a stochastic and a deterministic origin of the polarization variability. We have combined multiple data sets of R-band photometry and optical polarimetry measurements of 3C279, yielding exceptionally well-sampled flux density and polarization curves that cover a period of 2008-2012. Several large EVPA rotations are identified in the data. We introduce a quantitative measure for the EVPA curve smoothness, which is then used to test a set of simple random walk polarization variability models against the data. 3C279 shows different polarization variation characteristics during an optical low-flux state and a flaring state. The polarization variation during the flaring state, especially the smooth approx. 360 degrees rotation of the EVPA in mid-2011, is not consistent with the tested stochastic processes. We conclude that during the two different optical flux states, two different processes govern the polarization variation, possibly a stochastic process during the low-brightness state and a deterministic process during the flaring activity.


title: Dopant-controlled single-electron pumping through a metallic island

abstract: We investigate a hybrid metallic island / single dopant electron pump based on fully-depleted silicon on insulator technology. Electron transfer between the central metallic island and the leads is controlled by resonant tunneling through single phosphorus dopants in the barriers. Top gates above the barriers are used control the resonance conditions. Applying radio frequency signals to the gates, non-adiabatic quantized electron pumping is achieved. A simple deterministic model is presented and confirmed by comparing measurements with simulations.


title: An Introduction to L\'evy and Feller Processes. Advanced Courses in Mathematics - CRM Barcelona 2014

abstract: These lecture notes are an extended version of my lectures on L\'evy and L\'evy-type (Feller) processes given at the "Second Barcelona Summer School on Stochastic Analysis" 2014 organized by the Centre de Recerca Matemaatica (CRM). The lectures are aimed at advanced graduate and PhD students. In order to read these notes, one should have sound knowledge of measure theoretic probability theory and some background in stochastic processes, as it is covered in my books "Measures, Integals and Martingales" (Cambridge University Press) and "Brownian Motion" (de Gruyter). My purpose in these lectures is to give an introduction to Levy processes, and to show how one can extend this approach to space inhomogeneous processes which behave locally like L\'evy processes: L\'evy-type or Feller processes. These course notes will be published, together Davar Khoshnevisan's notes on "Invariance and Comparison Principles for Parabolic SPDEs" as "From L\'evy-Type Processes to Parabolic SPDEs" by the CRM, Barcelona and Birk\"auser, Cham 2017 (ISBN: 978-3-319-34119-4). The arXiv-version and the published version may differ in layout, pagination and wording, but not in content


title: Exact density profiles and symmetry classification for strongly interacting multi-component Fermi gases in tight waveguides

abstract: We consider a mixture of one-dimensional strongly interacting Fermi gases up to six components, subjected to a longitudinal harmonic confinement. In the limit of infinitely strong repulsions we provide an exact solution which generalizes the one for the two-component mixture. We show that an imbalanced mixture under harmonic confinement displays partial spatial separation among the components, with a structure which depends on the relative population of the various components. Furthermore, we provide a symmetry characterization of the ground and excited states of the mixture introducing and evaluating a suitable operator, namely the conjugacy class sum. We show that, even under external confinement, the gas has a definite symmetry which corresponds to the most symmetric one compatible with the imbalance among the components. This generalizes the predictions of the Lieb-Mattis theorem for a fermionic mixture with more than two components.


title: Magnetic properties of Fe$_5$SiB$_2$ and its alloys with P, S, and Co

abstract: Fe$_5$SiB$_2$ has been synthesized and magnetic measurements have been carried out, revealing that M$_{\text{sat}}$ = 0.92 MA/m at T = 300 K. The M vs T curve shows a broad peak around T = 160 K. The anisotropy constant, K$_1$, estimated at T = 300 K, is 0.25 MJ/m$^3$. Theoretical analysis of Fe$_5$SiB$_2$ system has been carried out and extended to the full range of Fe$_5$Si$_{1-x}$P$_x$B$_2$, Fe$_5$P$_{1-x}$S$_x$B$_2$, and (Fe$_{1-x}$Co$_x$)$_5$SiB$_2$ compositions. The electronic band structures have been calculated using the Full-Potential Local-Orbital Minimum-Basis Scheme (FPLO-14). The calculated total magnetic moments are 9.20, 9.15, 9.59 and 2.42$\mu_B$ per formula units of Fe$_5$SiB$_2$, Fe$_5$PB$_2$, Fe$_5$SB$_2$, and Co$_5$SiB$_2$, respectively. In agreement with experiment, magnetocrystalline anisotropy energies (MAE's) calculated for T = 0 K changes from a negative (easy-plane) anisotropy -0.28 MJ/m$^3$ for Fe$_5$SiB$_2$ to the positive (easy-axis) anisotropy 0.35 MJ/m$^3$ for Fe$_5$PB$_2$. Further increase of the number of p-electrons in Fe$_5$P$_{1-x}$S$_x$B$_2$ leads to an increase of MAE up to 0.77 MJ/m$^3$ for the hypothetical Fe$_5$P$_{0.4}$S$_{0.6}$B$_2$ composition. Volume variation and fixed spin moment calculations (FSM) performed for Fe$_5$SiB$_2$ show an inverse relation between MAE and magnetic moment in the region down to about 15\% reduction of the spin moment. The alloying of Fe$_5$SiB$_2$ with Co is proposed as a practical realization of magnetic moment reduction, which ought to increase MAE. MAE calculated in virtual crystal approximation (VCA) for a full range of (Fe$_{1-x}$Co$_x$)$_5$SiB$_2$ compositions reaches the maximum value of 1.16 MJ/m$^3$ at Co concentration x = 0.3, with the magnetic moment 7.75$\mu_B$ per formula unit. Thus, (Fe$_{0.7}$Co$_{0.3}$)$_5$SiB$_2$ is suggested as a candidate for a rare-earth free permanent magnet.


title: On the hardness of switching to a small number of edges

abstract: Seidel's switching is a graph operation which makes a given vertex adjacent to precisely those vertices to which it was non-adjacent before, while keeping the rest of the graph unchanged. Two graphs are called switching-equivalent if one can be made isomorphic to the other one by a sequence of switches. Jel\'inkov\'a et al. [DMTCS 13, no. 2, 2011] presented a proof that it is NP-complete to decide if the input graph can be switched to contain at most a given number of edges. There turns out to be a flaw in their proof. We present a correct proof. Furthermore, we prove that the problem remains NP-complete even when restricted to graphs whose density is bounded from above by an arbitrary fixed constant. This partially answers a question of Matou\v{s}ek and Wagner [Discrete Comput. Geom. 52, no. 1, 2014].


title: Revisiting the R$\nu$MDM Models

abstract: Combining neutrino mass generation and a dark matter candidate in a unified model has always been intriguing. We revisit the class of R$\nu$MDM models, which incorporate minimal dark matter in radiative neutrino mass models based on the one-loop ultraviolet completions of the Weinberg operator. The possibility of an exact accidental $Z_2$ is completely ruled out in this scenario. We study the phenomenology of one of the models with an approximate $Z_2$ symmetry. In addition to the Standard Model particles, it contains two real scalar quintuplets, one vector- like quadruplet fermion and a fermionic quintuplet. The neutral component of the fermionic quintuplet serves as a good dark matter candidate which can be tested by the future direct and indirect detection experiments. The constraints from flavor physics and electroweak-scale naturalness are also discussed.


title: Ordering phenomena in a heterostructure of frustrated and unfrustrated triangular-lattice Ising layers

abstract: We study critical and magnetic properties of a bilayer Ising system consisting of two triangular planes A and B, with the antiferromagnetic (AF) coupling $J_{\rm A}$ and the ferromagnetic (FM) one $J_{\rm B}$ for the respective layers, which are coupled by the interlayer interaction $J_{\rm AB}$ by using Monte Carlo simulations. When $J_{\rm A}$ and $J_{\rm B}$ are of the same order, the unfrustrated FM plane orders first at a high temperature $T_{c1} \sim J_{\rm B}$. The spontaneous FM order then exerts influence on the other frustrated AF plane as an effective magnetic field, which subsequently induces a ferrimagnetic order in this plane at low temperatures below $T_{c2}$. When short-range order is developed in the AF plane while the influence of the FM plane is still small, there appears a preemptive Berezinskii-Kosterlitz-Thouless-like pseudocritical crossover regime just above the ferrimagnetic phase transition point, where the short-distance behavior up to a rather large length scale exponentially diverging in $\propto \JA / T$ is controlled by a line of Gaussian fixed points at $T = 0$. In the crossover region, a continuous variation in the effective critical exponent $4/9 \lesssim \eta^{\rm eff} \lesssim 1/2$ is observed. The phase diagram by changing the ratio $J_{\rm A}/J_{\rm B}$ is also investigated.


title: On the effect of wear on asperity heigth distributions, and the corresponding effect in the mechanical response

abstract: Since the time of the original Greenwood & Williamson paper, it was noticed that abrasion and wear lead to possibly bimodal distribution of asperity height distribution, with the upper tail of asperities following from the characteristics of the process. Using a limit case solution due to Borucki for the wear of an originally Gaussian distribution, it is shown here that the tail is indeed always Gaussian, but with different equivalent parameters. Therefore, if the wear process is light, one obtains a bimodal distribution and both may affect the resulting contact mechanics behaviour. In this short note, we illustrate just the main features of the problem. We conclude that it is an oversimplification to consider surfaces Gaussian.


title: Giant spin Nernst effect induced by resonant scattering at surfaces of metallic films

abstract: A new concept realizing giant spin Nernst effect in nonmagnetic metallic films is introduced. It is based on the idea of engineering an asymmetric energy dependence of the longitudinal and transverse electrical conductivities, as well as a pronounced energy dependence of the spin Hall angle in the vicinity of the Fermi level by the resonant impurity states at the Fermi level. We employ an analytical model and demonstrate the emergence of a giant spin Nernst effect in Ag(111) films using {\it ab-initio} calculations combined with the Boltzmann approach for transport properties arising from skew scattering off impurities.


title: $L^p\ (p>1)$ solutions of BSDEs with generators satisfying some non-uniform conditions in $t$ and $\omega$

abstract: This paper is devoted to the $L^p$ ($p>1$) solutions of one-dimensional backward stochastic differential equations (BSDEs for short) with general time intervals and generators satisfying some non-uniform conditions in $t$ and $\omega$. An existence and uniqueness result, a comparison theorem and an existence result for the minimal solutions are respectively obtained, which considerably improve some known works. Some classical techniques used to deal with the existence and uniqueness of $L^p$ ($p>1$) solutions of BSDEs with Lipschitz or linear-growth generators are also developed in this paper.


title: Event Search and Analytics: Detecting Events in Semantically Annotated Corpora for Search and Analytics

abstract: In this article, I present the questions that I seek to answer in my PhD research. I posit to analyze natural language text with the help of semantic annotations and mine important events for navigating large text corpora. Semantic annotations such as named entities, geographic locations, and temporal expressions can help us mine events from the given corpora. These events thus provide us with useful means to discover the locked knowledge in them. I pose three problems that can help unlock this knowledge vault in semantically annotated text corpora: i. identifying important events; ii. semantic search; and iii. event analytics.


title: \textit{Ab initio} study of reactive collisions between Rb($^{2}S$) or Rb($^{2}P$) and OH$^{-}$($^{1}\Sigma^{+}$)

abstract: A theoretical rate constant for the associative detachment reaction Rb($^{2}S$)$+$OH$^{-}$($^{1}\Sigma^{+}$)$\rightarrow \,$ RbOH($^{1}\Sigma^{+}$)$+\,e^{-}$ of 4$\,\times$10$^{-10}$cm$^{3}$s$^{-1}$ at 300 K has been calculated. This result agrees with the experimental rate constant of 2$^{+2}_{-1}\,\times10^{-10}$cm$^{3}$s$^{-1}$ obtained by Deiglmayr \textit{et al.} (Phys. Rev. A 86, 2012) for a temperature between 200 K and 600 K. A Langevin-based dynamics which depends on the crossing point between the anion (RbOH$^{-}$) and neutral (RbOH) potential energy surfaces has been used. The calculation were performed using the ECP28MDF effective core potential to describe the rubidium atom at the CCSD(T) level of theory and extended basis sets. The effect of ECPs and basis set on the height of the crossing point, and hence the rate constant, has been investigated. The temperature dependence of the latter is also discussed. Preliminary work on the potential energy surface for the excited reaction channel Rb($^{2}P$)+OH$^{-}$($^{1}\Sigma^{+}$) calculated at the CASSF-icMRCI level of theory is shown. We qualitatively discuss the charge transfer and associative detachment reactions arising from this excited entrance channel.


title: Selecting M-giants with infra-red photometry: Distances, metallicities and the Sagittarius stream

abstract: Using a spectroscopically confirmed sample of M-giants, M-dwarfs and quasars from the LAMOST survey, we assess how well WISE $\&$ 2MASS color-cuts can be used to select M-giant stars. The WISE bands are very efficient at separating M-giants from M-dwarfs and we present a simple classification that can produce a clean and relatively complete sample of M-giants. We derive a new photometric relation to estimate the metallicity for M-giants, calibrated using data from the APOGEE survey. We find a strong correlation between the $(W1-W2)$ color and $\rm [M/H]$, where almost all of the scatter is due to photometric uncertainties. We show that previous photometric distance relations, which are mostly based on stellar models, may be biased and devise a new empirical distance relation, investigating trends with metallicity and star formation history. Given these relations, we investigate the properties of M-giants in the Sagittarius stream. The offset in the orbital plane between the leading and trailing tails is reproduced and, by identifying distant M-giants in the direction of the Galactic anti-center, we confirm that the previously detected debris in the outer halo is the apocenter of the trailing tail. We also find tentative evidence supporting an existing overdensity near the leading tail in the Northern Galactic hemisphere, possibly an extension to the trailing tail (so-called Branch C). We have measured the metallicity distribution along the stream, finding a clear metallicity offset between the leading and trailing tails, in agreement with models for the stream formation. We include an online table of M-giants to facilitate further studies.


title: Calibration of Radio Interferometers Using a Sparse DoA Estimation Framework

abstract: The calibration of modern radio interferometers is a significant challenge, specifically at low frequencies. In this perspective, we propose a novel iterative calibration algorithm, which employs the popular sparse representation framework, in the regime where the propagation conditions shift dissimilarly the directions of the sources. More precisely, our algorithm is designed to estimate the apparent directions of the calibration sources, their powers, the directional and undirectional complex gains of the array elements and their noise powers, with a reasonable computational complexity. Numerical simulations reveal that the proposed scheme is statistically efficient at low SNR and even with additional non-calibration sources at unknown directions.


title: Improved 2D Intelligent Driver Model simulating synchronized flow and evolution concavity in traffic flow

abstract: This paper firstly show that 2 Dimensional Intelligent Driver Model (Jiang et al., PloS one, 9(4), e94351, 2014) is not able to replicate the synchronized traffic flow. Then we propose an improved model by considering the difference between the driving behaviors at high speeds and that at low speeds. Simulations show that the improved model can reproduce the phase transition from synchronized flow to wide moving jams, the spatiotemporal patterns of traffic flow induced by traffic bottleneck, and the evolution concavity of traffic oscillations (i.e. the standard deviation of the velocities of vehicles increases in a concave/linear way along the platoon). Validating results show that the empirical time series of traffic speed obtained from Floating Car Data can be well simulated as well.


title: Spectral properties of unimodular lattice triangulations

abstract: Random unimodular lattice triangulations have been recently used as an embedded random graph model, which exhibit a crossover behaviour between an ordered, large-world and a disordered, small-world behaviour. Using the ergodic Pachner flips that transform such triangulations into another and an energy functional that corresponds to the degree distribution variance, Markov chain Monte-Carlo simulations can be applied to study these graphs. Here, we consider the spectra of the adja cency and the Laplacian matrix as well as the algebraic connectivity and the spectral radius. Power law dependencies on the system size can clearly be identified and compared to analytical solutions for periodic ground states. For random triangulations we find a qualitative agreement of the spectral properties with well-known random graph models. In the microcanonical ensemble analytical approximations agree with numerical simulations. In the canonical ensemble a crossover behaviour can be found for the algebraic connectivity and the spectral radius, thus combining large-world and small-world behavior in one model. The considered spectral properties can be applied to transport problems on triangulation graphs and the crossover behaviour allows a tuning of important transport quantities.


title: Can we close the Bohr-Einstein quantum debate?

abstract: Recent experiments allowed concluding that Bell-type inequalities are indeed violated thus it is important to understand what it means and how can we explain the existence of strong correlations between outcomes of distant measurements. Do we have to announce that: Einstein was wrong, Nature is nonlocal and nonlocal correlations are produced due to the quantum magic and emerge, somehow, from outside space time? Fortunately such conclusions are unfounded because if supplementary parameters describing measuring instruments are correctly incorporated in a theoretical model then Bell-type inequalities may not be proven .We construct a simple probabilistic model explaining these correlations in a locally causal way. In our model measurement outcomes are neither predetermined nor produced in irreducibly random way. We explain in detail why, contrary to the general belief; an introduction of setting dependent parameters does not restrict experimenters' freedom of choice. Since the violation of Bell-type inequalities does not allow concluding that Nature is nonlocal and that quantum theory is complete thus the Bohr-Einstein quantum debate may not be closed. The continuation of this debate is not only important for a better understanding of Nature but also for various practical applications of quantum phenomena.


title: Band-Passes and Long Virtual Knot Concordance

abstract: Every classical knot is band-pass equivalent to the unknot or the trefoil. The band-pass class of a knot is a concordance invariant. Every ribbon knot, for example, is band-pass equivalent to the unknot. Here we introduce the long virtual knot concordance group $\mathscr{VC}$. It is shown that for every concordance class $[K] \in \mathscr{VC}$, there is a $J \in [K]$ that is not band-pass equivalent to $K$ and an $L \in [K]$ that is not band-pass equivalent to either the long unknot or any long trefoil. This is accomplished by proving that $v_{2,1}+v_{2,2} \pmod 2$ is a band-pass invariant but not a concordance invariant of long virtual knots, where $v_{2,1}$ and $v_{2,2}$ generate the degree two Polyak group for long virtual knots.


title: Ergodic measures on compact metric spaces for isometric actions by inductively compact groups

abstract: We obtain a partial converse of Vershik's description of ergodic probability measures on a compact metric space with respect to an isometric action by an inductively compact group. This allows us to identify, in this setting, the set of ergodic probability measures with the set of weak limit points of orbital measures.


title: The Sobolev moment problem and Jordan dilations

abstract: Moment problems and orthogonal polynomials, both meant in a single real variable, belong to the oldest problems in Classical Analysis. They have been developing for over a century in two parallel, mostly independent streams. During the last 20 years a rapid advancement of polynomials orthogonal in Sobolev space has been noticed; their moment counterpart seems to be not paid enough attention as it deserves. In this paper we intend to resume the theme and open the door for further, deeper study of the problem.


title: On the origins of polarization holes in Bok globules

abstract: Context. Polarimetric observations of Bok globules frequently show a decrease in the degree of polarization towards their central dense regions (polarization holes). This behaviour is usually explained with increased disalignment owing to high density and temperature, or insufficient angular resolution of a possibly complex magnetic field structure. Aims. We investigate whether a significant decrease in polarized emission of dense regions in Bok globules is possible under certain physical conditions. For instance, we evaluate the impact of optical depth effects and various properties of the dust phase. Methods. We use radiative transfer modelling to calculate the temperature structure of an analytical Bok globule model and simulate the polarized thermal emission of elongated dust grains. For the alignment of the dust grains, we consider a magnetic field and include radiative torque and internal alignment. Results. Besides the usual explanations, selected conditions of the temperature and density distribution, the dust phase and the magnetic field are also able to significantly decrease the polarized emission of dense regions in Bok globules. Taking submm/mm grains and typical column densities of existing Bok globules into consideration, the optical depth is high enough to decrease the degree of polarization by up to {\Delta}P~10%. If limited to the densest regions, dust grain growth to submm/mm size and accumulated graphite grains decrease the degree of polarization by up to {\Delta}P~10% and {\Delta}P~5%, respectively. However, the effect of the graphite grains occurs only if they do not align with the magnetic field.


title: A canonical rate-independent model of geometrically linear isotropic gradient plasticity with isotropic hardening and plastic spin accounting for the Burgers vector

abstract: In this paper we propose a canonical variational framework for rate-independent phenomenological geometrically linear gradient plasticity with plastic spin. The model combines the additive decomposition of the total distortion into non-symmetric elastic and plastic distortions, with a defect energy contribution taking account of the Burgers vector through a dependence only on the dislocation density tensor Curl(p) giving rise to a non-symmetric nonlocal backstress, and isotropic hardening response only depending on the accumulated equivalent plastic strain. The model is fully isotropic and satisfies linearized gauge-invariance conditions, i.e., only true state-variables appear. The model satisfies also the principle of maximum dissipation which allows to show existence for the weak formulation. For this result, a recently introduced Korn's inequality for incompatible tensor fields is necessary. Uniqueness is shown in the class of strong solutions. For vanishing energetic length scale, the model reduces to classical elasto-plasticity with symmetric plastic strain sym(p) and standard isotropic hardening.


title: Stochastic systems with memory and jumps

abstract: Stochastic systems with memory naturally appear in life science, economy, and finance. We take the modelling point of view of stochastic functional delay equations and we study these structures when the driving noises admit jumps. Our results concern existence and uniqueness of strong solutions, estimates for the moments and the fundamental tools of calculus, such as the It\^o formula. We study the robustness of the solution to the change of noises. Specifically, we consider the noises with infinite activity jumps versus an adequately corrected Gaussian noise. The study is presented in two different frameworks: we work with random variables in infinite dimensions, where the values are considered either in an appropriate $L^p$-type space or in the space of c\`adl\`ag paths. The choice of the value space is crucial from the modelling point of view as the different settings allow for the treatment of different models of memory or delay. Our techniques involve tools of infinite dimensional calculus and the stochastic calculus via regularization.


title: Component Based Modeling of Ultrasound Signals

abstract: This work proposes a component based model for the raw ultrasound signals acquired by the transducer elements. Based on this approach, before undergoing the standard digital processing chain, every sampled raw signal is first decomposed into a smooth background signal and a strong reflectors component. The decomposition allows for a suited processing scheme to be adjusted for each component individually. We demonstrate the potential benefit of this approach in image enhancement by suppressing side lobe artifacts, and in improvement of digital data compression. Applying our proposed processing schemes to real cardiac ultrasound data, we show that by separating the two components and compressing them individually, over twenty-fold reduction of the data size is achieved while retaining the image contents.


title: WKB approach to pair creation in spacetime-dependent fields

abstract: Besides tunneling in static potential landscapes, for example, the Wentzel-Kramers-Brillouin (WKB) approach is a powerful nonperturbative approximation tool to study particle creation due to time-dependent background fields, such as cosmological particle production or the Sauter-Schwinger effect, i.e., electron-positron pair creation in a strong electric field. However, our understanding of particle creation processes in background fields depending on both space and time is rather incomplete. In order to venture into this direction, we propose a generalization of the WKB method to truly spacetime-dependent fields and apply it to the case of a spacetime-dependent mass.


title: Gland Segmentation in Colon Histology Images: The GlaS Challenge Contest

abstract: Colorectal adenocarcinoma originating in intestinal glandular structures is the most common form of colon cancer. In clinical practice, the morphology of intestinal glands, including architectural appearance and glandular formation, is used by pathologists to inform prognosis and plan the treatment of individual patients. However, achieving good inter-observer as well as intra-observer reproducibility of cancer grading is still a major challenge in modern pathology. An automated approach which quantifies the morphology of glands is a solution to the problem. This paper provides an overview to the Gland Segmentation in Colon Histology Images Challenge Contest (GlaS) held at MICCAI'2015. Details of the challenge, including organization, dataset and evaluation criteria, are presented, along with the method descriptions and evaluation results from the top performing methods.


title: Secondary fast reconnecting instability in the sawtooth crash

abstract: In this work we consider magnetic reconnection in thin current sheets with both resistive and electron inertia effects. When the current sheet is produced by a primary instability of the internal kink type, the analysis of secondary instabilities indicates that reconnection proceeds on a time scale much shorter than the primary instability characteristic time. In the case of a sawtooth crash, non-collisional physics becomes important above a value of the Lundquist number which scales like S ~ (R/d_e)^{12/5}, in terms of the tokamak major radius R and of the electron skin depth d_e. This value is commonly achieved in present day devices. As collisionality is further reduced, the characteristic rate increases, approaching Alfv\'enic values when the primary instability approaches the collisionless regime.


title: Exciton Radiative Lifetime in Transition Metal Dichalcogenide Monolayers

abstract: We have investigated the exciton dynamics in transition metal dichalcogenide mono-layers using time-resolved photoluminescence experiments performed with optimized time-resolution. For MoSe2 monolayers, we measure $\tau_{rad}=1.8\pm0.2$ ps that we interpret as the intrinsic radiative recombination time. Similar values are found for WSe2 mono-layers. Our detailed analysis suggests the following scenario: at low temperature (T $\leq$ 50 K), the exciton oscillator strength is so large that the entire light can be emitted before the time required for the establishment of a thermalized exciton distribution. For higher lattice temperatures, the photoluminescence dynamics is characterized by two regimes with very different characteristic times. First the PL intensity drops drastically with a decay time in the range of the picosecond driven by the escape of excitons from the radiative window due to exciton- phonon interactions. Following this first non-thermal regime, a thermalized exciton population is established gradually yielding longer photoluminescence decay times in the nanosecond range. Both the exciton effective radiative recombination and non-radiative recombination channels including exciton-exciton annihilation control the latter. Finally the temperature dependence of the measured exciton and trion dynamics indicates that the two populations are not in thermodynamical equilibrium.


title: Methods used by WHO to estimate the Global burden of TB disease

abstract: This paper describes methodological details used by WHO in 2015 to estimate TB incidence, prevalence and mortality. Incidence and mortality are disaggregated by HIV status, age and sex. Methods to derive MDR-TB burden indicators are detailed. Four main methods were used to derive incidence: (i) case notification data combined with expert opinion about case detection gaps (120 countries representing 51% of global incidence); (ii) results from national TB prevalence surveys (19 countries, 46% of global incidence); (iii) notifications in high-income countries adjusted by a standard factor to account for under-reporting and underdiagnosis (73 countries, 3% of global incidence) and (iv) capture-recapture modelling (5 countries, 0.5% of global incidence). Prevalence was obtained from results of national prevalence surveys in 21 countries, representing 69% of global prevalence). In other countries, prevalence was estimated from incidence and disease duration. Mortality was obtained from national vital registration systems of mortality surveys in 129 countries (43% of global HIV-negative TB mortality). In other countries, mortality was derived indirectly from incidence and case fatality ratio.


title: Fast iterative method with a second order implicit difference scheme for time-space fractional convection-diffusion equations

abstract: In this paper we want to propose practical numerical methods to solve a class of initial-boundary problem of time-space fractional convection-diffusion equations (TSFCDEs). To start with, an implicit difference method based on two-sided weighted shifted Gr\"{u}nwald formulae is proposed with a discussion of the stability and convergence. We construct an implicit difference scheme (IDS) and show that it converges with second order accuracy in both time and space. Then, we develop fast solution methods for handling the resulting system of linear equation with the Toeplitz matrix. The fast Krylov subspace solvers with suitable circulant preconditioners are designed to deal with the resulting Toeplitz linear systems. Each time level of these methods reduces the memory requirement of the proposed implicit difference scheme from $\mathcal{O}(N^2)$ to $\mathcal{O}(N)$ and the computational complexity from $O(N^3)$ to $O(N\log N)$ in each iterative step, where $N$ is the number of grid nodes. Extensive numerical example runs show the utility of these methods over the traditional direct solvers of the implicit difference methods, in terms of computational cost and memory requirements.


title: Smoothing and non-smoothing via a flow tangent to the Ricci flow

abstract: We study a transformation of metric measure spaces introduced by Gigli and Mantegazza consisting in replacing the original distance with the length distance induced by the transport distance between heat kernel measures. We study the smoothing effect of this procedure in two important examples. Firstly, we show that in the case of some Euclidean cones, a singularity persists at the apex. Secondly, we generalize the construction to a sub-Riemannian manifold, namely the Heisenberg group, and show that it regularizes the space instantaneously to a smooth Riemannian manifold.


title: Error analysis of linear optics measurements via turn-by-turn beam position data in circular accelerators

abstract: Many advanced techniques have been developed, tested and implemented in the last decades in almost all circular accelerators across the world to measure the linear optics. However, the greater availability and accuracy of beam diagnostics and the ever better correction of linear magnetic lattice imperfections (beta beating at 1% level and coupling at 0.1%) are reaching what seems to be the intrinsic accuracy and precision of different measurement techniques. This paper aims to highlight and quantify, when possible, the limitations of one standard method, the harmonic analysis of turn-by-turn beam position data. To this end, new analytic formulas for the evaluation of lattice parameters modified by focusing errors are derived. The unexpected conclusion of this study is that for the ESRF storage ring (and possibly for any third generation light source operating at ultra-low coupling and with similar diagnostics), measurement and correction of linear optics via orbit beam position data are to be preferred to the analysis of turn-by-turn data.


title: Cacti with Extremal PI Index

abstract: The vertex PI index $PI(G) = \sum_{xy \in E(G)} [n_{xy}(x) + n_{xy}(y)]$ is a distance-based molecular structure descriptor, where $n_{xy}(x)$ denotes the number of vertices which are closer to the vertex $x$ than to the vertex $y$ and which has been the considerable research in computational chemistry dating back to Harold Wiener in 1947. A connected graph is a cactus if any two of its cycles have at most one common vertex. In this paper, we completely determine the extremal graphs with the largest and smallest vertex PI indices among all the cacti. As a consequence, we obtain the sharp bounds with corresponding extremal cacti and extend a known result.


title: Visualizing "Fermi arcs" in the Weyl semimetal TaAs

abstract: One of the hallmarks of Weyl semi-metals is the existence of unusual topological surface states known as 'Fermi arcs' [1-3]. The formation of these states is guaranteed by the existence of bulk Weyl points with opposite chirality. Tantalum Arsenide (TaAs) [4-9], a member of the newly discovered family of Weyl semi-metals [4,5], harbors a host of non-topological ('trivial') surface states overlapping in energy with the predicted 12 'Fermi arcs'. This overlap poses a major challenge in identifying the signatures of the arcs [10]. Here we harness the inherently distinct spatial structure of trivial and Fermi arc states to visualize the Fermi arcs for the first time using scanning tunneling microscopy. We do so in four distinct ways, each of which highlights a different aspect of their unusual nature. We reveal their relatively isotropic scattering signature, their energy dispersion and its relation to the bulk Weyl points, their deep bulk penetration relative to that of non-topological surface states and their weak coupling to the atomic structure. The latter is obtained by accounting for the spatial structure of the Bloch wavefunction and its effect on the scattering properties of the electrons off lattice defects in general. It thus provides a novel analysis tool for the spectroscopic characterization of electronic wavefunctions using scanning tunneling microscopy.


title: Dual Smoothing and Level Set Techniques for Variational Matrix Decomposition

abstract: We focus on the robust principal component analysis (RPCA) problem, and review a range of old and new convex formulations for the problem and its variants. We then review dual smoothing and level set techniques in convex optimization, present several novel theoretical results, and apply the techniques on the RPCA problem. In the final sections, we show a range of numerical experiments for simulated and real-world problems.


title: Kernel-based Tests for Joint Independence

abstract: We investigate the problem of testing whether $d$ random variables, which may or may not be continuous, are jointly (or mutually) independent. Our method builds on ideas of the two variable Hilbert-Schmidt independence criterion (HSIC) but allows for an arbitrary number of variables. We embed the $d$-dimensional joint distribution and the product of the marginals into a reproducing kernel Hilbert space and define the $d$-variable Hilbert-Schmidt independence criterion (dHSIC) as the squared distance between the embeddings. In the population case, the value of dHSIC is zero if and only if the $d$ variables are jointly independent, as long as the kernel is characteristic. Based on an empirical estimate of dHSIC, we define three different non-parametric hypothesis tests: a permutation test, a bootstrap test and a test based on a Gamma approximation. We prove that the permutation test achieves the significance level and that the bootstrap test achieves pointwise asymptotic significance level as well as pointwise asymptotic consistency (i.e., it is able to detect any type of fixed dependence in the large sample limit). The Gamma approximation does not come with these guarantees; however, it is computationally very fast and for small $d$, it performs well in practice. Finally, we apply the test to a problem in causal discovery.


title: Redividing the Cake

abstract: The paper considers fair allocation of resources that are already allocated in an unfair way. This setting requires a careful balance between the fairness considerations and the rights of the present owners. The paper presents re-division algorithms that attain various trade-off points between fairness and ownership rights, in various settings differing in the geometric constraints on the allotments: (a) no geometric constraints; (b) connectivity -- the cake is a one-dimensional interval and each piece must be a contiguous interval; (c) rectangularity -- the cake is a two-dimensional rectangle or rectilinear polygon and the pieces should be rectangles; (d) convexity -- the cake is a two-dimensional convex polygon and the pieces should be convex. These re-division algorithms have implications on another problem: the price-of-fairness -- the loss of social welfare caused by fairness requirements. Each algorithm implies an upper bound on the price-of-fairness with the respective geometric constraints.


title: A $\gamma\gamma$ Collider for the 750 GeV Resonant State

abstract: Recent data collected by ATLAS and CMS at 13 TeV collision energy of the LHC indicate the existence of a new resonant state $\phi$ with a mass of 750 GeV decaying into two photons $\gamma\gamma$. The properties of $\phi$ should be studied further at the LHC and also future colliders. Since only $\phi \to \gamma\gamma$ decay channel has been measured, one of the best ways to extract more information about $\phi$ is to use a $\gamma\gamma$ collider to produce $\phi$ at the resonant energy. In this work we show how a $\gamma\gamma$ collider helps to verify the existence of $\phi$ and to provide some of the most important information about the properties of $\phi$, such as branching fractions of $\phi\to V_1V_2$. Here $V_i$ can be $\gamma$, $Z$, or $W^\pm$. We also show that by studying angular distributions of the final $\gamma$'s in $\gamma\gamma \to \phi \to \gamma\gamma$, one can obtain crucial information about whether this state is a spin-0 or a spin-2 state.


title: Effects of an impermeable wall in dissipative dynamics of saturated porous media

abstract: A phase transition model for porous media in consolidation is studied. The model is able to describe the phenomenon of fluid-segregation during the consolidation process, i.e., the coexistence of two phases differing from fluid content inside the porous medium under static load. Considering pure Darcy dissipation, the dynamics is described by a Cahn-Hilliard-like system of partial differential equations (PDE). The goal, here, is to study the dynamics of the formation of stationary fluid-rich bubbles. The evolution of the strain and fluid density profiles of the porous medium is analyzed in two physical situation: fluid free to flow through the boundaries of the medium and fluid flow prevented at one of the two boundaries. Morover, an analytic result on the position of the interface between the two phases is provided.


title: Boundary-Finite Element discretization of time dependent acoustic scattering by elastic obstacles with piezoelectric behavior

abstract: A coupled BEM/FEM formulation for the transient interaction between an acoustic field and a piezoelectric scatterer is proposed. The scattered part of the acoustic wave is represented in terms of retarded layer potentials while the elastic displacement and electric potential are treated variationally. This results in an integro-differential system. Well posedness of a general Galerkin semi-discretization in space of the problem is shown in the Laplace domain and translated into explicit stability bounds in the time domain. Trapezoidal-Rule and BDF2 Convolution Quadrature are used in combination with matching time stepping for time discretization. Second order convergence is proven for the BDF2-based method. Numerical experiments are provided for BDF2 and Trapezoidal Rule based time evolution.


title: Width of the exotic $X_b(5568)$ state through its strong decay to $B_s^{0} \pi^{+}$

abstract: The width of the newly observed exotic state $X_b(5568)$ is calculated via its dominant strong decay to $% B_s^{0} \pi^{+}$ using the QCD sum rule method on the light-cone in conjunction with the soft-meson approximation. To this end, the vertex $X_{b}B_{s}\pi$ is studied and the strong coupling $g_{X_{b}B_{s}\pi}$ is computed employing for $X_b(5568)$ state the interpolating diquark-antidiquark current of the $[su][\overline{b}\overline{d}]$ type. The obtained prediction for the decay width of $X_b(5568)$ is confronted and a nice agreement found with the experimental data of the D0 Collaboration.