title: Incorporating Views on Marginal Distributions in the Calibration of Risk Models

abstract: Entropy based ideas find wide-ranging applications in finance for calibrating models of portfolio risk as well as options pricing. The abstracted problem, extensively studied in the literature, corresponds to finding a probability measure that minimizes relative entropy with respect to a specified measure while satisfying constraints on moments of associated random variables. These moments may correspond to views held by experts in the portfolio risk setting and to market prices of liquid options for options pricing models. However, it is reasonable that in the former settings, the experts may have views on tails of risks of some securities. Similarly, in options pricing, significant literature focuses on arriving at the implied risk neutral density of benchmark instruments through observed market prices. With the intent of calibrating models to these more general stipulations, we develop a unified entropy based methodology to allow constraints on both moments as well as marginal distributions of functions of underlying securities. This is applied to Markowitz portfolio framework, where a view that a particular portfolio incurs heavy tailed losses is shown to lead to fatter and more reasonable tails for losses of component securities. We also use this methodology to price non-traded options using market information such as observed option prices and implied risk neutral densities of benchmark instruments.


title: Search for exotics at $BABAR$

abstract: One of the most intriguing puzzles in hadron spectroscopy are the numerous charmonium-like states observed in the last decade, including charged states that are manifestly exotic. The $BABAR$ experiment has extensively studied those in B meson decays, initial state radiation processes and two photon reactions. The study of the process $B \rightarrow J/\psi \phi K$, with a search for the resonant states X(4140) and X(4270) in their decays to $J/\psi \phi$, will be highlighted. The recent results of the Dalitz analysis of $\eta_c$ to 3 pseudoscalar mesons, via 2-photon interactions, will be presented in this report


title: Discovery of the VHE gamma-ray source HESS J1832-093 in the vicinity of SNR G22.7-0.2

abstract: The region around the supernova remnant (SNR) W41 contains several TeV sources and has prompted the H.E.S.S. Collaboration to perform deep observations of this field of view. This resulted in the discovery of the new very high energy (VHE) source HESS J1832-093, at the position $\rm RA=18^h 32^m 50^s \pm 3^s_{stat} \pm 2^s_{syst}, \rm Dec=-9^\circ 22' 36'' \pm 32''_{stat} \pm 20''_{syst} (J2000)$, spatially coincident with a part of the radio shell of the neighboring remnant G22.7-0.2. The photon spectrum is well described by a power-law of index $\Gamma = 2.6 \pm 0.3_{\rm stat} \pm 0.1_{\rm syst}$ and a normalization at 1 TeV of $\Phi_0=(4.8 \pm 0.8_{\rm stat}\pm 1.0_{\rm syst})\,\times\,10^{-13}\,\rm{cm} ^{-2}\,s^{-1}\,TeV^{-1}$. The location of the gamma-ray emission on the edge of the SNR rim first suggested a signature of escaping cosmic-rays illuminating a nearby molecular cloud. Then a dedicated XMM-Newton observation led to the discovery of a new X-ray point source spatially coincident with the TeV excess. Two other scenarios were hence proposed to identify the nature of HESS J1832-093. Gamma-rays from inverse Compton radiation in the framework of a pulsar wind nebula scenario or the possibility of gamma-ray production within a binary system are therefore also considered. Deeper multi-wavelength observations will help to shed new light on this intriguing VHE source.


abstract: Book review of Stochastic Chemical Kinetics. Theory and (Mostly) Systems Biological Applications, P. Erdi, G. Lente. Springer (2014)


title: Constraining neutrino electromagnetic properties by germanium detectors

abstract: The electromagnetic properties of neutrinos, which are either trivial or negligible in the context of the Standard Model, can probe new physics and have significant implications in astrophysics and cosmology. The current best direct limits on the neutrino millicharges and magnetic moments are both derived from data taken with germanium detectors with low thresholds at keV levels. In this paper, we discuss in detail a robust, ab initio method: the multiconfiguration relativistic random phase approximation, that enables us to reliably understand the germanium detector response at the sub-keV level, where atomic many-body physics matters. Using existing data with sub-keV thresholds, limits on reactor antineutrino's millicharge, magnetic moment, and charge radius squared are derived. The projected sensitivities for next generation experiments are also given and discussed.


title: KK-monopoles and G-structures in M-theory/type IIA reductions

abstract: We argue that M-theory/massive IIA backgrounds including KK-monopoles are suitably described in the language of G-structures and their intrinsic torsion. To this end, we study classes of minimal supergravity models that admit an interpretation as twisted reductions in which the twist parameters are not restricted to satisfy the Jacobi constraints $\omega\, \omega=0$ required by an ordinary Scherk-Schwarz reduction. We first derive the correspondence between four-dimensional data and torsion classes of the internal space and, then, check the one-to-one correspondence between higher-dimensional and four-dimensional equations of motion. Remarkably, the whole construction holds regardless of the Jacobi constraints, thus shedding light upon the string/M-theory interpretation of (smeared) KK-monopoles.


title: Ground states and concentration phenomena for the fractional Schr\"odinger equation

abstract: We consider here solutions of the nonlinear fractional Schr\"odinger equation $$\epsilon^{2s}(-\Delta)^s u+V(x)u=u^p.$$ We show that concentration points must be critical points for $V$. We also prove that, if the potential $V$ is coercive and has a unique global minimum, then ground states concentrate suitably at such minimal point as $\epsilon$ tends to zero. In addition, if the potential $V$ is radial, then the minimizer is unique provided $\epsilon$ is small.


title: The algebraic structure of quantum partial isometries

abstract: The partial isometries of $\mathbb R^N,\mathbb C^N$ form compact semigroups $\widetilde{O}_N,\widetilde{U}_N$. We discuss here the liberation question for these semigroups, and for their discrete versions $\widetilde{H}_N,\widetilde{K}_N$. Our main results concern the construction of half-liberations $\widetilde{H}_N^\times,\widetilde{K}_N^\times,\widetilde{O}_N^\times,\widetilde{U}_N^\times$ and of liberations $\widetilde{H}_N^+,\widetilde{K}_N^+,\widetilde{O}_N^+,\widetilde{U}_N^+$. We include a detailed algebraic and probabilistic study of all these objects, justifying our "half-liberation" and "liberation" claims.


title: Gaps problems and frequencies of patches in cut and project sets

abstract: We establish a connection between gaps problems in Diophantine approximation and the frequency spectrum of patches in cut and project sets with special windows. Our theorems provide bounds for the number of distinct frequencies of patches of size r, which depend on the precise cut and project sets being used, and which are almost always less than a power of log r. Furthermore, for a substantial collection of cut and project sets we show that the number of frequencies of patches of size r remains bounded as r tends to infinity. The latter result applies to a collection of cut and project sets of full Hausdorff dimension.


title: Variational Formulation of E & M Particle Simulation Algorithms in Cylindrical Geometry using an Angular Modal Decomposition

abstract: Taking advantage of the flexibility of the variational method with coordinate transformations, we derive a self-consistent set of equations of motion from a discretized Lagrangian to study kinetic plasmas using a Fourier decomposed cylindrical coordinate system. The phase-space distribution function was reduced to a collection of finite-sized macro-particles of arbitrary shape moving on a virtual Cartesian grid. However, the discretization of field quantities was performed in cylindrical coordinates and decomposed into a truncated Fourier series in angle. A straightforward finite element interpolation scheme is used to transform between the two grids. The equations of motion were then obtained by demanding the action be stationary. The primary advantage of the variational approach is preservation of Lagrangian symmetries. In the present case, this leads to exact energy conservation, thus avoiding possible difficulties with grid heating.


title: Extraction of the distribution function $h^{\perp}_{1T}$ from experimental data

abstract: We attempt an extraction of the pretzelosity distribution ($h^{\perp}_{1T}$) from preliminary COMPASS, HERMES, and JLAB experimental data on $\sin(3\phi_h - \phi_S)$ asymmetry on proton, and effective deuteron and neutron targets. The resulting distributions, albeit with big errors, for the first time show tendency for up-quark pretzelosity to be positive and down-quark pretzelosity to be negative. A model relation of pretzelosity distribution and orbital angular momentum of quarks is used to estimate contributions of up and down quarks.


title: Study of the hidden-order of URu2Si2 by point contact tunnel junctions

abstract: URu2Si2 presents superconductivity at temperatures below 1.5 K, and a hidden order (HO) at about 17.5 K, both electronic phenomena are in uenced by Fano and Kondo resonances. At 17.5 K the HO was related in the past to a Peierls distortion that produces an energy gap deformed by the resonances. This order has been studied by more than 20 years, still without a clear understanding. In this work we studied the electronic characteristics of URu2Si2 in a single crystal, with tunneling and metallic point contact spectroscopies. In the superconducting state, we determined the energy gap, which shows the in uence of the Fano and Kondo resonances. At temperatures where HO is observed, the tunnel junctions spectra show the in uence of the two resonances. Tunnel junctions characteristics show that the Fermi surface nesting depends on the crystallographic direction.


title: Affective Facial Expression Processing via Simulation: A Probabilistic Model

abstract: Understanding the mental state of other people is an important skill for intelligent agents and robots to operate within social environments. However, the mental processes involved in `mind-reading' are complex. One explanation of such processes is Simulation Theory - it is supported by a large body of neuropsychological research. Yet, determining the best computational model or theory to use in simulation-style emotion detection, is far from being understood. In this work, we use Simulation Theory and neuroscience findings on Mirror-Neuron Systems as the basis for a novel computational model, as a way to handle affective facial expressions. The model is based on a probabilistic mapping of observations from multiple identities onto a single fixed identity (`internal transcoding of external stimuli'), and then onto a latent space (`phenomenological response'). Together with the proposed architecture we present some promising preliminary results


title: A Hitchhiker's Guide to Automatic Differentiation

abstract: This article provides an overview of some of the mathematical principles of Automatic Differentiation (AD). In particular, we summarise different descriptions of the Forward Mode of AD, like the matrix-vector product based approach, the idea of lifting functions to the algebra of dual numbers, the method of Taylor series expansion on dual numbers and the application of the push-forward operator, and explain why they all reduce to the same actual chain of computations. We further give a short mathematical description of some methods of higher-order Forward AD and, at the end of this paper, briefly describe the Reverse Mode of Automatic Differentiation.


title: On the toplogical computation of K4 of the Gaussian and Eisenstein integers

abstract: In this paper we use topological tools to investigate the structure of the algebraic K-groups K_4 (Z[i]) and K_4 (Z[rho]), where i := sqrt{-1} and rho := (1+sqrt{-3})/2. We exploit the close connection between homology groups of GL_n(R) for n <= 5 and those of related classifying spaces, then compute the former using Voronoi's reduction theory of positive definite quadratic and Hermitian forms to produce a very large finite cell complex on which GL_n(R) acts. Our main result is that K_4 (Z[i]) and K_4 (Z[rho]) have no p-torsion for p >= 5.


title: Exotic hadron states

abstract: Many charmonium-like and bottomonium-like $XYZ$ resonances have been observed by the Belle, Babar, CLEO and BESIII collaborations in the past decade. They are difficult to fit in the conventional quark model and thus are considered as candidates of exotic hadrons, such as multi-quark states, meson molecules, and hybrids. In this talk, we first briefly introduce the method of QCD sum rules and then provide a short review of the mass spectra of the quarkonium-like tetraquark states and the heavy quarkonium hybrids in the QCD sum rules approach. Possible interpretations of the $XYZ$ resonances are briefly discussed.


title: Random walks in a one-dimensional L\'evy random environment

abstract: We consider a generalization of a one-dimensional stochastic process known in the physical literature as L\'evy-Lorentz gas. The process describes the motion of a particle on the real line in the presence of a random array of marked points, whose nearest-neighbor distances are i.i.d. and long-tailed (with finite mean but possibly infinite variance). The motion is a continuous-time, constant-speed interpolation of a symmetric random walk on the marked points. We first study the quenched random walk on the point process, proving the CLT and the convergence of all the accordingly rescaled moments. Then we derive the quenched and annealed CLTs for the continuous-time process.

categories: math.PR cond-mat.dis-nn cond-mat.stat-mech math-ph math.MP


title: Quantum Uncertainty and Error-Disturbance Tradeoff

abstract: The uncertainty principle is often interpreted by the tradeoff between the error of a measurement and the consequential disturbance to the followed ones, which originated long ago from Heisenberg himself but now falls into reexamination and even heated debate. Here we show that the tradeoff is switched on or off by the quantum uncertainties of two involved non-commuting observables: if one is more certain than the other, there is no tradeoff; otherwise, they do have tradeoff and the Jensen-Shannon divergence gives it a good characterization.


title: On Detecting Noun-Adjective Agreement Errors in Bulgarian Language Using GATE

abstract: In this article, we describe an approach for automatic detection of noun-adjective agreement errors in Bulgarian texts by explaining the necessary steps required to develop a simple Java-based language processing application. For this purpose, we use the GATE language processing framework, which is capable of analyzing texts in Bulgarian language and can be embedded in software applications, accessed through a set of Java APIs. In our example application we also demonstrate how to use the functionality of GATE to perform regular expressions over annotations for detecting agreement errors in simple noun phrases formed by two words - attributive adjective and a noun, where the attributive adjective precedes the noun. The provided code samples can also be used as a starting point for implementing natural language processing functionalities in software applications related to language processing tasks like detection, annotation and retrieval of word groups meeting a specific set of criteria.


title: Modular proximal optimization for multidimensional total-variation regularization

abstract: We study \emph{TV regularization}, a widely used technique for eliciting structured sparsity. In particular, we propose efficient algorithms for computing prox-operators for $\ell_p$-norm TV. The most important among these is $\ell_1$-norm TV, for whose prox-operator we present a new geometric analysis which unveils a hitherto unknown connection to taut-string methods. This connection turns out to be remarkably useful as it shows how our geometry guided implementation results in efficient weighted and unweighted 1D-TV solvers, surpassing state-of-the-art methods. Our 1D-TV solvers provide the backbone for building more complex (two or higher-dimensional) TV solvers within a modular proximal optimization approach. We review the literature for an array of methods exploiting this strategy, and illustrate the benefits of our modular design through extensive suite of experiments on (i) image denoising, (ii) image deconvolution, (iii) four variants of fused-lasso, and (iv) video denoising. To underscore our claims and permit easy reproducibility, we provide all the reviewed and our new TV solvers in an easy to use multi-threaded C++, Matlab and Python library.


title: Sparse matrices describing iterations of integer-valued functions

abstract: We consider iterations of integer-valued functions $\phi$, which have no fixed points in the domain of positive integers. We define a local function $\phi_n$, which is a sub-function of $\phi$ being restricted to the subdomain $\{0, ..., n \}$. The iterations of $\phi_n$ can be described by a certain $n \times n$ sparse matrix $M_n$ and its powers. The determinant of the related $n \times n$ matrix $\hat{M}_n = I - M_n$, where $I$ is the identity matrix, acts as an indicator, whether the iterations of the local function $\phi_n$ enter a cycle or not. If $\phi_n$ has no cycle, then $\det \hat{M}_n = 1$ and the structure of the inverse $\hat{M}_n^{-1}$ can be characterized. Subsequently, we give applications to compute the inverse $\hat{M}_n^{-1}$ for some special functions. At the end, we discuss the results in connection with the $3x+1$ and related problems.


title: Bayesian feature selection with strongly-regularizing priors maps to the Ising Model

abstract: Identifying small subsets of features that are relevant for prediction and/or classification tasks is a central problem in machine learning and statistics. The feature selection task is especially important, and computationally difficult, for modern datasets where the number of features can be comparable to, or even exceed, the number of samples. Here, we show that feature selection with Bayesian inference takes a universal form and reduces to calculating the magnetizations of an Ising model, under some mild conditions. Our results exploit the observation that the evidence takes a universal form for strongly-regularizing priors --- priors that have a large effect on the posterior probability even in the infinite data limit. We derive explicit expressions for feature selection for generalized linear models, a large class of statistical techniques that include linear and logistic regression. We illustrate the power of our approach by analyzing feature selection in a logistic regression-based classifier trained to distinguish between the letters B and D in the notMNIST dataset.


title: Two-band description of resonant superfluidity in atomic Fermi gases

abstract: Fermionic superfluidity in atomic Fermi gases across a Feshbach resonance is normally described by the atom-molecule theory, which treats the closed channel as a noninteracting point boson. In this work we present a theoretical description of the resonant superfluidity in analogy to the two-band superconductors. We employ the underlying two-channel scattering model of Feshbach resonance where the closed channel is treated as a composite boson with binding energy $\varepsilon_0$ and the resonance is triggered by the microscopic interchannel coupling $U_{12}$. The binding energy $\varepsilon_0$ naturally serves as an energy scale of the system, which has been sent to infinity in the atom-molecule theory. We show that the atom-molecule theory can be viewed as a leading-order low-energy effective theory of the underlying fermionic theory in the limit $\varepsilon_0\rightarrow\infty$ and $U_{12}\rightarrow0$, while keeping the phenomenological atom-molecule coupling finite. The resulting two-band description of the superfluid state is in analogy to the BCS theory of two-band superconductors. In the dilute limit $\varepsilon_0\rightarrow\infty$, the two-band description recovers precisely the atom-molecule theory. The two-band theory provides a natural approach to study the corrections because of a finite binding energy $\varepsilon_0$ in realistic experimental systems. For broad and moderate resonances, the correction is not important for current experimental densities. However, for extremely narrow resonance, we find that the correction becomes significant. The finite binding energy correction could be important for the stability of homogeneous polarized superfluid against phase separation in imbalanced Fermi gases across a narrow Feshbach resonance.


title: Two-Variable Ehrenfeucht-Fraisse Games over Omega-Terms

abstract: Fragments of first-order logic over words can often be characterized in terms of finite monoids, and identities of omega-terms are an effective mechanism for specifying classes of monoids. Huschenbett and the first author have shown how to use infinite Ehrenfeucht-Fraisse games on linear orders for showing that some given fragment satisfies an identity of omega-terms (STACS 2014). After revisiting this result, we show that for two-variable logic one can use simpler linear orders.


title: Multi-Cell Processing with Limited Cooperation: A Novel Framework to Timely Designs and Reduced CSI Feedback with General Inputs

abstract: We investigate the optimal power allocation and optimal precoding for a multi-cell-processing (MCP) framework with limited cooperation. In particular, we consider two base stations(BSs) which maximize the achievable rate for two users connecting to each BS and sharing channel state information (CSI). We propose a two way channel estimation or prediction process. Such framework has promising outcomes in terms of feedback reduction and acheivable rates moving the system from one with unkown CSI at the transmitter to a system with instantanous CSI at both sides of the communication. We derive new extentions of the fundamental relation between the gradient of the mutual information and the MMSE for the conditional and non-conditional mutual information. Capitalizing on such relations, we provide the optimal power allocation and optimal precoding designs with respect to the estimated channel and MMSE. The designs introduced are optimal for multiple access (MAC) Gaussian coherent time-varying fading channels with general inputs and can be specialized to multiple input multiple output (MIMO) channels by decoding interference. The impact of interference on the capacity is quantified by the gradient of the mutual information with respect to the power, channel, and error covariance of the interferer. We provide two novel distributed MCP algorithms that provide the solutions for the optimal power allocation and optimal precoding for the UL and DL with a two way channel estimation to keep track of the channel variations over blocks of data transmission. Therefore, we provide a novel solution that allows with limited cooperation: a significant reduction in the CSI feedback from the receiver to the transmitter, and timely optimal designs of the precoding and power allocation.


title: Reverse enGENEering of regulatory networks from Big Data: a guide for a biologist

abstract: Omics technologies enable unbiased investigation of biological systems through massively parallel sequence acquisition or molecular measurements, bringing the life sciences into the era of Big Data. A central challenge posed by such omics datasets is how to transform this data into biological knowledge. For example, how to use this data to answer questions such as: which functional pathways are involved in cell differentiation? Which genes should we target to stop cancer? Network analysis is a powerful and general approach to solve this problem consisting of two fundamental stages, network reconstruction and network interrogation. Herein, we provide an overview of network analysis including a step by step guide on how to perform and use this approach to investigate a biological question. In this guide, we also include the software packages that we and others employ for each of the steps of a network analysis workflow.


title: SNIa Host Galaxy Properties and the Dust Extinction Distribution

abstract: Supernovae Type Ia display a complex relation with their host galaxies. An important prior to the fit of the supernovae's lightcurve is the distribution of host galaxy extinction values that can be encountered. The SDSS-SN project has published light curve fits using both MLCS2k2 and SALT2. We use the former fits extinction parameter ($A_V$) to map this distribution of extinction values. We explore the dependence of this distribution on four observables; the inclination of the host galaxy disk, radial position of the supernova, redshift of the supernova and host, and the level of star-formation in the host galaxy. The distribution of $A_V$ values encountered by supernovae is typically characterised by: $\rm N_0 ~ e^{-A_V/\tau}$, with $\tau$= 0.4 or 0.33. We find that the inclination correction using an infinitely thin disk for the SNIa is sufficient, resulting in similar exponential $A_V$ distributions for high- and low-inclination disks. The $A_V$ distribution also depends on the radial position in the disk, consistent with previous results on the transparency of spiral disks. The distribution of $A_V$ values narrows with increased star-formation, possibly due to the destruction or dispersion of the dusty ISM by stellar winds prior to the ignition of the supernova. In future supernova searches, certainly the inclination of the host galaxy disk, should be considered in the construction of the \av \ prior with $\tau=0.4/cos(i)$ as the most likely prior in each individual host galaxy's case.


title: CzeV404 - an eclipsing dwarf nova in the period gap during its July 2014 superoutburst

abstract: Results of the CCD observations of CzeV404 Her are displayed. During the season of June-August 2014 we detected one outburst and one superoutburst of the star. Clear superhumps with the period of P_sh=0.10472(2) days were observed. The superhump period was decreasing with a high value of P_dot=-2.43(8) x 10^(-4). For 17 eclipses, we calculated an orbital period with the value of P_orb=0.0980203(6) days which indicates that CzeV404 Her belongs to period gap objects and it is the longest orbital period eclipsing SU UMa star. Based on superhump and orbital period determinations, the period excess 6.8 % +/- 0.02 % and the mass ratio q ~ 0.32 of the system were obtained.


title: A New Method of Accelerated Bayesian Inference for Comparable Mass Binaries in both Ground and Space-Based Gravitational Wave Astronomy

abstract: With the advance in computational resources, Bayesian inference is increasingly becoming the standard tool of practise in GW astronomy. However, algorithms such as Markov Chain Monte Carlo (MCMC) require a large number of iterations to guarantee convergence to the target density. Each chain demands a large number of evaluations of the likelihood function, and in the case of a Hessian MCMC, calculations of the Fisher information matrix for use as a proposal distribution. As each iteration requires the generation of at least one gravitational waveform, we very quickly reach a point of exclusion for current Bayesian algorithms, especially for low mass systems where the length of the waveforms is large and the waveform generation time is on the order of seconds. This suddenly demands a timescale of many weeks for a single MCMC. As each likelihood and Fisher information matrix calculation requires the evaluation of noise-weighted scalar products, we demonstrate that by using the linearity of integration, and the fact that more than 90% of the generation time is spent at frequencies less that one third of the maximum, we can construct composite integrals that speed up the MCMCs for comparable mass binaries by a factor of between 3.5 and 5.5, depending on the waveform length. This method is both source and detector type independent, and can be applied to any waveform that displays significant frequency evolution, such as stellar mass binaries with Advanced LIGO/Virgo, as well as supermassive black holes with eLISA


title: Dynamic spatial regression models for space-varying forest stand tables

abstract: Many forest management planning decisions are based on information about the number of trees by species and diameter per unit area. This information is commonly summarized in a stand table, where a stand is defined as a group of forest trees of sufficiently uniform species composition, age, condition, or productivity to be considered a homogeneous unit for planning purposes. Typically information used to construct stand tables is gleaned from observed subsets of the forest selected using a probability-based sampling design. Such sampling campaigns are expensive and hence only a small number of sample units are typically observed. This data paucity means that stand tables can only be estimated for relatively large areal units. Contemporary forest management planning and spatially explicit ecosystem models require stand table input at higher spatial resolution than can be affordably provided using traditional approaches. We propose a dynamic multivariate Poisson spatial regression model that accommodates both spatial correlation between observed diameter distributions and also correlation between tree counts across diameter classes within each location. To improve fit and prediction at unobserved locations, diameter specific intensities can be estimated using auxiliary data such as management history or remotely sensed information. The proposed model is used to analyze a diverse forest inventory dataset collected on the United States Forest Service Penobscot Experimental Forest in Bradley, Maine. Results demonstrate that explicitly modeling the residual spatial structure via a multivariate Gaussian process and incorporating information about forest structure from LiDAR covariates improve model fit and can provide high spatial resolution stand table maps with associated estimates of uncertainty.


title: A boundary Schwarz Lemma for holomorphic mappings between unit balls of different dimensions

abstract: In this paper, we give a general boundary Schwarz lemma for holomorphic mappings between unit balls in any dimensions. It is proved that if the mapping $f\in C^{1+\alpha}$ at $z_0\in \partial \mathbb B^n$ with $f(z_0)=w_0\in \partial \mathbb B^N$ for any $n,N\geq 1$, then the Jacobian matrix $J_f(z_0)$ maps the tangent space $T_{z_0}(\partial \mathbb B^n)$ to $T_{w_0}(\partial \mathbb B^N)$, and the holomorphic tangent space $T^{(1,0)}_{z_0}(\partial \mathbb B^n)$ to $T^{(1,0)}_{w_0}(\partial \mathbb B^N)$ as well.


title: Experimental Demonstration of the Co-existence of the Spin Hall and Rashba Effects in beta-Tantalum/Ferromagnet Bilayers

abstract: We have measured the spin torques of beta-Tantalum / Co20Fe60B20 bilayers versus Ta thickness at room temperature using an FMR technique. The spin Hall coefficient was calculated both from the observed change in damping coefficient of the ferromagnet with Ta thickness, and from the ratio of the symmetric and anti-symmetric components of the FMR signal. Results from these two methods yielded values for the spin Hall coefficient of -0.090+/-0.005 and -0.11+/-0.01, respectively. We have also identified a significant out-of-plane spin torque originating from Ta, which is constant with Ta thickness. We ascribe this to an interface spin orbit coupling, or Rashba effect, due to the strength and constancy of the torque with Ta thickness. From fitting measured data to a model including interface spin orbit coupling, we have determined the spin diffusion length for beta-Tantalum to be ~2.5 nm.


title: Factorbird - a Parameter Server Approach to Distributed Matrix Factorization

abstract: We present Factorbird, a prototype of a parameter server approach for factorizing large matrices with Stochastic Gradient Descent-based algorithms. We designed Factorbird to meet the following desiderata: (a) scalability to tall and wide matrices with dozens of billions of non-zeros, (b) extensibility to different kinds of models and loss functions as long as they can be optimized using Stochastic Gradient Descent (SGD), and (c) adaptability to both batch and streaming scenarios. Factorbird uses a parameter server in order to scale to models that exceed the memory of an individual machine, and employs lock-free Hogwild!-style learning with a special partitioning scheme to drastically reduce conflicting updates. We also discuss other aspects of the design of our system such as how to efficiently grid search for hyperparameters at scale. We present experiments of Factorbird on a matrix built from a subset of Twitter's interaction graph, consisting of more than 38 billion non-zeros and about 200 million rows and columns, which is to the best of our knowledge the largest matrix on which factorization results have been reported in the literature.


title: A boundary Schwarz lemma for holomorphic mappings from the polydisc to the unit ball

abstract: In this paper, we prove a general Schwarz lemma at the boundary for holomorphic mappings from the polydisc to the unit ball in any dimensions. For the special case of one complex variable, the obtained results give the classic boundary Schwarz lemma.


title: Correlation-driven topological Fermi surface transition in FeSe

abstract: The electronic structure and phase stability of paramagnetic FeSe is computed by using a combination of ab initio methods for calculating band structure and dynamical mean-field theory. Our results reveal a topological change (Lifshitz transition) of the Fermi surface upon a moderate expansion of the lattice. The Lifshitz transition is accompanied with a sharp increase of the local moments and results in an entire reconstruction of magnetic correlations from the in-plane magnetic wave vector $(\pi,\pi)$ to $(\pi,0)$. We attribute this behavior to a correlation-induced shift of the Van Hove singularity originating from the $d_{xy}$ and $d_{xz}/d_{yz}$ bands at the M-point across the Fermi level. We propose that superconductivity is strongly influenced, or even induced, by a Van Hove singularity.


title: The metaphysics of laws: dispositionalism vs. primitivism

abstract: The paper compares dispositionalism about laws of nature with primitivism. It argues that while the distinction between these two positions can be drawn in a clear-cut manner in classical mechanics, it is less clear in quantum mechanics, due to quantum non-locality. Nonetheless, the paper points out advantages for dispositionalism in comparison to primitivism also in the area of quantum mechanics, and of contemporary physics in general.


title: clustvarsel: A Package Implementing Variable Selection for Model-based Clustering in R

abstract: Finite mixture modelling provides a framework for cluster analysis based on parsimonious Gaussian mixture models. Variable or feature selection is of particular importance in situations where only a subset of the available variables provide clustering information. This enables the selection of a more parsimonious model, yielding more efficient estimates, a clearer interpretation and, often, improved clustering partitions. This paper describes the R package clustvarsel which performs subset selection for model-based clustering. An improved version of the methodology of Raftery and Dean (2006) is implemented in the new version 2 of the package to find the (locally) optimal subset of variables with group/cluster information in a dataset. Search over the solution space is performed using either a stepwise greedy search or a headlong algorithm. Adjustments for speeding up these algorithms are discussed, as well as a parallel implementation of the stepwise search. Usage of the package is presented through the discussion of several data examples.


title: Effective Polyakov line actions, and their solutions at finite chemical potential

abstract: I outline recent progress in the relative weights approach to deriving effective Polyakov line actions from an underlying lattice gauge theory, and compare mean field and complex Langevin methods for solving such theories at finite chemical potential.


title: Critical coupling to Tamm plasmons

abstract: The conditions of critical coupling of light to Tamm plasmons are investigated with comprehensive numerical simulations, highlighting the parameters that maximise absorption of incident light in the metal layer. The asymmetric response in reflection and absorption with respect to the direction of incidence is discussed, the two cases yielding different optimal coupling conditions. These findings are relevant for the design of optimised Tamm structures, particularly in applications such as narrow-band thermal emitters, field-enhanced spectroscopy and refractive-index sensing.


title: On a mass functional for initial data in 4+1 dimensional spacetime

abstract: We consider a broad class of asymptotically flat, maximal initial data sets satisfying the vacuum constraint equations, admitting two commuting rotational symmetries. We construct a mass functional for `$t-\phi^i$' symmetric data which evaluates to the ADM mass. We then show that $\mathbb{R} \times U(1)^2$-invariant solutions of the vacuum Einstein equations are critical points of this functional amongst this class of data. We demonstrate positivity of this functional for a class of rod structures which include the Myers-Perry initial data. The construction is a natural extension of Dain's mass functional to $D=5$, although several new features arise.


title: Planting colourings silently

abstract: Let $k\geq3$ be a fixed integer and let $Z_k(G)$ be the number of $k$-colourings of the graph $G$. For certain values of the average degree, the random variable $Z_k(G(n,m))$ is known to be concentrated in the sense that $\frac1n(\ln Z_k(G(n,m))-\ln E[Z_k(G(n,m))])$ converges to $0$ in probability [Achlioptas and Coja-Oghlan: FOCS 2008]. In the present paper we prove a significantly stronger concentration result. Namely, we show that for a wide range of average degrees, $\frac1\omega(\ln Z_k(G(n,m))-\ln E[Z_k(G(n,m))])$ converges to $0$ in probability for any diverging function $\omega=\omega(n)\to\infty$. For $k$ exceeding a certain constant $k_0$ this result covers all average degrees up to the so-called condensation phase transition, and this is best possible. As an application, we show that the experiment of choosing a $k$-colouring of the random graph $G(n,m)$ uniformly at random is contiguous with respect to the so-called "planted model".


title: Reaction rates for mesoscopic reaction-diffusion kinetics

abstract: The mesoscopic reaction-diffusion master equation (RDME) is a popular modeling framework, frequently applied to stochastic reaction-diffusion kinetics in systems biology. The RDME is derived from assumptions about the underlying physical properties of the system, and it may produce unphysical results for models where those assumptions fail. In that case, other more comprehensive models are better suited, such as hard-sphere Brownian dynamics (BD). Although the RDME is a model in its own right, and not inferred from any specific microscale model, it proves useful to attempt to approximate a microscale model by a specific choice of mesoscopic reaction rates. In this paper we derive mesoscopic reaction rates by matching certain statistics of the RDME solution to statistics of the solution of a widely used microscopic BD model: the Smoluchowski model with a mixed boundary condition at the reaction radius of two molecules. We also establish fundamental limits for the range of mesh resolutions for which this approach yields accurate results, and show both theoretically and in numerical examples that as we approach the lower fundamental limit, the mesoscopic dynamics approach the microscopic dynamics.


title: Radiation-induced re-emission in a 2D electron system

abstract: Recent experiments on re-emission of radiation and magnetotransport in photoexcited two-dimensional electron systems are theoretically analyzed. These experiments have concluded that there exists a strong correlation between the re-emitted radiation and the radiation-driven current through the sample. We study these remarkable experimental results using the radiationdriven electron orbit model. According to it, two-dimensional electrons under a magnetic field move back and forth driven by radiation. Therefore, they behave as oscillating dipoles re-emitting part of the radiation previously absorbed. Then, this theory would explain in a simple way the experimental results.


title: Topological Data Analysis and Cosheaves

abstract: This paper contains an expository account of persistent homology and its usefulness for topological data analysis. An alternative foundation for level-set persistence is presented using sheaves and cosheaves.


title: Time-Periodic Solutions of Driven-Damped Trimer Granular Crystals

abstract: In this work, we consider time-periodic structures of trimer granular crystals consisting of alternate chrome steel and tungsten carbide spherical particles yielding a spatial periodicity of three. The configuration at the left boundary is driven by a harmonic in-time actuation with given amplitude and frequency while the right one is a fixed wall. Similar to the case of a dimer chain, the combination of dissipation, driving of the boundary, and intrinsic nonlinearity leads to complex dynamics. For fixed driving frequencies in each of the spectral gaps, we find that the nonlinear surface modes and the states dictated by the linear drive collide in a saddle-node bifurcation as the driving amplitude is increased, beyond which the dynamics of the system become chaotic. While the bifurcation structure is similar for solutions within the first and second gap, those in the first gap appear to be less robust. We also conduct a continuation in driving frequency, where it is apparent that the nonlinearity of the system results in a complex bifurcation diagram, involving an intricate set of loops of branches, especially within the spectral gap. The theoretical findings are qualitatively corroborated by the experimental full-field visualization of the time-periodic structures.


title: Cheeger-Mueller Theorem on manifolds with cusps

abstract: We prove equality between the renormalized Ray-Singer analytic torsion and the intersection R-torsion on a Witt-manifold with cusps, up to an error term determined explicitly by the Betti numbers of the cross section of the cusp and the intersection R-torsion of a model cone. In the first step of the proof we compute explicitly the renormalized Ray-Singer analytic torsion of a model cusp in general dimension and without the Witt-condition. In the second step we establish a gluing formula for renormalized Ray-Singer analytic torsion on a general class of non-compact manifolds in any dimension that includes Witt-manifolds with cusps, but also scattering manifolds with asymptotically conical ends. In the final step, a Cheeger-Mueller theorem on cusps follows by a combination of the previous explicit computation and the gluing formula.


title: On the wellposedness of the exp-Rabelo equation

abstract: The exp-Rabelo equation describes pseudo-spherical surfaces. It is a nonlinear evolution equation. In this paper the wellposedness of bounded from above solutions for the initial value problem associated to this equation is studied.


title: Wellposedness of the Ostrovsky-Hunter equation Under the combined effects of dissipation and short wave dispersion

abstract: The Ostrovsky-Hunter equation provides a model for small-amplitude long waves in a rotating fluid of finite depth. It is a nonlinear evolution equation. In this paper we study the well-posedness for the Cauchy problem associated to this equation in presence of some weak dissipation effects.


title: Entropic destruction of heavy quarkonium in non-Abelian plasma from holography

abstract: Lattice QCD indicates a large amount of entropy associated with the heavy quark-antiquark pair immersed in the quark-gluon plasma. This entropy grows as a function of the inter-quark distance giving rise to an entropic force that can be very effective in dissociating the bound quarkonium states. In addition, the lattice data show a very sharp peak in the heavy quark-antiquark entropy at the deconfinement transition. Since the quark-gluon plasma around the deconfinement transition is strongly coupled, we employ the holographic correspondence to study the entropy associated with the heavy quark-antiquark pair in two theories: i) ${\cal{N}}=4$ supersymmetric Yang-Mills and ii) a confining Yang-Mills theory obtained by compactification on a Kaluza-Klein circle. In both cases we find the entropy growing with the inter-quark distance and evaluate the effect of the corresponding entropic forces. In the case ii), we find a sharp peak in the entropy near the deconfinement transition, in agreement with the lattice QCD results. This peak in our holographic description arises because the heavy quark pair acts as an eyewitness of the black hole formation in the bulk -- the process that describes the deconfinement transition. In terms of the boundary theory, this entropy likely emerges from the entanglement of a "long string" connecting the quark and antiquark with the rest of the system.


title: Gamma-Ray and Hard X-Ray Emission from Pulsar-Aided Supernovae as a Probe of Particle Acceleration in Embryonic Pulsar Wind Nebulae

abstract: It has been suggested that some classes of luminous supernovae (SNe) and gamma-ray bursts (GRBs) are driven by newborn magnetars. Fast-rotating proto-neutron stars have also been of interest as potential sources of gravitational waves (GWs). We show that for a range of rotation periods and magnetic fields, hard X rays and GeV gamma rays provide us with a promising probe of pulsar-aided SNe. It is observationally known that young pulsar wind nebulae (PWNe) in the Milky Way are very efficient lepton accelerators. We argue that, if embryonic PWNe satisfy similar conditions at early stages of SNe (in ~1-10 months after the explosion), external inverse-Compton emission via upscatterings of SN photons is naturally expected in the GeV range as well as broadband synchrotron emission. To fully take into account the Klein-Nishina effect and two-photon annihilation process that are important at early times, we perform detailed calculations including electromagnetic cascades. Our results suggest that hard X-ray telescopes such as NuSTAR can observe such early PWN emission by followup observations in months-to-years. GeV gamma rays may also be detected by Fermi for nearby SNe, which serve as counterparts of these GW sources. Detecting the signals will give us an interesting probe of particle acceleration at early times of PWNe, as well as clues to driving mechanisms of luminous SNe and GRBs. Since the Bethe-Heitler cross section is lower than the Thomson cross section, gamma rays would allow us to study subphotospheric dissipation. We encourage searches for high-energy emission from nearby SNe, especially Type Ibc SNe including super-luminous objects.


title: Classical, quantum, and phenomenological aspects of dark energy models

abstract: The origin of accelerating expansion of the Universe is one the biggest conundrum of fundamental physics. In this paper we review vacuum energy issues as the origin of accelerating expansion - generally called dark energy - and give an overview of alternatives, which a large number of them can be classified as interacting scalar field models. We review properties of these models both as classical field and as quantum condensates in the framework of non-equilibrium quantum field theory. Finally, we review phenomenology of models with the goal of discriminating between them.


title: Dealing in practice with hot-spots

abstract: The hot-spot phenomenon is a relatively frequent problem occurring in current photovoltaic generators. It entails both a risk for the photovoltaic module's lifetime and a decrease in its operational efficiency. Nevertheless, there is still a lack of widely accepted procedures for dealing with them in practice. This paper presents the IES-UPM observations on 200 affected modules. Visual and infrared inspection, electroluminescence, peak power and operating voltage tests have been accomplished. Hot-spot observation procedures and well defined acceptance and rejection criteria are proposed, addressing both the lifetime and the operational efficiency of the modules. The operating voltage has come out as the best parameter to control effective efficiency losses for the affected modules. This procedure is oriented to its possible application in contractual frameworks.


title: A Subspace Method for Array Covariance Matrix Estimation

abstract: This paper introduces a subspace method for the estimation of an array covariance matrix. It is shown that when the received signals are uncorrelated, the true array covariance matrices lie in a specific subspace whose dimension is typically much smaller than the dimension of the full space. Based on this idea, a subspace based covariance matrix estimator is proposed. The estimator is obtained as a solution to a semi-definite convex optimization problem. While the optimization problem has no closed-form solution, a nearly optimal closed-form solution is proposed making it easy to implement. In comparison to the conventional approaches, the proposed method yields higher estimation accuracy because it eliminates the estimation error which does not lie in the subspace of the true covariance matrices. The numerical examples indicate that the proposed covariance matrix estimator can significantly improve the estimation quality of the covariance matrix.


title: Search for CP violation in $D^0 \to h^+h^-$ decays at CDF

abstract: I report on a measurement of \CP-violating asymmetries ($A_{\Gamma}$) between effective lifetimes of $D^0$ or $\bar{D}^0$ in fully reconstructed $D^0\to K^+ K^-$ and $D^0\to \pi^+\pi^-$ decays collected in $p\bar{p}$ collisions by the Collider Detector at Fermilab experiment. The full CDF data set corresponding to $9.7$ fb$^{-1}$ of integrated luminosity is used. The flavor of the charm meson at production is determined by exploiting the strong-interaction decay $D^{*+} \to D^0 \pi^+$, while the contamination from mesons originated in $b$-hadron decays is evaluated and subtracted from the sample. Signal yields as functions of the observed decay-time distributions are extracted from maximum likelihood fits and used to measure the asymmetries. The results, $A_\Gamma (K^+K^-) = \bigl(-1.9 \pm 1.5 \mathrm{(stat)} \pm 0.4 \mathrm{(syst)} \bigr)\times10^{-3}$ and $A_\Gamma (\pi^+\pi^-)= \bigl(-0.1 \pm 1.8 \mathrm{(stat)} \pm 0.3 \mathrm{(syst)} \bigr)\times10^{-3}$, and their combination, $A_\Gamma = \bigl(-1.2 \pm 1.2\bigr)\times10^{-3}$, are consistent with the SM predictions and other experimental determinations.


title: On the Stanley depth of edge ideals of line and cyclic graphs

abstract: We prove that the edge ideals of line and cyclic graphs and their quotient rings satisfy the Stanley conjecture. We compute the Stanley depth for the quotient ring of the edge ideal associated to a cycle graph of length $n$, given a precise formula for $n\equiv 0,2 \pmod{3}$ and tight bounds for $n\equiv 1 \pmod{3}$. Also, we give bounds for the Stanley depth of a quotient of two monomial ideals, in combinatorial terms.


title: Schatten class Toeplitz operators acting on large weighted Bergman spaces

abstract: A full description of the membership in the Schatten ideal $S_ p(A^2_{\omega})$ for $0<p<\infty$ of the Toeplitz operator acting on large weighted Bergman spaces is obtained.


title: An Extremal Chiral Primary Three-Point Function at Two-loops in ABJ(M)

abstract: I compute the leading correction to the structure constant for the three-point function of two length-two and one length-four chiral primary operators in planar ABJ(M) theory at weak 't Hooft coupling. The computation is reduced to four-loop propagator type Feynman integrals via a manifestly finite integration over the position of the length-four operator.


title: On the structure of instability in moduli theory

abstract: We formulate a theory of instability and Harder-Narasimhan filtrations for an arbitrary moduli problem in algebraic geometry. We introduce the notion of a $\Theta$-stratification of a moduli problem, which generalizes the Kempf-Ness stratification in GIT as well as the Harder-Narasimhan stratification of the moduli of coherent sheaves on a projective scheme. Our main theorems establish necessary and sufficient conditions for the existence of these stratifications. We define a structure on an algebraic stack called a numerical invariant, and we show that in many situations a numerical invariant defines a $\Theta$-stratification on the stack, assuming a certain "HN boundedness" condition holds. We also discuss criteria under which the semistable locus has a moduli space. We apply our methods to an example that lies beyond the reach of classical methods: the stratification of the stack of objects in the heart of a Bridgeland stability condition.


title: PH = PSPACE

abstract: In this paper we show that PSPACE is equal to 4th level in the polynomial hierarchy. We also deduce a lot of important consequences. True quantified Boolean formula is a generalisation of the Boolean Satisfiability Problem, where determining of interpretation that satisfies a given Boolean formula is replaced by existence of Boolean functions that makes a given QBF to be tautology. Such functions are called the Skolem functions. The essential idea of the proof is to show that for any quantified Boolean formula $\phi$ we can obtain a formula $\phi'$ which is in the 4th level of the polynomial hierarchy, no more than polynomial in the size of a given $\phi$, such that the truth of $\phi$ can be determined from the truth of $\phi'$. The idea is to skolemize, and then use additional formulas from the 2nd level of the polynomial hierarchy inside the skolemized prefix to enforce that the skolem variables indeed depend only on the universally quantified variables they are supposed to. However, some dependence is lost when the quantification is reversed. It is called "XOR issue" because the functional dependence can be expressed by means of an XOR formula. Thus, it is needed to locate these XORs, but there is no need to locate all chains with XORs: any chain includes a XOR of only two variables. The last can be done locally in each iteration (keep in mind the algebraic normal form (ANF)), when all arguments are specified, i.e. as a polynomial subroutine. Relativization is defeated due to the well-known fact: PH = PSPACE iff second-order logic over finite structures gains no additional power from the addition of a transitive closure operator. The exchange is possible due to finite possibilities for arguments. So, the theorems with oracles are not applicable since a random oracle is an arbitrary set. And that's why PH is infinite relative to a random oracle with probability 1.


title: The Case for Muon-based Neutrino Beams

abstract: For the foreseeable future, high energy physics accelerator capabilities in the US will be deployed to study the physics of the neutrino sector. In this context, it is useful to explore the sensitivities and limiting systematic effects of the planned neutrino oscillation program, so that we can evaluate the issues that must be addressed in order to ensure the success of these efforts. It is only in this way that we will ultimately be able to elucidate the fundamental physics processes involved. We conclude that success can only be guaranteed by, at some point in the future, being able to deploy muon accelerator capabilities. Such capabilities provide the only route to precision neutrino beams with which to study and mitigate, at the sub-percent level, the limiting systematic issues of future oscillation measurements. Thus this analysis argues strongly for maintaining a viable accelerator research program towards future muon accelerator capabilities.


title: Active Inference for Binary Symmetric Hidden Markov Models

abstract: We consider active maximum a posteriori (MAP) inference problem for Hidden Markov Models (HMM), where, given an initial MAP estimate of the hidden sequence, we select to label certain states in the sequence to improve the estimation accuracy of the remaining states. We develop an analytical approach to this problem for the case of binary symmetric HMMs, and obtain a closed form solution that relates the expected error reduction to model parameters under the specified active inference scheme. We then use this solution to determine most optimal active inference scheme in terms of error reduction, and examine the relation of those schemes to heuristic principles of uncertainty reduction and solution unicity.

categories: stat.ML cond-mat.dis-nn cs.IT cs.LG math.IT


title: On the behavior of the bounds of the holographic theory for massive and massless particle systems

abstract: The aim of the present dissertation is to analyze the meaning of the entropy bounds of the holographic sector once tested for statistical ensembles of particles, in order to deeper investigate the nature of these constraints and their mutual links. From the Universal Energy Bound (UEB) simple time constraints can be argued, which are manifestations of the discrete nature of the space-time and of the presence of ultimate space-time scales. From the combined effort of the UEB and of the Holographic bound by 't Hooft and Susskind, an entropy density bound as a function of temperature is achieved.


title: Disappearing charged tracks in association with displaced leptons from supersymmetry

abstract: In this note we discuss a characteristic collider signature of models of gauge mediated supersymmetry breaking in which the selectron and smuon are mass-degenerate co-NLSP. In these models, all processes involving superpartners give rise to two NLSP selectrons or smuons, each of which subsequently decays to a nearly massless LSP gravitino and an electron or a muon. In a large region of the parameter space, the NLSPs travel macroscopic distances, of the order 0.1-1000 mm, before they decay. Hence, these models give rise to collider signatures involving charged tracks that end at vertices, which are inside the detector volume but displaced with respect to the original collision point, from which an electron or a muon is emitted. In order to probe this class of models we propose a search for disappearing charged tracks in association with displaced electrons or muons.


title: Measure of compactness for filters in product spaces: Kuratowski-Mr\`owka in CAP revisited

abstract: The first author introduced a measure of compactness for families of sets, relative to a class of filters, in the context of convergence approach spaces. We characterize a variety of maps (types of quotient maps, closed maps, and variants of perfect maps) as those respecting this measure of compactness under one form or another. We establish a product theorem for measure of compactness that yields as instances new product theorems for spaces and maps, and new product characterizations of spaces and maps, thus extending existing results from the category of convergence spaces to that of convergence approach spaces. In particular, results of the Mr\`owka-Kuratowski type are obtained, shedding new light on existing results for approach spaces.


title: A note on stability and optimal approximation estimates for symmetric saddle point systems

abstract: We establish sharp well-posedness and approximation estimates for variational saddle point systems at the continuous level. The main results of this note have been known to be true only in the finite dimensional case. Known spectral results from the discrete case are reformulated and proved using a functional analysis view, making the proofs in both cases, discrete and continuous, less technical than the known discrete approaches. We focus on analyzing the special case when the form $a(\cdot, \cdot)$ is bounded, symmetric, and coercive, and the mixed form $b(\cdot, \cdot)$ is bounded and satisfies a standard $\inf-\sup$ or LBB condition. We characterize the spectrum of the symmetric operators that describe the problem at the continuous level. For a particular choice of the inner product on the product space of $b(\cdot, \cdot)$, we prove that the spectrum of the operator representing the system at continuous level is $\left \{\frac{1-\sqrt{5}}{2}, 1, \frac{1+\sqrt{5}}{2} \right \}$. As consequences of the spectral description, we find the minimal length interval that contains the ratio between the norm of the data and the norm of the solution, and prove explicit approximation estimates that depend only on the continuity constant and the continuous and the discrete $\inf-\sup$ condition constants.


title: A Symmetry Approach to Geometric Phase for Quantum Ensembles

abstract: We use tools from the theory of dynamical systems with symmetries to stratify Uhlmann's standard purification bundle and derive a new connection for mixed quantum states. For unitarily evolving systems, this connection gives rise to the 'interferometric' geometric phase of Sj\"{o}qvist et al. [Phys. Rev. Lett. 85 2845 - 2849 (2000)], and for more generally evolving open systems it gives rise to the generalization of the interferometric geometric phase due to Tong et al. [Phys. Rev. Lett. 93 080405 (2004)].


title: Stellar magnetic activity and Star-Planet Interactions (invited review)

abstract: Stellar magnetic activity is an important factor in the formation and evolution of exoplanets. Magnetic phenomena like stellar flares, coronal mass ejections, and high-energy emission affect the exoplanetary atmosphere and its mass loss over time. One major question is whether the magnetic evolution of exoplanet host stars is the same as for stars without planets; tidal and magnetic interactions of a star and its close-in planets may play a role in this. Stellar magnetic activity also shapes our ability to detect exoplanets with different methods in the first place, and therefore we need to understand it properly to derive an accurate estimate of the existing exoplanet population. Here I review recent theoretical and observational results, as well as outline some avenues for future progress.


title: Tracing early stellar evolution with asteroseismology: pre-main sequence stars in NGC 2264

abstract: Asteroseismology has been proven to be a successful tool to unravel details of the internal structure for different types of stars in various stages of their main sequence and post-main sequence evolution. Recently, we found a relation between the detected pulsation properties in a sample of 34 pre-main sequence (pre-MS) delta Scuti stars and the relative phase in their pre-MS evolution. With this we are able to demonstrate that asteroseismology is similarly powerful if applied to stars in the earliest stages of evolution before the onset of hydrogen core burning.


title: Obliquities of "Top-Shaped" Asteroids May Not Imply Reshaping by YORP Spin-up

abstract: The timescales over which the YORP effect alters the rotation period and the obliquity of a small asteroid can be very different, because the corresponding torques couple to different aspects of the object's shape. For nearly axisymmetric, "top-shaped" near-Earth asteroids such as 101955 Bennu, spin timescales are an order of magnitude or more longer than obliquity timescales, which are of order 10^5 to 10^6 yr. The observed low obliquities (near 0 or 180 degrees) of top-shaped asteroids do not constitute evidence that they acquired their present shapes and spins through YORP spin-up, because low obliquities are expected regardless of the spin-up or reshaping mechanism.


title: The Feller Property for Graphs

abstract: The Feller property concerns the preservation of the space of functions vanishing at infinity by the semigroup generated by an operator. We study this property in the case of the Laplacian on infinite graphs with arbitrary edge weights and vertex measures. In particular, we give conditions for the Feller property involving curvature-type quantities for general graphs, characterize the property in the case of model graphs, and give some comparison results to the model case.


title: Electric field of a point charge in truncated hyperbolic motion

abstract: We find the electric field of a point charge in `truncated hyperbolic motion', in which the charge moves at a constant velocity followed by motion with a constant acceleration in its instantaneous rest frame. The same Lienard-Wiechert formula holds for the acceleration phase and the constant velocity phase of the charge's motion. The only modification is that the formula giving the retarded time is different for the two motions, and the acceleration is zero for the constant velocity motion. The electric field lines are continuous as the retarded time increases through the transition time between constant velocity and accelerated motion. As the transition time approaches negative infinity the electric field develops a delta function contribution that has been introduced by others as necessary to preserve Gauss's law for the electric field.


title: Astrometric Calibration of the Gemini NICI Planet-Finding Campaign

abstract: We describe the astrometric calibration of the Gemini NICI Planet-Finding Campaign. The Campaign requires a relative astrometric accuracy of $\approx$ 20 mas across multi-year timescales in order to distinguish true companions from background stars by verifying common proper motion and parallax with their parent stars. The calibration consists of a correction for instrumental optical image distortion, plus on-sky imaging of astrometric fields to determine the pixel scale and image orientation. We achieve an accuracy of $\lesssim 7$ mas between the center and edge of the 18$''$ NICI field, meeting the 20 mas requirement. Most of the Campaign data in the Gemini Science Archive are accurate to this level but we identify a number of anomalies and present methods to correct the errors.


title: An Alternative to Maxwell's Displacement Current

abstract: Though sufficient for local conservation of charge, we show that Maxwells displacement current is not necessary. An alternative to the Ampere Maxwell equation is exhibited and the alternative s electric and magnetic fields and scalar and vector potentials are expressed in terms of the charge and current densities. The alternative describes a theory in which action is instantaneous and so may provide a good approximation to Maxwells equations where and when the finite speed of light can be neglected. The result is reminiscent of the Darwin approximation which arose from the study classical charged point particles to order (v/c)2 in the Lagrangian. Unlike Darwin, this approach does not depend on the constitution of the electric current. Instead, this approach grows from a straightforward revision of the Ampere Equation which revision enforces the local conservation of charge.


title: Optical study of the hyper-luminous X-ray source 2XMM J011942.7+032421

abstract: We present the identification and characterization of the optical counterpart to 2XMM J011942.7+032421, one of the most luminous and distant ultra-luminous X-ray sources (ULXs). The counterpart is located near a star forming region in a spiral arm of the galaxy NGC 470 with u, g, and r magnitudes of 21.53, 21.69, and 21.71 mags, respectively. The luminosity of the counterpart is much larger than that of a single O-type star, indicating that it may be a stellar cluster. Our optical spectroscopic observations confirm the association of the X-ray source and the optical counterpart with its host galaxy NGC 470, which validates the high, > 10^41 erg/s, X-ray luminosity of the source. Its optical spectrum is embedded with numerous emission lines, including H recombination lines, metallic forbidden lines and more notably the high-ionization HeII (lambda 4686 A) line. This line shows a large velocity dispersion of $\simeq$ 410 \kms, consistent with the existence of a compact (< 5 AU) highly-ionized accretion disc rotating around the central X-ray source. The 1.4 x 10^37 erg/s luminosity of the HeII line emission makes the source one of the most luminous ULXs in the emission of that line. This, together with the high X-ray luminosity and the large velocity dispersion of the HeII emission, suggests that the source is an ideal candidate for more extensive follow-up observations for understanding the nature of hyper-luminous X-ray sources, a more luminous subgroup of ULXs and more likely candidates for intermediate-mass black holes.


title: Towards Tight Lower Bounds for Range Reporting on the RAM

abstract: In the orthogonal range reporting problem, we are to preprocess a set of $n$ points with integer coordinates on a $U \times U$ grid. The goal is to support reporting all $k$ points inside an axis-aligned query rectangle. This is one of the most fundamental data structure problems in databases and computational geometry. Despite the importance of the problem its complexity remains unresolved in the word-RAM. On the upper bound side, three best tradeoffs exists: (1.) Query time $O(\lg \lg n + k)$ with $O(nlg^{\varepsilon}n)$ words of space for any constant $\varepsilon>0$. (2.) Query time $O((1 + k) \lg \lg n)$ with $O(n \lg \lg n)$ words of space. (3.) Query time $O((1+k)\lg^{\varepsilon} n)$ with optimal $O(n)$ words of space. However, the only known query time lower bound is $\Omega(\log \log n +k)$, even for linear space data structures. All three current best upper bound tradeoffs are derived by reducing range reporting to a ball-inheritance problem. Ball-inheritance is a problem that essentially encapsulates all previous attempts at solving range reporting in the word-RAM. In this paper we make progress towards closing the gap between the upper and lower bounds for range reporting by proving cell probe lower bounds for ball-inheritance. Our lower bounds are tight for a large range of parameters, excluding any further progress for range reporting using the ball-inheritance reduction.


title: Reverse Hardy-type inequalities for supremal operators with measures

abstract: In this paper we characterize the validity of the inequalities $\|g\|_{p,(a,b),\lambda} \le c \|u(x) \|g\|_{\infty,(x,b),\mu}\|_{q,(a,b),\nu}$ and $\label{eq.0.1.2} \|g\|_{p,(a,b),\lambda} \le c \|u(x) \|g\|_{\infty,(a,x),\mu}\|_{q,(a,b),\nu} $ for all non-negative Borel measurable functions $g$ on the interval $(a,b)$, where $0 < p \le +\infty$, $0 < q \le +\infty$, $\lambda$, $\mu$ and $\nu$ are non-negative Borel measures on $(a,b)$, and $u$ is a weight function on $(a,b)$.


title: Curvature effects in statics and dynamics of low dimensional magnets

abstract: We develop an approach to treat magnetic energy of a ferromagnet for arbitrary curved wires and shells on the assumption that the anisotropy contribution much exceeds the dipolar and other weak interactions. We show that the curvature induces two effective magnetic interactions: effective magnetic anisotropy and effective Dzyaloshinskii-like interaction. We derive an equation of magnetisation dynamics and propose a general static solution for the limit case of strong anisotropy. To illustrate our approach we consider the magnetisation structure in a ring wire and a cone surface: ground states in both systems essentially depend on the curvature excluding strictly tangential solutions even in the case of strong anisotropy. We derive also the spectrum of spin waves in such systems.


title: Multiple Imputation Using Gaussian Copulas

abstract: Missing observations are pervasive throughout empirical research, especially in the social sciences. Despite multiple approaches to dealing adequately with missing data, many scholars still fail to address this vital issue. In this paper, we present a simple-to-use method for generating multiple imputations using a Gaussian copula. The Gaussian copula for multiple imputation (Hoff, 2007) allows scholars to attain estimation results that have good coverage and small bias. The use of copulas to model the dependence among variables will enable researchers to construct valid joint distributions of the data, even without knowledge of the actual underlying marginal distributions. Multiple imputations are then generated by drawing observations from the resulting posterior joint distribution and replacing the missing values. Using simulated and observational data from published social science research, we compare imputation via Gaussian copulas with two other widely used imputation methods: MICE and Amelia II. Our results suggest that the Gaussian copula approach has a slightly smaller bias, higher coverage rates, and narrower confidence intervals compared to the other methods. This is especially true when the variables with missing data are not normally distributed. These results, combined with theoretical guarantees and ease-of-use suggest that the approach examined provides an attractive alternative for applied researchers undertaking multiple imputations.


title: Random flights related to the Euler-Poisson-Darboux equation

abstract: This paper is devoted to the analysis of random motions on the line and in the space R^d (d > 1) performed at finite velocity and governed by a non-homogeneous Poisson process with rate \lambda(t). The explicit distributions p(x,t) of the position of the randomly moving particles are obtained solving initial-value problems for the Euler- Poisson-Darboux equation when \lambda(t) = \alpha/t, t > 0. We consider also the case where \lambda(t) = \lambda coth \lambda t and \lambda(t) = \lambda tanh \lambda t, where some Riccati differential equations emerge and the explicit distributions are obtained for d = 1. We also examine planar random motions with random velocities by projecting random flights in R^d onto the plane. Finally the case of planar motions with four orthogonal directions is considered and the corresponding higher-order equations with time-varying coefficients obtained.


title: A proto-cluster at z=2.45

abstract: We present the spectroscopic confirmation of a $z=2.45$ proto-cluster. Its member galaxies lie within a radius of 1.4Mpc (physical) on the sky and within $\Delta v \pm 700$km/s along the line of sight. We estimate an overdensity of 10, suggesting that the structure has made the turn-around but is not assembled yet. Comparison to the Millennium simulation suggests that analogous structures evolve into $10^{14}-10^{15}$M$_{\odot}$/h type dark matter haloes by $z=0$ qualifying the notion of "proto-cluster". The search for the complete census of mock progenitor galaxies at $z\sim2.5$ of these massive $z=0$ mock clusters reveals that they are widely spread over areas with a diameter of 3-20Mpc. This suggests that the optical selection of such proto-clusters can result in a rich diversity regarding their $z=0$ descendants. We also searched for signs of environmental differentiation in this proto-cluster. Whilst we see a weak trend for more massive and more quiescent galaxies in the proto-cluster, this is not statistically significant.


title: Proof of the satisfiability conjecture for large k

abstract: We establish the satisfiability threshold for random $k$-SAT for all $k\ge k_0$, with $k_0$ an absolute constant. That is, there exists a limiting density $\alpha_*(k)$ such that a random $k$-SAT formula of clause density $\alpha$ is with high probability satisfiable for $\alpha<\alpha_*$, and unsatisfiable for $\alpha>\alpha_*$. We show that the threshold $\alpha_*(k)$ is given explicitly by the one-step replica symmetry breaking prediction from statistical physics. The proof develops a new analytic method for moment calculations on random graphs, mapping a high-dimensional optimization problem to a more tractable problem of analyzing tree recursions. We believe that our method may apply to a range of random CSPs in the 1-RSB universality class.


title: Semistability and Simple Connectivity at Infinity of Finitely Generated Groups with a Finite Series of Commensurated Subgroups

abstract: A subgroup $H$ of a group $G$ is $commensurated$ in $G$ if for each $g\in G$, $gHg^{-1}\cap H$ has finite index in both $H$ and $gHg^{-1}$. If there is a sequence of subgroups $H=Q_0\prec Q_1\prec ...\prec Q_{k}\prec Q_{k+1}=G$ where $Q_i$ is commensurated in $Q_{i+1}$ for all $i$, then $Q_0$ is $subcommensurated$ in $G$. In this paper we introduce the notion of the simple connectivity at infinity of a finitely generated group (in analogy with that for finitely presented groups). Our main result is: If a finitely generated group $G$ contains an infinite, finitely generated, subcommensurated subgroup $H$, of infinite index in $G$, then $G$ is 1-ended and semistable at $\infty$. If additionally, $H$ is finitely presented and 1-ended, then $G$ is simply connected at $\infty$. A normal subgroup of a group is commensurated, so this result is a strict generalization of a number of results, including the main theorems of G. Conner and M. Mihalik \cite{CM}, B. Jackson \cite{J}, V. M. Lew \cite{L}, M. Mihalik \cite{M1}and \cite{M2}, and J. Profio \cite{P}.


title: Clustering memes in social media streams

abstract: The problem of clustering content in social media has pervasive applications, including the identification of discussion topics, event detection, and content recommendation. Here we describe a streaming framework for online detection and clustering of memes in social media, specifically Twitter. A pre-clustering procedure, namely protomeme detection, first isolates atomic tokens of information carried by the tweets. Protomemes are thereafter aggregated, based on multiple similarity measures, to obtain memes as cohesive groups of tweets reflecting actual concepts or topics of discussion. The clustering algorithm takes into account various dimensions of the data and metadata, including natural language, the social network, and the patterns of information diffusion. As a result, our system can build clusters of semantically, structurally, and topically related tweets. The clustering process is based on a variant of Online K-means that incorporates a memory mechanism, used to "forget" old memes and replace them over time with the new ones. The evaluation of our framework is carried out by using a dataset of Twitter trending topics. Over a one-week period, we systematically determined whether our algorithm was able to recover the trending hashtags. We show that the proposed method outperforms baseline algorithms that only use content features, as well as a state-of-the-art event detection method that assumes full knowledge of the underlying follower network. We finally show that our online learning framework is flexible, due to its independence of the adopted clustering algorithm, and best suited to work in a streaming scenario.


title: Afterpulsing studies of low noise InGaAs/InP single-photon negative feedback avalanche diodes

abstract: We characterize the temporal evolution of the afterpulse probability in a free-running negative feedback avalanche diode (NFAD) over an extended range, from $\sim$300 ns to $\sim$1 ms. This is possible thanks to an extremely low dark count rate on the order of 1 cps at 10% efficiency, achieved by operating the NFAD at a temperatures as low as 143 K. Experimental results in a large range of operating temperatures (223-143 K) are compared with a legacy afterpulsing model based on multiple trap families at discrete energy levels, which is found to be lacking in physical completeness. Subsequently, we expand on a recent proposal which considers a continuous spectrum of traps by introducing well defined edges to the spectrum, which are experimentally observed.


title: Combining Technical and Financial Impacts for Countermeasure Selection

abstract: Research in information security has generally focused on providing a comprehensive interpretation of threats, vulnerabilities, and attacks, in particular to evaluate their danger and prioritize responses accordingly. Most of the current approaches propose advanced techniques to detect intrusions and complex attacks but few of these approaches propose well defined methodologies to react against a given attack. In this paper, we propose a novel and systematic method to select security countermeasures from a pool of candidates, by ranking them based on the technical and financial impact associated to each alternative. The method includes industrial evaluation and simulations of the impact associated to a given security measure which allows to compute the return on response investment for different candidates. A simple case study is proposed at the end of the paper to show the applicability of the model.


title: The density of states from first principles

abstract: We present a novel algorithm to compute the density of states, which is proven to converge to the correct result. The algorithm is very general and can be applied to a wide range of models, in the frameworks of Statistical Mechanics and Lattice Gauge Theory. All the thermal or quantum expectation values can then be obtained by a simple integration of the density of states. As an application, a numerical study of 4d U(1) compact lattice gauge theory is presented.


title: Exotic G_2-manifolds

abstract: We exhibit the first examples of closed 7-dimensional Riemannian manifolds with holonomy G_2 that are homeomorphic but not diffeomorphic. These are also the first examples of closed Ricci-flat manifolds that are homeomorphic but not diffeomorphic. The examples are generated by applying the twisted connected sum construction to Fano 3-folds of Picard rank 1 and 2. The smooth structures are distinguished by the generalised Eells-Kuiper invariant introduced by the authors in arXiv:1406.2226.


title: Giant Gain Enhancement in Photonic Crystals with a Degenerate Band Edge

abstract: We propose a new approach leading to giant gain enhancement. It is based on unconventional slow wave resonance associated to a degenerate band edge (DBE) in the dispersion diagram for a special class of photonic crystals supporting two modes at each frequency. We show that the gain enhancement in a Fabry-Perot cavity (FPC) when operating at the DBE is several orders of magnitude stronger when compared to a cavity of the same length made of a standard photonic crystal with a regular band edge (RBE). The giant gain condition is explained by a significant increase in the photon lifetime and in the local density of states. We have demonstrated the existence of DBE operated special cavities that provide for superior gain conditions for solid-state lasers, quantum cascade lasers, traveling wave tubes, and distributed solid state amplifiers. We also report the possibility to achieve low-threshold lasing in FPC with DBE compared to RBE-based lasers.


title: Scaling Limit for the Kernel of the Spectral Projector and Remainder Estimates in the Pointwise Weyl Law

abstract: Let (M, g) be a compact smooth Riemannian manifold. We obtain new off-diagonal estimates as {\lambda} tend to infinity for the remainder in the pointwise Weyl Law for the kernel of the spectral projector of the Laplacian onto functions with frequency at most {\lambda}. A corollary is that, when rescaled around a non self-focal point, the kernel of the spectral projector onto the frequency interval (\lambda, \lambda + 1] has a universal scaling limit as {\lambda} goes to infinity (depending only on the dimension of M). Our results also imply that if M has no conjugate points, then immersions of M into Euclidean space by an orthonormal basis of eigenfunctions with frequencies in (\lambda, \lambda + 1] are embeddings for all {\lambda} sufficiently large.


title: Approximate Counting in SMT and Value Estimation for Probabilistic Programs

abstract: #SMT, or model counting for logical theories, is a well-known hard problem that generalizes such tasks as counting the number of satisfying assignments to a Boolean formula and computing the volume of a polytope. In the realm of satisfiability modulo theories (SMT) there is a growing need for model counting solvers, coming from several application domains (quantitative information flow, static analysis of probabilistic programs). In this paper, we show a reduction from an approximate version of #SMT to SMT. We focus on the theories of integer arithmetic and linear real arithmetic. We propose model counting algorithms that provide approximate solutions with formal bounds on the approximation error. They run in polynomial time and make a polynomial number of queries to the SMT solver for the underlying theory, exploiting "for free" the sophisticated heuristics implemented within modern SMT solvers. We have implemented the algorithms and used them to solve the value problem for a model of loop-free probabilistic programs with nondeterminism.


title: Many-body localization edge in the random-field Heisenberg chain

abstract: We present a large scale exact diagonalization study of the one dimensional spin $1/2$ Heisenberg model in a random magnetic field. In order to access properties at varying energy densities across the entire spectrum for system sizes up to $L=22$ spins, we use a spectral transformation which can be applied in a massively parallel fashion. Our results allow for an energy-resolved interpretation of the many body localization transition including the existence of an extensive many-body mobility edge. The ergodic phase is well characterized by Gaussian orthogonal ensemble statistics, volume-law entanglement, and a full delocalization in the Hilbert space. Conversely, the localized regime displays Poisson statistics, area-law entanglement and non ergodicity in the Hilbert space where a true localization never occurs. We perform finite size scaling to extract the critical edge and exponent of the localization length divergence.


title: ATCA-based ATLAS FTK input interface system

abstract: The first stage of the ATLAS Fast TracKer (FTK) is an ATCA-based input interface system, where hits from the entire silicon tracker are clustered and organized into overlapping eta-phi trigger towers before being sent to the tracking engines. First, FTK Input Mezzanine cards receive hit data and perform clustering to reduce data volume. Then, the ATCA-based Data Formatter system will organize the trigger tower data, sharing data among boards over full mesh backplanes and optic fibers. The board and system level design concepts and implementation details, as well as the operation experiences from the FTK full-chain testing, will be presented.


title: Broadband Energy Harvesting Using a Metamaterial Resonator Embedded With Non-Foster Impedance Circuitry

abstract: Radio Frequency Identification (RFID) and implantable biomedical devices need efficient power and data transfer with very low profile antennas. We propose a low profile electrically small antenna for near-field wireless power and data telemetry employing a metamaterial Split Ring Resonator (SRR) antenna. SRRs can be designed for operation over wide frequencies from RF to visible. However, they are inherently narrowband making them sensitive to component mismatch with respect to external transmit antenna. Here we propose an embedding of a non-foster impedance circuitry into the metamaterial SRR structure that imparts conjugate negative complex impedance to this resonator antenna thereby increasing the effective bandwidth and thus overcoming the fundamental limit for efficient signal coupling. We demonstrate the concept through extensive numerical simulations and a prototype system at the board level using discrete off-the-shelf components and printed circuit SRR antenna at 500 MHz. We show that the power transfer between SRR receive antenna and the external transmit loop antenna is broadened by almost 400 MHz which corresponds to increase in {\Delta}f/fC from 0.49 to 1.65, before and after non-foster circuit activation.


title: Realizing the Willis equations with pre-stresses

abstract: This paper proves that the linear elastic behavior of the material with inhomogeneous pre-stresses can be described by the Willis equations. In this case, the additional terms in the Willis equations, compared with the classical linear elastic equations for homogeneous media, are related to the gradient of pre-stresses. In this way, the material length scale is naturally incorporated in the framework of continuum mechanics. All these findings also coincide with the results of transformation elastodynamics, so that they can meet the requirement of the principle of material objectivity and the principle of general invariance.


title: Improving the sensitivity of stop searches with on-shell constrained invariant mass variables

abstract: The search for light stops is of paramount importance, both in general as a promising path to the discovery of beyond the standard model physics and more specifically as a way of evaluating the success of the naturalness paradigm. While the LHC experiments have ruled out much of the relevant parameter space, there are "stop gaps", i.e., values of sparticle masses for which existing LHC analyses have relatively little sensitivity to light stops. We point out that techniques involving on-shell constrained M_2 variables can do much to enhance sensitivity in this region and hence help close the stop gaps. We demonstrate the use of these variables for several benchmark points and describe the effect of realistic complications, such as detector effects and combinatorial backgrounds, in order to provide a useful toolkit for light stop searches in particular, and new physics searches at the LHC in general.


title: Building a Better Boosted Top Tagger

abstract: Distinguishing hadronically decaying boosted top quarks from massive QCD jets is an important challenge at the Large Hadron Collider. In this paper we use the power counting method to study jet substructure observables designed for top tagging, and gain insight into their performance. We introduce a powerful new family of discriminants formed from the energy correlation functions which outperform the widely used N-subjettiness. These observables take a highly non-trivial form, demonstrating the importance of a systematic approach to their construction.


title: Entropy Production in Relativistic Jet Boundary Layers

abstract: Hot relativistic jets, passing through a background medium with a pressure gradient p \propto r^{-\eta} where 2 < \eta <= 8/3, develop a shocked boundary layer containing a significant fraction of the jet power. In previous work, we developed a self-similar description of the boundary layer assuming isentropic flow, but we found that such models respect global energy conservation only for the special case \eta = 8/3. Here we demonstrate that models with \eta < 8/3 can be made self-consistent if we relax the assumption of constant specific entropy. Instead, the entropy must increase with increasing r along the boundary layer, presumably due to multiple shocks driven into the flow as it gradually collimates. The increase in specific entropy slows the acceleration rate of the flow and provides a source of internal energy that could be channeled into radiation. We suggest that this process may be important for determining the radiative characteristics of tidal disruption events and gamma-ray bursts from collapsars.


title: The Formation of Massive, Compact Galaxies at z=2 in the Illustris Simulation

abstract: Massive, quiescent galaxies at high redshift have been found to be considerably more compact than galaxies of similar mass in the local universe. How these compact galaxies formed has yet to be determined, though several progenitor populations have been proposed. Here we investigate the formation processes and quantify the assembly histories of such galaxies in Illustris, a suite of hydrodynamical cosmological simulations encompassing a sufficiently large volume to include rare objects, while simultaneously resolving the internal structure of galaxies. We select massive (~10^11 solar masses) and compact (stellar half-mass radius < 2 kpc) galaxies from the simulation at z=2. Within the Illustris suite, we find that these quantities are not perfectly converged, but are reasonably reliable for our purposes. The resulting population is composed primarily of quiescent galaxies, but we also find several star-forming compact galaxies. The simulated compact galaxies are similar to observed galaxies in star formation activity and appearance. We follow their evolution at high redshift in the simulation and find that there are multiple pathways to form these compact galaxies, dominated by two mechanisms: (i) intense, centrally concentrated starbursts generally triggered by gas-rich major mergers between z~2-4, reducing the galaxies' half-mass radii by a factor of a few to below 2 kpc, and (ii) assembly at very early times when the universe was much denser; the galaxies formed compact and remained so until z~2.


title: Invisible dark gauge boson search in top decays using a kinematic method

abstract: We discuss the discovery potential of a dark force carrier ($Z'$) of very light mass, $m_{Z'} \lesssim {\cal O}(1-10)$ GeV, at hadron colliders via rare top quark decays, especially when it decays invisibly in typical search schemes. We emphasize that the top sector is promising for the discovery of new particles because top quark pairs are copiously produced at the Large Hadron Collider. The signal process is initiated by a rare top decay into a bottom quark and a charged Higgs boson ($H^\pm$) decaying subsequently into a $W$ and one or multiple $Z'$s. The light $Z'$ can be invisible in collider searches in various scenarios, and it would be hard to distinguish the relevant collider signature from the regular $t\bar{t}$ process in the Standard Model. We suggest a search strategy using the recently proposed on-shell constrained $M_2$ variables. Our signal process is featured by an $\textit{asymmetric}$ event topology, while the $t\bar{t}$ is $\textit{symmetric}$. The essence behind the strategy is to evoke some contradiction in the relevant observables by applying the kinematic variables designed under the assumption of the $t\bar{t}$ event topology. To see the viability of the proposed technique, we perform Monte Carlo simulations including realistic effects such as cuts, backgrounds, detector resolution, and so on at the LHC of $\sqrt{s}=14$ TeV.


title: Model-independent precision constraints on dimension-6 operators

abstract: We discuss electroweak precision constraints on dimension-6 operators in the effective theory beyond the standard model. We identify the combinations of these operators that are constrained by the pole observables (the W and Z masses and on-shell decays) and by the W boson pair production. To this end, we define a set of effective couplings of W and Z bosons to fermions and to itself, which capture the effects of new physics corrections. This formalism clarifies which operators are constrained by which observable, independently of the adopted basis of operators. We obtain numerical constraints on the coefficients of dimension-6 operator in a form that can be easily adapted to any particular basis of operators, or any particular model with new heavy particles.


title: NuSTAR Spectroscopy of Multi-Component X-ray Reflection from NGC 1068

abstract: We report on observations of NGC1068 with NuSTAR, which provide the best constraints to date on its $>10$~keV spectral shape. We find no strong variability over the past two decades, consistent with its Compton-thick AGN classification. The combined NuSTAR, Chandra, XMM-Newton, and Swift-BAT spectral dataset offers new insights into the complex reflected emission. The critical combination of the high signal-to-noise NuSTAR data and a spatial decomposition with Chandra allow us to break several model degeneracies and greatly aid physical interpretation. When modeled as a monolithic (i.e., a single N_H) reflector, none of the common Compton-reflection models are able to match the neutral fluorescence lines and broad spectral shape of the Compton reflection. A multi-component reflector with three distinct column densities (e.g., N_H~1.5e23, 5e24, and 1e25 cm^{-2}) provides a more reasonable fit to the spectral lines and Compton hump, with near-solar Fe abundances. In this model, the higher N_H components provide the bulk of the Compton hump flux while the lower N_H component produces much of the line emission, effectively decoupling two key features of Compton reflection. We note that ~30% of the neutral Fe Kalpha line flux arises from >2" (~140 pc), implying that a significant fraction of the <10 keV reflected component arises from regions well outside of a parsec-scale torus. These results likely have ramifications for the interpretation of poorer signal-to-noise observations and/or more distant objects [Abridged].


title: SEEDS Adaptive Optics Imaging of the Asymmetric Transition Disk Oph IRS 48 in Scattered Light

abstract: We present the first resolved near infrared imagery of the transition disk Oph IRS 48 (WLY 2-48), which was recently observed with ALMA to have a strongly asymmetric sub-millimeter flux distribution. H-band polarized intensity images show a $\sim$60AU radius scattered light cavity with two pronounced arcs of emission, one from Northeast to Southeast and one smaller, fainter and more distant arc in the Northwest. K-band scattered light imagery reveals a similar morphology, but with a clear third arc along the Southwestern rim of the disk cavity. This arc meets the Northwestern arc at nearly a right angle, revealing the presence of a spiral arm or local surface brightness deficit in the disk, and explaining the East-West brightness asymmetry in the H-band data. We also present 0.8-5.4$\mu$m IRTF SpeX spectra of this object, which allow us to constrain the spectral class to A0$\pm$1 and measure a low mass accretion rate of 10$^{-8.5}$M$_{\odot}$/yr, both consistent with previous estimates. We investigate a variety of reddening laws in order to fit the mutliwavelength SED of Oph IRS 48 and find a best fit consistent with a younger, higher luminosity star than previous estimates.