title: A statistical interpretation of spectral embedding: the generalised random dot product graph

abstract: Spectral embedding is a procedure which can be used to obtain vector representations of the nodes of a graph. This paper proposes a generalisation of the latent position network model known as the random dot product graph, to allow interpretation of those vector representations as latent position estimates. The generalisation is needed to model heterophilic connectivity (e.g., `opposites attract') and to cope with negative eigenvalues more generally. We show that, whether the adjacency or normalised Laplacian matrix is used, spectral embedding produces uniformly consistent latent position estimates with asymptotically Gaussian error (up to identifiability). The standard and mixed membership stochastic block models are special cases in which the latent positions take only $K$ distinct vector values, representing communities, or live in the $(K-1)$-simplex with those vertices, respectively. Under the stochastic block model, our theory suggests spectral clustering using a Gaussian mixture model (rather than $K$-means) and, under mixed membership, fitting the minimum volume enclosing simplex, existing recommendations previously only supported under non-negative-definite assumptions. Empirical improvements in link prediction (over the random dot product graph), and the potential to uncover richer latent structure (than posited under the standard or mixed membership stochastic block models) are demonstrated in a cyber-security example.


title: Probability of Induced Emission in Atoms

abstract: Expression for the probability of induced emission of high-order harmonics is obtained in the region where the multiphoton approximation is applicable to the description of the ionization of an atom. The dependence of this probability on the main parameters of the pump wave and the atomic medium is established. Criteria for observing emission are formulated with the consideration of phase locking. The possibility of amplifying a UV probe wave aimed through the region where the atoms interact with the pump wave is considered.


title: On the nth record gap between primes in an arithmetic progression

abstract: Let $q>r\ge1$ be coprime integers. Let $R(n,q,r)$ be the $n$th record gap between primes in the arithmetic progression $r$, $r+q$, $r+2q,\ldots,$ and denote by $N_{q,r}(x)$ the number of such records observed below $x$. For $x\to\infty$, we heuristically argue that if the limit of $N_{q,r}(x)/\log x$ exists, then the limit is 2. We also conjecture that $R(n,q,r)=O_q(n^2)$. Numerical evidence supports the conjectural (a.s.) upper bound $$R(n,q,r)<\varphi(q)n^2+(n+2)q\log^2 q.$$ The median (over $r$) of $R(n,q,r)$ grows like a quadratic function of $n$; so do the mean and quartile points of $R(n,q,r)$. For fixed values of $q\gtrsim200$ and $n\approx10$, the distribution of $R(n,q,r)$ is skewed to the right and close to both Gumbel and lognormal distributions; however, the skewness appears to slowly decrease as $n$ increases. The existence of a limiting distribution of $R(n,q,r)$ is an open question.


title: Miquel dynamics for circle patterns

abstract: We study a new discrete-time dynamical system on circle patterns with the combinatorics of the square grid. This dynamics, called Miquel dynamics, relies on Miquel's six circles theorem. We provide a coordinatization of the appropriate space of circle patterns on which the dynamics acts and use it to derive local recurrence formulas. Isoradial circle patterns arise as periodic points of Miquel dynamics. Furthermore, we prove that certain signed sums of intersection angles are preserved by the dynamics. Finally, when the initial circle pattern is spatially biperiodic with a fundamental domain of size two by two, we show that the appropriately normalized motion of intersection points of circles takes place along an explicit quartic curve.


title: The Geometric Block Model

abstract: To capture the inherent geometric features of many community detection problems, we propose to use a new random graph model of communities that we call a Geometric Block Model. The geometric block model generalizes the random geometric graphs in the same way that the well-studied stochastic block model generalizes the Erdos-Renyi random graphs. It is also a natural extension of random community models inspired by the recent theoretical and practical advancement in community detection. While being a topic of fundamental theoretical interest, our main contribution is to show that many practical community structures are better explained by the geometric block model. We also show that a simple triangle-counting algorithm to detect communities in the geometric block model is near-optimal. Indeed, even in the regime where the average degree of the graph grows only logarithmically with the number of vertices (sparse-graph), we show that this algorithm performs extremely well, both theoretically and practically. In contrast, the triangle-counting algorithm is far from being optimum for the stochastic block model. We simulate our results on both real and synthetic datasets to show superior performance of both the new model as well as our algorithm.


title: Anosov diffeomorphisms of products I. Negative curvature and rational homology spheres

abstract: We show that various classes of products of manifolds do not support transitive Anosov diffeomorphisms. Exploiting the Ruelle-Sullivan cohomology class, we prove that the product of a negatively curved manifold with a rational homology sphere does not support transitive Anosov diffeomorphisms. We extend this result to products of finitely many negatively curved manifolds of dimensions at least three with a rational homology sphere that has vanishing simplicial volume. As an application of this study, we obtain new examples of manifolds that do not support transitive Anosov diffeomorphisms, including certain manifolds with non-trivial higher homotopy groups and certain products of aspherical manifolds.


title: Waterlike anomalies in a 2D core-softened potential

abstract: We investigate the structural, thermodynamic and dynamic behavior of of a two dimensional core-corona system using Langevin Dynamics simulations. The particles are modeled using a core-softened potential that exhibits waterlike anomalies for the 3D case. For quasi-2D systems we have observed previously a new region of structural anomaly. Now, our results show that a new region of structural, density and diffusion anomalies arises for the 2D system. Our findings indicates that while the anomalous region at lower densities observed is due the competition between the two length scales in the potential, the traditional mechanism, the higher densities anomalous region is related to changes in the particles conformation and a melting region.


title: Factorising Polynomials over Finite Fields

abstract: The aim of this paper is to show that there exists a deterministic algorithm that can be applied to compute the factors of a polynomial of degree 2, defined over a finite field, given certain conditions.


title: Regularization and Variable Selection with Copula Prior

abstract: In this work, we show that under specific choices of the copula, the lasso, elastic net, and $g$-prior are particular cases of `copula prior,' for regularization and variable selection method. We present `lasso with Gauss copula prior' and `lasso with t-copula prior.' The simulation study and real-world data for regression, classification, and large time-series data show that the `copula prior' often outperforms the lasso and elastic net while having a comparable sparsity of representation. Also, the copula prior encourages a grouping effect. The strongly correlated predictors tend to be in or out of the model collectively under the copula prior. The `copula prior' is a generic method, which can be used to define the new prior distribution. The application of copulas in modeling prior distribution for Bayesian methodology has not been explored much. We present the resampling-based optimization procedure to handle big data with copula prior.


title: Some variations on Ensembled Random Survival Forest with application to Cancer Research

abstract: In this paper we describe a novel implementation of adaboost for prediction of survival function. We take different variations of the algorithm and compare the algorithms based on system run time and root mean square error. Our construction includes right censoring data and competing risk data too. We take different data set to illustrate the performance of the algorithms.


title: Fluorescent Silicon Clusters and Nanoparticles

abstract: The fluorescence of silicon clusters is reviewed. Atomic clusters of silicon have been at the focus of research for several decades because of the relevance of size effects for material properties, the importance of silicon in electronics and the potential applications in bio-medicine. To date numerous examples of nanostructured forms of fluorescent silicon have been reported. This article introduces the principles and underlying concepts relevant for fluorescence of nanostructured silicon such as excitation, energy relaxation, radiative and non-radiative decay pathways and surface passivation. Experimental methods for the production of silicon clusters are presented. The geometric and electronic properties are reviewed and the implications for the ability to emit fluorescence are discussed. Free and pure silicon clusters produced in molecular beams appear to have properties that are unfavourable for light emission. However, when passivated or embedded in a suitable host, they may emit fluorescence. The current available data show that both quantum confinement and localised transitions, often at the surface, are responsible for fluorescence. By building silicon clusters atom by atom, and by embedding them in shells atom by atom, new insights into the microscopic origins of fluorescence from nanoscale silicon can be expected.


title: A superfield realization of the integrated vertex operator in an $AdS_5\times S^5$ background

abstract: The integrated massless vertex operator in an $AdS_5\times S^5$ background in the pure spinor formalism is constructed in terms of superfields.


title: Editorial (a preface to a Special Issue of journal Applied Sciences on "Guided-Wave Optics")

abstract: This is a preface introducing a Special Issue on the topic of "Guided-Wave Optics" published by Applied Sciences. The Special Issue is a collection of 20 original and review articles dealing with various fundamental and applied aspects of the topic.


title: Semi-Static and Sparse Variance-Optimal Hedging

abstract: We consider hedging of a contingent claim by a 'semi-static' strategy composed of a dynamic position in one asset and static (buy-and-hold) positions in other assets. We give general representations of the optimal strategy and the hedging error under the criterion of variance-optimality and provide tractable formulas using Fourier-integration in case of the Heston model. We also consider the problem of optimally selecting a sparse semi-static hedging strategy, i.e. a strategy which only uses a small subset of available hedging assets. The developed methods are illustrated in an extended numerical example where we compute a sparse semi-static hedge for a variance swap using European options as static hedging assets.


title: Remarks on entanglement and identical particles

abstract: We argue that in the case of identical particles the most natural identification of separability, that is of absence of non-classical correlations, is via the factorization of mean values of commuting observables. It thus follows that separability and entanglement depend both on the state and on the choice of observables and are not absolute notions. We compare this point of view with a recent novel approach to the entanglement of identical particles, which allows for the definition of an entanglement entropy from a suitably defined reduced particle density matrix, without the need of labeling the system constituents. We contrast this figure of merit with the aforementioned lack of an absolute notion of entanglement by considering few paradigmatic examples.


title: Ultra-broadband gradient-pitch Bragg-Berry mirrors

abstract: The realization of geometric phase optical device operating over a broad spectral range is usually confronted with intrinsic limitations depending of the physical process at play. Here we propose to use chiral nematic liquid crystal slabs with helical ordering that varies in three dimensions. Namely, gradient-pitch cholesterics endowed with in-plane space-variant angular positioning of the supramolecular helix. By doing so, we show that the recently introduced Bragg-Berry mirrors [Opt. Lett. {\bf 41}, 3972-3975 (2016)] can be endowed with ultra-broadband spectral range. Experimental demonstration is made in the case of ultra-broadband optical vortex generation in the visible domain. These results offer practical solution to the polychromatic management of the orbital angular momentum of light combining the circular Bragg reflection of chiral media with the Berry phase.


title: AISHELL-1: An Open-Source Mandarin Speech Corpus and A Speech Recognition Baseline

abstract: An open-source Mandarin speech corpus called AISHELL-1 is released. It is by far the largest corpus which is suitable for conducting the speech recognition research and building speech recognition systems for Mandarin. The recording procedure, including audio capturing devices and environments are presented in details. The preparation of the related resources, including transcriptions and lexicon are described. The corpus is released with a Kaldi recipe. Experimental results implies that the quality of audio recordings and transcriptions are promising.


title: Mixed eldfellite compounds \ce{Na(Fe_{1/2}M_{1/2})(SO4)2} (M = Mn, Co, Ni): A new family of high electrode potential cathodes for the sodium-ion battery

abstract: Natural abundance of sodium and its similar behavior to lithium triggered recent extensive studies of cost-effective sodium-ion batteries (SIBs) for large-scale energy storage systems. A challenge is to develop electrode materials with a high electrode potential, specific capacity and a good rate capability. In this work we propose mixed eldfellite compounds \ce{Na_x(Fe_{1/2}M_{1/2})(SO4)2} (M = Mn, Co, Ni) as a new family of high electrode potential cathodes of SIBs and present their material properties predicted by first-principles calculations. The structural optimizations show that these materials have significantly small volume expansion rates below 5\% upon Na insertion/desertion with negative Na binding energies. Through the electronic structure calculations, we find band insulating properties and hole (and/or electron) polaron hoping as a possible mechanism for the charge transfer. Especially we confirm the high electrode voltages over 4 V with reasonably high specific capacities. We also investigate the sodium ion mobility by estimating plausible diffusion pathways and calculating the corresponding activation barriers, demonstrating the reasonably fast migrations of sodium ions during the operation. Our calculation results indicate that these mixed eldfellite compounds can be suitable materials for high performance SIB cathodes.


title: A laser excitation scheme for $^{229\text{m}}$Th

abstract: Direct laser excitation of the lowest known nuclear excited state in $^{229}$Th has been a longstanding objective. It is generally assumed that reaching this goal would require a considerably reduced uncertainty of the isomer's excitation energy compared to the presently adopted value of $(7.8\pm 0.5)$ eV. Here we present a direct laser excitation scheme for $^{229\text{m}}$Th, which circumvents this requirement. The proposed excitation scheme makes use of already existing laser technology and therefore paves the way for nuclear laser spectroscopy. In this concept, the recently experimentally observed internal-conversion decay channel of the isomeric state is used for probing the isomeric population. A signal-to-background ratio of better than $10^4$ and a total measurement time of less than three days for laser scanning appear to be achievable.


title: Performance Evaluation of Spatial Complementary Code Keying Modulation in MIMO Systems

abstract: Spatial complementary code keying modulation (SCCKM) is proposed as a novel block coding modulation scheme. An input binary sequence is modulated based on the different lengths of complementary code keying (CCK) modulation and then spread across the transmit antennas (spatial domain) in a multiple input multiple output (MIMO) system exploiting orthogonal frequency division multiplexing (OFDM). At the receiver side, zero forcing equalization is applied to the OFDM modulated data to mitigate the effect of the multipath fast fading channel and then followed by maximum likelihood (ML) detection to retrieve the input sequence. The performance of SCCKM in different MIMO systems is compared to that of spatial modulation (SM) as a baseline scheme. Simulation results show that for the same spectral efficiency, SCCKM is able to substantially improve the bit error rate (BER).


title: Orientation-dependent deformation mechanisms of bcc niobium nanoparticles

abstract: Nanoparticles usually exhibit pronounced anisotropic properties, and a close insight into the atomic-scale deformation mechanisms is of great interest. In present study, atomic simulations are conducted to analyze the compression of bcc nanoparticles, and orientation-dependent features are addressed. It is revealed that surface morphology under indenter predominantly governs the initial elastic response. The loading curve follows the flat punch contact model in [110] compression, while it obeys the Hertzian contact model in [111] and [001] compressions. In plastic deformation regime, full dislocation gliding is dominated in [110] compression, while deformation twinning is prominent in [111] compression, and these two mechanisms coexist in [001] compression. Such deformation mechanisms are distinct from those in bulk crystals under nanoindentation and nanopillars under compression, and the major differences are also illuminated. Our results provide an atomic perspective on the mechanical behaviors of bcc nanoparticles and are helpful for the design of nanoparticle-based components and systems.


title: Semi-Static Variance-Optimal Hedging in Stochastic Volatility Models with Fourier Representation

abstract: In a financial market model, we consider the variance-optimal semi-static hedging of a given contingent claim, a generalization of the classic variance-optimal hedging. To obtain a tractable formula for the expected squared hedging error and the optimal hedging strategy, we use a Fourier approach in a general multidimensional semimartingale factor model. As a special case, we recover existing results for variance-optimal hedging in affine stochastic volatility models. We apply the theory to set up a variance-optimal semi-static hedging strategy for a variance swap in both the Heston and the 3/2-model, the latter of which is a non-affine stochastic volatility model.


title: Convergence Analysis and Design of Multi-block ADMM via Switched Control Theory

abstract: We consider three challenges in multi-block Alternating Direction Method of Multipliers (ADMM): building convergence conditions for ADMM with any block (variable) sequence, finding available block sequences to be fit for ADMM, and designing useful parameter controllers for ADMM with unfixed parameters. To address these challenges, we develop a switched control framework for studying multi-block ADMM. First, since ADMM recursively and alternately updates the block-variables, it is converted into a discrete-time switched dynamical system. Second, we study exponential stability and stabilizability of the switched system for linear convergence analysis and design of ADMM by employing switched Lyapunov functions. Moreover, linear matrix inequalities conditions are proposed to ensure convergence of ADMM under arbitrary sequence, to find convergent sequences, and to design the fixed parameters. These conditions are checked and solved by employing semidefinite programming. Numerical experiments further verify the effectiveness of our proposed theories.


title: Explicit Solution for Constrained Stochastic Linear-Quadratic Control with Multiplicative Noise

abstract: We study in this paper a class of constrained linear-quadratic (LQ) optimal control problem formulations for the scalar-state stochastic system with multiplicative noise, which has various applications, especially in the financial risk management. The linear constraint on both the control and state variables considered in our model destroys the elegant structure of the conventional LQ formulation and has blocked the derivation of an explicit control policy so far in the literature. We successfully derive in this paper the analytical control policy for such a class of problems by utilizing the state separation property induced from its structure. We reveal that the optimal control policy is a piece-wise affine function of the state and can be computed off-line efficiently by solving two coupled Riccati equations. Under some mild conditions, we also obtain the stationary control policy for infinite time horizon. We demonstrate the implementation of our method via some illustrative examples and show how to calibrate our model to solve dynamic constrained portfolio optimization problems.


title: Multiplicity of Solutions for Quasilinear Elliptic Problems

abstract: It is established existence, uniqueness and multiplicity of solutions for a quasilinear elliptic problem problems driven by $\Phi$-Laplacian operator. Here we consider the reflexive and nonreflexive cases using an auxiliary problem. In order to prove our main results we employ variational methods, regularity results and truncation techniques.


title: Information Geometry of Quantum Resources

abstract: I review recent works showing that information geometry is a useful framework to characterize quantum coherence and entanglement. Quantum systems exhibit peculiar properties which cannot be justified by classical physics, e.g. quantum coherence and quantum correlations. Once confined to thought experiments, they are nowadays created and manipulated by exerting an exquisite experimental control of atoms, molecules and photons. It is important to identify and quantify such quantum features, as they are deemed to be key resources to achieve supraclassical performances in computation and communication protocols. The information geometry viewpoint elucidates the advantage provided by quantum superpositions in phase estimation. Also, it enables to link measures of coherence and entanglement to observables, which can be evaluated in a laboratory by a limited number of measurements.


title: The Unknown Subgroup of $Aut(E_8)$

abstract: The $E_8$ lattice has been thoroughly studied for more than a century and nearly all the maximal subgroups of $W(E_8)$ have been described-all except $2A_9$. We will show that $2A_9$ has simple descriptions from three different perspectives: looking at $E_8/2E_8$; looking at the lattice's norm 2 vectors; and looking at its norm 4 vectors. Two of the three descriptions are especially simple: $2A_9$ stabilizes a partition of the 2160 norm 4 vectors into nine scale copies of $E_8$; and $2A_9$ stabilizes a partition of the 135 isotropic points of $E_8/2E_8$ into nine disjoint isotropic 4-spaces of 15 points each (so the well-known known isomorphism $A_8$ $\cong$ $L_4$(2) is visible within $E_8$.)


title: Topomap: Topological Mapping and Navigation Based on Visual SLAM Maps

abstract: Visual robot navigation within large-scale, semi-structured environments deals with various challenges such as computation intensive path planning algorithms or insufficient knowledge about traversable spaces. Moreover, many state-of-the-art navigation approaches only operate locally instead of gaining a more conceptual understanding of the planning objective. This limits the complexity of tasks a robot can accomplish and makes it harder to deal with uncertainties that are present in the context of real-time robotics applications. In this work, we present Topomap, a framework which simplifies the navigation task by providing a map to the robot which is tailored for path planning use. This novel approach transforms a sparse feature-based map from a visual Simultaneous Localization And Mapping (SLAM) system into a three-dimensional topological map. This is done in two steps. First, we extract occupancy information directly from the noisy sparse point cloud. Then, we create a set of convex free-space clusters, which are the vertices of the topological map. We show that this representation improves the efficiency of global planning, and we provide a complete derivation of our algorithm. Planning experiments on real world datasets demonstrate that we achieve similar performance as RRT* with significantly lower computation times and storage requirements. Finally, we test our algorithm on a mobile robotic platform to prove its advantages.


title: An alternative to continuous univariate distributions supported on a bounded interval: The BMT distribution

abstract: In this paper, we introduce the BMT distribution as an unimodal alternative to continuous univariate distributions supported on a bounded interval. The ideas behind the mathematical formulation of this new distribution come from computer aid geometric design, specifically from Bezier curves. First, we review general properties of a distribution given by parametric equations and extend the definition of a Bezier distribution. Then, after proposing the BMT cumulative distribution function, we derive its probability density function and a closed-form expression for quantile function, median, interquartile range, mode, and moments. The domain change from [0,1] to [c,d] is mentioned. Estimation of parameters is approached by the methods of maximum likelihood and maximum product of spacing. We test the numerical estimation procedures using some simulated data. Usefulness and flexibility of the new distribution are illustrated in three real data sets. The BMT distribution has a significant potential to estimate domain parameters and to model data outside the scope of the beta or similar distributions.


title: Towards a supercharacter theory of the parabolic subgroups

abstract: The supercharacter theory is constructed for the parabolic subgroups of $\mathrm{GL}(n,\Fq)$ with blocks of orders less or equal to two. The author formulated the hypotheses on construction of a supercharacter theory for an arbitrary parabolic subgroup in $\mathrm{GL}(n,\Fq)$.


title: Construction of full diversity $D_n$-lattices for all $n$

abstract: In this paper we construct some families of rotated $D_n$-lattices with full diversity for any $n$. These lattices can be good for signal transmission over both Gaussian and Rayleigh fading channels. In order to get bounds for their minimum product distances, we show that the $Z$-modules used in \cite{sethoggier} to obtain rotated $\mathbb{Z}^{n}$-lattices with $n$ odd are ideals and find a sufficient condition for such ideals being principal ideals.


title: A priori estimates for some elliptic equations involving the $p$-Laplacian

abstract: We consider the Dirichlet problem for positive solutions of the equation $-\Delta_p (u) = f(u)$ in a convex, bounded, smooth domain $\Omega \subset\R^N$, with $f$ locally Lipschitz continuous. \par We provide sufficient conditions guarantying $L^{\infty} $ a priori bounds for positive solutions of some elliptic equations involving the $p$-Laplacian and extend the class of known nonlinearities for which the solutions are $L^{\infty} $ a priori bounded. As a consequence we prove the existence of positive solutions in convex bounded domains.


title: DeepLung: 3D Deep Convolutional Nets for Automated Pulmonary Nodule Detection and Classification

abstract: In this work, we present a fully automated lung CT cancer diagnosis system, DeepLung. DeepLung contains two parts, nodule detection and classification. Considering the 3D nature of lung CT data, two 3D networks are designed for the nodule detection and classification respectively. Specifically, a 3D Faster R-CNN is designed for nodule detection with a U-net-like encoder-decoder structure to effectively learn nodule features. For nodule classification, gradient boosting machine (GBM) with 3D dual path network (DPN) features is proposed. The nodule classification subnetwork is validated on a public dataset from LIDC-IDRI, on which it achieves better performance than state-of-the-art approaches, and surpasses the average performance of four experienced doctors. For the DeepLung system, candidate nodules are detected first by the nodule detection subnetwork, and nodule diagnosis is conducted by the classification subnetwork. Extensive experimental results demonstrate the DeepLung is comparable to the experienced doctors both for the nodule-level and patient-level diagnosis on the LIDC-IDRI dataset.


title: On Isoperimetric Stability

abstract: We show that a non-empty subset of an abelian group with a small edge boundary must be large; in particular, if $A$ and $S$ are finite, non-empty subsets of an abelian group such that $S$ is independent, and the edge boundary of $A$ with respect to $S$ does not exceed $(1-\gamma)|S||A|$ with a real $\gamma\in(0,1]$, then $|A| \ge 4^{(1-1/d)\gamma |S|}$, where $d$ is the smallest order of an element of $S$. Here the constant $4$ is best possible. As a corollary, we derive an upper bound for the size of the largest independent subset of the set of popular differences of a finite subset of an abelian group. For groups of exponent $2$ and $3$, our bound translates into a sharp estimate for the additive dimension of the popular difference set. We also prove, as an auxiliary result, the following estimate of possible independent interest: if $A \subset \mathbb Z^n$ is a finite, non-empty downset then, denoting by $w(a)$ the number of non-zero components of the vector $a\in A$, we have \[\frac1{|A|} \sum_{a\in A} w(a) \le \frac12\, \log_2 |A|.\]


title: Primitive Element Pairs with One Prescribed Trace over a Finite Field

abstract: In this article, we establish a sufficient condition for the existence of a primitive element $\alpha \in {\mathbb{F}_{q^n}}$ such that the element $\alpha+\alpha^{-1}$ is also a primitive element of ${\mathbb{F}_{q^n}},$ and $Tr_{\mathbb{F}_{q^n}|\mathbb{F}_{q}}(\alpha)=a$ for any prescribed $a \in \mathbb{F}_q$, where $q=p^k$ for some prime $p$ and positive integer $k$. We prove that every finite field $\mathbb{F}_{q^n}~ (n \geq5),$ contains such primitive elements except for finitely many values of $q$ and $n$. Indeed, by computation, we conclude that there are no actual exceptional pairs $(q,n)$ for $n\geq5.$


title: First-order vortices in a gauged $CP(2)$ model with a Chern-Simons term

abstract: We consider a gauged $CP(2)$ theory in the presence of the Chern-Simons action, focusing our attention on those time-independent solutions possessing radial symmetry. In this context, we develop a coherent first-order framework via the Bogomol'nyi prescription, from which we obtain the corresponding energy lower-bound and the first-order equations the model supports. We use these expressions to introduce effective BPS scenarios, solving the resulting first-order equations by means of the finite-difference scheme, this way attaining genuine field solutions engendering topological configurations. We depict the new profiles, commenting on the main properties they engender.


title: CIRCE: The Canarias InfraRed Camera Experiment for the Gran Telescopio Canarias

abstract: The Canarias InfraRed Camera Experiment (CIRCE) is a near-infrared (1-2.5 micron) imager, polarimeter and low-resolution spectrograph operating as a visitor instrument for the Gran Telescopio Canarias 10.4-meter telescope. It was designed and built largely by graduate students and postdocs, with help from the UF astronomy engineering group, and is funded by the University of Florida and the U.S. National Science Foundation. CIRCE is intended to help fill the gap in near-infrared capabilities prior to the arrival of EMIR to the GTC, and will also provide the following scientific capabilities to compliment EMIR after its arrival: high-resolution imaging, narrowband imaging, high-time-resolution photometry, imaging polarimetry, low resolution spectroscopy. In this paper, we review the design, fabrication, integration, lab testing, and on-sky performance results for CIRCE. These include a novel approach to the opto-mechanical design, fabrication, and alignment.


title: Inclusive search for a highly boosted Higgs boson decaying to a bottom quark-antiquark pair

abstract: An inclusive search for the standard model Higgs boson ($\mathrm{H}$) produced with large transverse momentum ($p_\mathrm{T}$) and decaying to a bottom quark-antiquark pair ($\mathrm{b}\overline{\mathrm{b}}$) is performed using a data set of pp collisions at $\sqrt{s}=$ 13 TeV collected with the CMS experiment at the LHC. The data sample corresponds to an integrated luminosity of 35.9 fb$^{-1}$. A highly Lorentz-boosted Higgs boson decaying to $\mathrm{b}\overline{\mathrm{b}}$ is reconstructed as a single, large radius jet and is identified using jet substructure and dedicated $\mathrm{b}$ tagging techniques. The method is validated with $\mathrm{Z}\to\mathrm{b}\overline{\mathrm{b}}$ decays. The $\mathrm{Z}\to\mathrm{b}\overline{\mathrm{b}}$ process is observed for the first time in the single-jet topology with a local significance of 5.1 standard deviations (5.8 expected). For a Higgs boson mass of 125 GeV, an excess of events above the expected background is observed (expected) with a local significance of 1.5 (0.7) standard deviations. The measured cross section times branching fraction for production via gluon fusion of $\mathrm{H} \rightarrow \mathrm{b}\overline{\mathrm{b}}$ with $p_\mathrm{T} > $450 GeV and in the pseudorapidity range $-$2.5 $< \eta <$ 2.5 is 74 $\pm$ 48 (stat) $_{-10}^{+17}$ (syst) fb, which is consistent within uncertainties with the standard model prediction.


title: Existence result under flatness condition for a nonlinear elliptic equation with Sobolev exponent

abstract: In this paper, we consider the following nonlinear elliptic equation with Dirichlet boundary condition: $-\Delta u=K(x)u^{\frac{n+2}{n-2}},\, u>0$ in $\Omega,\, u=0$ on $\partial\Omega$, where $\Omega$ is a smooth bounded domain in $\mathbb{R}^n,$ $n\geqslant 4,$ and $K$ is a $\mathcal{C}^1$-positive function in $\bar{\Omega}$. Under the assumption that the order of flatness at each critical point of $K$ is $\beta \in ]\,n-2,\,n[,$ we give precise estimates on the looses of the compactness, and we prove an existence result through an Euler-Hopf type formula.


title: Generating Compact Tree Ensembles via Annealing

abstract: Tree ensembles are flexible predictive models that can capture relevant variables and to some extent their interactions in a compact and interpretable manner. Most algorithms for obtaining tree ensembles are based on versions of boosting or Random Forest. Previous work showed that boosting algorithms exhibit a cyclic behavior of selecting the same tree again and again due to the way the loss is optimized. At the same time, Random Forest is not based on loss optimization and obtains a more complex and less interpretable model. In this paper we present a novel method for obtaining compact tree ensembles by growing a large pool of trees in parallel with many independent boosting threads and then selecting a small subset and updating their leaf weights by loss optimization. We allow for the trees in the initial pool to have different depths which further helps with generalization. Experiments on real datasets show that the obtained model has usually a smaller loss than boosting, which is also reflected in a lower misclassification error on the test set.


title: A computational study of the mechanisms of growth-driven folding patterns on shells, with application to the developing brain

abstract: We consider the mechanisms by which folds, or sulci (troughs) and gyri (crests), develop in the brain. This feature, common to many gyrencephalic species including humans, has attracted recent attention from soft matter physicists. It occurs due to inhomogeneous, and predominantly tangential, growth of the cortex, which causes circumferential compression, leading to a bifurcation of the solution path to a folded configuration. The problem can be framed as one of buckling in the regime of linearized elasticity. However, the brain is a very soft solid, which is subject to large strains due to inhomogeneous growth. As a consequence, the morphomechanics of the developing brain demonstrates an extensive post-bifurcation regime. Nonlinear elasticity studies of growth-driven brain folding have established the conditions necessary for the onset of folding, and for its progression to configurations broadly resembling gyrencephalic brains. The reference, unfolded, configurations in these treatments have a high degree of symmetry--typically, ellipsoidal. Depending on the boundary conditions, the folded configurations have symmetric or anti-symmetric patterns. However, these configurations do not approximate the actual morphology of, e.g., human brains, which display unsymmetric folding. More importantly, from a neurodevelopmental standpoint, many of the unsymmetric sulci and gyri are notably robust in their locations. Here, we initiate studies on the physical mechanisms and geometry that control the development of primary sulci and gyri. In this preliminary communication we carry out computations with idealized geometries, boundary conditions and parameters, seeking a pattern resembling one of the first folds to form: the Central Sulcus.


title: Method for Mode Mixing Separation in Empirical Mode Decomposition

abstract: The Empirical Mode Decomposition (EMD) is a signal analysis method that separates multi-component signals into single oscillatory modes called intrinsic mode functions (IMFs), each of which can generally be associated to a physical meaning of the process from which the signal is obtained. When the phenomena of mode mixing occur, as a result of the EMD sifting process, the IMFs can lose their physical meaning hindering the interpretation of the results of the analysis. In the paper, "One or Two frequencies? The Empirical Mode Decomposition Answers", Gabriel Rilling and Patrick Flandrin [3] presented a rigorous mathematical analysis that explains how EMD behaves in the case of a composite two-tones signal and the amplitude and frequency ratios by which EMD will perform a good separation of tones. However, the authors did not propose a solution for separating the neighboring tones that will naturally remain mixed after an EMD. In this paper, based on the findings by Rilling and Flandrin, a method that can separate neighbouring spectral components, that will naturally remain within a single IMF, is presented. This method is based on reversing the conditions by which mode mixing occurs and that were presented in the map by Rilling and Flandrin in the above mentioned paper. Numerical experiments with signals containing closely spaced spectral components shows the effective separation of modes that EMD can perform after this principle is applied. The results verify also the regimes presented in the theoretical analysis by Rilling and Flandrin.


title: Forecasting of commercial sales with large scale Gaussian Processes

abstract: This paper argues that there has not been enough discussion in the field of applications of Gaussian Process for the fast moving consumer goods industry. Yet, this technique can be important as it e.g., can provide automatic feature relevance determination and the posterior mean can unlock insights on the data. Significant challenges are the large size and high dimensionality of commercial data at a point of sale. The study reviews approaches in the Gaussian Processes modeling for large data sets, evaluates their performance on commercial sales and shows value of this type of models as a decision-making tool for management.


title: The Main Belt Comets and Ice in the Solar System

abstract: We review the evidence for buried ice in the asteroid belt; specifically the questions around the so-called Main Belt Comets (MBCs). We summarise the evidence for water throughout the Solar System, and describe the various methods for detecting it, including remote sensing from ultraviolet to radio wavelengths. We review progress in the first decade of study of MBCs, including observations, modelling of ice survival, and discussion on their origins. We then look at which methods will likely be most effective for further progress, including the key challenge of direct detection of (escaping) water in these bodies.


title: Current reversals and metastable states in the infinite Bose-Hubbard chain with local particle loss

abstract: We present an algorithm which combines the quantum trajectory approach to open quantum systems with a density-matrix renormalization group scheme for infinite one-dimensional lattice systems. We apply this method to investigate the long-time dynamics in the Bose-Hubbard model with local particle loss starting from a Mott-insulating initial state with one boson per site. While the short-time dynamics can be described even quantitatively by an equation of motion (EOM) approach at the mean-field level, many-body interactions lead to unexpected effects at intermediate and long times: local particle currents far away from the dissipative site start to reverse direction ultimately leading to a metastable state with a total particle current pointing away from the lossy site. An alternative EOM approach based on an effective fermion model shows that the reversal of currents can be understood qualitatively by the creation of holon-doublon pairs at the edge of the region of reduced particle density. The doublons are then able to escape while the holes move towards the dissipative site, a process reminiscent---in a loose sense---of Hawking radiation.


title: Applying Machine Learning Methods to Enhance the Distribution of Social Services in Mexico

abstract: The Government of Mexico's social development agency, SEDESOL, is responsible for the administration of social services and has the mission of lifting Mexican families out of poverty. One key challenge they face is matching people who have social service needs with the services SEDESOL can provide accurately and efficiently. In this work we describe two specific applications implemented in collaboration with SEDESOL to enhance their distribution of social services. The first problem relates to systematic underreporting on applications for social services, which makes it difficult to identify where to prioritize outreach. Responding that five people reside in a home when only three do is a type of underreporting that could occur while a social worker conducts a home survey with a family to determine their eligibility for services. The second involves approximating multidimensional poverty profiles across households. That is, can we characterize different types of vulnerabilities -- for example, food insecurity and lack of health services -- faced by those in poverty? We detail the problem context, available data, our machine learning formulation, experimental results, and effective feature sets. As far as we are aware this is the first time government data of this scale has been used to combat poverty within Mexico. We found that survey data alone can suggest potential underreporting. Further, we found geographic features useful for housing and service related indicators and transactional data informative for other dimensions of poverty. The results from our machine learning system for estimating poverty profiles will directly help better match 7.4 million individuals to social programs.


title: Multivariate Gaussian Network Structure Learning

abstract: We consider a graphical model where a multivariate normal vector is associated with each node of the underlying graph and estimate the graphical structure. We minimize a loss function obtained by regressing the vector at each node on those at the remaining ones under a group penalty. We show that the proposed estimator can be computed by a fast convex optimization algorithm. We show that as the sample size increases, the estimated regression coefficients and the correct graphical structure are correctly estimated with probability tending to one. By extensive simulations, we show the superiority of the proposed method over comparable procedures. We apply the technique on two real datasets. The first one is to identify gene and protein networks showing up in cancer cell lines, and the second one is to reveal the connections among different industries in the US.


title: Accurate Anisotropic Gaussian Type Orbital Basis Sets for Atoms in Strong Magnetic Fields

abstract: In high magnetic field calculations, anisotropic Gaussian type orbital (AGTO) basis functions are capable of reconciling the competing demands of the spherically symmetric Coulombic interaction and cylindrical magnetic ($B$ field) confinement. However, the best available {\it a priori} procedure for composing highly accurate AGTO sets for atoms in a strong $B$ field [Phys.\ Rev. A {\bf 90}, 022504 (2014)] yields very large basis sets. Their size is problematical for use in any calculation with unfavorable computational cost scaling. Here we provide an alternative constructive procedure. It is based upon analysis of the underlying physics of atoms in $B$ fields that allows identification of several principles for the construction of AGTO basis sets. Aided by numerical optimization and parameter fitting, followed by fine tuning of fitting parameters, we devise formulae for generating accurate AGTO basis sets in an arbitrary $B$ field. For the hydrogen iso-electronic sequence, a set depends on $B$ field strength, nuclear charge, and upon orbital quantum numbers. For multi-electron systems, the basis set formulae also include adjustment to account for orbital occupations. Tests of the new basis sets for atoms H through C ($1 \le Z \le 6$), and ions Li$^+$, Be$^+$, and B$^+$, in a wide $B$ field range ($0 \le B \le 2000$ a.u.), show an accuracy better than a few $\mu$H for single-electron systems, and a few hundredths to a few mHs for multi-electron atoms. The relative errors are similar for different atoms and ions in a large $B$ field range, from a few to a couple of tens of millionths, thereby confirming rather uniform accuracy across the nuclear charge $Z$ and $B$ field strength values. Residual basis set errors are two to three orders of magnitude smaller than the electronic correlation energies in muti-electron atoms ...


title: Deep Automated Multi-task Learning

abstract: Multi-task learning (MTL) has recently contributed to learning better representations in service of various NLP tasks. MTL aims at improving the performance of a primary task, by jointly training on a secondary task. This paper introduces automated tasks, which exploit the sequential nature of the input data, as secondary tasks in an MTL model. We explore next word prediction, next character prediction, and missing word completion as potential automated tasks. Our results show that training on a primary task in parallel with a secondary automated task improves both the convergence speed and accuracy for the primary task. We suggest two methods for augmenting an existing network with automated tasks and establish better performance in topic prediction, sentiment analysis, and hashtag recommendation. Finally, we show that the MTL models can perform well on datasets that are small and colloquial by nature.


title: Moduli of semistable sheaves as quiver moduli

abstract: In the 1980s Dr\'ezet and Le Potier realized moduli spaces of Gieseker-semistable sheaves on $\mathbb{P}^2$ as what are now called quiver moduli spaces. We discuss how this construction can be understood using t-structures and exceptional collections on derived categories, and how it can be extended to a similar result on $\mathbb{P}^1\times\mathbb{P}^1$.


title: Anisotropic coalescence criterion for nanoporous materials

abstract: Ductile fracture through void growth and coalescence depends significantly on the plastic anisotropy of the material and on void size, as shown by experiments and/or numerical simulations through several studies. Macroscopic (homogenized) yield criteria aiming at modeling nanoporous materials have been proposed only for the growth regime, i.e. non-interacting voids. The aim of this study is thus to provide a yield criterion for nanoporous materials relevant for the coalescence regime, i.e. when plastic flow is localized between voids. Through homogenization and limit analysis, and accounting for interface stresses at the void-matrix interface, analytical coalescence criterion is derived under the following conditions: axisymmetric loading, orthotropic material obeying Hill's plasticity, cylindrical voids in cylindrical unit-cell. Incidentally, an orthotropic extension of the existing isotropic modeling of interface stresses through limit analysis is described and used. The proposed coalescence criterion is then extended to account for combined tension and shear loading conditions. Numerical limit analyses have been performed under specific conditions / materials parameters to get supposedly exact (up to numerical errors) results of coalescence stress. A good agreement between the analytical coalescence criteria derived in this study and numerical results is found for elongated spheroidal voids, making them usable to predict the onset of void coalescence in ductile fracture modeling of nanoporous materials.


title: Speech Dereverberation Using Nonnegative Convolutive Transfer Function and Spectro temporal Modeling

abstract: This paper presents two single channel speech dereverberation methods to enhance the quality of speech signals that have been recorded in an enclosed space. For both methods, the room acoustics are modeled using a nonnegative approximation of the convolutive transfer function (NCTF), and to additionally exploit the spectral properties of the speech signal, such as the low rank nature of the speech spectrogram, the speech spectrogram is modeled using nonnegative matrix factorization (NMF). Two methods are described to combine the NCTF and NMF models. In the first method, referred to as the integrated method, a cost function is constructed by directly integrating the speech NMF model into the NCTF model, while in the second method, referred to as the weighted method, the NCTF and NMF based cost functions are weighted and summed. Efficient update rules are derived to solve both optimization problems. In addition, an extension of the integrated method is presented, which exploits the temporal dependencies of the speech signal. Several experiments are performed on reverberant speech signals with and without background noise, where the integrated method yields a considerably higher speech quality than the baseline NCTF method and a state of the art spectral enhancement method. Moreover, the experimental results indicate that the weighted method can even lead to a better performance in terms of instrumental quality measures, but that the optimal weighting parameter depends on the room acoustics and the utilized NMF model. Modeling the temporal dependencies in the integrated method was found to be useful only for highly reverberant conditions.


title: Quantum phase transition of chiral Majorana fermion in the presence of disorder

abstract: We study the quantum phase transitions of a disordered two-dimensional quantum anomalous Hall insulator with $s$-wave superconducting proximity, which are governed by the percolation theory of chiral Majorana fermions. Based on symmetry arguments and a renormalization group analysis, we show there are generically two phase transitions from Bogoliubov-de Gennes Chern number $\mathcal{N}=0$ to $\mathcal{N}=1$ ($p+ip$ chiral topological superconductor) and then to $\mathcal{N}=2$, in agreement with the conclusion from the band theory without disorders. Further, we discuss the critical scaling behavior of the $e^2/2h$ conductance half plateau induced by $\mathcal{N}=1$ chiral topological superconductor recently observed in the experiment. In particular, we compare the critical behavior of the half plateau induced by topological superconductor with that predicted recently by alternative explanations of the half plateau, and show that they can be distinguished in experiments.


title: Nonnegative HMM for Babble Noise Derived from Speech HMM: Application to Speech Enhancement

abstract: Deriving a good model for multitalker babble noise can facilitate different speech processing algorithms, e.g. noise reduction, to reduce the so-called cocktail party difficulty. In the available systems, the fact that the babble waveform is generated as a sum of N different speech waveforms is not exploited explicitly. In this paper, first we develop a gamma hidden Markov model for power spectra of the speech signal, and then formulate it as a sparse nonnegative matrix factorization (NMF). Second, the sparse NMF is extended by relaxing the sparsity constraint, and a novel model for babble noise (gamma nonnegative HMM) is proposed in which the babble basis matrix is the same as the speech basis matrix, and only the activation factors (weights) of the basis vectors are different for the two signals over time. Finally, a noise reduction algorithm is proposed using the derived speech and babble models. All of the stationary model parameters are estimated using the expectation-maximization (EM) algorithm, whereas the time-varying parameters, i.e. the gain parameters of speech and babble signals, are estimated using a recursive EM algorithm. The objective and subjective listening evaluations show that the proposed babble model and the final noise reduction algorithm significantly outperform the conventional methods.


title: Twofold twist defect chains at criticality

abstract: The twofold twist defects in the $D(\mathbb{Z}_k)$ quantum double model (abelian topological phase) carry non-abelian fractional Majorana-like characteristics. We align these twist defects in a line and construct a one dimensional Hamiltonian which only includes the pairwise interaction. For the defect chain with even number of twist defects, it is equivalent to the $\mathbb{Z}_k$ clock model with periodic boundary condition (up to some phase factor for boundary term), while for odd number case, it maps to $\mathbb{Z}_k$ clock model with duality twisted boundary condition. At critical point, for both cases, the twist defect chain enjoys an additional translation symmetry, which corresponds to the Kramers-Wannier duality symmetry in the $\mathbb{Z}_k$ clock model and can be generated by a series of braiding operators for twist defects. We further numerically investigate the low energy excitation spectrum for $k=3,~4,~5$ and $6$. For even-defect chain, the critical points are the same as the $\mathbb{Z}_k$ clock conformal field theories (CFTs), while for odd-defect chain, when $k\neq 4$, the critical points correspond to orbifolding a $\mathbb{Z}_2$ symmetry of CFTs of the even-defect chain. For $k=4$ case, we numerically observe some similarity to the $\mathbb{Z}_4$ twist fields in $SU(2)_1/D_4$ orbifold CFT.


title: Exact combinatorial approach to finite coagulating systems

abstract: The paper outlines an exact combinatorial approach to finite coagulating systems. In this approach, cluster sizes and time are discrete, and the binary aggregation alone governs the time evolution of the systems. By considering the growth histories of all possible clusters, the exact expression is derived for the probability of a coagulating system with an arbitrary kernel being found in a given cluster configuration when monodisperse initial conditions are applied. Then, this probability is used to calculate the time-dependent distribution for the number of clusters of a given size, the average number of such clusters and that average's standard deviation. The correctness of our general expressions is proved based on the (analytical and numerical) results obtained for systems with the constant kernel. In addition, the results obtained are compared with the results arising from the solutions to the mean-field Smoluchowski coagulation equation, indicating its weak points. The paper closes with a brief discussion on the extensibility to other systems of the approach presented herein, emphasizing the issue of arbitrary initial conditions.


title: Efficient Statistically Accurate Algorithms for the Fokker-Planck Equation in Large Dimensions

abstract: Solving the Fokker-Planck equation for high-dimensional complex turbulent dynamical systems is an important and practical issue. However, most traditional methods suffer from the curse of dimensionality and have difficulties in capturing the fat tailed highly intermittent probability density functions (PDFs) of complex systems in turbulence, neuroscience and excitable media. In this article, efficient statistically accurate algorithms are developed for solving both the transient and the equilibrium solutions of Fokker-Planck equations associated with high-dimensional nonlinear turbulent dynamical systems with conditional Gaussian structures. The algorithms involve a hybrid strategy that requires only a small number of ensembles. Here, a conditional Gaussian mixture in a high-dimensional subspace via an extremely efficient parametric method is combined with a judicious non-parametric Gaussian kernel density estimation in the remaining low-dimensional subspace. Particularly, the parametric method provides closed analytical formulae for determining the conditional Gaussian distributions in the high-dimensional subspace and is therefore computationally efficient and accurate. The full non-Gaussian PDF of the system is then given by a Gaussian mixture. Different from the traditional particle methods, each conditional Gaussian distribution here covers a significant portion of the high-dimensional PDF. Therefore a small number of ensembles is sufficient to recover the full PDF, which overcomes the curse of dimensionality. Notably, the mixture distribution has a significant skill in capturing the transient behavior with fat tails of the high-dimensional non-Gaussian PDFs, and this facilitates the algorithms in accurately describing the intermittency and extreme events in complex turbulent systems.


title: Data Innovation for International Development: An overview of natural language processing for qualitative data analysis

abstract: Availability, collection and access to quantitative data, as well as its limitations, often make qualitative data the resource upon which development programs heavily rely. Both traditional interview data and social media analysis can provide rich contextual information and are essential for research, appraisal, monitoring and evaluation. These data may be difficult to process and analyze both systematically and at scale. This, in turn, limits the ability of timely data driven decision-making which is essential in fast evolving complex social systems. In this paper, we discuss the potential of using natural language processing to systematize analysis of qualitative data, and to inform quick decision-making in the development context. We illustrate this with interview data generated in a format of micro-narratives for the UNDP Fragments of Impact project.


title: Catalog of fine-structured electron velocity distribution functions $-$ Part 1: Antiparallel magnetic-field reconnection (Geospace Environmental Modeling case)

abstract: To understand the essential physics needed to reproduce magnetic reconnection events in 2.5-D particle-in-cell (PIC) simulations, we revisit the Geospace Environmental Modeling (GEM) setup. We set up a 2-D Harris current sheet (that also specifies the initial conditions) to evolve the reconnection of antiparallel magnetic fields. In contrast to the GEM setup, we use a much smaller initial perturbation to trigger the reconnection and evolve it more self-consistently. From PIC simulation data with high-quality particle statistics we study a symmetric reconnection site, including separatrix layers, as well as the inflow and the outflow regions. The velocity distribution functions (VDFs) of electrons have a fine structure and vary strongly depending on their location within the reconnection setup. The goal is to start cataloging multidimensional fine-structured electron velocity distributions showing different reconnection processes in the Earth's magnetotail under various conditions. This will enable a direct comparison with observations from, e.g., the NASA Magnetospheric MultiScale (MMS) mission, to identify reconnection-related events. We find regions with strong non-gyrotropy also near the separatrix layer and provide a refined criterion to identify an electron diffusion region in the magnetotail. The good statistical significance of this work for relatively small analysis areas reveals the gradual changes within the fine structure of electron VDFs depending on their sampling site.

categories: astro-ph.EP astro-ph.SR physics.plasm-ph physics.space-ph


title: Primordial black hole dark matter from single field inflation

abstract: We propose a model of inflation capable of generating a population of light black holes (about $10^{-16}$ - $10^{-14}$ solar masses) that might account for a significant fraction of the dark matter in the Universe. The effective potential of the model features an approximate inflection point arising from two-loop order logarithmic corrections in well-motivated and perturbative particle physics examples. This feature decelerates the inflaton before the end of inflation, enhancing the primordial spectrum of scalar fluctuations and triggering efficient black hole production with a peaked mass distribution. At larger field values, inflation occurs thanks to a generic small coupling between the inflaton and the curvature of spacetime. We compute accurately the peak mass and abundance of the primordial black holes using the Press-Schechter and Mukhanov-Sasaki formalisms, showing that the slow-roll approximation fails to reproduce the correct results by orders of magnitude. We study as well a qualitatively similar implementation of the idea, where the approximate inflection point is due to competing terms in a generic polynomial potential. In both models, requiring a significant part of the dark matter abundance to be in the form of black holes implies a small blue scalar tilt with a sizable negative running and a tensor spectrum that may be detected by the next-generation probes of the cosmic microwave background. We also comment on previous works on the topic.


title: Weibel instability in hot plasma flows with production of gamma-rays and electron-positron pairs

abstract: We present the results of theoretical analysis and numerical simulations of the Weibel instability in two counter-streaming hot relativistic plasma flows, e.g. flows of electron-proton plasma having rest-mass density $\rho \sim 10^{-4}\; \text{g}\, \text{cm}^{-3}$, Lorentz factors $\Gamma \sim 10$ and proper temperature $T \sim 10^{13}\; \text{K}$. The instability growth rate and the filament size at the linear stage are found analytically, and are in qualitative agreement with results of three-dimensional particle-in-cell simulations. In the simulations, incoherent synchrotron emission and pair photoproduction in electromagnetic fields are taken into account. If the plasma flows are dense, fast and/or hot enough, the overall energy of synchrotron photons can be much larger than the energy of generated electromagnetic fields. Furthermore, a sizable number of positrons can be produced due to the pair photoproduction in the generated magnetic field. We propose a rough criterion for judging copious pair production and synchrotron losses. By means of this criterion we conclude that incoherent synchrotron emission and pair production during the Weibel instability can have implications for the collapsar model of gamma-ray bursts.


title: The MEGaN project II. Gravitational waves from intermediate mass- and binary black holes around a supermassive black hole

abstract: We investigate the evolution of intermediate-mass (IMBHs), stellar (BHs) and binary black holes (BHBs), deposited near a supermassive black hole (SMBH) by a population of massive star clusters. Stellar BHs rapidly segregate around the SMBH, driving the formation of extreme mass-ratio inspirals that coalesce at a rate $\Gamma= 0.02-0.2$ yr$^{-1}$ Gpc$^{-3}$ at redshift $z=0$. A few IMBHs orbiting the SMBH favour the formation of massive pairs that coalescence within a Hubble time, being the merger rate for this channel $\Gamma =0.03$ yr$^{-1}$ Gpc$^{-3}$. Recoiling kicks post-merger can eject the remnant from the galaxy centre, especially in dwarf galaxies. Our results suggest that this mechanism can lead to up to $10^5$ ejected SMBH within 1 Gpc. An IMBH co-existing with a few single and binary BHs in the same cluster can affect significantly their evolution, either driving binary disruption, yielding to intermediate-mass ratio inspirals (merger rate $\Gamma =9.5$ yr$^{-1}$ Gpc$^{-3}$), or boosting BHBs coalescence ($\Gamma =2-8$ yr$^{-1}$ Gpc$^{-3}$). In a few simulations, the SMBH boosts BHBs coalescence, leading this process to a merger rate $\Gamma =1$ yr$^{-1}$ Gpc$^{-3}$. We note that BHBs experiencing a merger in a galactic nucleus can be erroneously estimated $\sim 30\%$ heavier than it really is because of the Doppler shift of the wave frequency as caused by the rapid motion around the SMBH. All our simulations are carried out using an $N$-body code tailored to treat close encounters and post-Newtonian dynamics, that includes also the galaxy field and dynamical friction in the particles' equation of motion.


title: Defects in the supersymmetric mKdV hierarchy via Backlund transformations

abstract: The integrability of the ${\cal N}=1$ supersymmetric modified Korteweg de-Vries (smKdV) hierarchy in the presence of defects is investigated through the construction of its super B\"acklund transformation. The construction of such transformation is performed by using essentially two methods: the B\"acklund-defect matrix approach and the superfield approach. Firstly, we employ the defect matrix associated to the hierarchy which turns out to be the same for the supersymmetric sinh-Gordon (sshG) model. The method is general for all flows and as an example we derive explicitly the B\"acklund equations in components for the first few flows of the hierarchy, namely $t_3$ and $t_5$. Secondly, the supersymmetric extension of the B\"acklund transformation in the superspace formalism is constructed for those flows. Finally, this super B\"acklund transformation is employed to introduce type I defects for the supersymmetric mKdV hierarchy. Further integrability aspects by considering modified conserved quantities are derived from the defect matrix.


title: Generalized eigenstate typicality in translation-invariant quasifree fermionic models

abstract: We demonstrate a generalized notion of eigenstate thermalization for translation-invariant quasifree fermionic models: the vast majority of eigenstates satisfying a finite number of suitable constraints (e.g. fixed energy and particle number) have the property that their reduced density matrix on small subsystems approximates the corresponding generalized Gibbs ensemble. To this end, we generalize analytic results by Lai and Yang (Phys. Rev. B 91, 081110 (2015)) and illustrate the claim numerically by example of the Jordan-Wigner transform of the XX spin chain.


title: A Well Defined Glass State Obtained by Oscillatory Shear

abstract: We investigate the process of shear melting and re-solidification of a colloidal glass, directly after loading (pre-yielding) and after a series of consecutive strain sweeps (post-yielding). The post-yielding glass shows a significant softening compared to the pre-yielding glass, together with the absence of history effects in successive shear melting protocols, indicating a reproducible process of fluidisation and re-solidification into a glass state unaffected by residual stresses. However, a significant hysteresis characterises strain sweeps with increasing or decreasing strain amplitude. The appearance of history and hysteresis effects coincides with the formation of a glass state, whereas it is not observed in the liquid. We can describe the onset of shear melting over a broad range of volume fractions and frequencies using a recently developed model which describes the yielding process in terms of loss of long-lived nearest neighbours.


title: Experimental Verification of Position-Dependent Angular-Momentum Selection Rules for Absorption of Twisted Light by a Bound Electron

abstract: We analyze the multipole excitation of atoms with twisted light, i.e., by a vortex light field that carries orbital angular momentum. A single trapped $^{40}$Ca$^+$ ion serves as a localized and positioned probe of the exciting field. We drive the $S_{1/2} \to D_{5/2}$ transition and observe the relative strengths of different transitions, depending on the ion's transversal position with respect to the center of the vortex light field. On the other hand, transition amplitudes are calculated for a twisted light field in form of a Bessel beam, a Bessel-Gauss and a Gauss-Laguerre mode. Analyzing experimental obtained transition amplitudes we find agreement with the theoretical predictions at a level of better than 3\%. Finally, we propose measurement schemes with two-ion crystals to enhance the sensing accuracy of vortex modes in future experiments.


title: Boundary Observer for Space and Time Dependent Reaction-Advection-Diffusion Equations

abstract: This paper presents boundary observer design for space and time dependent reaction-advection-diffusion equations using backstepping method. The method uses only a single measurement at the boundary of the systems. The existence of the observer kernel equation is proved using the method of successive approximation.


title: Luminous type IIP SN 2013ej with high-velocity Ni-56 ejecta

abstract: We explore the well-observed type IIP SN 2013ej with peculiar luminosity evolution. It is found that the hydrodynamic model cannot reproduce in detail the bolometric luminosity at both the plateau and the radioactive tail. Yet the ejecta mass of 23-26 Msun and the kinetic energy of (1.2-1.4)x10^51 erg are determined rather confidently. We suggest that the controversy revealed in hydrodynamic simulations stems from the strong asphericity of the Ni-56 ejecta. An analysis of the asymmetric nebular H-alpha line and of the peculiar radioactive tail made it possible to recover parameters of the asymmetric bipolar Ni-56 ejecta with the heavier jet residing in the rear hemisphere. The inferred Ni-56 mass is 0.039 Msun, twice as large compared to a straightforward estimate from the bolometric luminosity at the early radioactive tail. The bulk of ejected Ni-56 has velocities in the range of 4000-6500 km/s. The linear polarization predicted by the model with the asymmetric ionization produced by bipolar Ni-56 ejecta is consistent with the observational value.


title: The Effectve Potential in Massless Theories

abstract: The effective potential V in a massless self-coupled scalar theory and massless scalar electrodynamics is considered. Both the MS and Coleman-Weinberg renormalization schemes are examined. The renormalization scheme dependence of V is determined. Upon summing all of the logarithmic contributions to V, it is shown that the implicit and explicit dependence on the renormalization scale {\mu} cancels. In addition, if there is spontaneous symmetry breaking, then the dependence on the background field {\Phi} cancels, leaving V flat but with non-perturbative contributions. The quartic scalar coupling in the Coleman-Weinberg renormalization scheme consequently vanishes.


title: Wavepackets in inhomogeneous periodic media: propagation through a one-dimensional band crossing

abstract: We consider a model of an electron in a crystal moving under the influence of an external electric field: Schroedinger's equation in one spatial dimension with a potential which is the sum of a periodic function $V$ and a smooth function $W$. We assume that the period of $V$ is much shorter than the scale of variation of $W$ and denote the ratio of these scales by $\epsilon$. We consider the dynamics of $\textit{semiclassical wavepacket}$ asymptotic (in the limit $\epsilon \downarrow 0$) solutions which are spectrally localized near to a $\textit{crossing}$ of two Bloch band dispersion functions of the periodic operator $- \frac{1}{2} \partial_z^2 + V(z)$. We show that the dynamics is qualitatively different from the case where bands are well-separated: at the time the wavepacket is incident on the band crossing, a second wavepacket is `excited' which has $\textit{opposite}$ group velocity to the incident wavepacket. We then show that our result is consistent with the solution of a `Landau-Zener'-type model.

categories: math-ph cond-mat.mes-hall math.AP math.MP quant-ph


title: SKOS Concepts and Natural Language Concepts: an Analysis of Latent Relationships in KOSs

abstract: The vehicle to represent Knowledge Organization Systems (KOSs) in the environment of the Semantic Web and linked data is the Simple Knowledge Organization System (SKOS). SKOS provides a way to assign a URI to each concept, and this URI functions as a surrogate for the concept. This fact makes of main concern the need to clarify the URIs' ontological meaning. The aim of this study is to investigate the relation between the ontological substance of KOS concepts and concepts revealed through the grammatical and syntactic formalisms of natural language. For this purpose, we examined the dividableness of concepts in specific KOSs (i.e. a thesaurus, a subject headings system and a classification scheme) by applying Natural Language Processing (NLP) techniques (i.e. morphosyntactic analysis) to the lexical representations (i.e. RDF literals) of SKOS concepts. The results of the comparative analysis reveal that, despite the use of multi-word units, thesauri tend to represent concepts in a way that can hardly be further divided conceptually, while Subject Headings and Classification Schemes - to a certain extent - comprise terms that can be decomposed into more conceptual constituents. Consequently, SKOS concepts deriving from thesauri are more likely to represent atomic conceptual units and thus be more appropriate tools for inference and reasoning. Since identifiers represent the meaning of a concept, complex concepts are neither the most appropriate nor the most efficient way of modelling a KOS for the Semantic Web.


title: Vacuum thin shells in Einstein-Gauss-Bonnet brane-world cosmology

abstract: In this paper we construct new solutions of the Einstein-Gauss-Bonnet field equations in an isotropic Shiromizu-Maeda-Sasaki brane-world setting which represent a couple of $Z_2$-symmetric vacuum thin shells splitting from the central brane, and explore the main properties of the dynamics of the system. The matching of the separating vacuum shells with the brane-world is as smooth as possible and all matter fields are restricted to the brane. We prove the existence of these solutions, derive the criteria for their existence, analyse some fundamental aspects or their evolution and demonstrate the possibility of constructing cosmological examples that exhibit this feature at early times. We also comment on the possible implications for cosmology and the relation of this system with the thermodynamic instability of highly symmetric vacuum solutions of Lovelock theory.


title: Modified-Bloch Equation Based on Fractal Derivative for Analyzing PFG Anomalous Diffusion

abstract: A modified-Bloch equation based on the fractal derivative is proposed to analyze pulsed field gradient (PFG) anomalous diffusion. Anomalous diffusion exists in many systems such as in polymer or biological systems. PFG anomalous diffusion could be analyzed based on the fractal derivative or the fractional derivative. Compared to the fractional derivative, the fractal derivative is simpler, and it is faster in numerical evaluations. In this paper, the fractal derivative is employed to build the modified-Bloch equation that is a fundamental method to describe the spin magnetization evolution affected by fractional diffusion, Larmor precession, and relaxation. An equivalent form of the fractal derivative is proposed to convert the fractional diffusion equation, which can then be combined with the precession and relaxation equations to get the modified-Bloch equation. This modified-Bloch equation yields a general PFG signal attenuation expression that includes the finite gradient pulse width (FGPW) effect, namely, the signal attenuation during field gradient pulse. The FGPW effect needs to be considered in most clinical MRI applications, and including FGPW effect allows the detecting of slower diffusion that is often encountered in polymer systems. Additionally, the spin-spin relaxation effect can be analyzed, which provides a broad view of the dynamic process in materials. The modified-Bloch equation based on the fractal derivative could provide a fundamental theoretical model for PFG anomalous diffusion.


title: Cosmological scenarios in quantum cosmology of Brans-Dicke theory of gravity on Einstein and Jordan frames

abstract: The Hamiltonian of any quantum cosmological model must obey some conditions in order to be self-adjoint, or to admit self-adjoint extensions. We will consider in this paper the self-adjoint character of the Hamiltonian resulting from the Brans-Dicke theory, both in Einstein and Jordan frames. In both cases we find the wave-function solutions in order to obtain concrete predictions for the evolution of the Universe. We show that the specific cosmological scenarios obtained depend very weakly on the conditions of having a self-adjoint operator. The problem of equivalence between Einstein and Jordan frames is considered, and it is shown that this equivalence implies a specific ordering parameter, as well as a particular choice of the physical variables.


title: Natural extensions and Gauss measures for piecewise homographic continued fractions

abstract: We give a heuristic method to solve explicitly for an absolutely continuous invariant measure for a piecewise differentiable, expanding map of a compact subset $I$ of Euclidean space $R^d$. The method consists of constructing a skew product family of maps on $I\times R^d$, which has an attractor. Lebesgue measure is invariant for the skew product family restricted to this attractor. Under reasonable measure theoretic conditions, integration over the fibers gives the desired measure on $I$. Furthermore, the attractor system is then the natural extension of the original map with this measure. We illustrate this method by relating it to various results in the literature.


title: MultiNet: Multi-Modal Multi-Task Learning for Autonomous Driving

abstract: Autonomous driving requires operation in different behavioral modes ranging from lane following and intersection crossing to turning and stopping. However, most existing deep learning approaches to autonomous driving do not consider the behavioral mode in the training strategy. This paper describes a technique for learning multiple distinct behavioral modes in a single deep neural network through the use of multi-modal multi-task learning. We study the effectiveness of this approach, denoted MultiNet, using self-driving model cars for driving in unstructured environments such as sidewalks and unpaved roads. Using labeled data from over one hundred hours of driving our fleet of 1/10th scale model cars, we trained different neural networks to predict the steering angle and driving speed of the vehicle in different behavioral modes. We show that in each case, MultiNet networks outperform networks trained on individual modes while using a fraction of the total number of parameters.


title: Electronic properties of bilayer phosphorene quantum dots in the presence of perpendicular electric and magnetic fields

abstract: Using the tight-binding approach, we investigate the electronic properties of bilayer phosphorene (BLP) quantum dots (QDs) in the presence of perpendicular electric and magnetic fields. Since BLP consists of two coupled phosphorene layers, it is of interest to examine the layer-dependent electronic properties of BLP QDs, such as the electronic distributions over the two layers and the so-produced layer-polarization features, and to see how these properties are affected by the magnetic field and the bias potential. We find that in the absence of a bias potential only edge states are layer-polarized while the bulk states are not, and the layer-polarization degree (LPD) of the unbiased edge states increases with increasing magnetic field. However, in the presence of a bias potential both the edge and bulk states are layer-polarized, and the LPD of the bulk (edge) states depends strongly (weakly) on the interplay of the bias potential and the interlayer coupling. At high magnetic fields, applying a bias potential renders the bulk electrons in a BLP QD to be mainly distributed over the top or bottom layer, resulting in layer-polarized bulk Landau levels (LLs). In the presence of a large bias potential that can drive a semiconductor-to-semimetal transition in BLP, these bulk LLs exhibit different magnetic-field dependences, i.e., the zeroth LLs exhibit a linear-like dependence on the magnetic field while the other LLs exhibit a square-root-like dependence.


title: Mitigating Evasion Attacks to Deep Neural Networks via Region-based Classification

abstract: Deep neural networks (DNNs) have transformed several artificial intelligence research areas including computer vision, speech recognition, and natural language processing. However, recent studies demonstrated that DNNs are vulnerable to adversarial manipulations at testing time. Specifically, suppose we have a testing example, whose label can be correctly predicted by a DNN classifier. An attacker can add a small carefully crafted noise to the testing example such that the DNN classifier predicts an incorrect label, where the crafted testing example is called adversarial example. Such attacks are called evasion attacks. Evasion attacks are one of the biggest challenges for deploying DNNs in safety and security critical applications such as self-driving cars. In this work, we develop new methods to defend against evasion attacks. Our key observation is that adversarial examples are close to the classification boundary. Therefore, we propose region-based classification to be robust to adversarial examples. For a benign/adversarial testing example, we ensemble information in a hypercube centered at the example to predict its label. In contrast, traditional classifiers are point-based classification, i.e., given a testing example, the classifier predicts its label based on the testing example alone. Our evaluation results on MNIST and CIFAR-10 datasets demonstrate that our region-based classification can significantly mitigate evasion attacks without sacrificing classification accuracy on benign examples. Specifically, our region-based classification achieves the same classification accuracy on testing benign examples as point-based classification, but our region-based classification is significantly more robust than point-based classification to various evasion attacks.


title: Representation Learning on Graphs: Methods and Applications

abstract: Machine learning on graphs is an important and ubiquitous task with applications ranging from drug design to friendship recommendation in social networks. The primary challenge in this domain is finding a way to represent, or encode, graph structure so that it can be easily exploited by machine learning models. Traditionally, machine learning approaches relied on user-defined heuristics to extract features encoding structural information about a graph (e.g., degree statistics or kernel functions). However, recent years have seen a surge in approaches that automatically learn to encode graph structure into low-dimensional embeddings, using techniques based on deep learning and nonlinear dimensionality reduction. Here we provide a conceptual review of key advancements in this area of representation learning on graphs, including matrix factorization-based methods, random-walk based algorithms, and graph neural networks. We review methods to embed individual nodes as well as approaches to embed entire (sub)graphs. In doing so, we develop a unified framework to describe these recent approaches, and we highlight a number of important applications and directions for future work.


title: Rigorous Analysis for Efficient Statistically Accurate Algorithms for Solving Fokker-Planck Equations in Large Dimensions

abstract: This article presents a rigorous analysis for efficient statistically accurate algorithms for solving the Fokker-Planck equations associated with high-dimensional nonlinear turbulent dynamical systems with conditional Gaussian structures. Despite the conditional Gaussianity, these nonlinear systems contain many strong non-Gaussian features such as intermittency and fat-tailed probability density functions (PDFs). The algorithms involve a hybrid strategy that requires only a small number of samples $L$ to capture both the transient and the equilibrium non-Gaussian PDFs with high accuracy. Here, a conditional Gaussian mixture in a high-dimensional subspace via an extremely efficient parametric method is combined with a judicious Gaussian kernel density estimation in the remaining low-dimensional subspace. Rigorous analysis shows that the mean integrated squared error in the recovered PDFs in the high-dimensional subspace is bounded by the inverse square root of the determinant of the conditional covariance, where the conditional covariance is completely determined by the underlying dynamics and is independent of $L$. This is fundamentally different from a direct application of kernel methods to solve the full PDF, where $L$ needs to increase exponentially with the dimension of the system and the bandwidth shrinks. A detailed comparison between different methods justifies that the efficient statistically accurate algorithms are able to overcome the curse of dimensionality. It is also shown with mathematical rigour that these algorithms are robust in long time provided that the system is controllable and stochastically stable. Particularly, dynamical systems with energy-conserving quadratic nonlinearity as in many geophysical and engineering turbulence are proved to have these properties.


title: Generalized PMC model for the hybrid diagnosis of multiprocessor systems

abstract: Fault diagnosis is important to the design and maintenance of large multiprocessor systems. PMC model is the most famous diagnosis model in the system level diagnosis of multiprocessor systems. Under the PMC model, only node faults are allowed. But in real circumstances, link faults may occur. So based on the PMC model, we propose in this paper a diagnosis model called the generalized PMC(GPMC) model to adapt to the real circumstances. The foundation of GPMC model has been established. And to measure the fault diagnosis capability of multiprocessor systems under the GPMC model, the fault diagnosis capability measuring parameters: $h$-edge restricted diagnosability and $h$-vertex restricted edge diagnosability have been introduced. As an application, the $h$-edge restricted diagnosability and $h$-vertex restricted edge diagnosability of hypercubes are explored.


title: Character Distributions of Classical Chinese Literary Texts: Zipf's Law, Genres, and Epochs

abstract: We collect 14 representative corpora for major periods in Chinese history in this study. These corpora include poetic works produced in several dynasties, novels of the Ming and Qing dynasties, and essays and news reports written in modern Chinese. The time span of these corpora ranges between 1046 BCE and 2007 CE. We analyze their character and word distributions from the viewpoint of the Zipf's law, and look for factors that affect the deviations and similarities between their Zipfian curves. Genres and epochs demonstrated their influences in our analyses. Specifically, the character distributions for poetic works of between 618 CE and 1644 CE exhibit striking similarity. In addition, although texts of the same dynasty may tend to use the same set of characters, their character distributions still deviate from each other.


title: Hybrid Fault diagnosis capability analysis of Hypercubes under the PMC model and MM* model

abstract: System level diagnosis is an important approach for the fault diagnosis of multiprocessor systems. In system level diagnosis, diagnosability is an important measure of the diagnosis capability of interconnection networks. But as a measure, diagnosability can not reflect the diagnosis capability of multiprocessor systems to link faults which may occur in real circumstances. In this paper, we propose the definition of $h$-edge tolerable diagnosability to better measure the diagnosis capability of interconnection networks under hybrid fault circumstances. The $h$-edge tolerable diagnosability of a multiprocessor system $G$ is the maximum number of faulty nodes that the system can guarantee to locate when the number of faulty edges does not exceed $h$,denoted by $t_h^{e}(G)$. The PMC model and MM model are the two most widely studied diagnosis models for the system level diagnosis of multiprocessor systems. The hypercubes are the most well-known interconnection networks. In this paper, the $h$-edge tolerable diagnosability of $n$-dimensional hypercube under the PMC model and MM$^{*}$ is determined as follows: $t_h^{e}(Q_n)= n-h$, where $1\leq h<n$, $n\geq3$.


title: Two electrostatic gyrokinetic models derived by two different perturbative methods

abstract: This paper presents two different electrostatic gyrokinetic models derived through two different perturbative methods. One of the two models is just the conventional electrostatic gyrokinetic model, the derivation of which is repeated using the Lie transform perturbative method. One term is rectified in the derivation of the conventional model. To derive the other model, we use a new method, which is based on the covariant transform formula of the differential 1-form. The new method doesn't split the coordinate transform into the guiding-center transform and the gyrocenter transform. It carries out the coordinate transform up to the order equaling that of the amplitude of the perturbative wave. Compared with the conventional model, the finite Larmor radius terms are completely removed from the orbit equations of the new one, making it simpler for the numerical application.


title: Millimeter Wave Channel Measurements and Implications for PHY Layer Design

abstract: There has been an increasing interest in the millimeter wave (mmW) frequency regime in the design of next-generation wireless systems. The focus of this work is on understanding mmW channel properties that have an important bearing on the feasibility of mmW systems in practice and have a significant impact on physical (PHY) layer design. In this direction, simultaneous channel sounding measurements at 2.9, 29 and 61 GHz are performed at a number of transmit-receive location pairs in indoor office, shopping mall and outdoor environments. Based on these measurements, this paper first studies large-scale properties such as path loss and delay spread across different carrier frequencies in these scenarios. Towards the goal of understanding the feasibility of outdoor-to-indoor coverage, material measurements corresponding to mmW reflection and penetration are studied and significant notches in signal reception spread over a few GHz are reported. Finally, implications of these measurements on system design are discussed and multiple solutions are proposed to overcome these impairments.


title: On density of infinite subsets I

abstract: Let $Y$ be a compact metric space, $G$ be a group acting by transformations on $Y$. For any infinite subset $A\subset Y$, we study the density of $gA$ for $g\in G$ and quantitative density of the set $\displaystyle{\bigcup_{g\in G_n}gA}$ by the Hausdorff semimetric $d^H$. It is proven that for any integer $n\ge 2$, $\epsilon>0$, any infinite subset $A\subset \mathbb T^n$, there is a $g\in SL(n,\mathbb Z)$ such that $gA$ is $\epsilon$-dense. We also show that, for any infinite subset $A\subset [0,1]$, for generic rotation and generic 3-IET, $$\liminf_nn\cdot d^H\left(\bigcup_{k=0}^{n-1}T^kA,[0,1]\right)=0.$$


title: Computation of graphical derivatives of normal cone maps to a class of conic constraint sets

abstract: This paper concerns with the graphical derivative of the normals to the conic constraint $g(x)\in\!K$, where $g\!:\mathbb{X}\to\mathbb{Y}$ is a twice continuously differentiable mapping and $K\subseteq\mathbb{Y}$ is a nonempty closed convex set assumed to be $C^2$-cone reducible. Such a generalized derivative plays a crucial role in characterizing isolated calmness of the solution maps to generalized equations whose multivalued parts are modeled via the normals to the nonconvex set $\Gamma=g^{-1}(K)$. The main contribution of this paper is to provide an exact characterization for the graphical derivative of the normals to this class of nonconvex conic constraints under an assumption without requiring the nondegeneracy of the reference point as the papers \cite{Gfrerer17,Mordu15,Mordu151} do.


title: On density of infinite subsets II: dynamics on homogeneous spaces

abstract: Let $G$ be a noncompact semisimple Lie group, $\Gamma$ be an irreducible cocompact lattice in $G$, and $P<G$ be a minimal parabolic subgroup. We consider the dynamics of $P$ acting on $G/\Gamma$ by left translation. For any infinite subset $A\subset G/\Gamma$, we show that, for any $\epsilon>0$, there is a $g\in P$ such that $gA$ is $\epsilon$-dense. We also prove a similar result for certain discrete group actions on $\mathbb T^n$.


title: GDP growth rates as confined L\'evy flights

abstract: A new model that combines economic growth rate fluctuations at the microscopic and macroscopic level is presented. At the microscopic level, firms are growing at different rates while also being exposed to idiosyncratic shocks at the firm and sector level. We describe such fluctuations as independent L\'evy-stable fluctuations, varying over multiple orders of magnitude. These fluctuations are aggregated and measured at the macroscopic level in averaged economic output quantities such as GDP. A fundamental question is thereby to what extend individual firm size fluctuations can have a noticeable impact on the overall economy. We argue that this question can be answered by considering the L\'evy fluctuations as embedded in a steep confining potential well, ensuring nonlinear mean-reversal behavior, without having to rely on microscopic details of the system. The steepness of the potential well directly controls the extend towards which idiosyncratic shocks to firms and sectors are damped at the level of the economy. Additionally, the theory naturally accounts for business cycles, represented in terms of a bimodal economic output distribution, and thus connects two so far unrelated fields in economics. By analyzing 200 years of US GDP growth rates, we find that the model is in good agreement with the data.


title: Chromatic numbers and a Lov\'asz type inequality for non-commutative graphs

abstract: Non-commutative graph theory is an operator space generalization of graph theory. Well known graph parameters such as the independence number and Lov\'asz theta function were first generalized to this setting by Duan, Severini, and Winter. We introduce two new generalizations of the chromatic number to non-commutative graphs and provide a generalization of the Lov\'asz sandwich inequality. In particular, we show the chromatic number of the orthogonal complement of a non-commutative graph is bounded below by its theta number. We also provide a generalization of both Sabadussi's Theorem and Hedetniemi's conjecture to non-commutative graphs.


title: Further insights into the damping-induced self-recovery phenomenon

abstract: In a series of papers, D. E. Chang, et al., proved and experimentally demonstrated a phenomenon they termed "damping-induced self-recovery". However, these papers left a few questions concerning the observed phenomenon unanswered - in particular, the effect of the intervening lubricant-fluid and its viscosity on the recovery, the abrupt change in behaviour with the introduction of damping, a description of the energy dynamics, and the curious occurrence of overshoots and oscillations and its dependence on the control law. In this paper we attempt to answer these questions through theory. In particular, we derive an expression for the infinite-dimensional fluid-stool-wheel system, that approximates its dynamics to that of the better understood finite-dimensional case.


title: On the Application of Optimal Control Theory to Climate Engineering

abstract: The new method is considered for identifying the perfect (optimal) scenario of climate engineering operations based on the optimal control theory. The application of this approach is demonstrated using zero-dimensional energy-balance climate model. The global average surface temperature represents the state variable and the reflective power of aerosols injected into the stratosphere is designated as the control variable. Maximum principle is used to find the best possible control and the associated climate system trajectory for a specific objective function (performance measure). Discussed illustrative results were analytically and numerically calculated using the four greenhouse gas concentration scenarios based on the Representative Concentration Pathways.