title: Monte Carlo integration with a growing number of control variates

abstract: It is well known that Monte Carlo integration with variance reduction by means of control variates can be implemented by the ordinary least squares estimator for the intercept in a multiple linear regression model. A central limit theorem is established for the integration error if the number of control variates tends to infinity. The integration error is scaled by the standard deviation of the error term in the regression model. If the linear span of the control variates is dense in a function space that contains the integrand, the integration error tends to zero at a rate which is faster than the square root of the number of Monte Carlo replicates. Depending on the situation, increasing the number of control variates may or may not be computationally more efficient than increasing the Monte Carlo sample size.


title: On the existence of an ultra central approximate identity for certain semigroup algebras

abstract: In this paper we characterize the existance of an ultra central approximate identity for $\ell^{1}(S)$, where $S$ is a uniformly locally finite inverse semigroup. As an application, for the Brandt semigroup $S=M^{0}(G,I)$ over a non-empty set $I$, we show that $\ell^{1}(S)$ has an ultra central approximate identity if and only if $I$ is finite.


title: Closed-form Marginal Likelihood in Gamma-Poisson Matrix Factorization

abstract: We present novel understandings of the Gamma-Poisson (GaP) model, a probabilistic matrix factorization model for count data. We show that GaP can be rewritten free of the score/activation matrix. This gives us new insights about the estimation of the topic/dictionary matrix by maximum marginal likelihood estimation. In particular, this explains the robustness of this estimator to over-specified values of the factorization rank, especially its ability to automatically prune irrelevant dictionary columns, as empirically observed in previous work. The marginalization of the activation matrix leads in turn to a new Monte Carlo Expectation-Maximization algorithm with favorable properties.


title: Momentum-Field Interactions Beyond Standard Quadratic Optomechanics

abstract: This chapter summarizes the recent progress in the theory and analytical tools of quadratic optomechanical interactions, as one of the prominent domains of contemporary nonlinear quantum optics. Emphasis has been put here first to show what types of nonlinear interactions do exist, and what physical interpretations follow each. The standard quadratic interactions between light and mechanical motion is expressed as the product of cavity light intensity and squared mirror position. However, there exists a non-standard quadratic optomechanical interaction as well, which assumes a mathematically different form and appears as the squared product of field and mirror momenta. This non-standard type of quadratic interaction originates from two corrections: the momentum exchange and conservation among mirror and field, as well as relativistic corrections due to different mechanisms. Both these types of non-standard interactions become relevant when the ratio of mechanical to optical frequency is no longer negligible. Next, we turn to the solution technique of such interactions, and introduce a formal higher-order operator method to tackle the nonlinear evolution of quantum systems. This enables one to accurately study any type of quantum nonlinear interaction using the analysis tools of linear algebra. In order to employ the analytical power of higher-order operator method, one first needs to identify a closed Lie algebra, which should satisfy closedness property under commutation either exactly or approximately, and is referred to as the basis. Having the basis of higher-order operators known, one may proceed to construct the corresponding Langevin equations, which can be now conveniently analyzed using the existing mathematical toolbox of linear algebra to yield the spectral densities, moments, and expectation values.


title: On a connectedness principle of Shokurov-Koll\'ar type

abstract: Let $(X,\Delta)$ be a log pair over $S$, such that $-(K_X+\Delta)$ is nef over $S$. It is conjectured that the intersection of the non-klt (non Kawamata log terminal) locus of $(X,\Delta)$ with any fiber $X_s$ has at most two connected components. We prove this conjecture in dimension $\leq 4$ and in arbitrary dimension assuming the termination of klt flips.


title: Neighborhood-Prime Labelings of Trees and Other Classes of Graphs

abstract: A neighborhood-prime labeling of a graph is a variation of a prime labeling in which the vertices are assigned labels from $1$ to $|V(G)|$ such that the $\gcd$ of the labels in the neighborhood of each non-degree $1$ vertex is equal to $1$. In this paper, we examine neighborhood-prime labelings of several classes of graphs such as polygonal snakes and books, with a focus on trees including caterpillars, spiders, and firecrackers.


title: Scheduling Policies for Minimizing Age of Information in Broadcast Wireless Networks

abstract: We consider a wireless broadcast network with a base station sending time-sensitive information to a number of clients through unreliable channels. The Age of Information (AoI), namely the amount of time that elapsed since the most recently delivered packet was generated, captures the freshness of the information. We formulate a discrete-time decision problem to find a transmission scheduling policy that minimizes the expected weighted sum AoI of the clients in the network. We first show that in symmetric networks a Greedy policy, which transmits the packet with highest current age, is optimal. For general networks, we develop three low-complexity scheduling policies: a randomized policy, a Max-Weight policy and a Whittle's Index policy, and derive performance guarantees as a function of the network configuration. To the best of our knowledge, this is the first work to derive performance guarantees for scheduling policies that attempt to minimize AoI in wireless networks with unreliable channels. Numerical results show that both Max-Weight and Whittle's Index policies outperform the other scheduling policies in every configuration simulated, and achieve near optimal performance.


title: Compressibility in turbulent MHD and passive scalar transport: mean-field theory

abstract: We develop a mean-field theory of compressibility effects in turbulent magnetohydrodynamics and passive scalar transport using the quasi-linear approximation and the spectral $\tau$-approach. We find that compressibility decreases the $\alpha$ effect and the turbulent magnetic diffusivity both at small and large magnetic Reynolds numbers, Rm. Similarly, compressibility decreases the turbulent diffusivity for passive scalars both at small and large P\'eclet numbers, Pe. On the other hand, compressibility does not affect the effective pumping velocity of the magnetic field for large Rm, but it decreases it for small Rm. Density stratification causes turbulent pumping of passive scalars, but it is found to become weaker with increasing compressibility. No such pumping effect exists for magnetic fields. However, compressibility results in a new passive scalar pumping effect from regions of low to high turbulent intensity both for small and large P\'eclet numbers. It can be interpreted as compressible turbophoresis of noninertial particles and gaseous admixtures, while the classical turbophoresis effect exists only for inertial particles and causes them to be pumped to regions with lower turbulent intensity.

categories: physics.flu-dyn astro-ph.SR physics.plasm-ph physics.space-ph


title: Methodology for the optimization of detection systems of Plasma corpuscular radiation, working under intensive conditions of n-g background radiation

abstract: This Thesis present the process of development very low sensitivity to background radiation n-g scintillation detectors with close to 100% efficiency to corpuscular diagnostic of thermonuclear plasma. The detectors developed using this methodology were installed in two major Neutral Particle Analyzers for the JET tokamak (Gemma-2M and ISEP) and they operated successfully for many years. The configuration for such detectors is based on the growing of very thin scintillators of CsI(Tl) on top of quartz supporters and then glued to photomultipliers by keeping the thickness of the supporter and the entrance window of the photocamera of the photomultiplier very thin low sensitivity to background radiation can be achieved.


title: An elliptic boundary value problem for $G_{2}$ structures

abstract: We show that the $G_{2}$ holonomy equation on a manifold with boundary, with prescribed 3-form on the boundary, is elliptic. The main point is to set up a suitable linear elliptic boundary value problem. This result leads to a deformation theory. In particular we establish the existence of certain $G_{2}$ cobordisms between two small deformations of a Calabi-Yau 3-fold.


title: A Greedy Search Tree Heuristic for Symbolic Regression

abstract: Symbolic Regression tries to find a mathematical expression that describes the relationship of a set of explanatory variables to a measured variable. The main objective is to find a model that minimizes the error and, optionally, that also minimizes the expression size. A smaller expression can be seen as an interpretable model considered a reliable decision model. This is often performed with Genetic Programming which represents their solution as expression trees. The shortcoming of this algorithm lies on this representation that defines a rugged search space and contains expressions of any size and difficulty. These pose as a challenge to find the optimal solution under computational constraints. This paper introduces a new data structure, called Interaction-Transformation (IT), that constrains the search space in order to exclude a region of larger and more complicated expressions. In order to test this data structure, it was also introduced an heuristic called SymTree. The obtained results show evidence that SymTree are capable of obtaining the optimal solution whenever the target function is within the search space of the IT data structure and competitive results when it is not. Overall, the algorithm found a good compromise between accuracy and simplicity for all the generated models.


title: On behavior of homeomorphisms with inverse modulus conditions

abstract: We consider some class of homeomorphisms of domains of Euclidean space, which are more general than quasiconformal mappings. For these homeomorphisms, we have obtained theorems on local behavior of it's inverse mappings in a given domain. Under some additional conditions, we proved results about behavior of mappings mentioned above in the closure of the domain.


title: Centrifugally forced Rayleigh-Taylor Instability

abstract: The effect of rotation on the classical gravity-driven Rayleigh-Taylor instability has been shown to influence the scale of the perturbations that develop at the unstable interface and consequently alter the speed of propagation of the front. The present authors argued that this is as a result of a competition between the destabilizing effect of gravity and the stabilizing effect of the rotation. In the present paper we consider the extreme limit of high rotation rates in which rotational forces dominate and gravitational forces may be ignored. The two liquid layers initially form concentric cylinders, centred on the axis of rotation. The configuration may be thought of as a fluid-fluid centrifuge. There are two types of perturbation to the interface that may be considered, an azimuthal perturbation around the circumference of the interface and a varicose perturbation in the axial direction along the length of the interface. It is the first of these types of perturbation that we consider here, and so the flow may be considered essentially two-dimensional, taking place in a circular domain. We carry out a linear stability analysis on a perturbation to the hydrostatic background state and derive a fourth order Orr-Sommerfeld-like equation that governs the system. We consider the dynamics of systems of stable and unstable configurations, inviscid and viscous fluids, immiscible fluid layers with surface tension, and miscible fluid layers that may have some initial diffusion of density. Theoretical predictions are compared with numerical experiments and the agreement is shown to be good. We do not restrict our analysis to equal volume fluid layers and so our results also have applications in coating and lubrication problems in rapidly rotating systems and machinery.


title: Bayesian calibration of a numerical code for prediction

abstract: Field experiments are often difficult and expensive to make. To bypass these issues, industrial companies have developed computational codes. These codes intend to be representative of the physical system, but come with a certain amount of problems. The code intends to be as close as possible to the physical system. It turns out that, despite continuous code development, the difference between the code outputs and experiments can remain significant. Two kinds of uncertainties are observed. The first one comes from the difference between the physical phenomenon and the values recorded experimentally. The second concerns the gap between the code and the physical system. To reduce this difference, often named model bias, discrepancy, or model error, computer codes are generally complexified in order to make them more realistic. These improvements lead to time consuming codes. Moreover, a code often depends on parameters to be set by the user to make the code as close as possible to field data. This estimation task is called calibration. This paper proposes a review of Bayesian calibration methods and is based on an application case which makes it possible to discuss the various methodological choices and to illustrate their divergences. This example is based on a code used to predict the power of a photovoltaic plant.


title: SABCEMM-A Simulator for Agent-Based Computational Economic Market Models

abstract: We introduce the simulation tool SABCEMM (Simulator for Agent-Based Computational Economic Market Models) for agent-based computational economic market (ABCEM) models. Our simulation tool is implemented in C++ and we can easily run ABCEM models with several million agents. The object-oriented software design enables the isolated implementation of building blocks for ABCEM models, such as agent types and market mechanisms. The user can design and compare ABCEM models in a unified environment by recombining existing building blocks using the XML-based SABCEMM configuration file. We introduce an abstract ABCEM model class which our simulation tool is built upon. Furthermore, we present the software architecture as well as computational aspects of SABCEMM. Here, we focus on the efficiency of SABCEMM with respect to the run time of our simulations. We show the great impact of different random number generators on the run time of ABCEM models. The code and documentation is published on GitHub at https://github.com/SABCEMM/SABCEMM, such that all results can be reproduced by the reader.


title: Horospheres in Teichm\"uller space and mapping class group

abstract: We study the geometry of horospheres in Teichm\"uller space of Riemann surfaces of genus g with n punctures, where $3g-3+n\geq 2$. We show that every $C^1$-diffeomorphism of Teichm\"uller space to itself that preserves horospheres is an element of the extended mapping class group. Using the relation between horospheres and metric balls, we obtain a new proof of Royden's Theorem that the isometry group of the Teichm\"uller metric is the extended mapping class group.


title: On a numerical upper bound for the extended Goldbach conjecture

abstract: The Goldbach conjecture states that every even number can be decomposed as the sum of two primes. Let $D(N)$ denote the number of such prime decompositions for an even $N$. It is known that $D(N)$ can be bounded above by $$ D(N) \leq C^* \Theta(N), \quad \Theta(N):= \frac{N}{\log^2 N}\prod_{\substack{p|N p>2}} \left( 1 + \frac{1}{p-2}\right)\prod_{p>2}\left(1-\frac{1}{(p-1)^2}\right) $$ where $C^*$ denotes Chen's constant. It is conjectured that $C^*=2$. In 2004, Wu showed that $C^* \leq 7.8209$. We attempted to replicate his work in computing Chen's constant, and in doing so we provide an improved approximation of the Buchstab function $\omega(u)$, \begin{align*} \omega(u)=1/u, & \quad (1\leq u\leq 2), (u \omega(u))'=\omega(u-1), & \quad (u\geq 2). \end{align*} based on work done by Cheer and Goldston. For each interval $[j,j+1]$, they expressed $\omega(u)$ as a Taylor expansion about $u=j+1$. We expanded about the point $u=j+0.5$, so $\omega(u)$ was never evaluated more than $0.5$ away from the center of the Taylor expansion, which gave much stronger error bounds. Issues arose while using this Taylor expansion to compute the required integrals for Chen's constant, so we proceeded with solving the above differential equation to obtain $\omega(u)$, and then integrating the result. Although the values that were obtained undershot Wu's results, we pressed on and refined Wu's work by discretising his integrals with finer granularity. The improvements to Chen's constant were negligible (as predicted by Wu). This provides experimental evidence, but not a proof, that were Wu's integrals computed on smaller intervals in exact form, the improvement to Chen's constant would be negligible. Thus, any substantial improvement on Chen's constant likely requires a radically different method to what Wu provided.


title: Automated Conjecturing VII: The Graph Brain Project & Big Mathematics

abstract: The Graph Brain Project is an experiment in how the use of automated mathematical discovery software, databases, large collaboration, and systematic investigation provide a model for how mathematical research might proceed in the future. Our Project began with the development of a program that can be used to generate invariant-relation and property-relation conjectures in many areas of mathematics. This program can produce conjectures which are not implied by existing (published) theorems. Here we propose a new approach to push forward existing mathematical research goals---using automated mathematical discovery software. We suggest how to initiate and harness large-scale collaborative mathematics. We envision mathematical research labs similar to what exist in other sciences, new avenues for funding, new opportunities for training students, and a more efficient and effective use of published mathematical research. And our experiment in graph theory can be imitated in many other areas of mathematical research. Big Mathematics is the idea of large, systematic, collaborative research on problems of existing mathematical interest. What is possible when we put our skills, tools, and results together systematically?


title: Once More on Coulomb-Nuclear Interference

abstract: This is a critical reconsideration of the standard way of account for Coulomb-nuclear interference in elastic scattering amplitude.


title: Finding the seed of uniform attachment trees

abstract: A uniform attachment tree is a random tree that is generated dynamically. Starting from a fixed "seed" tree, vertices are added sequentially by attaching each vertex to an existing vertex chosen uniformly at random. Upon observing a large (unlabeled) tree, one wishes to find the initial seed. We investigate to what extent seed trees can be recovered, at least partially. We consider three types of seeds: a path, a star, and a random uniform attachment tree. We propose and analyze seed-finding algorithms for all three types of seed trees.


title: Scaling Behavior in the 3D Random Field $XY$ Model

abstract: We have performed studies of the 3D random field $XY$ model on $L \times L \times L$ simple cubic lattices with periodic boundary conditions, with a random field strength of $h_r$ = 1.875, for $L = 64$ and $L = 96$, using a parallelized Monte Carlo algorithm. We present results for the angle-averaged magnetic structure factor, $S ( k )$ at $T = 1.00$, which appears to be the temperature at which small jumps in the magnetization per spin and the energy per spin occur. The results indicate the existence of an approximately logarithmic divergence of $S ( k )$ as $k \to 0$. This suggests that the lower critical dimension for long range order in this model is three.


title: Towards a quantum time mirror for nonrelativistic wave packets

abstract: We propose a method -- a quantum time mirror (QTM) -- for simulating a partial time-reversal of the free-space motion of a nonrelativistic quantum wave packet. The method is based on a short-time spatially-homogeneous perturbation to the wave packet dynamics, achieved by adding a nonlinear time-dependent term to the underlying Schr\"odinger equation. Numerical calculations, supporting our analytical considerations, demonstrate the effectiveness of the proposed QTM for generating a time-reversed echo image of initially localized matter-wave packets in one and two spatial dimensions. We also discuss possible experimental realizations of the proposed QTM.


title: Twisted arithmetic Siegel Weil formula on X0(N)

abstract: In this paper, we study twisted arithmetic divisors on the modular curve X_0(N) with N square-free. For each pair (\Delta, r) where \Delta >0 and \Delta \equiv r^2 \mod 4N, we constructed a twisted arithmetic theta function \phi_{\Delta, r}(\tau) which is a generating function of arithmetic twisted Heegner divisors. We prove the modularity of \phi_{\Delta, r}(\tau), along the way, we also identify the arithmetic pairing \langle \phi_{\Delta, r}(\tau),\widehat{\omega}_N \rangle with special value of some Eisenstein series, where \widehat{\omega}_N is a normalized metric Hodge line bundle.


title: Design and Implementation of a Polar Codes Blind Detection Scheme

abstract: In blind detection, a set of candidates has to be decoded within a strict time constraint, to identify which transmissions are directed at the user equipment. Blind detection is required by the 3GPP LTE/LTE-Advanced standard, and it will be required in the 5th generation wireless communication standard (5G) as well. Polar codes have been selected for use in 5G: thus, the issue of blind detection of polar codes must be addressed. We propose a polar code blind detection scheme where the user ID is transmitted instead of some of the frozen bits. A first, coarse decoding phase helps selecting a subset of candidates that is decoded by a more powerful algorithm: an early stopping criterion is also introduced for the second decoding phase. Simulations results show good missed detection and false alarm rates, along with substantial latency gains thanks to early stopping. We then propose an architecture to implement the devised blind detection scheme, based on a tunable decoder that can be used for both phases. The architecture is synthesized and implementation results are reported for various system parameters. The reported area occupation and latency, obtained in 65 nm CMOS technology, are able to meet 5G requirements, and are guaranteed to meet them with even less resource usage in the latest technology nodes.


title: Mean-field modeling of the basal ganglia-thalamocortical system. II. Dynamics of parkinsonian oscillations

abstract: Neuronal correlates of Parkinson's disease (PD) include a slowing of the electroencephalogram (EEG) and enhanced synchrony at 3-7 and 7-30 Hz in the basal ganglia, thalamus, and cortex. This study describes the dynamics of a physiologically based mean-field model of the basal ganglia-thalamocortical system, and shows how it accounts for key electrophysiological correlates of PD. Its connectivity comprises partially segregated direct and indirect pathways through the striatum, a hyperdirect pathway involving a corticosubthalamic projection, thalamostriatal feedback, and local inhibition in striatum and external pallidum (GPe). In a companion paper, realistic steady-state firing rates were obtained for the healthy state, and after dopamine loss modeled by weaker direct and stronger indirect pathways, reduced intrapallidal inhibition, lower firing thresholds of the GPe and subthalamic nucleus (STN), a stronger striato-GPe projection, and weaker cortical interactions. Here we show that oscillations around 5 and 20 Hz can arise with a strong indirect pathway, which also increases synchrony throughout the basal ganglia. Further, increased theta power with nigrostriatal degeneration correlates with reduced alpha power and peak frequency, matching experiments. Unlike the hyperdirect pathway, the indirect pathway sustains oscillations with realistic phase relationships. Changes in basal ganglia responses to transient stimuli accord with experimental data. Reduced cortical gains due to both nigrostriatal and mesocortical dopamine loss lead to slower cortical activity changes and may be related to bradykinesia. Finally, increased EEG power found in some studies may be partly explained by a lower effective GPe firing threshold, reduced GPe-GPe inhibition, and/or weaker intracortical connections in PD. Strict separation of the direct and indirect pathways is not necessary for these results.


title: Mean-field modeling of the basal ganglia-thalamocortical system. I. Firing rates in healthy and parkinsonian states

abstract: Parkinsonism leads to various electrophysiological changes in the basal ganglia-thalamocortical system (BGTCS), often including elevated discharge rates of the subthalamic nucleus (STN) and the output nuclei, and reduced activity of the globus pallidus external segment (GPe). These rate changes have been explained qualitatively in terms of the direct/indirect pathway model, involving projections of distinct striatal populations to the output nuclei and GPe. Although these populations partly overlap, evidence suggests dopamine depletion differentially affects cortico-striato-pallidal connection strengths to the two pallidal segments. Dopamine loss may also decrease the striatal signal-to-noise ratio, reducing both corticostriatal coupling and striatal firing thresholds. Here we present a mean-field model of the BGTCS with structure and parameter estimates closely based on physiology and anatomy. Changes in firing rates due to modeled dopamine loss are compared with experiment. Our results suggest that a stronger indirect pathway, possibly combined with a weakened direct pathway, is compatible with empirical evidence. However, altered corticostriatal connection strengths are probably not solely responsible for substantially increased STN activity often found. A lower STN firing threshold, weaker intracortical inhibition, and stronger striato-GPe inhibition help explain the relatively large increase in STN rate. Changes in cortex, GPe, and STN help normalize the cortical rate, also in accord with experiments. The model integrates the basal ganglia into a unified framework along with an existing thalamocortical model that accounts for a wide range of electrophysiological phenomena. A companion paper discusses the dynamics and oscillations of this combined system.


title: Separating intrinsic interactions from extrinsic correlations in a network of sensory neurons

abstract: Correlations in sensory neural networks have both extrinsic and intrinsic origins. Extrinsic or stimulus correlations arise from shared inputs to the network, and thus depend strongly on the stimulus ensemble. Intrinsic or noise correlations reflect biophysical mechanisms of interactions between neurons, which are expected to be robust to changes of the stimulus ensemble. Despite the importance of this distinction for understanding how sensory networks encode information collectively, no method exists to reliably separate intrinsic interactions from extrinsic correlations in neural activity data, limiting our ability to build predictive models of the network response. In this paper we introduce a general strategy to infer {population models of interacting neurons that collectively encode stimulus information}. The key to disentangling intrinsic from extrinsic correlations is to infer the {couplings between neurons} separately from the encoding model, and to combine the two using corrections calculated in a mean-field approximation. We demonstrate the effectiveness of this approach on retinal recordings. The same coupling network is inferred from responses to radically different stimulus ensembles, showing that these couplings indeed reflect stimulus-independent interactions between neurons. The inferred model predicts accurately the collective response of retinal ganglion cell populations as a function of the stimulus.


title: Bringing Information Centric Networking to Challenged Environments: An overview of the Second Workshop on Future Internet Architecture for Developing Regions

abstract: The 2nd workshop on Future Internet Architecture for Developing Regions (FI4D) took place on January 8th, 2017. The workshop was hosted by the 14th Annual IEEE Consumer Communications $\&$ Networking Conference in Las Vegas, USA. This report summarizes the motivation, discussions and perspectives held during the half-day event.


title: Towards Understanding and Answering Multi-Sentence Recommendation Questions on Tourism

abstract: We introduce the first system towards the novel task of answering complex multisentence recommendation questions in the tourism domain. Our solution uses a pipeline of two modules: question understanding and answering. For question understanding, we define an SQL-like query language that captures the semantic intent of a question; it supports operators like subset, negation, preference and similarity, which are often found in recommendation questions. We train and compare traditional CRFs as well as bidirectional LSTM-based models for converting a question to its semantic representation. We extend these models to a semisupervised setting with partially labeled sequences gathered through crowdsourcing. We find that our best model performs semi-supervised training of BiDiLSTM+CRF with hand-designed features and CCM(Chang et al., 2007) constraints. Finally, in an end to end QA system, our answering component converts our question representation into queries fired on underlying knowledge sources. Our experiments on two different answer corpora demonstrate that our system can significantly outperform baselines with up to 20 pt higher accuracy and 17 pt higher recall.


title: Compactifiable classes of compacta

abstract: We introduce the notion of compactifiable classes -- these are classes of metrizable compact spaces that can be up to homeomorphic copies ``disjointly combined'' into one metrizable compact space. This is witnessed by so-called compact composition of the class. Analogously, we consider Polishable classes and Polish compositions. The question of compactifiability or Polishability of a class is related to hyperspaces. Strongly compactifiable and strongly Polishable classes may be characterized by the existence of a corresponding family in the hyperspace of all metrizable compacta. We systematically study the introduced notions -- we give several characterizations, consider preservation under various constructions, and raise several questions.


title: Role of kinetic instability in runaway electron avalanche and elevated critical electric fields

abstract: The effects of kinetic whistler-wave instabilities on the runaway-electron (RE) avalanche is investigated. With parameters from DIII-D experiments, we show that RE scattering from excited whistler waves can explain several poorly understood experimental results seen in a variety of tokamaks. We find an increase of the avalanche growth rate and threshold electric field, bringing the present model much closer to observations than previous results. The excitation of kinetic instabilities and the scattering of resonant electrons are calculated self-consistently using a quasilinear model. We also explain the observed fast growth of electron cyclotron emission (ECE) signals and excitation of very low-frequency whistler modes observed in the quiescent RE experiments at DIII-D [D. A. Spong et al., submitted to Phys. Rev. Lett.]. These results indicate that by controlling the background thermal plasma temperature, the plasma wave can be excited spontaneously in tokamak disruptions and the avalanche generation of runaway electrons may be suppressed.


title: Shielding Google's language toxicity model against adversarial attacks

abstract: Lack of moderation in online communities enables participants to incur in personal aggression, harassment or cyberbullying, issues that have been accentuated by extremist radicalisation in the contemporary post-truth politics scenario. This kind of hostility is usually expressed by means of toxic language, profanity or abusive statements. Recently Google has developed a machine-learning-based toxicity model in an attempt to assess the hostility of a comment; unfortunately, it has been suggested that said model can be deceived by adversarial attacks that manipulate the text sequence of the comment. In this paper we firstly characterise such adversarial attacks as using obfuscation and polarity transformations. The former deceives by corrupting toxic trigger content with typographic edits, whereas the latter deceives by grammatical negation of the toxic content. Then, we propose a two--stage approach to counter--attack these anomalies, bulding upon a recently proposed text deobfuscation method and the toxicity scoring model. Lastly, we conducted an experiment with approximately 24000 distorted comments, showing how in this way it is feasible to restore toxicity of the adversarial variants, while incurring roughly on a twofold increase in processing time. Even though novel adversary challenges would keep coming up derived from the versatile nature of written language, we anticipate that techniques combining machine learning and text pattern recognition methods, each one targeting different layers of linguistic features, would be needed to achieve robust detection of toxic language, thus fostering aggression--free digital interaction.


title: On the Logic (plus some history and philosophy) of Statistical Tests and Scientific Investigation

abstract: Every scientific endeavour consists of (at least) two components: A hypothesis on the one hand and data on the other. There is always a more or less abstract level - some theory, a set of concepts, certain relations of ideas - and a concrete level, i.e., empirical evidence, experiments or some observations which constitute matters of fact. The focus of this contribution is on elementary models connecting both levels that have been very popular in psychological research - statistical tests. Going from simple to complex we will examine four paradigms of statistical testing (Fisher, Likelihood, Bayes, Neyman & Pearson) and an elegant contemporary treatment. In a nutshell, testing is an easy problem that has a straightforward mathematical solution. However, it is rather surprising that the statistical mainstream has pursued a different line of argument. The application of the latter theory in psychology and other fields has brought some progress but has also impaired scientific thinking.


title: Multipliers in Bessel potential spaces. The case of different sign smooth indices

abstract: The objective of this paper is to describe the space of multipliers acting from a Bessel potential space $H^s_p(\mathbb R^n)$ into another space $H^{-t}_q(\mathbb R^n)$, provided that the smooth indices of these spaces have different signs, i.e. $s, t \geqslant 0$. This space of multipliers consists of distributions $u$, such that for all $\varphi \in H^s_p(\mathbb R^n)$ the product $\varphi \cdot u$ is well-defined and belongs to the space $H^{-t}_q(\mathbb R^n)$. We succeed to describe this space explicitly, provided that $p \leqslant q$ and one of the following conditions $$ s \geqslant t \geqslant 0, \ s > n/p \ \ \, \text{or} \ \ \, t \geqslant s \geqslant 0, \ t > n/q' \quad (\: \text{where} \; 1/q +1/q' = 1), $$ holds. In this case one has $$ M[H^s_p(\mathbb{R}^n) \to H^{-t}_{q}(\mathbb{R}^n)] = H^{-t}_{q, \: unif}(\mathbb{R}^n) \cap H^{-s}_{p', \: unif}(\mathbb{R}^n), $$ where $H^\gamma_{r, \: unif}(\mathbb{R}^n), \: \gamma \in \mathbb{R}, \: r > 1$ is the scale of uniformly localized Bessel potential spaces. In particular but important case $s = t < n/\max (p,q')$ we prove two-sided continuous embeddings $$ H^{-s}_{r_1, \: unif}(\mathbb{R}^n) \subset M[H^s_p(\mathbb{R}^n) \to H^{-s}_q(\mathbb{R}^n)] \subset H^{-s}_{r_2, \: unif}(\mathbb{R}^n), $$ where $r_2 = \max (p', q), \ r_1 =[s/n-(1/p -1/q)]^{-1}$.


title: The Herschel-ATLAS: magnifications and physical sizes of $500\,\mu$m-selected strongly lensed galaxies

abstract: We perform lens modelling and source reconstruction of Submillimeter Array (SMA) data for a sample of 12 strongly lensed galaxies selected at 500$\mu$m in the Herschel Astrophysical Terahertz Large Area Survey H-ATLAS. A previous analysis of the same dataset used a single S\`ersic profile to model the light distribution of each background galaxy. Here we model the source brightness distribution with an adaptive pixel scale scheme, extended to work in the Fourier visibility space of interferometry. We also present new SMA observations for seven other candidate lensed galaxies from the H-ATLAS sample. Our derived lens model parameters are in general consistent with previous findings. However, our estimated magnification factors, ranging from 3 to 10, are lower. The discrepancies are observed in particular where the reconstructed source hints at the presence of multiple knots of emission. We define an effective radius of the reconstructed sources based on the area in the source plane where emission is detected above 5$\sigma$. We also fit the reconstructed source surface brightness with an elliptical Gaussian model. We derive a median value $r_{eff}\,\sim 1.77\,$kpc and a median Gaussian full width at half maximum $\sim1.47\,$kpc. After correction for magnification, our sources have intrinsic star formation rates SFR$\,\sim900-3500\,M_{\odot}yr^{-1}$, resulting in a median star formation rate surface density $\Sigma_{SFR}\sim132\,M_{\odot}$ yr$^{-1}$ kpc$^{-2}$ (or $\sim 218\,M_{\odot}$ yr$^{-1}$ kpc$^{-2}$ for the Gaussian fit). This is consistent with what observed for other star forming galaxies at similar redshifts, and is significantly below the Eddington limit for a radiation pressure regulated starburst.


title: Magnetic-oscillation Mechanism for Understanding Timing Noise of Pulsar Observation

abstract: Highly precise pulsar timing is very important for understanding the nature of a neutron star, and it can even be used to detect gravitational waves. Unfortunately, the accuracy of the pulsar timing is seriously affected by the spin-down irregularities of pulsars.


title: Seymour's conjecture on 2-connected graphs of large pathwidth

abstract: We prove the conjecture of Seymour (1993) that for every apex-forest $H_1$ and outerplanar graph $H_2$ there is an integer $p$ such that every 2-connected graph of pathwidth at least $p$ contains $H_1$ or $H_2$ as a minor. An independent proof was recently obtained by Dang and Thomas.


title: Understanding Type Ia supernovae through their U-band spectra

abstract: Context. Observations of Type Ia supernovae (SNe Ia) can be used to derive accurate cosmological distances through empirical standardization techniques. Despite this success neither the progenitors of SNe Ia nor the explosion process are fully understood. The U-band region has been less well observed for nearby SNe, due to technical challenges, but is the most readily accessible band for high-redshift SNe. Aims. Using spectrophotometry from the Nearby Supernova Factory, we study the origin and extent of U-band spectroscopic variations in SNe Ia and explore consequences for their standardization and the potential for providing new insights into the explosion process. Methods. We divide the U-band spectrum into four wavelength regions {\lambda}(uNi), {\lambda}(uTi), {\lambda}(uSi) and {\lambda}(uCa). Two of these span the Ca H&K {\lambda}{\lambda} 3934, 3969 complex. We employ spectral synthesis using SYNAPPS to associate the two bluer regions with Ni/Co and Ti. Results. (1) The flux of the uTi feature is an extremely sensitive temperature/luminosity indicator, standardizing the SN peak luminosity to 0.116 $\pm$ 0.011 mag RMS. A traditional SALT2.4 fit on the same sample yields a 0.135 mag RMS. Standardization using uTi also reduces the difference in corrected magnitude between SNe originating from different host galaxy environments. (2) Early U-band spectra can be used to probe the Ni+Co distribution in the ejecta, thus offering a rare window into the source of lightcurve power. (3) The uCa flux further improves standardization, yielding a 0.086 $\pm$ 0.010 mag RMS without the need to include an additional intrinsic dispersion to reach {\chi}$^2$/dof $\sim$ 1. This reduction in RMS is partially driven by an improved standardization of Shallow Silicon and 91T-like SNe.


title: Elastic-Viscoplastic Model for Clays: Development, Validation, and Application

abstract: This paper presents an elastic-viscoplastic (EVP) constitutive model in triaxial space and general stress space for isotropic clays. The EVP model is anchored in the bounding surface theory along with the mapping rule and adopts a critical state soil mechanics framework. It incorporates creep effects, and a non-linear creep function is used in the model. The EVP deformation of clay is integrated considering a reference surface and loading surface. An image parameter is deduced to establish the image surface. The strain rate tensor of the model comprises elastic-strain-rate tensor and viscoplastic-strain-rate tensor. The model formulation is capable of accounting for composite as well as single surface ellipses. Parameters of the model can be extracted from conventional oedometer and triaxial tests. The model performance is validated by capturing the behaviours in creep test, relaxation test, strain-rate effect test, and over consolidation ratio effect test of Kaolin clay, Hong Kong Marine Deposit clay, and Fukakusa clay. The model is also implemented in a Finite Element (FE) code and used to predict the long-term performance of the Nerang Broadbeach Roadway embankment constructed in Australia. The long-term settlement prediction of this embankment is also compared with that obtained with the Modified Cam Clay (MCC) model. Pertinent details of the theoretical framework of the proposed EVP model along with its validation, FE implementation and field application are discussed in this paper.


title: Dynamical phase transition in neutron stars

abstract: In this work, we have studied the dynamical evolution of the shock front in a neutron star. The shock wave is expected to possess enough strength to ignite the nuclear matter thereby converting it to quark matter. The conversion of nuclear to quark matter is assumed to take place at the shock discontinuity. The density and pressure discontinuity is studied both spatially and temporally as it starts near the center of the star and moves towards the surface. Polytropic equations of state which mimics real original nuclear matter and quark matter equations of state are used to study such dynamical phase transition. Solving relativistic hydrodynamic equations for a spherically symmetric star we have studied the phase transition assuming a considerable density discontinuity near the center. We find that as the shock wave propagates outwards, its intensity decreases with time, however the shock velocity peaks up and reaches a value close to that of light. Such fast shock velocity indicates rapid phase transition in neutron star taking place on a timescale of some tens of microseconds. Such a result is quite an interesting one, and it differs from previous calculations that the phase transition in neutron stars takes at least some tens of milliseconds. Rapid phase transition can have significant observational significance because such fast phase transition would imply quite strong gravitational wave signals but very short lived. Such short-lived gravitational wave signals would be accompanied with short-lived gamma-ray bursts and neutrino signals originating from the neutrino and gamma-ray generation from the phase transition of nuclear matter to quark matter.


title: The NANOGrav 11-year Data Set: High-precision timing of 45 Millisecond Pulsars

abstract: We present high-precision timing data over time spans of up to 11 years for 45 millisecond pulsars observed as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project, aimed at detecting and characterizing low-frequency gravitational waves. The pulsars were observed with the Arecibo Observatory and/or the Green Bank Telescope at frequencies ranging from 327 MHz to 2.3 GHz. Most pulsars were observed with approximately monthly cadence, with six high--timing-precision pulsars observed weekly, and all were observed at widely separated frequencies at each observing epoch in order to fit for time-variable dispersion delays. We describe our methods for data processing, time-of-arrival (TOA) calculation, and the implementation of a new, automated method for removing outlier TOAs. We fit a timing model for each pulsar that includes spin, astrometric, and, if necessary, binary parameters, in addition to time-variable dispersion delays and parameters that quantify pulse-profile evolution with frequency. The new timing solutions provide three new parallax measurements, two new Shapiro delay measurements, and two new measurements of large orbital-period variations. We fit models that characterize sources of noise for each pulsar. We find that 11 pulsars show significant red noise, with generally smaller spectral indices than typically measured for non-recycled pulsars, possibly suggesting a different origin. Future papers will use these data to constrain or detect the signatures of gravitational-wave signals.


title: A Bramble-Pasciak conjugate gradient method for discrete Stokes equations with random viscosity

abstract: We study the iterative solution of linear systems of equations arising from stochastic Galerkin finite element discretizations of saddle point problems. We focus on the Stokes model with random data parametrized by uniformly distributed random variables and discuss well-posedness of the variational formulations. We introduce a Bramble-Pasciak conjugate gradient method as a linear solver. It builds on a non-standard inner product associated with a block triangular preconditioner. The block triangular structure enables more sophisticated preconditioners than the block diagonal structure usually applied in MINRES methods. We show how the existence requirements of a conjugate gradient method can be met in our setting. We analyze the performance of the solvers depending on relevant physical and numerical parameters by means of eigenvalue estimates. For this purpose, we derive bounds for the eigenvalues of the relevant preconditioned sub-matrices. We illustrate our findings using the flow in a driven cavity as a numerical test case, where the viscosity is given by a truncated Karhunen-Lo\`eve expansion of a random field. In this example, a Bramble-Pasciak conjugate gradient method with block triangular preconditioner outperforms a MINRES method with block diagonal preconditioner in terms of iteration numbers.


title: Characterization of rational matrices that admit finite digit representations

abstract: Let $A$ be an $n \times n$ matrix with rational entries and let \[ \mathbb{Z}^n[A] := \bigcup_{k=1}^{\infty} \left( \mathbb{Z}^n + A\mathbb{Z}^n + \dots + A^{k-1}\mathbb{Z}^n\right) \] be the minimal $A$-invariant $\mathbb{Z}$-module containing the lattice $\mathbb{Z}^n$. If $\mathcal{D}\subset\mathbb{Z}^n[A]$ is a finite set we call the pair $(A,\mathcal{D})$ a digit system. We say that $(A,\mathcal{D})$ has the finiteness property if each $\mathbf{z} \in \mathbb{Z}^n[A]$ can be written in the form \[ \mathbf{z} = \mathbf{d}_0 + A\mathbf{d}_1 + \dots + A^k\mathbf{d}_k, \] with $k\in\mathbb{N}$ and digits $\mathbf{d}_j \in \mathcal{D}$ for $0\le j\le k$. We prove that for a given matrix $A \in M_n(\mathbb{Q})$ there is a finite set $\mathcal{D}\subset\mathbb{Z}^n[A]$ such that $(A, \mathcal{D})$ has the finiteness property if and only if $A$ has no eigenvalue of absolute value $< 1$. This result is the matrix analogue of the height reducing property of algebraic numbers. In proving this result we also characterize integer polynomials $P \in \mathbb{Z}[x]$ that admit digit systems having the finiteness property in the quotient ring $\mathbb{Z}[x]/(P)$.


title: The Simplest Form of the Lorentz Transformations

abstract: We report the simplest possible form to compute rotations around arbitrary axis and boosts in arbitrary directions for 4-vectors (space-time points, energy-momentum) and bi-vectors (electric and magnetic field vectors) by symplectic similarity transformations. The Lorentz transformations are based exclusively on real $4\times 4$-matrices and require neither complex numbers nor special implementations of abstract entities like quaternions or Clifford numbers. No raising or lowering of indices is necessary. It is explained how the Lorentz transformations can be derived from the most simple second order Hamiltonian of general significance. Since this approach exclusively uses the real Clifford algebra $Cl(3,1)$, all calculations are based on real $4\times 4$ matrix algebra.


title: Two-Tier Mapper: a user-independent clustering method for global gene expression analysis based on topology

abstract: There is a growing need for unbiased clustering methods, ideally automated. We have developed a topology-based analysis tool called Two-Tier Mapper (TTMap) to detect subgroups in global gene expression datasets and identify their distinguishing features. First, TTMap discerns and adjusts for highly variable features in the control group and identifies outliers. Second, the deviation of each test sample from the control group in a high-dimensional space is computed and the test samples are clustered in a global and local network using a new topological algorithm based on Mapper. Validation of TTMap on both synthetic and biological datasets shows that it outperforms current clustering methods in sensitivity and stability; clustering is not affected by removal of samples from the control group, choice of normalization nor subselection of data. There is no user induced bias because all parameters are data-driven. Datasets can readily be combined into one analysis. TTMap reveals hitherto undetected gene expression changes in mouse mammary glands related to hormonal changes during the estrous cycle. This illustrates the ability to extract information from highly variable biological samples and its potential for personalized medicine.


title: Modeling of Negative Capacitance of Ferroelectric Capacitors as a Non-Quasi Static Effect

abstract: Pulse-based studies of ferroelectric capacitor systems have been used by several groups to experimentally probe the mechanisms of apparent negative capacitance. In this paper, the behavior of such systems is modeled through SPICE simulation with a delayed-response Preisach model, and the results are compared to available data. It is found that a simple ferroelectric domain delay model can explain much of the observed behavior, capturing the qualitatively different effects of bipolar and unipolar switching, as well as the voltage dependence of said switching. The observed behavior and its modeling suggests that the observed negative capacitance is in fact associated with ferroelectric switching, and its presence is highly sensitive to the switching frequency and other details of the ferroelectric system.


title: Design and Performance Characterization of RADICAL-Pilot on Titan

abstract: Many extreme scale scientific applications have workloads comprised of a large number of individual high-performance tasks. The Pilot abstraction decouples workload specification, resource management, and task execution via job placeholders and late-binding. As such, suitable implementations of the Pilot abstraction can support the collective execution of large number of tasks on supercomputers. We introduce RADICAL-Pilot (RP) as a portable, modular and extensible Python-based Pilot system. We describe RP's design, architecture and implementation. We characterize its performance and show its ability to scalably execute workloads comprised of thousands of MPI tasks on Titan--a DOE leadership class facility. Specifically, we investigate RP's weak (strong) scaling properties up to 131K (65K) cores and 4096 (16384) 32 core tasks. RADICAL-Pilot can be used stand-alone, as well as integrated with other tools as a runtime system.


title: Oscillating Entropy and Spin Precession in the Ensemble of Qubits Interacting with Thermal Systems

abstract: We present a simple model which allows us to explain the physical nature of the oscillating entropy. We consider an ensemble of qubits interacting with thermal two-level systems. The entropy of the qubits oscillates between zero and the value of entropy of the thermal systems. We show that the oscillations of the entropy can be clearly explained by the precession of the real or effective spins of the qubits.


title: Precise measurements of torque in von Karman swirling flow driven by a bladed disk

abstract: Scrupulous measurements and detailed data analysis of the torque in a swirling turbulent flow driven by counter-rotating bladed disks reveals an apparent breaking of the law of similarity. Potentially, such breakdown could arise from several possible factors, including dependence on dimensionless numbers other that $Re$ or velocity coupling to other fields such as temperature. However, careful redesign and calibration of the experiment showed that this unexpected result was due to background errorscaused by minute misalignments which lead to a noisy and irreproducible torque signal at low rotation speeds and prevented correct background subtraction normally ascribed to frictional losses. An important lesson to be learnt is that multiple minute misalignments can nonlinearly couple to the torque signal and provide a dc offset that cannot be removed by averaging. That offset can cause the observed divergence of the friction coefficient C_f from its constant value observed in the turbulent regime. To minimize the friction and misalignments, we significantly modified the experimental setup and carried out the experiment with one bladed disk where the disk, torque meter and motor shaft axes can be aligned with significantly smaller error, close to the torque meter resolution. As a result we made precise measurements with high resolution and sensitivity of the small torques produced for low rotation speeds for several water-glycerin solutions of different viscosities and confirmed the similarity law in a wide range of Re in particular in low viscosity fluids.


title: Search for new physics in events with two soft oppositely charged leptons and missing transverse momentum in proton-proton collisions at $\sqrt{s}=$ 13 TeV

abstract: A search is presented for new physics in events with two low-momentum, oppositely charged leptons (electrons or muons) and missing transverse momentum in proton-proton collisions at a centre-of-mass energy of 13 TeV. The data collected using the CMS detector at the LHC correspond to an integrated luminosity of 35.9 fb$^{-1}$. The observed event yields are consistent with the expectations from the standard model. The results are interpreted in terms of pair production of charginos and neutralinos ($\widetilde{\chi}^\pm_1$ and $\widetilde{\chi}^0_2$) with nearly degenerate masses, as expected in natural supersymmetry models with light higgsinos, as well as in terms of the pair production of top squarks ($\widetilde{\mathrm{t}}$), when the lightest neutralino and the top squark have similar masses. At 95% confidence level, wino-like $\widetilde{\chi}^\pm_1$/$\widetilde{\chi}^0_2$ masses are excluded up to 230 GeV for a mass difference of 20 GeV relative to the lightest neutralino. In the higgsino-like model, masses are excluded up to 168 GeV for the same mass difference. For $\widetilde{\mathrm{t}}$ pair production, top squark masses up to 450 GeV are excluded for a mass difference of 40 GeV relative to the lightest neutralino.


title: Learning Implicit Brain MRI Manifolds with Deep Learning

abstract: An important task in image processing and neuroimaging is to extract quantitative information from the acquired images in order to make observations about the presence of disease or markers of development in populations. Having a lowdimensional manifold of an image allows for easier statistical comparisons between groups and the synthesis of group representatives. Previous studies have sought to identify the best mapping of brain MRI to a low-dimensional manifold, but have been limited by assumptions of explicit similarity measures. In this work, we use deep learning techniques to investigate implicit manifolds of normal brains and generate new, high-quality images. We explore implicit manifolds by addressing the problems of image synthesis and image denoising as important tools in manifold learning. First, we propose the unsupervised synthesis of T1-weighted brain MRI using a Generative Adversarial Network (GAN) by learning from 528 examples of 2D axial slices of brain MRI. Synthesized images were first shown to be unique by performing a crosscorrelation with the training set. Real and synthesized images were then assessed in a blinded manner by two imaging experts providing an image quality score of 1-5. The quality score of the synthetic image showed substantial overlap with that of the real images. Moreover, we use an autoencoder with skip connections for image denoising, showing that the proposed method results in higher PSNR than FSL SUSAN after denoising. This work shows the power of artificial networks to synthesize realistic imaging data, which can be used to improve image processing techniques and provide a quantitative framework to structural changes in the brain.


title: Gamma-ray burst afterglow blast waves

abstract: The various stages of baryonic gamma-ray burst afterglow blast waves are reviewed. These are responsible for the afterglow emission from which much of our understanding of gamma-ray bursts derives. Initially, the blast waves are confined to the dense medium surrounding the burster (stellar envelope or dense wind), giving rise to a jet-cocoon structure. A massive ejecta is released and potentially fed by ongoing energy release from the burster and a forward-reverse shock system is set up between ejecta and ambient density. Ultimately the blast wave spreads sideways and slows down, and the dominant afterglow emission shifts from X-rays down to radio. Over the past years significant progress has been made both observationally and theoretically/numerically in our understanding of these blast waves, unique in the universe due to their often incredibly high initial Lorentz factors of 100-1000. The recent discovery of a short gamma-ray burst counterpart to a gravitational wave detection (GW 170817) brings the promise of a completely new avenue to explore and constrain the dynamics of gamma-ray burst blast waves.


title: Hi-Fi: Hierarchical Feature Integration for Skeleton Detection

abstract: In natural images, the scales (thickness) of object skeletons may dramatically vary among objects and object parts, making object skeleton detection a challenging problem. We present a new convolutional neural network (CNN) architecture by introducing a novel hierarchical feature integration mechanism, named Hi-Fi, to address the skeleton detection problem. The proposed CNN-based approach has a powerful multi-scale feature integration ability that intrinsically captures high-level semantics from deeper layers as well as low-level details from shallower layers. % By hierarchically integrating different CNN feature levels with bidirectional guidance, our approach (1) enables mutual refinement across features of different levels, and (2) possesses the strong ability to capture both rich object context and high-resolution details. Experimental results show that our method significantly outperforms the state-of-the-art methods in terms of effectively fusing features from very different scales, as evidenced by a considerable performance improvement on several benchmarks.


title: Hyperbolic HHS II: Graphs of hierarchically hyperbolic groups

abstract: In this paper we consider a large family of graphs of hierarchically hyperbolic groups (HHG) and show that their fundamental groups admit HHG structures. To do that, we will investigate the notion of hierarchical quasi convexity and show that for a hyperbolic HHS it coincides with the notion of quasi-convexity. The main technical result, for which we expect further applications, is that it is possible to incorporate the HHG structure of a hierarchically hyperbolically embedded subgroup into the HHG structure of the ambient group. This generalizes and provides some additional details to a procedure described in [DHS17].


title: A closer look at the quadruply lensed quasar PSOJ0147: spectroscopic redshifts and microlensing effect

abstract: I present timely spectroscopic follow-up of the newly discovered, quadruply lensed quasar PSOJ0147 from the Pan-STARRS 1 survey. The newly acquired optical spectra with GMOS onboard the Gemini north telescope allow us to pin down the redshifts of both the foreground lensing galaxy and the background lensed quasar to be z=0.572 and 2.341, providing firm basis for cosmography with future high cadence photometric monitoring. I also inspect difference spectra from two of the quasar images, revealing the microlensing effect. Long-term spectroscopic follow-ups will shed lights on the structure of the AGN and its environment.


title: The triply-ionized carbon forest from eBOSS: cosmological correlations with quasars in SDSS-IV DR14

abstract: We present measurements of the cross-correlation of the triply-ionized carbon (CIV) forest with quasars using Sloan Digital Sky Survey Data Release 14. The study exploits a large sample of new quasars from the first two years of observations by the Extended Baryon Oscillation Spectroscopic Survey (eBOSS). The CIV forest is a weaker tracer of large-scale structure than the Ly$\alpha$ forest, but benefits from being accessible at redshifts $z<2$ where the quasar number density from eBOSS is high. Our data sample consists of 287,651 CIV forest quasars in the redshift range $1.4<z<3.5$ and 387,315 tracer quasars with $1.2<z<3.5$. We measure large-scale correlations from CIV absorption occuring in three distinct quasar rest-frame wavelength bands of the spectra referred to as the CIV forest, the SiIV forest and the Ly$\alpha$ forest. From the combined fit to the quasar-CIV cross-correlations for the CIV forest and the SiIV forest, the CIV redshift-space distortion parameter is $\beta_{\rm CIV}=0.27_{\ -0.14}^{\ +0.16}$ and its combination with the CIV linear transmission bias parameter is $b_{\rm CIV}(1+\beta_{\rm CIV})=-0.0183_{\ -0.0014}^{\ +0.0013}$ ($1\sigma$ statistical error) at the mean redshift $z=2.00$. Splitting the sample at $z=2.2$ to constrain the bias evolution with redshift yields the power-law exponent $\gamma=0.60\pm0.63$, indicating a significantly weaker redshift-evolution than for the Ly$\alpha$ forest linear transmission bias. We demonstrate that CIV absorption has the potential to be used as a probe of baryon acoustic oscillations (BAO). While the current data set is insufficient for a detection of the BAO peak feature, the final quasar samples for redshifts $1.4<z<2.2$ from eBOSS and the Dark Energy Spectroscopic Instrument (DESI) are expected to provide measurements of the isotropic BAO scale to $\sim7\%$ and $\sim3\%$ precision, respectively, at $z\simeq1.6$.


title: Maximum entropy models reveal the excitatory and inhibitory correlation structures in cortical neuronal activity

abstract: Maximum Entropy models can be inferred from large data-sets to uncover how collective dynamics emerge from local interactions. Here, such models are employed to investigate neurons recorded by multielectrode arrays in the human and monkey cortex. Taking advantage of the separation of excitatory and inhibitory neuron types, we construct a model including this distinction. This approach allows to shed light upon differences between excitatory and inhibitory activity across different brain states such as wakefulness and deep sleep, in agreement with previous findings. Additionally, Maximum Entropy models can also unveil novel features of neuronal interactions, which are found to be dominated by pairwise interactions during wakefulness, but are population-wide during deep sleep. In particular, inhibitory neurons are observed to be strongly tuned to the inhibitory population. Overall, we demonstrate Maximum Entropy models can be useful to analyze data-sets with classified neuron types, and to reveal the respective roles of excitatory and inhibitory neurons in organizing coherent dynamics in the cerebral cortex.


title: The IceCube Neutrino Observatory, the Pierre Auger Observatory and the Telescope Array: Joint Contribution to the 35th International Cosmic Ray Conference (ICRC 2017)

abstract: Joint contributions of the IceCube Collaboration, the Telescope Array Collaboration, and the Pierre Auger Collaboration to the 35th International Cosmic Ray Conference (ICRC 2017), 12-20 July 2017, Bexco, Busan, Korea.


title: Supercurrent through a spin-split quasi-ballistic point contact in an InAs nanowire

abstract: We study the superconducting proximity effect in an InAs nanowire contacted by Ta-based superconducting electrodes. Using local bottom gates, we control the potential landscape along the nanowire, tuning its conductance to a quasi-ballistic regime. At high magnetic field ($B$), we observe approximately quantized conductance plateaus associated with the first two spin-polarized one-dimensional modes. For $B < 1$ T, the onset of superconductivity occurs in concomitance with the development of sizeable charge localization leading to a 0.7-type conductance anomaly. In this regime, the proximity supercurrent exhibits an unusual, non-monotonic $B$ dependence. We interpret this finding in terms of a competition between the Kondo effect, dominating near $B=0$, and the Zeeman effect, enforcing spin polarization and the emergence of a $\pi$ phase shift in the Josephson relation at higher $B$.


title: Spatial Spread of Epidemic Diseases in Geographical Settings: Seasonal Influenza Epidemics in Puerto Rico

abstract: Deterministic models are developed for the spatial spread of epidemic diseases in geographical settings. The models are focused on outbreaks that arise from a small number of infected hosts imported into sub-regions of the geographical settings. The goal is to understand how spatial heterogeneity influences the transmission dynamics of the susceptible and infected populations. The models consist of systems of partial differential equations with diffusion terms describing the spatial spread of the underlying microbial infectious agents. The model is compared with real data from seasonal influenza epidemics in Puerto Rico.


title: On the Kahler-Einstein metric at strictly pseudoconvex points

abstract: We prove a local boundary regularity result for the complete Kahler-Einstein metrics of negative Ricci curvature near strictly pseudoconvex boundary point. We also study the asymptotic behaviour of their holomorphic bisectional curvatures near such points.


title: Special moves for open book decompositions of 3-manifolds

abstract: We provide a complete set of two moves that suffice to relate any two open book decompositions of a given 3-manifold. One of the moves is the usual plumbing with a positive or negative Hopf band, while the other one is a special local version of Harer's twisting, which is presented in two different (but stably equivalent) forms. Our approach relies on 4-dimensional Lefschetz fibrations, and on 3-dimensional contact topology, via the Giroux-Goodman stable equivalence theorem for open book decompositions representing homologous contact structures.


title: The Arbitrarily Varying Relay Channel

abstract: We study the arbitrarily varying relay channel, and establish the cutset bound and partial decode-forward bound on the random code capacity. We further determine the random code capacity for special cases. Then, we consider conditions under which the deterministic code capacity is determined as well.


title: Controllability of the Navier-Stokes equation in a rectangle with a little help of a distributed phantom force

abstract: We consider the 2D incompressible Navier-Stokes equation in a rectangle with the usual no-slip boundary condition prescribed on the upper and lower boundaries. We prove that for any positive time, for any finite energy initial data, there exist controls on the left and right boundaries and a distributed force, which can be chosen arbitrarily small in any Sobolev norm in space, such that the corresponding solution is at rest at the given final time. Our work improves earlier results where the distributed force is small only in a negative Sobolev space. It is a further step towards an answer to Jacques-Louis Lions' question about the small-time global exact boundary controllability of the Navier-Stokes equation with the no-slip boundary condition, for which no distributed force is allowed. Our analysis relies on the well-prepared dissipation method already used for Burgers and for Navier-Stokes in the case of the Navier slip-with-friction boundary condition. In order to handle the larger boundary layers associated with the no-slip boundary condition, we perform a preliminary regularization into analytic functions with arbitrarily large analytic radius and prove a long-time nonlinear Cauchy-Kovalevskaya estimate relying only on horizontal analyticity.


title: Complex Reaction Kinetics in Chemistry: A unified picture suggested by Mechanics in Physics

abstract: Complex biochemical pathways or regulatory enzyme kinetics can be reduced to chains of elementary reactions, which can be described in terms of chemical kinetics. This discipline provides a set of tools for quantifying and understanding the dialogue between reactants, whose framing into a solid and consistent mathematical description is of pivotal importance in the growing field of biotechnology. Among the elementary reactions so far extensively investigated, we recall the socalled Michaelis-Menten scheme and the Hill positive-cooperative kinetics, which apply to molecular binding and are characterized by the absence and the presence, respectively, of cooperative interactions between binding sites, giving rise to qualitative different phenomenologies. However, there is evidence of reactions displaying a more complex, and by far less understood, pattern: these follow the positive-cooperative scenario at small substrate concentration, yet negative-cooperative effects emerge and get stronger as the substrate concentration is increased. In this paper we analyze the structural analogy between the mathematical backbone of (classical) reaction kinetics in Chemistry and that of (classical) mechanics in Physics: techniques and results from the latter shall be used to infer properties on the former.


title: High-order symbolic strong-coupling expansion for the Bose-Hubbard model

abstract: Combining the process-chain method with a symbolized evaluation we work out in detail a high-order symbolic strong-coupling expansion (HSSCE) for determining the quantum phase boundaries between the Mott insulator and the superfluid phase of the Bose-Hubbard model for different fillings in hypercubic lattices of different dimensions. With a subsequent Pad{\'e} approximation we achieve for the quantum phase boundaries a high accuracy, which is comparable to high-precision quantum Monte-Carlo simulations, and show that all the Mott lobes can be rescaled to a single one. As a further cross-check, we find that the HSSCE results coincide with a hopping expansion of the quantum phase boundaries, which follow from the effective potential Landau theory (EPLT).


title: Simulating quantum dynamical phenomena using classical oscillators: Landau-Zener-Stuckelberg-Majorana interferometry, latching modulation, and motional averaging

abstract: A quantum system can be driven by either sinusoidal, rectangular, or noisy signals. In the literature, these regimes are referred to as Landau-Zener-Stuckelberg-Majorana (LZSM) interferometry, latching modulation, and motional averaging, respectively. We demonstrate that these pronounced and interesting effects are also inherent in the dynamics of classical two-state systems. We discuss how such classical systems are realized using either mechanical, electrical, or optical resonators. In addition to the fundamental interest of such dynamical phenomena linking classical and quantum physics, we believe that these are attractive for the classical analogue simulation of quantum systems.


title: An asymptotic bound for Castelnuovo-Mumford regularity of certain Ext modules over graded complete intersection rings

abstract: Set $ A := Q/({\bf z}) $, where $ Q $ is a polynomial ring over a field, and $ {\bf z} = z_1,\ldots,z_c $ is a homogeneous $ Q $-regular sequence. Let $ M $ and $ N $ be finitely generated graded $ A $-modules, and $ I $ be a homogeneous ideal of $ A $. We show that (1) $ \mathrm{reg}\left( \mathrm{Ext}_A^{i}(M, I^nN) \right) \le \rho_N(I) \cdot n - f \cdot \left\lfloor \frac{i}{2} \right\rfloor + b \mbox{ for all } i, n \ge 0 $, (2) $ \mathrm{reg}\left( \mathrm{Ext}_A^{i}(M,N/I^nN) \right) \le \rho_N(I) \cdot n - f \cdot \left\lfloor \frac{i}{2} \right\rfloor + b' \mbox{ for all } i, n \ge 0 $, where $ b $ and $ b' $ are some constants, $ f := \mathrm{min}\{ \mathrm{deg}(z_j) : 1 \le j \le c \} $, and $ \rho_N(I) $ is an invariant defined in terms of reduction ideals of $ I $ with respect to $ N $. There are explicit examples which show that these inequalities are sharp.


title: Gravitational fermion creation during an anisotropic phase of cosmological expansion

abstract: The free Dirac equation is solved in a Bianchi Type I space-time, which represents a homogeneous but anisotropic universe, to show the creation of fermionic particles. It is found that unlike in the isotropic case, massless fermion production is possible. An estimate of the energy density of massless particles created during an early anisotropic phase of cosmological expansion is shown to cause substantial back-reaction on the gravitational field. The potential relevance to dark matter particle production, primordial magnetogenesis, and early universe cosmology is discussed briefly.


title: The Reeb Graph Edit Distance is Universal

abstract: We consider the setting of Reeb graphs of piecewise linear functions and study distances between them that are stable, meaning that functions which are similar in the supremum norm ought to have similar Reeb graphs. We define an edit distance for Reeb graphs and prove that it is stable and universal, meaning that it provides an upper bound to any other stable distance. In contrast, via a specific construction, we show that the interleaving distance and the functional distortion distance on Reeb graphs are not universal.


title: Limit-Periodic Schr\"odinger Operators With Lipschitz Continuous IDS

abstract: We show that there exist limit-periodic Schr\"odinger operators such that the associated integrated density of states is Lipschitz continuous. These operators arise in the inverse spectral theoretic KAM approach of P\"oschel.


title: Multiplicity results and qualitative properties for semilinear elliptic problems with Neumann boundary conditions

abstract: In this paper we study multiplicity and qualitative behavior of solutions for semilinear elliptic problems with neumann boundary condition and asymptotically linear smooth nonlinearity. We provide sufficient conditions on the number of eigenvalues the derivative of the nonlinearity crosses to guarantee existence of at least five nontrivial solutions. The techniques we use are a combination of minimization, Leray-Schauder degree, Morse Theory and Reduction method a la Castro-Lazer.


title: The dynamical structure of political corruption networks

abstract: Corruptive behaviour in politics limits economic growth, embezzles public funds, and promotes socio-economic inequality in modern democracies. We analyse well-documented political corruption scandals in Brazil over the past 27 years, focusing on the dynamical structure of networks where two individuals are connected if they were involved in the same scandal. Our research reveals that corruption runs in small groups that rarely comprise more than eight people, in networks that have hubs and a modular structure that encompasses more than one corruption scandal. We observe abrupt changes in the size of the largest connected component and in the degree distribution, which are due to the coalescence of different modules when new scandals come to light or when governments change. We show further that the dynamical structure of political corruption networks can be used for successfully predicting partners in future scandals. We discuss the important role of network science in detecting and mitigating political corruption.


title: Collective motion in triaxial nuclei within minimal length concept

abstract: The concept of minimal length, inspired by Heisenberg algebra, is applied to the geometrical collective Bohr- Mottelson model (BMM) of nuclei. With the deformed canonical commutation relation and the Pauli-Podolsky prescription, we have derived the quantized Hamiltonian operator for triaxial nuclei as we have previously done for axial prolate $\gamma$-rigid ones (M. Chabab et al., Phys. Lett. B 758 (2016) 212-218). By considering an infinite square well like potential in $\beta$ collective shape variable, the eigenvalues of the Hamiltonian are obtained in terms of zeros of Bessel functions of irrational order with an explicit dependence on the minimal length parameter. Moreover, the associated symmetry with the model that we have constructed here can be considered as a new quasi-dynamical critical point symmetries (CPSs) in nuclear structure. The theoretical results indicate a dramatic contribution (Low effect) of the minimal length to energy levels for lower values of the angular momentum and regular (significant effect) for higher values. In fact, these features show that the minimal length scenario which is useful in recognizing the properties of real deformed nuclei having high spins or strong rotations. Finally, numerical calculations are performed for some nuclei such as ${}^{124,128,130}$Xe and ${}^{114}$Pd revealing a qualitative agreement with the experimental data.


title: Exotic energy injection with ExoCLASS: Application to the Higgs portal model and evaporating black holes

abstract: We devise a new user-friendly tool interfaced with the Boltzmann code CLASS to deal with any kind of exotic electromagnetic energy injection in the universe and its impact on anisotropies of the Cosmic Microwave Background. It makes use of the results from standard electromagnetic cascade calculations develop in the context of WIMP annihilation, generalized to incorporate any injection history. We first validate it on a specific WIMP scenario, the Higgs Portal model, confirming that the standard effective on-the-spot treatment is accurate enough. We then analyze the more involved example of evaporating Primordial Black Holes (PBHs) with masses in the range $[3\times10^{13},5\times10^{16}]$g, for which the standard approximations break down. We derive robust CMB bounds on the relic density of evaporating PBHs, ruling out the possibility for PBHs with a monochromatic distribution of masses in the range $[3\times10^{13},2.5\times10^{16}]$g to represent all of the Dark Matter in our Universe. Remarkably, we confirm with an accurate study that the CMB bounds are several orders of magnitude stronger than those from the galactic gamma-ray background in the range $[3\times10^{13},3\times10^{14}]$g. A future CMB experiment like CORE+, or an experiment attempting at measuring the 21 cm signal from the Dark Ages could greatly improve the sensitivity to these models.


title: Prospects of relativistic astrophysics: Viewpoints on the most influential journey

abstract: Based on several recent developments in relativistic astrophysics and observational cosmology, this article discusses prospects of (1) theoretical structure of General Relativity (GR) in the new era and (2) submission of GR to new and previously unanticipated tests. It explores the internal theoretical challenges experienced by the theory and challenges produced by it in the otherwise successful journey of 100 years. Gravitational collapse, geometrization of laws of physics, absence of strong field astronomical tests, the dark components of the universe (Cold Dark Matter (CDM) + Cosmological Constant ($\Lambda$ ) or something similar) and cosmological inflation have not only created serious challenges for the theory but also have ignited contemplation for alternatives and attempts for constraining them through observations. The article first explains the nature of the theory and how it was perceived by scientific community a century ago. It then highlights certain perspectives from cosmology and quantum gravity, the ultimate applications of the theory. With the advent of modern observing techniques, the theory has been experiencing tug of war with several tests. Our true journey with the theory has started only recently, thanks to the development of gravitational wave astronomy and extremely large telescopes some of which are and will be dedicated to test the theory by exploring previously uncharted regions of the universe. Incorporating several aspects of observational astronomy and cosmology, gravitational wave research and theoretical physics the article presents the new challenges created by the theory and also their prospects.


title: When is $X\times Y$ homeomorphic to $X\times_l Y$?

abstract: We identify a class of linearly ordered topological spaces $X$ that may satisfy the property that $X\times X$ is homeomorphic to $X\times_l X$ or can be embedded into a linearly ordered space with the stated property. We justify the conjectures by partial results.


title: RobustGaSP: Robust Gaussian Stochastic Process Emulation in R

abstract: Gaussian stochastic process emulation is a powerful tool for approximating computationally intensive computer models. However, estimation of parameters in the GaSP emulator is a challenging task. No closed-form estimator is available and many numerical problems arise with standard estimates, e.g., the maximum likelihood estimator. In this package, we implement a marginal posterior mode estimator, for special priors and parameterizations, an estimation method that meets the robust parameter estimation criteria discussed in \cite{Gu2018robustness}; mathematical reasons are provided therein to explain why robust parameter estimation can greatly improve predictive performance of the emulator. In addition, inert inputs (inputs that almost have no effect on the variability of a function) can be identified from the marginal posterior mode estimation, at no extra computational cost. The package also implements the parallel partial Gaussian stochastic process (PP GaSP) emulator (\cite{gu2016parallel}) for the scenario where the computer model has multiple outputs on e.g. spatial-temporal coordinates. The package can be operated in a default mode, but also allows numerous user specifications, such as the capability of specifying trend functions and noise terms. Examples are studied herein to highlight the performance of the package in terms of out-of-sample prediction.}


title: Near Optimal Coded Data Shuffling for Distributed Learning

abstract: Data shuffling between distributed cluster of nodes is one of the critical steps in implementing large-scale learning algorithms. Randomly shuffling the data-set among a cluster of workers allows different nodes to obtain fresh data assignments at each learning epoch. This process has been shown to provide improvements in the learning process. However, the statistical benefits of distributed data shuffling come at the cost of extra communication overhead from the master node to worker nodes, and can act as one of the major bottlenecks in the overall time for computation. There has been significant recent interest in devising approaches to minimize this communication overhead. One approach is to provision for extra storage at the computing nodes. The other emerging approach is to leverage coded communication to minimize the overall communication overhead. The focus of this work is to understand the fundamental trade-off between the amount of storage and the communication overhead for distributed data shuffling. In this work, we first present an information theoretic formulation for the data shuffling problem, accounting for the underlying problem parameters (number of workers, $K$, number of data points, $N$, and the available storage, $S$ per node). We then present an information theoretic lower bound on the communication overhead for data shuffling as a function of these parameters. We next present a novel coded communication scheme and show that the resulting communication overhead of the proposed scheme is within a multiplicative factor of at most $\frac{K}{K-1}$ from the information-theoretic lower bound. Furthermore, we present the aligned coded shuffling scheme for some storage values, which achieves the optimal storage vs communication trade-off for $K<5$, and further reduces the maximum multiplicative gap down to $\frac{K-\frac{1}{3}}{K-1}$, for $K\geq 5$.


title: Root Mean Square Minimum Distance as a Quality Metric for Localization Nanoscopy Images

abstract: A localization algorithm in stochastic optical localization nanoscopy plays an important role in obtaining a high-quality image. A universal and objective metric is crucial and necessary to evaluate qualities of nanoscopy images and performances of localization algorithms. In this paper, we propose root mean square minimum distance (RMSMD) as a quality metric for localization nanoscopy images. RMSMD measures an average, local, and mutual fitness between two sets of points. Its properties common to a distance metric as well as unique to itself are presented. The ambiguity, discontinuity, and inappropriateness of the metrics of accuracy, precision, recall, and Jaccard index, which are currently used in the literature, are analyzed. A numerical example demonstrates the advantages of RMSMD over the four existing metrics that fail to distinguish qualities of different nanoscopy images in certain conditions. The unbiased Gaussian estimator that achieves the Fisher information and Cramer-Rao lower bound (CRLB) of a single data frame is proposed to benchmark the quality of localization nanoscopy images and the performance of localization algorithms. The information-achieving estimator is simulated in an example and the result demonstrates the superior sensitivity of RMSMD over the other four metrics. As a universal and objective metric, RMSMD can be broadly employed in various applications to measure the mutual fitness of two sets of points.


title: Prediction of flow dynamics using point processes

abstract: Describing a time series parsimoniously is the first step to study the underlying dynamics. For a time-discrete system, a generating partition provides a compact description such that a time series and a symbolic sequence are one-to-one. But, for a time-continuous system, such a compact description does not have a solid basis. Here, we propose to describe a time-continuous time series using a local cross section and the times when the orbit crosses the local cross section. We show that if such a series of crossing times and some past observations are given, we can predict the system's dynamics with fine accuracy. This reconstructability neither depends strongly on the size nor the placement of the local cross section if we have a sufficiently long database. We demonstrate the proposed method using the Lorenz model as well as the actual measurement of wind speed.


title: A rank 2 Dijkgraaf-Moore-Verlinde-Verlinde formula

abstract: We conjecture a formula for the virtual elliptic genera of moduli spaces of rank 2 sheaves on minimal surfaces $S$ of general type. We express our conjecture in terms of the Igusa cusp form $\chi_{10}$ and Borcherds type lifts of three quasi-Jacobi forms which are all related to the Weierstrass elliptic function. We also conjecture that the generating function of virtual cobordism classes of these moduli spaces depends only on $\chi(\mathcal{O}_S)$ and $K_S^2$ via two universal functions, one of which is determined by the cobordism classes of Hilbert schemes of points on $K3$. We present generalizations of these conjectures, e.g. to arbitrary surfaces with $p_g>0$ and $b_1=0$. We use a result of J. Shen to express the virtual cobordism class in terms of descendent Donaldson invariants. In a prequel we used T. Mochizuki's formula, universality, and toric calculations to compute such Donaldson invariants in the setting of virtual $\chi_y$-genera. Similar techniques allow us to verify our new conjectures in many cases.


title: Linear-time general decoding algorithm for the surface code

abstract: A quantum error correcting protocol can be substantially improved by taking into account features of the physical noise process. We present an efficient decoder for the surface code which can account for general noise features, including coherences and correlations. We demonstrate that the decoder significantly outperforms the conventional matching algorithm on a variety of noise models, including non-Pauli noise and spatially correlated noise. The algorithm is based on an approximate calculation of the logical channel using a tensor-network description of the noisy state.


title: Suzaku Observations of Heavily Obscured (Compton-thick) Active Galactic Nuclei selected by Swift/BAT Hard X-ray Survey

abstract: We present a uniform broadband X-ray (0.5-100.0 keV) spectral analysis of 12 Swift/Burst Alert Telescope (BAT) selected Compton-thick ($\log N_{\mathrm{H}}/\mathrm{cm}^{-2} \geq 24$) Active Galactic Nuclei (CTAGNs) observed with Suzaku. The Suzaku data of 3 objects are published here for the first time. We fit the Suzaku and Swift spectra with models utilizing an analytic reflection code and those utilizing the Monte Carlo based model from an AGN torus by Ikeda et al. 2009. The main results are as follows. (1) The estimated intrinsic luminosity of a CTAGN strongly depends on the model; applying Compton scattering to the transmitted component in an analytic model may largely overestimates the intrinsic luminosity at large column densities. (2) Unabsorbed reflection components are commonly observed, suggesting that the tori are clumpy. (3) Most of CTAGNs show small scattering fractions (<0.5%) implying a buried AGN nature. (4) Comparison with the results obtained for Compton-thin AGNs (Kawamuro et al. 2016) suggests that the properties of these CTAGNs can be understood as a smooth extension from Compton-thin AGNs with heavier obscuration; we find no evidence that the bulk of the population of hard X-ray selected CTAGN is different from less obscured objects.


title: The NuSTAR Extragalactic Surveys: source catalog and the Compton-thick fraction in the UDS field

abstract: We present the results and the source catalog of the NuSTAR survey in the UKIDSS Ultra Deep Survey (UDS) field, bridging the gap in depth and area between NuSTAR's ECDFS and COSMOS surveys. The survey covers a $\sim 0.6$ deg$^2$ area of the field for a total observing time of $\sim$ 1.75 Ms, to a half-area depth of $\sim$ 155 ks corrected for vignetting at $3-24$ keV, and reaching sensitivity limits at half-area in the full ($3-24$ keV), soft ($3-8$ keV) and hard ($8-24$ keV) bands of $2.2 \times 10^{-14}$ erg cm$^{-2}$ s$^{-1}$, $1.0 \times 10^{-14}$ erg cm$^{-2}$ s$^{-1}$, and $2.7 \times 10^{-14}$ erg cm$^{-2}$ s$^{-1}$, respectively. A total of 67 sources are detected in at least one of the three bands, 56 of which have a robust optical redshift with a median of $\langle z\rangle \sim 1.1$. Through a broadband ($0.5-24$ keV) spectral analysis of the whole sample combined with the NuSTAR hardness ratios, we compute the observed Compton-thick (CT; $N_{\rm H} > 10^{24}$ cm$^{-2}$) fraction. Taking into account the uncertainties on each $N_{\rm H}$ measurement, the final number of CT sources is $6.8\pm1.2$. This corresponds to an observed CT fraction of $(11.5\pm2.0)\%$, providing a robust lower limit to the intrinsic fraction of CT AGN and placing constraints on cosmic X-ray background synthesis models.


title: The 105 month Swift-BAT all-sky hard X-ray survey

abstract: We present a catalog of hard X-ray sources detected in the first 105 months of observations with the Burst Alert Telescope (BAT) coded mask imager on board the Swift observatory. The 105 month Swift-BAT survey is a uniform hard X-ray all-sky survey with a sensitivity of $8.40\times 10^{-12}\ {\rm erg\ s^{-1}\ cm^{-2}}$ over 90% of the sky and $7.24\times 10^{-12}\ {\rm erg\ s^{-1}\ cm^{-2}}$ over 50% of the sky in the 14-195 keV band. The Swift-BAT 105 month catalog provides 1632 (422 new detections) hard X-ray sources in the 14-195 keV band above the 4.8{\sigma} significance level. Adding to the previously known hard X-ray sources, 34% (144/422) of the new detections are identified as Seyfert AGN in nearby galaxies (z<0.2). The majority of the remaining identified sources are X-ray binaries (7%, 31) and blazars/BL Lac objects (10%, 43). As part of this new edition of the Swift-BAT catalog, we release eight-channel spectra and monthly sampled light curves for each object in the online journal and at the Swift-BAT 105 month Web site.


title: Evidence for radial variations in the stellar mass-to-light ratio of massive galaxies from weak and strong lensing

abstract: The Initial Mass Function (IMF) for massive galaxies can be constrained by combining stellar dynamics with strong gravitational lensing. However, this method is limited by degeneracies between the density profile of dark matter and the stellar mass-to-light ratio. In this work we reduce this degeneracy by combining weak lensing together with strong lensing and stellar kinematics. Our analysis is based on two galaxy samples: 45 strong lenses from the SLACS survey and 1,700 massive quiescent galaxies from the SDSS main spectroscopic sample with weak lensing measurements from the Hyper Suprime-Cam survey. We use a Bayesian hierarchical approach to jointly model all three observables. We fit the data with models of varying complexity and show that a model with a radial gradient in the stellar mass-to-light ratio is required to simultaneously describe both galaxy samples. This result is driven by a subset of strong lenses with very steep total density profile, that cannot be fitted by models with no gradient. Our measurements are unable to determine whether $M_*/L$ gradients are due to variations in stellar population parameters at fixed IMF, or to gradients in the IMF itself. The inclusion of $M_*/L$ gradients decreases dramatically the inferred IMF normalisation, compared to previous lensing-based studies, with the exact value depending on the assumed dark matter profile. The main effect of strong lensing selection is to shift the stellar mass distribution towards the high mass end, while the halo mass and stellar IMF distribution at fixed stellar mass are not significantly affected.


title: Unsupervised Low-Dimensional Vector Representations for Words, Phrases and Text that are Transparent, Scalable, and produce Similarity Metrics that are Complementary to Neural Embeddings

abstract: Neural embeddings are a popular set of methods for representing words, phrases or text as a low dimensional vector (typically 50-500 dimensions). However, it is difficult to interpret these dimensions in a meaningful manner, and creating neural embeddings requires extensive training and tuning of multiple parameters and hyperparameters. We present here a simple unsupervised method for representing words, phrases or text as a low dimensional vector, in which the meaning and relative importance of dimensions is transparent to inspection. We have created a near-comprehensive vector representation of words, and selected bigrams, trigrams and abbreviations, using the set of titles and abstracts in PubMed as a corpus. This vector is used to create several novel implicit word-word and text-text similarity metrics. The implicit word-word similarity metrics correlate well with human judgement of word pair similarity and relatedness, and outperform or equal all other reported methods on a variety of biomedical benchmarks, including several implementations of neural embeddings trained on PubMed corpora. Our implicit word-word metrics capture different aspects of word-word relatedness than word2vec-based metrics and are only partially correlated (rho = ~0.5-0.8 depending on task and corpus). The vector representations of words, bigrams, trigrams, abbreviations, and PubMed title+abstracts are all publicly available from http://arrowsmith.psych.uic.edu for release under CC-BY-NC license. Several public web query interfaces are also available at the same site, including one which allows the user to specify a given word and view its most closely related terms according to direct co-occurrence as well as different implicit similarity metrics.