title: Quantitative Constraints on the Reionization History from the IGM Damping Wing Signature in Two Quasars at z > 7

abstract: During reionization, neutral hydrogen in the intergalactic medium (IGM) imprints a damping wing absorption feature on the spectrum of high-redshift quasars. A detection of this signature provides compelling evidence for a significantly neutral Universe, and enables measurements of the hydrogen neutral fraction $x_{\rm HI}(z)$ at that epoch. Obtaining reliable quantitative constraints from this technique, however, is challenging due to stochasticity induced by the patchy inside-out topology of reionization, degeneracies with quasar lifetime, and the unknown unabsorbed quasar spectrum close to rest-frame Ly$\alpha$. We combine a large-volume semi-numerical simulation of reionization topology with 1D radiative transfer through high-resolution hydrodynamical simulations of the high-redshift Universe to construct models of quasar transmission spectra during reionization. Our state-of-the-art approach captures the distribution of damping wing strengths in biased quasar halos that should have reionized earlier, as well as the erosion of neutral gas in the quasar environment caused by its own ionizing radiation. Combining this detailed model with our new technique for predicting the quasar continuum and its associated uncertainty, we introduce a Bayesian statistical method to jointly constrain the neutral fraction of the Universe and the quasar lifetime from individual quasar spectra. We apply this methodology to the spectra of the two highest redshift quasars known, ULAS J1120+0641 and ULAS J1342+0928, and measured volume-averaged neutral fractions $\langle x_{\rm HI} \rangle(z=7.09)=0.48^{+0.26}_{-0.26}$ and $\langle x_{\rm HI} \rangle(z=7.54)=0.60^{+0.20}_{-0.23}$ (posterior medians and 68% credible intervals) when marginalized over quasar lifetimes of $10^3 \leq t_{\rm q} \leq 10^8$ years.


title: Algorithmic improvements for the CIECAM02 and CAM16 color appearance models

abstract: This note is concerned with the CIECAM02 color appearance model and its successor, the CAM16 color appearance model. Several algorithmic flaws are pointed out and remedies are suggested. The resulting color model is algebraically equivalent to CIECAM02/CAM16, but shorter, more efficient, and works correctly for all edge cases.


title: Artificial intelligence and pediatrics: A synthetic mini review

abstract: The use of artificial intelligence intelligencein medicine can be traced back to 1968 when Paycha published his paper Le diagnostic a l'aide d'intelligences artificielle, presentation de la premiere machine diagnostri. Few years later Shortliffe et al. presented an expert system named Mycin which was able to identify bacteria causing severe blood infections and to recommend antibiotics. Despite the fact that Mycin outperformed members of the Stanford medical school in the reliability of diagnosis it was never used in practice due to a legal issue who do you sue if it gives a wrong diagnosis?. However only in 2016 when the artificial intelligence software built into the IBM Watson AI platform correctly diagnosed and proposed an effective treatment for a 60-year-old womans rare form of leukemia the AI use in medicine become really popular.On of first papers presenting the use of AI in paediatrics was published in 1984. The paper introduced a computer-assisted medical decision making system called SHELP.


title: Inflation in a scale invariant universe

abstract: A scale-invariant universe can have a period of accelerated expansion at early times: inflation. We use a frame-invariant approach to calculate inflationary observables in a scale invariant theory of gravity involving two scalar fields - the spectral indices, the tensor to scalar ratio, the level of isocurvature modes and non-Gaussianity. We show that scale symmetry leads to an exact cancellation of isocurvature modes and that, in the scale-symmetry broken phase, this theory is well described by a single scalar field theory. We find the predictions of this theory strongly compatible with current observations.


title: Diversity is All You Need: Learning Skills without a Reward Function

abstract: Intelligent creatures can explore their environments and learn useful skills without supervision. In this paper, we propose DIAYN ('Diversity is All You Need'), a method for learning useful skills without a reward function. Our proposed method learns skills by maximizing an information theoretic objective using a maximum entropy policy. On a variety of simulated robotic tasks, we show that this simple objective results in the unsupervised emergence of diverse skills, such as walking and jumping. In a number of reinforcement learning benchmark environments, our method is able to learn a skill that solves the benchmark task despite never receiving the true task reward. We show how pretrained skills can provide a good parameter initialization for downstream tasks, and can be composed hierarchically to solve complex, sparse reward tasks. Our results suggest that unsupervised discovery of skills can serve as an effective pretraining mechanism for overcoming challenges of exploration and data efficiency in reinforcement learning.


title: The origins and legacy of Kolmogorov's Grundbegriffe

abstract: April 25, 2003, marked the 100th anniversary of the birth of Andrei Nikolaevich Kolmogorov, the twentieth century's foremost contributor to the mathematical and philosophical foundations of probability. The year 2003 was also the 70th anniversary of the publication of Kolmogorov's Grundbegriffe der Wahrscheinlichkeitsrechnung. Kolmogorov's Grundbegriffe put probability's modern mathematical formalism in place. It also provided a philosophy of probability - an explanation of how the formalism can be connected to the world of experience. In this article, we examine the sources of these two aspects of the Grundbegriffe - the work of the earlier scholars whose ideas Kolmogorov synthesized.


title: 3-D Volumetric Gamma-ray Imaging and Source Localization with a Mobile Robot

abstract: Radiation detection has largely been a manual inspection process with point sensors such as Geiger-Muller counters and scintillation spectrometers to date. While their observations of source proximity prove useful, they lack the directional information necessary for efficient source localization and characterization in cluttered environments with multiple radiation sources. The recent commercialization of Compton gamma cameras provides directional information to the broader radiation detection community for the first time. This paper presents the integration of a Compton gamma camera with a self-localizing ground robot for accurate 3D radiation mapping. Using the position and orientation of the robot, radiation images from the gamma camera are accumulated over a traversed path in a shared frame of reference to construct a consistent voxel grid-based radiation map. The peaks of the map at pre-specified energy windows are selected as the source location estimates, which are compared to the ground truth source locations. The proposed approach localizes multiple sources to within an average of 0.2 m in two 5 x 4 m^2 and 14 x 6 m^2 laboratory environments.


title: An Introduction to Schur Polynomials

abstract: Notes from a course at the ATM Workshop on Schubert Varieties, held at The Institute of Mathematical Sciences, Chennai, in November 2017. Various expansions of Schur functions, the Lindstr\"om-Gessel-Viennot lemma, semistandard Young tableaux, Schensted's insertion algorithm, the plactic monoid, the RSK correspondence, and the Littlewood-Richardson rule are discussed.


title: Direct Calculation of the Transfer Map of Electrostatic Deflectors, and Comparison with the Codes COSY INFINITY and GIOS

abstract: COSY INFINITY uses a beamline coordinate system with a Frenet-Serret frame relative to the reference particle, and calculates differential algebra-valued transfer maps by integrating the ODEs of motion in the respective vector space over a differential algebra (DA). We will describe and perform computation of the DA transfer map of an electrostatic spherical deflector in a laboratory coordinate system using two conventional methods: (1) by integrating the ODEs of motion using a 4th order Runge-Kutta integrator and (2) by computing analytically and in closed form the properties of the respective elliptical orbits from Kepler theory. We will compare the resulting transfer maps with (3) the DA transfer map of COSY INFINITY's built-in electrostatic spherical deflector element ESP and (4) the transfer map of the electrostatic spherical deflector computed using the program GIOS. In addition to the electrostatic spherical deflector, we study an electrostatic cylindrical deflector, where the Kepler theory is not applicable. We compute the DA transfer map by the ODE integration method (1), and compare it with the transfer maps by (3) COSY INFINITY's built-in electrostatic cylindrical deflector element ECL, and (4) GIOS. In addition to the code listings in the appendices, the codes to run the test cases are available at http://bt.pa.msu.edu/cgi-bin/display.pl?name=ELSPHTM17


title: Resurrecting the Power-law, Intermediate, and Logamediate Inflations in the DBI Scenario with Constant Sound Speed

abstract: We investigate the power-law, intermediate, and logamediate inflationary models in the framework of DBI non-canonical scalar field with constant sound speed. In the DBI setting, we first represent the power spectrum of both scalar density and tensor gravitational perturbations. Then, we derive different inflationary observables including the scalar spectral index $n_s$, the running of the scalar spectral index $dn_s/d\ln k$, and the tensor-to-scalar ratio $r$. We show that the 95\% CL constraint of the Planck 2015 T+E data on the non-Gaussianity parameter $f_{{\rm NL}}^{{\rm DBI}}$ leads to the sound speed bound $c_{s}\geq0.087$ in the DBI inflation. Moreover, our results imply that, although the predictions of the power-law, intermediate, and logamediate inflations in the standard canonical framework ($c_s=1$) are not consistent with the Planck 2015 data, in the DBI scenario with constant sound speed $c_s<1$, the result of the $r-n_s$ diagram for these models can lie inside the 68\% CL region favored by the Planck 2015 TT,TE,EE+lowP data. We also specify the parameter space of the power-law, intermediate, and logamediate inflations for which our models are compatible with the 68\% or 95\% CL regions of the Planck 2015 TT,TE,EE+lowP data. Using the allowed ranges of the parameter space of the intermediate and logamediate inflationary models, we estimate the running of the scalar spectral index and find that it is compatible with the 95\% CL constraint from the Planck 2015 TT,TE,EE+lowP data.


title: Gamma-ray intensities in multi-gated spectra

abstract: The level structure of nuclei offers a large amount and variety of information to improve our knowledge of the strong interaction and of mesoscopic quantum systems. Gamma spectroscopy is a powerful tool to perform such studies: modern gamma multi-detectors present increasing performances in terms of sensitivity and efficiency, allowing to extend ever more our ability to observe and characterize abundant nuclear states. For instance, the high-spin part of level schemes often reflects intriguing nuclear shape phenomena: this behaviour is unveiled by high-fold experimental data analysed through multi-coincidence spectra, in which long deexcitation cascades become observable. Determining the intensity of newly discovered transitions is important to characterize the nuclear structure and formation mechanism of the emitting levels. However, it is not trivial to relate the apparent intensity observed in multi-gated spectra to the actual transition intensity. In this work, we introduce the basis of a formalism affiliated with graph theory: we have obtained analytic expressions from which data-analysis methods can eventually be derived to recover this link in a rigorous way.


title: A sampling-based approximation of the complex error function and its implementation without poles

abstract: Recently we developed a new sampling methodology based on incomplete cosine expansion of the sinc function and applied it in numerical integration in order to obtain a rational approximation for the complex error function $w\left(z \right) = e^{- {z^2}}\left(1 + \frac{2i}{\sqrt \pi}\int_0^z e^{t^2}dt\right),$ where $z = x + iy$. As a further development, in this work we show how this sampling-based rational approximation can be transformed into alternative form for efficient computation of the complex error function $w\left(z \right)$ at smaller values of the imaginary argument $y=\operatorname{Im}\left[z \right]$. Such an approach enables us to avoid poles in implementation and to cover the entire complex plain with high accuracy in a rapid algorithm. An optimized Matlab code utilizing only three rapid approximations is presented.


title: Decay Modes of the Hoyle State in $^{12}C$

abstract: Recent experimental results give an upper limit less than 0.043\% (95\% C.L.) to the direct decay of the Hoyle state into 3$\alpha$ respect to the sequential decay into $^8${Be}+$\alpha$. We performed one and two-dimensional tunneling calculations to estimate such a ratio and found it to be more than one order of magnitude smaller than experiment depending on the range of the nuclear force. This is within high statistics experimental capabilities. Our results can also be tested by measuring the decay modes of high excitation energy states of $^{12}$C where the ratio of direct to sequential decay might reach 10\% at $E^*$($^{12}$C)=10.3 MeV. The link between a Bose Einstein Condensate (BEC) and the direct decay of the Hoyle state is also addressed. We discuss a hypothetical `Efimov state' at $E^*$($^{12}$C)=7.458 MeV, which would mainly {\it sequentially} decay with 3$\alpha$ of {\it equal energies}: a counterintuitive result of tunneling. Such a state, if it would exist, is at least 8 orders of magnitude less probable than the Hoyle's, thus below the sensitivity of recent and past experiments.


title: Global-scale phylogenetic linguistic inference from lexical resources

abstract: Automatic phylogenetic inference plays an increasingly important role in computational historical linguistics. Most pertinent work is currently based on expert cognate judgments. This limits the scope of this approach to a small number of well-studied language families. We used machine learning techniques to compile data suitable for phylogenetic inference from the ASJP database, a collection of almost 7,000 phonetically transcribed word lists over 40 concepts, covering two third of the extant world-wide linguistic diversity. First, we estimated Pointwise Mutual Information scores between sound classes using weighted sequence alignment and general-purpose optimization. From this we computed a dissimilarity matrix over all ASJP word lists. This matrix is suitable for distance-based phylogenetic inference. Second, we applied cognate clustering to the ASJP data, using supervised training of an SVM classifier on expert cognacy judgments. Third, we defined two types of binary characters, based on automatically inferred cognate classes and on sound-class occurrences. Several tests are reported demonstrating the suitability of these characters for character-based phylogenetic inference.


title: Possibility of Entanglement at LIGO!

abstract: It is shown that a linearized classical gravity wave $\hat{a}$ {\em la} Einstein can get entangled with an array of test masses in a plane perpendicular to its direction of propagation. A Bell-CHSH inequality based on the requirement of noncontextuality for classical realism is derived, and it is shown that the entangled state produced violates this inequality.


title: A cooling neutron star crust after recurrent outbursts: Modelling the accretion outburst history of Aql X-1

abstract: With our neutron star crust cooling code {\tt NSCool} we track the thermal evolution of the neutron star in Aql X-1 over the full accretion outburst history from 1996 until 2015. For the first time, we model many outbursts (23 outbursts were detected) collectively and in great detail. This allows us to investigate the influence of previous outbursts on the internal temperature evolution and to test different neutron star crust cooling scenarios. Aql X-1 is an ideal test source for this purpose, because it shows frequent, short outbursts and thermally dominated quiescence spectra. The source goes into outburst roughly once a year for a few months. Assuming that the quiescent {\it Swift}/XRT observations of Aql X-1 can be explained within the crust cooling scenario (Waterhouse et al. 2016), we find three main conclusions. Firstly, the data are well reproduced by our model if the envelope composition and shallow heating parameters are allowed to change between outbursts. This is not the case if both shallow heating parameters (strength and depth) are tied throughout all accretion episodes, supporting earlier results that the properties of the shallow heating mechanism are not constant between outbursts. Second, from our models shallow heating could not be connected to one specific spectral state during outburst. Third, and most importantly, we find that the neutron star in Aql X-1 does not have enough time between outbursts to cool down to crust-core equilibrium and that heating during one outburst influences the cooling curves of the next.


title: Studies of Lepton Flavor Violation at the LHC

abstract: We examine the charged lepton flavor violating process $gg \rightarrow \mu^\pm \tau^\mp$ at the $\sqrt{s} = 13$ TeV LHC. Operators generating this process can be induced by new physics (NP) at dimension 8. Despite the power suppression associated with dimension 8 operators, we show that the LHC's large gluon luminosity makes it possible to probe this channel. For an integrated luminosity of 100 fb$^{-1}$ at the LHC, we predict a constraint on the NP scale $\Lambda \gtrsim 3$ TeV. In addition, we point out that such operators can be induced through top quark loops in models that generate dimension 6 operators of the form $\overline{t} t \, \mu \tau$. We find that the NP scale of these dimension 6 operators can be constrained to be $\Lambda \gtrsim 3.4-4.1$ TeV with 100 fb$^{-1}$ of data.


title: ALMA view of a massive spheroid progenitor: a compact rotating core of molecular gas in an AGN host at z=2.226

abstract: We present ALMA observations at 107.291 GHz (band 3) and 214.532 GHz (band 6) of GMASS 0953, a star-forming galaxy at z=2.226 hosting an obscured AGN that has been proposed as a progenitor of compact quiescent galaxies (QG). We measure for the first time the size of the dust and molecular gas emission of GMASS 0953 that we find to be extremely compact ($\sim$1 kpc). This result, coupled with a very high ISM density (n$\sim$10$^{5.5}$ cm$^{-3}$), a low gas mass fraction ($\sim$0.2) and a short gas depletion timescale ($\sim$150 Myr) imply that GMASS 0953 is experiencing an episode of intense star-formation in its central region that will rapidly exhaust its gas reservoirs, likely aided by AGN-induced feedback, confirming its fate as a compact QG. Kinematic analysis of the CO(6-5) line shows evidence of rapidly-rotating gas ($V_{rot}$=320$^{+92}_{-53}$ km s$^{-1}$), as observed also in a handful of similar sources at the same redshift. On-going quenching mechanisms could either destroy the rotation or leave it intact leading the galaxy to evolve into a rotating QG.


title: The Influence of Galaxy Environment on the Stellar Initial Mass Function of Early-Type Galaxies

abstract: In this paper we investigate whether the stellar initial mass function of early-type galaxies depends on their host environment. To this purpose, we have selected a sample of early-type galaxies from the SPIDER catalogue, characterized their environment through the group catalogue of Wang et al. and used their optical SDSS spectra to constrain the IMF slope, through the analysis of IMF-sensitive spectral indices. To reach a high enough signal-to-noise ratio, we have stacked spectra in velocity dispersion ($\sigma_0$) bins, on top of separating the sample by galaxy hierarchy and host halo mass, as proxies for galaxy environment. In order to constrain the IMF, we have compared observed line strengths to predictions of MIUSCAT/EMILES synthetic stellar population models, with varying age, metallicity, and "bimodal" (low-mass tapered) IMF slope ($\rm \Gamma_b$). Consistent with previous studies, we find that $\rm \Gamma_b$ increases with $\sigma_0$, becoming bottom-heavy (i.e. an excess of low-mass stars with respect to the Milky-Way-like IMF) at high $\sigma_0$. We find that this result is robust against the set of isochrones used in the stellar population models, as well as the way the effect of elemental abundance ratios is taken into account. We thus conclude that it is possible to use currently state-of-the-art stellar population models and intermediate resolution spectra to consistently probe IMF variations. For the first time, we show that there is no dependence of $\Gamma_b$ on environment or galaxy hierarchy, as measured within the $3"$ SDSS fibre, thus leaving the IMF as an intrinsic galaxy property, possibly set already at high redshift.


title: Kaluza-Klein reductions and AdS/Ricci-flat correspondence

abstract: The AdS/Ricci-flat (AdS/RF) correspondence is a map between families of asymptotically locally AdS solutions on a torus and families of asymptotically flat spacetimes on a sphere. The aim of this work is to perturbatively extend this map to general AdS and asymptotically flat solutions. A prime application for such map would be the development of holography for Minkowski spacetime. In this paper we perform a Kaluza-Klein (KK) reduction of AdS on a torus and of Minkowski on a sphere, keeping all massive KK modes. Such computation is interesting on its own, as there are relatively few examples of such explicit KK reductions in the literature. We perform both KK reductions in parallel to illustrate their similarity. In particular, we show how to construct gauge invariant variables, find the field equations they satisfy, and construct a corresponding effective action. We further diagonalize all equations and find their general solution in closed form. Surprisingly, in the limit of large dimension of the compact manifolds (torus and sphere), the AdS/RF correspondence maps individual KK modes from one side to the other. In a sequel of this paper we will discuss how the AdS/RF maps acts when the dimension of the compact space is finite.


title: The convective photosphere of the red supergiant CE Tau. I. VLTI/PIONIER H-band interferometric imaging

abstract: Context. Red supergiant stars are one of the latest stages in the evolution of massive stars. Their photospheric convection may play an important role in the launching mechanism of their mass loss; however, its characteristics and dynamics are still poorly constrained. Aims. By observing red supergiant stars with near infrared interferometry at different epochs, we expect to reveal the evolution of bright convective features on their stellar surface. Methods. We observed the M2Iab-Ib red supergiant star CE Tau with the VLTI/PIONIER instrument in the H band at two different epochs separated by one month. Results. We derive the angular diameter of the star and basic stellar parameters, and reconstruct two reliable images of its H-band photosphere. The contrast of the convective pattern of the reconstructed images is $5 \pm 1 \%$ and $6 \pm 1 \%$ for our two epochs of observation. Conclusions. The stellar photosphere shows few changes between the two epochs. The contrast of the convective pattern is below the average contrast variations obtained on 30 randomly chosen snapshots of the best matching 3D radiative hydrodynamics simulation: $23 \pm 1~\%$ for the original simulation images and $16 \pm 1~\%$ for the maps degraded to the reconstruction resolution. We offer two hypotheses to explain this observation. CE Tau may be experiencing a quiet convective activity episode or it could be a consequence of its warmer effective temperature (hence its smaller radius) compared to the simulation.


title: Boson-fermion duality in a gravitational background

abstract: We study the $2+1$ dimensional boson-fermion duality in the presence of background curvature and electromagnetic fields. The main players are, on the one hand, a massive complex $|\phi|^4$ scalar field coupled to a $U(1)$ Maxwell-Chern-Simons gauge field at level $1$, representing a relativistic composite boson with one unit of attached flux, and on the other hand, a massive Dirac fermion. We show that, in a curved background and at the level of the partition function, the relativistic composite boson, in the infinite coupling limit, is dual to a short-range interacting Dirac fermion. The coupling to the gravitational spin connection arises naturally from the spin factors of the Wilson loop in the Chern-Simons theory. A non-minimal coupling to the scalar curvature is included on the bosonic side in order to obtain agreement between partition functions. Although an explicit Lagrangian expression for the fermionic interactions is not obtained, their short-range nature constrains them to be irrelevant, which protects the duality in its strong interpretation as an exact mapping at the IR fixed point between a Wilson-Fischer-Chern-Simons complex scalar and a free Dirac fermion. We also show that, even away from the IR, keeping the $|\phi|^4$ term is of key importance as it provides the short-range bosonic interactions necessary to prevent intersections of worldlines in the path integral, thus forbidding unknotting of knots and ensuring preservation of the worldline topologies.


title: Reducing and Analyzing the PHAT Survey with the Cloud

abstract: We discuss the technical challenges we faced and the techniques we used to overcome them when reducing the PHAT photometric data set on the Amazon Elastic Compute Cloud (EC2). We first describe the architecture of our photometry pipeline, which we found particularly efficient for reducing the data in multiple ways for different purposes. We then describe the features of EC2 that make this architecture both efficient to use and challenging to implement. We describe the techniques we adopted to process our data, and suggest ways these techniques may be improved for those interested in trying such reductions in the future. Finally, we summarize the output photometry data products, which are now hosted publicly in two places in two formats. They are in simple fits tables in the high-level science products on MAST, and on a queryable database available through the NOAO Data Lab.


title: Self-organization on Riemannian manifolds

abstract: We consider an aggregation model that consists of an active transport equation for the macroscopic population density, where the velocity has a nonlocal functional dependence on the density, modelled via an interaction potential. We set up the model on general Riemannian manifolds and provide a framework for constructing interaction potentials which lead to equilibria that are constant on their supports. We consider such potentials for two specific cases (the two-dimensional sphere and the two-dimensional hyperbolic space) and investigate analytically and numerically the long-time behaviour and equilibrium solutions of the aggregation model on these manifolds. Equilibria obtained numerically with other interaction potentials are also presented.


title: Unveiled electric profiles within hydrogen bonds suggest DNA base pairs with similar bond strengths

abstract: Electrical forces are the background of all the interactions occurring in biochemical systems. From here and by using a combination of ab-initio and ad-hoc models, we introduce the first description of electric field profiles with intrabond resolution to support a characterization of single bond forces attending to its electrical origin. This fundamental issue has eluded a physical description so far. Our method is applied to describe hydrogen bonds (HB) in DNA base pairs. Numerical results reveal that base pairs in DNA could be equivalent considering HB strength contributions, which challenges previous interpretations of thermodynamic properties of DNA based on the assumption that Adenine/Thymine pairs are weaker than Guanine/Cytosine pairs due to the sole difference in the number of HB. Thus, our methodology provides solid foundations to support the development of extended models intended to go deeper into the molecular mechanisms of DNA functioning.


title: Bridging Cognitive Programs and Machine Learning

abstract: While great advances are made in pattern recognition and machine learning, the successes of such fields remain restricted to narrow applications and seem to break down when training data is scarce, a shift in domain occurs, or when intelligent reasoning is required for rapid adaptation to new environments. In this work, we list several of the shortcomings of modern machine-learning solutions, specifically in the contexts of computer vision and in reinforcement learning and suggest directions to explore in order to try to ameliorate these weaknesses.


title: Four moments theorems on Markov chaos

abstract: We obtain quantitative Four Moments Theorems establishing convergence of the laws of elements of a Markov chaos to a Pearson distribution, where the only assumption we make on the Pearson distribution is that it admits four moments. While in general one cannot use moments to establish convergence to a heavy-tailed distributions, we provide a context in which only the first four moments suffices. These results are obtained by proving a general carr\'e du champ bound on the distance between laws of random variables in the domain of a Markov diffusion generator and invariant measures of diffusions. For elements of a Markov chaos, this bound can be reduced to just the first four moments.


title: Gradient descent with identity initialization efficiently learns positive definite linear transformations by deep residual networks

abstract: We analyze algorithms for approximating a function $f(x) = \Phi x$ mapping $\Re^d$ to $\Re^d$ using deep linear neural networks, i.e. that learn a function $h$ parameterized by matrices $\Theta_1,...,\Theta_L$ and defined by $h(x) = \Theta_L \Theta_{L-1} ... \Theta_1 x$. We focus on algorithms that learn through gradient descent on the population quadratic loss in the case that the distribution over the inputs is isotropic. We provide polynomial bounds on the number of iterations for gradient descent to approximate the least squares matrix $\Phi$, in the case where the initial hypothesis $\Theta_1 = ... = \Theta_L = I$ has excess loss bounded by a small enough constant. On the other hand, we show that gradient descent fails to converge for $\Phi$ whose distance from the identity is a larger constant, and we show that some forms of regularization toward the identity in each layer do not help. If $\Phi$ is symmetric positive definite, we show that an algorithm that initializes $\Theta_i = I$ learns an $\epsilon$-approximation of $f$ using a number of updates polynomial in $L$, the condition number of $\Phi$, and $\log(d/\epsilon)$. In contrast, we show that if the least squares matrix $\Phi$ is symmetric and has a negative eigenvalue, then all members of a class of algorithms that perform gradient descent with identity initialization, and optionally regularize toward the identity in each layer, fail to converge. We analyze an algorithm for the case that $\Phi$ satisfies $u^{\top} \Phi u > 0$ for all $u$, but may not be symmetric. This algorithm uses two regularizers: one that maintains the invariant $u^{\top} \Theta_L \Theta_{L-1} ... \Theta_1 u > 0$ for all $u$, and another that "balances" $\Theta_1, ..., \Theta_L$ so that they have the same singular values.

categories: cs.LG cs.NE math.OC math.ST stat.ML stat.TH


title: Examining the rank of Semi-definite Programming for Power System State Estimation

abstract: In the power system, state estimation (SE) is important monitoring task for the reliable operation of the system. The optimal estimate from the SE is delivered to all EMS application such as fault analysis, automatic generation control. Hence, it is crucial to have good estimation before taking any critical actions. However, the SE problem is challenging problem due to nonconvexity of power flow equations in the nonlinear AC power flow model, which give us a usually local solution. To deal with this nonconvexity, some recent literatures applied the convex semi-definite (SDP) relaxation technique to relax the SE problem attaining or approximating a global solution. In this paper, we investigate the rank of this technique, which is critical to yield a physically meaningful solution with the five-bus test system and propose new approach to possibly reduce the rank by complementing the traditional set of measurement with PMU data. Numerical tests on the standard IEEE 14, 30, 57, 118-bus test system are presented for the demonstration.


title: Mining Sub-Interval Relationships In Time Series Data

abstract: Time-series data is being increasingly collected and stud- ied in several areas such as neuroscience, climate science, transportation, and social media. Discovery of complex patterns of relationships between individual time-series, using data-driven approaches can improve our understanding of real-world systems. While traditional approaches typically study relationships between two entire time series, many interesting relationships in real-world applications exist in small sub-intervals of time while remaining absent or feeble during other sub-intervals. In this paper, we define the notion of a sub-interval relationship (SIR) to capture inter- actions between two time series that are prominent only in certain sub-intervals of time. We propose a novel and efficient approach to find most interesting SIR in a pair of time series. We evaluate our proposed approach on two real-world datasets from climate science and neuroscience domain and demonstrated the scalability and computational efficiency of our proposed approach. We further evaluated our discovered SIRs based on a randomization based procedure. Our results indicated the existence of several such relationships that are statistically significant, some of which were also found to have physical interpretation.


title: Rho-meson properties in medium

abstract: Properties of $\rho$-meson in symmetric nuclear matter are investigated in a light-front constituent quark model (LFCQM), using the in-medium inputs calculated by the quark-meson coupling (QMC) model. The LFCQM used in this study was already applied for the studies of the electromagnetic properties of $\rho$-meson in vacuum, namely, the charge~$G_0$, magnetic~$G_1$, and quadrupole~$G_2$ form factors, electromagnetic charge radius, and electromagnetic decay constant. We predict that the electromagnetic decay constant, charge radius, and quadrupole moment are enhanced as increasing the nuclear matter density, while the magnetic moment is slightly quenched. Furthermore, we predict that the value $Q^2_{\rm zero}$, which crosses zero of the charge form factor, $G_0(Q^2_{\rm zero})=0$ ($Q^2 = -q^2 > 0$ with $q$ being the four-momentum transfer), decreases as increasing the nuclear matter density.


title: Characterization of finite metric space by their isometric sequences

abstract: Let $(X,d)$ be a finite metric space with $|X|=n$. For a positive integer $k$ we define $A_k(X)$ to be the quotient set of all $k$-subsets of $X$ by isometry, and we denote $|A_k(X)|$ by $a_k$. The sequence $(a_1,a_2,\ldots,a_{n})$ is called the isometric sequence of $(X,d)$. In this article we aim to characterize finite metric spaces by their isometric sequences under one of the following assumptions: (i) $a_k=1$ for some $k$ with $2\leq k\leq n-2$; (ii) $a_k=2$ for some $k$ with $4\leq k\leq \frac{1+\sqrt{1+4n}}{2}$; (iii) $a_3=2$; (iv) $a_2=a_3=3$. Furthermore, we give some criterion on how to embed such finite metric spaces to Euclidean spaces. We give some maximum cardinalities of subsets in the $d$-dimensional Euclidean space with small $a_3$, which are analogue problems on a sets with few distinct triangles discussed by Epstein, Lott, Miller and Palsson.


title: Characterization of finite colored spaces with certain conditions

abstract: A colored space is the pair $(X,r)$ of a set $X$ and a function $r$ whose domain is $\binom{X}{2}$. Let $(X,r)$ be a finite colored space and $Y,Z\subseteq X$. We shall write $Y\simeq_r Z$ if there exists a bijection $f:Y\to Z$ such that $r(U)=r(f(U))$ for each $U\in\binom{Y}{2}$. We denote the numbers of equivalence classes with respect to $\simeq_r$ contained in $\binom{X}{2}$ and $\binom{X}{3}$ by $a_2(r)$ and $a_3(r)$, respectively. In this paper we prove that $a_2(r)\leq a_3(r)$ when $5\leq |X|$, and show what happens when the equality holds.


title: A problem in control of elastodynamics with piezoelectric effects

abstract: We consider an optimal control problem where the state equations are a coupled hyperbolic-elliptic system. This system arises in elastodynamics with piezoelectric effects -- the elastic stress tensor is a function of elastic displacement and electric potential. The electric flux acts as the control variable and, in addition to the state constraints, the bound constraints on the control are considered. We develop a complete analysis for the state equations and the control problem. The requisite regularity on the control, to show the well-posedness of state equations, is enforced using the cost functional. We rigorously derive the first order necessary and sufficient conditions using adjoint equations and further study their well-posedness. For spatially discrete (time continuous) problem, we show the convergence of our numerical scheme. Three dimensional numerical experiments are provided showing convergence properties of a fully discrete method and the practical applicability of our approach.


title: Extreme Value Analysis of Solar Flare Events

abstract: Space weather events such as solar flares can be harmful for life and infrastructure on earth or in near-earth orbit. In this paper we employ extreme value theory (EVT) to model extreme solar flare events; EVT offers the appropriate tools for the study and estimation of probabilities for extrapolation to ranges outside of those that have already been observed. In the past such phenomena have been modelled as following a power law which may gives poor estimates of such events due to overestimation. The data used in the study were X-ray fluxes from NOAA/GOES and the expected return levels for Carrington or Halloween like events were calculated with the outcome that the existing data predict similar events happening in 110 and 38 years respectively.

categories: stat.AP astro-ph.SR physics.data-an physics.space-ph


title: Market Impact in a Latent Order Book

abstract: The latent order book of \cite{donier2015fully} is one of the most promising agent-based models for market impact. This work extends the minimal model by allowing agents to exhibit mean-reversion, a commonly observed pattern in real markets. This modification leads to new order book dynamics, which we explicitly study and analyze. Underlying our analysis is a mean-field assumption that views the order book through its \textit{average} density. We show how price impact develops in this new model, providing a flexible family of solutions that can potentially improve calibration to real data. While no closed-form solution is provided, we complement our theoretical investigation with extensive numerical results, including a simulation scheme for the entire order book.


title: Mapping fast evolution of transient surface photovoltage dynamics using G-Mode Kelvin probe force microscopy

abstract: Optoelectronic phenomena in materials such as organic/inorganic hybrid perovskites depend on a complex interplay between light induced carrier generation and fast (electronic) and slower (ionic) processes, all of which are known to be strongly affected by structural inhomogeneities such as interfaces and grain boundaries. Here, we develop a time resolved Kelvin probe force microscopy (KPFM) approach, based on the G-Mode SPM platform, allowing quantification of surface photovoltage (SPV) with microsecond temporal and nanoscale spatial resolution. We demonstrate the approach on methylammonium lead bromide (MAPbBr3) thin films and further highlight the usefulness of unsupervised clustering methods to quickly discern spatial variability in the information rich SPV dataset. Using this technique, we observe concurrent spatial and ultra-fast temporal variations in the SPV generated across the thin film, indicating that structure is likely responsible for the heterogenous behavior.


title: Counting Homomorphisms to Trees Modulo a Prime

abstract: Many important graph theoretic notions can be encoded as counting graph homomorphism problems, such as partition functions in statistical physics, in particular, independent sets and colourings. In this article we study the complexity of $\#_p\mathsf{H}$OMS$\mathsf{T}$O$H$, the problem of counting graph homomorphisms from an input graph to a graph $H$ modulo a prime number $p$. Dyer and Greenhill proved a dichotomy stating that the tractability of non-modular counting graph homomorphisms depends on the structure of the input graph. Many intractable cases in non-modular counting become tractable in modular counting due to the common phenomenon of cancellation. However, in subsequent studies on counting modulo $2$ the influence, the structure of $H$ has on the tractability, was shown to persist, yielding similar dichotomies. Our main result shows that for every tree $H$ and every prime $p$ the problem $\#_p\mathsf{H}$OMS$\mathsf{T}$O$H$ is either polynomial time computable or $\#_p\mathsf{P}$-complete. This addresses the conjecture of Faben and Jerrum stating this dichotomy for every graph $H$ when counting modulo 2. In order to prove this result, we study the structural properties of a homomorphism. As an important interim, this study yields a dichotomy for the problem of weighted counting independent sets in a bipartite graph modulo some prime $p$. Our results are the first suggesting that such dichotomies hold not only for the one-bit functions of the modulo~2 case but for the modular counting functions of all primes $p$.


title: Information-theoretic Limits for Community Detection in Network Models

abstract: We analyze the information-theoretic limits for the recovery of node labels in several network models. This includes the Stochastic Block Model, the Exponential Random Graph Model, the Latent Space Model, the Directed Preferential Attachment Model, and the Directed Small-world Model. For the Stochastic Block Model, the non-recoverability condition depends on the probabilities of having edges inside a community, and between different communities. For the Latent Space Model, the non-recoverability condition depends on the dimension of the latent space, and how far and spread are the communities in the latent space. For the Directed Preferential Attachment Model and the Directed Small-world Model, the non-recoverability condition depends on the ratio between homophily and neighborhood size. We also consider dynamic versions of the Stochastic Block Model and the Latent Space Model.


title: Sub-tree counts on hyperbolic random geometric graphs

abstract: We study the hyperbolic random geometric graph introduced in Krioukov et al. For a sequence $R_n \to \infty$, we define these graphs to have the vertex set as Poisson points distributed uniformly in balls $B(0,R_n) \subset B_d^{\alpha}$, the $d$-dimensional Poincar\'e ball (unit d-ball with the Poincar\'e metric $d_{\alpha}$ corresponding to negative curvature $-\alpha^2, \alpha > 0$) by connecting any two points within a distance $R_n$ according to the metric $d_{\zeta}, \zeta > 0$. Denoting these graphs by $HG_n(R_n ; \alpha, \zeta)$, we study asymptotic counts of copies of a fixed tree $\Gamma_k$ (with the ordered degree sequence $d_{(1)} \leq \ldots \leq d_{(k)}$) in $HG_n(R_n ; \alpha, \zeta)$. Unlike earlier works, we count more involved structures, allowing for $d > 2$, and in many places, more general choices of $R_n$ rather than $R_n = 2[\zeta (d-1)]^{-1}\log (n/ \nu), \nu \in (0,\infty)$. The latter choice of $R_n$ for $\alpha / \zeta > 1/2$ corresponds to the thermodynamic regime. We show multiple phase transitions in $HG_n(R_n ; \alpha, \zeta)$ as $\alpha / \zeta$ increases, i.e., the space $B_d^{\alpha}$ becomes more hyperbolic. In particular, our analyses reveal that the sub-tree counts exhibit an intricate dependence on the degree sequence $d_{(1)},\ldots,d_{(k)}$ of $\Gamma_k$ as well as the ratio $\alpha/\zeta$. Under a more general radius regime $R_n$ than that described above, we investigate the asymptotics of the expectation and variance of sub-tree counts. Moreover, we prove the corresponding central limit theorem as well. Our proofs rely crucially on a careful analysis of the sub-tree counts near the boundary using Palm calculus for Poisson point processes along with estimates for the hyperbolic metric and measure. For the central limit theorem, we use the abstract normal approximation result from Last et al. derived using the Malliavin-Stein method.


title: A deep learning framework for turbulence modeling using data assimilation and feature extraction

abstract: Turbulent problems in industrial applications are predominantly solved using Reynolds Averaged Navier Stokes (RANS) turbulence models. The accuracy of the RANS models is limited due to closure assumptions that induce uncertainty into the RANS modeling. We propose the use of deep learning algorithms via convolution neural networks along with data from direct numerical simulations to extract the optimal set of features that explain the evolution of turbulent flow statistics. Statistical tests are used to determine the correlation of these features with the variation in the quantities of interest that are to be predicted. These features are then used to develop improved partial differential equations that can replace classical Reynolds Averaged Navier Stokes models and show improvement in the accuracy of the predictions.


title: Stable rationality in smooth families of threefolds

abstract: We exhibit families of smooth projective threefolds with both stably rational and non stably rational fibers.


title: Modeling the Formation of Social Conventions from Embodied Real-Time Interactions

abstract: What is the role of real-time control and learning in the formation of social conventions? To answer this question, we propose a computational model that matches human behavioral data in a social decision-making game that was analyzed both in discrete-time and continuous-time setups. Furthermore, unlike previous approaches, our model takes into account the role of sensorimotor control loops in embodied decision-making scenarios. For this purpose, we introduce the Control-based Reinforcement Learning (CRL) model. CRL is grounded in the Distributed Adaptive Control (DAC) theory of mind and brain, where low-level sensorimotor control is modulated through perceptual and behavioral learning in a layered structure. CRL follows these principles by implementing a feedback control loop handling the agent's reactive behaviors (pre-wired reflexes), along with an adaptive layer that uses reinforcement learning to maximize long-term reward. We test our model in a multi-agent game-theoretic task in which coordination must be achieved to find an optimal solution. We show that CRL is able to reach human-level performance on standard game-theoretic metrics such as efficiency in acquiring rewards and fairness in reward distribution.


title: Toeplitz algebras in quantum Hopf fibrations

abstract: The paper presents applications of Toeplitz algebras in Noncommutative Geometry. As an example, a quantum Hopf fibration is given by gluing trivial U(1) bundles over quantum discs (or, synonymously, Toeplitz algebras) along their boundaries. The construction yields associated quantum line bundles over the generic Podles spheres which are isomorphic to those from the well-known Hopf fibration of quantum SU(2). The relation between these two versions of quantum Hopf fibrations is made precise by giving an isomorphism in the category of right U(1)-comodules and left modules over the C*-algebra of the generic Podles spheres. It is argued that the gluing construction yields a significant simplification of index computations by obtaining elementary projections as representatives of K-theory classes.


title: Light, the universe, and everything -- 12 Herculean tasks for quantum cowboys and black diamond skiers

abstract: The Winter Colloquium on the Physics of Quantum Electronics (PQE) has been a seminal force in quantum optics and related areas since 1971. It is rather mindboggling to recognize how the concepts presented at these conferences have transformed scientific understanding and human society. In January, 2017, the participants of PQE were asked to consider the equally important prospects for the future, and to formulate a set of questions representing some of the greatest aspirations in this broad field. The result is this multi-authored paper, in which many of the world's leading experts address the following fundamental questions: (1) What is the future of gravitational wave astronomy? (2) Are there new quantum phases of matter away from equilibrium that can be found and exploited - such as the time crystal? (3) Quantum theory in uncharted territory: What can we learn? (4) What are the ultimate limits for laser photon energies? (5) What are the ultimate limits to temporal, spatial, and optical resolution? (6) What novel roles will atoms play in technology? (7) What applications lie ahead for nitrogen-vacancy centers in diamond? (8) What is the future of quantum coherence, squeezing, and entanglement for enhanced superresolution and sensing? (9) How can we solve (some of) humanity's biggest problems through new quantum technologies? (10) What new understanding of materials and biological molecules will result from their dynamical characterization with free electron lasers? (11) What new technologies and fundamental discoveries might quantum optics achieve by the end of this century? (12) What novel topological structures can be created and employed in quantum optics?


title: The small-x gluon distribution in centrality biased pA and pp collisions

abstract: The nuclear modification factor $R_{pA}(p_T)$ provides information on the small-$x$ gluon distribution of a nucleus at hadron colliders. Several experiments have recently measured the nuclear modification factor not only in minimum bias but also for central $pA$ collisions. In this paper we analyze the bias on the configurations of soft gluon fields introduced by a centrality selection via the number of hard particles. Such bias can be viewed as reweighting of configurations of small-$x$ gluons. We find that the biased nuclear modification factor ${\cal Q}_{pA}(p_T)$ for central collisions is above $R_{pA}(p_T)$ for minimum bias events, and that it may redevelop a "Cronin peak" even at small $x$. The magnitude of the peak is predicted to increase approximately like $1/{A_{\perp}}^\nu$, $\nu\sim0.6\pm0.1$, if one is able to select more compact configurations of the projectile proton where its gluons occupy a smaller transverse area $A_\perp$. We predict an enhanced ${\cal Q}_{pp}(p_T)-1 \sim 1/(p_T^2)^\nu$ and a Cronin peak even for central $pp$ collisions.


title: Affine quadrics and the Picard group of the motivic category

abstract: In this article we study the subgroup of the Picard group of Voevodsky's category of geometric motives generated by the reduced motives of affine quadrics. Our main tools here are the functors of Bachmann, but we also provide an alternative method. We show that the group in question can be described in terms of indecomposable direct summands in the motives of projective quadrics. In particular, we describe all the relations among the reduced motives of affine quadrics. We also extend the Criterion of motivic equivalence of projective quadrics.


title: Optics and bremsstrahlung estimates for channeling radiation experiments at FAST

abstract: This paper presents X-ray spectra of channeling radiation (CR) expected at the FAST (Fermi Accelerator Science and Technology) facility in Fermilab. Our purpose is to produce high brightness quasi-monochromatic X-rays in an energy range from 40 keV to 110 keV. We will use a diamond crystal and low emittance electrons with an energy of around 43 MeV. The quality of emitted X-rays strongly depends on parameters of the electron beam at the crystal. We present simulations of the beam optics for high brightness and high yield operations with bunch charges of 1 pC, 20 pC, and 200 pC. We estimate the X-ray spectra including bremsstrahlung background for a charge of 20 pC. The electron beam distributions with and without channeling in the diamond crystal are calculated. We discuss an X-ray detector system to avoid pile-up effect during high charge operations.


title: Restricting the bi-equivariant spectral triple on quantum SU(2) to the Podles spheres

abstract: It is shown that the isospectral bi-equivariant spectral triple on quantum SU(2) and the isospectral equivariant spectral triples on the Podles spheres are related by restriction. In this approach, the equatorial Podles sphere is distinguished because only in this case the restricted spectral triple admits an equivariant grading operator together with a real structure (up to infinitesimals of arbitrary high order). The real structure is expressed by the Tomita operator on quantum SU(2) and it is shown that the failure of the real structure to satisfy the commutant property is related to the failure of the universal R-matrix operator to be unitary.


title: On the dissociation between potential vorticity conservation and symmetries

abstract: Using a four-dimensional manifestly covariant formalism suitable for classical fluid dynamics, it is shown that the conservation of potential vorticity is not associated with any symmetry of the equations of motion but is instead a trivial conservation law of the second kind. The demonstration is provided in arbitrary coordinates and therefore applies to comoving (or label) coordinates. Since this is at odds with previous studies, which claimed that potential vorticity conservation is associated with a symmetry under particle-relabeling, a detailed discussion on relabeling transformations is also presented.


title: Evolution of H$_2$O Production in Comet C/2012 S1 (ISON) as Inferred from Forbidden Oxygen and OH Emission

abstract: We present H$_2$O production rates for comet C/2012 S1 (ISON) derived from observations of [OI] and OH emission during its inbound leg, covering a heliocentric distance range of 1.8-0.44 AU. Our production rates are in agreement with previous measurements using a variety of instruments and techniques and with data from the various observatories greatly differing in their projected fields of view. The consistent results across all data suggest the absence of an extended source of H$_2$O production, for example sublimation of icy grains in the coma, or a source with spatial extent confined to the dimensions of the smallest projected field of view (in this case $<$ 1,000 km). We find that ISON had an active area of around 10 km$^2$ for heliocentric distances R$_h$ > 1.2 AU, which then decreased to about half this value from R$_h$=1.2-0.9 AU. This was followed by a rapid increase in active area at about R$_h$=0.6 AU, corresponding to the first of three major outbursts ISON experienced inside of 1 AU. The combination of a detected outburst in the light curve and rapid increase in active area likely indicates a major nucleus fragmentation event. The 5-10 km$^2$ active area observed outside of R$_h$=0.6 AU is consistent with a 50-100% active fraction for the nucleus, larger than typically observed for cometary nuclei. Although the absolute value of the active area is somewhat dependent on the thermal model employed, the changes in active area observed are consistent among models. The conclusion of a 50-100+% active fraction is robust for realistic thermal models of the nucleus. However the possibility of a contribution of a spatially unresolved distribution of icy grains cannot be discounted.


title: Scenarios: A New Representation for Complex Scene Understanding

abstract: The ability for computational agents to reason about the high-level content of real world scene images is important for many applications. Existing attempts at addressing the problem of complex scene understanding lack representational power, efficiency, and the ability to create robust meta-knowledge about scenes. In this paper, we introduce scenarios as a new way of representing scenes. The scenario is a simple, low-dimensional, data-driven representation consisting of sets of frequently co-occurring objects and is useful for a wide range of scene understanding tasks. We learn scenarios from data using a novel matrix factorization method which we integrate into a new neural network architecture, the ScenarioNet. Using ScenarioNet, we can recover semantic information about real world scene images at three levels of granularity: 1) scene categories, 2) scenarios, and 3) objects. Training a single ScenarioNet model enables us to perform scene classification, scenario recognition, multi-object recognition, content-based scene image retrieval, and content-based image comparison. In addition to solving many tasks in a single, unified framework, ScenarioNet is more computationally efficient than other CNNs because it requires significantly fewer parameters while achieving similar performance on benchmark tasks and is more interpretable because it produces explanations when making decisions. We validate the utility of scenarios and ScenarioNet on a diverse set of scene understanding tasks on several benchmark datasets.


title: Tracking critical points on evolving curves and surfaces

abstract: In recent years it became apparent that geophysical abrasion can be well characterized by the time evolution $N(t)$ of the number $N$ of static balance points of the abrading particle. Static balance points correspond to the critical points of the particle's surface represented as a scalar distance function $r$, measured from the center of mass, so their time evolution can be expressed as $N(r(t))$. The mathematical model of the particle can be constructed on two scales: on the macro scale the particle may be viewed as a smooth, convex manifold described by the smooth distance function $r$ with $N=N(r)$ equilibria, while on the micro scale the particle's natural model is a finely discretized, convex polyhedral approximation $r^{\Delta}$ of $r$, with $N^{\Delta}=N(r^{\Delta})$ equilibria. There is strong intuitive evidence suggesting that under some particular evolution models $N(t)$ and $N^{\Delta}(t)$ primarily evolve in the opposite manner. Here we create the mathematical framework necessary to understand these phenomenon more broadly, regardless of the particular evolution equation. We study micro and macro events in one-parameter families of curves and surfaces, corresponding to bifurcations triggering the jumps in $N(t)$ and $N^{\Delta}(t)$. We show that the intuitive picture developed for curvature-driven flows is not only correct, it has universal validity, as long as the evolving surface $r$ is smooth. In this case, bifurcations associated with $r$ and $r^{\Delta}$ are coupled to some extent: resonance-like phenomena in $N^{\Delta}(t)$ can be used to forecast downward jumps in $N(t)$ (but not upward jumps). Beyond proving rigorous results for the $\Delta \to 0$ limit on the nontrivial interplay between singularities in the discrete and continuum approximations we also show that our mathematical model is structurally stable, i.e. it may be verified by computer simulations.


title: Superconductivity in Cage Compounds La$Tr_{2}$Al$_{20}$ with $Tr$ = Ti, V, Nb, and Ta

abstract: Electrical resistivity, magnetic susceptibility, and specific heat measurements on single crystals of La$Tr_{2}$Al$_{20}$ ($Tr$ = Ti, V, Nb, and Ta) revealed that these four compounds exhibit weak-coupling superconductivity with transition temperatures $T_{\rm c}$ = 0.46, 0.15, 1.05, and 1.03 K, respectively. LaTi$_{2}$Al$_{20}$ is most probably a type-I superconductor, which is quite rare among intermetallic compounds. Single-crystal X-ray diffraction suggests "rattling" anharmonic large-amplitude oscillations of Al ions (16$c$ site) on the Al$_{16}$ cage, while no such feature is suggested for the cage-center La ion. Using a parameter $d_{\rm GFS}$ quantifying the "guest free space" of the cage-center ion, we demonstrate that nonmagnetic $RTr_{2}$Al$_{20}$ superconductors are classified into two groups, i.e., (A) $d_{\rm GFS} \ne 0$ and $T_{\rm c}$ correlates with $d_{\rm GFS}$, and (B) $d_{\rm GFS} \simeq 0$ and $T_{\rm c}$ seems to be governed by other factors.


title: Simple Bounds for Utility Maximization with Small Transaction Costs

abstract: Using elementary arguments, we show how to derive $\mathbf{L}_p$-error bounds for the approximation of frictionless wealth process in markets with proportional transaction costs. For utilities with bounded risk aversion, these estimates yield lower bounds for the frictional value function, which pave the way for its asymptotic analysis using stability results for viscosity solutions. Using tools from Malliavin calculus, we also derive simple sufficient conditions for the regularity of frictionless optimal trading strategies, the second main ingredient for the asymptotic analysis of small transaction costs.


title: Contextuality and The Single-Qubit Stabilizer Formalism

abstract: Contextuality is a fundamental non-classical property of quantum theory, which has recently been proven to be a key resource for achieving quantum speed-ups in some leading models of quantum computation. However, which of the forms of contextuality, and how much thereof, are required to obtain a speed-up in an arbitrary model of quantum computation remains unclear. In this paper, we show that the relation between contextuality and a compuational advantage is more complicated than previously thought. We achieve this by proving that generalized contextuality is present even within the simplest subset of quantum operations, the so-called single-qubit stabilizer theory, which offers no computational advantage and was previously believed to be completely non-contextual. However, the contextuality of the single-qubit stabilizer theory can be confined to transformations. Therefore our result also demonstrates that the commonly considered prepare-and-measure scenarios (which ignore transformations) do not fully capture the contextuality of quantum theory.


title: The operator approach to the truncated multidimensional moment problem

abstract: We study the truncated multidimensional moment problem with a general type of truncations. The operator approach to the moment problem is presented. A way to construct atomic solutions of the moment problem is indicated.


title: SBP-SAT finite difference discretization of acoustic wave equations on staggered block-wise uniform grids

abstract: We consider the numerical simulation of the acoustic wave equations arising from seismic applications, for which staggered grid finite difference methods are popular choices due to their simplicity and efficiency. We relax the uniform grid restriction on finite difference methods and allow the grids to be block-wise uniform with nonconforming interfaces. In doing so, variations in the wave speeds of the subterranean media can be accounted for more efficiently. Staggered grid finite difference operators satisfying the summation-by-parts (SBP) property are devised to approximate the spatial derivatives appearing in the acoustic wave equation. These operators are applied within each block independently. The coupling between blocks is achieved through simultaneous approximation terms (SATs), which impose the interface condition weakly, i.e., by penalty. Ratio of the grid spacing of neighboring blocks is allowed to be rational number, for which specially designed interpolation formulas are presented. These interpolation formulas constitute key pieces of the simultaneous approximation terms. The overall discretization is shown to be energy-conserving and examined on test cases of both theoretical and practical interests, delivering accurate and stable simulation results.


title: Dirac operator on a noncommutative Toeplitz torus

abstract: We construct a 1+ summable regular even spectral triple for a noncommutative torus defined by a C*-subalgebra of the Toeplitz algebra.


title: Characterization of Type Ia Supernova Light Curves Using Principal Component Analysis of Sparse Functional Data

abstract: With growing data from ongoing and future supernova surveys it is possible to empirically quantify the shapes of SNIa light curves in more detail, and to quantitatively relate the shape parameters with the intrinsic properties of SNIa. Building such relationship is critical in controlling systematic errors associated with supernova cosmology. Based on a collection of well-observed SNIa samples accumulated in the past years, we construct an empirical SNIa light curve model using a statistical method called the functional principal component analysis (FPCA) for sparse and irregularly sampled functional data. Using this method, the entire light curve of an SNIa is represented by a linear combination of principal component functions, and the SNIa is represented by a few numbers called principal component scores. These scores are used to establish relations between light curve shapes and physical quantities such as intrinsic color, interstellar dust reddening, spectral line strength, and spectral classes. These relations allow for descriptions of some critical physical quantities based purely on light curve shape parameters. Our study shows that some important spectral feature information is being encoded in the broad band light curves, for instance, we find that the light curve shapes are correlated with the velocity and velocity gradient of the Si II $\lambda$6355 line. This is important for supernova surveys, e.g., LSST and WFIRST. Moreover, the FPCA light curve model is used to construct the entire light curve shape, which in turn is used in a functional linear form to adjust intrinsic luminosity when fitting distance models.


title: The Mean-Field Approximation: Information Inequalities, Algorithms, and Complexity

abstract: The mean field approximation to the Ising model is a canonical variational tool that is used for analysis and inference in Ising models. We provide a simple and optimal bound for the KL error of the mean field approximation for Ising models on general graphs, and extend it to higher order Markov random fields. Our bound improves on previous bounds obtained in work in the graph limit literature by Borgs, Chayes, Lov\'asz, S\'os, and Vesztergombi and another recent work by Basak and Mukherjee. Our bound is tight up to lower order terms. Building on the methods used to prove the bound, along with techniques from combinatorics and optimization, we study the algorithmic problem of estimating the (variational) free energy for Ising models and general Markov random fields. For a graph $G$ on $n$ vertices and interaction matrix $J$ with Frobenius norm $\| J \|_F$, we provide algorithms that approximate the free energy within an additive error of $\epsilon n \|J\|_F$ in time $\exp(poly(1/\epsilon))$. We also show that approximation within $(n \|J\|_F)^{1-\delta}$ is NP-hard for every $\delta > 0$. Finally, we provide more efficient approximation algorithms, which find the optimal mean field approximation, for ferromagnetic Ising models and for Ising models satisfying Dobrushin's condition.

categories: cs.LG cond-mat.stat-mech math-ph math.CO math.MP


title: K-theory and index pairings for C*-algebras generated by q-normal operators

abstract: The paper presents a detailed description of the K-theory and K-homology of C*-algebras generated by q-normal operators including generators and the index pairing. The C*-algebras generated by q-normal operators can be viewed as a q-deformation of the quantum complex plane. In this sense, we find deformations of the classical Bott projections describing complex line bundles over the 2-sphere, but there are also simpler generators for the K_0-groups, for instance 1-dimensional Powers-Rieffel type projections and elementary projections belonging to the C*-algebra. The index pairing between these projections and generators of the even K-homology group is computed, and the result is used to express the K_0-classes of the quantized line bundles of any winding number in terms of the other projections.


title: Topological edge states in the Su-Schrieffer-Heeger model subject to balanced particle gain and loss

abstract: We investigate the Su-Schrieffer-Heeger model in presence of an injection and removal of particles, introduced via a master equation in Lindblad form. It is shown that the dynamics of the density matrix follows the predictions of calculations in which the gain and loss are modeled by complex $\mathcal{PT}$-symmetric potentials. In particular it is found that there is a clear distinction in the dynamics between the topologically different cases known from the stationary eigenstates.


title: The Vertex Sample Complexity of Free Energy is Polynomial

abstract: We study the following question: given a massive Markov random field on $n$ nodes, can a small sample from it provide a rough approximation to the free energy $\mathcal{F}_n = \log{Z_n}$? Results in graph limit literature by Borgs, Chayes, Lov\'asz, S\'os, and Vesztergombi show that for Ising models on $n$ nodes and interactions of strength $\Theta(1/n)$, an $\epsilon$ approximation to $\log Z_n / n$ can be achieved by sampling a randomly induced model on $2^{O(1/\epsilon^2)}$ nodes. We show that the sampling complexity of this problem is {\em polynomial in} $1/\epsilon$. We further show a polynomial dependence on $\epsilon$ cannot be avoided. Our results are very general as they apply to higher order Markov random fields. For Markov random fields of order $r$, we obtain an algorithm that achieves $\epsilon$ approximation using a number of samples polynomial in $r$ and $1/\epsilon$ and running time that is $2^{O(1/\epsilon^2)}$ up to polynomial factors in $r$ and $\epsilon$. For ferromagnetic Ising models, the running time is polynomial in $1/\epsilon$. Our results are intimately connected to recent research on the regularity lemma and property testing, where the interest is in finding which properties can tested within $\epsilon$ error in time polynomial in $1/\epsilon$. In particular, our proofs build on results from a recent work by Alon, de la Vega, Kannan and Karpinski, who also introduced the notion of polynomial vertex sample complexity. Another critical ingredient of the proof is an effective bound by the authors of the paper relating the variational free energy and the free energy.


title: Fast, Trainable, Multiscale Denoising

abstract: Denoising is a fundamental imaging problem. Versatile but fast filtering has been demanded for mobile camera systems. We present an approach to multiscale filtering which allows real-time applications on low-powered devices. The key idea is to learn a set of kernels that upscales, filters, and blends patches of different scales guided by local structure analysis. This approach is trainable so that learned filters are capable of treating diverse noise patterns and artifacts. Experimental results show that the presented approach produces comparable results to state-of-the-art algorithms while processing time is orders of magnitude faster.


title: Design and Development of Effective Transmission Mechanisms on a Tendon Driven Hand Orthosis for Stroke Patients

abstract: Tendon-driven hand orthoses have advantages over exoskeletons with respect to wearability and safety because of their low-profile design and ability to fit a range of patients without requiring custom joint alignment. However, no existing study on a wearable tendon-driven hand orthosis for stroke patients presents evidence that such devices can overcome spasticity given repeated use and fatigue, or discusses transmission efficiency. In this study, we propose two designs that provide effective force transmission by increasing moment arms around finger joints. We evaluate the designs with geometric models and experiment using a 3D-printed artificial finger to find force and joint angle characteristics of the suggested structures. We also perform clinical tests with stroke patients to demonstrate the feasibility of the designs. The testing supports the hypothesis that the proposed designs efficiently elicit extension of the digits in patients with spasticity as compared to existing baselines.


title: Interaction Matters: A Note on Non-asymptotic Local Convergence of Generative Adversarial Networks

abstract: Motivated by the pursuit of a systematic computational and algorithmic understanding of Generative Adversarial Networks (GANs), we present a simple yet unified non-asymptotic local convergence theory for smooth two-player games, which subsumes several discrete-time gradient-based saddle point dynamics. The analysis reveals the surprising nature of the off-diagonal interaction term as both a blessing and a curse. On the one hand, this interaction term explains the origin of the slow-down effect in the convergence of Simultaneous Gradient Ascent (SGA) to stable Nash equilibria. On the other hand, for the unstable equilibria, exponential convergence can be proved thanks to the interaction term, for four modified dynamics proposed to stabilize GAN training: Optimistic Mirror Descent (OMD), Consensus Optimization (CO), Implicit Updates (IU) and Predictive Method (PM). The analysis uncovers the intimate connections among these stabilizing techniques, and provides detailed characterization on the choice of learning rate. As a by-product, we present a new analysis for OMD proposed in Daskalakis, Ilyas, Syrgkanis, and Zeng [2017] with improved rates.


title: Dual phase-space cascades in 3D hybrid-Vlasov-Maxwell turbulence

abstract: To explain energy dissipation via turbulence in collisionless, magnetized plasmas, the existence of a dual real- and velocity-space cascade of ion-entropy fluctuations below the ion gyroradius has been proposed. Such a dual cascade, predicted by the gyrokinetic theory, has previously been observed in gyrokinetic simulations of two-dimensional, electrostatic turbulence. For the first time we show evidence for a dual phase-space cascade of ion-entropy fluctuations in a three-dimensional simulation of hybrid-kinetic, electromagnetic turbulence. Some of the scalings observed in the energy spectra are consistent with a generalized theory for the cascade that accounts for the spectral anisotropy of critically balanced, intermittent, sub-ion-Larmor-scale fluctuations. The observed velocity-space cascade is also anisotropic with respect to the magnetic-field direction, with linear phase mixing along magnetic-field lines proceeding mainly at spatial scales above the ion gyroradius and nonlinear phase mixing across magnetic-field lines proceeding at perpendicular scales below the ion gyroradius. Such phase-space anisotropy could be sought in heliospheric and magnetospheric data of solar-wind turbulence and has far-reaching implications for the dissipation of turbulence in weakly collisional astrophysical plasmas.


title: Towards commissioning the Fermilab Muon G-2 Experiment

abstract: Starting this summer, Fermilab will host a key experiment dedicated to the search for signals of new physics: The Fermilab Muon g-2 Experiment. Its aim is to precisely measure the anomalous magnetic moment of the muon. In full operation, in order to avoid contamination, the newly born secondary beam is injected into a 505 m long Delivery Ring (DR) wherein it makes several revolutions before being sent to the experiment. Part of the commissioning scenario will execute a running mode wherein the passage from the DR will be skipped. With the aid of numerical simulations, we provide estimates of the expected performance.


title: Null Conservation Laws for Gravity

abstract: We give a full analysis of the conservation along null surfaces of generalized energy and super-momenta, for gravitational systems enclosed by a finite boundary. In particular we interpret the conservation equations in a canonical manner, revealing a notion of symplectic potential and a boundary current intrinsic to null surfaces. This generalizes similar analyses done at asymptotic infinity or on horizons.


title: Self-organized oscillations of Leidenfrost drops

abstract: In the Leidenfrost effect, a thin layer of evaporated vapor forms between a liquid and a hot solid. The complex interactions between the solid, liquid, and vapor phases can lead to rich dynamics even in a single Leidenfrost drop. Here we investigate the self-organized oscillations of Leidenfrost drops that are excited by a constant flow of evaporated vapor beneath the drop. We show that for small Leidenfrost drops, the frequency of a recently reported "breathing mode" can be explained by a simple balance of gravitational and surface tension forces. For large Leidenfrost drops, azimuthal star-shaped oscillations are observed. Our previous work showed how the coupling between the rapid evaporated vapor flow and the vapor-liquid interface excites the star oscillations (Ma \textit{et al., Phys. Rev. Fluids}, 2, 2017, 031602). In our experiments, star-shaped oscillation modes of $n=2$ to 13 are observed in different liquids, the number of observed modes depends sensitively on the liquid viscosity. Here we expand on this work by directly comparing the oscillations with theoretical predictions, and show how the oscillations are initiated by a parametric forcing mechanism through pressure oscillations in the vapor layer. The pressure oscillations are driven by the capillary waves of a characteristic wavelength beneath the drop. These capillary waves can be generated by a large shear stress at the liquid-vapor interface due to the rapid flow of evaporated vapor. We also explore potential effects of thermal convection in the liquid. Although the measured Rayleigh number is significantly larger than the critical value, the frequency (wavelength) of the oscillations depends only on the capillary length of the liquid, and is independent of the drop radius and substrate temperature. Thus convection seems to play a minor role in Leidenfrost drop oscillations, which are mostly hydrodynamic in origin.


title: A Unified Framework for Planning in Adversarial and Cooperative Environments

abstract: Users of AI systems may rely upon them to produce plans for achieving desired objectives. Such AI systems should be able to compute obfuscated plans whose execution in adversarial situations protects privacy, as well as legible plans which are easy for team members to understand in cooperative situations. We develop a unified framework that addresses these dual problems by computing plans with a desired level of comprehensibility from the point of view of a partially informed observer. For adversarial settings, our approach produces obfuscated plans with observations that are consistent with at least k goals from a set of decoy goals. By slightly varying our framework, we present an approach for goal legibility in cooperative settings which produces plans that achieve a goal while being consistent with at most j goals from a set of confounding goals. In addition, we show how the observability of the observer can be controlled to either obfuscate or clarify the next actions in a plan when the goal is known to the observer. We present theoretical results on the complexity analysis of our problems. We demonstrate the execution of obfuscated and legible plans in a cooking domain using a physical robot Fetch. We also provide an empirical evaluation to show the feasibility and usefulness of our approaches using IPC domains.


title: Detecting Social Influence in Event Cascades by Comparing Discriminative Rankers

abstract: The global dynamics of event cascades are often governed by the local dynamics of peer influence. However, detecting social influence from observational data is challenging due to confounds like homophily and practical issues like missing data. We propose a simple discriminative method to detect influence from observational data. The core of the approach is to train a ranking algorithm to predict the source of the next event in a cascade, and compare its out-of-sample accuracy against a competitive baseline which lacks access to features corresponding to social influence. We analyze synthetically generated data to show that this method correctly identifies influence in the presence of confounds, and is robust to both missing data and misspecification --- unlike well-known alternatives. We apply the method to two real-world datasets: (1) the co-sponsorship of legislation in the U.S. House of Representatives on a social network of shared campaign donors; (2) rumors about the Higgs boson discovery on a follower network of $10^5$ Twitter accounts. Our model identifies the role of social influence in these scenarios and uses it to make more accurate predictions about the future trajectory of cascades.


title: Reactive Reinforcement Learning in Asynchronous Environments

abstract: The relationship between a reinforcement learning (RL) agent and an asynchronous environment is often ignored. Frequently used models of the interaction between an agent and its environment, such as Markov Decision Processes (MDP) or Semi-Markov Decision Processes (SMDP), do not capture the fact that, in an asynchronous environment, the state of the environment may change during computation performed by the agent. In an asynchronous environment, minimizing reaction time---the time it takes for an agent to react to an observation---also minimizes the time in which the state of the environment may change following observation. In many environments, the reaction time of an agent directly impacts task performance by permitting the environment to transition into either an undesirable terminal state or a state where performing the chosen action is inappropriate. We propose a class of reactive reinforcement learning algorithms that address this problem of asynchronous environments by immediately acting after observing new state information. We compare a reactive SARSA learning algorithm with the conventional SARSA learning algorithm on two asynchronous robotic tasks (emergency stopping and impact prevention), and show that the reactive RL algorithm reduces the reaction time of the agent by approximately the duration of the algorithm's learning update. This new class of reactive algorithms may facilitate safer control and faster decision making without any change to standard learning guarantees.


title: Real-Time 3D Shape of Micro-Details

abstract: Motivated by the growing demand for interactive environments, we propose an accurate real-time 3D shape reconstruction technique. To provide a reliable 3D reconstruction which is still a challenging task when dealing with real-world applications, we integrate several components including (i) Photometric Stereo (PS), (ii) perspective Cook-Torrance reflectance model that enables PS to deal with a broad range of possible real-world object reflections, (iii) realistic lightening situation, (iv) a Recurrent Optimization Network (RON) and finally (v) heuristic Dijkstra Gaussian Mean Curvature (DGMC) initialization approach. We demonstrate the potential benefits of our hybrid model by providing 3D shape with highly-detailed information from micro-prints for the first time. All real-world images are taken by a mobile phone camera under a simple setup as a consumer-level equipment. In addition, complementary synthetic experiments confirm the beneficial properties of our novel method and its superiority over the state-of-the-art approaches.


title: A generic characterization of Pol(C)

abstract: We investigate the polynomial closure operation (C -> Pol(C)) defined on classes of regular languages. We present an interesting and useful connection relating the separation problem for the class C and the membership problem for it polynomial closure Pol(C). This connection is formulated as an algebraic characterization of Pol(C) which holds when C is an arbitrary \pvari of regular languages and whose statement is parameterized by C-separation. Its main application is an effective reduction from Pol(C)-membership to C-separation. Thus, as soon as one designs a C-separation algorithm, this yields "for free" a membership algorithm for the more complex class Pol(C).


title: A noncommutative 2-sphere generated by the quantum complex plane

abstract: S. L. Woronowicz's theory of introducing C*-algebras generated by unbounded elements is applied to q-normal operators satisfying the defining relation of the quantum complex plane. The unique non-degenerate C*-algebra of bounded operators generated by a q-normal operator is computed and an abstract description is given by using crossed product algebras. If the spectrum of the modulus of the q-normal operator is the positive half line, this C*-algebra will be considered as the algebra of continuous functions on the quantum complex plane vanishing at infinity, and its unitization will be viewed as the algebra of continuous functions on a quantum 2-sphere.


title: On the Tur\'an density of $\{1, 3\}$-Hypergraphs

abstract: In this paper, we consider the Tur\'an problems on $\{1,3\}$-hypergraphs. We prove that a $\{1, 3\}$-hypergraph is degenerate if and only if it's $H^{\{1, 3\}}_5$-colorable, where $H^{\{1, 3\}}_5$ is a hypergraph with vertex set $V=[5]$ and edge set $E=\{\{2\}, \{3\}, \{1, 2, 4\}, \{1, 3, 5\}, \{1, 4, 5\}\}.$ Using this result, we further prove that for any finite set $R$ of distinct positive integers, except the case $R=\{1, 2\}$, there always exist non-trivial degenerate $R$-graphs. We also compute the Tur\'an densities of some small $\{1,3\}$-hypergraphs.


title: Function algebras on a 2-dimensional quantum complex plane

abstract: The well-behaved representations of the coordinate algebra of a 2-dimensional quantum complex plane are classified and a C*-algebra is defined which can be viewed as the algebra of continuous functions on the 2-dimensional quantum complex plane vanishing at infinity.


title: Selmer groups are intersection of two direct summands of the adelic cohomology

abstract: We give a positive answer to a Conjecture by Manjul Bhargava, Daniel M. Kane, Hendrik W. Lenstra Jr., Bjorn Poonen and Eric Rains, concerning the cohomology of torsion subgroups of elliptic curves over global fields. This implies that, given a global field $k$ and an integer $n$, for $100\%$ of elliptic curves $E$ defined over $k$, the $n$-th Selmer group of $E$ is the intersection of two direct summands of the adelic cohomology group $H^1(\mathbf{A},E[n])$. We also give examples of elliptic curves for which the conclusion of this conjecture does not hold.


title: Representations of quantum SU(2) operators on a local chart

abstract: Hilbert space representations of quantum SU(2) by multiplication operators on a local chart are constructed, where the local chart is given by tensor products of square integrable functions on a quantum disc and on the classical unit circle. The actions of generators of quantum SU(2), generators of the opposite algebra, and noncommutative partial derivatives are computed on a Hilbert space basis.


title: Doubling inequality and nodal sets for solutions of bi-Laplace equations

abstract: We investigate the doubling inequality and nodal sets for the solutions of bi-Laplace equations. A polynomial upper bound for the nodal sets of solutions and their gradient is obtained based on the recent development of nodal sets for Laplace eigenfunctions by Logunov. In addition, we derive an implicit upper bound for the nodal sets of solutions. We show two types of doubling inequalities for the solutions of bi-Laplace equations. As a consequence, the rate of vanishing is given for the solutions.


title: Multi-user Beam-Alignment for Millimeter-Wave Networks

abstract: Millimeter-wave communications is the most promising technology for next-generation cellular wireless systems, thanks to the large bandwidth available compared to sub-6 GHz networks. Nevertheless, communication at these frequencies requires narrow beams via massive MIMO and beamforming to overcome the strong signal attenuation, and thus precise beam-alignment between transmitter and receiver is needed. The resulting signaling overhead may become a severe impairment, especially in mobile networks with high users density. Therefore, it is imperative to optimize the beam-alignment protocol to minimize the signaling overhead. In this paper, the design of energy efficient joint beam-alignment protocols for two users is addressed, with the goal to minimize the power consumption during data transmission, subject to rate constraints for both users, under analog beamforming constraints. It is proved that a bisection search algorithm is optimal. Additionally, the optimal scheduling strategy of the two users in the data communication phase is optimized based on the outcome of beam-alignment, according to a time division multiplexing scheme. The numerical results show significant decrease in the power consumption for the proposed joint beam-alignment scheme compared to exhaustive search and a single-user beam-alignment scheme taking place separately for each user.


title: Search for narrow resonances in the b-tagged dijet mass spectrum in proton-proton collisions at $\sqrt{s} =$ 8 TeV

abstract: A search for narrow resonances decaying to bottom quark-antiquark pairs is presented, using a data sample of proton-proton collisions at $\sqrt{s} =$ 8 TeV corresponding to an integrated luminosity of 19.7 fb$^{-1}$. The search is extended to masses lower than those reached in typical searches for resonances decaying into jet pairs at the LHC, by taking advantage of triggers that identify jets originating from bottom quarks. No significant excess of events is observed above the background predictions. Limits are set on the product of cross section and branching fraction to bottom quarks for spin 0, 1, and 2 resonances in the mass range of 325-1200 GeV. These results significantly improve on the limits for resonances decaying into jet pairs in the 325-500 GeV mass range.


title: Robust s+/- pairing in CaK[Fe(1-x)Ni(x)]4As4$ (x = 0 and 0.05) from the response to electron irradiation

abstract: Controlled point-like disorder introduced by 2.5 MeV electron irradiation was used to probe the superconducting state of single crystals of \CaKx\ superconductor at $x = 0$ and 0.05 doping levels. Both compositions show an increase of the residual resistivity and a decrease of the superconducting transition temperature, $T_c$ at the rate of $dT_c/d\rho(T_c) \approx$ 0.19 K(\textmu$\Omega$cm)$^{-1}$ for $x=0$ and 0.38 K(\textmu$\Omega$cm)$^{-1}$ for $x=\:$0.05, respectively. In Ni - doped, $x = 0.05$, compound the coexisting spin-vortex crystal (SVC) magnetic phase is suppressed at the rate of $dT_N/d\rho(T_N)\approx$ 0.16 K(\textmu$\Omega$cm)$^{-1}$. Low - temperature variation of London penetration depth is well approximated by the power law, $\Delta \lambda (T) = AT^n$ with $n\approx\,$2.5 for $x=0$ and $n\approx\,$1.9 for $x=0.05$ in the pristine state. Electron irradiation leads to the exponent $n$ increase above 2 in $x=0.05$ suggesting superconducting gap with significant anisotropy that is smeared by the disorder scattering. Detailed analysis of $\lambda (T)$ and \(T_{c}\) evolution with disorder is consistent with two effective nodeless superconducting energy gaps due to robust s$_{\pm}$ pairing. Overall the behavior of \CaKx\ at $x = 0$ is similar to a slightly overdoped \BaK\ at $y \approx$ 0.5 and at $x= 0.05$ to an underdoped composition at $y \approx$ 0.2.


title: Scalable Inference for Space-Time Gaussian Cox Processes

abstract: The log-Gaussian Cox process is a flexible and popular class of point pattern models for capturing spatial and space-time dependence for point patterns. Model fitting requires approximation of stochastic integrals which is implemented through discretization over the domain of interest. With fine scale discretization, inference based on Markov chain Monte Carlo is computationally burdensome because of the cost of matrix decompositions and storage, such as the Cholesky, for high dimensional covariance matrices associated with latent Gaussian variables. This article addresses these computational bottlenecks by combining two recent developments: (i) a data augmentation strategy that has been proposed for space-time Gaussian Cox processes that is based on exact Bayesian inference and does not require fine grid approximations for infinite dimensional integrals, and (ii) a recently developed family of sparsity-inducing Gaussian processes, called nearest-neighbor Gaussian processes, to avoid expensive matrix computations. Our inference is delivered within the fully model-based Bayesian paradigm and does not sacrifice the richness of traditional log-Gaussian Cox processes. We apply our method to crime event data in San Francisco and investigate the recovery of the intensity surface.


title: Plasma diagnostics of coronal dimming events

abstract: Coronal mass ejections (CMEs) are often associated with coronal dimmings, i.e. transient dark regions that are most distinctly observed in Extreme Ultra-violet (EUV) wavelengths. Using Atmospheric Imaging Assembly (AIA) data, we apply Differential Emission Measure (DEM) diagnostics to study the plasma characteristics of six coronal dimming events. In the core dimming region, we find a steep and impulsive decrease of density with values up to 50-70%. Five of the events also reveal an associated drop in temperature of 5-25%. The secondary dimming regions also show a distinct decrease in density, but less strong, decreasing by 10-45%. In both the core and the secondary dimming the density changes are much larger than the temperature changes, confirming that the dimming regions are mainly caused by plasma evacuation. In the core dimming, the plasma density reduces rapidly within the first 20-30 min after the flare start, and does not recover for at least 10 hrs later, whereas the secondary dimming tends to be more gradual and starts to replenish after 1-2 hrs. The pre-event temperatures are higher in the core dimming (1.7-2.6 MK) than in the secondary dimming regions (1.6-2.0 MK). Both core and secondary dimmings are best observed in the AIA 211 \AA\ and 193 \AA\ filters. These findings suggest that the core dimming corresponds to the footpoints of the erupting flux rope rooted in the AR, while the secondary dimming represents plasma from overlying coronal structures that expand during the CME eruption.


title: A Likelihood-Free Inference Framework for Population Genetic Data using Exchangeable Neural Networks

abstract: An explosion of high-throughput DNA sequencing in the past decade has led to a surge of interest in population-scale inference with whole-genome data. Recent work in population genetics has centered on designing inference methods for relatively simple model classes, and few scalable general-purpose inference techniques exist for more realistic, complex models. To achieve this, two inferential challenges need to be addressed: (1) population data are exchangeable, calling for methods that efficiently exploit the symmetries of the data, and (2) computing likelihoods is intractable as it requires integrating over a set of correlated, extremely high-dimensional latent variables. These challenges are traditionally tackled by likelihood-free methods that use scientific simulators to generate datasets and reduce them to hand-designed, permutation-invariant summary statistics, often leading to inaccurate inference. In this work, we develop an exchangeable neural network that performs summary statistic-free, likelihood-free inference. Our framework can be applied in a black-box fashion across a variety of simulation-based tasks, both within and outside biology. We demonstrate the power of our approach on the recombination hotspot testing problem, outperforming the state-of-the-art.


title: InP femtosecond mode-locked laser in a compound feedback cavity with a switchable repetition rate

abstract: A monolithically integrated mode-locked semiconductor laser is proposed. The compound ring cavity is composed of a colliding pulse mode-locking (ML) subcavity and a passive Fabry-Perot feedback subcavity. These two 1.6 mm long subcavities are coupled by using on-chip reflectors at both ends, enabling harmonic mode locking. By changing DC-bias conditions, optical mode spacing from 50 to 450 GHz is experimentally demonstrated. Ultrafast pulses shorter than 0.3 ps emitted from this laser diode are shown in autocorrelation traces.


title: Oriented Borel-Moore homologies of toric varieties

abstract: We generalize the K\"unneth formula for Chow groups to an arbitrary OBM-homology theory satisfying descent (e.g. algebraic cobordism) when taking a product with a toric variety. As a corollary we obtain a universal coefficient theorem for the operational cohomology rings. We also give a description for the homology groups of smooth toric varieties, that allows the calculation of algebraic cobordism groups of singular toric varieties. Some computations are carried out.